title: Optimising a magnitude-limited spectroscopic training sample for photometric classification of supernovae

abstract: In preparation for photometric classification of transients from the Legacy Survey of Space and Time (LSST) we run tests with different training data sets. Using estimates of the depth to which the 4-metre Multi-Object Spectroscopic Telescope (4MOST) Time Domain Extragalactic Survey (TiDES) can classify transients, we simulate a magnitude-limited sample reaching $r_{\textrm{AB}} \approx$ 22.5 mag. We run our simulations with the software snmachine, a photometric classification pipeline using machine learning. The machine-learning algorithms struggle to classify supernovae when the training sample is magnitude-limited, in contrast to representative training samples. Classification performance noticeably improves when we combine the magnitude-limited training sample with a simulated realistic sample of faint, high-redshift supernovae observed from larger spectroscopic facilities; the algorithms' range of average area under ROC curve (AUC) scores over 10 runs increases from 0.547-0.628 to 0.946-0.969 and purity of the classified sample reaches 95 per cent in all runs for 2 of the 4 algorithms. By creating new, artificial light curves using the augmentation software avocado, we achieve a purity in our classified sample of 95 per cent in all 10 runs performed for all machine-learning algorithms considered. We also reach a highest average AUC score of 0.986 with the artificial neural network algorithm. Having `true' faint supernovae to complement our magnitude-limited sample is a crucial requirement in optimisation of a 4MOST spectroscopic sample. However, our results are a proof of concept that augmentation is also necessary to achieve the best classification results.


title: Machine Learning Algorithm for NLOS Millimeter Wave in 5G V2X Communication

abstract: The 5G vehicle-to-everything (V2X) communication for autonomous and semi-autonomous driving utilizes the wireless technology for communication and the Millimeter Wave bands are widely implemented in this kind of vehicular network application. The main purpose of this paper is to broadcast the messages from the mmWave Base Station to vehicles at LOS (Line-of-sight) and NLOS (Non-LOS). Relay using Machine Learning (RML) algorithm is formulated to train the mmBS for identifying the blockages within its coverage area and broadcast the messages to the vehicles at NLOS using a LOS nodes as a relay. The transmission of information is faster with higher throughput and it covers a wider bandwidth which is reused, therefore when performing machine learning within the coverage area of mmBS most of the vehicles in NLOS can be benefited. A unique method of relay mechanism combined with machine learning is proposed to communicate with mobile nodes at NLOS.


title: Non-autoregressive electron flow generation for reaction prediction

abstract: Reaction prediction is a fundamental problem in computational chemistry. Existing approaches typically generate a chemical reaction by sampling tokens or graph edits sequentially, conditioning on previously generated outputs. These autoregressive generating methods impose an arbitrary ordering of outputs and prevent parallel decoding during inference. We devise a novel decoder that avoids such sequential generating and predicts the reaction in a Non-Autoregressive manner. Inspired by physical-chemistry insights, we represent edge edits in a molecule graph as electron flows, which can then be predicted in parallel. To capture the uncertainty of reactions, we introduce latent variables to generate multi-modal outputs. Following previous works, we evaluate our model on USPTO MIT dataset. Our model achieves both an order of magnitude lower inference latency, with state-of-the-art top-1 accuracy and comparable performance on Top-K sampling.


title: Deep Learning of Cell Classification using Microscope Images of Intracellular Microtubule Networks

abstract: Microtubule networks (MTs) are a component of a cell that may indicate the presence of various chemical compounds and can be used to recognize properties such as treatment resistance. Therefore, the classification of MT images is of great relevance for cell diagnostics. Human experts find it particularly difficult to recognize the levels of chemical compound exposure of a cell. Improving the accuracy with automated techniques would have a significant impact on cell therapy. In this paper we present the application of Deep Learning to MT image classification and evaluate it on a large MT image dataset of animal cells with three degrees of exposure to a chemical agent. The results demonstrate that the learned deep network performs on par or better at the corresponding cell classification task than human experts. Specifically, we show that the task of recognizing different levels of chemical agent exposure can be handled significantly better by the neural network than by human experts.


title: Structure and Complexity of Bag Consistency

abstract: Since the early days of relational databases, it was realized that acyclic hypergraphs give rise to database schemas with desirable structural and algorithmic properties. In a by-now classical paper, Beeri, Fagin, Maier, and Yannakakis established several different equivalent characterizations of acyclicity; in particular, they showed that the sets of attributes of a schema form an acyclic hypergraph if and only if the local-to-global consistency property for relations over that schema holds, which means that every collection of pairwise consistent relations over the schema is globally consistent. Even though real-life databases consist of bags (multisets), there has not been a study of the interplay between local consistency and global consistency for bags. We embark on such a study here and we first show that the sets of attributes of a schema form an acyclic hypergraph if and only if the local-to global consistency property for bags over that schema holds. After this, we explore algorithmic aspects of global consistency for bags by analyzing the computational complexity of the global consistency problem for bags: given a collection of bags, are these bags globally consistent? We show that this problem is in NP, even when the schema is part of the input. We then establish the following dichotomy theorem for fixed schemas: if the schema is acyclic, then the global consistency problem for bags is solvable in polynomial time, while if the schema is cyclic, then the global consistency problem for bags is NP-complete. The latter result contrasts sharply with the state of affairs for relations, where, for each fixed schema, the global consistency problem for relations is solvable in polynomial time.


title: The Control of Lung Branching Morphogenesis

abstract: Branching morphogenesis generates epithelial trees which facilitate gas exchange, filtering, as well as secretion processes with their large surface to volume ratio. In this review, we focus on the developmental mechanisms that control the early stages of lung branching morphogenesis. Lung branching morphogenesis involves the stereotypic, recurrent definition of new branch points, subsequent epithelial budding, and lung tube elongation. We discuss current models and experimental evidence for each of these steps. Finally, we discuss the role of the mesenchyme in determining the organ-specific shape.


title: All-sky search in early O3 LIGO data for continuous gravitational-wave signals from unknown neutron stars in binary systems

abstract: Rapidly spinning neutron stars are promising sources of persistent, continuous gravitational waves. Detecting such a signal would allow probing of the physical properties of matter under extreme conditions. A significant fraction of the known pulsar population belongs to binary systems. Searching for unknown neutron stars in binary systems requires specialized algorithms to address unknown orbital frequency modulations. We present a search for continuous gravitational waves emitted by neutron stars in binary systems in early data from the third observing run of the Advanced LIGO and Advanced Virgo detectors using the semicoherent, GPU-accelerated, BinarySkyHough pipeline. The search analyzes the most sensitive frequency band of the LIGO detectors, 50 - 300 Hz. Binary orbital parameters are split into four regions, comprising orbital periods of 3 - 45 days and projected semimajor axes of 2 - 40 light-seconds. No detections are reported. We estimate the sensitivity of the search using simulated continuous wave signals, achieving the most sensitive results to date across the analyzed parameter space.


title: Asymptotic optimality of the triangular lattice for a class of optimal location problems

abstract: We prove an asymptotic crystallization result in two dimensions for a class of nonlocal particle systems. To be precise, we consider the best approximation with respect to the 2-Wasserstein metric of a given absolutely continuous probability measure $f \mathrm{d}x$ by a discrete probability measure $\sum_i m_i \delta_{z_i}$, subject to a constraint on the particle sizes $m_i$. The locations $z_i$ of the particles, their sizes $m_i$, and the number of particles are all unknowns of the problem. We study a one-parameter family of constraints. This is an example of an optimal location problem (or an optimal sampling or quantization problem) and it has applications in economics, signal compression, and numerical integration. We establish the asymptotic minimum value of the (rescaled) approximation error as the number of particles goes to infinity. In particular, we show that for the constrained best approximation of the Lebesgue measure by a discrete measure, the discrete measure whose support is a triangular lattice is asymptotically optimal. In addition, we prove an analogous result for a problem where the constraint is replaced by a penalization. These results can also be viewed as the asymptotic optimality of the hexagonal tiling for an optimal partitioning problem. They generalise the crystallization result of Bourne, Peletier and Theil (Communications in Mathematical Physics, 2014) from a single particle system to a class of particle systems, and prove a case of a conjecture by Bouchitt\'{e}, Jimenez and Mahadevan (Journal de Math\'ematiques Pures et Appliqu\'ees, 2011). Finally, we prove a crystallization result which states that optimal configurations with energy close to that of a triangular lattice are geometrically close to a triangular lattice.


title: A Computational Framework for Solving Nonlinear Binary OptimizationProblems in Robust Causal Inference

abstract: Identifying cause-effect relations among variables is a key step in the decision-making process. While causal inference requires randomized experiments, researchers and policymakers are increasingly using observational studies to test causal hypotheses due to the wide availability of observational data and the infeasibility of experiments. The matching method is the most used technique to make causal inference from observational data. However, the pair assignment process in one-to-one matching creates uncertainty in the inference because of different choices made by the experimenter. Recently, discrete optimization models are proposed to tackle such uncertainty. Although a robust inference is possible with discrete optimization models, they produce nonlinear problems and lack scalability. In this work, we propose greedy algorithms to solve the robust causal inference test instances from observational data with continuous outcomes. We propose a unique framework to reformulate the nonlinear binary optimization problems as feasibility problems. By leveraging the structure of the feasibility formulation, we develop greedy schemes that are efficient in solving robust test problems. In many cases, the proposed algorithms achieve global optimal solutions. We perform experiments on three real-world datasets to demonstrate the effectiveness of the proposed algorithms and compare our result with the state-of-the-art solver. Our experiments show that the proposed algorithms significantly outperform the exact method in terms of computation time while achieving the same conclusion for causal tests. Both numerical experiments and complexity analysis demonstrate that the proposed algorithms ensure the scalability required for harnessing the power of big data in the decision-making process.


title: The compression semigroup of the dual Vinberg cone

