title: Multivariate binary probability distribution in the Grassmann formalism

abstract: We propose a probability distribution for multivariate binary random variables. For this purpose, we use the Grassmann number, an anti-commuting number. In our model, the partition function, the central moment, and the marginal and conditional distributions are expressed analytically by the matrix of the parameters analogous to the covariance matrix in the multivariate Gaussian distribution. That is, summation over all possible states is not necessary for obtaining the partition function and various expected values, which is a problem with the conventional multivariate Bernoulli distribution. The proposed model has many similarities to the multivariate Gaussian distribution. For example, the marginal and conditional distributions are expressed by the parameter matrix and its inverse matrix, respectively. That is, the inverse matrix expresses a sort of partial correlation. Analytical expressions for the marginal and conditional distributions are also useful in generating random numbers for multivariate binary variables. Hence, we validated the proposed method using synthetic datasets. We observed that the sampling distributions of various statistics are consistent with the theoretical predictions and estimates are consistent and asymptotically normal.


title: Using angular momentum maps to detect kinematically distinct galactic components

abstract: In this work we introduce a physically motivated method of performing disc/spheroid decomposition of simulated galaxies, which we apply to the Eagle sample. We make use of the HEALPix package to create Mollweide projections of the angular momentum map of each galaxy's stellar particles. A number of features arise on the angular momentum space which allows us to decompose galaxies and classify them into different morphological types. We assign stellar particles with angular separation of less/greater than 30 degrees from the densest grid cell on the angular momentum sphere to the disc/spheroid components, respectively. We analyse the spatial distribution for a subsample of galaxies and show that the surface density profiles of the disc and spheroid closely follow an exponential and a Sersic profile, respectively. In addition discs rotate faster, have smaller velocity dispersions, are younger and are more metal rich than spheroids. Thus our morphological classification reproduces the observed properties of such systems. Finally, we demonstrate that our method is able to identify a significant population of galaxies with counter-rotating discs and provide a more realistic classification of such systems compared to previous methods.


title: On the influence of a Chern-Simons term in the dynamical fermion masses in Reduced or Pseudo QED

abstract: Mixed dimensional theories have been used to describe condensed matter systems where fermions are constrained to a plane while the gauge fields they interact with remain four-dimensional. Here we investigate dynamical breaking of chiral symmetry in the framework of a mixed-dimensional theory that has been successful in describing graphene and other planar Dirac materials and has been dubbed as Pseudo or Reduced Quantum Electrodynamics. We explore the interplay between the gauge and fermion sector when a Chern-Simons term is considered and, by exploring the tensor structure of the fermion self energy, we show that the radiative corrections may induce both, a Dirac and a Haldane mass term. Furthermore, by solving the corresponding Schwinger-Dyson equation in a suitable truncation, we explore non-perturbatively the dynamical generation of fermion masses for different values of the electromagnetic coupling and the Chern-Simons coefficient. We also show that for definite parameter values, the contributions from each chirality cancel out and the chiral symmetry may be restored. Possible implications of this result for physical systems, in particular on what concerns the chiral magnetic effect, are discussed.


title: Congruences on K-theoretic Gromov--Witten invariants

abstract: We study K-theoretic Gromov--Witten invariants of projective hypersurfaces using a virtual localization formula under finite group actions. In particular, it provides all K-theoretic Gromov--Witten invariants of the quintic threefold modulo 41, up to genus 19 and degree 40. As an illustration, we give an instance in genus one and degree one. Applying the same idea to a K-theoretic version of FJRW theory, we determine it modulo 41 for the quintic polynomial with minimal group and narrow insertions, in every genus.


title: An improved existence criterion and an optimal result

abstract: We are concerned with a semi-linear elliptic equation on a smooth bounded domain $\Omega$ of $\mathbb{R}^n,\,n\geq 5,$ which involves a critical nonlinearity and a linear term of the form $K(x)u^{(n+2)/(n-2)}$ and $\mu u,$ respectively. By using a test function procedure, we give an existence criterion involving the parameter $\mu$ and the function $K(x).$ For a particular case of $\Omega,\,K(x)$ and $n,$ we prove its optimality through a Pohozaev type identity.


title: Symmetry mediated tunable molecular magnetism on a 2D material

abstract: The induction of unconventional superconductivity by twisting two layers of graphene a small angle was groundbreaking1, and since then has attracted widespread attention to novel phenomena caused by lattice or angle mismatch between two-dimensional (2D) materials2. While many studies address the influence of angle mismatch between layered 2D materials3-5 , the impact of the absorption alignment on the physical properties of planar molecules on 2D substrates has not been studied in detail. Using scanning probe microscopy (SPM) we show that individual cobalt phthalocyanine (CoPc) molecules adsorbed on the layered superconductor 2H-NbSe2 change drastically their charge and spin state when the symmetry axes of the molecule and the substrate are twisted with respect to each other. The CoPc changes from an effective spin-1/2 as found in gas-phase6 to a molecule with non-magnetic ground-state. On the latter we observe a singlet-triplet transition originating from an antiferromagnetic interaction between the central-ion spin and a distributed magnetic moment on the molecular ligands. Because the Ising superconductor 2H-NbSe2 lacks inversion symmetry and has large spin-orbit coupling7 this intramolecular magnetic exchange has significant non-collinear Dzyaloshinskii-Moriya (DM)8, 9 contribution.


title: $C^*$-correspondence functoriality of Cuntz-Pimsner algebras

abstract: We construct a functor that maps $C^*$-correspondences to their Cuntz-Pimsner algebras. The objects in our domain category are $C^*$-correspondences, and the morphisms are the isomorphism classes of $C^*$-correspondences satisfying certain conditions. As an application, we recover a well-known result of Muhly and Solel. In fact, we show that functoriality leads us to a more generalized result: strongly Morita equivalent $C^*$-correspondences have Morita equivalent Cuntz-Pimsner algebras.


title: Quantum Probability's Algebraic Origin

abstract: Max Born's statistical interpretation made probabilities play a major role in quantum theory. Here we show that these quantum probabilities and the classical probabilities have very different origins. While the latter always result from an assumed probability measure, the first include transition probabilities with a purely algebraic origin. Moreover, the general definition of transition probability introduced here comprises not only the well-known quantum mechanical transition probabilities between pure states or wave functions, but further physically meaningful and experimentally verifiable novel cases. A transition probability that differs from 0 and 1 manifests the typical quantum indeterminacy in a similar way as Heisenberg's and others' uncertainty relations and, furthermore, rules out deterministic states in the same way as the Bell-Kochen-Specker theorem. However, the transition probability defined here achieves a lot more beyond that: it demonstrates that the algebraic structure of the Hilbert space quantum logic dictates the precise values of certain probabilities and it provides an unexpected access to these quantum probabilities that does not rely on states or wave functions.


title: Spatio-Temporal Probabilistic Voltage Sensitivity Analysis - A Novel Framework for Hosting Capacity Analysis

abstract: Smart grids are envisioned to accommodate high penetration of distributed photovoltaic (PV) generation, which may cause adverse grid impacts in terms of voltage violations. Therefore, PV Hosting capacity (HC) is being used as a planning tool to determine the maximum PV installation capacity that causes the first voltage violation and above which would require infrastructure upgrades. Traditional methods of HC analysis are computationally complex as they are based on iterative load flow algorithms that require investigation of a large number of scenarios for accurate assessment of PV impacts. This paper first presents a computationally efficient analytical approach to compute the probability distribution of voltage change at a particular node due to random behavior of randomly located multiple distributed PVs. Next, the derived distribution is used to identify voltage violations for various PV penetration levels and subsequently determine the HC of the system without the need to examine multiple scenarios. Results from the proposed spatio-temporal probabilistic voltage sensitivity analysis and the HC are validated via conventional load flow based simulation approach on the IEEE 37 and IEEE 123 node test systems.


title: Globular structure of the hypermineralized tissue in human femoral neck

abstract: Bone becomes more fragile with ageing. Among many structural changes, a thin layer of highly mineralized and brittle tissue covers part of the external surface of the thin femoral neck cortex in older people and has been proposed to increase hip fragility. However, there have been very limited reports on this hypermineralized tissue in the femoral neck, especially on its ultrastructure. Such information is critical to understanding both the mineralization process and its contributions to hip fracture. Here, we use multiple advanced techniques to characterize the ultrastructure of the hypermineralized tissue in the neck across various length scales. Synchrotron radiation micro-CT found larger but less densely distributed cellular lacunae in hypermineralized tissue than in lamellar bone. When examined under FIB-SEM, the hypermineralized tissue was mainly composed of mineral globules with sizes varying from submicron to a few microns. Nano-sized channels were present within the mineral globules and oriented with the surrounding organic matrix. Transmission electron microscopy showed the apatite inside globules were poorly crystalline, while those at the boundaries between the globules had well-defined lattice structure with crystallinity similar to the apatite mineral in lamellar bone. No preferred mineral orientation was observed both inside each globule and at the boundaries. Collectively, we conclude based on these new observations that the hypermineralized tissue is non-lamellar and has less organized mineral, which may contribute to the high brittleness of the tissue.


title: Designing knowledge plane to optimize leaf and spine data center

abstract: In the last few decades, data center architecture evolved from the traditional client-server to access-aggregation-core architectures. Recently there is a new shift in the data center architecture due to the increasing need for low latency and high throughput between server-to-server communications, load balancing and, loop-free environment. This new architecture, known as leaf and spine architecture, provides low latency and minimum packet loss by enabling the addition and deletion of network nodes on demand. Network nodes can be added or deleted from the network based on network statistics like link speed, packet loss, latency, and throughput. With the maturity of Open Virtual Switch (OvS) and OpenFlow based Software Defined Network (SDN) controllers, network automation through programmatic extensions has become possible based on network statistics. The separation of the control plane and data plane has enabled automated management of network and Machine Learning (ML) can be applied to learn and optimize the network. In this publication, we propose the design of an ML-based approach to gather network statistics and build a knowledge plane. We demonstrate that this knowledge plane enables data center optimization using southbound APIs and SDN controllers. We describe the design components of this approach - using a network simulator and show that it can maintain the historical patterns of network statistics to predict future growth or decline. We also provide an open-source software that can be utilized in a leaf and spine data center to provide elastic capacity based on load forecasts.


