title: Magnetic and charge orders in the ground state of the Emery model -- accurate numerical results

abstract: We perform extensive auxiliary-field quantum Monte Carlo (AFQMC) calculations for the three-band Hubbard (Emery) model in order to study the ground-state properties of Copper-Oxygen planes in the cuprates. Employing cutting-edge AFQMC techniques with a self-consistent gauge constraint in auxiliary-field space to control the sign problem, we reach supercells containing around 500 atoms to capture collective modes in the charge and spin orders and characterize the behavior in the thermodynamic limit. The self-consistency scheme interfacing with generalized Hartree-Fock calculations allows high accuracy in AFQMC to resolve small energy scales, which is crucial for determining the complex candidate orders in such a system. We present detailed information on the charge order, spin order, momentum distribution, and localization properties as a function of charge-transfer energy for the the under-doped regime. In contrast with the stripe and spiral orders under hole-doping, we find that the corresponding 1/8 electron-doped system exhibits purely antiferromagnetic order in the three-band model, consistent with the asymmetry between electron and hole-doping in the phase diagram of cuprates.


title: Flat band topology of magic angle graphene on a transition metal dichalcogenide

abstract: We consider twisted bilayer graphene on a transition metal dichalcogenide substrate, where proximity-induced spin-orbit coupling significantly alters the eight flat bands which occur near the magic angle. The resulting band structure features a pair of extremely flat bands across most of the mini-Brillouin zone. Further details depend sensitively on the symmetries of the heterostructure; we find semiconducting band structures when all two-fold rotations around in-plane axis are broken, and semi-metallic band structures otherwise. We calculate the Chern numbers of the different isolated bands, and identify the parameter regimes and filling factors where valley Chern insulators and topological insulators are realized. Interestingly, we find that for realistic values of the proximity-induced terms, it is possible to realize a topological insulator protected by time-reversal symmetry by doping two holes or two electrons per superlattice unit cell into the system.


title: Three Variants of Differential Privacy: Lossless Conversion and Applications

abstract: We consider three different variants of differential privacy (DP), namely approximate DP, R\'enyi DP (RDP), and hypothesis test DP. In the first part, we develop a machinery for optimally relating approximate DP to RDP based on the joint range of two $f$-divergences that underlie the approximate DP and RDP. In particular, this enables us to derive the optimal approximate DP parameters of a mechanism that satisfies a given level of RDP. As an application, we apply our result to the moments accountant framework for characterizing privacy guarantees of noisy stochastic gradient descent (SGD). When compared to the state-of-the-art, our bounds may lead to about 100 more stochastic gradient descent iterations for training deep learning models for the same privacy budget. In the second part, we establish a relationship between RDP and hypothesis test DP which allows us to translate the RDP constraint into a tradeoff between type I and type II error probabilities of a certain binary hypothesis test. We then demonstrate that for noisy SGD our result leads to tighter privacy guarantees compared to the recently proposed $f$-DP framework for some range of parameters.


title: On Explaining the Surprising Success of Reservoir Computing Forecaster of Chaos? The Universal Machine Learning Dynamical System with Contrasts to VAR and DMD

abstract: Machine learning has become a widely popular and successful paradigm, including in data-driven science and engineering. A major application problem is data-driven forecasting of future states from a complex dynamical. Artificial neural networks (ANN) have evolved as a clear leader amongst many machine learning approaches, and recurrent neural networks (RNN) are considered to be especially well suited for forecasting dynamical systems. In this setting, the echo state networks (ESN) or reservoir computer (RC) have emerged for their simplicity and computational complexity advantages. Instead of a fully trained network, an RC trains only read-out weights by a simple, efficient least squares method. What is perhaps quite surprising is that nonetheless an RC succeeds to make high quality forecasts, competitively with more intensively trained methods, even if not the leader. There remains an unanswered question as to why and how an RC works at all, despite randomly selected weights. We explicitly connect the RC with linear activation and linear read-out to well developed time-series literature on vector autoregressive averages (VAR) that includes theorems on representability through the WOLD theorem, which already perform reasonably for short term forecasts. In the case of a linear activation and now popular quadratic read-out RC, we explicitly connect to a nonlinear VAR (NVAR), which performs quite well. Further, we associate this paradigm to the now widely popular dynamic mode decomposition (DMD), and thus these three are in a sense different faces of the same thing. We illustrate our observations in terms of popular benchmark examples including Mackey-Glass differential delay equations and the Lorenz63 system.


title: The Average Order of Dominating Sets of a Graph

abstract: This papers focuses on the average order of dominating sets of a graph. We find the extremal graphs for the maximum and minimum value over all graphs on $n$ vertices, while for trees we prove that the star minimizes the average order of dominating sets. We prove the average order of dominating sets in graphs without isolated vertices is at most $3n/4$, but provide evidence that the actual upper bound is $2n/3$. Finally, we show that the normalized average, while dense in $[1/2,1]$, tends to $\frac{1}{2}$ for almost all graphs.


title: Equilibrium frame reveals hidden PT symmetry of passive systems

abstract: We discuss how introducing an equilibrium frame, in which a given Hamiltonian has balanced loss and gain terms, can reveal PT symmetry hidden in non-Hermitian Hamiltonians of dissipative systems. Passive PT-symmetric Hamiltonians, in which only loss is present and gain is absent, can also display exceptional points, just like PT-symmetric systems, and therefore are extensively investigated. We demonstrate that non-Hermitian Hamiltonians, which can be divided into a PT-symmetric term and a term commuting with the Hamiltonian, possess hidden PT symmetries. These symmetries become apparent in the equilibrium frame. We also show that the number of eigenstates having the same value in an exceptional point is usually smaller in the initial frame than in the equilibrium frame. This property is associated with the second part of the Hamiltonian.


title: Photonic simulation of giant atom decay

abstract: Spontaneous emission of an excited atom in a featureless continuum of electromagnetic modes is a fundamental process in quantum electrodynamics associated with an exponential decay of the quantum emitter to its ground state accompanied by an irreversible emission of a photon. However, such a simple scenario is deeply modified when considering a $^{\prime}$giant$^{\prime}$ atom, i.e an atom whose dimension is larger than the wavelength of the emitted photon. In such an unconventional regime, non-Markovian effects and strong deviations from an exponential decay are observed owing to interference effects arising from non-local light-atom coupling. Here we suggest a photonic simulation of non-Markovian giant atom decay, based on light escape dynamics in an optical waveguide non-locally-coupled to a waveguide lattice. Major effects such as non-exponential decay, enhancement or slowing down of the decay, and formation of atom-field dark states can be emulated in this system


title: MatryODShka: Real-time 6DoF Video View Synthesis using Multi-Sphere Images

abstract: We introduce a method to convert stereo 360{\deg} (omnidirectional stereo) imagery into a layered, multi-sphere image representation for six degree-of-freedom (6DoF) rendering. Stereo 360{\deg} imagery can be captured from multi-camera systems for virtual reality (VR), but lacks motion parallax and correct-in-all-directions disparity cues. Together, these can quickly lead to VR sickness when viewing content. One solution is to try and generate a format suitable for 6DoF rendering, such as by estimating depth. However, this raises questions as to how to handle disoccluded regions in dynamic scenes. Our approach is to simultaneously learn depth and disocclusions via a multi-sphere image representation, which can be rendered with correct 6DoF disparity and motion parallax in VR. This significantly improves comfort for the viewer, and can be inferred and rendered in real time on modern GPU hardware. Together, these move towards making VR video a more comfortable immersive medium.


title: Infrared color separation between thin-shelled oxygen-rich and carbon-rich AGB stars

abstract: We present 43 GHz VLA spectra for 51 AGB sources with the goal of verifying an infrared (IR) color cut intended to separate Carbon-rich (C) and Oxygen-rich (O) AGB sources throughout the Galaxy. The color cut is a simple line in the $[K_s]-[A]$ versus $[A]-[E]$ color-color diagram based on 2MASS and MSX photometry, and was originally derived from SiO detection rates in the Bulge Asymmetries and Dynamical Evolution (BAaDE) sample. The division is fully supported by the spectra presented here, which show that SiO maser detections lie on the O-rich side, and SiO nondetections and a single HC$_3$N detection are found on the C-rich side of the division. We further compare the color cut with classifications of the sources based on Low-Resolution Spectra (LRS) from the Infrared Astronomical Satellite (IRAS), and find good agreement, verifying that the division is a reliable and efficient method for differentiating O- and C-rich AGB sources. These observations also demonstrate that single lines detected in the BAaDE survey around 42.9 GHz are almost certainly the $^{29}$SiO v=0 line. SiO maser sources where this rare isotopologue transition is brighter than the dominant $^{28}$SiO lines have not been reported before, and our observations show that these sources can reverse their behavior such that the typical ratios of $^{28}$SiO and $^{29}$SiO are restored within a few years.


title: Making Distributed Mobile Applications SAFE: Enforcing User Privacy Policies on Untrusted Applications with Secure Application Flow Enforcement

abstract: Today's mobile devices sense, collect, and store huge amounts of personal information, which users share with family and friends through a wide range of applications. Once users give applications access to their data, they must implicitly trust that the apps correctly maintain data privacy. As we know from both experience and all-too-frequent press articles, that trust is often misplaced. While users do not trust applications, they do trust their mobile devices and operating systems. Unfortunately, sharing applications are not limited to mobile clients but must also run on cloud services to share data between users. In this paper, we leverage the trust that users have in their mobile OSes to vet cloud services. To do so, we define a new Secure Application Flow Enforcement (SAFE) framework, which requires cloud services to attest to a system stack that will enforce policies provided by the mobile OS for user data. We implement a mobile OS that enforces SAFE policies on unmodified mobile apps and two systems for enforcing policies on untrusted cloud services. Using these prototypes, we demonstrate that it is possible to enforce existing user privacy policies on unmodified applications.


title: The Relevance of Classic Fuzz Testing: Have We Solved This One?