abstract: We investigate the semigroup associated to the dual Vinberg cone and prove its triple and Ol'shanski\u{\i} polar decompositions. Moreover, we show that the semigroup does not have the contraction property with respect to the canonical Riemannian metric on the cone.


title: MOVES IV. Modelling the influence of stellar XUV-flux, cosmic rays, and stellar energetic particles on the atmospheric composition of the hot Jupiter HD 189733b

abstract: Hot Jupiters provide valuable natural laboratories for studying potential contributions of high-energy radiation to prebiotic synthesis in the atmospheres of exoplanets. In this fourth paper of the MOVES (Multiwavelength Observations of an eVaporating Exoplanet and its Star) programme, we study the effect of different types of high-energy radiation on the production of organic and prebiotic molecules in the atmosphere of the hot Jupiter HD 189733b. Our model combines X-ray and UV observations from the MOVES programme and 3D climate simulations from the 3D Met Office Unified Model to simulate the atmospheric composition and kinetic chemistry with the STAND2019 network. Also, the effects of galactic cosmic rays and stellar energetic particles are included. We find that the differences in the radiation field between the irradiated dayside and the shadowed nightside lead to stronger changes in the chemical abundances than the variability of the host star's XUV emission. We identify ammonium (NH4+) and oxonium (H3O+) as fingerprint ions for the ionization of the atmosphere by both galactic cosmic rays and stellar particles. All considered types of high-energy radiation have an enhancing effect on the abundance of key organic molecules such as hydrogen cyanide (HCN), formaldehyde (CH2O), and ethylene (C2H4). The latter two are intermediates in the production pathway of the amino acid glycine (C2H5NO2) and abundant enough to be potentially detectable by JWST.

categories: astro-ph.EP astro-ph.HE astro-ph.SR physics.chem-ph


title: Finitely-valued propositional dynamic logic

abstract: We study a many-valued generalization of Propositional Dynamic Logic where formulas in states and accessibility relations between states of a Kripke model are evaluated in a finite FL-algebra. One natural interpretation of this framework is related to reasoning about costs of performing structured actions. We prove that PDL over any finite FL-algebra is decidable. We also establish a general completeness result for a class of PDLs based on commutative integral FL-algebras with canonical constants.


title: The effect of doping on the lattice parameter and properties of cubic boron nitride

abstract: The effect of doping of cubic boron nitride with beryllium, silicon, sulfur and magnesium on the lattice parameters, electrical conductivity and EPR spectra has been studied. It is established that the degree of doping increases significantly in the case of crystallization of cubic boron nitride from BN solutions in supercritical ammonia at 3.9-4.2 GPa and 1100{\deg}C in comparison with the conventional synthesis from melts of the Mg-B-N system at 4.2 GPa and 1400{\deg}C. Doping with silicon and beryllium results in semiconductor properties of cubic boron nitride.


title: COVID-19: Optimal Design of Serosurveys for Disease Burden Estimation

abstract: We provide a methodology by which an epidemiologist may arrive at an optimal design for a survey whose goal is to estimate the disease burden in a population. For serosurveys with a given budget of $C$ rupees, a specified set of tests with costs, sensitivities, and specificities, we show the existence of optimal designs in four different contexts, including the well known c-optimal design. Usefulness of the results are illustrated via numerical examples. Our results are applicable to a wide range of epidemiological surveys under the assumptions that the estimate's Fisher-information matrix satisfies a uniform positive definite criterion.


title: Highly tunable charge transport in defective graphene nanoribbons under external local forces and constraints: A hybrid computational study

abstract: In this paper, we propose a combined modeling of molecular mechanics (MM) and the tight-binding (TB) approach, which enables us to study the effect of factors such as external local forces, constraints, and vacancy defects on electronic transport properties of nanomaterials. Nanostructures selected in this work are armchair graphene nanoribbons (AGNRs). According to this method, the nanostructure is modeled as a frame, and the beam element is applied for illustrating the covalent interatomic interactions in bonds. In our calculations, the terms of torsional, stretching, and bending energies are considered. The selected pristine nanoribbon is a metal, and the purpose of this study is to find the effects of mechanical loading, the vacancy defects and their positions on the electrical conductance of the structure. We observe that the presence of vacancy defects in the structure leads to the opening of an energy gap, which changes the phase of the nanostructure from metal to the semiconductor. We find that with increasing the number of point defects, the energy gap size of the strained system grows. Besides, increasing the magnitude of the local force reduces the conductance, and the energy gap of the system. By changing parameters such as the number of point defects and magnitude local forces, the transport gap of the system can be controlled. The results of this research may be useful in the design of nanoelectromechanical systems.


title: Projected Stochastic Gradient Langevin Algorithms for Constrained Sampling and Non-Convex Learning

abstract: Langevin algorithms are gradient descent methods with additive noise. They have been used for decades in Markov chain Monte Carlo (MCMC) sampling, optimization, and learning. Their convergence properties for unconstrained non-convex optimization and learning problems have been studied widely in the last few years. Other work has examined projected Langevin algorithms for sampling from log-concave distributions restricted to convex compact sets. For learning and optimization, log-concave distributions correspond to convex losses. In this paper, we analyze the case of non-convex losses with compact convex constraint sets and IID external data variables. We term the resulting method the projected stochastic gradient Langevin algorithm (PSGLA). We show the algorithm achieves a deviation of $O(T^{-1/4}(\log T)^{1/2})$ from its target distribution in 1-Wasserstein distance. For optimization and learning, we show that the algorithm achieves $\epsilon$-suboptimal solutions, on average, provided that it is run for a time that is polynomial in $\epsilon^{-1}$ and slightly super-exponential in the problem dimension.


title: Projection-Free Bandit Optimization with Privacy Guarantees

abstract: We design differentially private algorithms for the bandit convex optimization problem in the projection-free setting. This setting is important whenever the decision set has a complex geometry, and access to it is done efficiently only through a linear optimization oracle, hence Euclidean projections are unavailable (e.g. matroid polytope, submodular base polytope). This is the first differentially-private algorithm for projection-free bandit optimization, and in fact our bound of $\widetilde{O}(T^{3/4})$ matches the best known non-private projection-free algorithm (Garber-Kretzu, AISTATS `20) and the best known private algorithm, even for the weaker setting when projections are available (Smith-Thakurta, NeurIPS `13).


title: Image to Bengali Caption Generation Using Deep CNN and Bidirectional Gated Recurrent Unit

abstract: There is very little notable research on generating descriptions of the Bengali language. About 243 million people speak in Bengali, and it is the 7th most spoken language on the planet. The purpose of this research is to propose a CNN and Bidirectional GRU based architecture model that generates natural language captions in the Bengali language from an image. Bengali people can use this research to break the language barrier and better understand each other's perspectives. It will also help many blind people with their everyday lives. This paper used an encoder-decoder approach to generate captions. We used a pre-trained Deep convolutional neural network (DCNN) called InceptonV3image embedding model as the encoder for analysis, classification, and annotation of the dataset's images Bidirectional Gated Recurrent unit (BGRU) layer as the decoder to generate captions. Argmax and Beam search is used to produce the highest possible quality of the captions. A new dataset called BNATURE is used, which comprises 8000 images with five captions per image. It is used for training and testing the proposed model. We obtained BLEU-1, BLEU-2, BLEU-3, BLEU-4 and Meteor is 42.6, 27.95, 23, 66, 16.41, 28.7 respectively.


title: Betti structures of hypergeometric equations

abstract: We study Betti structures in the solution complexes of confluent hypergeometric equations. We use the framework of enhanced ind-sheaves and the irregular Riemann-Hilbert correspondence of D'Agnolo-Kashiwara. The main result is a group theoretic criterion that ensures that enhanced solutions of such systems are defined over certain subfields of the complex numbers. The proof uses a description of the hypergeometric systems as exponentially twisted Gauss-Manin systems of certain Laurent polynomials.


title: Learning to Initialize Gradient Descent Using Gradient Descent

abstract: Non-convex optimization problems are challenging to solve; the success and computational expense of a gradient descent algorithm or variant depend heavily on the initialization strategy. Often, either random initialization is used or initialization rules are carefully designed by exploiting the nature of the problem class. As a simple alternative to hand-crafted initialization rules, we propose an approach for learning "good" initialization rules from previous solutions. We provide theoretical guarantees that establish conditions that are sufficient in all cases and also necessary in some under which our approach performs better than random initialization. We apply our methodology to various non-convex problems such as generating adversarial examples, generating post hoc explanations for black-box machine learning models, and allocating communication spectrum, and show consistent gains over other initialization techniques.


title: High-Speed Robot Navigation using Predicted Occupancy Maps

abstract: Safe and high-speed navigation is a key enabling capability for real world deployment of robotic systems. A significant limitation of existing approaches is the computational bottleneck associated with explicit mapping and the limited field of view (FOV) of existing sensor technologies. In this paper, we study algorithmic approaches that allow the robot to predict spaces extending beyond the sensor horizon for robust planning at high speeds. We accomplish this using a generative neural network trained from real-world data without requiring human annotated labels. Further, we extend our existing control algorithms to support leveraging the predicted spaces to improve collision-free planning and navigation at high speeds. Our experiments are conducted on a physical robot based on the MIT race car using an RGBD sensor where were able to demonstrate improved performance at 4 m/s compared to a controller not operating on predicted regions of the map.


title: Introducing CPL

abstract: CPL here stands for a computer programming language conceived and developed by the author since 1993, but published for the first time in 2020. It was born as a Compiled Programming Language, designed together with its compiler and therefore suitable for computationally intensive numerical applications, although some years later an interpreter was also provided for interactive usage. CPL's distinctive features are Concealed Pointer Lookup, the ability to implicitly dereference pointers based on the type of operands involved, Consistent Procedure Linkage, the enforcement of function prototypes without dedicated header or interface files, and Coactive Parameter Lists, the ability to overload function names which are then distinguished by the type of their parameters and/or parameter separators. Perhaps even more distinctly, CPL's syntax can be extended on the fly by the program being compiled; library modules tap this feature to seamlessly add real and complex matrix operations, graphics, parallel-computing extensions, and symbolic differentiation. The CPL coding software is available for free download at http://CPLcode.net .


