title: Electronic structure of water from Koopmans-compliant functionals

abstract: Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW+f$_{xc}$). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that the different MD approaches lead to up to 1 eV change in the average band gap, thus, we focused on the band gap dependence with the geometrical properties of the system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance, and HOH angle, as well as the intermolecular hydrogen bond O$\cdot\cdot\cdot$O distance, and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O$\cdot\cdot\cdot$O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows a good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, a broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.


title: Deep Learning for direct Dark Matter search with nuclear emulsions

abstract: We propose a new method for the discrimination of sub-micron nuclear recoil tracks from an instrumental background in fine-grain nuclear emulsions used in the directional dark matter search. The proposed method uses a 3D Convolutional Neural Network, whose parameters are optimised by Bayesian search. Unlike previous studies focused on extracting the directional information, we focus on the signal/background separation exploiting the polarisation dependence of the Localised Surface Plasmon Resonance phenomenon. Comparing the proposed method with the conventional cut-based approach shows a significant boost in the reduction factor for given signal efficiency.


title: Introducing endogenous transport provision in a LUTI model to explore polycentric governance systems

abstract: Models focusing on interactions between land-use and transportation mostly assume an exogenous provision of transportation infrastructures. We investigate here co-evolutionary processes between land-use and transportation, at the scale of Mega-City Regions, by introducing a toy model of corresponding processes. In particular, our model is specifically tailored to include governance processes ruling the growth of transportation infrastructure. We show through stylised numerical simulations the potentialities of our model to reproduce a variety of dynamics when co-evolution is taken into account. We then apply the model to a case study, by calibrating it for the Pearl River Delta Mega-city Region (China, 1990-2010). To go beyond this first modelling step, we elaborate on the challenges to overcome to go further towards more complex models integrating co-evolution between transportation networks and territories in urban systems.


title: Dissipative quasi-particle picture for quadratic Markovian open quantum systems

abstract: Correlations between different regions of a quantum many-body system can be quantified through measures based on entropies of (reduced) subsystem states. For closed systems, several analytical and numerical tools, e.g., hydrodynamic theories or tensor networks, can accurately capture the time-evolution of subsystem entropies, thus allowing for a profound understanding of the unitary dynamics of quantum correlations. However, so far, these methods either cannot be applied to open quantum systems or do not permit an efficient computation of quantum entropies for mixed states. Here, we make progress in solving this issue by formulating a dissipative quasi-particle picture -- describing the dynamics of quantum entropies in the hydrodynamic limit -- for a general class of noninteracting open quantum systems. Our results show that also in dissipative many-body systems, correlations are generically established through the propagation of quasi-particles.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el quant-ph


title: Towards an independent determination of muon g-2 from muonium spectroscopy

abstract: We show that muonium spectroscopy in the coming years can reach a precision high enough to determine the anomalous magnetic moment of the muon below one part per million (ppm). Such an independent determination of muon g-2, which is not limited by hadronic uncertainties, would certainly shed light on the 2ppm difference currently observed between spin-precession measurements and (R-ratio based) Standard Model predictions. The magnetic dipole interaction between electrons and (anti)muons bound in muonium gives rise to a hyperfine splitting (HFS) of the ground state which is sensitive to the muon anomalous magnetic moment. A direct comparison of the muonium frequency measurements of the HFS at J-PARC and the 1S-2S transition at PSI with theory predictions will allow to extract muon g-2 with high precision. Improving the accuracy of QED calculations of these transitions by about one order of magnitude is also required. Moreover, the good agreement between theory and experiment for the electron g-2 indicates that new physics interactions are unlikely to affect muonium spectroscopy down to the envisaged precision.


title: On Factorizable S-matrices, Generalized TTbar, and the Hagedorn Transition

abstract: We study solutions of the Thermodynamic Bethe Ansatz equations for relativistic theories defined by the factorizable $S$-matrix of an integrable QFT deformed by CDD factors. Such $S$-matrices appear under generalized TTbar deformations of integrable QFT by special irrelevant operators. The TBA equations, of course, determine the ground state energy $E(R)$ of the finite-size system, with the spatial coordinate compactified on a circle of circumference $R$. We limit attention to theories involving just one kind of stable particles, and consider deformations of the trivial (free fermion or boson) $S$-matrix by CDD factors with two elementary poles and regular high energy asymptotics -- the "2CDD model". We find that for all values of the parameters (positions of the CDD poles) the TBA equations exhibit two real solutions at $R$ greater than a certain parameter-dependent value $R_*$, which we refer to as the primary and secondary branches. The primary branch is identified with the standard iterative solution, while the secondary one is unstable against iterations and needs to be accessed through an alternative numerical method known as pseudo-arc-length continuation. The two branches merge at the "turning point" $R_*$ (a square-root branching point). The singularity signals a Hagedorn behavior of the density of high energy states of the deformed theories, a feature incompatible with the Wilsonian notion of a local QFT originating from a UV fixed point, but typical for string theories. This behavior of $E(R)$ is qualitatively the same as the one for standard TTbar deformations of local QFT.

categories: hep-th cond-mat.stat-mech math-ph math.MP nlin.SI


title: Absence of precursor incommensurate charge order in electronic nematic Ba0.35Sr0.65Ni2As2

abstract: Recent discoveries of charge order and electronic nematic order in the iron-based superconductors and cuprates have pointed towards the possibility of nematic and charge fluctuations playing a role in the enhancement of superconductivity. The Ba1-xSrxNi2As2 system, closely related in structure to the BaFe2As2 system, has recently been shown to exhibit both types of ordering without the presence of any magnetic order. We report single crystal X-ray diffraction, resistance transport measurements, and magnetization of \BaSrLate, providing evidence that the previously reported incommensurate charge order with wavevector $(0,0.28,0)_{tet}$ in the tetragonal state of \BaNi~vanishes by this concentration of Sr substitution together with nematic order. Our measurements suggest that the nematic and incommensurate charge orders are closely tied in the tetragonal state, and show that the $(0,0.33,0)_{tri}$ charge ordering in the triclinic phase of BaNi2As2 evolves to become $(0,0.5,0)_{tri}$ charge ordering at $x$=0.65 before vanishing at $x$=0.71.


title: Inference in High-dimensional Linear Regression

abstract: This paper develops an approach to inference in a linear regression model when the number of potential explanatory variables is larger than the sample size. The approach treats each regression coefficient in turn as the interest parameter, the remaining coefficients being nuisance parameters, and seeks an optimal interest-respecting transformation, inducing sparsity on the relevant blocks of the notional Fisher information matrix. The induced sparsity is exploited through a marginal least squares analysis for each variable, as in a factorial experiment, thereby avoiding penalization. One parameterization of the problem is found to be particularly convenient, both computationally and mathematically. In particular, it permits an analytic solution to the optimal transformation problem, facilitating theoretical analysis and comparison to other work. In contrast to regularized regression such as the lasso and its extensions, neither adjustment for selection nor rescaling of the explanatory variables is needed, ensuring the physical interpretation of regression coefficients is retained. Recommended usage is within a broader set of inferential statements, so as to reflect uncertainty over the model as well as over the parameters. The considerations involved in extending the work to other regression models are briefly discussed.


title: Coordinates for non-integrable Lie algebroids

abstract: We construct local coordinates for the Weinstein groupoid of a non-integrable Lie algebroid. To this end, we reformulate the notion of bi-submersion in a completely algebraic way and prove the existence of bi-submersions as such for the Weinstein groupoid. This implies that a $C^*$-algebra can be attached to every Lie algebroid.


title: Global 3D Radiation Hydrodynamic Simulations of Proto-Jupiter's Convective Envelope

abstract: The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter's envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated spherical envelopes and envelopes embedded in discs. The envelopes are heated at given rates to achieve steady states, enabling comparisons with 1D models. We vary envelope properties to obtain both radiative and convective solutions. Using a passive scalar, we observe significant mass recycling on the orbital timescale. For a radiative envelope, recycling can only penetrate from the disc surface until $\sim$0.1-0.2 planetary Hill radii, while for a convective envelope, the convective motion can "dredge up" the deeper part of the envelope so that the entire convective envelope is recycled efficiently. This recycling, however, has only limited effects on the envelopes' thermal structure. The radiative envelope embedded in the disc has identical structure as the isolated envelope. The convective envelope has a slightly higher density when it is embedded in the disc. We introduce a modified 1D approach which can fully reproduce our 3D simulations. With our updated opacity and 1D model, we recompute Jupiter's envelope accretion with a 10 $M_{\oplus}$ core, and the timescale to runaway accretion is shorter than the disc lifetime as in prior studies. Finally, we discuss the implications of the efficient recycling on the observed chemical abundances of the planetary atmosphere (especially for super-Earths and mini-Neptunes).


title: Modelling lipid-coated microbubbles in focused ultrasound applications at subresonance frequencies

abstract: We present a computational study of the behaviour of a lipid-coated SonoVue microbubble with initial radius $1 \, \mu \text{m} \leq R_0 \leq 2 \, \mu \text{m}$, excited at frequencies (200-1500 kHz) significantly below the linear resonance frequency and pressure amplitudes of up to 1500 kPa, an excitation regime used in many focused ultrasound applications. The bubble dynamics are simulated using the Rayleigh-Plesset equation and the Gilmore equation, in conjunction with the Marmottant model for the lipid monolayer coating. Also, a new continuously differentiable variant of the Marmottant model is introduced. Below the onset of inertial cavitation, a linear regime is identified in which the maximum pressure at the bubble wall is linearly proportional to the excitation pressure amplitude and, likewise, the mechanical index. This linear regime is bounded by the Blake pressure and, in line with recent in vitro experiments, the onset of inertial cavitation is found to occur approximately at an excitation pressure amplitude of 130-190 kPa, dependent on the initial bubble size. In the nonlinear regime the maximum pressure at the bubble wall is found to be readily predicted by the maximum bubble radius and both the Rayleigh-Plesset and Gilmore equations are shown to predict the onset of sub- and ultraharmonic frequencies of the acoustic emissions compared to in vitro experiments. Neither the surface dilatational viscosity of the lipid monolayer nor the compressibility of the liquid have a discernible influence on the studied quantities, yet accounting for the lipid coating is critical for the accurate prediction of the bubble behaviour. The Gilmore equation is shown to be valid for the considered bubbles and excitation regime, and the Rayleigh-Plesset equation also provides accurate qualitative predictions, even though it is outside its range of validity for many of the considered cases.