abstract: As fuzz testing has passed its 30th anniversary, and in the face of the incredible progress in fuzz testing techniques and tools, the question arises if the classic, basic fuzz technique is still useful and applicable? In that tradition, we have updated the basic fuzz tools and testing scripts and applied them to a large collection of Unix utilities on Linux, FreeBSD, and MacOS. As before, our failure criteria was whether the program crashed or hung. We found that 9 crash or hang out of 74 utilities on Linux, 15 out of 78 utilities on FreeBSD, and 12 out of 76 utilities on MacOS. A total of 24 different utilities failed across the three platforms. We note that these failure rates are somewhat higher than our in previous 1995, 2000, and 2006 studies of the reliability of command line utilities. In the basic fuzz tradition, we debugged each failed utility and categorized the causes the failures. Classic categories of failures, such as pointer and array errors and not checking return codes, were still broadly present in the current results. In addition, we found a couple of new categories of failures appearing. We present examples of these failures to illustrate the programming practices that allowed them to happen. As a side note, we tested the limited number of utilities available in a modern programming language (Rust) and found them to be of no better reliability than the standard ones.


title: Uncertainty Estimates of Solar Wind Prediction using HMI Photospheric Vector and Spatial Standard Deviation Synoptic Maps

abstract: Current solar wind prediction is based on the Wang & Sheeley empirical relationship between the solar wind speed observed at 1 AU and the rate of magnetic flux tube expansion (FTE) between the photosphere and the inner corona, where FTE is computed by coronal models that take the photospheric flux density synoptic maps as their inner boundary conditions to extrapolate the photospheric magnetic fields to deduce the coronal and the heliospheric magnetic field configuration. Since these synoptic maps are among the most widely-used of all solar magnetic data products, the uncertainties in the model predictions that are caused by the uncertainties in the synoptic maps are worthy of study. However, such an estimate related to synoptic map construction was not available until Bertello et al. (Solar Physics, 289, 2014) obtained the spatial standard deviation synoptic maps; 98 Monte-Carlo realizations of the spatial standard deviation maps for each photospheric synoptic maps. In this paper, we present an estimate of uncertainties in the solar wind speed predicted at 1 AU by the CSSS model due to the uncertainties in the photospheric synoptic maps. We also present a comparison of the coronal hole locations predicted by the models with the STEREO/SECCHI EUV synoptic maps. In order to quantify the extent of the uncertainties involved, we compared the predicted speeds with the OMNI solar wind data during the same period (taking the solar wind transit time into account) and obtained the root mean square error between them. To illustrate the significance of the uncertainty estimate in the solar wind prediction, we carried out the analysis for three Carrington rotations, CR 2102, CR 2137 and CR 2160 at different phases of the solar cycle. The uncertainty estimate is critical information necessary for the current and future efforts of improving the solar wind prediction accuracies.


title: Rotation-time symmetry in bosonic systems and the existence of exceptional points in the absence of $\mathcal{PT}$ symmetry

abstract: We study symmetries of open bosonic systems in the presence of laser pumping. Non-Hermitian Hamiltonians describing these systems can be parity-time (${\cal{PT}}$) symmetric in special cases only. Systems exhibiting this symmetry are characterised by real-valued energy spectra and can display exceptional points, where a symmetry-breaking transition occurs. We demonstrate that there is a more general type of symmetry, i.e., rotation-time (${\cal{RT}}$) symmetry. We observe that ${\cal{RT}}$-symmetric non-Hermitian Hamiltonians exhibit real-valued energy spectra which can be made singular by symmetry breaking. To calculate the spectra of the studied bosonic non-diagonalisable Hamiltonians we apply diagonalisation methods based on bosonic algebra. Finally, we list a versatile set rules allowing to immediately identifying or constructing ${\cal{RT}}$-symmetric Hamiltonians. We believe that our results on the ${\cal{RT}}$-symmetric class of bosonic systems and their spectral singularities can lead to new applications inspired by those of the ${\cal{PT}}$-symmetric systems.


title: Superradiance paradox in waveguide lattices

abstract: Recently, it has been suggested that the collective radiative decay of two point-like quantum emitters coupled to a waveguide, separated by a distance comparable to the coherence length of a spontaneously emitted photon, leads to an apparent $^{\prime}$superradiance paradox$^{\prime}$ by which one cannot decide whether independent or collective emission occurs. The resolution of the paradox stems from the strong non-Markovian dynamics arising from the delayed field-mediated atom interaction. Here we suggest an integrated optics platform to emulate the superradiance paradox, based on photon escape dynamics in waveguide lattices. Remarkably, Markovian decay dynamics and independent photon emission can be restored by frequent (Zeno-like) observation of the system.


title: Biasing relation, environmental dependencies and estimation of the growth rate from star forming galaxies

abstract: The connection between galaxy star formation rate (SFR) and dark matter (DM) is of paramount importance for the extraction of cosmological information from next generation spectroscopic surveys that will target emission line star forming galaxies. Using publicly available mock galaxy catalogs obtained from various semi-analytic models (SAMs) we explore the SFR-DM connection in relation to the speed-from-light method (Feix et al. 2016) for inferring the growth rate, $f$, from luminosity/SFR shifts. Emphasis is given to the dependence of the SFR distribution on the environment density on scales of 10s-100s Mpc. We show that the application of the speed-from-light method to an Euclid-like survey is not biased by environmental effects. In all models, the precision on the measured $\beta=f/b$ parameter is $\sigma_\beta < 0.17$ at $z=1$. This translates into errors of $\sigma_f \sim 0.22$ and $\sigma_{(f\sigma_8)}\sim 0.1$, without invoking assumptions on the mass power spectrum. These errors are in the same ballpark as recent analyses of the redshift space distortions in galaxy clustering. In agreement with previous studies, the bias factor, $b$ is roughly a scale-independent, constant function of the SFR for star forming galaxies. Its value at $z=1$ ranges from $1.2$ to $1.5$ depending on the SAM recipe. Although in all SAMs denser environments host galaxies with higher stellar masses, the dependence of the SFR on the environment is more involved. In most models the SFR probability distribution is skewed to larger values in denser regions. One model exhibits an inverted trend where high SFR is suppressed in dense environment.


title: A Scalable, Adaptive and Sound Nonconvex Regularizer for Low-rank Matrix Completion

abstract: Matrix learning is at the core of many machine learning problems. A number of real-world applications such as collaborative filtering and text mining can be formulated as a low-rank matrix completion problem, which recovers incomplete matrix using low-rank assumptions. To ensure that the matrix solution has a low rank, a recent trend is to use nonconvex regularizers that adaptively penalize singular values. They offer good recovery performance and have nice theoretical properties, but are computationally expensive due to repeated access to individual singular values. In this paper, based on the key insight that adaptive shrinkage on singular values improve empirical performance, we propose a new nonconvex low-rank regularizer called "nuclear norm minus Frobenius norm" regularizer, which is scalable, adaptive and sound. We first show it provably holds the adaptive shrinkage property. Further, we discover its factored form which bypasses the computation of singular values and allows fast optimization by general optimization algorithms. Stable recovery and convergence are guaranteed. Extensive low-rank matrix completion experiments on a number of synthetic and real-world data sets show that the proposed method obtains state-of-the-art recovery performance while being the fastest in comparison to existing low-rank matrix learning methods.


title: AntiDote: Attention-based Dynamic Optimization for Neural Network Runtime Efficiency

abstract: Convolutional Neural Networks (CNNs) achieved great cognitive performance at the expense of considerable computation load. To relieve the computation load, many optimization works are developed to reduce the model redundancy by identifying and removing insignificant model components, such as weight sparsity and filter pruning. However, these works only evaluate model components' static significance with internal parameter information, ignoring their dynamic interaction with external inputs. With per-input feature activation, the model component significance can dynamically change, and thus the static methods can only achieve sub-optimal results. Therefore, we propose a dynamic CNN optimization framework in this work. Based on the neural network attention mechanism, we propose a comprehensive dynamic optimization framework including (1) testing-phase channel and column feature map pruning, as well as (2) training-phase optimization by targeted dropout. Such a dynamic optimization framework has several benefits: (1) First, it can accurately identify and aggressively remove per-input feature redundancy with considering the model-input interaction; (2) Meanwhile, it can maximally remove the feature map redundancy in various dimensions thanks to the multi-dimension flexibility; (3) The training-testing co-optimization favors the dynamic pruning and helps maintain the model accuracy even with very high feature pruning ratio. Extensive experiments show that our method could bring 37.4% to 54.5% FLOPs reduction with negligible accuracy drop on various of test networks.


title: Reinvestigation of the homogeneous spin model in YbMgGaO$_4$

abstract: Motivated by a recent inelastic neutron scattering experiment on $\mathrm{YbMgGaO}_4$ \cite{William2019}, we reinvestigate the homogeneous spin model on the triangular lattice. Using the cluster mean-field theory, we study the phase diagram and the magnetic-field-induced phase transition. We find that the phase boundary between the stripe state and the $120^{\circ}$ antiferromagnetic state is broadened by the magnetic field, leading to a field-induced phase transition. This phase transition is suppressed by the next-nearest neighbor exchange interaction $J_2/J_1$ and vanishes as $J_2/J_1>0.13$. We find a parameter space at $J_2/J_1=0.1$, in which the field-induce transition can be achieved and the deviation of theoretical spin excitation energies from experimental data is only $5.4\%$. Our results imply that an effective homogeneous spin model still works in $\mathrm{YbMgGaO}_4$.


title: GANplifying Event Samples

abstract: A critical question concerning generative networks applied to event generation in particle physics is if the generated events add statistical precision beyond the training sample. We show for a simple example with increasing dimensionality how generative networks indeed amplify the training statistics. We quantify their impact through an amplification factor or equivalent numbers of sampled events.


title: Learning Lyapunov Functions for Piecewise Affine Systems with Neural Network Controllers