title: Integrating computing in the statistics and data science curriculum: Creative structures, novel skills and habits, and ways to teach computational thinking

abstract: Nolan and Temple Lang (2010) argued for the fundamental role of computing in the statistics curriculum. In the intervening decade the statistics education community has acknowledged that computational skills are as important to statistics and data science practice as mathematics. There remains a notable gap, however, between our intentions and our actions. In this special issue of the *Journal of Statistics and Data Science Education* we have assembled a collection of papers that (1) suggest creative structures to integrate computing, (2) describe novel data science skills and habits, and (3) propose ways to teach computational thinking. We believe that it is critical for the community to redouble our efforts to embrace sophisticated computing in the statistics and data science curriculum. We hope that these papers provide useful guidance for the community to move these efforts forward.


title: A COLREGs-Compliant Motion Planner for Autonomous Maneuvering of Marine Vessels in Complex Environments

abstract: An enabling technology for future sea transports is safe and energy-efficient autonomous maritime navigation in narrow environments with other marine vessels present. This requires that the algorithm controlling the ship is able to account for the vessel's dynamics, and obeys the rules specified in the international regulations for preventing collision at sea (COLREGs). To account for these properties, this work proposes a two-step optimization-based and COLREGs-compliant motion planner tailored for large vessels. In the first step, a lattice-based motion planner is used to compute a suboptimal trajectory based on a library of precomputed motion primitives. To comply with the rules specified in COLREGs, the lattice-based planner is augmented with discrete states that represent what type of COLREGs situations that are active with respect to other nearby vessels. The trajectory computed by the lattice-based planner is then used as warm-start for the second receding-horizon step based on direct optimal control techniques. The aim of the second step is to solve the problem to local optimality with the discrete COLREGs states remained fixed. The applicability of the proposed motion planner is demonstrated in simulations, where it computes energy-efficient and COLREGs-compliant trajectories. The results also show that the motion planner is able to prevent complex collision situations from occurring.


title: Mutually interacting superprocesses with migration

abstract: A system of mutually interacting superprocesses with migration is constructed as the limit of a sequence of branching particle systems arising from population models. The uniqueness in law of the superprocesses is established using the pathwise uniqueness of a system of stochastic partial differential equations with non-Lipschitz coefficients, which is satisfied by the corresponding system of distribution-function-valued processes.


title: Another presentation of orthogonal Steinberg groups

abstract: We use the pro-group approach to show that $\mathrm{StO}(M, q)$ admits van der Kallen's "another presentation", where $M$ is a module over a commutative ring with sufficiently isotropic quadratic form $q$. Moreover, we construct an analog of ESD-transvections in orthogonal Steinberg pro-groups under some assumptions on their parameters.


title: Cabling Legendrian and transverse knots

abstract: In this paper we will show how to classify Legendrian and transverse knots in the knot type of "sufficiently positive" cables of a knot in terms of the classification of the underlying knot. We will also completely explain the phenomena of "Legendrian large" cables. These are Legendrian representatives of cables that have Thurston-Bennequin invariant larger that the framing coming from the cabling torus. Such examples have only recently, and unexpectedly, been found. We will also give criteria that determines the classification of Legendrian and transverse knots the the knot type of negative cables.


title: Hierarchically nanostructured thermoelectric materials: Challenges and opportunities for improved power factors

abstract: The field of thermoelectric materials has undergone a revolutionary transformation over the last couple of decades as a result of the ability to nanostructure and synthesize myriads of materials and their alloys. The ZT figure of merit, which quantifies the performance of a thermoelectric material has more than doubled after decades of inactivity, reaching values larger than two, consistently across materials and temperatures. Central to this ZT improvement is the drastic reduction in the material thermal conductivity due to the scattering of phonons on the numerous interfaces, boundaries, dislocations, point defects, phases, etc., which are purposely included. In these new generation of nanostructured materials, phonon scattering centers of different sizes and geometrical configurations (atomic, nano- and macro-scale) are formed, which are able to scatter phonons of mean-free-paths across the spectrum. Beyond thermal conductivity reductions, ideas are beginning to emerge on how to use similar hierarchical nanostructuring to achieve power factor improvements. Ways that relax the adverse interdependence of the electrical conductivity and Seebeck coefficient are targeted, which allows power factor improvements. For this, elegant designs are required, that utilize for instance non-uniformities in the underlying nanostructured geometry, non-uniformities in the dopant distribution, or potential barriers that form at boundaries between materials. A few recent reports, both theoretical and experimental, indicate that extremely high power factor values can be achieved, even for the same geometries that also provide ultra-low thermal conductivities. Despite the experimental complications that can arise in having the required control in nanostructure realization, in this colloquium, we aim to demonstrate, mostly theoretically, that it is a very promising path worth exploring.


title: Numerical approximation of singular-degenerate parabolic stochastic PDEs

abstract: We study a general class of singular degenerate parabolic stochastic partial differential equations (SPDEs) which include, in particular, the stochastic porous medium equations and the stochastic fast diffusion equation. We propose a fully discrete numerical approximation of the considered SPDEs based on the very weak formulation. By exploiting the monotonicity properties of the proposed formulation we prove the convergence of the numerical approximation towards the unique solution. Furthermore, we construct an implementable finite element scheme for the spatial discretization of the very weak formulation and provide numerical simulations to demonstrate the practicability of the proposed discretization.


title: Updated universal relations for tidal deformabilities of neutron stars from phenomenological equations of state

abstract: Equation of state (EOS) insensitive relations, so-called universal relations, between the neutron star (NS) compactness, its multipolar tidal deformability coefficients, and between the tidal parameters for binary systems are essential to break degeneracies in gravitational wave data analysis. Here, we validate and recalibrate these universal relations using a large set of almost 2 million phenomenological EOSs that are consistent with current observations. In doing so, we extend universal relations to a larger region of the EOS parameter space, most notably to softer EOSs and larger compactnesses. We show that waveform models that neglect higher-than-leading-order tidal deformations of the NSs accumulate as much as $3.5$ radians of dephasing from $20\,{\rm Hz}$ to merger. We also perform a full Bayesian parameter estimation of the GW170817 data, and we compare the NS radius constraints produced using universal relations from the literature and the updated fits we propose here. We find that the new fits yield a NS radius that is smaller by about 500 meters. This difference is less than the statistical uncertainty on the radius at the signal-to-noise-ratio of GW170817, but it is significantly larger than the precision anticipated for next-generation detectors.


title: Photometric analysis of three totally eclipsing W UMa stars with increasing periods: TYC 3700-1384-1, V1511 Her and V1179 Her

abstract: The first multi-colour light curve models and period studies for the totally eclipsing W UMa stars TYC 3700-1384-1, V1511 Her and V1179 Her are presented. All three stars are A-subtype W UMa stars of spectral type F. The light curve solutions show that TYC 3700-1384-1 has a moderately low mass ratio of q = 0.182 +/- 0.001 and a degree of overcontact of f = 49%. For V1179 Her a mass ratio q = 0.153 +/- 0.001 and a degree of overcontact of f = 48% is derived. The solution for V1511 Her is inconclusive, however the mass ratio is expected to be between 0.13 < q < 0.15. The evolutionary status is compared with zero-age main sequence stars taking into account energy transfer from the primary to the secondary component. The primary component of TYC 3700-1384-1 fits well in the main-sequence, while V1179 Her is more evolved. The period study reveals for all three stars a continuously increasing period at a rate of dP/dt = 6.1x10-7 d/yr, dP/dt = 5.0x10-7 d/yr and dP/dt = 9.6x10-7 d/yr for TYC 3700-1384-1, V1511 Her and V1179 Her respectively. The estimated mass transfer rates derived from these period changes are 1.6x10-7 solMass/yr for TYC 3700-1384-1 and 1.9x10-7 solMass/yr for V1179 Her.


title: An Improved Algorithm for Coarse-Graining Cellular Automata

abstract: In studying the predictability of emergent phenomena in complex systems, Israeli & Goldenfeld (Phys. Rev. Lett., 2004; Phys. Rev. E, 2006) showed how to coarse-grain (elementary) cellular automata (CA). Their algorithm for finding coarse-grainings of supercell size $N$ took doubly-exponential $2^{2^N}$-time, and thus only allowed them to explore supercell sizes $N \leq 4$. Here we introduce a new, more efficient algorithm for finding coarse-grainings between any two given CA that allows us to systematically explore all elementary CA with supercell sizes up to $N=7$, and to explore individual examples of even larger supercell size. Our algorithm is based on a backtracking search, similar to the DPLL algorithm with unit propagation for the NP-complete problem of Boolean Satisfiability.

categories: nlin.CG cond-mat.stat-mech cs.DS nlin.PS


title: How dark is the dark side of diversification?

abstract: Against the widely held belief that diversification at banking institutions contributes to the stability of the financial system, Wagner (2010) found that diversification actually makes systemic crisis more likely. While it is true, as Wagner asserts, that the probability of joint default of the diversified portfolios is larger; we contend that, as common practice, the effect of diversification is examined with respect to a risk measure like VaR. We find that when banks use VaR, diversification does reduce individual and systemic risk. This, in turn, generates a different set of incentives for banks and regulators.


title: Estimation of discrete choice models with hybrid stochastic adaptive batch size algorithms