title: Finite elements for Helmholtz equations with a nonlocal boundary condition

abstract: Numerical resolution of exterior Helmholtz problems requires some approach to domain truncation. As an alternative to approximate nonreflecting boundary conditions and invocation of the Dirichlet-to-Neumann map, we introduce a new, nonlocal boundary condition. This condition is exact and requires the evaluation of layer potentials involving the free space Green's function. However, it seems to work in general unstructured geometry, and Galerkin finite element discretization leads to convergence under the usual mesh constraints imposed by G{\aa}rding-type inequalities. The nonlocal boundary conditions are readily approximated by fast multipole methods, and the resulting linear system can be preconditioned by the purely local operator involving transmission boundary conditions.


title: Practical Link Adaptation Algorithm with Power Density Offsets for 5G Uplink Channels

abstract: This letter proposes a pragmatic link adaptation algorithm considering power density offsets (PDOs) for next-generation uplink wireless channels. The proposed algorithm consists of PDO calculation between a physical uplink shared channel and its associated sounding reference signal, key channel state metric generation, and modulation and coding scheme (MCS) adaptation with respect to the PDO. Scaling is applied to estimated channel matrices based on multiple reference PDO points to generate corresponding reference mutual information (MI) values, followed by interpolation or extrapolation to obtain the adapted MI and ultimately MCS. The proposed algorithm has low complexity in terms of hardware implementation, while yielding satisfactory block error rates and throughput for a wide range of PDOs as shown by simulation results.


title: Building Containerized Environments for Reproducibility and Traceability of Scientific Workflows

abstract: Scientists rely on simulations to study natural phenomena. Trusting the simulation results is vital to develop sciences in any field. One approach to build trust is to ensure the reproducibility and traceability of the simulations through the annotation of executions at the system-level; by the generation of record trails of data moving through the simulation workflow. In this work, we present a system-level solution that leverages the intrinsic characteristics of containers (i.e., portability, isolation, encapsulation, and unique identifiers). Our solution consists of a containerized environment capable to annotate workflows, capture provenance metadata, and build record trails. We assess our environment on four different workflows and measure containerization costs in terms of time and space. Our solution, built with a tolerable time and space overhead, enables transparent and automatic provenance metadata collection and access, an easy-to-read record trail, and tight connections between data and metadata.


title: A Fast and Robust Method for Global Topological Functional Optimization

abstract: Topological statistics, in the form of persistence diagrams, are a class of shape descriptors that capture global structural information in data. The mapping from data structures to persistence diagrams is almost everywhere differentiable, allowing for topological gradients to be backpropagated to ordinary gradients. However, as a method for optimizing a topological functional, this backpropagation method is expensive, unstable, and produces very fragile optima. Our contribution is to introduce a novel backpropagation scheme that is significantly faster, more stable, and produces more robust optima. Moreover, this scheme can also be used to produce a stable visualization of dots in a persistence diagram as a distribution over critical, and near-critical, simplices in the data structure.


title: The Next Big Thing(s) in Unsupervised Machine Learning: Five Lessons from Infant Learning

abstract: After a surge in popularity of supervised Deep Learning, the desire to reduce the dependence on curated, labelled data sets and to leverage the vast quantities of unlabelled data available recently triggered renewed interest in unsupervised learning algorithms. Despite a significantly improved performance due to approaches such as the identification of disentangled latent representations, contrastive learning, and clustering optimisations, the performance of unsupervised machine learning still falls short of its hypothesised potential. Machine learning has previously taken inspiration from neuroscience and cognitive science with great success. However, this has mostly been based on adult learners with access to labels and a vast amount of prior knowledge. In order to push unsupervised machine learning forward, we argue that developmental science of infant cognition might hold the key to unlocking the next generation of unsupervised learning approaches. Conceptually, human infant learning is the closest biological parallel to artificial unsupervised learning, as infants too must learn useful representations from unlabelled data. In contrast to machine learning, these new representations are learned rapidly and from relatively few examples. Moreover, infants learn robust representations that can be used flexibly and efficiently in a number of different tasks and contexts. We identify five crucial factors enabling infants' quality and speed of learning, assess the extent to which these have already been exploited in machine learning, and propose how further adoption of these factors can give rise to previously unseen performance levels in unsupervised learning.


title: Accelerated paths and Unruh effect I: scalars and fermions in Anti De Sitter spacetime

abstract: We have investigated the Unruh effect in Anti de-Sitter (AdS) spacetime by examining the response function of an Unruh-DeWitt particle detector with uniform constant acceleration. An exact expression of the detector response function for the scalar field has been obtained with different levels of non-linearity in even dimensional AdS spacetime. We also showed how the response of the accelerated Unruh detector coupled quadratically to massless Dirac field in $D$ dimensional $(D \geq 2)$ AdS spacetime is proportional to that of a detector linearly coupled to a massless scalar field in $2D$ dimensional AdS spacetime. Here, the fermionic and scalar matter field is coupled minimally and conformally to the background AdS metric, respectively. Finally, we discuss about the extension of the results for more general stationary motion.


title: Consistency of nucleon-transfer sum rules in well-deformed nuclei

abstract: Nucleon-transfer sum rules have been assessed via a consistent reanalysis of cross-section data from neutron-adding ($d$,$p$) and -removing ($d$,$t$) reactions on well-deformed isotopes of Gd, Dy, Er, Yb, and W, with $92\leq N\leq108$, studied at the Niels Bohr Institute in the 1960s and 1970s. These are complemented by new measurements of cross sections using the ($d$,$p$), ($d$,$t$), and ($p$,$d$) reactions on a subset of these nuclei. The sum rules, defined in a Nilsson-model framework, are remarkably consistent. A single overall normalization is used in the analysis, which appears to be sensitive to assumptions about the reaction mechanism, and in the case of sums using the ($d$,$t$) reaction, differs from values determined from reactions on spherical systems.


title: RoadRunner: a fast and flexible exoplanet transit model

abstract: I present RoadRunner, a fast exoplanet transit model that can use any radially symmetric function to model stellar limb darkening while still being faster to evaluate than the analytical transit model for quadratic limb darkening by Mandel & Agol (2002). CPU and GPU implementations of the model are available in the PyTransit transit modelling package, and come with platform-independent parallelisation, supersampling, and support for modelling complex heterogeneous time series. The code is written in numba-accelerated Python (and the GPU model in OpenCL) without C or Fortran dependencies, which allows for the limb darkening model to be given as any Python-callable function. Finally, as an example of the flexibility of the approach, the latest version of PyTransit comes with a numerical limb darkening model that uses LDTk-generated limb darkening profiles directly without approximating them by analytical models.


title: Separating physically distinct mechanisms in complex infrared plasmonic nanostructures via machine learning enhanced electron energy loss spectroscopy

abstract: Low-loss electron energy loss spectroscopy (EELS) has emerged as a technique of choice for exploring the localization of plasmonic phenomena at the nanometer level, necessitating analysis of physical behaviors from 3D spectral data sets. For systems with high localization, linear unmixing methods provide an excellent basis for exploratory analysis, while in more complex systems large numbers of components are needed to accurately capture the true plasmonic response and the physical interpretability of the components becomes uncertain. Here, we explore machine learning based analysis of low-loss EELS data on heterogeneous self-assembled monolayer films of doped-semiconductor nanoparticles, which support infrared resonances. We propose a pathway for supervised analysis of EELS datasets that separate and classify regions of the films with physically distinct spectral responses. The classifications are shown to be robust, to accurately capture the common spatiospectral tropes of the complex nanostructures, and to be transferable between different datasets to allow high-throughput analysis of large areas of the sample. As such, it can be used as a basis for automated experiment workflows based on Bayesian optimization, as demonstrated on the ex situ data. We further demonstrate the use of non-linear autoencoders (AE) combined with clustering in the latent space of the AE yields highly reduced representations of the system response that yield insight into the relevant physics that do not depend on operator input and bias. The combination of these supervised and unsupervised tools provides complementary insight into the nanoscale plasmonic phenomena.


title: Trade-offs in absorption and scattering by nanophotonic structures

abstract: Trade-offs between feasible absorption and scattering cross sections of obstacles confined to an arbitrarily shaped volume are formulated as a multi-objective optimization problem solvable by Lagrangian-dual methods. Solutions to this optimization problem yield a Pareto-optimal set, the shape of which reveals the feasibility of achieving simultaneously extremal absorption and scattering. Two forms of the trade-off problems are considered involving both loss and reactive material parameters. Numerical comparisons between the derived multi-objective bounds and several classes of realized structures are made. Additionally, low-frequency (electrically small, long wavelength) limits are examined for certain special cases.


title: Near-field simulations of pellet ablation for disruptions mitigation in tokamaks

abstract: Detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic and low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons passing through the ablation cloud, expansion of the ablation cloud along magnetic field lines and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich-Kwong corrections to the ideal gas equation of state for cold and dense gases have been used near the pellet surface. The FronTier pellet code is the next generation of the code described in [R. Samulyak, T. Lu, P. Parks, Nuclear Fusion, (47) 2007, 103--118]. It has been validated against the semi-analytic improved Neutral Gas Shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich-Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.


title: Investigation of entanglement measures across the magnetization process of a highly frustrated spin-1/2 Heisenberg octahedral chain as a new paradigm of the localized-magnon approach