abstract: We propose a learning-based method for Lyapunov stability analysis of piecewise affine dynamical systems in feedback with piecewise affine neural network controllers. The proposed method consists of an iterative interaction between a learner and a verifier, where in each iteration, the learner uses a collection of samples of the closed-loop system to propose a Lyapunov function candidate as the solution to a convex program. The learner then queries the verifier, which solves a mixed-integer program to either validate the proposed Lyapunov function candidate or reject it with a counterexample, i.e., a state where the stability condition fails. This counterexample is then added to the sample set of the learner to refine the set of Lyapunov function candidates. We design the learner and the verifier based on the analytic center cutting-plane method, in which the verifier acts as the cutting-plane oracle to refine the set of Lyapunov function candidates. We show that when the set of Lyapunov functions is full-dimensional in the parameter space, the overall procedure finds a Lyapunov function in a finite number of iterations. We demonstrate the utility of the proposed method in searching for quadratic and piecewise quadratic Lyapunov functions.


title: Collapsed Anosov flows and self orbit equivalences

abstract: We propose a generalization of the concept of discretized Anosov flows that covers a wide class of partially hyperbolic diffeomorphisms in 3-manifolds, and that we call collapsed Anosov flows. They are related with Anosov flows via a self orbit equivalence of the flow. We show that all the examples constructed in a paper by Bonatti, Gogolev, Hammerlindl and the third author belong to this class, and that it is an open and closed class among partially hyperbolic diffeomorphisms. We provide some equivalent definitions which may be more amenable to analysis and are useful in different situations. Conversely we describe the class of partially hyperbolic diffeomorphisms that are collapsed Anosov flows associated with certain types of Anosov flows.


title: Metallic Diluted Dimerization in VO2 Tweeds

abstract: The observation of electronic phase separation textures in vanadium dioxide (VO2), a prototypical electron-correlated oxide, has recently added new perspectives on the long standing debate about its metal-insulator transition and its applications. Yet, the lack of atomically-resolved information on phases accompanying such complex patterns still hinders a comprehensive understanding of the transition and its implementation in practical devices. In this work, atomic resolution imaging and spectroscopy unveils the existence of ferroelastic tweed structures on ~5 nm length scales, well below the resolution limit of currently employed spectroscopic imaging techniques. Moreover, density functional theory calculations show that such fine tweeds are formed by a metallic structure, formed by partially dimerized V-chains, that appears tetragonal (rutile) on average, but is locally monoclinic. These observations suggest that the metal/insulator coexistence among low-symmetry forms of VO2 is driven by lattice degrees of freedom, and provide a multiscale perspective for the interpretation of existing data, whereby phase coexistence and structural intermixing can occur all the way down to the atomic scale.


title: Heterogeneous Interventions Reduce the Spread of COVID-19 in Simulations on Real Mobility Data

abstract: Major interventions have been introduced worldwide to slow down the spread of the SARS-CoV-2 virus. Large-scale lockdowns of human movements are effective in reducing the spread, but they come at a cost of significantly limited societal functions. We show that natural human movements are statistically diverse, and the spread of the disease is significantly influenced by a small group of active individuals and gathering venues. We find that interventions focused on these most mobile individuals and popular venues reduce both the peak infection rate and the total infected population while retaining high social activity levels. These trends are seen consistently in simulations with real human mobility data of different scales, resolutions, and modalities from multiple cities across the world. The observation implies that compared to broad sweeping interventions, more heterogeneous strategies that are targeted based on the network effects in human mobility provide a better balance between pandemic control and regular social activities.


title: Deep Learning for estimating parameters of Gravitational Waves

abstract: In recent years, improvements in Deep Learning (DL) techniques towards Gravitational Wave (GW) astronomy have led to a significant rise in the development of various classification algorithms that have been successfully employed to extract GWs of binary blackhole merger events from noisy time-series data. However, the success of these models is constrained by the length of time-sample and the class of GW source: binary blackhole and neutron star binaries to some extent. In this work, we intended to advance the boundaries of DL techniques using Convolutional Neural Networks, to go beyond binary classification and predict the physical parameters of the events. We aim to propose an alternative method that can be employed for realtime detection and parameter prediction. The DL model we present has been trained on 12s of data to predict the GW source parameters if detected. During training, the maximum accuracy attained was 90.93%, with a validation accuracy of 89.97%.


title: Sketch-Guided Object Localization in Natural Images