title: Exploring the Representational Power of Graph Autoencoder

abstract: While representation learning has yielded a great success on many graph learning tasks, there is little understanding behind the structures that are being captured by these embeddings. For example, we wonder if the topological features, such as the Triangle Count, the Degree of the node and other centrality measures are concretely encoded in the embeddings. Furthermore, we ask if the presence of these structures in the embeddings is necessary for a better performance on the downstream tasks, such as clustering and classification. To address these questions, we conduct an extensive empirical study over three classes of unsupervised graph embedding models and seven different variants of Graph Autoencoders. Our results show that five topological features: the Degree, the Local Clustering Score, the Betweenness Centrality, the Eigenvector Centrality, and Triangle Count are concretely preserved in the first layer of the graph autoencoder that employs the SUM aggregation rule, under the condition that the model preserves the second-order proximity. We supplement further evidence for the presence of these features by revealing a hierarchy in the distribution of the topological features in the embeddings of the aforementioned model. We also show that a model with such properties can outperform other models on certain downstream tasks, especially when the preserved features are relevant to the task at hand. Finally, we evaluate the suitability of our findings through a test case study related to social influence prediction.


title: Conformal TBA for resolved conifolds

abstract: We revisit the Riemann-Hilbert problem determined by Donaldson-Thomas invariants for the resolved conifold and for other small crepant resolutions. While this problem can be recast as a system of TBA-type equations in the conformal limit, solutions are ill-defined due to divergences in the sum over infinite trajectories in the spectrum of D2-D0-brane bound states. We explore various prescriptions to make the sum well-defined, show that one of them reproduces the existing solution in the literature, and identify an alternative solution which is better behaved in a certain limit. Furthermore, we show that a suitable asymptotic expansion of the $\tau$ function reproduces the genus expansion of the topological string partition function for any small crepant resolution. As a by-product, we conjecture new integral representations for the triple sine function, similar to Woronowicz' integral representation for Faddeev's quantum dilogarithm.


title: Energy hardware and workload aware job scheduling towards interconnected HPC environments

abstract: New HPC machines are getting close to the exascale. Power consumption for those machines has been increasing, and researchers are studying ways to reduce it. A second trend is HPC machines' growing complexity, with increasing heterogeneous hardware components and different clusters architectures cooperating in the same machine. We refer to these environments with the term heterogeneous multi-cluster environments. With the aim of optimizing performance and energy consumption in these environments, this paper proposes an Energy-Aware-Multi-Cluster (EAMC) job scheduling policy. EAMC-policy is able to optimize the scheduling and placement of jobs by predicting performance and energy consumption of arriving jobs for different hardware architectures and processor frequencies, reducing workload's energy consumption, makespan, and response time. The policy assigns a different priority to each job-resource combination so that the most efficient ones are favored, while less efficient ones are still considered on a variable degree, reducing response time and increasing cluster utilization. We implemented EAMC-policy in Slurm, and we evaluated a scenario in which two CPU clusters collaborate in the same machine. Simulations of workloads running applications modeled from real-world show a reduction of response time and makespan by up to 25% and 6% while saving up to 20% of total energy consumed when compared to policies minimizing runtime, and by 49%, 26%, and 6% compared to policies minimizing energy.


title: Transparent planar solar absorber for winter thermal management

abstract: Indoor heating during winters accounts for a significant portion of energy consumed by buildings in regions of cold climate. Development of transparent coatings for windows that efficiently harvest solar energy can play a major role in reducing energy consumption and fuel costs incurred for winter heating. In recent years, there has been a great research effort towards designing transparent solar absorber coatings using nanophotonic structures. The potential of coatings based on planar multilayer structures, however, has received very little attention. In this work we investigate the performance of planar multilayer thin films using low cost materials for design of transparent solar absorber window coatings. Our study led to the proposal of two planar multilayer designs. Simulation results predict that an increase in surface temperature by 21 K and 25 K, while maintaining mean visible transmittance of over 60% is possible using these designs. These results illustrate the great promise planar multilayer structures hold for winter thermal management of buildings.


title: ALMA CO Observations of the Mixed-Morphology Supernova Remnant W49B: Efficient Production of Recombining Plasma and Hadronic Gamma-rays via Shock-Cloud Interactions

abstract: We carried out new CO($J$ = 2-1) observations toward the mixed-morphology supernova remnant (SNR) W49B with the Atacama Large Millimeter/submillimeter Array (ALMA). We found that CO clouds at $\sim$10 km s$^{-1}$ show a good spatial correspondence with synchrotron radio continuum as well as an X-ray deformed shell. The bulk mass of molecular clouds accounts for the western part of the shell, not for the eastern shell where near-infrared H$_2$ emission is detected. The molecular clouds at $\sim$10 km s$^{-1}$ show higher kinetic temperature of $\sim$20-60 K, suggesting that modest shock-heating occurred. The expanding motion of the clouds with $\Delta V \sim$6 km s$^{-1}$ was formed by strong winds from the progenitor system. We argue that the barrel-like structure of Fe rich ejecta was possibly formed not only by an asymmetric explosion, but also by interactions with dense molecular clouds. We also found a negative correlation between the CO intensity and the electron temperature of recombining plasma, implying that the origin of the high-temperature recombining plasma in W49B can be understood as the thermal conduction model. The total energy of accelerated cosmic-ray protons $W_\mathrm{p}$ is estimated to be $\sim$$2\times 10^{49}$ erg by adopting an averaged gas density of $\sim$$650\pm200$ cm$^{-3}$. The SNR age-$W_\mathrm{p}$ diagram indicates that W49B shows one of the highest in-situ values of $W_\mathrm{p}$ in the gamma-ray bright SNRs.


title: Space-time Non-local multi-continua upscaling for parabolic equations with moving channelized media

abstract: In this paper, we consider a parabolic problem with time-dependent heterogeneous coefficients. Many applied problems have coupled space and time heterogeneities. Their homogenization or upscaling requires cell problems that are formulated in space-time representative volumes for problems with scale separation. In problems without scale separation, local problems include multiple macroscopic variables and oversampled local problems, where these macroscopic parameters are computed. These approaches, called Non-local multi-continua, are proposed for problems with complex spatial heterogeneities in a number of previous papers. In this paper, we extend this approach for space-time heterogeneities, by identifying macroscopic parameters in space-time regions. Our proposed method space-time Non-local multi-continua (space-time NLMC) is an efficient numerical solver to deal with time-dependent heterogeneous coefficients. It provides a flexible and systematic way to construct multiscale basis functions to approximate the solution. These multiscale basis functions are constructed by solving a local energy minimization problems in the oversampled space-time regions such that these multiscale basis functions decay exponentially outside the oversampled domain. Unlike the classical time-stepping methods combined with full-discretization technique, our space-time NLMC efficiently constructs the multiscale basis functions in a space-time domain and can provide a computational savings compared to space-only approaches as we discuss in the paper. We present two numerical experiments, which show that the proposed approach can provide a good accuracy.


title: P2T: Pyramid Pooling Transformer for Scene Understanding

abstract: Recently, the vision transformer has achieved great success by pushing the state-of-the-art of various vision tasks. One of the most challenging problems in the vision transformer is that the large sequence length of image tokens leads to high computational cost (quadratic complexity). A popular solution to this problem is to use a single pooling operation to reduce the sequence length. This paper considers how to improve existing vision transformers, where the pooled feature extracted by a single pooling operation seems less powerful. To this end, we note that pyramid pooling has been demonstrated to be effective in various vision tasks owing to its powerful ability in context abstraction. However, pyramid pooling has not been explored in backbone network design. To bridge this gap, we propose to adapt pyramid pooling to Multi-Head Self-Attention (MHSA) in the vision transformer, simultaneously reducing the sequence length and capturing powerful contextual features. Plugged with our pooling-based MHSA, we build a universal vision transformer backbone, dubbed Pyramid Pooling Transformer (P2T). Extensive experiments demonstrate that, when applied P2T as the backbone network, it shows substantial superiority in various vision tasks such as image classification, semantic segmentation, object detection, and instance segmentation, compared to previous CNN- and transformer-based networks. The code will be released at https://github.com/yuhuan-wu/P2T.


title: Test-time Collective Prediction

abstract: An increasingly common setting in machine learning involves multiple parties, each with their own data, who want to jointly make predictions on future test points. Agents wish to benefit from the collective expertise of the full set of agents to make better predictions than they would individually, but may not be willing to release their data or model parameters. In this work, we explore a decentralized mechanism to make collective predictions at test time, leveraging each agent's pre-trained model without relying on external validation, model retraining, or data pooling. Our approach takes inspiration from the literature in social science on human consensus-making. We analyze our mechanism theoretically, showing that it converges to inverse meansquared-error (MSE) weighting in the large-sample limit. To compute error bars on the collective predictions we propose a decentralized Jackknife procedure that evaluates the sensitivity of our mechanism to a single agent's prediction. Empirically, we demonstrate that our scheme effectively combines models with differing quality across the input space. The proposed consensus prediction achieves significant gains over classical model averaging, and even outperforms weighted averaging schemes that have access to additional validation data.