abstract: The emergence of Big Data has enabled new research perspectives in the discrete choice community. While the techniques to estimate Machine Learning models on a massive amount of data are well established, these have not yet been fully explored for the estimation of statistical Discrete Choice Models based on the random utility framework. In this article, we provide new ways of dealing with large datasets in the context of Discrete Choice Models. We achieve this by proposing new efficient stochastic optimization algorithms and extensively testing them alongside existing approaches. We develop these algorithms based on three main contributions: the use of a stochastic Hessian, the modification of the batch size, and a change of optimization algorithm depending on the batch size. A comprehensive experimental comparison of fifteen optimization algorithms is conducted across ten benchmark Discrete Choice Model cases. The results indicate that the HAMABS algorithm, a hybrid adaptive batch size stochastic method, is the best performing algorithm across the optimization benchmarks. This algorithm speeds up the optimization time by a factor of 23 on the largest model compared to existing algorithms used in practice. The integration of the new algorithms in Discrete Choice Models estimation software will significantly reduce the time required for model estimation and therefore enable researchers and practitioners to explore new approaches for the specification of choice models.


title: Bosonic impurity in a one-dimensional quantum droplet in the Bose-Bose mixture

abstract: We study an impurity immersed in the mixture of Bose ultracold gases in the regime where a quantum droplet exists. The quasi-one-dimensional geometry is considered. We find an effective attractive potential that acts by the quantum droplet onto the impurity. The bound states of the impurity in this potential are investigated. These impurity bound states can provide potential probes for the presence of quantum fluctuations effects on the droplet properties. In the case of strong impurity-BEC coupling, we study the properties of the nonlinear local modes on the impurity induced by the quantum fluctuations.


title: On Symmetry and the Reality of Holomorphic Hartree--Fock Wavefunctions

abstract: The coalescence and disappearance of Hartree--Fock (HF) solutions as the molecular structure varies have been a common source of criticism for the breakdown of the HF approximation to the potential energy surfaces. However, recent developments in holomorphic HF theory show that this disappearing behavior is only a manifestation of the way conventional HF equations prevent solutions from being analytically continued, but it is unclear what factors govern the existence and the locations of these disappearances. In this work, we explore some of these factors from the perspective of spatial symmetry by introducing a classification for symmetry constraints on electronic-structure calculations. This forms a framework for us to systematically investigate several analytic holomorphic HF solutions of a model $\textrm{[H}_4\textrm{]}^{2+}$ system in STO-3G and demonstrate that, under appropriate conditions, spatial symmetry imposes strict requirements on the reality of certain solutions. The implications for self-consistent-field HF search algorithms are then discussed. Throughout this article, the term reality means the quality of a holomorphic HF solution having real molecular orbitals.


title: Mechanical parametric feedback-cooling for pendulum-based gravity experiments

abstract: Gravitational forces that oscillate at audio-band frequencies are measured with masses suspended as pendulums that have resonance frequencies even lower. If the pendulum is excited by thermal energy or by seismic motion of the environment, the measurement sensitivity is reduced. Conventionally, this problem is mitigated by seismic isolation and linear damping, potentially combined with cryogenic cooling. Here, we propose mechanical parametric cooling of the pendulum motion during the gravitational field measurement. We report a proof of principle demonstration in the seismic noise dominated regime and achieve a damping factor of the pendulum motion of 5.7. We find a model system for which mechanical parametric feedback cooling reaches the quantum mechanical regime near the ground state. More feasible applications we anticipate in gravitational-wave detectors.


title: A Funk perspective on billiards, projective geometry and Mahler volume

abstract: We explore connections furnished by the Funk metric, a relative of the Hilbert metric, between projective geometry, billiards, convex geometry and affine inequalities. We first show that many metric invariants of the Funk metric are invariant under projective transformations as well as projective duality. These include the Holmes-Thompson volume and surface area of convex subsets, and the length spectrum of their boundary, extending results of Holmes-Thompson and \'Alvarez Paiva on Sch\"affer's dual girth conjecture. We explore in particular Funk billiards, which generalize hyperbolic billiards in the same way that Minkowski billiards generalize Euclidean ones, and extend a result of Gutkin-Tabachnikov on the duality of Minkowski billiards. We next consider the volume of outward balls in Funk geometry. We conjecture a general affine inequality corresponding to the volume maximizers, which includes the Blaschke-Santal\'o and centro-affine isoperimetric inequalities as limit cases, and prove it for unconditional bodies, yielding a new proof of the volume entropy conjecture for the Hilbert metric for unconditional bodies. As a by-product, we obtain generalizations to higher moments of inequalities of Ball and Huang-Li, which in turn strengthen the Blaschke-Santal\'o inequality for unconditional bodies. Lastly, we introduce a regularization of the total volume of a smooth strictly convex 2-dimensional set equipped with the Funk metric, resembling the O'Hara M\"obius energy of a knot, and show that it is a projective invariant of the convex body.


title: Convolutional Recurrent Network for Road Boundary Extraction

abstract: Creating high definition maps that contain precise information of static elements of the scene is of utmost importance for enabling self driving cars to drive safely. In this paper, we tackle the problem of drivable road boundary extraction from LiDAR and camera imagery. Towards this goal, we design a structured model where a fully convolutional network obtains deep features encoding the location and direction of road boundaries and then, a convolutional recurrent network outputs a polyline representation for each one of them. Importantly, our method is fully automatic and does not require a user in the loop. We showcase the effectiveness of our method on a large North American city where we obtain perfect topology of road boundaries 99.3% of the time at a high precision and recall.


title: Tropical compactification via Ganter's algorithm

abstract: We describe a canonical compactification of a polyhedral complex in Euclidean space. When the recession cones of the polyhedral complex form a fan, the compactified polyhedral complex is a subspace of a tropical toric variety. In this case, the procedure is analogous to the tropical compactifications of subvarieties of tori. We give an analysis of the combinatorial structure of the compactification and show that its Hasse diagram can be computed via Ganter's algorithm. Our algorithm is implemented in and shipped with polymake.


title: Generalization of group-theoretic coherent states for variational calculations

abstract: We introduce new families of pure quantum states that are constructed on top of the well-known Gilmore-Perelomov group-theoretic coherent states. We do this by constructing unitaries as the exponential of operators quadratic in Cartan subalgebra elements and by applying these unitaries to regular group-theoretic coherent states. This enables us to generate entanglement not found in the coherent states themselves, while retaining many of their desirable properties. Most importantly, we explain how the expectation values of physical observables can be evaluated efficiently. Examples include generalized spin-coherent states and generalized Gaussian states, but our construction can be applied to any Lie group represented on the Hilbert space of a quantum system. We comment on their applicability as variational families in condensed matter physics and quantum information.

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


title: Fast response photogating in monolayer MoS2 phototransistors

abstract: Two-dimensional transition metal dichalcogenide (TMD) phototransistors have been object of intensive research during the last years due to their potential for photodetection. Photoresponse in these devices is typically caused by a combination of two physical mechanisms: photoconductive effect (PCE) and photogating effect (PGE). In earlier literature for monolayer (1L) MoS2 phototransistors PGE is generally attributed to charge trapping by polar molecules adsorbed to the semiconductor channel, giving rise to a very slow photoresponse. Thus, the photoresponse of 1L-MoS2 phototransistors at high-frequency light modulation is assigned to PCE alone. Here we investigate the photoresponse of a fully h-BN encapsulated monolayer (1L) MoS2 phototransistor. In contrast with previous understanding, we identify a rapidly responding PGE mechanism that becomes the dominant contribution to photoresponse under high-frequency light modulation. Using a Hornbeck-Haynes model for the photocarrier dynamics, we fit the illumination power dependence of this PGE and estimate the energy level of the involved traps. The resulting energies are compatible with shallow traps in MoS2 caused by the presence of sulfur vacancies.


title: A numerical study of an Heaviside function driven degenerate diffusion equation

abstract: We analyze a nonlinear degenerate parabolic problem whose diffusion coefficient is the Heaviside function of the distance of the solution itself from a given target function. We show that this model behaves as an evolutive variational inequality having the target as an obstacle: under suitable hypotheses, starting from an initial state above the target the solution evolves in time towards an asymptotic solution, eventually getting in contact with part of the target itself. We also study a finite difference approach to the solution of this problem, using the exact Heaviside function or a regular approximation of it, showing the results of some numerical tests.


title: QCD sum rules for parameters of the $B$-meson distribution amplitudes

abstract: We obtain new estimates for the parameters $\lambda_{E}^2$, $\lambda_H^2$ and their ratio $\mathcal{R} = \lambda_{E}^2/\lambda_H^2$, which appear in the second moments of the $B$-meson light-cone distribution amplitudes defined in the heavy-quark effective field theory. The computation is based on two-point QCD sum rules for the diagonal correlation function and includes all contributions up to mass dimension seven in the operator-product expansion. For the ratio we get $\mathcal{R} = (0.1 \pm 0.1)$ with $\lambda_H^2 = (0.15 \pm 0.05) \, \text{GeV}^2$ and $\lambda_E^2 = (0.01 \pm 0.01) \, \text{GeV}^2$.


title: Unifying the Hyperbolic and Spherical 2-Body Problem with Biquaternions

abstract: The 2-body problem on the sphere and hyperbolic space are both real forms of holomorphic Hamiltonian systems defined on the complex sphere. This admits a natural description in terms of biquaternions and allows us to address questions concerning the hyperbolic system by complexifying it and treating it as the complexification of a spherical system. In this way, results for the 2-body problem on the sphere are readily translated to the hyperbolic case. For instance, we implement this idea to completely classify the relative equilibria for the 2-body problem on hyperbolic 3-space for a strictly attractive potential.


title: Sensitivity analysis in the infinite dimensional Heston model

abstract: We consider the infinite dimensional Heston stochastic volatility model proposed in \arXiv:1706:03500. The price of a forward contract on a non-storable commodity is modelled by a generalized Ornstein-Uhlenbeck process in the Filipovi\'{c} space with this volatility. We prove different representation formulas for the forward price. Then we consider prices of options written on these forward contracts and we study sensitivity analysis with computation of the Greeks with respect to different parameters in the model. Since these parameters are infinite dimensional, we need to reinterpret the meaning of the Greeks. For this we use infinite dimensional Malliavin calculus and a randomization technique.