abstract: We introduce the novel problem of localizing all the instances of an object (seen or unseen during training) in a natural image via sketch query. We refer to this problem as sketch-guided object localization. This problem is distinctively different from the traditional sketch-based image retrieval task where the gallery set often contains images with only one object. The sketch-guided object localization proves to be more challenging when we consider the following: (i) the sketches used as queries are abstract representations with little information on the shape and salient attributes of the object, (ii) the sketches have significant variability as they are hand-drawn by a diverse set of untrained human subjects, and (iii) there exists a domain gap between sketch queries and target natural images as these are sampled from very different data distributions. To address the problem of sketch-guided object localization, we propose a novel cross-modal attention scheme that guides the region proposal network (RPN) to generate object proposals relevant to the sketch query. These object proposals are later scored against the query to obtain final localization. Our method is effective with as little as a single sketch query. Moreover, it also generalizes well to object categories not seen during training and is effective in localizing multiple object instances present in the image. Furthermore, we extend our framework to a multi-query setting using novel feature fusion and attention fusion strategies introduced in this paper. The localization performance is evaluated on publicly available object detection benchmarks, viz. MS-COCO and PASCAL-VOC, with sketch queries obtained from `Quick, Draw!'. The proposed method significantly outperforms related baselines on both single-query and multi-query localization tasks.


title: On leaky forcing and resilience

abstract: A leak is a vertex that is not allowed to perform a force during the zero forcing process. Leaky forcing was recently introduced as a new variation of zero forcing in order to analyze how leaks in a network disrupt the zero forcing process. The $\ell$-leaky forcing number of a graph is the size of the smallest zero forcing set that can force a graph despite $\ell$ leaks. A graph $G$ is $\ell$-resilient if its zero forcing number is the same as its $\ell$-leaky forcing number. In this paper, we analyze $\ell$-leaky forcing and show that if an $(\ell-1)$-leaky forcing set $B$ is robust enough, then $B$ is an $\ell$-leaky forcing set. This provides the framework for characterizing $\ell$-leaky forcing sets. Furthermore, we consider structural implications of $\ell$-resilient graphs. We apply these results to bound the $\ell$-leaky forcing number of several graph families including trees, supertriangles, and grid graphs. In particular, we resolve a question posed by Dillman and Kenter concerning the upper bound on the $1$-leaky forcing number of grid graphs.


title: P-Wave Two-Particle Bound and Scattering States in a Finite Volume including QED

abstract: The mass shifts for two-fermion bound and scattering P-wave states subject to the long-range interactions due to QED in the non-relativistic regime are derived. Introducing a short range force coupling the spinless fermions to one unit of angular momentum in the framework of pionless EFT, we first calculate both perturbatively and non-perturbatively the Coulomb corrections to fermion-fermion scattering in the continuum and infinite volume context. Motivated by the research on particle-antiparticle bound states, we extend the results to fermions of identical mass and opposite charge. Second, we transpose the system onto a cubic lattice with periodic boundary conditions and we calculate the finite volume corrections to the energy of the lowest bound and unbound $T_1^{\pm}$ eigenstates. In particular, power law corrections proportional to the fine structure constant and resembling the recent results for S-wave states are found. Higher order contributions in $\alpha$ are neglected, since the gapped nature of the momentum operator in the lattice environnement allows for a perturbative treatment of the QED interactions.


title: On the Hardness of Massively Parallel Computation

abstract: We investigate whether there are inherent limits of parallelization in the (randomized) massively parallel computation (MPC) model by comparing it with the (sequential) RAM model. As our main result, we show the existence of hard functions that are essentially not parallelizable in the MPC model. Based on the widely-used random oracle methodology in cryptography with a cryptographic hash function $h:\{0,1\}^n \rightarrow \{0,1\}^n$ computable in time $t_h$, we show that there exists a function that can be computed in time $O(T\cdot t_h)$ and space $S$ by a RAM algorithm, but any MPC algorithm with local memory size $s < S/c$ for some $c>1$ requires at least $\tilde{\Omega}(T)$ rounds to compute the function, even in the average case, for a wide range of parameters $n \leq S \leq T \leq 2^{n^{1/4}}$. Our result is almost optimal in the sense that by taking $T$ to be much larger than $t_h$, \textit{e.g.}, $T$ to be sub-exponential in $t_h$, to compute the function, the round complexity of any MPC algorithm with small local memory size is asymptotically the same (up to a polylogarithmic factor) as the time complexity of the RAM algorithm. Our result is obtained by adapting the so-called compression argument from the data structure lower bounds and cryptography literature to the context of massively parallel computation.


title: A New Perspective on Pool-Based Active Classification and False-Discovery Control

abstract: In many scientific settings there is a need for adaptive experimental design to guide the process of identifying regions of the search space that contain as many true positives as possible subject to a low rate of false discoveries (i.e. false alarms). Such regions of the search space could differ drastically from a predicted set that minimizes 0/1 error and accurate identification could require very different sampling strategies. Like active learning for binary classification, this experimental design cannot be optimally chosen a priori, but rather the data must be taken sequentially and adaptively. However, unlike classification with 0/1 error, collecting data adaptively to find a set with high true positive rate and low false discovery rate (FDR) is not as well understood. In this paper we provide the first provably sample efficient adaptive algorithm for this problem. Along the way we highlight connections between classification, combinatorial bandits, and FDR control making contributions to each.


title: Data-Informed Decomposition for Localized Uncertainty Quantification of Dynamical Systems

abstract: Industrial dynamical systems often exhibit multi-scale response due to material heterogeneities, operation conditions and complex environmental loadings. In such problems, it is the case that the smallest length-scale of the systems dynamics controls the numerical resolution required to effectively resolve the embedded physics. In practice however, high numerical resolutions is only required in a confined region of the system where fast dynamics or localized material variability are exhibited, whereas a coarser discretization can be sufficient in the rest majority of the system. To this end, a unified computational scheme with uniform spatio-temporal resolutions for uncertainty quantification can be very computationally demanding. Partitioning the complex dynamical system into smaller easier-to-solve problems based of the localized dynamics and material variability can reduce the overall computational cost. However, identifying the region of interest for high-resolution and intensive uncertainty quantification can be a problem dependent. The region of interest can be specified based on the localization features of the solution, user interest, and correlation length of the random material properties. For problems where a region of interest is not evident, Bayesian inference can provide a feasible solution. In this work, we employ a Bayesian framework to update our prior knowledge on the localized region of interest using measurements and system response. To address the computational cost of the Bayesian inference, we construct a Gaussian process surrogate for the forward model. Once, the localized region of interest is identified, we use polynomial chaos expansion to propagate the localization uncertainty. We demonstrate our framework through numerical experiments on a three-dimensional elastodynamic problem.


title: On the globalization of Riemannian Newton method

abstract: In the present paper, in order to fnd a singularity of a vector field defined on Riemannian manifolds, we present a new globalization strategy of Newton method and establish its global convergence with superlinear rate. In particular, this globalization generalizes for a general retraction the existing damped Newton's method. The presented global convergence analysis does not require any hypotesesis on singularity of the vector field. We applied the proposed method to solve the truncated singular value problem on the product of two Stiefel manifolds, the dextrous hand grasping problem on the cone of symmetric positive definite matrices and the Rayleigh quotient on the sphere. Moreover, some academic problems are solved. Numerical experiments are presented showing that the proposed algorithm has better robustness compared with the aforementioned method.


title: Donkin-Koppinen filtration for GL(m|n) and generalized Schur superalgebras

abstract: The paper contains results that characterize the Donkin-Koppinen filtration of the coordinate superalgebra $K[G]$ of the general linear supergroup $G=GL(m|n)$ by its subsupermodules $C_{\Gamma}=O_{\Gamma}(K[G])$. Here, the supermodule $C_{\Gamma}$ is the largest subsupermodule of $K[G]$ whose composition factors are irreducible supermodules of highest weight $\lambda$, where $\lambda$ belongs to a finitely-generated ideal $\Gamma$ of the poset $X(T)^+$ of dominant weights of $G$. A decomposition of $G$ as a product of subsuperschemes $U^-\times G_{ev}\times U^+$ induces a superalgebra isomorphism $\phi^* : K[U^-]\otimes K[G_{ev}]\otimes K[U^+]\simeq K[G]$. We show that $C_{\Gamma}=\phi^*(K[U^-]\otimes M_{\Gamma}\otimes K[U^+])$, where $M_{\Gamma}=O_{\Gamma}(K[G_{ev}])$. Using the basis of the module $M_{\Gamma}$, given by generalized bideterminants, we describe a basis of $C_{\Gamma}$. Since each $C_{\Gamma}$ is a subsupercoalgebra of $K[G]$, its dual $C_{\Gamma}^*=S_{\Gamma}$ is a (pseudocompact) superalgebra, called the generalized Schur superalgebra. There is a natural superalgebra morphism $\pi_{\Gamma}:Dist(G)\to S_{\Gamma}$ such that the image of the distribution algebra $Dist(G)$ is dense in $S_{\Gamma}$. For the ideal $X(T)^+_{l}$, of all weights of fixed length $l$, the generators of the kernel of $\pi_{X(T)^+_{l}}$ are described.


title: Performance characterization of a novel deep learning-based MR image reconstruction pipeline

abstract: A novel deep learning-based magnetic resonance imaging reconstruction pipeline was designed to address fundamental image quality limitations of conventional reconstruction to provide high-resolution, low-noise MR images. This pipeline's unique aims were to convert truncation artifact into improved image sharpness while jointly denoising images to improve image quality. This new approach, now commercially available at AIR Recon DL (GE Healthcare, Waukesha, WI), includes a deep convolutional neural network (CNN) to aid in the reconstruction of raw data, ultimately producing clean, sharp images. Here we describe key features of this pipeline and its CNN, characterize its performance in digital reference objects, phantoms, and in-vivo, and present sample images and protocol optimization strategies that leverage image quality improvement for reduced scan time. This new deep learning-based reconstruction pipeline represents a powerful new tool to increase the diagnostic and operational performance of an MRI scanner.


title: Dust content solutions for the Alcubierre warp drive spacetime

abstract: The Alcubierre metric is a spacetime geometry where a massive particle inside a spacetime distortion, called warp bubble, is able to travel at velocities arbitrarily higher than the velocity of light, a feature known as the warp drive. This is a consequence of general relativity, which allows global superluminal velocities but restricts local speeds to subluminal ones as required by special relativity. In this work we solved the Einstein equations for the Alcubierre warp drive spacetime geometry considering the dust matter distribution as source, since the Alcubierre metric was not originally advanced as a solution of the Einstein equations, but as a spacetime geometry proposed without a source gravity field. We found out that all Einstein equations solutions of this geometry containing pressureless dust lead to vacuum solutions. We also concluded that these solutions connect the Alcubierre metric to the Burgers equation, which describes shock waves moving through an inviscid fluid. Our results also indicated that these shock waves behave as plane waves.


title: Rapid charge redistribution leading to core hollowing in a high-intensity ion beam

abstract: Recently, the first direct measurement of a full 6D accelerator beam distribution was reported [1]. That work observed a correlation between energy and transverse coordinates, for which the energy distribution becomes hollowed and double-peaked near the transverse core. In this article, a similar structure is shown to emerge in expansion of an initially uncorrelated, high density bunched beam as the result of velocity perturbation from nonlinear space charge forces. This hollowing is obscured when the 6D phase space is projected onto one- and two-dimensional axes. This phenomenon has not been widely recognized in accelerator systems, but parallels can be drawn to observations of laser-ionized nanoclusters and electron sources for diffraction. While this effect provides insight into the origin of the measured core correlation, it does not provide a complete description. A better reproduction of the measured structure can be obtained via self-consistent simulation through the radio-frequency quadrupole. [1] B. Cathey, S. Cousineau, A. Aleksandrov, and A. Zhukov, PRL 121, 064804 (2018).


title: Ab Initio Molecular Dynamics on Quantum Computers

abstract: Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, a quantum computer-based AIMD method is presented. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. We compute the energy gradients numerically using the Hellmann-Feynman theorem, finite differences, and a correlated sampling technique. Our method only requires additional classical calculations of electron integrals for each degree of freedom, without any additional computations on a quantum computer beyond the initial VQE run. To achieve comparable accuracy, our gradient calculation method requires three to five orders of magnitude fewer measurements than other brute force methods without correlated sampling. As a proof of concept, AIMD dynamics simulations are demonstrated for the H2 molecule on IBM quantum devices. To the best of our knowledge, it is the first successful attempt to run AIMD on quantum devices for a chemical system. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction (FCI) wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular and material systems.


title: Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds from ASAS-SN: Implications for Thorne-\.Zytkow Objects and Super-Asymptotic Giant Branch Stars

abstract: Stars with unusual properties can provide a wealth of information about rare stages of stellar evolution and exotic physics. However, determining the true nature of peculiar stars is often difficult. In this work, we conduct a systematic search for cool and luminous stars in the Magellanic Clouds with extreme variability, motivated by the properties of the unusual SMC star and Thorne-\.Zytkow Object (T\.ZO) candidate HV2112. Using light curves from ASAS-SN we identify 38 stars with surface temperatures T $<$ 4800K, luminosities $\log$(L/L$_\odot$) $>$ 4.3, variability periods $>$ 400 days, and variability amplitudes $\Delta$V $>$ 2.5 mag. Eleven of these stars possess the distinctive double-peaked light curve morphology of HV2112. We use the pulsation properties and derived occurrence rates for these 12 objects to constrain their nature. From comparisons to stellar populations and models, we find that one star may be a red supergiant with large amplitude pulsations. For the other 11 stars we derive current masses of $\sim$5-10 M$_{\odot}$, below the theoretical minimum mass of $\sim$15 M$_{\odot}$ for T\.ZOs to be stable, casting doubt on this interpretation. Instead, we find that the temperatures, luminosities, mass-loss rates, and periods of these stars are consistent with predictions for super-Asymptotic Giant Branch (s-AGB) stars that have begun carbon burning but have not reached the superwind phase. We infer lifetimes in this phase of $\sim($1$-$7) $\times$ 10$^{4}$ years, also consistent with an s-AGB interpretation. If confirmed, these objects would represent the first identified population of s-AGB stars, illuminating the transition between low- and high-mass stellar evolution.


title: Optimal selection of the number of control units in kNN algorithm to estimate average treatment effects

abstract: We propose a simple approach to optimally select the number of control units in k nearest neighbors (kNN) algorithm focusing in minimizing the mean squared error for the average treatment effects. Our approach is non-parametric where confidence intervals for the treatment effects were calculated using asymptotic results with bias correction. Simulation exercises show that our approach gets relative small mean squared errors, and a balance between confidence intervals length and type I error. We analyzed the average treatment effects on treated (ATET) of participation in 401(k) plans on accumulated net financial assets confirming significant effects on amount and positive probability of net asset. Our optimal k selection produces significant narrower ATET confidence intervals compared with common practice of using k=1.


title: High dimensional characterization of the longitudinal phase space formed in a radio frequency quadrupole

abstract: Modern accelerator front ends almost exclusively include radio-frequency quadrupoles for initial capture and focusing of the low-energy beam. Dynamics in the RFQ define the longitudinal bunch parameters. Simulation of the SNS RFQ with PARMTEQ seeded with a realistic LEBT distribution produces a 2.5 MeV, 40 mA H- beam with root-mean-square emittance of 130 deg-keV. In measurement, a detailed characterization of the longitudinal phase space is made, including a novel study of the dependence of longitudinal emittance on transverse coordinates. This work introduces a new virtual slit technique that provides sub-slit resolution in an energy spectrometer as well as an approach for visualizing 4D phase space data. Through simulation and measurement, the RFQ-formed bunch is confirmed to have significant internal correlated structure. The high-dimensional features are shown to be in qualitative agreement. However, the measured rms emittances are up to 30% lower than predicted, closer to the design value of 95 deg-keV.


title: Summary of Workshop on Common Neutrino Event Generator Tools

abstract: A neutrino community workshop was held at Fermilab in Jan 2020, with the aim of developing an implementation plan for a set of common interfaces to Neutrino Event Generators. This white paper summarizes discussions at the workshop and the resulting plan.


title: On the fully nonlinear Alt-Phillips equation

abstract: For a parameter $\gamma\in(1,2)$, we study the fully nonlinear version of the Alt-Phillips equation, $F(D^2u)=u^{\gamma-1}$, for $u\ge 0.$ We establish the optimal regularity of the solution, as well as the $C^1$ regularity of the regular part of the free boundary.


title: An Evaluation Framework of End-to-End 5G Millimeter Wave Communication for Connected Vehicle Applications

abstract: The internet-of-things (IoT) environment connects different intelligent devices together and enables seamless data communication between the connected devices. Connected vehicles (CVs) are one of the primary example of the IoT, and the efficient, reliable, and safe operation of CVs demands a reliable wireless communication system, which can ensure high throughput and low communication latency. The 5G millimeter wave (5G mmWave) wireless communication network offers such benefits, which can be the enabler of CV applications, especially for dense urban areas with high number of CVs. In this study, we present a simulation-based evaluation framework of end-to-end 5G mmWave communication for CV applications. In addition, we compare the 5G mmWave with the Dedicated Short Range Communication (DSRC) technology for a CV application. The simulation framework is developed using two simulators, a network simulator and a traffic simulator. In order to develop the framework in this study, we have used Network Simulator 3 (ns-3) and SUMO, an open-source microscopic roadway traffic simulator. We have used end-to-end latency, packet loss and throughput as the performance evaluation metrics. We have found that for dense urban areas, 5G mmWave can achieve higher throughput, lower latency and lower data loss compared to DSRC. 5G mmWave can support CV applications with high throughput requirement on the downlink data flow. Through further investigation, we have found that the performance of 5G mmWave is significantly impacted by the penetration level of CVs, maximum CV speed, and CV application requirements.


title: New blow-up result for the weakly coupled wave equations with a scale-invariant damping and time derivative nonlinearity

abstract: We consider in this article the weakly coupled system of wave equations in the \textit{scale-invariant case} and with time-derivative nonlinearities. Under the usual assumption of small initial data, we obtain an improvement of the delimitation of the blow-up region by obtaining a new candidate for the critical curve. More precisely, we enhance the results obtained in \cite{Palmieri} for the system under consideration in the present work. We believe that our result is optimal in the sense that beyond the blow-up region obtained here we may conjecture the global existence of the solution.


title: Fast Dimension Independent Private AdaGrad on Publicly Estimated Subspaces

abstract: We revisit the problem of empirical risk minimziation (ERM) with differential privacy. We show that noisy AdaGrad, given appropriate knowledge and conditions on the subspace from which gradients can be drawn, achieves a regret comparable to traditional AdaGrad plus a well-controlled term due to noise. We show a convergence rate of $O(\text{Tr}(G_T)/T)$, where $G_T$ captures the geometry of the gradient subspace. Since $\text{Tr}(G_T)=O(\sqrt{T})$ we can obtain faster rates for convex and Lipschitz functions, compared to the $O(1/\sqrt{T})$ rate achieved by known versions of noisy (stochastic) gradient descent with comparable noise variance. In particular, we show that if the gradients lie in a known constant rank subspace, and assuming algorithmic access to an envelope which bounds decaying sensitivity, one can achieve faster convergence to an excess empirical risk of $\tilde O(1/\epsilon n)$, where $\epsilon$ is the privacy budget and $n$ the number of samples. Letting $p$ be the problem dimension, this result implies that, by running noisy Adagrad, we can bypass the DP-SGD bound $\tilde O(\sqrt{p}/\epsilon n)$ in $T=(\epsilon n)^{2/(1+2\alpha)}$ iterations, where $\alpha \geq 0$ is a parameter controlling gradient norm decay, instead of the rate achieved by SGD of $T=\epsilon^2n^2$. Our results operate with general convex functions in both constrained and unconstrained minimization. Along the way, we do a perturbation analysis of noisy AdaGrad of independent interest. Our utility guarantee for the private ERM problem follows as a corollary to the regret guarantee of noisy AdaGrad.


title: Toward an End-to-End Auto-tuning Framework in HPC PowerStack

abstract: Efficiently utilizing procured power and optimizing performance of scientific applications under power and energy constraints are challenging. The HPC PowerStack defines a software stack to manage power and energy of high-performance computing systems and standardizes the interfaces between different components of the stack. This survey paper presents the findings of a working group focused on the end-to-end tuning of the PowerStack. First, we provide a background on the PowerStack layer-specific tuning efforts in terms of their high-level objectives, the constraints and optimization goals, layer-specific telemetry, and control parameters, and we list the existing software solutions that address those challenges. Second, we propose the PowerStack end-to-end auto-tuning framework, identify the opportunities in co-tuning different layers in the PowerStack, and present specific use cases and solutions. Third, we discuss the research opportunities and challenges for collective auto-tuning of two or more management layers (or domains) in the PowerStack. This paper takes the first steps in identifying and aggregating the important R&D challenges in streamlining the optimization efforts across the layers of the PowerStack.


title: Modeling of relative permeabilities including dynamic wettability transition zones

abstract: Wettability is a pore-scale property that impacts the relative movement and distribution of fluids in a porous medium. There are reservoir fluids that provoke the surface within pores to undergo a wettability change. This wettability change, in turn, alters the dynamics of relative permeabilities at the Darcy scale. Thus, modeling the impact of wettability change in relative permeabilities is essential to understand fluids interaction in porous media. In this study, we include time-dependent wettability change into the relative permeability--saturation relation by modifying the existing relative permeability function. To do so, we assume the wettability change is represented by the sorption-based model that is exposure time and chemistry dependent. This pore-scale model is then coupled with a triangular bundle-of-tubes model to simulate exposure time-dependent relative permeabilities data. The simulated data is used to characterize and quantify the wettability dynamics in the relative permeability--saturation curves. This study further shows the importance of accurate prediction of the relative permeability in a dynamically altering porous medium.


title: Interaction of a wave packet with potential structures moving with acceleration

abstract: The paper is devoted to a numerical study of the problem of interaction of the wave packet with potential structures moving with constant acceleration. In all the cases considered the result of the interaction is a change in the velocity spectrum. In the first approximation the magnitude of the shift in the spectrum is determined by the product of acceleration by group delay time. Also, as the direction of acceleration reverses the effect changes its sign. The results are completely consistent with the idea of universality of the Effect of Acceleration which consists in a change in the frequency of the wave at scattering on an object moving with acceleration.


title: Rotating Black Holes in Three-Dimensional Horava Gravity Revisited

abstract: I revisit rotating black hole solutions in three-dimensional Horava gravity with z = 2 as a simpler set-up of the renormalizable quantum gravity `a la Lifshitz and DeWitt. The solutions have a curvature singularity at the origin for a non-vanishing rotation parameter J, unlike the black holes in three-dimensional Einstein gravity. For anti-de Sitter space, there are black hole event horizons as usual and the singularity is not naked, in agreement with the cosmic censorship. On the other hand, for flat or de Sitter space, the earlier solution has also a cosmic-censorship problem because there are no conventional black hole horizons as in Einstein gravity, other than the usual cosmological horizon for the latter case, so that the singularity could be naked in Horava gravity. However, with the help of recent corrections, I show that the solutions have a peculiar black hole horizon at the origin so that the singularity is not naked even without the conventional black hole horizons in flat or de Sitter case, due to the Lorentz-violating higher-derivative terms. On the other hand, I note also that a new ``cosmological" horizon exists even for the flat case, contrary to the usual wisdom, due to combined effects of the higher derivatives and the angular-momentum barrier. I study an unified treatment of their unusual black hole thermodynamics for the flat and de Sitter spaces, as well as the anti-de Sitter space, which might be due to lack of the absolute horizons in the Lorentz-violating gravity.