title: The ASKAP-EMU Early Science Project: 888 MHz Radio Continuum Survey of the Large Magellanic Cloud

abstract: We present an analysis of a new 120 deg$^{2}$ radio continuum image of the Large Magellanic Cloud (LMC) at 888 MHz with a bandwidth of 288 MHz and beam size of $13\rlap{.}^{\prime\prime}9\times12\rlap{.}^{\prime\prime}1$, from the Australian Square Kilometre Array Pathfinder (ASKAP) processed as part of the Evolutionary Map of the Universe (EMU) survey. The median Root Mean Squared noise is 58 $\mu$Jy beam$^{-1}$. We present a catalogue of 54,612 sources, divided over a GOLD list (30,866 sources) complete down to 0.5 mJy uniformly across the field, a SILVER list (22,080 sources) reaching down to $<$ 0.2 mJy and a BRONZE list (1,666 sources) of visually inspected sources in areas of high noise and/or near bright complex emission. We discuss detections of planetary nebulae and their radio luminosity function, young stellar objects showing a correlation between radio luminosity and gas temperature, novae and X-ray binaries in the LMC, and active stars in the Galactic foreground that may become a significant population below this flux level. We present examples of diffuse emission in the LMC (H II regions, supernova remnants, bubbles) and distant galaxies showcasing spectacular interaction between jets and intracluster medium. Among 14,333 infrared counterparts of the predominantly background radio source population we find that star-forming galaxies become more prominent below 3 mJy compared to active galactic nuclei. We combine the new 888 MHz data with archival Australia Telescope Compact Array data at 1.4 GHz to determine spectral indices; the vast majority display synchrotron emission but flatter spectra occur too. We argue that the most extreme spectral index values are due to variability.


title: Transfer Learning of Deep Spatiotemporal Networks to Model Arbitrarily Long Videos of Seizures

abstract: Detailed analysis of seizure semiology, the symptoms and signs which occur during a seizure, is critical for management of epilepsy patients. Inter-rater reliability using qualitative visual analysis is often poor for semiological features. Therefore, automatic and quantitative analysis of video-recorded seizures is needed for objective assessment. We present GESTURES, a novel architecture combining convolutional neural networks (CNNs) and recurrent neural networks (RNNs) to learn deep representations of arbitrarily long videos of epileptic seizures. We use a spatiotemporal CNN (STCNN) pre-trained on large human action recognition (HAR) datasets to extract features from short snippets (approx. 0.5 s) sampled from seizure videos. We then train an RNN to learn seizure-level representations from the sequence of features. We curated a dataset of seizure videos from 68 patients and evaluated GESTURES on its ability to classify seizures into focal onset seizures (FOSs) (N = 106) vs. focal to bilateral tonic-clonic seizures (TCSs) (N = 77), obtaining an accuracy of 98.9% using bidirectional long short-term memory (BLSTM) units. We demonstrate that an STCNN trained on a HAR dataset can be used in combination with an RNN to accurately represent arbitrarily long videos of seizures. GESTURES can provide accurate seizure classification by modeling sequences of semiologies.


title: Lattice points problem, equidistribution and ergodic theorems for certain arithmetic spheres