title: Hahn polynomials for hypergeometric distribution

abstract: Orthogonal polynomials for the multivariate hypergeometric distribution are defined on lattices in polyhedral domains in $\RR^d$. Their structures are studied through a detailed analysis of classical Hahn polynomials with negative integer parameters. Factorization of the Hahn polynomials is explored and used to explain the relation between the index set of orthogonal polynomials and the lattice set in polyhedral domain. In the multivariate case, these constructions lead to nontrivial families of hypergeometric polynomials vanishing on lattice polyhedra. The generating functions and bispectral properties of the orthogonal polynomials are also discussed.


title: C-Watcher: A Framework for Early Detection of High-Risk Neighborhoods Ahead of COVID-19 Outbreak

abstract: The novel coronavirus disease (COVID-19) has crushed daily routines and is still rampaging through the world. Existing solution for nonpharmaceutical interventions usually needs to timely and precisely select a subset of residential urban areas for containment or even quarantine, where the spatial distribution of confirmed cases has been considered as a key criterion for the subset selection. While such containment measure has successfully stopped or slowed down the spread of COVID-19 in some countries, it is criticized for being inefficient or ineffective, as the statistics of confirmed cases are usually time-delayed and coarse-grained. To tackle the issues, we propose C-Watcher, a novel data-driven framework that aims at screening every neighborhood in a target city and predicting infection risks, prior to the spread of COVID-19 from epicenters to the city. In terms of design, C-Watcher collects large-scale long-term human mobility data from Baidu Maps, then characterizes every residential neighborhood in the city using a set of features based on urban mobility patterns. Furthermore, to transfer the firsthand knowledge (witted in epicenters) to the target city before local outbreaks, we adopt a novel adversarial encoder framework to learn "city-invariant" representations from the mobility-related features for precise early detection of high-risk neighborhoods, even before any confirmed cases known, in the target city. We carried out extensive experiments on C-Watcher using the real-data records in the early stage of COVID-19 outbreaks, where the results demonstrate the efficiency and effectiveness of C-Watcher for early detection of high-risk neighborhoods from a large number of cities.


title: Characteristic classes for families of bundles

abstract: The generalized Miller-Morita-Mumford classes of a manifold bundle with fiber $M$ depend only on the underlying $\tau_M$-fibration, meaning the family of vector bundles formed by the tangent bundles of the fibers. This motivates a closer study of the classifying space for $\tau_M$-fibrations, $Baut(\tau_M)$, and its cohomology ring, i.e., the ring of characteristic classes of $\tau_M$-fibrations. For a bundle $\xi$ over a simply connected Poincar\'e duality space, we construct a relative Sullivan model for the universal orientable $\xi$-fibration together with explicit cocycle representatives for the characteristic classes of the canonical bundle over its total space. This yields tools for computing the rational cohomology ring of $Baut(\xi)$ as well as the subring generated by the generalized Miller-Morita-Mumford classes. To illustrate, we carry out sample computations for spheres and complex projective spaces. We discuss applications to tautological rings of simply connected manifolds and to the problem of deciding whether a given $\tau_M$-fibration comes from a manifold bundle.


title: (3+1)-Formulation for Gravity with Torsion and Non-Metricity: The Stress-Energy-Momentum Equation

abstract: We derive the generalized Gauss-Codazzi-Mainardi (GCM) equation for a general affine connection with torsion and non-metricity. Moreover, we show that the metric compatibility and torsionless condition of a connection on a manifold are inherited to the connection of its hypersurface. As a physical application to these results, we derive the (3+1)-Einstein Field Equation (EFE) for a special case of Metric-Affine f(R)-gravity when f(R)=R, the Metric-Affine General Relativity (MAGR). Motivated by the concept of geometrodynamics, we introduce additional variables on the hypersurface as a consequence of non-vanishing torsion and non-metricity. With these additional variables, we show that for MAGR, the energy, momentum, and the stress-energy part of the EFE are dynamical, i.e., all of them contain the derivative of a quantity with respect to the time coordinate. For the Levi-Civita connection, one could recover the Hamiltonian and the momentum (diffeomorphism) constraint, and obtain the standard dynamics of GR.


title: A Photometric and Dynamic Study of Comet C/2013 A1 (Siding Spring) from Observations at a Heliocentric Distance of ~4.1 AU

abstract: An analysis is presented for the photometric data on comet C/2013 A1 (Siding Spring) from observations at a large heliocentric distance (near 4.1 AU). Comet C/2013 A1 (Siding Spring) displays intense activity despite the relatively large heliocentric distance. The morphology of the comet's coma is analyzed. The following parameters are measured: the color indices V-R, the normalized spectral gradient of the reflectivity of the comet's dust S', and the dust production rate Af\r{ho}. A numerical simulation is performed for the evolution of the comet's orbit after a close encounter with Mars. The most probable values are obtained for the Keplerian orbital elements of the comet over a hundred-year period. The comet's orbit remains nearly parabolic after passing the orbits of all the Solar System planets.


title: Constructive and Destructive Interference of Kerker-type Scattering in an Ultra-thin Silicon Huygens Metasurface

abstract: High refractive index dielectric nanoparticles have provided a new platform for exotic light manipulation through the interference of multipole modes. The Kerker effect is one example of a Huygens source design. Rather than exploiting interference between the electric dipole and magnetic dipole, as in many conventional Huygens source designs, we explore Kerker-type suppressed backward scattering mediated by the dominant electric dipole, toroidal dipole and magnetic quadrupole. These modes are provided by a designed and fabricated CMOS compatible ultra-thin Silicon nanodisk metasurface with a suppressed magnetic dipole contribution, and verified through multipole decomposition. The non-trivial substrate effect is considered using a semi-analytical transfer matrix model. The model successfully predicts the observed reflection dip. By applying a general criterion for constructive and destructive interference, it is shown that while constructive interference occurs between the electric and toroidal dipole contributions, the experimentally observed suppressed backward Kerker-type scattering arises from the destructive interference between backward scattered contributions due to the total electric dipole and the magnetic quadrupole. Our study paves the way towards new types of Huygens sources or metasurface design, such as for peculiar transverse Kerker scattering.


title: Fundamental Limits of Controlled Stochastic Dynamical Systems: An Information-Theoretic Approach

abstract: In this paper, we examine the fundamental performance limitations in the control of stochastic dynamical systems; more specifically, we derive generic $\mathcal{L}_p$ bounds that hold for any causal (stabilizing) controllers and any stochastic disturbances, by an information-theoretic analysis. We first consider the scenario where the plant (i.e., the dynamical system to be controlled) is linear time-invariant, and it is seen in general that the lower bounds are characterized by the unstable poles (or nonminimum-phase zeros) of the plant as well as the conditional entropy of the disturbance. We then analyze the setting where the plant is assumed to be (strictly) causal, for which case the lower bounds are determined by the conditional entropy of the disturbance. We also discuss the special cases of $p = 2$ and $p = \infty$, which correspond to minimum-variance control and controlling the maximum deviations, respectively. In addition, we investigate the power-spectral characterization of the lower bounds as well as its relation to the Kolmogorov-Szeg\"o formula.

categories: eess.SY cs.IT cs.LG cs.RO cs.SY math.IT math.OC


title: Latent Feature Representation via Unsupervised Learning for Pattern Discovery in Massive Electron Microscopy Image Volumes

abstract: We propose a method to facilitate exploration and analysis of new large data sets. In particular, we give an unsupervised deep learning approach to learning a latent representation that captures semantic similarity in the data set. The core idea is to use data augmentations that preserve semantic meaning to generate synthetic examples of elements whose feature representations should be close to one another. We demonstrate the utility of our method applied to nano-scale electron microscopy data, where even relatively small portions of animal brains can require terabytes of image data. Although supervised methods can be used to predict and identify known patterns of interest, the scale of the data makes it difficult to mine and analyze patterns that are not known a priori. We show the ability of our learned representation to enable query by example, so that if a scientist notices an interesting pattern in the data, they can be presented with other locations with matching patterns. We also demonstrate that clustering of data in the learned space correlates with biologically-meaningful distinctions. Finally, we introduce a visualization tool and software ecosystem to facilitate user-friendly interactive analysis and uncover interesting biological patterns. In short, our work opens possible new avenues in understanding of and discovery in large data sets, arising in domains such as EM analysis.


title: Statistically significant tests of multiparticle quantum correlations based on randomized measurements

abstract: We consider statistical methods based on finite samples of locally randomized measurements in order to certify different degrees of multiparticle entanglement in intermediate-scale quantum systems. We first introduce hierarchies of multi-qubit criteria, satisfied by states which are separable with respect to partitions of different size, involving only second moments of the underlying probability distribution. Then, we analyze in detail the statistical error of the estimation in experiments and present several approaches for estimating the statistical significance based on large deviation bounds. The latter allows us to characterize the measurement resources required for the certification of multiparticle correlations, as well as to analyze given experimental data in detail.


title: Quantum Convolutional Neural Networks for High Energy Physics Data Analysis

abstract: This work presents a quantum convolutional neural network (QCNN) for the classification of high energy physics events. The proposed model is tested using a simulated dataset from the Deep Underground Neutrino Experiment. The proposed architecture demonstrates the quantum advantage of learning faster than the classical convolutional neural networks (CNNs) under a similar number of parameters. In addition to faster convergence, the QCNN achieves greater test accuracy compared to CNNs. Based on experimental results, it is a promising direction to study the application of QCNN and other quantum machine learning models in high energy physics and additional scientific fields.

categories: cs.LG cs.AI hep-ex physics.data-an quant-ph


title: Reducing Solid-State Drive Read Latency by Optimizing Read-Retry (Extended Abstract)