title: CHESS: an Innovative Concept for High-Resolution, Far-UV Spectroscopy

abstract: The space ultraviolet (UV) is a critical astronomical observing window, where a multitude of atomic, ionic, and molecular signatures provide crucial insight into planetary, interstellar, stellar, intergalactic, and extragalactic objects. The next generation of large space telescopes require highly sensitive, moderate-to-high resolution UV spectrograph. However, sensitive observations in the UV are difficult, as UV optical performance and imaging efficiencies have lagged behind counterparts in the visible and infrared regimes. This has historically resulted in simple, low-bounce instruments to increase sensitivity. In this study, we present the design, fabrication, and calibration of a simple, high resolution, high throughput far-UV spectrograph - the Colorado High-resolution Echelle Stellar Spectrograph (CHESS). CHESS is a sounding rocket payload to demonstrate the instrument design for the next-generation UV space telescopes. We present tests and results on the performance of several state-of-the-art diffraction grating and detector technologies for far-UV astronomical applications that were flown aboard the first two iterations of CHESS. The CHESS spectrograph was used to study the atomic-to-molecular transitions within translucent cloud regions in the interstellar medium (ISM) through absorption spectroscopy. The first two flights looked at the sightlines towards alpha Virgo and epsilon Persei, and flight results are presented.


title: Single Board Computers (SBC): The Future of Next Generation Pedagogies in Pakistan