abstract: We establish an asymptotic formula for the number of lattice points in the sets \[ \mathbf S_{h_1, h_2, h_3}(\lambda): =\{x\in\mathbb Z_+^3:\lfloor h_1(x_1)\rfloor+\lfloor h_2(x_2)\rfloor+\lfloor h_3(x_3)\rfloor=\lambda\} \quad \text{with}\quad \lambda\in\mathbb Z_+; \] where functions $h_1, h_2, h_3$ are constant multiples of regularly varying functions of the form $h(x):=x^c\ell_h(x)$, where the exponent $c>1$ (but close to $1$) and a function $\ell_h(x)$ is taken from a certain wide class of slowly varying functions. Taking $h_1(x)=h_2(x)=h_3(x)=x^c$ we will also derive an asymptotic formula for the number of lattice points in the sets \[ \mathbf S_{c}^3(\lambda) := \{x \in \mathbb Z^3 : \lfloor |x_1|^c \rfloor + \lfloor |x_2|^c \rfloor + \lfloor |x_3|^c \rfloor= \lambda \} \quad \text{with}\quad \lambda\in\mathbb Z_+; \] which can be thought of as a perturbation of the classical Waring problem in three variables. We will use the latter asymptotic formula to study, the main results of this paper, norm and pointwise convergence of the ergodic averages \[ \frac{1}{\#\mathbf S_{c}^3(\lambda)}\sum_{n\in \mathbf S_{c}^3(\lambda)}f(T_1^{n_1}T_2^{n_2}T_3^{n_3}x) \quad \text{as}\quad \lambda\to\infty; \] where $T_1, T_2, T_3:X\to X$ are commuting invertible and measure-preserving transformations of a $\sigma$-finite measure space $(X, \nu)$ for any function $f\in L^p(X)$ with $p>\frac{11-4c}{11-7c}$. Finally, we will study the equidistribution problem corresponding to the spheres $\mathbf S_{c}^3(\lambda)$.


title: On Matrix Factorizations in Subspace Clustering

abstract: This article explores subspace clustering algorithms using CUR decompositions, and examines the effect of various hyperparameters in these algorithms on clustering performance on two real-world benchmark datasets, the Hopkins155 motion segmentation dataset and the Yale face dataset. Extensive experiments are done for a variety of sampling methods and oversampling parameters for these datasets, and some guidelines for parameter choices are given for practical applications.


title: Plurisupported Currents on Compact K\"ahler Manifolds

abstract: Let $X$ be a compact K\"ahler manifold. We study plurisupported currents on $X$, i.e. closed, positive $(1,1)$-currents which are supported on a pluripolar set. In particular, we are able present a technical generalization of Witt-Nystr\"om's proof of the BDPP conjecture on projective manifolds, showing that this conjecture holds on $X$ admitting at least one plurisupported current $T$ such that $[T]$ is K\"ahler. One of the steps in our proof is to show an upper-bound for the pluripolar mass of certain envelopes of quasi-psh functions when the cohomology class is shifted, a result of independent interest. Using this, we are able to generalize an inequality of McKinnon and Roth to arbitrary pseudoeffective classes on compact K\"ahler manifolds.


title: Back to quiescence: post-outburst evolution of the pulsar J1119-6127 and its wind nebula

abstract: We report on the analysis of a deep Chandra observation of the high-magnetic field pulsar (PSR) J1119-6127 and its compact pulsar wind nebula (PWN) taken in October 2019, three years after the source went into outburst. The 0.5-7 keV post-outburst (2019) spectrum of the pulsar is best described by a two-component blackbody plus powerlaw model with a temperature of 0.2\pm0.1 keV, photon index of 1.8\pm0.4 and X-ray luminosity of ~1.9e33 erg s^{-1}, consistent with its pre-burst quiescent phase. We find that the pulsar has gone back to quiescence. The compact nebula shows a jet-like morphology elongated in the north-south direction, similar to the pre-burst phase. The post-outburst PWN spectrum is best fit by an absorbed powerlaw with a photon index of 2.3\pm0.5 and flux of ~3.2e-14 erg cm^{-2} s^{-1} (0.5-7 keV). The PWN spectrum shows evidence of spectral softening in the post-outburst phase, with the pre-burst photon index of 1.2\pm0.4 changing to 2.3\pm0.5, and pre-burst luminosity of ~1.5e32 erg s^{-1} changing to 2.7e32 erg s^{-1} in the 0.5-7 keV band, suggesting magnetar outbursts can impact PWNe. The observed timescale for returning to quiescence, of just a few years, implies a rather fast cooling process and favors a scenario where J1119 is temporarily powered by magnetic energy following the magnetar outburst, in addition to its spin-down energy.


title: Length-Preserving Directions and Some Diophantine Equations

abstract: We study directions along which the norms of vectors are preserved under a linear map. In particular, we find families of matrices for which these directions are determined by integer vectors. We consider the two-dimensional case in detail, and also discuss the extension to the three-dimensional case.


title: Flutter instability and Ziegler destabilization paradox for elastic rods subject to non-holonomic constraints

abstract: Two types of non-holonomic constraints (imposing a prescription on velocity) are analyzed, connected to an end of a (visco)elastic rod, straight in its undeformed configuration. The equations governing the nonlinear dynamics are obtained and then linearized near the trivial equilibrium configuration. The two constraints are shown to lead to the same equations governing the linearized dynamics of the Beck (or Pfluger) column in one case and of the Reut column in the other. Therefore, although the structural systems are fully conservative (when viscosity is set to zero), they exhibit flutter and divergence instability. In addition, the Ziegler's destabilization paradox is found when dissipation sources are introduced. It follows that these features are proven to be not only a consequence of 'unrealistic non-conservative loads' (as often stated in the literature), rather, the models proposed by Beck, Reut, and Ziegler can exactly describe the linearized dynamics of structures subject to non-holonomic constraints, which are made now fully accessible to experiments.


title: DetectX -- Adversarial Input Detection using Current Signatures in Memristive XBar Arrays

abstract: Adversarial input detection has emerged as a prominent technique to harden Deep Neural Networks(DNNs) against adversarial attacks. Most prior works use neural network-based detectors or complex statistical analysis for adversarial detection. These approaches are computationally intensive and vulnerable to adversarial attacks. To this end, we propose DetectX - a hardware friendly adversarial detection mechanism using hardware signatures like Sum of column Currents (SoI) in memristive crossbars (XBar). We show that adversarial inputs have higher SoI compared to clean inputs. However, the difference is too small for reliable adversarial detection. Hence, we propose a dual-phase training methodology: Phase1 training is geared towards increasing the separation between clean and adversarial SoIs; Phase2 training improves the overall robustness against different strengths of adversarial attacks. For hardware-based adversarial detection, we implement the DetectX module using 32nm CMOS circuits and integrate it with a Neurosim-like analog crossbar architecture. We perform hardware evaluation of the Neurosim+DetectX system on the Neurosim platform using datasets-CIFAR10(VGG8), CIFAR100(VGG16) and TinyImagenet(ResNet18). Our experiments show that DetectX is 10x-25x more energy efficient and immune to dynamic adversarial attacks compared to previous state-of-the-art works. Moreover, we achieve high detection performance (ROC-AUC > 0.95) for strong white-box and black-box attacks. The code has been released at https://github.com/Intelligent-Computing-Lab-Yale/DetectX


title: Promoter methylation in a mixed feedback loop circadian clock model

abstract: We introduce and analyze an extension of the mixed feedback loop model of Fran\c{c}ois and Hakim. Our extension includes an additional promoter state and allows for reversible protein sequestration, which was absent from the original studies of the mixed feedback loop model. Motivated by experimental observations that link DNA methylation with circadian gene expression, we use our extended model to investigate the role of DNA methylation(s) in the mammalian circadian clock. We extend the perturbation analysis of Fran\c{c}ois and Hakim to determine how methylation affects the presence and the periodicity of oscillations. We derive a modified Goodwin oscillator model as an approximation to show that although methylation contributes to period control, excessive methylation can abolish rhythmicity.


title: Team PyKale (xy9) Submission to the EPIC-Kitchens 2021 Unsupervised Domain Adaptation Challenge for Action Recognition

abstract: This report describes the technical details of our submission to the EPIC-Kitchens 2021 Unsupervised Domain Adaptation Challenge for Action Recognition. The EPIC-Kitchens dataset is more difficult than other video domain adaptation datasets due to multi-tasks with more modalities. Firstly, to participate in the challenge, we employ a transformer to capture the spatial information from each modality. Secondly, we employ a temporal attention module to model temporal-wise inter-dependency. Thirdly, we employ the adversarial domain adaptation network to learn the general features between labeled source and unlabeled target domain. Finally, we incorporate multiple modalities to improve the performance by a three-stream network with late fusion. Our network achieves the comparable performance with the state-of-the-art baseline T$A^3$N and outperforms the baseline on top-1 accuracy for verb class and top-5 accuracies for all three tasks which are verb, noun and action. Under the team name xy9, our submission achieved 5th place in terms of top-1 accuracy for verb class and all top-5 accuracies.


title: Q-Learning Lagrange Policies for Multi-Action Restless Bandits

abstract: Multi-action restless multi-armed bandits (RMABs) are a powerful framework for constrained resource allocation in which $N$ independent processes are managed. However, previous work only study the offline setting where problem dynamics are known. We address this restrictive assumption, designing the first algorithms for learning good policies for Multi-action RMABs online using combinations of Lagrangian relaxation and Q-learning. Our first approach, MAIQL, extends a method for Q-learning the Whittle index in binary-action RMABs to the multi-action setting. We derive a generalized update rule and convergence proof and establish that, under standard assumptions, MAIQL converges to the asymptotically optimal multi-action RMAB policy as $t\rightarrow{}\infty$. However, MAIQL relies on learning Q-functions and indexes on two timescales which leads to slow convergence and requires problem structure to perform well. Thus, we design a second algorithm, LPQL, which learns the well-performing and more general Lagrange policy for multi-action RMABs by learning to minimize the Lagrange bound through a variant of Q-learning. To ensure fast convergence, we take an approximation strategy that enables learning on a single timescale, then give a guarantee relating the approximation's precision to an upper bound of LPQL's return as $t\rightarrow{}\infty$. Finally, we show that our approaches always outperform baselines across multiple settings, including one derived from real-world medication adherence data.


title: Rogue waves on an elliptic function background in complex modified Korteweg-de Vries equation

abstract: With the assistance of one fold Darboux transformation formula, we derive rogue wave solutions of the complex modified Korteweg-de Vries equation on an elliptic function background. We employ an algebraic method to find the necessary squared eigenfunctions and eigenvalues. To begin we construct the elliptic function background. Then, on top of this background, we create a rogue wave. We demonstrate the outcome for three distinct elliptic modulus values. We find that when we increase the modulus value the amplitude of rogue waves on the dn-periodic background decreases whereas it increases in the case of cn-periodic background.


title: The Neurally-Guided Shape Parser: Grammar-based Labeling of 3D Shape Regions with Approximate Inference

abstract: We propose the Neurally-Guided Shape Parser (NGSP), a method that learns how to assign fine-grained semantic labels to regions of a 3D shape. NGSP solves this problem via MAP inference, modeling the posterior probability of a label assignment conditioned on an input shape with a learned likelihood function. To make this search tractable, NGSP employs a neural guide network that learns to approximate the posterior. NGSP finds high-probability label assignments by first sampling proposals with the guide network and then evaluating each proposal under the full likelihood. We evaluate NGSP on the task of fine-grained semantic segmentation of manufactured 3D shapes from PartNet, where shapes have been decomposed into regions that correspond to part instance over-segmentations. We find that NGSP delivers significant performance improvements over comparison methods that (i) use regions to group per-point predictions, (ii) use regions as a self-supervisory signal or (iii) assign labels to regions under alternative formulations. Further, we show that NGSP maintains strong performance even with limited labeled data or noisy input shape regions. Finally, we demonstrate that NGSP can be directly applied to CAD shapes found in online repositories and validate its effectiveness with a perceptual study.


title: ABCD: A Graph Framework to Convert Complex Sentences to a Covering Set of Simple Sentences

abstract: Atomic clauses are fundamental text units for understanding complex sentences. Identifying the atomic sentences within complex sentences is important for applications such as summarization, argument mining, discourse analysis, discourse parsing, and question answering. Previous work mainly relies on rule-based methods dependent on parsing. We propose a new task to decompose each complex sentence into simple sentences derived from the tensed clauses in the source, and a novel problem formulation as a graph edit task. Our neural model learns to Accept, Break, Copy or Drop elements of a graph that combines word adjacency and grammatical dependencies. The full processing pipeline includes modules for graph construction, graph editing, and sentence generation from the output graph. We introduce DeSSE, a new dataset designed to train and evaluate complex sentence decomposition, and MinWiki, a subset of MinWikiSplit. ABCD achieves comparable performance as two parsing baselines on MinWiki. On DeSSE, which has a more even balance of complex sentence types, our model achieves higher accuracy on the number of atomic sentences than an encoder-decoder baseline. Results include a detailed error analysis.


title: A bijection between the indecomposable summands of two multiplicity free tilting modules

abstract: We show that there is a reflection type bijection between the indecomposable summands of two multiplicity free tilting modules $X$ and $Y$. This bijection fixes the common indecomposable summands of $X$ and $Y$ and sends indecomposable projective (resp., injective) summands of exactly one module to non-projective (resp., non-injective) summands of the other. Moreover, this bijection interchanges the two possible non-isomorphic complements of an almost complete tilting module.


title: MIMHD: Accurate and Efficient Hyperdimensional Inference Using Multi-Bit In-Memory Computing

abstract: Hyperdimensional Computing (HDC) is an emerging computational framework that mimics important brain functions by operating over high-dimensional vectors, called hypervectors (HVs). In-memory computing implementations of HDC are desirable since they can significantly reduce data transfer overheads. All existing in-memory HDC platforms consider binary HVs where each dimension is represented with a single bit. However, utilizing multi-bit HVs allows HDC to achieve acceptable accuracies in lower dimensions which in turn leads to higher energy efficiencies. Thus, we propose a highly accurate and efficient multi-bit in-memory HDC inference platform called MIMHD. MIMHD supports multi-bit operations using ferroelectric field-effect transistor (FeFET) crossbar arrays for multiply-and-add and FeFET multi-bit content-addressable memories for associative search. We also introduce a novel hardware-aware retraining framework (HWART) that trains the HDC model to learn to work with MIMHD. For six popular datasets and 4000 dimension HVs, MIMHD using 3-bit (2-bit) precision HVs achieves (i) average accuracies of 92.6% (88.9%) which is 8.5% (4.8%) higher than binary implementations; (ii) 84.1x (78.6x) energy improvement over a GPU, and (iii) 38.4x (34.3x) speedup over a GPU, respectively. The 3-bit $\times$ is 4.3x and 13x faster and more energy-efficient than binary HDC accelerators while achieving similar accuracies.


title: A Simple and Practical Approach to Improve Misspellings in OCR Text

abstract: The focus of our paper is the identification and correction of non-word errors in OCR text. Such errors may be the result of incorrect insertion, deletion, or substitution of a character, or the transposition of two adjacent characters within a single word. Or, it can be the result of word boundary problems that lead to run-on errors and incorrect-split errors. The traditional N-gram correction methods can handle single-word errors effectively. However, they show limitations when dealing with split and merge errors. In this paper, we develop an unsupervised method that can handle both errors. The method we develop leads to a sizable improvement in the correction rates. This tutorial paper addresses very difficult word correction problems - namely incorrect run-on and split errors - and illustrates what needs to be considered when addressing such problems. We outline a possible approach and assess its success on a limited study.


title: Cancellation properties of graded and nonunital rings. Graded clean and graded exchange Leavitt path algebras

abstract: Various authors have been generalizing some unital ring properties to nonunital rings. We consider properties related to cancellation of modules (being unit-regular, having stable range one, being directly finite, exchange, or clean) and their "local" versions. We explore their relationships and extend the defined concepts to graded rings. With graded clean and graded exchange rings suitably defined, we study how these properties behave under the formation of graded matrix rings. We exhibit properties of a graph $E$ which are equivalent to the unital Leavitt path algebra $L_K(E)$ being graded clean. We also exhibit some graph properties which are necessary and some which are sufficient for $L_K(E)$ to be graded exchange.


title: Experimental Quantum Computing to Solve Network DC Power Flow Problem

abstract: Practical quantum computing applications to power grids are nonexistent at the moment. This paper investigates how a fundamental grid problem, namely DC power flow, can be solved using quantum computing. Power flow is the most widely used power system analysis technique, either as a stand-alone application or embedded in other applications; therefore, its fast and accurate solution is of utmost significance for grid operators. We base our studies on the Harrow-Hassidim-Lloyd (HHL) quantum algorithm, which has a proven theoretical speedup over classical algorithms in solving a system of linear equations. Practical studies on a quantum computer are conducted using the WSCC 9-bus system.


title: Strategic Liquidity Provision in Uniswap v3

abstract: Uniswap v3 is the largest decentralized exchange for digital currencies. A novelty of its design is that it allows a liquidity provider (LP) to allocate liquidity to one or more closed intervals of the price of an asset instead of the full range of possible prices. An LP earns fee rewards proportional to the amount of its liquidity allocation when prices move in this interval. This induces the problem of {\em strategic liquidity provision}: smaller intervals result in higher concentration of liquidity and correspondingly larger fees when the price remains in the interval, but with higher risk as prices may exit the interval leaving the LP with no fee rewards. Although reallocating liquidity to new intervals can mitigate this loss, it comes at a cost, as LPs must expend gas fees to do so. We formalize the dynamic liquidity provision problem and focus on a general class of strategies for which we provide a neural network-based optimization framework for maximizing LP earnings. We model a single LP that faces an exogenous sequence of price changes that arise from arbitrage and non-arbitrage trades in the decentralized exchange. We present experimental results informed by historical price data that demonstrate large improvements in LP earnings over existing allocation strategy baselines. Moreover we provide insight into qualitative differences in optimal LP behaviour in different economic environments.


title: Pure Exploration in Kernel and Neural Bandits

abstract: We study pure exploration in bandits, where the dimension of the feature representation can be much larger than the number of arms. To overcome the curse of dimensionality, we propose to adaptively embed the feature representation of each arm into a lower-dimensional space and carefully deal with the induced model misspecification. Our approach is conceptually very different from existing works that can either only handle low-dimensional linear bandits or passively deal with model misspecification. We showcase the application of our approach to two pure exploration settings that were previously under-studied: (1) the reward function belongs to a possibly infinite-dimensional Reproducing Kernel Hilbert Space, and (2) the reward function is nonlinear and can be approximated by neural networks. Our main results provide sample complexity guarantees that only depend on the effective dimension of the feature spaces in the kernel or neural representations. Extensive experiments conducted on both synthetic and real-world datasets demonstrate the efficacy of our methods.


title: GRB afterglow parameters in the era of TeV observations: the case of GRB~190114C

abstract: The afterglow of GRB~190114C has been observed at {60-1200~s} after the burst in the {sub-TeV range} by the {MAGIC Cherenkov} telescope. The simultaneous observations in X-ray range, which is presumed to be of synchrotron origin, and in sub-TeV range, where the emission is presumed to be inverse Compton, provide new stringent constraints on the conditions within the emitting regions {and their evolution in time}. While the additional data contain a lot of new information, it turns out that fitting both the X-ray and the TeV emission, is much more complicated than what was originally anticipated. We find that optical flux measurements provide important complementary information that in combination with TeV measurements breaks degeneracy in the parameter space. We present here a numerical fit to the multi-wavelength observed spectrum using a new code that calculates the single-zone synchrotron including self-Compton emission taking into account the exact Klein-Nishina cross-sections as well as pair production via absorption of the high-energy photons inside the emitting zone and the emission from the resulting secondary pairs. We also present a revised set of single zone parameters and a method for fitting the data to the observations. Our model for GRB~190114C that fits all the observations, from the optical data point to the sub-TeV range, suggests that it is in the fast cooling regime. The inferred parameters for observations {at two separate moments of time} show significant deviations from some of the common expectations in afterglow modeling but are all consistent with the predictions of the {pair-balance} model.


title: Effect of forward speed on the level-crossing distribution of kinematic variables in multidirectional ocean waves

abstract: The influence of forward speed on stochastic free-surface crossing, in a Gaussian wave field, is investigated. The case of a material point moving with a constant forward speed is considered; the wave field is assumed stationary in time, and homogeneous in space. The focus is on up-crossing events, which are defined as the material point crossing the free surface, into the water domain. The effect of the Doppler shift (induced by the forward speed) on the up-crossing frequency, and the related conditional joint distribution of wave kinematic variables is analytically investigated. Some general trends are illustrated through different examples, where three kinds of wave direction distribution are considered: unidirectional, short-crested anisotropic, and isotropic. The way the developed approach may be used in the context of slamming on marine structures is briefly discussed.


title: Listen to Your Favorite Melodies with img2Mxml, Producing MusicXML from Sheet Music Image by Measure-based Multimodal Deep Learning-driven Assembly

abstract: Deep learning has recently been applied to optical music recognition (OMR). However, currently OMR processing from various sheet music images still lacks precision to be widely applicable. Here, we present an MMdA (Measure-based Multimodal deep learning (DL)-driven Assembly) method allowing for end-to-end OMR processing from various images including inclined photo images. Using this method, measures are extracted by a deep learning model, aligned, and resized to be used for inference of given musical symbol components by using multiple deep learning models in sequence or in parallel. Use of each standardized measure enables efficient training of the models and accurate adjustment of five staff lines in each measure. Multiple musical symbol component category models with a small number of feature types can represent a diverse set of notes and other musical symbols including chords. This MMdA method provides a solution to end-to-end OMR processing with precision.


title: Parallel Exponential Time Differencing Methods for Geophysical Flow Simulations

abstract: Two ocean models are considered for geophysical flow simulations: the multi-layer shallow water equations and the multi-layer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock-Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. These methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrate a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

categories: physics.ao-ph cs.NA math.NA physics.flu-dyn


title: Clustering of check-in sequences using the mixture Markov chain process

abstract: This work is devoted to the clustering of check-in sequences from a geosocial network. We used the mixture Markov chain process as a mathematical model for time-dependent types of data. For clustering, we adjusted the Expectation-Maximization (EM) algorithm. As a result, we obtained highly detailed communities (clusters) of users of the now defunct geosocial network, Weeplaces.


title: Hairy Black holes in General Minimal Massive Gravity

abstract: In this work, we investigate the near horizon and asymptotic symmetries of static and rotating hairy$-$AdS black hole in the framework of general minimal massive gravity. We obtain energy, angular momentum and entropy of the solutions. Then we show that our results for these quantities are consistent with the first law of black hole thermodynamics. By considering the near horizon geometry of black hole, we find near horizon conserved charges and their algebra. By writing the algebra of conserved charges in terms of Fourier modes we have obtained the conserved charges in terms of zero modes.


title: Analysis of the Evolution of Parametric Drivers of High-End Sea-Level Hazards

abstract: Climate models are critical tools for developing strategies to manage the risks posed by sea-level rise to coastal communities. While these models are necessary for understanding climate risks, there is a level of uncertainty inherent in each parameter in the models. This model parametric uncertainty leads to uncertainty in future climate risks. Consequently, there is a need to understand how those parameter uncertainties impact our assessment of future climate risks and the efficacy of strategies to manage them. Here, we use random forests to examine the parametric drivers of future climate risk and how the relative importances of those drivers change over time. We find that the equilibrium climate sensitivity and a factor that scales the effect of aerosols on radiative forcing are consistently the most important climate model parametric uncertainties throughout the 2020 to 2150 interval for both low and high radiative forcing scenarios. The near-term hazards of high-end sea-level rise are driven primarily by thermal expansion, while the longer-term hazards are associated with mass loss from the Antarctic and Greenland ice sheets. Our results highlight the practical importance of considering time-evolving parametric uncertainties when developing strategies to manage future climate risks.


title: HydroPower Plant Planning for Resilience Improvement of Power Systems using Fuzzy-Neural based Genetic Algorithm

abstract: This paper will propose a novel technique for optimize hydropower plant in small scale based on load frequency control (LFC) which use self-tuning fuzzy Proportional- Derivative (PD) method for estimation and prediction of planning. Due to frequency is not controlled by any dump load or something else, so this power plant is under dynamic frequency variations that will use PD controller which optimize by fuzzy rules and then with neural deep learning techniques and Genetic Algorithm optimization. The main purpose of this work is because to maintain frequency in small-hydropower plant at nominal value. So, proposed controller means Fuzzy PD optimization with Genetic Algorithm will be used for LFC in small scale of hydropower system. The proposed schema can be used in different designation of both diesel generator and mini-hydropower system at low stream flow. It is also possible to use diesel generator at the hydropower system which can be turn off when Consumer demand is higher than electricity generation. The simulation will be done in MATLAB/Simulink to represent and evaluate the performance of this control schema under dynamic frequency variations. Spiking Neural Network (SNN) used as the main deep learning techniques to optimizing this load frequency control which turns into Deep Spiking Neural Network (DSNN). Obtained results represented that the proposed schema has robust and high-performance frequency control in comparison to other methods.


title: Introducing Quantum Entanglement to First-Year Students: Resolving the Trilemma

abstract: While quantum mechanics (QM) is covered at length in introductory physics textbooks, the concept of quantum entanglement is typically not covered at all, despite its importance in the rapidly growing area of quantum information science and its extensive experimental confirmation. Thus, physics educators are left to their own devices as to how to introduce this important concept. Regardless of how a physics educator chooses to introduce quantum entanglement, they face a trilemma involving its mysterious Bell-inequality-violating correlations. They can compromise on the the completeness of their introduction and simply choose not to share that fact, totally ignoring the 2022 Nobel Prize in Physics. They can frustrate their more curious students by introducing the mystery and simply telling them that the QM formalism with its associated (equally mysterious) conservation law maps beautifully to the experiments, so there is nothing else that needs to be said. Or, they can compromise the rigor of their presentation and attempt to resolve the mystery by venturing into the metaphysical quagmire of competing QM interpretations. Herein, we resolve this trilemma in precisely the same way that Einstein resolved the mysteries of time dilation and length contraction that existed in the late nineteenth century. That is, we resort to "principle" explanation based on the mathematical consequences of "empirically discovered" facts. Indeed, our principle account of quantum entanglement is even based on the same principle Einstein used, i.e., the relativity principle or "no preferred reference frame." Thus, this principle resolution of the trilemma is as complete, satisfying, analytically rigorous, and accessible as the standard introduction of special relativity for first-year physics students.


title: Empathy and Hope: Resource Transfer to Model Inter-country Social Media Dynamics

abstract: The ongoing COVID-19 pandemic resulted in significant ramifications for international relations ranging from travel restrictions, global ceasefires, and international vaccine production and sharing agreements. Amidst a wave of infections in India that resulted in a systemic breakdown of healthcare infrastructure, a social welfare organization based in Pakistan offered to procure medical-grade oxygen to assist India -- a nation which was involved in four wars with Pakistan in the past few decades. In this paper, we focus on Pakistani Twitter users' response to the ongoing healthcare crisis in India. While #IndiaNeedsOxygen and #PakistanStandsWithIndia featured among the top-trending hashtags in Pakistan, divisive hashtags such as #EndiaSaySorryToKashmir simultaneously started trending. Against the backdrop of a contentious history including four wars, divisive content of this nature, especially when a country is facing an unprecedented healthcare crisis, fuels further deterioration of relations. In this paper, we define a new task of detecting \emph{supportive} content and demonstrate that existing \emph{NLP for social impact} tools can be effectively harnessed for such tasks within a quick turnaround time. We also release the first publicly available data set at the intersection of geopolitical relations and a raging pandemic in the context of India and Pakistan.


title: Deep learning for improved global precipitation in numerical weather prediction systems

abstract: The formation of precipitation in state-of-the-art weather and climate models is an important process. The understanding of its relationship with other variables can lead to endless benefits, particularly for the world's monsoon regions dependent on rainfall as a support for livelihood. Various factors play a crucial role in the formation of rainfall, and those physical processes are leading to significant biases in the operational weather forecasts. We use the UNET architecture of a deep convolutional neural network with residual learning as a proof of concept to learn global data-driven models of precipitation. The models are trained on reanalysis datasets projected on the cubed-sphere projection to minimize errors due to spherical distortion. The results are compared with the operational dynamical model used by the India Meteorological Department. The theoretical deep learning-based model shows doubling of the grid point, as well as area averaged skill measured in Pearson correlation coefficients relative to operational system. This study is a proof-of-concept showing that residual learning-based UNET can unravel physical relationships to target precipitation, and those physical constraints can be used in the dynamical operational models towards improved precipitation forecasts. Our results pave the way for the development of online, hybrid models in the future.


title: KrigR -- A tool for downloading and statistically downscaling climate reanalysis data

abstract: Advances in climate science have rendered obsolete gridded observation data sets commonly used in macroecological analyses. Novel climate reanalysis products outperform legacy data products in accuracy, temporal resolution, and provision of uncertainty metrics. Consequently, there is an urgent need to develop a workflow through which to integrate these improved data into analyses. The ERA5 product family are the latest and most advanced global reanalysis products created by the ECMWF. These data products offer up to 83 essential climate variables at hourly intervals for the time-period of 1981 to today with preliminary back-extensions being available for 1950-1981. Spatial resolutions range from 30km (ERA5) to 11km (ERA5-Land) and can be statistically downscaled to study-requirements at finer spatial resolutions. Kriging is one such method to interpolate data to finer resolutions and has the advantages that one can leverage additional covariate information and obtain the uncertainty associated with the downscaling. The KrigR R Package enables users to (1) download ERA5(-Land) climate reanalysis data for a user-specified region, and time-period, (2) aggregate these climate products to desired temporal resolutions and metrics, (3) acquire topographical co-variates, and (4) statistically downscale spatial data to a user-specified resolution using co-variate data via kriging. This R package provides a workflow for implementation of state-of-the-art climate data into biological analyses while avoiding issues of storage limitations at high temporal and spatial resolutions by providing data according to user-needs rather than in global data sets. Consequently, KrigR provides a toolbox to obtain a wide range of tailored climate data at unprecedented combinations of high temporal and spatial resolutions thus enabling the use of world-leading climate data in the R-environment.


title: Fluctuations of water quality time series in rivers follow superstatistics

abstract: Superstatistics is a general method from nonequilibrium statistical physics which has been applied to a variety of complex systems, ranging from hydrodynamic turbulence to traffic delays and air pollution dynamics. Here, we investigate water quality time series (such as dissolved oxygen concentrations and electrical conductivity) as measured in rivers, and provide evidence that they exhibit superstatistical behaviour. Our main example are time series as recorded in the river Chess in South East England. Specifically, we use seasonal detrending and empirical mode decomposition (EMD) to separate trends from fluctuations for the measured data. With either detrending method, we observe heavy-tailed fluctuation distributions, which are well described by a log-normal superstatistics for dissolved oxygen. Contrarily, we find a double peaked non-standard superstatistics for the electrical conductivity data, which we model using two combined $\chi^2$-distributions.

categories: physics.ao-ph physics.bio-ph physics.data-an stat.AP


title: The Dirichlet problem for a family of totally degenerate differential operators

abstract: In the framework of the Potential Theory we prove existence and uniqueness for the Perron-Wiener-Brelot solution to the Dirichlet problem related to a family of totally degenerate, in the sense of Bony, differential operators. We also state and prove a Wiener-type criterium and an exterior cone condition for boundary regularity. Our results apply to a wide family of strongly degenerate operators that includes the following example~$\l = t^2\Delta_x + \langle x, \nabla_y \rangle -\partial_t$, for~$(x,y,t) \in \R^N \times \R^{N} \times \R$.


title: Pricing American options with the Runge-Kutta-Legendre finite difference scheme

abstract: This paper presents the Runge-Kutta-Legendre finite difference scheme, allowing for an additional shift in its polynomial representation. A short presentation of the stability region, comparatively to the Runge-Kutta-Chebyshev scheme follows. We then explore the problem of pricing American options with the Runge-Kutta-Legendre scheme under the one factor Black-Scholes and the two factor Heston stochastic volatility models, as well as the pricing of butterfly spread and digital options under the uncertain volatility model, where a Hamilton-Jacobi-Bellman partial differential equation needs to be solved. We explore the order of convergence in these problems, as well as the option greeks stability, compared to the literature and popular schemes such as Crank-Nicolson, with Rannacher time-stepping.


title: Classical and quantum $f(R)$ cosmology: The big rip, the little rip and the little sibling of the big rip

abstract: The big rip, the little rip and the little sibling of the big rip are cosmological doomsdays predicted by some phantom dark energy models that could describe the future evolution of our Universe. When the universe evolves towards either of these future cosmic events all bounded structures and, ultimately, space-time itself are ripped apart. Nevertheless, it is commonly belief that quantum gravity effects may smooth or even avoid these classically predicted singularities. In this review, we discuss the classical and quantum occurrence of these riplike events in the scheme of metric $f(R)$ theories of gravity. The quantum analysis is performed in the framework of $f(R)$ quantum geometrodynamics. In this context, we analyse the fulfilment of the DeWitt criterion for the avoidance of these singular fates.


title: More stochastic expansions for the pricing of vanilla options with cash dividends

abstract: There is no exact closed form formula for pricing of European options with discrete cash dividends under the model where the underlying asset price follows a piecewise lognormal process with jumps at dividend ex-dates. This paper presents alternative expansions based on the technique of Etore and Gobet, leading to more robust first, second and third-order expansions across the range of strikes and the range of dividend dates.


title: Volume Rendering of Neural Implicit Surfaces

abstract: Neural volume rendering became increasingly popular recently due to its success in synthesizing novel views of a scene from a sparse set of input images. So far, the geometry learned by neural volume rendering techniques was modeled using a generic density function. Furthermore, the geometry itself was extracted using an arbitrary level set of the density function leading to a noisy, often low fidelity reconstruction. The goal of this paper is to improve geometry representation and reconstruction in neural volume rendering. We achieve that by modeling the volume density as a function of the geometry. This is in contrast to previous work modeling the geometry as a function of the volume density. In more detail, we define the volume density function as Laplace's cumulative distribution function (CDF) applied to a signed distance function (SDF) representation. This simple density representation has three benefits: (i) it provides a useful inductive bias to the geometry learned in the neural volume rendering process; (ii) it facilitates a bound on the opacity approximation error, leading to an accurate sampling of the viewing ray. Accurate sampling is important to provide a precise coupling of geometry and radiance; and (iii) it allows efficient unsupervised disentanglement of shape and appearance in volume rendering. Applying this new density representation to challenging scene multiview datasets produced high quality geometry reconstructions, outperforming relevant baselines. Furthermore, switching shape and appearance between scenes is possible due to the disentanglement of the two.


title: Geometric Horizons: A Frame Approach

abstract: In the numerical investigation of the physical merger of two black holes, it is crucial to locate a black hole locally. This is usually done utilizing an apparent horizon. An alternative proposal is to identify a geometric horizon (GH), which is characterized by a surface in the spacetime on which the curvature tensor or its covariant derivatives are algebraically special. This necessitates the choice of a special null frame, which we shall refer to as an algebraically preferred null frame (APNF). The GH is then identified by surfaces of vanishing scalar curvature invariants but, unfortunately, these are difficult to compute. However, the algebraic nature of a GH means that the APNF plays a central role and suggests a null frame approach to characterizing the GH. Indeed, if we employ the Cartan-Karlhede algorithm to completely fix the null frame invariantly, then all of the remaining non vanishing components of the curvature tensor and its covariant derivatives are Cartan scalars. Hence the GH is characterized by the vanishing of certain Cartan scalars. A null frame approach is useful in the numerical investigation of the merger of two black holes in general, but we will focus on the application to identifying a GH. We begin with a review of the use of APNF and GH in previous work. The APNF is then defined and chosen so that the Weyl tensor is algebraically special, and we must examine the covariant derivatives of the Weyl tensor in this frame. We show how to invariantly fix the null frame, and hence characterize the APNF, and describe how to then identify the GH using the zero-set of certain Cartan scalars. Our ultimate aim is to apply this frame formalism to the numerical collapse of two black holes. As an example, we investigate the axisymmetric evolution of a two black hole Kastor-Traschen spacetime.


title: Automatic Head Overcoat Thickness Measure with NASNet-Large-Decoder Net

abstract: Transmission electron microscopy (TEM) is one of the primary tools to show microstructural characterization of materials as well as film thickness. However, manual determination of film thickness from TEM images is time-consuming as well as subjective, especially when the films in question are very thin and the need for measurement precision is very high. Such is the case for head overcoat (HOC) thickness measurements in the magnetic hard disk drive industry. It is therefore necessary to develop software to automatically measure HOC thickness. In this paper, for the first time, we propose a HOC layer segmentation method using NASNet-Large as an encoder and then followed by a decoder architecture, which is one of the most commonly used architectures in deep learning for image segmentation. To further improve segmentation results, we are the first to propose a post-processing layer to remove irrelevant portions in the segmentation result. To measure the thickness of the segmented HOC layer, we propose a regressive convolutional neural network (RCNN) model as well as orthogonal thickness calculation methods. Experimental results demonstrate a higher dice score for our model which has lower mean squared error and outperforms current state-of-the-art manual measurement.


title: Dominance-based linear formulation for the Anchor-Robust Project Scheduling Problem

abstract: In project scheduling under processing times uncertainty, the Anchor-Robust Project Scheduling Problem is to find a baseline schedule of bounded makespan and a max-weight subset of jobs whose starting times are guaranteed. The problem was proven NP-hard even for budgeted uncertainty. In the present work we design mixed-integer programming (MIP) formulations that are valid for a variety of uncertainty sets encompassing budgeted uncertainty. A new dominance among solutions is proposed, resulting into an MIP formulation. We further study the combinatorial structure of the problem. Non-trivial polynomial cases under budgeted uncertainty are exhibited, where the dominance-based formulation yields a polyhedral characterization of integer solutions. In more general cases, the dominance-based formulation is shown to be tighter than all previously known formulations. In numerical experiments we investigate how the formulation performs on instances around the polynomial cases, for both budgeted uncertainty sets and more elaborate uncertainty sets involving several budgets.


title: On Positivity Bias in Negative Reviews

abstract: Prior work has revealed that positive words occur more frequently than negative words in human expressions, which is typically attributed to positivity bias, a tendency for people to report positive views of reality. But what about the language used in negative reviews? Consistent with prior work, we show that English negative reviews tend to contain more positive words than negative words, using a variety of datasets. We reconcile this observation with prior findings on the pragmatics of negation, and show that negations are commonly associated with positive words in negative reviews. Furthermore, in negative reviews, the majority of sentences with positive words express negative opinions based on sentiment classifiers, indicating some form of negation.


title: Multivariate Generating Functions for Information Spread on Multi-Type Random Graphs

abstract: We study the spread of information on multi-type directed random graphs. In such graphs the vertices are partitioned into distinct types (communities) that have different transmission rates between themselves and with other types. We construct multivariate generating functions and use multi-type branching processes to derive an equation for the size of the large out-components in multi-type random graphs with a general class of degree distributions. We use our methods to analyse the spread of epidemics and verify the results with population based simulations


title: Dubrovnik Skein Theory and Power Sum Elements

abstract: In this work, we extend some results from the Kauffman bracket and HOMFLYPT skein theories to the Kauffman (Dubrovnik) skein theory. A definition is given for ``power sum" type elements $\widetilde{P}_k$ in the Dubrovnik skein algebra of the annulus $\mathcal{D}(A)$. These elements generalize the Chebyshev polynomials often used when studying Kauffman bracket skein algebras. Threadings of the $\widetilde{P}_k$ are used as generators in a presentation of the Dubrovnik skein algebra of the torus $\mathcal{D}(T^2)$, where they are shown to satisfy simple relations. This description of $\mathcal{D}(T^2)$ is used to describe the natural action of this algebra on the skein module of the solid torus. We give evidence that the universal character rings for the orthogonal and symplectic Lie groups correspond to the skein algebra $\mathcal{D}(A)$ such that the Schur functions of type either B, C or D correspond to annular closures $\widetilde{Q}_\lambda$ of minimal idempotents of the Birman-Murakami-Wenzl algebras $BMW_n$. We also record some miscellaneous applications of the $\widetilde{P}_k$, such as commutation relations for the annular closures of BMW symmetrizers $\widetilde{Q}_{(n)}$ and an expression of central elements of $BMW_n$ in terms of Jucys-Murphy elements.


title: Stochastic Batch Acquisition: A Simple Baseline for Deep Active Learning

abstract: We examine a simple stochastic strategy for adapting well-known single-point acquisition functions to allow batch active learning. Unlike acquiring the top-K points from the pool set, score- or rank-based sampling takes into account that acquisition scores change as new data are acquired. This simple strategy for adapting standard single-sample acquisition strategies can even perform just as well as compute-intensive state-of-the-art batch acquisition functions, like BatchBALD or BADGE, while using orders of magnitude less compute. In addition to providing a practical option for machine learning practitioners, the surprising success of the proposed method in a wide range of experimental settings raises a difficult question for the field: when are these expensive batch acquisition methods pulling their weight?


title: Faster Randomized Methods for Orthogonality Constrained Problems

abstract: Recent literature has advocated the use of randomized methods for accelerating the solution of various matrix problems arising throughout data science and computational science. One popular strategy for leveraging randomization is to use it as a way to reduce problem size. However, methods based on this strategy lack sufficient accuracy for some applications. Randomized preconditioning is another approach for leveraging randomization, which provides higher accuracy. The main challenge in using randomized preconditioning is the need for an underlying iterative method, thus randomized preconditioning so far have been applied almost exclusively to solving regression problems and linear systems. In this article, we show how to expand the application of randomized preconditioning to another important set of problems prevalent across data science: optimization problems with (generalized) orthogonality constraints. We demonstrate our approach, which is based on the framework of Riemannian optimization and Riemannian preconditioning, on the problem of computing the dominant canonical correlations and on the Fisher linear discriminant analysis problem. For both problems, we evaluate the effect of preconditioning on the computational costs and asymptotic convergence, and demonstrate empirically the utility of our approach.


title: Intrinsic versus extrinsic orbital and electronic reconstructions at complex oxide interfaces

abstract: The interface between the insulators LaAlO$_3$ and SrTiO$_3$ accommodates a two-dimensional electron liquid (2DEL) -- a high mobility electron system exhibiting superconductivity as well as indications of magnetism and correlations. While this flagship oxide heterostructure shows promise for electronics applications, the origin and microscopic properties of the 2DEL remain unclear. The uncertainty remains in part because the electronic structures of such nanoscale buried interfaces are difficult to probe, and is compounded by the variable presence of oxygen vacancies and coexistence of both localized and delocalized charges. These various complications have precluded decisive tests of intrinsic electronic and orbital reconstruction at this interface. Here we overcome prior difficulties by developing an interface analysis based on the inherently interface-sensitive resonant x-ray reflectometry. We discover a high charge density of 0.5 electrons per interfacial unit cell for samples above the critical LaAlO$_3$ thickness, and extract the depth dependence of both the orbital and electronic reconstructions near the buried interface. We find that the majority of the reconstruction phenomena are confined to within 2 unit cells of the interface, and we quantify how oxygen vacancies significantly affect the electronic system. Our results provide strong support for the existence of polarity induced electronic reconstruction, clearly separating its effects from those of oxygen vacancies.


title: A Practical & Unified Notation for Information-Theoretic Quantities in ML

abstract: A practical notation can convey valuable intuitions and concisely express new ideas. Information theory is of importance to machine learning, but the notation for information-theoretic quantities is sometimes opaque. We propose a practical and unified notation and extend it to include information-theoretic quantities between observed outcomes (events) and random variables. This includes the point-wise mutual information known in NLP and mixed quantities such as specific surprise and specific information in the cognitive sciences and information gain in Bayesian optimal experimental design. We apply our notation to prove a version of Stirling's approximation for binomial coefficients mentioned by MacKa (2003) using new intuitions. We also concisely rederive the evidence lower bound for variational auto-encoders and variational inference in approximate Bayesian neural networks. Furthermore, we apply the notation to a popular information-theoretic acquisition function in Bayesian active learning which selects the most informative (unlabelled) samples to be labelled by an expert and extend this acquisition function to the core-set problem with the goal of selecting the most informative samples given the labels.


title: Families of similar simplices inscribed in most smoothly embedded spheres

abstract: Let $\Delta$ denote a non-degenerate $k$-simplex in $\mathbb{R}^k$. The set $\text{Sim}(\Delta)$ of simplices in $\mathbb{R}^k$ similar to $\Delta$ is diffeomorphic to $O(k)\times [0,\infty)\times \mathbb{R}^k$, where the factor in $O(k)$ is a matrix called the {\em pose}. Among $(k-1)$-spheres smoothly embedded in $\mathbb{R}^k$ and isotopic to the identity, there is a dense family of spheres, for which the subset of $\text{Sim}(\Delta)$ of simplices inscribed in each embedded sphere contains a similar simplex of every pose $U\in O(k)$. Further, the intersection of $\text{Sim}(\Delta)$ with the configuration space of $k+1$ distinct points on an embedded sphere is a manifold whose top homology class maps to the top class in $O(k)$ via the pose map. This gives a high dimensional generalization of classical results on inscribing families of triangles in plane curves. We use techniques established in our previous paper on the square-peg problem where we viewed inscribed simplices in spheres as transverse intersections of submanifolds of compactified configuration spaces.


title: Assessing epidemic curves for evidence of superspreading

abstract: The expected number of secondary infections arising from each index case, referred to as the reproduction or $R$ number, is a vital summary statistic for understanding and managing epidemic diseases. There are many methods for estimating $R$; however, few explicitly model heterogeneous disease reproduction, which gives rise to superspreading within the population. We propose a parsimonious discrete-time branching process model for epidemic curves that incorporates heterogeneous individual reproduction numbers. Our Bayesian approach to inference illustrates that this heterogeneity results in less certainty on estimates of the time-varying cohort reproduction number $R_t$. We apply these methods to a COVID-19 epidemic curve for the Republic of Ireland and find support for heterogeneous disease reproduction. Our analysis allows us to estimate the expected proportion of secondary infections attributable to the most infectious proportion of the population. For example, we estimate that the 20% most infectious index cases account for approximately 75-98% of the expected secondary infections with 95% posterior probability. In addition, we highlight that heterogeneity is a vital consideration when estimating $R_t$.


title: Polycyclic Codes Associated with Trinomials: Good Codes and Open Questions

abstract: Polycyclic codes are a generalization of cyclic and constacyclic codes. Even though they have been known since 1972 and received some attention more recently, there have not been many studies on polycyclic codes. This paper presents an in-depth investigation of polycyclic codes associated with trinomials. Our results include a number of facts about trinomials, some properties of polycyclic codes, and many new quantum codes derived from polycyclic codes. We also state several conjectures about polynomials and polycyclic codes. Hence, we show useful features of polycyclic codes and present some open problems related to them.


title: It's All in the Heads: Using Attention Heads as a Baseline for Cross-Lingual Transfer in Commonsense Reasoning

abstract: Commonsense reasoning is one of the key problems in natural language processing, but the relative scarcity of labeled data holds back the progress for languages other than English. Pretrained cross-lingual models are a source of powerful language-agnostic representations, yet their inherent reasoning capabilities are still actively studied. In this work, we design a simple approach to commonsense reasoning which trains a linear classifier with weights of multi-head attention as features. To evaluate this approach, we create a multilingual Winograd Schema corpus by processing several datasets from prior work within a standardized pipeline and measure cross-lingual generalization ability in terms of out-of-sample performance. The method performs competitively with recent supervised and unsupervised approaches for commonsense reasoning, even when applied to other languages in a zero-shot manner. Also, we demonstrate that most of the performance is given by the same small subset of attention heads for all studied languages, which provides evidence of universal reasoning capabilities in multilingual encoders.


title: Nature of the orbital angular momentum (OAM) fields in a multilayered fiber

abstract: We provide a theoretical analysis of the nature of the orbital angular momentum (OAM) modal fields in a multilayered fiber, such as the step-index fiber and the ring-core fiber. In a detailed study of the vector field solutions of the step-index fiber (in the exponential basis), we discover that the polarization-induced field component is a modified scalar OAM field (as opposed to a standard OAM scalar field) with a shifted intensity pattern in the weakly guiding approximation (WGA); the familiar intensity donut pattern is reduced or increased in radius depending upon whether it is a case of spin-alignment or anti-alignment with the OAM. Such a shift in the intensity pattern appears to be a general feature of the field of a multilayered fiber as seen from an extension to the ring-core fiber. Additionally, we derive a general expression for the polarization-correction to the scalar propagation constant, which includes, for the first time, the contribution of the polarization-induced field. All the analytic expressions are illustrated and validated numerically with application to a step-index fiber, whose analytic solutions are well-known.


title: The Rate of Convergence of Variation-Constrained Deep Neural Networks

abstract: Multi-layer feedforward networks have been used to approximate a wide range of nonlinear functions. An important and fundamental problem is to understand the learnability of a network model through its statistical risk, or the expected prediction error on future data. To the best of our knowledge, the rate of convergence of neural networks shown by existing works is bounded by at most the order of $n^{-1/4}$ for a sample size of $n$. In this paper, we show that a class of variation-constrained neural networks, with arbitrary width, can achieve near-parametric rate $n^{-1/2+\delta}$ for an arbitrarily small positive constant $\delta$. It is equivalent to $n^{-1 +2\delta}$ under the mean squared error. This rate is also observed by numerical experiments. The result indicates that the neural function space needed for approximating smooth functions may not be as large as what is often perceived. Our result also provides insight to the phenomena that deep neural networks do not easily suffer from overfitting when the number of neurons and learning parameters rapidly grow with $n$ or even surpass $n$. We also discuss the rate of convergence regarding other network parameters, including the input dimension, network layer, and coefficient norm.


title: Heat transport in an ordered harmonic chain in presence of a uniform magnetic field

abstract: We consider heat transport across a harmonic chain of charged particles, with transverse degrees of freedom, in the presence of a uniform magnetic field. For an open chain connected to heat baths at the two ends we obtain the nonequilibrium Green's function expression for the heat current. This expression involves two different Green's functions which can be identified as corresponding respectively to scattering processes within or between the two transverse waves. The presence of the magnetic field leads to two phonon bands of the isolated system and we show that the net transmission can be written as a sum of two distinct terms attributable to the two bands. Exact expressions are obtained for the current in the thermodynamic limit, for the the cases of free and fixed boundary conditions. In this limit, we find that at small frequency $\omega$, the effective transmission has the frequency-dependence $\omega^{3/2}$ and $\omega^{1/2}$ for fixed and free boundary conditions respectively. This is in contrast to the zero magnetic field case where the transmission has the dependence $\omega^2$ and $\omega^0$ for the two boundary conditions respectively, and can be understood as arising from the quadratic low frequency phonon dispersion.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: Towards Consistent Predictive Confidence through Fitted Ensembles

abstract: Deep neural networks are behind many of the recent successes in machine learning applications. However, these models can produce overconfident decisions while encountering out-of-distribution (OOD) examples or making a wrong prediction. This inconsistent predictive confidence limits the integration of independently-trained learning models into a larger system. This paper introduces separable concept learning framework to realistically measure the performance of classifiers in presence of OOD examples. In this setup, several instances of a classifier are trained on different parts of a partition of the set of classes. Later, the performance of the combination of these models is evaluated on a separate test set. Unlike current OOD detection techniques, this framework does not require auxiliary OOD datasets and does not separate classification from detection performance. Furthermore, we present a new strong baseline for more consistent predictive confidence in deep models, called fitted ensembles, where overconfident predictions are rectified by transformed versions of the original classification task. Fitted ensembles can naturally detect OOD examples without requiring auxiliary data by observing contradicting predictions among its components. Experiments on MNIST, SVHN, CIFAR-10/100, and ImageNet show fitted ensemble significantly outperform conventional ensembles on OOD examples and are possible to scale.


title: New Techniques for Worldline Integration

abstract: The worldline formalism provides an alternative to Feynman diagrams in the construction of amplitudes and effective actions that shares some of the superior properties of the organization of amplitudes in string theory. In particular, it allows one to write down integral representations combining the contributions of large classes of Feynman diagrams of different topologies. However, calculating these integrals analytically without splitting them into sectors corresponding to individual diagrams poses a formidable mathematical challenge. We summarize the history and state of the art of this problem, including some natural connections to the theory of Bernoulli numbers and polynomials and multiple zeta values.


title: Bayesian estimation for collisional thermometry

abstract: Quantum thermometry exploits the high level of control in coherent devices to offer enhanced precision for temperature estimation. This highlights the need for constructing concrete estimation strategies. Of particular importance is collisional thermometry, where a series of ancillas are sent sequentially to probe the system's temperature. In this paper we put forth a complete framework for analyzing collisional thermometry using Bayesian inference. The approach is easily implementable and experimentally friendly. Moreover, it is guaranteed to always saturate the Cramer-Rao bound in the long-time limit. Subtleties concerning the prior information about the system's temperature are also discussed, and analyzed in terms of a modified Cramer-Rao bound associated to Van Trees and Sch\"utzenberger.


title: Noncommutative Differential K-theory

abstract: We introduce a differential extension of algebraic K-theory of an algebra using Karoubi's Chern character. In doing so, we develop a necessary theory of secondary transgression forms as well as a differential refinement of the smooth Serre--Swan correspondence. Our construction subsumes the differential K-theory of a smooth manifold when the algebra is complex-valued smooth functions. Furthermore, our construction fits into a noncommutative differential cohomology hexagon diagram.


title: Reachability Analysis of Convolutional Neural Networks

abstract: Deep convolutional neural networks have been widely employed as an effective technique to handle complex and practical problems. However, one of the fundamental problems is the lack of formal methods to analyze their behavior. To address this challenge, we propose an approach to compute the exact reachable sets of a network given an input domain, where the reachable set is represented by the face lattice structure. Besides the computation of reachable sets, our approach is also capable of backtracking to the input domain given an output reachable set. Therefore, a full analysis of a network's behavior can be realized. In addition, an approach for fast analysis is also introduced, which conducts fast computation of reachable sets by considering selected sensitive neurons in each layer. The exact pixel-level reachability analysis method is evaluated on a CNN for the CIFAR10 dataset and compared to related works. The fast analysis method is evaluated over a CNN CIFAR10 dataset and VGG16 architecture for the ImageNet dataset.