abstract: 3D NAND flash memory with advanced multi-level cell techniques provides high storage density, but suffers from significant performance degradation due to a large number of read-retry operations. Although the read-retry mechanism is essential to ensuring the reliability of modern NAND flash memory, it can significantly increase the read latency of an SSD by introducing multiple retry steps that read the target page again with adjusted read-reference voltage values. Through a detailed analysis of the read mechanism and rigorous characterization of 160 real 3D NAND flash memory chips, we find new opportunities to reduce the read-retry latency by exploiting two advanced features widely adopted in modern NAND flash-based SSDs: 1) the CACHE READ command and 2) strong ECC engine. First, we can reduce the read-retry latency using the advanced CACHE READ command that allows a NAND flash chip to perform consecutive reads in a pipelined manner. Second, there exists a large ECC-capability margin in the final retry step that can be used for reducing the chip-level read latency. Based on our new findings, we develop two new techniques that effectively reduce the read-retry latency: 1) Pipelined Read-Retry (PR$^2$) and 2) Adaptive Read-Retry (AR$^2$). PR$^2$ reduces the latency of a read-retry operation by pipelining consecutive retry steps using the CACHE READ command. AR$^2$ shortens the latency of each retry step by dynamically reducing the chip-level read latency depending on the current operating conditions that determine the ECC-capability margin. Our evaluation using twelve real-world workloads shows that our proposal improves SSD response time by up to 31.5% (17% on average) over a state-of-the-art baseline with only small changes to the SSD controller.


title: Timely Monitoring of Dynamic Sources with Observations from Multiple Wireless Sensors

abstract: Age of Information (AoI) has recently received much attention due to its relevance in IoT sensing and monitoring applications. In this paper, we consider the problem of minimizing the AoI in a system in which a set of sources are observed by multiple sensors in a many-to-many relationship, and the probability that a sensor observes a source depends on the state of the source. This model represents many practical scenarios, such as the ones in which multiple cameras or microphones are deployed to monitor objects moving in certain areas. We formulate the scheduling problem as a Markov Decision Process, and show how the age-optimal scheduling policy can be obtained. We further consider partially observable variants of the problem, and devise approximate policies for large state spaces. Our evaluations show that the approximate policies work well in the considered scenarios, and that the fact that sensors can observe multiple sources is beneficial, especially when there is high uncertainty of the source states.


title: Cloud removal in remote sensing images using generative adversarial networks and SAR-to-optical image translation

abstract: Satellite images are often contaminated by clouds. Cloud removal has received much attention due to the wide range of satellite image applications. As the clouds thicken, the process of removing the clouds becomes more challenging. In such cases, using auxiliary images such as near-infrared or synthetic aperture radar (SAR) for reconstructing is common. In this study, we attempt to solve the problem using two generative adversarial networks (GANs). The first translates SAR images into optical images, and the second removes clouds using the translated images of prior GAN. Also, we propose dilated residual inception blocks (DRIBs) instead of vanilla U-net in the generator networks and use structural similarity index measure (SSIM) in addition to the L1 Loss function. Reducing the number of downsamplings and expanding receptive fields by dilated convolutions increase the quality of output images. We used the SEN1-2 dataset to train and test both GANs, and we made cloudy images by adding synthetic clouds to optical images. The restored images are evaluated with PSNR and SSIM. We compare the proposed method with state-of-the-art deep learning models and achieve more accurate results in both SAR-to-optical translation and cloud removal parts.


title: Designing an Interactive Visualization System for Monitoring Participant Compliance in a Large-Scale, Longitudinal Study

abstract: Frequent monitoring of participant compliance is necessary when conducting large-scale, longitudinal studies to ensure that the collected data is of sufficiently high quality. While the need for achieving high compliance has been underscored and there are discussions on incentives and factors affecting compliance, little is shared about the actual processes and tools used for monitoring compliance in such studies. Monitoring participant compliance with respect to multi-modal data can be a tedious process, especially if there are only a few personnel involved. In this case study, we describe the iterative design of an interactive visualization system we developed for monitoring compliance and refined based on changing requirements in an ongoing study. We find that the visualization system, leveraging the digital medium, both facilitates the exploratory tasks of monitoring participant compliance and supports asynchronous collaboration among non-co-located researchers. Our documented requirements for checking participant compliance as well as the design of the visualization system can help inform the compliance-monitoring process in future studies.


title: Reduced Complexity Simulation of Wireless Sensor Networks for Application Development

abstract: This paper presents an approach for low-cost simulation modeling for application development for wireless sensor networks. Computational complexity of simulating wireless sensor networks can be very high and as such must be carefully managed. Application-level code prototyping with reasonable accuracy and fidelity can be accomplished through simulation that models only the effects of the wireless and distributed computations which materialize mainly as delay and drop for the messages being exchanged among the motes. This approach employs the abstraction that all physical or communication and protocol level operations can be represented in terms of their effects as message delay and drop at the application level for a wireless sensor network. This study proposes that idea of empirical modeling of delay and drop and employing those models to affect the reception times of wirelessly communicated messages. It further proposes the delay and drop to be modeled as random variables with probability distributions empirically approximated based on the data reported in the literature. The proposed approach is demonstrated through development of a neural network application with neurons distributed across the motes of a wireless sensor network. Delay and drop are incorporated into wireless communications, which carry neuron output values among motes. A set of classification data sets from the Machine Learning Repository are employed to demonstrate the performance of the proposed system in a comparative context with the similar studies in the literature. Results and findings indicate that the proposed approach of abstracting wireless sensor network operation in terms of message delay and drop at the application level is feasible to facilitate development of applications with competitive performance profiles while minimizing the spatio-temporal cost of simulation.


title: A PMU-Based Machine Learning Application for Fast Detection of Forced Oscillations from Wind Farms

abstract: Today's evolving power system contains an increasing amount of power electronic interfaced energy sources and loads that require a paradigm shift in utility operations. Sub-synchronous oscillations at frequencies around 13-15 Hz, for instance, have been reported by utilities due to wind farm controller interactions with the grid. Dynamics at such frequencies are unobservable by most SCADA tools due to low sampling frequencies and lack of synchronization. Real-time or off-line frequency domain analysis of phasor measurement unit (PMU) data has become a valuable method to identify such phenomena, at the expense of costly power system data and communication infrastructure. This article proposes an alternative machine learning (ML) based application for sub-synchronous oscillation detection in wind farm applications. The application is targeted for real-time implementation at the edge, resulting in significant savings in terms of data and communication requirements. Validation is performed using data from a North American wind farm operator.


title: COVID-19 Emotion Monitoring as a Tool to Increase Preparedness for Disease Outbreaks in Developing Regions

abstract: The COVID-19 pandemic brought many challenges, from hospital-occupation management to lock-down mental-health repercussions such as anxiety or depression. In this work, we present a solution for the later problem by developing a Twitter emotion-monitor system based on a state-of-the-art natural-language processing model. The system monitors six different emotions on accounts in cities, as well as politicians and health-authorities Twitter accounts. With an anonymous use of the emotion monitor, health authorities and private health-insurance companies can develop strategies to tackle problems such as suicide and clinical depression. The model chosen for such a task is a Bidirectional-Encoder Representations from Transformers (BERT) pre-trained on a Spanish corpus (BETO). The model performed well on a validation dataset. The system is deployed online as part of a web application for simulation and data analysis of COVID-19, in Colombia, available at https://epidemiologia-matematica.org.


title: A Shell Bonded to an Elastic Foundation and the Existence of Optimal Elastic and Geometric Conditions for the Membrane Case