abstract: ARM processors have taken over the mobile industry from a long time now. Future of data centers and the IT industry is estimated to make use of ARM Processors. Projects like Openstack on ARM are enabling use of ARM in data centers . Single board computers (SBCs) based on ARM processors have become the norm these days. Reason for their popularity lies in their cost effective and power efficient nature. There are hundreds of them available in the market having different sizes, compute power and prices. The reason for their popularity is largely due to the rise of new technology called IoT (Internet of Things) but there is another perspective where they can become handy. Low Price and Power Usage of single board computers makes them top candidate to be used for teaching many courses with hands-on experience in developing countries like Pakistan. Many boards support full Linux distributions and can be used as general-purpose computers while many of them are open hardware based. In this paper, we have reviewed the famous options available and tried to figure out which of them are better for teaching what kind of courses.


title: Cycles of a given length in tournaments

abstract: We study the asymptotic behavior of the maximum number of directed cycles of a given length in a tournament: let $c(\ell)$ be the limit of the ratio of the maximum number of cycles of length $\ell$ in an $n$-vertex tournament and the expected number of cycles of length $\ell$ in the random $n$-vertex tournament, when $n$ tends to infinity. It is well-known that $c(3)=1$ and $c(4)=4/3$. We show that $c(\ell)=1$ if and only if $\ell$ is not divisible by four, which settles a conjecture of Bartley and Day. If $\ell$ is divisible by four, we show that $1+2\cdot\left(2/\pi\right)^{\ell}\le c(\ell)\le 1+\left(2/\pi+o(1)\right)^{\ell}$ and determine the value $c(\ell)$ exactly for $\ell = 8$. We also give a full description of the asymptotic structure of tournaments with the maximum number of cycles of length $\ell$ when $\ell$ is not divisible by four or $\ell\in\{4,8\}$.


title: Vortex dynamics and transport phenomena in stenotic aortic models using Echo-PIV

abstract: In this work, we propose a novel approach which combines ultrasound with Eulerian and Lagrangian descriptors, to analyse blood flow dynamics and fluid transport in stenotic aortic models with morphology, mechanical and optical properties close to those of real arteries. To this end, vorticity, particle residence time (PRT), particle's final position (FP) and finite time Lyapunov's exponents (FTLE) were computed from the experimental fluid velocity fields acquired using ultrasonic particle imaging velocimetry (Echo-PIV). For the experiments, CT-images were used to create morphological realistic models of the descending aorta with 0%, 35% and 50% occlusion degree with same mechanical properties as real arteries. Each model was connected to a circuit with a pulsatile programmable pump which mimics physiological flow and pressure conditions. The pulsatile frequency was set to 0.9 Hz (55 bpm) and the upstream peak Reynolds number (Re) was changed from 1100 to 2000. Flow in the post-stenotic region was composed of two main structures: a high velocity jet over the stenosis throat and a recirculation region behind the stenosis where vortex form and shed. We characterized vortex kinematics showing that vortex propagation velocity increases with Re. Moreover, from the FTLE field we identified Lagrangian Coherent Structures (i.e. material barriers) that dictate transport behind the stenosis. The size and strength of those barriers increased with Re and the occlusion degree. Finally, from the PRT and FP, we showed that independently of Re, the same amount of fluid remains on the stenosis over more than a pulsatile period, which combined with large FTLE values may provide an alternative way to understand stenosis growth.


title: Degree for weakly upper semicontinuous perturbations of quasi-$m$-accretive operators

abstract: In the paper we provide the construction of a coincidence degree being a homotopy invariant detecting the existence of solutions of equations or inclusions of the form $Ax\in F(x)$, $x\in U$, where $A\colon D(A)\multimap E$ is an $m$-accretive operator in a Banach space $ E$, $F\colon K\multimap E$ is a weakly upper semicontinuous set-valued map constrained to an open subset $U$ of a closed set $K\subset E$. Two different approaches will be presented. The theory is applied to show the existence of nontrivial positive solutions of some nonlinear second order partial differential equations with discontinuities.


title: Adaptation Algorithms for Neural Network-Based Speech Recognition: An Overview

abstract: We present a structured overview of adaptation algorithms for neural network-based speech recognition, considering both hybrid hidden Markov model / neural network systems and end-to-end neural network systems, with a focus on speaker adaptation, domain adaptation, and accent adaptation. The overview characterizes adaptation algorithms as based on embeddings, model parameter adaptation, or data augmentation. We present a meta-analysis of the performance of speech recognition adaptation algorithms, based on relative error rate reductions as reported in the literature.


title: Audio-Visual Event Localization via Recursive Fusion by Joint Co-Attention

abstract: The major challenge in audio-visual event localization task lies in how to fuse information from multiple modalities effectively. Recent works have shown that attention mechanism is beneficial to the fusion process. In this paper, we propose a novel joint attention mechanism with multimodal fusion methods for audio-visual event localization. Particularly, we present a concise yet valid architecture that effectively learns representations from multiple modalities in a joint manner. Initially, visual features are combined with auditory features and then turned into joint representations. Next, we make use of the joint representations to attend to visual features and auditory features, respectively. With the help of this joint co-attention, new visual and auditory features are produced, and thus both features can enjoy the mutually improved benefits from each other. It is worth noting that the joint co-attention unit is recursive meaning that it can be performed multiple times for obtaining better joint representations progressively. Extensive experiments on the public AVE dataset have shown that the proposed method achieves significantly better results than the state-of-the-art methods.


title: Stable Polar Oxynitrides through Epitaxial Strain

abstract: We investigate energetically favorable structures of ABO$_2$N oxynitrides as functions of pressure and strain via swarm-intelligence-based structure prediction methods, DFT lattice dynamics and first-principles molecular dynamics. We predict several thermodynamically stable polar oxynitride perovskites under high pressures. In addition, we find that ferroelectric polar phases of perovskite-structured oxynitrides can be thermodynamically stable and synthesized at high pressure on appropriate substrates. The dynamical stability of the ferroelectric oxynitrides under epitaxial strain at ambient pressure also imply the possibility to synthesize them using pulsed laser deposition or other atomic layer deposition methods. Our results have broad implications for further exploration of other oxynitride materials as well. We performed first-principles molecular dynamics and find that the polar perovskite of YSiO$_2$N is metastable up to at least 600 K under compressive epitaxial strain before converting to the stable wollastonite-like structures. YSiO$_2$N is stabilized under pressure with extensional epitaxial strain. We predict that LaSi$_2$N, LaGeO$_2$N, BiSiO$_2$N, and BiGeO$_2$N are metastable as ferroelectric perovskites at zero pressure even without epitaxial strain.


title: Minimal Surfaces in $\mathbb{H}_2\times \mathbb{R}$: Nonfillable Curves

abstract: We study the asymptotic Plateau problem in $\mathbb{H}_2\times \mathbb{R}$. We give the first examples of non-fillable finite curves with no thin tail in the asymptotic cylinder. Furthermore, we study the fillability question for infinite curves in the asymptotic boundary.


title: Convergence of exclusion processes and KPZ equation to the KPZ fixed point

abstract: We show that under the 1:2:3 scaling, critically probing large space and time, the height function of finite range asymmetric exclusion processes and the KPZ equation converge to the KPZ fixed point, constructed earlier as a limit of the totally asymmetric simple exclusion process through exact formulas. Consequently, based on recent results of \cite{wu},\cite{DM20}, the KPZ line ensemble converges to the Airy line ensemble. For the KPZ equation we are able to start from a continuous function plus a finite collection of narrow wedges. For nearest neighbour exclusions, we can take (discretizations) of height functions with $h(x)\le C(1+|x|)$. For non-nearest neighbour exclusions, we are restricted at the present time to a class of (random) initial data, dense in continuous functions in the topology of uniform convergence on compacts. The method is by comparison of the transition probabilities of finite range exclusion processes and the totally asymmetric simple exclusion processes using energy estimates.


title: COVID-Robot: Monitoring Social Distancing Constraints in Crowded Scenarios

abstract: Maintaining social distancing norms between humans has become an indispensable precaution to slow down the transmission of COVID-19. We present a novel method to automatically detect pairs of humans in a crowded scenario who are not adhering to the social distance constraint, i.e. about 6 feet of space between them. Our approach makes no assumption about the crowd density or pedestrian walking directions. We use a mobile robot with commodity sensors, namely an RGB-D camera and a 2-D lidar to perform collision-free navigation in a crowd and estimate the distance between all detected individuals in the camera's field of view. In addition, we also equip the robot with a thermal camera that wirelessly transmits thermal images to a security/healthcare personnel who monitors if any individual exhibits a higher than normal temperature. In indoor scenarios, our mobile robot can also be combined with static mounted CCTV cameras to further improve the performance in terms of number of social distancing breaches detected, accurately pursuing walking pedestrians etc. We highlight the performance benefits of our approach in different static and dynamic indoor scenarios.


title: An Approximate Maximum Likelihood Time Synchronization Algorithm for Zero-padded OFDM in Channels with Impulsive Gaussian Noise

abstract: In wireless communication systems, Orthogonal Frequency-Division Multiplexing (OFDM) includes variants using either a cyclic prefix (CP) or a zero padding (ZP) as the guard interval to avoid inter-symbol interference. OFDM is ideally suited to deal with frequency-selective channels and additive white Gaussian noise (AWGN); however, its performance may be dramatically degraded in the presence of impulse noise. While the ZP variants of OFDM exhibit lower bit error rate (BER)and higher energy efficiency compared to their CP counterparts,they demand strict time synchronization, which is challenging in the absence of pilot and CP. Moreover, on the contrary to AWGN, impulse noise severely corrupts data. In this paper, a new low-complexity timing offset (TO) estimator for ZP-OFDM for practical impulsive-noise environments is proposed, where relies on the second-other statistics of the multipath fading channel and noise. Performance comparison with existing TO estimators demonstrates either a superior performance in terms of lock-in probability or a significantly lower complexity over a wide range of signal-to-noise ratio (SNR) for various practical scenarios.