abstract: The bipartite entanglement across the magnetization process of a highly frustrated spin-1/2 Heisenberg octahedral chain is examined within the concept of localized magnons, which enables a simple calculation of the concurrence measuring a strength of the pairwise entanglement between nearest-neighbor and next-nearest-neighbor spins from square plaquettes. A full exact diagonalization of the finite-size Heisenberg octahedral chain with up to 4 unit cells (20 spins) evidences an extraordinary high precision of the localized-magnon theory in predicting measures of the bipartite entanglement at sufficiently low temperatures. While the monomer-tetramer phase emergent at low enough magnetic fields exhibits presence (absence) of the bipartite entanglement between the nearest-neighbor (next-nearest-neighbor) spins, the magnon-crystal phase emergent below the saturation field contrarily displays identical bipartite entanglement between the nearest-neighbor and next-nearest-neighbor spins. The presented results verify a new paradigm of the localized-magnon approach concerned with a simple calculation of entanglement measures.


title: A consequence of the relative Bogomolov conjecture

abstract: We propose a formulation of the relative Bogomolov conjecture and show that it gives an affirmative answer to a question of Mazur's concerning the uniformity of the Mordell-Lang conjecture for curves. In particular we show that the relative Bogomolov conjecture implies the uniform Manin-Mumford conjecture for curves. The proof is built up on our previous work "Uniformity in Mordell-Lang for curves".


title: Influence of extended defects on the formation energy, the hyperfine structure, and the zero-field splitting of NV centers in diamond

abstract: We present a density functional theory analysis of nitrogen-vacancy (NV) centers in diamond which are located in the vicinity of extended defects, namely intrinsic stacking faults (ISF), extrinsic stacking faults (ESF), and coherent twin boundaries (CTB) on \{111\} planes in diamond crystals. Several sites for NV centers close to the extended defects are energetically preferred with respect to the bulk crystal. This indicates that NV centers may be enriched at extended defects. We report the hyperfine structure (HFS) and zero-field splitting (ZFS) parameters of the NV centers at the extended defects which typically deviate by about 10\% but in some cases up to 90\% from their bulk values. Furthermore, we find that the influence of the extended defects on the NV centers is of short range: NV centers that are about three double layers (corresponding to $\sim6$ $\mathring{A}$ away from defect planes already show bulk-like behavior.


title: Reconstructing Actions To Explain Deep Reinforcement Learning

abstract: Feature attribution has been a foundational building block for explaining the input feature importance in supervised learning with Deep Neural Network (DNNs), but face new challenges when applied to deep Reinforcement Learning (RL).We propose a new approach to explaining deep RL actions by defining a class of \emph{action reconstruction} functions that mimic the behavior of a network in deep RL. This approach allows us to answer more complex explainability questions than direct application of DNN attribution methods, which we adapt to \emph{behavior-level attributions} in building our action reconstructions. It also allows us to define \emph{agreement}, a metric for quantitatively evaluating the explainability of our methods. Our experiments on a variety of Atari games suggest that perturbation-based attribution methods are significantly more suitable in reconstructing actions to explain the deep RL agent than alternative attribution methods, and show greater \emph{agreement} than existing explainability work utilizing attention. We further show that action reconstruction allows us to demonstrate how a deep agent learns to play Pac-Man game.


title: Boundary effects on symmetry resolved entanglement

abstract: We study the symmetry resolved entanglement entropies in one-dimensional systems with boundaries. We provide some general results for conformal invariant theories and then move to a semi-infinite chain of free fermions. We consider both an interval starting from the boundary and away from it. We derive exact formulas for the charged and symmetry resolved entropies based on theorems and conjectures about the spectra of Toeplitz+Hankel matrices. En route to characterise the interval away from the boundary, we prove a general relation between the eigenvalues of Toeplitz+Hankel matrices and block Toeplitz ones. An important aspect is that the saddle-point approximation from charged to symmetry resolved entropies introduces algebraic corrections to the scaling that are much more severe than in systems without boundaries.


title: Testing eigenstate decoherence hypothesis in a model of collisional decoherence

abstract: The eigenstate decoherence hypothesis (EDH) asserts that each individual eigenstate of a large closed system is locally classical-like. We test this hypothesis for a heavy particle interacting with a gas of light particles. This system is paradigmatic for studies of the quantum-to-classical transition: The reduced state of the heavy particle is widely believed to rapidly loose any nonclassical features due to the interaction with the gas. Yet, we find numerical evidence that the EDH is violated: certain eigenstates of this model are manifestly non-classical. Only the weak version of EDH referring to the majority (instead of the totality) of eigenstates holds.


title: Automatic deep learning for trend prediction in time series data

abstract: Recently, Deep Neural Network (DNN) algorithms have been explored for predicting trends in time series data. In many real world applications, time series data are captured from dynamic systems. DNN models must provide stable performance when they are updated and retrained as new observations becomes available. In this work we explore the use of automatic machine learning techniques to automate the algorithm selection and hyperparameter optimisation process for trend prediction. We demonstrate how a recent AutoML tool, specifically the HpBandSter framework, can be effectively used to automate DNN model development. Our AutoML experiments found optimal configurations that produced models that compared well against the average performance and stability levels of configurations found during the manual experiments across four data sets.


title: Smartphone Camera De-identification while Preserving Biometric Utility

abstract: The principle of Photo Response Non Uniformity (PRNU) is often exploited to deduce the identity of the smartphone device whose camera or sensor was used to acquire a certain image. In this work, we design an algorithm that perturbs a face image acquired using a smartphone camera such that (a) sensor-specific details pertaining to the smartphone camera are suppressed (sensor anonymization); (b) the sensor pattern of a different device is incorporated (sensor spoofing); and (c) biometric matching using the perturbed image is not affected (biometric utility). We employ a simple approach utilizing Discrete Cosine Transform to achieve the aforementioned objectives. Experiments conducted on the MICHE-I and OULU-NPU datasets, which contain periocular and facial data acquired using 12 smartphone cameras, demonstrate the efficacy of the proposed de-identification algorithm on three different PRNU-based sensor identification schemes. This work has application in sensor forensics and personal privacy.


title: Intelligent Reflecting Surface Aided Pilot Contamination Attack and Its Countermeasure

abstract: Pilot contamination attack (PCA) in a time division duplex wireless communication system is considered, where an eavesdropper (Eve) attacks the reverse pilot transmission phase in order to wiretap the data transmitted from a transmitter, Alice, to a receiver, Bob. We propose a new PCA scheme for Eve, wherein Eve does not emit any signal by itself but uses an intelligent reflecting surface (IRS) to reflect the pilot sent by Bob to Alice. The proposed new PCA scheme, referred to as IRS-PCA, increases the signal leakage from Alice to the IRS during the data transmission phase, which is then reflected by the IRS to Eve in order to improve the wiretapping capability of Eve. The proposed IRS-PCA scheme disables many existing countermeasures on PCA due to the fact that with IRS-PCA, Eve no longer needs to know the pilot sequence of Bob, and therefore, poses severe threat to the security of the legitimate wireless communication system. In view of this, the problems of 1) IRS-PCA detection and 2) secure transmission under IRSPCA are considered in this paper. For IRS-PCA detection, a generalized cumulative sum (GCUSUM) detection procedure is proposed based on the framework of quickest detection, aiming at detecting the occurrence of IRS-PCA as soon as possible once it occurs. For secure transmission under IRS-PCA, a cooperative channel estimation scheme is proposed to estimate the channel of the IRS, based on which zero-forcing beamforming is designed to reduce signal leakage.


title: Practical Quantum Computing: The value of local computation

abstract: As we enter the era of useful quantum computers we need to better understand the limitations of classical support hardware, and develop mitigation techniques to ensure effective qubit utilisation. In this paper we discuss three key bottlenecks in near-term quantum computers: bandwidth restrictions arising from data transfer between central processing units (CPUs) and quantum processing units (QPUs), latency delays in the hardware for round-trip communication, and timing restrictions driven by high error rates. In each case we consider a near-term quantum algorithm to highlight the bottleneck: randomised benchmarking to showcase bandwidth limitations, adaptive noisy, intermediate scale quantum (NISQ)-era algorithms for the latency bottleneck and quantum error correction techniques to highlight the restrictions imposed by qubit error rates. In all three cases we discuss how these bottlenecks arise in the current paradigm of executing all the classical computation on the CPU, and how these can be mitigated by providing access to local classical computational resources in the QPU.


title: Active particles in linear viscoelastic fluids and the scallop theorem

abstract: We derive a general formula for the inertialess dynamics of active particles in linear viscoelastic fluids by means of a modified reciprocal theorem. We then demonstrate that force-free active particles in Maxwell-like linear viscoelastic fluids with no retardation have exactly the same dynamics as in Newtonian fluids. In contrast, active particles in Jeffreys-like fluids with retardation can display markedly different dynamics, including net motion under reciprocal actuation thereby breaking the scallop theorem.


title: A computational framework for evaluating the role of mobility on the propagation of epidemics on point processes

abstract: This paper is focused on SIS (Susceptible-Infected-Susceptible) epidemic dynamics (also known as the contact process) on populations modelled by homogeneous Poisson point processes of the Euclidean plane, where the infection rate of a susceptible individual is proportional to the number of infected individuals in a disc around it. The main focus of the paper is a model where points are also subject to some random motion. Conservation equations for moment measures are leveraged to analyze the stationary regime of the point processes of infected and susceptible individuals. A heuristic factorization of the third moment measure is then proposed to obtain simple polynomial equations allowing one to derive closed form approximations for the fraction of infected individuals in the steady state. These polynomial equations also lead to a phase diagram which tentatively delineates the regions of the space of parameters (population density, infection radius, infection and recovery rate, and motion rate) where the epidemic survives and those where there is extinction. A key take-away from this phase diagram is that the extinction of the epidemic is not always aided by a decrease in the motion rate. These results are substantiated by simulations on large two dimensional tori. These simulations show that the polynomial equations accurately predict the fraction of infected individuals when the epidemic survives. The simulations also show that the proposed phase diagram accurately predicts the parameter regions where the mean survival time of the epidemic increases (resp. decreases) with motion rate.