abstract: In this article, we derive a mathematical model for a shell (i.e. a thin elastic body) bonded to an elastic foundation by modifying Koiter's linear shell equations. We prove the existence and the uniqueness of solutions, and we explicitly derive the governing equations and the boundary conditions for the general case. Finally, with numerical modelling and asymptotic analysis, we show that there exist optimal values for the Young's modulus, the Poisson's ratio and the thickness of the shell (relative to the elastic foundation), and the radius of curvature of the contact region such that the planar solution derived by the shell model (i.e. the membrane case, where stretching effects are dominant) results in a good approximation of the thin body. It is often regarded in the field of stretchable and flexible electronics that the planar solution is mostly accurate when the thin body (i.e. the shell or the membrane) is relatively stiffer (i.e. has a high Young's modulus) than the thicker foundation. As far as we are aware, this is the first analysis showing that for a membrane (or the planar solution of a shell) bonded to an elastic foundation, higher stiffness (relative to the elastic foundation's stiffness) alone would not guarantee a more accurate solution.


title: Learning to Play Imperfect-Information Games by Imitating an Oracle Planner

abstract: We consider learning to play multiplayer imperfect-information games with simultaneous moves and large state-action spaces. Previous attempts to tackle such challenging games have largely focused on model-free learning methods, often requiring hundreds of years of experience to produce competitive agents. Our approach is based on model-based planning. We tackle the problem of partial observability by first building an (oracle) planner that has access to the full state of the environment and then distilling the knowledge of the oracle to a (follower) agent which is trained to play the imperfect-information game by imitating the oracle's choices. We experimentally show that planning with naive Monte Carlo tree search does not perform very well in large combinatorial action spaces. We therefore propose planning with a fixed-depth tree search and decoupled Thompson sampling for action selection. We show that the planner is able to discover efficient playing strategies in the games of Clash Royale and Pommerman and the follower policy successfully learns to implement them by training on a few hundred battles.


title: Near-Field Angular Scan Enhancement of Antenna Arrays Using Metasurfaces

abstract: In this publication we continue our previous work on extending the angular scan range of phased arrays using metasurfaces. We consider in detail scan enhancement using a single metasurface lens in the near-field of a source and provide analytical expressions for source excitation to obtain a desired beam. We show that such a device suffers from the same directivity degradation as in the case of far-field lens placement. Moreover, we discuss in detail that this approach has limitations, and one cannot place the lens arbitrarily close to the source array. We then propose a two-lens near-field scan enhancer which is not subject to the same limitations, although still subject to the same directivity degradation. Finally we propose a binary metasurface concept which theoretically achieves desired scan enhancement without directivity degradation. Some theoretical claims are then verified via full-field simulations.


title: Automated Multi-Channel Segmentation for the 4D Myocardial Velocity Mapping Cardiac MR

abstract: Four-dimensional (4D) left ventricular myocardial velocity mapping (MVM) is a cardiac magnetic resonance (CMR) technique that allows assessment of cardiac motion in three orthogonal directions. Accurate and reproducible delineation of the myocardium is crucial for accurate analysis of peak systolic and diastolic myocardial velocities. In addition to the conventionally available magnitude CMR data, 4D MVM also acquires three velocity-encoded phase datasets which are used to generate velocity maps. These can be used to facilitate and improve myocardial delineation. Based on the success of deep learning in medical image processing, we propose a novel automated framework that improves the standard U-Net based methods on these CMR multi-channel data (magnitude and phase) by cross-channel fusion with attention module and shape information based post-processing to achieve accurate delineation of both epicardium and endocardium contours. To evaluate the results, we employ the widely used Dice scores and the quantification of myocardial longitudinal peak velocities. Our proposed network trained with multi-channel data shows enhanced performance compared to standard U-Net based networks trained with single-channel data. Based on the results, our method provides compelling evidence for the design and application for the multi-channel image analysis of the 4D MVM CMR data.


title: Detecting Planetary-mass Primordial Black Holes with Resonant Electromagnetic Gravitational Wave Detectors

abstract: The possibility to detect gravitational waves (GW) from planetary-mass primordial black hole (PBH) binaries with electromagnetic (EM) detectors of high-frequency GWs is investigated. We consider two patented experimental designs, based on the inverse Gertsenshtein effect, in which incoming GWs passing through a static magnetic field induce EM excitations inside either a TM cavity or a TEM waveguide. The frequency response of the detectors is computed for post-newtonian GW waveforms. We find that such EM detectors based on current technology may achieve a strain sensitivity down to $h \sim 10^{-30}$, which generates an EM power variation of $10^{-10}$ W. This allows the detection of PBH binary mergers of mass around $10^{-5} M_\odot$ if they constitute more than $0.01$ percent of the dark matter, as suggested by recent microlensing observations. We envision that this class of detectors could also be used to detect cosmological GW backgrounds and probe sources in the early Universe at energies up to the GUT scale.


title: Identification of Additive Link Metrics: Proof of Selected Theorems

abstract: This is a technical report, containing all the theorem proofs in the following two papers: (1) Liang Ma, Ting He, Kin K. Leung, Ananthram Swami, and Don Towsley, "Identifiability of Link Metrics Based on End-to-end Path Measurements," in ACM IMC, 2013. (2) Liang Ma, Ting He, Kin K. Leung, Ananthram Swami, and Don Towsley, "Inferring Link Metrics from End-to-end Path Measurements: Identifiability and Monitor Placement," IEEE/ACM Transactions on Networking, vol. 22, no. 4, pp. 1351-1368, 2014.


title: Efficient Identification of Additive Link Metrics: Theorem Proof and Evaluations

abstract: This is a technical report, containing all the theorem proofs and additional evaluations in paper "Efficient Identification of Additive Link Metrics via Network Tomography" by Liang Ma, Ting He, Kin K. Leung, Don Towsley, and Ananthram Swami, published in IEEE ICDCS, 2013.


title: Query Answering via Decentralized Search

abstract: Expert networks are formed by a group of expert-professionals with different specialties to collaboratively resolve specific queries posted to the network. In such networks, when a query reaches an expert who does not have sufficient expertise, this query needs to be routed to other experts for further processing until it is completely solved; therefore, query answering efficiency is sensitive to the underlying query routing mechanism being used. Among all possible query routing mechanisms, decentralized search, operating purely on each expert's local information without any knowledge of network global structure, represents the most basic and scalable routing mechanism, which is applicable to any network scenarios even in dynamic networks. However, there is still a lack of fundamental understanding of the efficiency of decentralized search in expert networks. In this regard, we investigate decentralized search by quantifying its performance under a variety of network settings. Our key findings reveal the existence of network conditions, under which decentralized search can achieve significantly short query routing paths (i.e., between $O(\log n)$ and $O(\log^2 n)$ hops, $n$: total number of experts in the network). Based on such theoretical foundation, we further study how the unique properties of decentralized search in expert networks is related to the anecdotal small-world phenomenon. In addition, we demonstrate that decentralized search is robust against estimation errors introduced by misinterpreting the required expertise levels. To the best of our knowledge, this is the first work studying fundamental behaviors of decentralized search in expert networks. The developed performance bounds, confirmed by real datasets, are able to assist in predicting network performance and designing complex expert networks.


title: Bridgeland Moduli spaces for Gushel-Mukai threefolds and Kuznetsov's Fano threefold conjecture

abstract: We study the Hilbert scheme $\mathcal{H}$ of twisted cubics on a special smooth Gushel-Mukai threefolds $X_{10}$. We show that it is a smooth irreducible projective threefold if $X_{10}$ is general among special Gushel-Mukai threefolds, while it is singular if $X_{10}$ is not general. We construct an irreducible component of a moduli space of Bridgeland stable objects in the Kuznetsov component of $X_{10}$ as a divisorial contraction of $\mathcal{H}$. We also identify the minimal model of Fano surface $\mathcal{C}(X_{10}')$ of conics on a smooth ordinary Gushel-Mukai threefold with moduli space of Bridgeland stable objects in the Kuznetsov component of $X_{10}'$. As a result, we show that the Kuznetsov's Fano threefold conjecture is not true


title: Derivation of an extended Bjerrum equation for the activity coefficient of ions based on an analysis of Coulombic forces

abstract: The activity coefficient of ions in solution was proposed by Bjerrum (1916) to depend on the cube root of concentration, because of a good fit with data for low and moderate salt concentration. However, Debye and H\"uckel (DH) later developed a theory that prevailed, and this theory has a square root dependence. The derivation of the DH equation is not easy or intuitive, and has various uncertain elements, and the fit to data is not very good for salts with non-unity valencies. We develop a model for the activity coefficient based on the Coulombic forces between an ion and the nearest ions of opposite charge, including the distribution of separations between these ions. For a symmetric salt, we only have to analyse the interactions between one anion and one cation, and we can derive an expression for the dilute limit that depends on the Bjerrum length and the cube root of salt concentration, the same as an expression put forward by Bjerrum. Our theory also analytically describes higher salt concentrations and a non-zero ion size, not only for 1:1 but also for 2:2 and 3:3 salts. For a 1:1 salt, the cube root scaling law describes data of the activity coefficient very well, both at low and intermediate salt concentrations, while the extended Bjerrum equation describes data even better and up to higher concentrations, and also includes the effect of ion size. For asymmetric salts (2:1 and 3:1 salts), a numerical procedure can be used which considers all possible orbits of two or three monovalent ions around a central multivalent ion of opposite charge, which also fits data accurately.


title: Labels Are Not Perfect: Inferring Spatial Uncertainty in Object Detection

abstract: The availability of many real-world driving datasets is a key reason behind the recent progress of object detection algorithms in autonomous driving. However, there exist ambiguity or even failures in object labels due to error-prone annotation process or sensor observation noise. Current public object detection datasets only provide deterministic object labels without considering their inherent uncertainty, as does the common training process or evaluation metrics for object detectors. As a result, an in-depth evaluation among different object detection methods remains challenging, and the training process of object detectors is sub-optimal, especially in probabilistic object detection. In this work, we infer the uncertainty in bounding box labels from LiDAR point clouds based on a generative model, and define a new representation of the probabilistic bounding box through a spatial uncertainty distribution. Comprehensive experiments show that the proposed model reflects complex environmental noises in LiDAR perception and the label quality. Furthermore, we propose Jaccard IoU (JIoU) as a new evaluation metric that extends IoU by incorporating label uncertainty. We conduct an in-depth comparison among several LiDAR-based object detectors using the JIoU metric. Finally, we incorporate the proposed label uncertainty in a loss function to train a probabilistic object detector and to improve its detection accuracy. We verify our proposed methods on two public datasets (KITTI, Waymo), as well as on simulation data. Code is released at https://bit.ly/2W534yo.


title: Identifiability of Bifactor Models

abstract: The bifactor model and its extensions are multidimensional latent variable models, under which each item measures up to one subdimension on top of the primary dimension(s). Despite their wide applications to educational and psychological assessments, this type of multidimensional latent variable models may suffer from non-identifiability, which can further lead to inconsistent parameter estimation and invalid inference. The current work provides a relatively complete characterization of identifiability for the linear and dichotomous bifactor models and the linear extended bifactor model with correlated subdimensions. In addition, similar results for the two-tier models are also developed. Illustrative examples are provided on checking model identifiability through inspecting the factor loading structure. Simulation studies are reported that examine estimation consistency when the identifiability conditions are/are not satisfied.


title: On the spread of infinite groups

abstract: A group is $\frac32$-generated if every non-trivial element is part of a generating pair. In 2019, Donoven and Harper showed that many Thompson groups are $\frac32$-generated and posed five questions. The first of these is whether there exists a 2-generated group with every proper quotient cyclic that is not $\frac32$-generated. This is a natural question given the significant work in proving that no finite group has this property, but we show that there is such an infinite group. The groups we consider are a family of finite index subgroups of the Houghton group $\text{FSym}(\mathbb{Z})\rtimes\mathbb{Z}$. We then show that the first two groups in our family are $\frac32$-generated, and investigate the related notion of spread for these groups. We are able to show that they have finite spread which is greater than 2. These are therefore the first infinite groups to be shown to have finite positive spread, and the first to be shown to have spread greater than 2 (other than $\mathbb{Z}$ and the Tarski monsters, which have infinite spread).


title: Forensic strength of evidence statements should preferably be likelihood ratios calculated using relevant data, quantitative measurements, and statistical models

abstract: Lennard (2013) [Fingerprint identification: how far have we come? Aus J Forensic Sci. doi:10.1080/00450618.2012.752037] proposes that the numeric output of statistical models should not be presented in court (except "if necessary" / "if required"). Instead he argues in favour of an "expert opinion" which may be informed by a statistical model but which is not itself the output of a statistical model. We argue that his proposed procedure lacks the transparency, the ease of testing of validity and reliability, and the relative robustness to cognitive bias that are the strengths of a likelihood-ratio approach based on relevant data, quantitative measurements, and statistical models, and that the latter is therefore preferable.


title: Involuntary unemployment in overlapping generations model due to instability of the economy

abstract: The existence of involuntary unemployment advocated by J. M. Keynes is a very important problem of the modern economic theory. Using a three-generations overlapping generations model, we show that the existence of involuntary unemployment is due to the instability of the economy. Instability of the economy is the instability of the difference equation about the equilibrium price around the full-employment equilibrium, which means that a fall in the nominal wage rate caused by the presence of involuntary unemployment further reduces employment. This instability is due to the negative real balance effect that occurs when consumers' net savings (the difference between savings and pensions) are smaller than their debt multiplied by the marginal propensity to consume from childhood consumption.


title: Robustness of Higher Dimensional Nonlocality against dual noise and sequential measurements

abstract: Robustness in the violation of Collins-Linden-Gisin-Masser-Popescu (CGLMP) inequality is investigated from the dual perspective of noise in measurements as well as in states. To quantify it, we introduce a quantity called area of nonlocal region which reveals a dimensional advantage. Specifically, we report that with the increase of dimension, the maximally violating states (MVS) show a greater enhancement in the area of nonlocal region in comparison to the maximally entangled states (MES) and the scaling of the increment, in this case, grows faster than visibility. Moreover, we examine the robustness in the sequential violation of CGLMP inequality using weak measurements, and find that even for higher dimensions, two observers showing a simultaneous violation of the CGLMP inequality as obtained for two-qubit states persists. We notice that the complementarity between information gain and disturbance by measurements is manifested by the decrease of the visibility in the first round and the increase of the same in the second round with dimensions. Furthermore, the amount of white noise that can be added to an MES so that it gives two rounds of the violation, decreases with the dimension, while the same does not appreciably change for the MVS.


title: The GALAH+ Survey: A New Library of Observed Stellar Spectra Improves Radial Velocities and Hints at Motions within M67

abstract: GALAH+ is a magnitude-limited survey of high resolution stellar spectra obtained by the HERMES spectrograph at the Australian Astronomical Observatory. Its third data release provides reduced spectra with new derivations of stellar parameters and abundances of 30 chemical elements for 584,015 dwarfs and giants, 88% of them in the Gaia magnitude range 11 < G < 14. Here we use these improved values of stellar parameters to build a library of observed spectra which is useful to study variations of individual spectral lines with stellar parameters. This and other improvements are used to derive radial velocities with uncertainties which are generally within 0.1 km/s or ~25% smaller than in the previous release. Median differences in radial velocities measured here and by the Gaia DR2 or APOGEE DR16 surveys are smaller than 30 m/s, a larger offset is present only for Gaia measurements of giant stars. We identify 4483 stars with intrinsically variable velocities and 225 stars for which the velocity stays constant over >=3 visits spanning more than a year. The combination of radial velocities from GALAH+ with distances and sky plane motions from Gaia enables studies of dynamics within streams and clusters. For example, we estimate that the open cluster M67 has a total mass of ~3300 Msun and its outer parts seem to be expanding, though astrometry with a larger time-span than currently available from Gaia eDR3 is needed to judge if the latter result is real.


title: Machine Learning and cosmographic reconstructions of quintessence and the Swampland conjectures

abstract: We present model independent reconstructions of quintessence and the Swampland conjectures (SC) using both Machine Learning (ML) and cosmography. In particular, we demonstrate how the synergies between theoretical analyses and ML can provide key insights on the nature of dark energy and modified gravity. Using the Hubble parameter $H(z)$ data from the cosmic chronometers we find that the ML and cosmography reconstructions of the SC are compatible with observations at low redshifts. Finally, including the growth rate data $f\sigma_8(z)$ we perform a model independent test of modified gravity cosmologies through two phase diagrams, namely $H-f\sigma_8$ and $\eta-f\sigma_8$, where the anisotropic stress parameter $\eta$ is obtained via the $E_g$ statistics, which is related to gravitational lensing data. While the first diagram is consistent within the errors with the $\Lambda$CDM model, the second one has a $\sim 2\sigma$ deviation of the anisotropic stress from unity at $z\sim 0.3$ and a $\sim 4\sigma$ deviation at $z\sim 0.9$, thus pointing toward mild deviations from General Relativity, which could be further tested with upcoming large-scale structure surveys.


title: Autonomous sPOMDP Environment Modeling With Partial Model Exploitation

abstract: A state space representation of an environment is a classic and yet powerful tool used by many autonomous robotic systems for efficient and often optimal solution planning. However, designing these representations with high performance is laborious and costly, necessitating an effective and versatile tool for autonomous generation of state spaces for autonomous robots. We present a novel state space exploration algorithm by extending the original surprise-based partially-observable Markov Decision Processes (sPOMDP), and demonstrate its effective long-term exploration planning performance in various environments. Through extensive simulation experiments, we show the proposed model significantly increases efficiency and scalability of the original sPOMDP learning techniques with a range of 31-63% gain in training speed while improving robustness in environments with less deterministic transitions. Our results pave the way for extending sPOMDP solutions to a broader set of environments.


title: Velocity-dependent interacting dark energy and dark matter with a Lagrangian description of perfect fluids

abstract: We consider a cosmological scenario where the dark sector is described by two perfect fluids that interact through a velocity-dependent coupling. This coupling gives rise to an interaction in the dark sector driven by the relative velocity of the components, thus making the background evolution oblivious to the interaction and only the perturbed Euler equations are affected at first order. We obtain the equations governing this system with the Schutz-Sorkin Lagrangian formulation for perfect fluids and derive the corresponding stability conditions to avoid ghosts and Laplacian instabilities. As a particular example, we study a model where dark energy behaves as a radiation fluid at high redshift while it effectively becomes a cosmological constant in the late Universe. Within this scenario, we show that the interaction of both dark components leads to a suppression of the dark matter clustering at late times. We also argue the possibility that this suppression of clustering together with the additional dark radiation at early times can simultaneously alleviate the $\sigma_8$ and $H_0$ tensions.


title: NDRIO White Paper: Envisioning Digital Research Infrastructure for the Simons Observatory

abstract: Observations of the cosmic microwave background (CMB) are an incredibly fertile source of information for studying the origins and evolution of the Universe. Canadian digital research infrastructure (DRI) has played a key role in reducing ever-larger quantities of raw data into maps of the CMB suitable for scientific analysis, as exemplified by the many scientific results produced by the Atacama Cosmology Telescope (ACT) over the past decade. The Simons Observatory (SO), due to start observing in 2023, will be able to measure the CMB with about an order of magnitude more sensitivity than ACT and other current telescopes. In this White Paper we outline how Canadian DRI under the New Digital Research Infrastructure Organization (NDRIO) could build upon the legacy of ACT and play a pivotal role in processing SO data, helping to produce data products that will be central to the cosmology community for years to come. We include estimates of DRI resources required for this work to indicate what kind of advanced research computing (ARC) would best support an SO-like project. Finally, we comment on how ARC allocations could be structured for large collaborations like SO and propose a research data management (RDM) system that makes public data releases available not only for download but also for direct analysis on Canadian DRI.


title: FracBNN: Accurate and FPGA-Efficient Binary Neural Networks with Fractional Activations

abstract: Binary neural networks (BNNs) have 1-bit weights and activations. Such networks are well suited for FPGAs, as their dominant computations are bitwise arithmetic and the memory requirement is also significantly reduced. However, compared to start-of-the-art compact convolutional neural network (CNN) models, BNNs tend to produce a much lower accuracy on realistic datasets such as ImageNet. In addition, the input layer of BNNs has gradually become a major compute bottleneck, because it is conventionally excluded from binarization to avoid a large accuracy loss. This work proposes FracBNN, which exploits fractional activations to substantially improve the accuracy of BNNs. Specifically, our approach employs a dual-precision activation scheme to compute features with up to two bits, using an additional sparse binary convolution. We further binarize the input layer using a novel thermometer encoding. Overall, FracBNN preserves the key benefits of conventional BNNs, where all convolutional layers are computed in pure binary MAC operations (BMACs). We design an efficient FPGA-based accelerator for our novel BNN model that supports the fractional activations. To evaluate the performance of FracBNN under a resource-constrained scenario, we implement the entire optimized network architecture on an embedded FPGA (Xilinx Ultra96v2). Our experiments on ImageNet show that FracBNN achieves an accuracy comparable to MobileNetV2, surpassing the best-known BNN design on FPGAs with an increase of 28.9% in top-1 accuracy and a 2.5x reduction in model size. FracBNN also outperforms a recently introduced BNN model with an increase of 2.4% in top-1 accuracy while using the same model size. On the embedded FPGA device, FracBNN demonstrates the ability of real-time image classification.


title: Model independent analysis of the angular observables in $B^{0} \to K^{*0} \mu^+ \mu^-$ and $B^{+} \to K^{*+} \mu^+ \mu^-$

abstract: We analyse the results recently presented on the $B^{+} \to K^{*+} \mu^+ \mu^-$ angular observables by the LHCb Collaboration which show indications for New Physics beyond the Standard Model. Within a model-independent analysis, we compare the fit results with the corresponding results for the angular observables in $B^{0} \to K^{*0} \mu^+ \mu^-$.