title: An Object-Oriented Framework for Designing Reusable and Maintainable DEVS Models using Design Patterns

abstract: Design patterns are well practices to share software development experiences. These patterns allow enhancing reusability, readability and maintainability of architecture and code of software applications. As simulation applies computerized models to produce traces in order to obtain results and conclusions, designers of simulation explored design patterns to make the simulation code more reusable, more readable and easy to maintain, in addition to design complex software oriented simulation modeling. In DEVS (Discrete Event System specification), the designers have successfully designed simulations, frameworks, tools, etc. However, some issues remain still open and should be explored like how a piece of code that implements a set of states, events and transitions may be reused to design a new DEVS model? How may a DEVS model be extended to a new formalism? Etc. In this paper, we address these issues and we propose a set of patterns that may serve as guidelines to designers of DEVS models and its extensions and may contribute to the design of an operational simulation framework. These patterns are inspired partly by the available designs of DEVS community and software engineering developers.


title: Ab initio neutrinoless double-beta decay matrix elements for 48Ca, 76Ge, and 82Se

abstract: We calculate basis-space converged neutrinoless $\beta \beta$ decay nuclear matrix elements for the lightest candidates: 48Ca, 76Ge and 82Se. Starting from initial two- and three-nucleon forces, we apply the ab initio in-medium similarity renormalization group to construct valence-space Hamiltonians and consistently transformed $\beta \beta$-decay operators. We find that the tensor component is non-negligible in 76Ge and 82Se, and resulting nuclear matrix elements are overall 25-45% smaller than those obtained from the phenomenological shell model. While a final matrix element with uncertainties still requires substantial developments, this work nevertheless opens a path toward a true first-principles calculation of neutrinoless $\beta \beta$ decay in all nuclei relevant for ongoing large-scale searches.


title: Technical Report: Property-Directed Verified Monitoring of Signal Temporal Logic

abstract: Signal Temporal Logic monitoring over numerical simulation traces has emerged as an effective approach to approximate verification of continuous and hybrid systems. In this report we explore an exact verification procedure for STL properties based on monitoring verified traces in the form of Taylor model flowpipes as produced by the Flow* verified integrator. We explore how tight integration with Flow*'s symbolic flowpipe representation can lead to more precise and more efficient monitoring. We then show how the performance of monitoring can be increased substantially by introducing masks, a property-directed refinement of our method which restricts flowpipe monitoring to the time regions relevant to the overall truth of a complex proposition. Finally, we apply our implementation of these methods to verifying properties of a challenging continuous system, evaluating the impact of each aspect of our procedure on monitoring performance.


title: Periodic Solutions to Reversible Second Order Autonomous DDEs in Prescribed Symmetric Nonconvex Domains

abstract: The existence and spatio-temporal patterns of $2\pi$-periodic solutions to second order reversible equivariant autonomous systems with commensurate delays are studied using the Brouwer $O(2) \times \Gamma \times \mathbb Z_2$-equivariant degree theory. The solutions are supposed to take their values in a prescribed symmetric domain $D$, while $O(2)$ is related to the reversal symmetry combined with the autonomous form of the system. The group $\Gamma$ reflects symmetries of $D$ and/or possible coupling in the corresponding network of identical oscillators, and $\mathbb Z_2$ is related to the oddness of the right-hand side. Abstract results, based on the use of Gauss curvature of $\partial D$, Hartman-Nagumo type {\it a priori bounds} and Brouwer equivariant degree techniques, are supported by a concrete example with $\Gamma = D_8$ -- the dihedral group of order $16$.


title: New Techniques for Proving Fine-Grained Average-Case Hardness

abstract: The recent emergence of fine-grained cryptography strongly motivates developing an average-case analogue of Fine-Grained Complexity (FGC). This paper defines new versions of OV, $k$SUM and zero-$k$-clique that are both worst-case and average-case fine-grained hard assuming the core hypotheses of FGC. We then use these as a basis for fine-grained hardness and average-case hardness of other problems. The new problems represent their inputs in a certain ``factored'' form. We call them ``factored''-OV, ``factored''-zero-$k$-clique and ``factored''-$3$SUM. We show that factored-$k$-OV and factored $k$SUM are equivalent and are complete for a class of problems defined over Boolean functions. Factored zero-$k$-clique is also complete, for a different class of problems. Our hard factored problems are also simple enough that we can reduce them to many other problems, e.g.~to edit distance, $k$-LCS and versions of Max-Flow. We further consider counting variants of the factored problems and give WCtoACFG reductions for them for a natural distribution. Through FGC reductions we then get average-case hardness for well-studied problems like regular expression matching from standard worst-case FGC assumptions. To obtain our WCtoACFG reductions, we formalize the framework of [Boix-Adsera et al. 2019] that was used to give a WCtoACFG reduction for counting $k$-cliques. We define an explicit property of problems such that if a problem has that property one can use the framework on the problem to get a WCtoACFG self reduction. We then use the framework to slightly extend Boix-Adsera et al.'s average-case counting $k$-cliques result to average-case hardness for counting arbitrary subgraph patterns of constant size in $k$-partite graphs...


title: The rising of collective forgetting and cultural selectivity in inventors and physicists communities

abstract: How long until this paper is forgotten? Collective forgetting is the process by which the attention received by cultural pieces decays as time passes. Recent work modeled this decay as the result of two different processes, one linked to communicative memory --memories sustained by human communication-- and cultural memory --memories sustained by the physical recording of content. Yet, little is known on how the collective forgetting dynamic changes over time. Are older cultural pieces forgotten at a lower rate than newer ones? Here, we study the temporal changes of collective memory and attention by focusing on two knowledge communities: inventors and physicists. We use data on patents from the United States Patent and Trademark Office (USPTO) and physics papers published in the American Physical Society (APS) to quantify how collective forgetting has changed over time. The model enables us to distinguish between two branches of forgetting. One branch is short-lived, going directly from communicative memory to oblivion. The other one is long-lived going from communicative to cultural memory and then to oblivion. The data analysis shows an increasing forgetting rate for both communities as the information grows. Furthermore, these knowledge communities seem to be increasing their selectivity at storing valuable cultural pieces in their cultural memory. These findings provide empirical confirmation on the forgetting as an annulment hypothesis and show that knowledge communities can effectively slow down the rising of collective forgetting at improving their cultural selectivity.


title: Isometric flows of $G_2$-structures

abstract: We survey recent progress in the study of flows of isometric $G_2$-structures on 7-dimensional manifolds, that is, flows that preserve the metric, while modifying the $G_2$-structure. In particular, heat flows of isometric $G_2$-structures have been recently studied from several different perspectives, in particular in terms of $3$-forms, octonions, vector fields, and geometric structures. We will give an overview of each approach, the results obtained, and compare the different perspectives.


title: Inefficient magnetic-field amplification in supersonic laser-plasma turbulence

abstract: We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic plasma with a large magnetic Reynolds number ($\mathrm{Rm} \approx 45$) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller than the outer scale of the turbulent cascade was also observed. Our results support the notion that moderately supersonic, low-magnetic-Prandtl-number plasma turbulence is inefficient at amplifying magnetic fields.


title: Decision-making at Unsignalized Intersection for Autonomous Vehicles: Left-turn Maneuver with Deep Reinforcement Learning

abstract: Decision-making module enables autonomous vehicles to reach appropriate maneuvers in the complex urban environments, especially the intersection situations. This work proposes a deep reinforcement learning (DRL) based left-turn decision-making framework at unsignalized intersection for autonomous vehicles. The objective of the studied automated vehicle is to make an efficient and safe left-turn maneuver at a four-way unsignalized intersection. The exploited DRL methods include deep Q-learning (DQL) and double DQL. Simulation results indicate that the presented decision-making strategy could efficaciously reduce the collision rate and improve transport efficiency. This work also reveals that the constructed left-turn control structure has a great potential to be applied in real-time.


title: A Note on the Likelihood Ratio Test in High-Dimensional Exploratory Factor Analysis

abstract: The likelihood ratio test is widely used in exploratory factor analysis to assess the model fit and determine the number of latent factors. Despite its popularity and clear statistical rationale, researchers have found that when the dimension of the response data is large compared to the sample size, the classical chi-square approximation of the likelihood ratio test statistic often fails. Theoretically, it has been an open problem when such a phenomenon happens as the dimension of data increases; practically, the effect of high dimensionality is less examined in exploratory factor analysis, and there lacks a clear statistical guideline on the validity of the conventional chi-square approximation. To address this problem, we investigate the failure of the chi-square approximation of the likelihood ratio test in high-dimensional exploratory factor analysis, and derive the necessary and sufficient condition to ensure the validity of the chi-square approximation. The results yield simple quantitative guidelines to check in practice and would also provide useful statistical insights into the practice of exploratory factor analysis.


title: Weakly supervised cross-domain alignment with optimal transport

abstract: Cross-domain alignment between image objects and text sequences is key to many visual-language tasks, and it poses a fundamental challenge to both computer vision and natural language processing. This paper investigates a novel approach for the identification and optimization of fine-grained semantic similarities between image and text entities, under a weakly-supervised setup, improving performance over state-of-the-art solutions. Our method builds upon recent advances in optimal transport (OT) to resolve the cross-domain matching problem in a principled manner. Formulated as a drop-in regularizer, the proposed OT solution can be efficiently computed and used in combination with other existing approaches. We present empirical evidence to demonstrate the effectiveness of our approach, showing how it enables simpler model architectures to outperform or be comparable with more sophisticated designs on a range of vision-language tasks.


title: A Stochastic Control Approach to Defined Contribution Plan Decumulation: "The Nastiest, Hardest Problem in Finance"