title: Particle acceleration by relativistic magnetic reconnection driven by kink instability turbulence in Poynting flux dominated jets

abstract: Particle acceleration in magnetized relativistic jets still puzzles theorists, specially when one tries to explain the highly variable emission observed in blazar jets or gamma-ray bursts putting severe constraints on current models. In this work we investigate the acceleration of particles injected in a three-dimensional relativistic magnetohydrodynamical jet subject to current driven kink instability (CDKI), which drives turbulence and fast magnetic reconnection. Test protons injected in the nearly stationary snapshots of the jet, experience an exponential acceleration up to a maximum energy. For a background magnetic field of $B \sim 0.1$ G, this saturation energy is $\sim 10^{16}$ eV, while for $B \sim 10$ G it is $\sim 10^{18}$ eV. The simulations also reveal a clear association of the accelerated particles with the regions of fast reconnection. In the early stages of the development of the non-linear growth of CDKI in the jet, when there are still no sites of fast reconnection, injected particles are also efficiently accelerated, but by magnetic curvature drift in the wiggling jet spine. However, they have to be injected with an initial energy much larger than that required for particles to accelerate in reconnection sites. Finally, we have also obtained from the simulations an acceleration time due to reconnection with a weak dependence on the particles energy $E$, $t_A \propto E^{0.1}$. The energy spectrum of the accelerated particles develops a high energy tail with a power law index $p \sim$ -1.2 in the beginning of the acceleration, in agreement with earlier works. Our results provide an appropriate multi-dimensional framework for exploring this process in real systems and explain their complex emission patterns, specially in the very high energy bands and the associated neutrino emission recently detected in some blazars.


title: Extracting chromatic properties of electron beams from spectral analysis of turn-by-turn beam position data

abstract: A method to estimate linear chromaticity, RMS energy spread, and chromatic beta-beating, directly from turn-by-turn beam position data in a circular electron accelerator, is presented. This technique is based on frequency analysis of a transversely excited beam, in the presence of finite chromaticity. Due to the turn-by-turn chromatic modulation of the beam's envelope, betatron sidebands appear around the main frequency of the Fourier spectra. By determining the amplitude of both sidebands, chromatic properties of the beam can be estimated. In this paper, analytical derivations justifying the proposed method are given, along with results from tracking simulations. To this end, results from practical applications of this technique at the KARA electron ring are demonstrated.


title: Slow scrambling in extremal BTZ and microstate geometries

abstract: Out-of-time-order correlators (OTOCs) that capture maximally chaotic properties of a black hole are determined by scattering processes near the horizon. This prompts the question to what extent OTOCs display chaotic behaviour in horizonless microstate geometries. This question is complicated by the fact that Lyapunov growth of OTOCs requires nonzero temperature, whereas constructions of microstate geometries have been mostly restricted to extremal black holes. In this paper, we compute OTOCs for a class of extremal black holes, namely maximally rotating BTZ black holes, and show that on average they display "slow scrambling", characterized by cubic (rather than exponential) growth. Superposed on this average power-law growth is a sawtooth pattern, whose steep parts correspond to brief periods of Lyapunov growth associated to the nonzero temperature of the right-moving degrees of freedom in a dual conformal field theory. Next we study the extent to which these OTOCs are modified in certain "superstrata", horizonless microstate geometries corresponding to these black holes. Rather than an infinite throat ending on a horizon, these geometries have a very deep but finite throat ending in a cap. We find that the superstrata display the same slow scrambling as maximally rotating BTZ black holes, except that for large enough time intervals the growth of the OTOC is cut off by effects related to the cap region, some of which we evaluate explicitly.


title: Properties of the density of shear transformations in driven amorphous solids

abstract: The strain load $\Delta\gamma$ that triggers consecutive avalanches is a key observable in the slow deformation of amorphous solids. Its temporally averaged value $\langle \Delta\gamma \rangle$ displays a non-trivial system-size dependence that constitutes one of the distinguishing features of the yielding transition. Details of this dependence are not yet fully understood. We address this problem by means of theoretical analysis and simulations of elastoplastic models for amorphous solids. An accurate determination of the size dependence of $\langle \Delta\gamma \rangle$ leads to a precise evaluation of the steady-state distribution of local distances to instability $x$. We find that the usually assumed form $P(x)\sim x^\theta$ (with $\theta$ being the so-called pseudo-gap exponent) is not accurate at low $x$ and that in general $P(x)$ tends to a system-size-dependent \textit{finite} limit as $x\to 0$. We work out the consequences of this finite-size dependence standing on exact results for random-walks and disclosing an alternative interpretation of the mechanical noise felt by a reference site. We test our predictions in two- and three-dimensional elastoplastic models, showing the crucial influence of the saturation of $P(x)$ at small $x$ on the size dependence of $\langle \Delta\gamma \rangle$ and related scalings.


title: Skein lasagna modules for 2-handlebodies

abstract: Morrison, Walker, and Wedrich used the blob complex to construct a generalization of Khovanov-Rozansky homology to links in the boundary of a 4-manifold. The degree zero part of their theory, called the skein lasagna module, admits an elementary definition in terms of certain diagrams in the 4-manifold. We give a description of the skein lasagna module for 4-manifolds without 1- and 3-handles, and present some explicit calculations for disk bundles over $S^2$.


title: Riesz bases of port-Hamiltonian systems

abstract: The location of the spectrum and the Riesz basis property of well-posed homogeneous infinite-dimensional linear port-Hamiltonian systems on a 1D spatial domain are studied. It is shown that the Riesz basis property is equivalent to the fact that system operator generates a strongly continuous group. Moreover, in this situation the spectrum consists of eigenvalues only, located in a strip parallel to the imaginary axis and they can decomposed into finitely many sets having each a uniform gap.


title: A categorical sl_2 action on some moduli spaces of sheaves

abstract: We study certain sequences of moduli spaces of sheaves on K3 surfaces, building on work of Markman, Yoshioka, and Nakajima. We show that these sequences can be given the structure of a geometric categorical sl_2 action in the sense of Cautis, Kamnitzer, and Licata. As a corollary, we get an equivalence between derived categories of some moduli spaces that are birational via stratified Mukai flops.


title: Quantum scarring in a spin-boson system: fundamental families of periodic orbits

abstract: As the name indicates, a periodic orbit is a solution for a dynamical system that repeats itself in time. In the regular regime, periodic orbits are stable, while in the chaotic regime, they become unstable. The presence of unstable periodic orbits is directly associated with the phenomenon of quantum scarring, which restricts the degree of delocalization of the eigenstates and leads to revivals in the dynamics. Here, we study the Dicke model in the superradiant phase and identify two sets of fundamental periodic orbits. This experimentally realizable atom-photon model is regular at low energies and chaotic at high energies. We study the effects of the periodic orbits in the structure of the eigenstates in both regular and chaotic regimes and obtain their quantized energies. We also introduce a measure to quantify how much scarred an eigenstate gets by each family of periodic orbits and compare the dynamics of initial coherent states close and away from those orbits.


title: Emergence of multistability and strongly asymmetric collective modes in two quorum sensing coupled identical ring oscillators

abstract: The simplest ring oscillator is made from three strongly nonlinear elements repressing each other unidirectionally resulting in the emergence of a limit cycle. A popular implementation of this scheme uses repressive genes in bacteria creating the synthetic genetic oscillator known as the Repressilator. Here, we consider the main collective modes produced when two identical Repressilators are mean-field coupled via the quorum sensing (QS) mechanism which is realized via production of diffusive signal molecules. Using the rate of the repressor's production and the value of coupling strength as the bifurcation parameters, we performed analysis of dynamical regimes starting from the two Andronov-Hopf bifurcations of unstable homogeneous steady state, which generate in-phase and anti-phase limit cycles. Pitchfork bifurcation of the unstable in-phase cycle leads to creation of inhomogeneous limit cycles with very different amplitudes in contrast to well-known asymmetrical limit cycles arising from oscillation death. Neimark-Sacker bifurcation of the anti-phase cycle determines the border of an island in two-parameter space containing almost all the interesting regimes including the set of resonant limit cycles, the area with stable inhomogenous cycle, and very large areas with chaotic regimes resulting from torus destruction, period doubling of resonant cycles and inhomogenous cycles. We discuss the structure of chaos skeleton to show the role of inhomogeneous cycles in its formation. Many regions of multistability and transitions between regimes are presented. These results provide new insights into the coupling-dependent mechanisms of multistability and collective regime symmetry breaking in populations of identical multidimensional oscillators.


title: Deep Learning & Software Engineering: State of Research and Future Directions