abstract: We pose the decumulation strategy for a Defined Contribution (DC) pension plan as a problem in optimal stochastic control. The controls are the withdrawal amounts and the asset allocation strategy. We impose maximum and minimum constraints on the withdrawal amounts, and impose no-shorting no-leverage constraints on the asset allocation strategy. Our objective function measures reward as the expected total withdrawals over the decumulation horizon, and risk is measured by Expected Shortfall (ES) at the end of the decumulation period. We solve the stochastic control problem numerically, based on a parametric model of market stochastic processes. We find that, compared to a fixed constant withdrawal strategy, with minimum withdrawal set to the constant withdrawal amount, the optimal strategy has a significantly higher expected average withdrawal, at the cost of a very small increase in ES risk. Tests on bootstrapped resampled historical market data indicate that this strategy is robust to parametric model misspecification.


title: Wikidata on MARS

abstract: Multi-attributed relational structures (MARSs) have been proposed as a formal data model for generalized property graphs, along with multi-attributed rule-based predicate logic (MARPL) as a useful rule-based logic in which to write inference rules over property graphs. Wikidata can be modelled in an extended MARS that adds the (imprecise) datatypes of Wikidata. The rules of inference for the Wikidata ontology can be modelled as a MARPL ontology, with extensions to handle the Wikidata datatypes and functions over these datatypes. Because many Wikidata qualifiers should participate in most inference rules in Wikidata a method of implicitly handling qualifier values on a per-qualifier basis is needed to make this modelling useful. The meaning of Wikidata is then the extended MARS that is the closure of running these rules on the Wikidata data model. Wikidata constraints can be modelled as multi-attributed predicate logic (MAPL) formulae, again extended with datatypes, that are evaluated over this extended MARS. The result models Wikidata in a way that fixes several of its major problems.


title: Impact of a course transformation on students' reasoning about measurement uncertainty

abstract: Physics lab courses are integral parts of an undergraduate physics education, and offer a variety of opportunities for learning. Many of these opportunities center around a common learning goal in introductory physics lab courses: measurement uncertainty. Accordingly, when the stand-alone introductory lab course at the University of Colorado Boulder (CU) was recently transformed, measurement uncertainty was the focus of a learning goal of that transformation. The Physics Measurement Questionnaire (PMQ), a research-based assessment of student understanding around statistical measurement uncertainty, was used to measure the effectiveness of that transformation. Here, we analyze student responses to the PMQ at the beginning and end of the CU course. We also compare such responses from two semesters: one before and one after the transformation. We present evidence that students in both semesters shifted their reasoning in ways aligned with the measurement uncertainty learning goal. Furthermore, we show that more students in the transformed semester shifted in ways aligned with the learning goal, and that those students tended to communicate their reasoning with greater sophistication than students in the original course. These findings provide evidence that even a traditional lab course can support valuable learning, and that transforming such a course to align with well-defined learning goals can result in even more effective learning experiences.


title: Intelligent Service Selection in a Multi-dimensional Environment of Cloud Providers for IoT stream Data through cloudlets

abstract: The expansion of the Internet of Things(IoT) services and a huge amount of data generated by different sensors, signify the importance of cloud computing services like Storage as a Service more than ever. IoT traffic imposes such extra constraints on the cloud storage service as sensor data preprocessing capability and load-balancing between data centers and servers in each data center. Also, it should be allegiant to the Quality of Service (QoS). The hybrid MWG algorithm has been proposed in this work, which considers different objectives such as energy, processing time, transmission time, and load balancing in both Fog and Cloud Layer. The MATLAB script is used to simulate and implement our algorithms, and services of different servers, e.g. Amazon, Dropbox, Google Drive, etc. have been considered. The MWG has 7%, 13%, and 25% improvement in comparison with MOWCA, KGA, and NSGAII in metric of spacing, respectively. Moreover, the MWG has 4%, 4.7%, and 7.3% optimization in metric of quality in comparison to MOWCA, KGA, and NSGAII, respectively. The overall optimization shows that the MWG algorithm has 7.8%, 17%, and 21.6% better performance in comparison with MOWCA, KGA, and NSGAII in the obtained best result by considering different objectives, respectively.


title: Probing the Dark Matter of Three-loop Radiative Neutrino Mass Generation Model with the Cherenkov Telescope Array

abstract: We investigate the prospect of detecting the Dark Matter (DM) candidate in the three-loop radiative neutrino mass generation model extended with large electroweak multiplets of the Standard Model (SM) gauge group, at the future imaging atmospheric Cherenkov telescope known as the Cherenkov Telescope Array (CTA). We find that the addition of such large electroweak multiplets leads to a sizable Sommerfeld enhanced annihilation of the DM with $O(\text{TeV})$ mass, into the SM gauge bosons which results in continuum and line-like spectra of very high energy (VHE) gamma-rays, and therefore becomes observable for the CTA. We determine the viable models by setting the upper limit on the $SU(2)_{L}$ isospin of the multiplets from the partial-wave unitarity constraints and the appearance of low-scale Landau pole in the gauge coupling. Afterwards, by considering the continuum VHE gamma-rays produced from the DM annihilation at the Galactic center, we probe the parameter space of the model using the sensitivity reach of the CTA.


title: X-ray photolysis of CH$_3$COCH$_3$ ice: Implications for the radiation effects of compact objects towards astrophysical ices

abstract: In this study, we employed broadband X-rays ($6-2000$ eV) to irradiate the frozen acetone CH$_3$COCH$_3$, at the temperature of 12 K, with different photon fluences up to $2.7\times 10^{18}$ photons cm$^{-2}$. Here, we consider acetone as a representative complex organic molecule (COM) present on interstellar ice grains. The experiments were conduced at the Brazilian synchrotron facility (LNLS/CNPEN) employing infrared spectroscopy (FTIR) to monitor chemical changes induced by radiation in the ice sample. We determined the effective destruction cross-section of the acetone molecule and the effective formation cross-section for daughter species. Chemical equilibrium, obtained for fluence $2\times 10^{18}$ photons cm$^{-2}$, and molecular abundances at this stage were determined, which also includes the estimates for the abundance of unknown molecules, produced but not detected, in the ice. Timescales for ices, at hypothetical snow line distances, to reach chemical equilibrium around several compact and main-sequence X-ray sources are given. We estimate timescales of 18 days, 3.6 and 1.8 months, $1.4\times 10^9-6\times 10^{11}$ years, 600 and $1.2\times 10^7$ years, and $10^7$ years, for the Sun at 5 AU, for O/B stars at 5 AU, for white dwarfs at 1 LY, for the Crab pulsar at 2.25 LY, for Vela pulsar at 2.25 LY, and for Sagittarius A* at 3 LY, respectively. This study improves our current understanding about radiation effects on the chemistry of frozen material, in particular, focusing for the first time, the effects of X-rays produced by compact objects in their eventual surrounding ices.

categories: astro-ph.SR astro-ph.EP astro-ph.IM physics.space-ph


title: Model-Free Optimal Control of Linear Multi-Agent Systems via Decomposition and Hierarchical Approximation

abstract: Designing the optimal linear quadratic regulator (LQR) for a large-scale multi-agent system (MAS) is time-consuming since it involves solving a large-size matrix Riccati equation. The situation is further exasperated when the design needs to be done in a model-free way using schemes such as reinforcement learning (RL). To reduce this computational complexity, we decompose the large-scale LQR design problem into multiple smaller-size LQR design problems. We consider the objective function to be specified over an undirected graph, and cast the decomposition as a graph clustering problem. The graph is decomposed into two parts, one consisting of independent clusters of connected components, and the other containing edges that connect different clusters. Accordingly, the resulting controller has a hierarchical structure, consisting of two components. The first component optimizes the performance of each independent cluster by solving the smaller-size LQR design problem in a model-free way using an RL algorithm. The second component accounts for the objective coupling different clusters, which is achieved by solving a least squares problem in one shot. Although suboptimal, the hierarchical controller adheres to a particular structure as specified by inter-agent couplings in the objective function and by the decomposition strategy. Mathematical formulations are established to find a decomposition that minimizes the number of required communication links or reduces the optimality gap. Numerical simulations are provided to highlight the pros and cons of the proposed designs.


title: On the importance of being primitive

abstract: We give a brief survey of primitivity in ring theory and in particular look at characterizations of primitive ideals in the prime spectrum for various classes of rings.


title: Quantification of BERT Diagnosis Generalizability Across Medical Specialties Using Semantic Dataset Distance

abstract: Deep learning models in healthcare may fail to generalize on data from unseen corpora. Additionally, no quantitative metric exists to tell how existing models will perform on new data. Previous studies demonstrated that NLP models of medical notes generalize variably between institutions, but ignored other levels of healthcare organization. We measured SciBERT diagnosis sentiment classifier generalizability between medical specialties using EHR sentences from MIMIC-III. Models trained on one specialty performed better on internal test sets than mixed or external test sets (mean AUCs 0.92, 0.87, and 0.83, respectively; p = 0.016). When models are trained on more specialties, they have better test performances (p < 1e-4). Model performance on new corpora is directly correlated to the similarity between train and test sentence content (p < 1e-4). Future studies should assess additional axes of generalization to ensure deep learning models fulfil their intended purpose across institutions, specialties, and practices.


title: Self-supervised Contrastive Video-Speech Representation Learning for Ultrasound

abstract: In medical imaging, manual annotations can be expensive to acquire and sometimes infeasible to access, making conventional deep learning-based models difficult to scale. As a result, it would be beneficial if useful representations could be derived from raw data without the need for manual annotations. In this paper, we propose to address the problem of self-supervised representation learning with multi-modal ultrasound video-speech raw data. For this case, we assume that there is a high correlation between the ultrasound video and the corresponding narrative speech audio of the sonographer. In order to learn meaningful representations, the model needs to identify such correlation and at the same time understand the underlying anatomical features. We designed a framework to model the correspondence between video and audio without any kind of human annotations. Within this framework, we introduce cross-modal contrastive learning and an affinity-aware self-paced learning scheme to enhance correlation modelling. Experimental evaluations on multi-modal fetal ultrasound video and audio show that the proposed approach is able to learn strong representations and transfers well to downstream tasks of standard plane detection and eye-gaze prediction.