abstract: Given the current transformative potential of research that sits at the intersection of Deep Learning (DL) and Software Engineering (SE), an NSF-sponsored community workshop was conducted in co-location with the 34th IEEE/ACM International Conference on Automated Software Engineering (ASE'19) in San Diego, California. The goal of this workshop was to outline high priority areas for cross-cutting research. While a multitude of exciting directions for future work were identified, this report provides a general summary of the research areas representing the areas of highest priority which were discussed at the workshop. The intent of this report is to serve as a potential roadmap to guide future work that sits at the intersection of SE & DL.


title: The equivariant cohomology for semidirect product actions

abstract: The rational Borel equivariant cohomology for actions of a compact connected Lie group is determined by restriction of the action to a maximal torus. We show that a similar reduction holds for any compact Lie group $G$ when there is a closed subgroup $K$ such that the cohomology of the classifying space $BK$ is free over the cohomology of $BG$ for field coefficients. We study the particular case when $G$ is a semi-direct product and $K$ is its maximal elementary abelian 2-subgroup for cohomology with coefficients in a field of characteristic two. This provides a different approach to investigate the syzygy order of the equivariant cohomology of a space with a torus action and a compatible involution, and we relate this description with results for 2-torus actions.


title: Realizations of non-commutative rational functions around a matrix centre, II: The lost-abbey conditions

abstract: In a previous paper the authors generalized classical results of minimal realizations of non-commutative (nc) rational functions, using nc Fornasini--Marchesini realizations which are centred at an arbitrary matrix point. In particular, it was proved that the domain of regularity of a nc rational function is contained in the invertibility set of a corresponding pencil of any minimal realization of the function. In this paper we prove an equality between the domain of a nc rational function and the domain of any of its minimal realizations. As for evaluations over stably finite algebras, we show that the domain of the realization w.r.t any such algebra coincides with the so called matrix domain of the function w.r.t the algebra. As a corollary we show that the domain of regularity and the stable extended domain coincide. In contrary to both the classical case and the scalar case -- where every matrix coefficients which satisfy the controllability and observability conditions can appear in a minimal realization of a nc rational function -- the matrix coefficients in our case have to satisfy certain equations, called linearized lost-abbey conditions, which are related to Taylor--Taylor expansions in nc function theory.


title: SREC: Proactive Self-Remedy of Energy-Constrained UAV-Based Networks via Deep Reinforcement Learning

abstract: Energy-aware control for multiple unmanned aerial vehicles (UAVs) is one of the major research interests in UAV based networking. Yet few existing works have focused on how the network should react around the timing when the UAV lineup is changed. In this work, we study proactive self-remedy of energy-constrained UAV networks when one or more UAVs are short of energy and about to quit for charging. We target at an energy-aware optimal UAV control policy which proactively relocates the UAVs when any UAV is about to quit the network, rather than passively dispatches the remaining UAVs after the quit. Specifically, a deep reinforcement learning (DRL)-based self remedy approach, named SREC-DRL, is proposed to maximize the accumulated user satisfaction scores for a certain period within which at least one UAV will quit the network. To handle the continuous state and action space in the problem, the state-of-the-art algorithm of the actor-critic DRL, i.e., deep deterministic policy gradient (DDPG), is applied with better convergence stability. Numerical results demonstrate that compared with the passive reaction method, the proposed SREC-DRL approach shows a $12.12\%$ gain in accumulative user satisfaction score during the remedy period.


title: Complex solid solution electrocatalyst discovery by prediction and high-throughput experimentation

abstract: Efficient discovery of electrocatalysts for electrochemical energy conversion reactions is of utmost importance to combat climate change. With the example of the oxygen reduction reaction we show that by utilising a data-driven discovery cycle, the multidimensionality challenge offered by compositionally complex solid solution (high entropy alloy) electrocatalysts can be mastered. Iteratively refined computational models predict activity trends for quinary target compositions, around which continuous composition spread thin-film libraries are synthesized. High-throughput characterisation datasets are then input for refinement of the model. The refined model correctly predicts activity maxima of the exemplary model system Ag-Ir-Pd-Pt-Ru for the oxygen reduction reaction. The method can identify optimal complex solid solutions for electrochemical reactions in an unprecedented manner.


title: The quotient criterion for syzygies in equivariant cohomology for elementary abelian $2$-group actions

abstract: Let $G$ be a elementary abelian $2$-group and $X$ be a manifold with a locally standard action of $G$. We provide a criterion to determine the syzygy order of the $G$-equivariant cohomology of $X$ with coefficients over a field of characteristic two using a complex associated to the cohomology of the face filtration of the manifold with corners $X/G$. This result is the real version of the quotient criterion for locally standard torus actions developed by M. Franz in arXiv:1205.4462.


title: Positive Solutions For a Schr\"odinger-Bopp-Podolsky system in $\mathbb R^{3}$

abstract: We consider the following Schr\"odinger-Bopp-Podolsky system in $\mathbb R^{3}$ $$\left\{ \begin{array}{c} -\varepsilon^{2} \Delta u + V(x)u + \phi u = f(u)\\ -\varepsilon^{2} \Delta \phi + \varepsilon^{4} \Delta^{2}\phi = 4\pi\varepsilon u^{2}\\ \end{array} \right.$$ where $\varepsilon > 0$ with $ V:\mathbb{R}^{3} \rightarrow \mathbb{R}, f:\mathbb{R} \rightarrow \mathbb{R}$ satisfy suitable assumptions. By using variational methods, we prove that the number of positive solutions is estimated below by the Ljusternick-Schnirelmann category of $M$, the set of minima of the potential $V$.


title: Radio number of Hamming graphs of diameter 3

abstract: For $G$ a simple, connected graph, a vertex labeling $f:V(G)\rightarrow \mathbb{Z}_+$ is called a $\textit{radio labeling of}$ $G$ if it satisfies $|f(u)-f(v)|\geq \operatorname{diam}(G) + 1 - d(u,v)$ for all distinct vertices $u,v\in V(G)$. The $\textit{radio number}$ of $G$ is the minimal span over all radio labelings of $G$. If a bijective radio labeling onto $\{1,2,...,|V(G)|\}$ exists, $G$ is called a $\textit{radio graceful graph}$. We determine the radio number of all diameter $3$ Hamming graphs and show that an infinite subset of them is radio graceful.


title: Robust Asymptotic Growth in Stochastic Portfolio Theory under Long-Only Constraints

abstract: We consider the problem of maximizing the asymptotic growth rate of an investor under drift uncertainty in the setting of stochastic portfolio theory (SPT). As in the work of Kardaras and Robertson we take as inputs (i) a Markovian volatility matrix $c(x)$ and (ii) an invariant density $p(x)$ for the market weights, but we additionally impose long-only constraints on the investor. Our principal contribution is proving a uniqueness and existence result for the class of concave functionally generated portfolios and developing a finite dimensional approximation, which can be used to numerically find the optimum. In addition to the general results outlined above, we propose the use of a broad class of models for the volatility matrix $c(x)$, which can be calibrated to data and, under which, we obtain explicit formulas of the optimal unconstrained portfolio for any invariant density.


title: Assessment of Efficiency, Impact Factor, Impact of Probe Mass, Probe Life Expectancy, and Reliability of Mars Missions

abstract: Mars is the next frontier after Moon for space explorers to demonstrate the extent of human expedition and technology beyond low-earth orbit. Government space agencies as well as private space sectors are extensively endeavouring for a better space enterprise. Focusing on the inspiration to reach Mars by robotic satellite, we have interpreted some of the significant mission parameters like proportionality of mission attempts, efficiency and reliability of Mars probes, Impact and Impact Factor of mass on mission duration, Time lag between consecutive mission attempts, interpretation of probe life and transitional region employing various mathematical analysis. And we have discussed the importance of these parameters for a prospective mission accomplishment. Our novelty in this paper is we have found a deep relation describing that the probe mass adversely affects the mission duration. Applying this relation, we also interpreted the duration of probe life expectancy for upcoming missions.


title: Numerical Testing of a New Positivity-Preserving Interpolation Algorithm

abstract: An important component of a number of computational modeling algorithms is an interpolation method that preserves the positivity of the function being interpolated. This report describes the numerical testing of a new positivity-preserving algorithm that is designed to be used when interpolating from a solution defined on one grid to different spatial grid. The motivating application for this work was a numerical weather prediction (NWP) code that uses a spectral element mesh discretization for its dynamics core and a cartesian tensor product mesh for the evaluation of its physics routines. This coupling of spectral element mesh, which uses nonuniformly spaced quadrature/collocation points, and uniformly-spaced cartesian mesh combined with the desire to maintain positivity when moving between these meshes necessitates our work. This new approach is evaluated against several typical algorithms in use on a range of test problems in one or more space dimensions. The results obtained show that the new method is competitive in terms of observed accuracy while at the same time preserving the underlying positivity of the functions being interpolated.


title: A stabilizer free weak Galerkin finite element method on polytopal mesh: Part III

abstract: A weak Galerkin (WG) finite element method without stabilizers was introduced in [J. Comput. Appl. Math., 371 (2020). arXiv:1906.06634] on polytopal mesh. Then it was improved in [arXiv:2008.13631] with order one superconvergence. The goal of this paper is to develop a new stabilizer free WG method on polytopal mesh. This method has convergence rates two orders higher than the optimal convergence rates for the corresponding WG solution in both an energy norm and the $L^2$ norm. The numerical examples are tested for low and high order elements in two and three dimensional spaces.


title: Direct Imaging Discovery of a Young Brown Dwarf Companion to an A2V Star

abstract: We present the discovery and spectroscopy of HIP 75056Ab, a companion directly imaged at a very small separation of 0.125 arcsec to an A2V star in the Scorpius-Centaurus OB2 association. Our observations utilized VLT/SPHERE between 2015$-$2019, enabling low-resolution spectroscopy (0.95$-$1.65 $\mu m$), dual-band imaging (2.1$-$2.25 $\mu m$), and relative astrometry over a four-year baseline. HIP 75056Ab is consistent with spectral types in the range of M6$-$L2 and $T_{\rm eff}\sim$ 2000$-$2600 K. A comparison of the companion's brightness to evolutionary tracks suggests a mass of $\sim$20$-$30 M$_{Jup}$. The astrometric measurements are consistent with an orbital semi-major axis of $\sim$15$-$45 au and an inclination close to face-on (i$\lesssim$35$^o$). In this range of mass and orbital separation, HIP 75056Ab is likely at the low-mass end of the distribution of companions formed via disk instability, although a formation of the companion via core accretion cannot be excluded. The orbital constraints are consistent with the modest eccentricity values predicted by disk instability, a scenario that can be confirmed by further astrometric monitoring. HIP 75056Ab may be utilized as a low-mass atmospheric comparison to older, higher-mass brown dwarfs, and also to young giant planets. Finally, the detection of HIP 75056Ab at 0.125 arcsec represents a milestone in detecting low-mass companions at separations corresponding to the habitable zones of nearby Sun-like stars.


title: Resonant scattering of electromagnetic waves by small metal particles

abstract: The review is devoted to a discussion of new (and often unexpected) aspects of the old problem of elastic light scattering by small metal particles, whose size is comparable to or smaller than the thickness of the skin layer. The main focus is put on elucidating the physical grounds for these new aspects. It is shown that, in many practically important cases, the scattering of light by such particles, despite their smallness, may have almost nothing in common with the Rayleigh one. The so-called, anomalous scattering and absorption, as well as Fano resonances, including unconventional (associated with the excitation of longitudinal electromagnetic oscillations) and directional Fano resonances, observed only in a small solid angle, are discussed in detail. The review contains a Mathematical Supplement, which includes a summary of the main results of the Mie theory and a discussion of some general properties of the scattering coefficients. In addition to purely academic interest, the phenomena considered in this review can find wide applications in biology, medicine, pharmacology, genetic engineering, imaging of ultra-small objects, ultra-high-resolution spectroscopy, information transmission, recording, and processing, and many other applications and technologies. The reported study was funded by RFBR, project number 19-11-00001 and the project of the Russian Science Foundation No. 19-72-30012, within the framework of which all the original calculations given in this publication were performed.

categories: physics.optics physics.app-ph physics.med-ph physics.plasm-ph


title: Objective, Probabilistic, and Generalized Noise Level Dependent Classifications of sets of more or less 2D Periodic Images into Plane Symmetry Groups

abstract: Crystallographic symmetry classifications from real-world images with periodicities in two dimensions (2D) are of interest to crystallographers and practitioners of computer vision studies alike. Currently, these classifications are typically made by both communities in a subjective manner that relies on arbitrary thresholds for judgments, and are reported under the pretense of being definitive, which is impossible. Moreover, the computer vision community tends to use direct space methods to make such classifications instead of more powerful and computationally efficient Fourier space methods. This is because the proper functioning of those methods requires more periodic repeats of a unit cell motif than are commonly present in images analyzed by the computer vision community. We demonstrate a novel approach to plane symmetry group classifications that is enabled by Kenichi Kanatani's Geometric Akaike Information Criterion and associated Geometric Akaike weights. Our approach leverages the advantages of working in Fourier space, is well suited for handling the hierarchic nature of crystallographic symmetries, and yields probabilistic results that are generalized noise level dependent. The latter feature means crystallographic symmetry classifications can be updated when less noisy image data and more accurate processing algorithms become available. We demonstrate the ability of our approach to objectively estimate the plane symmetry and pseudosymmetries of sets of synthetic 2D-periodic images with varying amounts of red-green-blue and spread noise. Additionally, we suggest a simple solution to the problem of too few periodic repeats in an input image for practical application of Fourier space methods. In doing so, we effectively solve the decades-old and heretofore intractable problem from computer vision of symmetry detection and classification from images in the presence of noise.


title: Partial Actions of a Hopf algebra on its base field and the corresponding partial smash product algebra

abstract: We introduce the concept of a $\lambda$-Hopf algebra as a Hopf algebra obtained as the partial smash product algebra of a Hopf algebra and its base field, and show that every Hopf algebra is a $\lambda$-Hopf algebra. Moreover, a method to compute partial actions of a given Hopf algebra on its base field is developed and, as an application, we exhibit all partial actions of such type for some families of Hopf algebras.


title: Variational Disentanglement for Rare Event Modeling

abstract: Combining the increasing availability and abundance of healthcare data and the current advances in machine learning methods have created renewed opportunities to improve clinical decision support systems. However, in healthcare risk prediction applications, the proportion of cases with the condition (label) of interest is often very low relative to the available sample size. Though very prevalent in healthcare, such imbalanced classification settings are also common and challenging in many other scenarios. So motivated, we propose a variational disentanglement approach to semi-parametrically learn from rare events in heavily imbalanced classification problems. Specifically, we leverage the imposed extreme-distribution behavior on a latent space to extract information from low-prevalence events, and develop a robust prediction arm that joins the merits of the generalized additive model and isotonic neural nets. Results on synthetic studies and diverse real-world datasets, including mortality prediction on a COVID-19 cohort, demonstrate that the proposed approach outperforms existing alternatives.


title: The Poincare lemma for codifferential, anticoexact forms, and applications to physics

abstract: The linear homotopy theory for codifferential operator on Riemannian manifolds is developed in analogy to a similar idea for exterior derivative. The main object is the cohomotopy operator, which singles out a module of anticoexact forms from the module of differential forms defined on a star-shaped open subset of a manifold. It is shown that there is a direct sum decomposition of a differential form into coexact and anticoexat parts. This decomposition gives a new way of solving exterior differential systems. The method is applied to equations of fundamental physics, including vacuum Dirac-K\"{a}hler equation, coupled Maxwell-Kalb-Ramond system of equations occurring in a bosonic string theory and its reduction to the Dirac equation.


title: Estimation of Health and Demographic Indicators with Incomplete Geographic Information

abstract: In low and middle income countries, household surveys are a valuable source of information for a range of health and demographic indicators. Increasingly, subnational estimates are required for targeting interventions and evaluating progress towards targets. In the majority of cases, stratified cluster sampling is used, with clusters corresponding to enumeration areas. The reported geographical information varies. A common procedure, to preserve confidentiality, is to give a jittered location with the true centroid of the cluster is displaced under a known algorithm. An alternative situation, which was used for older surveys in particular, is to report the geographical region within the cluster lies. In this paper, we describe a spatial hierarchical model in which we account for inaccuracies in the cluster locations. The computational algorithm we develop is fast and avoids the heavy computation of a pure MCMC approach. We illustrate by simulation the benefits of the model, over naive alternatives.


title: Charge-regulation effects in nanoparticle self-assembly

abstract: Nanoparticles in solution acquire charge through dissociation or association of surface groups. Thus, a proper description of their electrostatic interactions requires the use of charge-regulating boundary conditions rather than the commonly employed constant-charge approximation. We implement a hybrid Monte Carlo/Molecular Dynamics scheme that dynamically adjusts the charges of individual surface groups of objects while evolving their trajectories. Charge-regulation effects are shown to qualitatively change self-assembled structures due to global charge redistribution, stabilizing asymmetric constructs. We delineate under which conditions the conventional constant-charge approximation may be employed and clarify the interplay between charge regulation and dielectric polarization.


title: Asymptotic Performance Prediction for ADMM-Based Compressed Sensing

abstract: In this paper, we propose a method to predict the asymptotic performance of the alternating direction method of multipliers (ADMM) for compressed sensing, where we reconstruct an unknown structured signal from its underdetermined linear measurements. The derivation of the proposed method is based on the recently developed convex Gaussian min-max theorem (CGMT), which can be applied to various convex optimization problems to obtain its asymptotic error performance. Our main idea is to analyze the convex subproblem in the update of ADMM iteratively and characterize the asymptotic distribution of the tentative estimate obtained at each iteration. However, since the original CGMT cannot be used directly for the analysis of the iterative updates, we intuitively assume an extended version of CGMT in the derivation of the proposed method. Under the assumption, the result shows that the update equations in ADMM can be decoupled into a scalar-valued stochastic process in the asymptotic regime with the large system limit. From the asymptotic result, we can predict the evolution of the error (e.g., mean-square-error (MSE) and symbol error rate (SER)) in ADMM for large-scale compressed sensing problems. Simulation results show that the empirical performance of ADMM and its prediction are close to each other in sparse vector reconstruction and binary vector reconstruction.


title: Ankara \"{U}niversitesi Kreiken Rasathanesinde Bulunan Radyo Teleskobun Y\"{o}nlendirme Mekanizmas{\i}n{\i}n \.{I}ncelenip Yenilenmesi, Anten Benzetimleri ve Uyumlama Tasar{\i}m{\i}

abstract: To make the historical radio telescope of A.\"{U}. Kreiken Observatory that includes full-wavelength VHF and UHF dipole arrays operational, the antenna rotator mechanism is renovated and the antenna matching structure is designed in a simulation environment. The matching structure is designed by using short-circuited parallel transmission line stubs. The use of high-gain full-wavelength elements, as opposed to the conventional half-wavelength elements, increases the importance of the structure and the matching method. -- A.\"{U}. Kreiken Rasathanesinde bulunan, tam dalgaboyunda dipol elemanlara sahip VHF ve UHF anten dizilerinden olu\c{s}an tarihi radyo teleskobun faal hale getirilmesi i\c{c}in anten y\"{o}nlendirme mekanizmasi incelenip tamir edilmi\c{s}, teleskop anteninin uyumlama tasarimi benzetim ortaminda ger\c{c}ekle\c{s}tirilmi\c{s}tir. Kisa devre ile sonlandirilmi\c{s} paralel iletim hatti saplamalari kullanarak uyumlama tasarlanmi\c{s}tir. Dipol anten elemanlarinin boyunun yaygin tercih edilen yarim dalgaboyu yerine y\"{u}ksek kazan\c{c}li tam dalgaboyu olmasi, antenin ve uyumlama y\"{o}nteminin \"{o}nemini artirmaktadir.


title: Practical Dynamic SC-Flip Polar Decoders: Algorithm and Implementation

abstract: SC-Flip (SCF) is a low-complexity polar code decoding algorithm with improved performance, and is an alternative to high-complexity (CRC)-aided SC-List (CA-SCL) decoding. However, the performance improvement of SCF is limited since it can correct up to only one channel error ($\omega=1$). Dynamic SCF (DSCF) algorithm tackles this problem by tackling multiple errors ($\omega \geq 1$), but it requires logarithmic and exponential computations, which make it infeasible for practical applications. In this work, we propose simplifications and approximations to make DSCF practically feasible. First, we reduce the transcendental computations of DSCF decoding to a constant approximation. Then, we show how to incorporate special node decoding techniques into DSCF algorithm, creating the Fast-DSCF decoding. Next, we reduce the search span within the special nodes to further reduce the computational complexity. Following, we describe a hardware architecture for the Fast-DSCF decoder, in which we introduce additional simplifications such as metric normalization and sorter length reduction. All the simplifications and approximations are shown to have minimal impact on the error-correction performance, and the reported Fast-DSCF decoder is the only SCF-based architecture that can correct multiple errors. The Fast-DSCF decoders synthesized using TSMC $65$nm CMOS technology can achieve a $1.25$, $1.06$ and $0.93$ Gbps throughput for $\omega \in \{1,2,3\}$, respectively. Compared to the state-of-the-art fast CA-SCL decoders with equivalent FER performance, the proposed decoders are up to $5.8\times$ more area-efficient. Finally, observations at energy dissipation indicate that the Fast-DSCF is more energy-efficient than its CA-SCL-based counterparts.


title: Constructing highly regular expanders from hyperbolic Coxeter groups

abstract: A graph $X$ is defined inductively to be $(a_0,\dots,a_{n-1})$-regular if $X$ is $a_0$-regular and for every vertex $v$ of $X$, the sphere of radius $1$ around $v$ is an $(a_1,\dots,a_{n-1})$-regular graph. Such a graph $X$ is said to be highly regular (HR) of level $n$ if $a_{n-1}\neq 0$. Chapman, Linial and Peled studied HR-graphs of level 2 and provided several methods to construct families of graphs which are expanders "globally and locally". They ask whether such HR-graphs of level 3 exist. In this paper we show how the theory of Coxeter groups, and abstract regular polytopes and their generalisations, can lead to such graphs. Given a Coxeter system $(W,S)$ and a subset $M$ of $S$, we construct highly regular quotients of the 1-skeleton of the associated Wythoffian polytope $\mathcal{P}_{W,M}$, which form an infinite family of expander graphs when $(W,S)$ is indefinite and $\mathcal{P}_{W,M}$ has finite vertex links. The regularity of the graphs in this family can be deduced from the Coxeter diagram of $(W,S)$. The expansion stems from applying superapproximation to the congruence subgroups of the linear group $W$. This machinery gives a rich collection of families of HR-graphs, with various interesting properties, and in particular answers affirmatively the question asked by Chapman, Linial and Peled.


title: Bounds on Sweep-Covers by Raney Numbers

abstract: In this work, we introduce a vertex separator in trees known as a sweep-cover that is defined by an ancestor-descendent relationship with all nodes in the tree. We prove the recurrence relation of sweep-covers with $n$ subcovers $P_{\Delta, \gamma}(n)$ on a class of infinite $\Delta$-ary trees with constant path lengths $\gamma$ between the $\Delta$-star internal nodes. Then, we provide recurrence relations for Raney numbers over integer compositions and show that they provide a lower-bound for sweep-covers such that $P_{\Delta, \gamma}(n) = \Omega\left( \frac{\sqrt{2 \pi} n^{\Delta n + \Delta + \frac{3}{2}}}{e^n ((\Delta-1)n+\Delta+1)!(n+1)!} \gamma \right)$.


title: Discovery of an Ultra-Faint Stellar System near the Magellanic Clouds with the DECam Local Volume Exploration (DELVE) Survey

abstract: We report the discovery of a new ultra-faint stellar system found near the Magellanic Clouds in the DECam Local Volume Exploration (DELVE) Survey. This new system, DELVE J0155$-$6815 (DELVE 2), is located at a heliocentric distance of $D_{\odot} = 71 \pm 4\text{ kpc}$, which places it at a 3D physical separation of 12 kpc from the center of Small Magellanic Cloud (SMC) and 28 kpc from the center of the Large Magellanic Cloud (LMC). DELVE 2 is identified as a resolved overdensity of old ($\tau > 13.3\text{ Gyr}$) and metal-poor (${\rm [Fe/H]} = -2.0_{-0.5}^{+0.2}$ dex) stars with a projected half-light radius of $r_{1/2} = 21^{+4}_{-3}\text{ pc}$ and an absolute magnitude of $M_V = -2.1^{+0.4}_{-0.5}\text{ mag}$. The size and luminosity of DELVE 2 are consistent with both the population of recently discovered ultra-faint globular clusters and the smallest ultra-faint dwarf galaxies. However, its age and metallicity would place it among the oldest and most metal-poor globular clusters in the Magellanic system. DELVE 2 is detected in Gaia DR2 with a clear proper motion signal, with multiple blue horizontal branch stars near the centroid of the system with proper motions consistent with the systemic mean. We measure the system proper motion to be $(\mu_{\alpha} \cos \delta, \mu_{\delta})= (1.02_{-0.25}^{+0.24}, -0.85_{-0.19}^{+0.18})$ mas yr$^{-1}$. We compare the spatial position and proper motion of DELVE 2 with simulations of the accreted satellite population of the LMC and find that it is very likely to be associated with the LMC.


title: Kohn-Sham equations as regularizer: building prior knowledge into machine-learned physics

abstract: Including prior knowledge is important for effective machine learning models in physics, and is usually achieved by explicitly adding loss terms or constraints on model architectures. Prior knowledge embedded in the physics computation itself rarely draws attention. We show that solving the Kohn-Sham equations when training neural networks for the exchange-correlation functional provides an implicit regularization that greatly improves generalization. Two separations suffice for learning the entire one-dimensional H$_2$ dissociation curve within chemical accuracy, including the strongly correlated region. Our models also generalize to unseen types of molecules and overcome self-interaction error.


title: Structured Attention for Unsupervised Dialogue Structure Induction

abstract: Inducing a meaningful structural representation from one or a set of dialogues is a crucial but challenging task in computational linguistics. Advancement made in this area is critical for dialogue system design and discourse analysis. It can also be extended to solve grammatical inference. In this work, we propose to incorporate structured attention layers into a Variational Recurrent Neural Network (VRNN) model with discrete latent states to learn dialogue structure in an unsupervised fashion. Compared to a vanilla VRNN, structured attention enables a model to focus on different parts of the source sentence embeddings while enforcing a structural inductive bias. Experiments show that on two-party dialogue datasets, VRNN with structured attention learns semantic structures that are similar to templates used to generate this dialogue corpus. While on multi-party dialogue datasets, our model learns an interactive structure demonstrating its capability of distinguishing speakers or addresses, automatically disentangling dialogues without explicit human annotation.


title: Generation-Augmented Retrieval for Open-domain Question Answering

abstract: We propose Generation-Augmented Retrieval (GAR) for answering open-domain questions, which augments a query through text generation of heuristically discovered relevant contexts without external resources as supervision. We demonstrate that the generated contexts substantially enrich the semantics of the queries and GAR with sparse representations (BM25) achieves comparable or better performance than state-of-the-art dense retrieval methods such as DPR. We show that generating diverse contexts for a query is beneficial as fusing their results consistently yields better retrieval accuracy. Moreover, as sparse and dense representations are often complementary, GAR can be easily combined with DPR to achieve even better performance. GAR achieves state-of-the-art performance on Natural Questions and TriviaQA datasets under the extractive QA setup when equipped with an extractive reader, and consistently outperforms other retrieval methods when the same generative reader is used.


title: Long-Distance Coupling and Energy Transfer between Exciton States in Magnetically Controlled Microcavities

abstract: Coupling of quantum emitters in a semiconductor relies, generally, on short-range dipole-dipole or electronic exchange type interactions. Consistently, energy transfer between exciton states, that is, electron-hole pairs bound by Coulomb interaction, is limited to distances of the order of 10~nm. Here, we demonstrate polariton-mediated coupling and energy transfer between excitonic states over a distance exceeding 2~$\mu$m. We accomplish this by coupling quantum well-confined excitons through the delocalized mode of two coupled optical microcavities. Use of magnetically doped quantum wells enables us to tune the confined exciton energy by the magnetic field and in this way to control the spatial direction of the transfer. Such controlled, long-distance interaction between coherently coupled quantum emitters opens possibilities of a scalable implementation of quantum networks and quantum simulators based on solid-state, multi-cavity systems.


title: Improved recovery guarantees and sampling strategies for TV minimization in compressive imaging

abstract: In this paper, we consider the use of Total Variation (TV) minimization for compressive imaging; that is, image reconstruction from subsampled measurements. Focusing on two important imaging modalities -- namely, Fourier imaging and structured binary imaging via the Walsh--Hadamard transform -- we derive uniform recovery guarantees asserting stable and robust recovery for arbitrary random sampling strategies. Using this, we then derive a class of theoretically-optimal sampling strategies. For Fourier sampling, we show recovery of an image with approximately $s$-sparse gradient from $m \gtrsim_d s \cdot \log^2(s) \cdot \log^4(N)$ measurements, in $d \geq 1$ dimensions. When $d = 2$, this improves the current state-of-the-art result by a factor of $\log(s) \cdot \log(N)$. It also extends it to arbitrary dimensions $d \geq 2$. For Walsh sampling, we prove that $m \gtrsim_d s \cdot \log^2(s) \cdot \log^2(N/s) \cdot \log^3(N) $ measurements suffice in $d \geq 2$ dimensions. To the best of our knowledge, this is the first recovery guarantee for structured binary sampling with TV minimization.


title: Thermal dimensional reduction and black hole evaporation

abstract: We analyze how a quantum-gravity-induced change in the number of thermal dimensions (through a modified dispersion relation) affects the geometry and the thermodynamics of a charged black hole. To that end we resort to Kiselev's solution as the impact such modifications have on the evaporation rate of the black hole becomes more clear. As an application, we study the case for which the thermal dimension is reduced to two.


title: Ever-present Majorana bound state in a generic one-dimensional superconductor with odd number of Fermi surfaces

abstract: A quasi-1D superconductor with odd number of Fermi surfaces is expected to exhibit a nondegenerate Majorana bound state at the Fermi level at its boundary with an insulator (where the latter could be an actual insulator material or vacuum, for a terminated sample). Previous explicit theoretical demonstrations of this property were done for specific microscopic models of the bulk Hamiltonian and, most importantly, of the boundary. In this work, we theoretically demonstrate that this property holds for the whole class of systems, using the symmetry-based formalism of low-energy continuum models and general boundary conditions. We derive the general form of the Bogoliubov-de Gennes low-energy Hamiltonian that is subject only to charge-conjugation symmetry $\mathcal{C}_+$ of the type $\mathcal{C}_+^2=+1$ and a few minimal assumptions. Crucially, we also derive the most general form of the boundary conditions describing the boundary with an insulator, subject only to the fundamental principle of the probability-current conservation and $\mathcal{C}_+$ symmetry. Such {\em normal-reflection} boundary conditions do not contain scattering between electrons and holes. We find that for odd number of Fermi surfaces a Majorana bound state always exists as long as the bulk is in the gapped superconducting state, irrespective of the parameters of the bulk Hamiltonian and boundary conditions. Importantly, our general model includes a possible {\em Fermi-point mismatch}, when the two Fermi points are not at exactly opposite momenta, which disfavors superconductivity. We find that the Fermi-point mismatch does {\em not} have a direct destructive effect on the Majorana bound state, in the sense that once the bulk gap is opened the bound state is always present.


title: The Ruelle zeta function at zero for nearly hyperbolic 3-manifolds

abstract: We show that for a generic conformal metric perturbation of a compact hyperbolic 3-manifold $\Sigma$ with Betti number $b_1$, the order of vanishing of the Ruelle zeta function at zero equals $4-b_1$, while in the hyperbolic case it is equal to $4-2b_1$. This is in contrast to the 2-dimensional case where the order of vanishing is a topological invariant. The proof uses the microlocal approach to dynamical zeta functions, giving a geometric description of generalized Pollicott-Ruelle resonant differential forms at 0 in the hyperbolic case and using first variation for the perturbation. To show that the first variation is generically nonzero we introduce a new identity relating pushforwards of products of resonant and coresonant 2-forms on the sphere bundle $S\Sigma$ with harmonic 1-forms on $\Sigma$.


title: On Primes, Log-Loss Scores and (No) Privacy

abstract: Membership Inference Attacks exploit the vulnerabilities of exposing models trained on customer data to queries by an adversary. In a recently proposed implementation of an auditing tool for measuring privacy leakage from sensitive datasets, more refined aggregates like the Log-Loss scores are exposed for simulating inference attacks as well as to assess the total privacy leakage based on the adversary's predictions. In this paper, we prove that this additional information enables the adversary to infer the membership of any number of datapoints with full accuracy in a single query, causing complete membership privacy breach. Our approach obviates any attack model training or access to side knowledge with the adversary. Moreover, our algorithms are agnostic to the model under attack and hence, enable perfect membership inference even for models that do not memorize or overfit. In particular, our observations provide insight into the extent of information leakage from statistical aggregates and how they can be exploited.


title: Small but Mighty: New Benchmarks for Split and Rephrase

abstract: Split and Rephrase is a text simplification task of rewriting a complex sentence into simpler ones. As a relatively new task, it is paramount to ensure the soundness of its evaluation benchmark and metric. We find that the widely used benchmark dataset universally contains easily exploitable syntactic cues caused by its automatic generation process. Taking advantage of such cues, we show that even a simple rule-based model can perform on par with the state-of-the-art model. To remedy such limitations, we collect and release two crowdsourced benchmark datasets. We not only make sure that they contain significantly more diverse syntax, but also carefully control for their quality according to a well-defined set of criteria. While no satisfactory automatic metric exists, we apply fine-grained manual evaluation based on these criteria using crowdsourcing, showing that our datasets better represent the task and are significantly more challenging for the models.


title: Loci of the Brocard Points over Selected Triangle Families

abstract: We study the loci of the Brocard points over two selected families of triangles: (i) 2 vertices fixed on a circumference and a third one which sweeps it, (ii) Poncelet 3-periodics in the homothetic ellipse pair. Loci obtained include circles, ellipses, and teardrop-like curves. We derive expressions for both curves and their areas. We also study the locus of the vertices of Brocard triangles over the homothetic and Brocard-poristic Poncelet families.


title: Bounds for Learning Lossless Source Coding

abstract: This paper asks a basic question: how much training is required to beat a universal source coder? Traditionally, there have been two types of source coders: fixed, optimum coders such as Huffman coders; and universal source coders, such as Lempel-Ziv The paper considers a third type of source coders: learned coders. These are coders that are trained on data of a particular type, and then used to encode new data of that type. This is a type of coder that has recently become very popular for (lossy) image and video coding. The paper consider two criteria for performance of learned coders: the average performance over training data, and a guaranteed performance over all training except for some error probability $P_e$. In both cases the coders are evaluated with respect to redundancy. The paper considers the IID binary case and binary Markov chains. In both cases it is shown that the amount of training data required is very moderate: to code sequences of length $l$ the amount of training data required to beat a universal source coder is $m=K\frac{l}{\log l}$, where the constant in front depends the case considered.


title: SCREENet: A Multi-view Deep Convolutional Neural Network for Classification of High-resolution Synthetic Mammographic Screening Scans

abstract: Purpose: To develop and evaluate the accuracy of a multi-view deep learning approach to the analysis of high-resolution synthetic mammograms from digital breast tomosynthesis screening cases, and to assess the effect on accuracy of image resolution and training set size. Materials and Methods: In a retrospective study, 21,264 screening digital breast tomosynthesis (DBT) exams obtained at our institution were collected along with associated radiology reports. The 2D synthetic mammographic images from these exams, with varying resolutions and data set sizes, were used to train a multi-view deep convolutional neural network (MV-CNN) to classify screening images into BI-RADS classes (BI-RADS 0, 1 and 2) before evaluation on a held-out set of exams. Results: Area under the receiver operating characteristic curve (AUC) for BI-RADS 0 vs non-BI-RADS 0 class was 0.912 for the MV-CNN trained on the full dataset. The model obtained accuracy of 84.8%, recall of 95.9% and precision of 95.0%. This AUC value decreased when the same model was trained with 50% and 25% of images (AUC = 0.877, P=0.010 and 0.834, P=0.009 respectively). Also, the performance dropped when the same model was trained using images that were under-sampled by 1/2 and 1/4 (AUC = 0.870, P=0.011 and 0.813, P=0.009 respectively). Conclusion: This deep learning model classified high-resolution synthetic mammography scans into normal vs needing further workup using tens of thousands of high-resolution images. Smaller training data sets and lower resolution images both caused significant decrease in performance.


title: SISTA: learning optimal transport costs under sparsity constraints

abstract: In this paper, we describe a novel iterative procedure called SISTA to learn the underlying cost in optimal transport problems. SISTA is a hybrid between two classical methods, coordinate descent ("S"-inkhorn) and proximal gradient descent ("ISTA"). It alternates between a phase of exact minimization over the transport potentials and a phase of proximal gradient descent over the parameters of the transport cost. We prove that this method converges linearly, and we illustrate on simulated examples that it is significantly faster than both coordinate descent and ISTA. We apply it to estimating a model of migration, which predicts the flow of migrants using country-specific characteristics and pairwise measures of dissimilarity between countries. This application demonstrates the effectiveness of machine learning in quantitative social sciences.


title: Almost certain loss from black holes: critical comments on the black hole final state proposal

abstract: In this paper, we critically revisit the Horowitz-Maldacena proposal and its generalization by Lloyd. In the original proposal, as well as in Lloyd's generalization, Hawking radiation involves a pair of maximally entangled quantum states in which the ingoing partner state and the collapsed matter form either a maximally entangled pair or a Schmidt decomposed random state near the singularity. We point out that the unitary matrix introduced in Lloyd's fidelity calculation depends on initial matter states; hence, his result on the high average fidelity may not represent an almost unitary evolution. In opposition to Lloyd's conclusion, when we do not include the state-dependent unitary matrix for the fidelity computation, we analytically and numerically confirm that information will almost certainly be lost because the fidelity will approach zero as the degrees of freedom increase.


title: MUTANT: A Training Paradigm for Out-of-Distribution Generalization in Visual Question Answering

abstract: While progress has been made on the visual question answering leaderboards, models often utilize spurious correlations and priors in datasets under the i.i.d. setting. As such, evaluation on out-of-distribution (OOD) test samples has emerged as a proxy for generalization. In this paper, we present MUTANT, a training paradigm that exposes the model to perceptually similar, yet semantically distinct mutations of the input, to improve OOD generalization, such as the VQA-CP challenge. Under this paradigm, models utilize a consistency-constrained training objective to understand the effect of semantic changes in input (question-image pair) on the output (answer). Unlike existing methods on VQA-CP, MUTANT does not rely on the knowledge about the nature of train and test answer distributions. MUTANT establishes a new state-of-the-art accuracy on VQA-CP with a $10.57\%$ improvement. Our work opens up avenues for the use of semantic input mutations for OOD generalization in question answering.


title: The Multipolar Magnetic Field of Millisecond Pulsar PSR J0030+0451

abstract: Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. (2019) and Riley et al. (2019), we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions to the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous gamma-ray light curves following Kalapotharakos et al. (2014) these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the gamma-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and gamma-ray light curves. We conclude with a discussion of the broader implications of our study.


title: Inference for Large-Scale Linear Systems with Known Coefficients

abstract: This paper considers the problem of testing whether there exists a non-negative solution to a possibly under-determined system of linear equations with known coefficients. This hypothesis testing problem arises naturally in a number of settings, including random coefficient, treatment effect, and discrete choice models, as well as a class of linear programming problems. As a first contribution, we obtain a novel geometric characterization of the null hypothesis in terms of identified parameters satisfying an infinite set of inequality restrictions. Using this characterization, we devise a test that requires solving only linear programs for its implementation, and thus remains computationally feasible in the high-dimensional applications that motivate our analysis. The asymptotic size of the proposed test is shown to equal at most the nominal level uniformly over a large class of distributions that permits the number of linear equations to grow with the sample size.