title: Inference on quantile processes with a finite number of clusters

abstract: I introduce a generic method for inference on entire quantile and regression quantile processes in the presence of a finite number of large and arbitrarily heterogeneous clusters. The method asymptotically controls size by generating statistics that exhibit enough distributional symmetry such that randomization tests can be applied. The randomization test does not require ex-ante matching of clusters, is free of user-chosen parameters, and performs well at conventional significance levels with as few as five clusters. The method tests standard (non-sharp) hypotheses and can even be asymptotically similar in empirically relevant situations. The main focus of the paper is inference on quantile treatment effects but the method applies more broadly. Numerical and empirical examples are provided.


title: Robust construction of differential emission measure profiles using a regularized maximum likelihood method

abstract: Extreme-ultraviolet (EUV) observations provide considerable insight into evolving physical conditions in the active solar atmosphere. For a prescribed density and temperature structure, it is straightforward to construct the corresponding differential emission measure profile $\xi(T)$, such that $\xi(T) \, dT$ is proportional to the emissivity from plasma in the temperature range $[T, T + dT]$. Here we study the inverse problem of obtaining a valid $\xi(T)$ profile from a set of EUV spectral line intensities observed at a pixel within a solar image. Our goal is to introduce and develop a regularized maximum likelihood (RML) algorithm designed to address the mathematically ill-posed problem of constructing differential emission measure profiles from a discrete set of EUV intensities in specified wavelength bands, specifically those observed by the Atmospheric Imaging Assembly (AIA) on the NASA Solar Dynamics Observatory. The RML method combines features of maximum likelihood and regularized approaches used by other authors. It is also guaranteed to produce a positive definite differential emission measure profile. We evaluate and compare the effectiveness of the method against other published algorithms, using both simulated data generated from parametric differential emission profile forms, and AIA data from a solar eruptive event on 2010 November 3. Similarities and differences between the differential emission measure profiles and maps reconstructed by the various algorithms are discussed. The RML inversion method is mathematically rigorous, computationally efficient, and produces acceptable measures of performance in the following three key areas: fidelity to the data, accuracy in the reconstruction, and robustness in the presence of data noise. As such, it shows considerable promise for computing differential emission measure profiles from datasets of discrete spectral lines.


title: Weakly asymmetric facilitated exclusion process

abstract: We consider the facilitated exclusion process, an interacting particle system on the integer line where particles hop to one of their left or right neighbouring site only when the other neighbouring site is occupied by a particle. A peculiarity of this system is that, starting from the step initial condition, the density profile develops a downward jump discontinuity around the position of the first particle, unlike other exclusion processes such as the asymmetric simple exclusion process (ASEP). In the weakly asymmetric regime, we show that the field of particle positions around the jump discontinuity converges to the solution of the multiplicative noise stochastic heat equation (i.e. the exponential of a solution to the KPZ equation) on a half-line subject to Dirichlet boundary condition, with initial condition given by the derivative of a Dirac delta function. We prove this result by reformulating the problem in terms of ASEP on a half-line with a boundary reservoir, for which we extend known proofs of convergence to deal with Dirichlet boundary condition and the very singular type of initial condition that arises in our case.


title: A Functorial Perspective on (Multi)computational Irreducibility

abstract: This article aims to provide a novel formalization of the concept of computational irreducibility in terms of the exactness of functorial correspondence between a category of data structures and elementary computations and a corresponding category of (1-dimensional) cobordisms. We proceed to demonstrate that, by equipping both categories with a symmetric monoidal structure and considering the case of higher-dimensional cobordism categories, we obtain a natural extension of this formalism that serves also to encompass non-deterministic or ``multiway'' computations, in which one quantifies not only the irreducibility in the behavior of a single (deterministic) computation path, but in the branching and merging behavior of an entire ``multiway system'' of such paths too. We finally outline how, in the most general case, the resulting symmetric monoidal functor may be considered to be adjoint to the functor characterizing the Atiyah-Segal axiomatization of a functorial quantum field theory. Thus, we conclude by arguing that the irreducibility of (multi)computations may be thought of as being dual to the locality of time evolution in functorial approaches to quantum mechanics and quantum field theory. In the process, we propose an extension of the methods of standard (monoidal) category theory, in which morphisms are effectively equipped with intrinsic computational complexity data, together with an algebra for how those complexities compose (both in sequence and in parallel, subject to the monoidal structure).


title: QoS Based Contract Design for Profit Maximization in IoT-Enabled Data Markets

abstract: The massive deployment of Internet of Things (IoT) devices, including sensors and actuators, is ushering in smart and connected communities of the future. The massive deployment of Internet of Things (IoT) devices, including sensors and actuators, is ushering in smart and connected communities of the future. The availability of real-time and high-quality sensor data is crucial for various IoT applications, particularly in healthcare, energy, transportation, etc. However, data collection may have to be outsourced to external service providers (SPs) due to cost considerations or lack of specialized equipment. Hence, the data market plays a critical role in such scenarios where SPs have different quality levels of available data, and IoT users have different application-specific data needs. The pairing between data available to the SP and users in the data market requires an effective mechanism design that considers the SPs' profitability and the quality-of-service (QoS) needs of the users. We develop a generic framework to analyze and enable such interactions efficiently, leveraging tools from contract theory and mechanism design theory. It can enable and empower emerging data sharing paradigms such as Sensing-as-a-Service (SaaS). The contract design creates a pricing structure for on-demand sensing data for IoT users. By considering a continuum of user types, we capture a diverse range of application requirements and propose optimal pricing and allocation rules that ensure QoS provisioning and maximum profitability for the SP. Furthermore, we provide analytical solutions for fixed distributions of user types to analyze the developed approach. For comparison, we consider the benchmark case assuming complete information of the user types and obtain optimal contract solutions. Finally, a case study is presented to demonstrate the efficacy of the proposed contract design framework.


title: Rotation Periods, Inclinations, and Obliquities of Cool Stars Hosting Directly Imaged Substellar Companions: Spin-Orbit Misalignments are Common

abstract: The orientation between a star's spin axis and a planet's orbital plane provides valuable information about the system's formation and dynamical history. For non-transiting planets at wide separations, true stellar obliquities are challenging to measure, but lower limits on spin-orbit orientations can be determined from the difference between the inclination of the star's rotational axis and the companion's orbital plane ($\Delta i$). We present results of a uniform analysis of rotation periods, stellar inclinations, and obliquities of cool stars (SpT $\gtrsim$ F5) hosting directly imaged planets and brown dwarf companions. As part of this effort, we have acquired new $v \sin i_*$ values for 22 host stars with the high-resolution Tull spectrograph at the Harlan J. Smith telescope. Altogether our sample contains 62 host stars with rotation periods, most of which are newly measured using light curves from the Transiting Exoplanet Survey Satellite. Among these, 53 stars have inclinations determined from projected rotational and equatorial velocities, and 21 stars predominantly hosting brown dwarfs have constraints on $\Delta i$. Eleven of these (52$^{+10}_{-11}$% of the sample) are likely misaligned, while the remaining ten host stars are consistent with spin-orbit alignment. As an ensemble, the minimum obliquity distribution between 10-250 AU is more consistent with a mixture of isotropic and aligned systems than either extreme scenario alone--pointing to direct cloud collapse, formation within disks bearing primordial alignments and misalignments, or architectures processed by dynamical evolution. This contrasts with stars hosting directly imaged planets, which show a preference for low obliquities. These results reinforce an emerging distinction between the orbits of long-period brown dwarfs and giant planets in terms of their stellar obliquities and orbital eccentricities.


title: Why the observed spin evolution of older-than-solar like stars might not require a dynamo mode change

abstract: The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, stellar activity, and consequent influence on companion planets. Observations of older than solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the sun. Here we take a different approach which considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it.

categories: astro-ph.SR astro-ph.EP astro-ph.HE astro-ph.IM


title: Science Priorities for the Extraction of the Solid MSR Samples from their Sample Tubes

abstract: Preservation of the chemical and structural integrity of samples that will be brought back from Mars is paramount to achieving the scientific objectives of MSR. Given our knowledge of the nature of the samples retrieved at Jezero by Perseverance, at least two options need to be tested for opening the sample tubes: (1) One or two radial cuts at the end of the tube to slide the sample out. (2) Two radial cuts at the ends of the tube and two longitudinal cuts to lift the upper half of the tube and access the sample. Strategy 1 will likely minimize contamination but incurs the risk of affecting the physical integrity of weakly consolidated samples. Strategy 2 will be optimal for preserving the physical integrity of the samples but increases the risk of contamination and mishandling of the sample as more manipulations and additional equipment will be needed. A flexible approach to opening the sample tubes is therefore required, and several options need to be available, depending on the nature of the rock samples returned. Both opening strategies 1 and 2 may need to be available when the samples are returned to handle different sample types (e.g., loosely bound sediments vs. indurated magmatic rocks). This question should be revisited after engineering tests are performed on analogue samples. The MSR sample tubes will have to be opened under stringent BSL4 conditions and this aspect needs to be integrated into the planning.


title: Learning Continuous Mesh Representation with Spherical Implicit Surface

abstract: As the most common representation for 3D shapes, mesh is often stored discretely with arrays of vertices and faces. However, 3D shapes in the real world are presented continuously. In this paper, we propose to learn a continuous representation for meshes with fixed topology, a common and practical setting in many faces-, hand-, and body-related applications. First, we split the template into multiple closed manifold genus-0 meshes so that each genus-0 mesh can be parameterized onto the unit sphere. Then we learn spherical implicit surface (SIS), which takes a spherical coordinate and a global feature or a set of local features around the coordinate as inputs, predicting the vertex corresponding to the coordinate as an output. Since the spherical coordinates are continuous, SIS can depict a mesh in an arbitrary resolution. SIS representation builds a bridge between discrete and continuous representation in 3D shapes. Specifically, we train SIS networks in a self-supervised manner for two tasks: a reconstruction task and a super-resolution task. Experiments show that our SIS representation is comparable with state-of-the-art methods that are specifically designed for meshes with a fixed resolution and significantly outperforms methods that work in arbitrary resolutions.


title: On Modeling Network Slicing Communication Resources with SARSA Optimization

abstract: Network slicing is a crucial enabler to support the composition and deployment of virtual network infrastructures required by the dynamic behavior of networks like 5G/6G mobile networks, IoT-aware networks, e-health systems, and industry verticals like the internet of vehicles (IoV) and industry 4.0. The communication slices and their allocated communication resources are essential in slicing architectures for resource orchestration and allocation, virtual network function (VNF) deployment, and slice operation functionalities. The communication slices provide the communications capabilities required to support slice operation, SLA guarantees, and QoS/ QoE application requirements. Therefore, this contribution proposes a networking slicing conceptual model to formulate the optimization problem related to the sharing of communication resources among communication slices. First, we present a conceptual model of network slicing, we then formulate analytically some aspects of the model and the optimization problem to address. Next, we proposed to use a SARSA agent to solve the problem and implement a proof of concept prototype. Finally, we present the obtained results and discuss them.


title: Regular Black Holes and Horizonless Ultra-Compact Objects in Lorentz-Violating Gravity

abstract: There is growing evidence that Ho\v{r}ava gravity may be a viable quantum theory of gravity. It is thus legitimate to expect that gravitational collapse in the full, non-projectable version of the theory should result in geometries that are free of spacetime singularities. Previous analyses have shown that such geometries must belong to one of the following classes: simply connected regular black holes with inner horizons; non-connected black holes "hiding" a wormhole mouth (black bounces); simply connected or non-connected horizonless compact objects. Here, we consider a singular black hole in the low-energy limit of non-projectable Ho\v{r}ava gravity, i.e. khronometric theory, and describe examples of its possible "regularisations", covering all of the viable classes. To our knowledge, these examples constitute the first instances of black holes with inner universal horizons, of black bounces and of stars with a de Sitter core in the context of Lorentz-violating theories of gravity.


title: Accelerating large-eddy simulations of clouds with Tensor Processing Units

abstract: Clouds, especially low clouds, are crucial for regulating Earth's energy balance and mediating the response of the climate system to changes in greenhouse gas concentrations. Despite their importance for climate, they remain relatively poorly understood and are inaccurately represented in climate models. A principal reason is that the high computational expense of simulating them with large-eddy simulations (LES) has inhibited broad and systematic numerical experimentation and the generation of large datasets for training parametrization schemes for climate models. Here we demonstrate LES of low clouds on Tensor Processing Units (TPUs), application-specific integrated circuits that were originally developed for machine learning applications. We show that TPUs in conjunction with tailored software implementations can be used to simulate computationally challenging stratocumulus clouds in conditions observed during the Dynamics and Chemistry of Marine Stratocumulus (DYCOMS) field study. The TPU-based LES code successfully reproduces clouds during DYCOMS and opens up the large computational resources available on TPUs to cloud simulations. The code enables unprecedented weak and strong scaling of LES, making it possible, for example, to simulate stratocumulus with $10\times$ speedup over real-time evolution in domains with a $34.7~\mathrm{km} \times 53.8~\mathrm{km}$ horizontal cross section. The results open up new avenues for computational experiments and for substantially enlarging the sample of LES available to train parameterizations of low clouds.

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


title: CMZoom III: Spectral Line Data Release

abstract: We present an overview and data release of the spectral line component of the SMA Large Program, \textit{CMZoom}. \textit{CMZoom} observed $^{12}$CO(2-1), $^{13}$CO(2-1) and C$^{18}$O(2-1), three transitions of H$_{2}$CO, several transitions of CH$_{3}$OH, two transitions of OCS and single transitions of SiO and SO, within gas above a column density of N(H$_2$)$\ge 10^{23}$\,cm$^{-2}$ in the Central Molecular Zone (CMZ; inner few hundred pc of the Galaxy). We extract spectra from all compact 1.3\,mm \emph{CMZoom} continuum sources and fit line profiles to the spectra. We use the fit results from the H$_{2}$CO 3(0,3)-2(0,2) transition to determine the source kinematic properties. We find $\sim 90$\% of the total mass of \emph{CMZoom} sources have reliable kinematics. Only four compact continuum sources are formally self-gravitating. The remainder are consistent with being in hydrostatic equilibrium assuming that they are confined by the high external pressure in the CMZ. Based on the mass and density of virially bound sources, and assuming star formation occurs within one free-fall time with a star formation efficiency of $10\% - 75\%$, we place a lower limit on the future embedded star-formation rate of $0.008 - 0.06$\,M$_{\odot}$\,yr$^{-1}$. We find only two convincing proto-stellar outflows, ruling out a previously undetected population of very massive, actively accreting YSOs with strong outflows. Finally, despite having sufficient sensitivity and resolution to detect high-velocity compact clouds (HVCCs), which have been claimed as evidence for intermediate mass black holes interacting with molecular gas clouds, we find no such objects across the large survey area.


title: Continuous-Time Formulations for Multi-Mode Project Scheduling

abstract: This paper reviews compact continuous-time formulations for the multi-mode resource-constrained project scheduling problem. Specifically, we first point out a serious flaw in an existing start-end-event-based formulation owing to inconsistent mode choices. We propose two options to formulate the missing constraints and we consider an equivalent reformulation with sparser constraint matrix. Second, we formulate an aggregate variant of an existing model that relies on on-off-events and clarify the role of mode consistency issues in such models. Third, we suggest two variants of an existing network flow formulation. We enhance our models by adapting several techniques that have been used previously, e.g., in cases with only a single mode. A large set of benchmark instances from the literature provides the basis for an up-to-date and fair computational study with an out-of-the-box solver package. We compare our models against two models from the literature. Our experiments assert confidently that network flow formulations prevail in the test bed, and they provide a hint on why event-based models become less competitive in multi-mode settings.


title: Indices of diagonalizable and universal realizability of spectra

abstract: A list $\Lambda =\{\lambda _{1},\ldots ,\lambda _{n}\}$ of complex numbers (repeats allowed) is said to be \textit{realizable} if it is the spectrum of an entrywise nonnegative matrix $A$. $\Lambda $ is \textit{diagonalizably realizable} if the realizing matrix $A$ is diagonalizable. $\Lambda $ is said to be \textit{universally realizable} if it is \textit{\ realizable} for each possible Jordan canonical form allowed by $\Lambda .$ Here, we study the connection between diagonalizable realizability and universal realizability of spectra. In particular, we establish \textit{\ indices of realizability} for diagonalizable and universal realizability. We also define the merge of two spectra and we prove a result that allow us to easily decide, in many cases, about the universal realizability of spectra.


title: Physics Simulation Via Quantum Graph Neural Network

abstract: We develop and implement two realizations of quantum graph neural networks (QGNN), applied to the task of particle interaction simulation. The first QGNN is a speculative quantum-classical hybrid learning model that relies on the ability to directly utilize superposition states as classical information to propagate information between particles. The second is an implementable quantum-classical hybrid learning model that propagates particle information directly through the parameters of $RX$ rotation gates. A classical graph neural network (CGNN) is also trained in the same task. Both the Speculative QGNN and CGNN act as controls against the Implementable QGNN. Comparison between classical and quantum models is based on the loss value and accuracy of each model. Overall, each model had a high learning efficiency, in which the loss value rapidly approached zero during training; however, each model was moderately inaccurate. Comparing performances, our results show that the Implementable QGNN has a potential advantage over the CGNN. Additionally, we show that a slight alteration in hyperparameters in the CGNN notably improves accuracy, suggesting that further fine tuning could mitigate the issue of moderate inaccuracy in each model.


title: Anomalous flux state in higher-order topological superconductors

abstract: We investigate the anomalous flux state of interacting higher-order topological superconductors (HOTSC) protected by rotation symmetries. By introducing a $\pi$ superconducting flux in 2D HOTSC, we demonstrate the existence of a robust zero mode trapped at the flux center. Remarkably, the rotation symmetry and fermion parity display projective representation inside the $\pi$ flux with $N=2$ supersymmetry algebra. A similar gapless flux pattern also exists in 3D HOTSCs, the flux lines of which carry anomalous helical modes that cannot be realized on purely one-dimensional lattice models. Notably, these exotic phenomena can be manifested in a 2D frustrated quantum magnet whose low energy excitation characterizes emergent Majoranas with HOTSC band structure and the $Z_2$ flux exhibiting supersymmetries.


title: SensePOLAR: Word sense aware interpretability for pre-trained contextual word embeddings

abstract: Adding interpretability to word embeddings represents an area of active research in text representation. Recent work has explored thepotential of embedding words via so-called polar dimensions (e.g. good vs. bad, correct vs. wrong). Examples of such recent approaches include SemAxis, POLAR, FrameAxis, and BiImp. Although these approaches provide interpretable dimensions for words, they have not been designed to deal with polysemy, i.e. they can not easily distinguish between different senses of words. To address this limitation, we present SensePOLAR, an extension of the original POLAR framework that enables word-sense aware interpretability for pre-trained contextual word embeddings. The resulting interpretable word embeddings achieve a level of performance that is comparable to original contextual word embeddings across a variety of natural language processing tasks including the GLUE and SQuAD benchmarks. Our work removes a fundamental limitation of existing approaches by offering users sense aware interpretations for contextual word embeddings.


title: Inverse Quantum Fourier Transform Inspired Algorithm for Unsupervised Image Segmentation

abstract: Image segmentation is a very popular and important task in computer vision. In this paper, inverse quantum Fourier transform (IQFT) for image segmentation has been explored and a novel IQFT-inspired algorithm is proposed and implemented by leveraging the underlying mathematical structure of the IQFT. Specifically, the proposed method takes advantage of the phase information of the pixels in the image by encoding the pixels' intensity into qubit relative phases and applying IQFT to classify the pixels into different segments automatically and efficiently. To the best of our knowledge, this is the first attempt of using IQFT for unsupervised image segmentation. The proposed method has low computational cost comparing to the deep learning-based methods and more importantly it does not require training, thus make it suitable for real-time applications. The performance of the proposed method is compared with K-means and Otsu-thresholding. The proposed method outperforms both of them on the PASCAL VOC 2012 segmentation benchmark and the xVIEW2 challenge dataset by as much as 50% in terms of mean Intersection-Over-Union (mIOU).


title: Aspects of $\mathbb{Z}_N$ rank-2 gauge theory in $(2+1)$ dimensions: construction schemes, holonomies, and sublattice one-form symmetries

abstract: Rank-2 toric code (R2TC), a prototypical archetype of the discrete rank-2 symmetric gauge theory, has properties that differ from those of the standard toric code. Specifically, it features a blending of UV and IR in its ground state, restricted mobility of its quasiparticles, and variations in the braiding statistics of its quasiparticles based on their position. In this paper, we investigate various aspects of $\mathbb{Z}_N$ rank-2 gauge theory in ${(2+1)}$-dimensional spacetime. Firstly, we demonstrate that $U(1)$ rank-2 gauge theory can arise from ${U(1)\times U(1)}$ rank-1 gauge theory after condensing the gauge charges in a specific way. This construction scheme of $U(1)$ rank-2 gauge theory carries over to the $\mathbb{Z}_N$ case simply by Higgsing $U(1)$ to $\mathbb{Z}_N$, after which the resulting rank-2 gauge theory can be tuned to the R2TC. The holonomy operators of R2TC are readily identified using this scheme and are given clear physical interpretation as the pair creation/annihilation of various monopoles and dipoles. Explicit tensor network construction of the ground states of R2TC are given as two copies of the ground states of Kitaev's toric code that are `sewn together' according to the condensation scheme. In addition, through a similar anyon condensation protocol, we present a double semion version of rank-2 toric code whose flux excitations exhibit restricted mobility and semionic statistics. Finally, we identify the generalized discrete symmetries of the R2TC, which are much more complex than typical 1-form symmetries. They include conventional and unconventional 1-form symmetries, such as framed 1-form symmetries and what we call sublattice 1-form symmetries. Using these, we interpret the R2TC's unique properties (UV/IR mixing, position-dependent braiding, etc.) from the modern perspective of generalized spontaneous symmetry breaking and 't Hooft anomalies.


title: Optimal coverage-based placement of static leak detection devices for pipeline water supply networks

abstract: In this paper, we provide a mathematical optimization-based framework to determine the location of leak detection devices along a network. Assuming that the devices are endowed with a known coverage area, we analyze two different models. The first model aims to minimize the number of devices to be located in order to (fully or partially) cover the volume of the network. In the second model, the number of devices is given, and the goal is to locate them to provide a coverage volume as broad as possible. Unlike other approaches in the literature, in our models, it is not assumed that the devices are located on the network (nodes or edges) but in the whole space and that the different segments in the networks may be partially covered, which allows for more flexible coverage. We also derive a method to construct initial solutions as well as a math-heuristic approach for solving the problem for larger instances. We report the results of a series of experiments on real-world water supply pipeline networks, supporting the validity of our models.


title: Two-Photon EXchange -- TPEX

abstract: We propose a new measurement of the ratio of positron-proton to electron-proton, elastic scattering at DESY to determine the contributions beyond single-photon exchange, which are essential to the QED description of the most fundamental process in hadronic physics. A 20~cm long liquid hydrogen target together with the extracted beam from the DESY synchrotron would yield an average luminosity of $2.12\times10^{35}$~cm$^{-2}\cdot$s$^{-1}\cdot$sr$^{-1}$ ($\sim200$ times the luminosity achieved by OLYMPUS). A commissioning run at 2 GeV followed by measurements at 3 GeV would provide new data up to $Q^2=4.6$~(GeV/$c$)$^2$ (twice the range of current measurements). Lead tungstate calorimeters would be used to detect the scattered leptons at polar angles of $30^\circ$, $50^\circ$, $70^\circ$, $90^\circ$, and $110^\circ$. The measurements could be scheduled to not interfere with the operation of PETRA. We present rate estimates and simulations for the planned measurements including background considerations. Initial measurements at the DESY test beam facility using prototype lead tungstate calorimeters in 2019, 2021, and 2022 were made to check the Monte Carlo simulations and the performance of the calorimeters. These tests also investigated different readout schemes (triggered and streaming). Various upgrades are possible to shorten the running time and to make higher beam energies and thus greater $Q^2$ ranges accessible.


title: Causal Abstraction for Faithful Model Interpretation

abstract: A faithful and interpretable explanation of an AI model's behavior and internal structure is a high-level explanation that is human-intelligible but also consistent with the known, but often opaque low-level causal details of the model. We argue that the theory of causal abstraction provides the mathematical foundations for the desired kinds of model explanations. In causal abstraction analysis, we use interventions on model-internal states to rigorously assess whether an interpretable high-level causal model is a faithful description of an AI model. Our contributions in this area are: (1) We generalize causal abstraction to cyclic causal structures and typed high-level variables. (2) We show how multi-source interchange interventions can be used to conduct causal abstraction analyses. (3) We define a notion of approximate causal abstraction that allows us to assess the degree to which a high-level causal model is a causal abstraction of a lower-level one. (4) We prove constructive causal abstraction can be decomposed into three operations we refer to as marginalization, variable-merge, and value-merge. (5) We formalize the XAI methods of LIME, causal effect estimation, causal mediation analysis, iterated nullspace projection, and circuit-based explanations as special cases of causal abstraction analysis.


title: Density Mediated Spin Correlations Drive Edge to Bulk Flow Transition in Active Chiral Matter

abstract: We demonstrate that edge currents develop in active chiral matter -- composed of spinning disk-shaped grains with chirally arranged tilted legs confined in a circular vibrating chamber -- due to boundary shielding over a wide range of densities corresponding to a gas, fluid, and crystal. The edge currents are then shown to increasingly drive circulating bulk flows with area fraction $\phi$ due to increasing spin-coupling between neighbors mediated by frictional contacts, as percolating clusters develop. Edge currents are observed even in the dilute limit. While, at low $\phi$, the average flux vanishes except within a distance of a single particle diameter of the boundary, the penetration depth grows with increasing $\phi$ till a solid body rotation is achieved corresponding to the highest packing, where the particles are fully caged with hexagonal order and spin in phase with the entire packing. A coarse-grained model, based on the increased collisional interlocking of the particles with $\phi$ and the emergence of order, captures the observed flow fields.


title: Refined theoretical values of field and mass isotope shifts in thallium to extract charge radii of Tl isotopes

abstract: Electronic factors for the field and mass isotope shifts in the $6p\ ^{2}P_{3/2} \to 7s\ ^{2}S_{1/2}$ (535 nm), $6p\ ^{2}P_{1/2} \to 6d\ ^{2}D_{3/2}$ (277 nm) and $6p\ ^{2}P_{1/2} \to 7s\ ^{2}S_{1/2}$ (378~nm) transitions in the neutral thallium were calculated within the high-order relativistic coupled cluster approach. These factors were used to reinterpret previous experimental isotope shift measurements in terms of charge radii of a wide range of Tl isotopes. Good agreement between theoretical and experimental King-plot parameters was found for the $6p\ ^{2}P_{3/2} \to 7s\ ^{2}S_{1/2}$ and $6p\ ^{2}P_{1/2} \to 6d\ ^{2}D_{3/2}$ transitions. It was shown that the value of the specific mass shift factor for the $6p\ ^{2}P_{3/2} \to 7s\ ^{2}S_{1/2}$ transition is not negligible compared to the value of normal mass shift in contrary to what had been suggested previously. Theoretical uncertainties in mean square charge radii were estimated. They were substantially reduced compared to the previously ascribed ones and amounted to less than 2.6%. The achieved accuracy paves the way to a more reliable comparison of the charge radii trends in the lead region.


title: An analytic formula for entraining CAPE in mid-latitude storm environments

abstract: This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from the storm environment. Extending previous formulas derived in Peters et al. (2020), entrainment is connected to the background environment via an analytic manipulation of the equations of motion that yields a direct correspondence between the storm relative flow and the updraft radius, and an inverse scaling between the updraft radius squared and entrainment rate. These concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, a second formula is derived that accounts for updraft enhancement via kinetic energy extracted from the cloud's background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity $w_{max}$ in an updraft. Hence, these formulas are evaluated using $w_{max}$ from past numerical modeling studies. Both of the new formulas improve predictions of $w_{max}$ substantially over undiluted CAPE, ECAPE with a prescribed entrainment rate, and the ECAPE formula from Peters et al. (2020). The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of $\sqrt{2\text{CAPE}}$ by $w_{max}$ in simulations, and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies.


title: How do "technical" design-choices made when building algorithmic decision-making tools for criminal justice authorities create constitutional dangers?

abstract: This two part paper argues that seemingly "technical" choices made by developers of machine-learning based algorithmic tools used to inform decisions by criminal justice authorities can create serious constitutional dangers, enhancing the likelihood of abuse of decision-making power and the scope and magnitude of injustice. Drawing on three algorithmic tools in use, or recently used, to assess the "risk" posed by individuals to inform how they should be treated by criminal justice authorities, we integrate insights from data science and public law scholarship to show how public law principles and more specific legal duties that are rooted in these principles, are routinely overlooked in algorithmic tool-building and implementation. We argue that technical developers must collaborate closely with public law experts to ensure that if algorithmic decision-support tools are to inform criminal justice decisions, those tools are configured and implemented in a manner that is demonstrably compliant with public law principles and doctrine, including respect for human rights, throughout the tool-building process.


title: $RO(C_2)$-graded Cohomology of $C_2$-equivariant Eilenberg-Mac Lane spaces

abstract: In this paper, we calculate $RO(C_2)$-graded cohomology of $C_2$-equivariant Eilenberg-Mac Lane spaces $K(\underline{Z/2}, n+\sigma)$ for $n\geq 0$. These can be used to give the relation between equivariant lambda algebra and equivariant Adams resolution and equivariant unstable Adams spectral sequence, which are defined in author`s dissertation.


title: How do "technical" design-choices made when building algorithmic decision-making tools for criminal justice authorities create constitutional dangers? Part II

abstract: This two-part paper argues that seemingly "technical" choices made by developers of machine-learning based algorithmic tools used to inform decisions by criminal justice authorities can create serious constitutional dangers, enhancing the likelihood of abuse of decision-making power and the scope and magnitude of injustice. Drawing on three algorithmic tools in use, or recently used, to assess the "risk" posed by individuals to inform how they should be treated by criminal justice authorities, we integrate insights from data science and public law scholarship to show how public law principles and more specific legal duties that are rooted in these principles, are routinely overlooked in algorithmic tool-building and implementation. We argue that technical developers must collaborate closely with public law experts to ensure that if algorithmic decision-support tools are to inform criminal justice decisions, those tools are configured and implemented in a manner that is demonstrably compliant with public law principles and doctrine, including respect for human rights, throughout the tool-building process.


title: Measurement of the $B^{0} \rightarrow D^{*-} \ell^{+} \nu_{\ell}$ branching ratio and $|V_{cb}|$ with a fully reconstructed accompanying $B$ meson in 2019-2021 Belle II data

abstract: We present a measurement of the $B^{0} \rightarrow D^{*-} \ell^{+} \nu_{\ell}$ ($\ell=e,\mu$) branching ratio and of the CKM parameter $|V_{cb}|$ using signal decays accompanied by a fully reconstructed $B$ meson. The Belle II data set of electron-positron collisions at the $\Upsilon(4S)$ resonance, corresponding to 189.3$\,$fb$^{-1}$ of integrated luminosity, is analyzed. With the Caprini-Lellouch-Neubert form factor parameterization, the parameters $\eta_{\rm EW} F(1) |V_{cb}|$ and $\rho^{2}$ are extracted, where $\eta_{\rm EW}$ is an electroweak correction, $F(1)$ is a normalization factor and $\rho^{2}$ is a form factor shape parameter. We reconstruct 516 signal decays and thereby obtain $\mathcal{B} (B^{0} \rightarrow D^{*-} \ell^{+} \nu_{\ell} ) = \left(5.27 \pm 0.22~\rm{\left(stat\right)} \pm 0.38~\rm{\left(syst\right)}\right) \%$, $\eta_{EW} F(1) |V_{cb}| \times 10^{3} = 34.6 \pm 1.8~\rm{\left(stat\right)} \pm 1.7~\rm{\left(syst\right)}$, and $\rho^{2} = 0.94 \pm 0.18~\rm{\left(stat\right)} \pm 0.11~\rm{\left(syst\right)}$.


title: Fitting Bell Curves to Data Distributions using Visualization

abstract: Idealized probability distributions, such as normal or other curves, lie at the root of confirmatory statistical tests. But how well do people understand these idealized curves? In practical terms, does the human visual system allow us to match sample data distributions with hypothesized population distributions from which those samples might have been drawn? And how do different visualization techniques impact this capability? This paper shares the results of a crowdsourced experiment that tested the ability of respondents to fit normal curves to four different data distribution visualizations: bar histograms, dotplot histograms, strip plots, and boxplots. We find that the crowd can estimate the center (mean) of a distribution with some success and little bias. We also find that people generally overestimate the standard deviation, which we dub the "umbrella effect" because people tend to want to cover the whole distribution using the curve, as if sheltering it from the heavens above, and that strip plots yield the best accuracy.


title: About the connection of the electron binding energy of a single carbon anion with binding energies of an electron attached to carbon molecules

abstract: We demonstrate that the model of zero-range potentials can be successfully employed for the description of attached electrons in atomic and molecular anions, for example, negatively charged carbon clusters. To illustrate the capability of the model we calculate the energies of the attached electron for the family of carbon cluster anions consisting of two-, three- (equilateral triangle), and four (tetrahedron) carbon atoms equidistant from each other as well as for a C3 molecule having a chain structure. The considered approach can be easily extended to carbon clusters containing an arbitrary number of atoms arranged in an arbitrary configuration.


title: How To Optimize My Blockchain? A Multi-Level Recommendation Approach

abstract: Aside from the conception of new blockchain architectures, existing blockchain optimizations in the literature primarily focus on system or data-oriented optimizations within prevailing blockchains. However, since blockchains handle multiple aspects ranging from organizational governance to smart contract design, a holistic approach that encompasses all the different layers of a given blockchain system is required to ensure that all optimization opportunities are taken into consideration. In this vein, we define a multi-level optimization recommendation approach that identifies optimization opportunities within a blockchain at the system, data, and user level. Multiple metrics and attributes are derived from a blockchain log and nine optimization recommendations are formalized. We implement an automated optimization recommendation tool, BlockOptR, based on these concepts. The system is extensively evaluated with a wide range of workloads covering multiple real-world scenarios. After implementing the recommended optimizations, we observe an average of 20% improvement in the success rate of transactions and an average of 40% improvement in latency.


title: An Architecture For Cooperative Mobile Health Applications

abstract: Mobile health applications are steadily gaining momentum in the modern world given the omnipresence of various mobile or WiFi connections. Given that the bandwidth of these connections increases over time, especially in conjunction with advanced modulation and error-correction codes, whereas the latency drops, the cooperation between mobile applications becomes gradually easier. This translates to reduced computational burden and heat dissipation for each isolated device at the expense of increased privacy risks. This chapter presents a configurable and scalable edge computing architecture for cooperative digital health mobile applications.


title: The Core Normal Type Ia Supernova 2019np: An Overall Spherical Explosion with an Aspherical Surface Layer and an Aspherical 56Ni Core

abstract: Optical spectropolarimetry of the normal thermonuclear supernova SN2019np from -14.5 to +14.5 days relative to B-band maximum detected an intrinsic continuum polarization, p(cont), of 0.21+-0.09% at the first epoch. Between days -11.5 to +05, p(cont) remained about 0 and by day +14.5 was again significant at 0.19+-0.10%. Not considering the first epoch, the dominant axis of SiII(6355A) was roughly constant, staying close to the continuum until both rotated in opposite directions on day +14.5. Detailed radiation-hydrodynamical simulations produce a very steep density slope in the outermost ejecta so that the low first-epoch p(const) of about 0.2% nevertheless suggests a separate structure with an axis ratio of about 2 in the outer carbon-rich 3...5E-3 Mo. Large-amplitude fluctuations in the polarization profiles and a flocculent appearance of the polar diagram for the CaII near-infrared triplet (NIR3) may be related by a common origin. The temporal evolution of the polarization spectra agrees with an off-center delayed detonation. The late-time increase in polarization and the possible change in position angle are also consistent with an aspherical 56Ni core. The p(cont) and the absorptions due to Si II(6355A) and the CaII NIR3 form in the same region of the extended photosphere, with an interplay between line occultation and thermalization producing p. Small-scale polarization features may be due to small-scale structures, but many could be related to atomic patterns of the quasi-continuum; the hardly have an equivalent in the total-flux spectrum. We compare SN2019np to other SNe and develop future objectives and strategies for SNIa spectropolarimetry.


title: Rate of Convergence in Multiple SLE using Random Matrix Theory

abstract: We provide an order of convergence for a version of the Carath\'eodory convergence for the multiple SLE model with a Dyson Brownian motion driver towards its hydrodynamic limit, for $\beta=1$ and $\beta=2$. The result is obtained by combining techniques from the field of Schramm-Loewner Evolutions with modern techniques from random matrices. Our approach shows how one can apply modern tools used in the proof of universality in random matrix theory, in the field of Schramm-Loewner Evolutions.


title: Millimetre-waves to Terahertz SISO and MIMO Continuous Variable Quantum Key Distribution

abstract: With the exponentially increased demands for large bandwidth, it is important to think about the best network platform as well as the security and privacy of the information in communication networks. Millimetre (mm)-waves and terahertz (THz) with high carrier frequency are proposed as the enabling technologies to overcome Shannons channel capacity limit of existing communication systems by providing ultrawide bandwidth signals. Mm-waves and THz are also able to build wireless links compatible with optical communication systems. However, most solid-state components that can operate reasonably efficiently at these frequency ranges (100GHz-10THz), especially sources and detectors, require cryogenic cooling, as is a requirement for most quantum systems. Here, we show that secure mm-waves and THz QKD can be achieved when the sources and detectors operate at cryogenic temperatures down to T= 4K. We compare single-input single-output (SISO) and multiple-input multiple-output (MIMO) Continuous Variable THz Quantum Key Distribution (CVQKD) schemes and find the positive secret key rate in the frequency ranges between f=100 GHz and 1 THz. Moreover, we find that the maximum transmission distance could be extended, the secret key rate could be improved in lower temperatures, and achieve a maximum secrete communication distance of more than 5 km at f=100GHz and T=4K by using 1024*1024 antennas. Our results may contribute to the efforts to develop next-generation secure wireless communication systems and quantum internet for applications from inter-satellite and deep space, to indoor and short-distance communications.


title: Decoding Structure-Spectrum Relationships with Physically Organized Latent Spaces

abstract: A new semi-supervised machine learning method for the discovery of structure-spectrum relationships is developed and demonstrated using the specific example of interpreting X-ray absorption near-edge structure (XANES) spectra. This method constructs a one-to-one mapping between individual structure descriptors and spectral trends. Specifically, an adversarial autoencoder is augmented with a novel rank constraint (RankAAE). The RankAAE methodology produces a continuous and interpretable latent space, where each dimension can track an individual structure descriptor. As a part of this process, the model provides a robust and quantitative measure of the structure-spectrum relationship by decoupling intertwined spectral contributions from multiple structural characteristics. This makes it ideal for spectral interpretation and the discovery of new descriptors. The capability of this procedure is showcased by considering five local structure descriptors and a database of over fifty thousand simulated XANES spectra across eight first-row transition metal oxide families. The resulting structure-spectrum relationships not only reproduce known trends in the literature, but also reveal unintuitive ones that are visually indiscernible in large data sets. The results suggest that the RankAAE methodology has great potential to assist researchers to interpret complex scientific data, test physical hypotheses, and reveal new patterns that extend scientific insight.


title: Blockchain For Mobile Health Applications: Acceleration With GPU Computing

abstract: Blockchain is a linearly linked, distributed, and very robust data structure. Originally proposed as part of the Bitcoin distributed stack, it found a number of applications in a number of fields, most notably in smart contracts, social media, secure IoT, and cryptocurrency mining. It ensures data integrity by distributing strongly encrypted data in widely redundant segments. Each new insertion requires verification and approval by the majority of the users of the blockchain. Both encryption and verification are computationally intensive tasks which cannot be solved with ordinary off-the-shelf CPUs. This has resulted in a renewed scientific interest in secure distributed communication and coordination protocols. Mobile health applications are growing progressively popular and have the enormous advantage of timely diagnosis of certain conditions. However, privacy concerns have been raised as mobile health application by default have access to highly sensitive personal data. This chapter presents concisely how blockchain can be applied to mobile health applications in order to enhance privacy.


title: On the zeros of certain Sheffer sequences and their cognate sequences

abstract: Given a Sheffer sequence of polynomials, we introduce the notion of an associated sequence called the cognate sequence. We study the relationship between the zeros of this pair of associated sequences and show that in case of an Appell sequence, as well as a more general family of Sheffer sequences, the zeros of the members of each sequence (for large n) are either real, or lie on a line $\Re (z)=c$. In addition to finding the zero locus, we also find the limiting probability distribution function of such sequences.


title: A Quantum Algorithm for Shapley Value Estimation

abstract: The introduction of the European Union's (EU) set of comprehensive regulations relating to technology, the General Data Protection Regulation, grants EU citizens the right to explanations for automated decisions that have significant effects on their life. This poses a substantial challenge, as many of today's state-of-the-art algorithms are generally unexplainable black boxes. Simultaneously, we have seen an emergence of the fields of quantum computation and quantum AI. Due to the fickle nature of quantum information, the problem of explainability is amplified, as measuring a quantum system destroys the information. As a result, there is a need for post-hoc explanations for quantum AI algorithms. In the classical context, the cooperative game theory concept of the Shapley value has been adapted for post-hoc explanations. However, this approach does not translate to use in quantum computing trivially and can be exponentially difficult to implement if not handled with care. We propose a novel algorithm which reduces the problem of accurately estimating the Shapley values of a quantum algorithm into a far simpler problem of estimating the true average of a binomial distribution in polynomial time.


title: Patch Locale of a Spectral Locale in Univalent Type Theory

abstract: Stone locales together with continuous maps form a coreflective subcategory of spectral locales and perfect maps. A proof in the internal language of an elementary topos was previously given by the second-named author. This proof can be easily translated to univalent type theory using resizing axioms. In this work, we show how to achieve such a translation without resizing axioms, by working with large, locally small, and small complete frames with small bases. This turns out to be nontrivial and involves predicative reformulations of several fundamental concepts of locale theory.


title: PL-genus of surfaces in homology balls

abstract: We consider manifold-knot pairs $(Y,K)$ where $Y$ is a homology sphere that bounds a homology ball. We show that the minimum genus of a PL surface $\Sigma$ in a homology ball $X$ such that $\partial (X, \Sigma) = (Y, K)$ can be arbitrarily large. Equivalently, the minimum genus of a surface cobordism in a homology cobordism from $(Y, K)$ to any knot in $S^3$ can be arbitrarily large. The proof relies on Heegaard Floer homology.


title: Graphene amplifier reaches the quantum-noise limit

abstract: News & Views on the graphene-based Josephson parametric amplifiers.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.supr-con quant-ph


title: Extended Corona Models of X-ray Reverberation in the AGN 1H~0707-495 and IRAS 13224-3809

abstract: We fit a new vertically extended corona model to previously measured reverberation time lags observed by \emph{XMM-Newton} in two extremely variable Narrow Line Seyfert 1 Active Galactic Nuclei (AGN), 1H~0707-495 and IRAS~13224-3809, in a variety of similarly observed flux groups and explore the model in all observations over a 16 year period. The model employs two X-ray sources located along the black hole rotational axis at height, $h_1$ and $h_2$ respectively. These sources have their associated photon indices $\Gamma_1$ and $\Gamma_2$ which respond to fluctuations in the disc with a maximum response duration of $t_\text{max}$ and a propagation delay between the response of the two of $t_\text{shift}$. We find that for 1H 0707-495, $h_2$ is significantly correlated with $\Gamma_1$ and anti-correlated with ionisation $\xi$. Whilst the 1H 0707-495 corona extends upwards, the emission appears softer and the disc is less ionised. We find similarities in IRAS 13224-3809, but significant anti-correlation between $\Gamma_2$ and both $t_\text{max}$ and $t_\text{shift}$. This suggests that when the IRAS 13224-3809 corona becomes softer while extending vertically upwards, the overall corona response occurs faster. This may also suggest that the inner disc also becomes more active. In addition, $\Gamma_1$ and $\Gamma_2$ are extreme, relatively less variable, but more separate in IRAS 13224-3809 than in 1H 0707-495. This suggests that the IRAS 13224-3809 corona may be more patchy in the sense that it has two more clear distinct spectral zones of $\Gamma_1$ and $\Gamma_2$ (possibly relating to two distinct zones of coronal temperature) when compared to 1H 0707-495.


title: Semi-infinite construction for the double Yangian of type $A_1^{(1)}$

abstract: We consider certain infinite dimensional modules of level 1 for the double Yangian $\text{DY}(\mathfrak{gl}_2)$ which are based on the Iohara-Kohno realization. We show that they possess topological bases of Feigin-Stoyanovsky-type, i.e. the bases expressed in terms of semi-infinite monomials of certain integrable operators which stabilize and satisfy the difference two condition. Finally, we give some applications of these bases to the representation theory of the corresponding quantum affine vertex algebra.


title: AGMN: Association Graph-based Graph Matching Network for Coronary Artery Semantic Labeling on Invasive Coronary Angiograms

abstract: Semantic labeling of coronary arterial segments in invasive coronary angiography (ICA) is important for automated assessment and report generation of coronary artery stenosis in the computer-aided diagnosis of coronary artery disease (CAD). Inspired by the training procedure of interventional cardiologists for interpreting the structure of coronary arteries, we propose an association graph-based graph matching network (AGMN) for coronary arterial semantic labeling. We first extract the vascular tree from invasive coronary angiography (ICA) and convert it into multiple individual graphs. Then, an association graph is constructed from two individual graphs where each vertex represents the relationship between two arterial segments. Using the association graph, the AGMN extracts the vertex features by the embedding module, aggregates the features from adjacent vertices and edges by graph convolution network, and decodes the features to generate the semantic mappings between arteries. By learning the mapping of arterial branches between two individual graphs, the unlabeled arterial segments are classified by the labeled segments to achieve semantic labeling. A dataset containing 263 ICAs was employed to train and validate the proposed model, and a five-fold cross-validation scheme was performed. Our AGMN model achieved an average accuracy of 0.8264, an average precision of 0.8276, an average recall of 0.8264, and an average F1-score of 0.8262, which significantly outperformed existing coronary artery semantic labeling methods. In conclusion, we have developed and validated a new algorithm with high accuracy, interpretability, and robustness for coronary artery semantic labeling on ICAs.


title: Construction of Wannier functions from the spectral moments of correlated electron systems

abstract: When the first four spectral moments are considered, spectral features missing in standard Kohn-Sham (KS) density-functional theory (DFT), such as upper and lower Hubbard bands, as well as spectral satellite peaks, can be described, and the bandwidths can be corrected. Therefore, we have devised a \textit{moment-functional based spectral density functional theory} (MFbSDFT) recently. However, many computational tools in theoretical solid state physics, such as the construction of maximally localized Wannier functions (MLWFs), have been developed for KS-DFT and require modifications if they are supposed to be used in MFbSDFT. Here, we show how generalized Wannier functions may be constructed from the first four spectral moment matrices. We call these functions \textit{maximally localized spectral moment Wannier functions} (MLSMWFs). We demonstrate how MLSMWFs may be used to compute the anomalous Hall effect (AHE) in fcc Ni by Wannier interpolation. More generally, MLSMWFs may be computed from the first $2P$ moments ($P=1,2,3,\dots$). Using more than 4 moments opens the perspective of reproducing all spectral features accurately in MFbSDFT.


title: Revisiting Pure State Transformations with Zero Communication

abstract: It is known that general convertibility of bipartite entangled states is not possible to arbitrary error without some classical communication. While some trade-offs between communication cost and conversion error have been proven, these bounds can be very loose. In particular, there are many cases in which tolerable error might be achievable using zero-communication protocols. In this work we address these cases by deriving the optimal fidelity of pure state conversions under local unitaries as well as local operations and shared randomness (LOSR). We also uses these results to explore catalytic conversions between pure states using zero communication.


title: Twisted recurrence for dynamical systems with exponential decay of correlations

abstract: We study the set of points returning infinitely often to a sequence of targets dependent on the starting points. With an assumption of decay of correlations for L1 against bounded variations, we prove a generalized quantitative recurrence result under Lipschitz twists.


title: Regularity results for a model in magnetohydrodynamics with imposed pressure

abstract: The magnetohydrodynamics (MHD) problem is most often studied in a framework where Dirichlet type boundary conditions on the velocity field is imposed. In this Note, we study the (MHD) system with pressure boundary condition, together with zero tangential trace for the velocity and the magnetic field. In a three-dimensional bounded possibly multiply connected domain, we first prove the existence of weak solutions in the Hilbert case, and later, the regularity in $W^{1,p}(\Omega)$ for $p\geq 2$ and in $W^{2,p}(\Omega)$ for $p\geq 6/5$ using the regularity results for some Stokes and elliptic problems with this type of boundary conditions. Furthermore, under the condition of small data, we obtain the existence and uniqueness of solutions in $W^{1,p}(\Omega)$ for $3/2<p<2$ by using a fixed-point technique over a linearized (MHD) problem.


title: Upper bounds for the Hausdorff dimension of Weierstrass curves

abstract: We produce an upper bound for the Hausdorff dimension of the graph of a Weierstrass-type function. Whilst strictly weaker than existing results, it has the advantage of being directly computable from the theory of hyperbolic iterated function systems.


title: Making sense of noise: introducing students to stochastic processes in order to better understand biological behaviors

abstract: Biological systems are characterized by the ubiquitous roles of weak, that is, non-covalent molecular interactions, small, often very small, numbers of specific molecules per cell, and Brownian motion. These combine to produce stochastic behaviors at all levels from the molecular and cellular to the behavioral. That said, students are rarely introduced to the ubiquitous role of stochastic processes in biological systems, and how they produce unpredictable behaviors. Here I present the case that they need to be and provide some suggestions as to how it might be approached.


title: The Berkelmans-Pries Feature Importance Method: A Generic Measure of Informativeness of Features

abstract: Over the past few years, the use of machine learning models has emerged as a generic and powerful means for prediction purposes. At the same time, there is a growing demand for interpretability of prediction models. To determine which features of a dataset are important to predict a target variable $Y$, a Feature Importance (FI) method can be used. By quantifying how important each feature is for predicting $Y$, irrelevant features can be identified and removed, which could increase the speed and accuracy of a model, and moreover, important features can be discovered, which could lead to valuable insights. A major problem with evaluating FI methods, is that the ground truth FI is often unknown. As a consequence, existing FI methods do not give the exact correct FI values. This is one of the many reasons why it can be hard to properly interpret the results of an FI method. Motivated by this, we introduce a new global approach named the Berkelmans-Pries FI method, which is based on a combination of Shapley values and the Berkelmans-Pries dependency function. We prove that our method has many useful properties, and accurately predicts the correct FI values for several cases where the ground truth FI can be derived in an exact manner. We experimentally show for a large collection of FI methods (468) that existing methods do not have the same useful properties. This shows that the Berkelmans-Pries FI method is a highly valuable tool for analyzing datasets with complex interdependencies.

categories: cs.LG cs.GT math.ST stat.AP stat.ML stat.TH


title: Efficient Preference-Based Reinforcement Learning Using Learned Dynamics Models

abstract: Preference-based reinforcement learning (PbRL) can enable robots to learn to perform tasks based on an individual's preferences without requiring a hand-crafted reward function. However, existing approaches either assume access to a high-fidelity simulator or analytic model or take a model-free approach that requires extensive, possibly unsafe online environment interactions. In this paper, we study the benefits and challenges of using a learned dynamics model when performing PbRL. In particular, we provide evidence that a learned dynamics model offers the following benefits when performing PbRL: (1) preference elicitation and policy optimization require significantly fewer environment interactions than model-free PbRL, (2) diverse preference queries can be synthesized safely and efficiently as a byproduct of standard model-based RL, and (3) reward pre-training based on suboptimal demonstrations can be performed without any environmental interaction. Our paper provides empirical evidence that learned dynamics models enable robots to learn customized policies based on user preferences in ways that are safer and more sample efficient than prior preference learning approaches. Supplementary materials and code are available at https://sites.google.com/berkeley.edu/mop-rl.


title: HADA: A Graph-based Amalgamation Framework in Image-text Retrieval

abstract: Many models have been proposed for vision and language tasks, especially the image-text retrieval task. All state-of-the-art (SOTA) models in this challenge contained hundreds of millions of parameters. They also were pretrained on a large external dataset that has been proven to make a big improvement in overall performance. It is not easy to propose a new model with a novel architecture and intensively train it on a massive dataset with many GPUs to surpass many SOTA models, which are already available to use on the Internet. In this paper, we proposed a compact graph-based framework, named HADA, which can combine pretrained models to produce a better result, rather than building from scratch. First, we created a graph structure in which the nodes were the features extracted from the pretrained models and the edges connecting them. The graph structure was employed to capture and fuse the information from every pretrained model with each other. Then a graph neural network was applied to update the connection between the nodes to get the representative embedding vector for an image and text. Finally, we used the cosine similarity to match images with their relevant texts and vice versa to ensure a low inference time. Our experiments showed that, although HADA contained a tiny number of trainable parameters, it could increase baseline performance by more than 3.6% in terms of evaluation metrics in the Flickr30k dataset. Additionally, the proposed model did not train on any external dataset and did not require many GPUs but only 1 to train due to its small number of parameters. The source code is available at https://github.com/m2man/HADA.


title: Analysis of Interior Rubble Void Spaces at Champlain Towers South Collapse

abstract: The 2021 Champlain Towers South Condominiums collapse in Surfside, Florida, resulted 98 deaths. Nine people are thought to have survived the initial collapse, and might have been rescued if rescue workers could have located them. Perhaps, if rescue workers had been able to use robots to search the interior of the rubble pile, outcomes might have been better. An improved understanding of the environment in which a robot would have to operate to be able to search the interior of a rubble pile would help roboticists develop better suited robotic platforms and control strategies. To this end, this work offers an approach to characterize and visualize the interior of a rubble pile and conduct a preliminary analysis of the occurrence of voids. Specifically, the analysis makes opportunistic use of four days of aerial imagery gathered from responders at Surfside to create a 3D volumetric aggregated model of the collapse in order to identify and characterize void spaces in the interior of the rubble. The preliminary results confirm expectations of small number and scale of these interior voids. The results can inform better selection and control of existing robots for disaster response, aid in determining the design specifications (specifically scale and form factor), and improve control of future robotic platforms developed for search operations in rubble.


title: An Oxygen Target for (Anti)neutrinos

abstract: We discuss a method to obtain an effective oxygen target within a low-density detector allowing an accurate characterization of the various event topologies in $\nu (\bar \nu)$-oxygen interactions. Results can be of interest for long-baseline neutrino oscillation experiments utilizing water targets. In particular, the combination of both oxygen and hydrogen targets within the same detector can provide in-situ measurements of nuclear effects and of the (anti)neutrino flux, which are the leading sources of systematic uncertainties in long-baseline oscillation analyses. These measurements can also provide useful information about the nuclear modifications of bound nucleons, as well as about the isospin symmetry in nucleons and nuclei.


title: Fast Persistent Homology Computation for Functions on $\mathbb{R}$

abstract: 0-dimensional persistent homology is known, from a computational point of view, as the easy case. Indeed, given a list of $n$ edges in non-decreasing order of filtration value, one only needs a union-find data structure to keep track of the connected components and we get the persistence diagram in time $O(n\alpha(n))$. The running time is thus usually dominated by sorting the edges in $\Theta(n\log(n))$. A little-known fact is that, in the particularly simple case of studying the sublevel sets of a piecewise-linear function on $\mathbb{R}$ or $\mathbb{S}^1$, persistence can actually be computed in linear time. This note presents a simple algorithm that achieves this complexity and an extension to image persistence. An implementation is available in Gudhi.


title: Switchable Lightweight Anti-symmetric Processing (SLAP) with CNN Outspeeds Data Augmentation by Smaller Sample -- Application in Gomoku Reinforcement Learning

abstract: To replace data augmentation, this paper proposed a method called SLAP to intensify experience to speed up machine learning and reduce the sample size. SLAP is a model-independent protocol/function to produce the same output given different transformation variants. SLAP improved the convergence speed of convolutional neural network learning by 83% in the experiments with Gomoku game states, with only one eighth of the sample size compared with data augmentation. In reinforcement learning for Gomoku, using AlphaGo Zero/AlphaZero algorithm with data augmentation as baseline, SLAP reduced the number of training samples by a factor of 8 and achieved similar winning rate against the same evaluator, but it was not yet evident that it could speed up reinforcement learning. The benefits should at least apply to domains that are invariant to symmetry or certain transformations. As future work, SLAP may aid more explainable learning and transfer learning for domains that are not invariant to symmetry, as a small step towards artificial general intelligence.


title: Equitable Network-Aware Decarbonization of Residential Heating at City Scale

abstract: Residential heating, primarily powered by natural gas, accounts for a significant portion of residential sector energy use and carbon emissions in many parts of the world. Hence, there is a push towards decarbonizing residential heating by transitioning to energy-efficient heat pumps powered by an increasingly greener and less carbon-intensive electric grid. However, such a transition will add additional load to the electric grid triggering infrastructure upgrades, and subsequently erode the customer base using the gas distribution network. Utilities want to guide these transition efforts to ensure a phased decommissioning of the gas network and deferred electric grid infrastructure upgrades while achieving carbon reduction goals. To facilitate such a transition, we present a network-aware optimization framework for decarbonizing residential heating at city scale with an objective to maximize carbon reduction under budgetary constraints. Our approach operates on a graph representation of the gas network topology to compute the cost of transitioning and select neighborhoods for transition. We further extend our approach to explicitly incorporate equity and ensure an equitable distribution of benefits across different socioeconomic groups. We apply our framework to a city in the New England region of the U.S., using real-world gas usage, electric usage, and grid infrastructure data. We show that our network-aware strategy achieves 55% higher carbon reductions than prior network-oblivious work under the same budget. Our equity-aware strategy achieves an equitable outcome while preserving the carbon reduction benefits of the network-aware strategy.


title: LSDM: Long-Short Diffeomorphic Motion for Weakly-Supervised Ultrasound Landmark Tracking

abstract: Accurate tracking of an anatomical landmark over time has been of high interests for disease assessment such as minimally invasive surgery and tumor radiation therapy. Ultrasound imaging is a promising modality benefiting from low-cost and real-time acquisition. However, generating a precise landmark tracklet is very challenging, as attempts can be easily distorted by different interference such as landmark deformation, visual ambiguity and partial observation. In this paper, we propose a long-short diffeomorphic motion network, which is a multi-task framework with a learnable deformation prior to search for the plausible deformation of landmark. Specifically, we design a novel diffeomorphism representation in both long and short temporal domains for delineating motion margins and reducing long-term cumulative tracking errors. To further mitigate local anatomical ambiguity, we propose an expectation maximisation motion alignment module to iteratively optimize both long and short deformation, aligning to the same directional and spatial representation. The proposed multi-task system can be trained in a weakly-supervised manner, which only requires few landmark annotations for tracking and zero annotation for long-short deformation learning. We conduct extensive experiments on two ultrasound landmark tracking datasets. Experimental results show that our proposed method can achieve better or competitive landmark tracking performance compared with other state-of-the-art tracking methods, with a strong generalization capability across different scanner types and different ultrasound modalities.


title: Orthogonal polynomials in weighted Bergman spaces

abstract: Let $w$ be a weight on the unit disk $\mathbb{D}$ having the form \[w(z)=|v(z)|^2\prod_{k=1}^s\left|\frac{z-a_k}{1-z\overline{a}_k}\right|^{m_k}\,,\quad m_k>-2,\ |a_k|<1,\] where $v$ is analytic and free of zeros in $\overline{\mathbb{D}}$, and let $(p_n)_{n=0}^\infty$ be the sequence of polynomials ($p_n$ of degree $n$) orthonormal over $\mathbb{D}$ with respect to $w$. We give an integral representation for $p_n$ from which it is in principle possible to derive its asymptotic behavior as $n\to\infty$ at every point $z$ of the complex plane, the asymptotic analysis of the integral being primarily dependent on the nature of the first singularities encountered by the function $v(z)^{-1}\prod_{k=1}^s(1-z\overline{a}_k)^{-1}$.


title: Charge transfer due to defects in hexagonal boron nitride/graphene heterostructures: an ab initio study

abstract: Using density functional theory (DFT), we study charge transfer between hexagonal boron nitride (h-BN) point defects and graphene in h-BN/graphene heterostructures for a range of intrinsic defects -- nitrogen vacancy, boron vacancy, nitrogen antisite and boron antisite. We show that traditional methods that calculate charge transfer by spatial discrimination of charge to different atoms suffer from the misallocation of charge and introduce an alternative method that relies on the integration of the density of states. We also show that DFT calculations of charge transfer have cell size dependencies due to a change in the density of states in the vicinity of the defect levels. Our results indicate that the nitrogen and boron anitsites do not participate in charge transfer, whereas the nitrogen and boron vacancies experience the transfer of a whole electron. Additionally, we show that a change in the geometry of a defect corresponds to a change in the charge state of the defect. The results of our study will be invaluable for a wide variety of device applications that involve charge transfer between h-BN defects and graphene in h-BN/graphene heterostructures, while our methodology can be feasibly extended to a wide range of point defects and heterostructures.


title: Artificial Intelligence Generated Coins for Size Comparison

abstract: Authors of scientific articles use coins in photographs as a size reference for objects. For this purpose, coins are placed next to objects when taking the photo. In this letter we propose a novel method that uses artificial intelligence (AI) generated images of coins to provide a size reference in photos. The newest generation is able to quickly generate realistic high-quality images from textual descriptions. With the proposed method no physical coin is required while taking photos. Coins can be added to photos that contain none. Furthermore, we show how the coin motif can be matched to the object.


title: Semantic Web Enabled Geographic Question Answering Framework: GeoTR

abstract: With the considerable growth of linked data, researchers have focused on how to increase the availability of semantic web technologies to provide practical usages for real life systems. Question answering systems are an example of real-life systems that communicate directly with end users, understand user intention and generate answers. End users do not care about the structural query language or the vocabulary of the knowledge base where the point of a problem arises. In this study, a question answering framework that converts Turkish natural language input into SPARQL queries in the geographical domain is proposed. Additionally, a novel Turkish ontology, which covers a 10th grade geography lesson named Spatial Synthesis Turkey, has been developed to be used as a linked data provider. Moreover, a gap in the literature on Turkish question answering systems, which utilizes linked data in the geographical domain, is addressed. A hybrid system architecture that combines natural language processing techniques with linked data technologies to generate answers is also proposed. Further related research areas are suggested.


title: Cache-Aided $K$-User Broadcast Channels with State Information at Receivers

abstract: We study a $K$-user coded-caching broadcast problem in a joint source-channel coding framework. The transmitter observes a database of files that are being generated at a certain rate per channel use, and each user has a cache, which can store a fixed fraction of the generated symbols. In the delivery phase, the transmitter broadcasts a message so that the users can decode their desired files using the received signal and their cache content. The communication between the transmitter and the receivers happens over a (deterministic) \textit{time-varying} erasure broadcast channel, and the channel state information is only available to the users. We characterize the maximum achievable source rate for the $2$-user and the degraded $K$-user problems. We provide an upper bound for any caching strategy's achievable source rates. Finally, we present a linear programming formulation to show that the upper bound is not a sharp characterization. Closing the gap between the achievable rate and the optimum rate remains open.


title: A Point-of-Care Biosensor for Rapid Detection and Differentiation of COVID-19 Virus (SARS-CoV-2) and Influenza Virus Using Subwavelength Grating Micro-ring Resonator

abstract: In the context of continued spread of coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 and the emergence of new variants, the demand for rapid, accurate, and frequent detection is increasing. Besides, the new predominant strain, Omicron variant, manifests more similar clinical features to those of other common respiratory infections. The concurrent detection of multiple potential pathogens helps distinguish SARS-CoV-2 infection from other diseases with overlapping symptoms, which is significant for patients to receive tailored treatment and containing the outbreak. Here, we report a lab-on-a-chip biosensing platform for SARS-CoV-2 detection based on subwavelength grating micro-ring resonator. The sensing surface is functionalized by specific antibody against SARS-CoV-2 spike protein, which could produce redshifts of resonant peaks by antigen-antibody combination, thus achieving quantitative detection. Additionally, the sensor chip is integrated with a microfluidic chip with an anti-backflow Y-shaped structure that enables the concurrent detection of two analytes. In this study, we realized the detection and differentiation of COVID-19 and influenza A H1N1. Experimental results show that the limit of detection of our device reaches 100 fg/mL (1.31 fM) within 15 min detecting time, and cross-reactivity tests manifest the specificity of the optical diagnostic assay. Further, the integrated packaging and streamlined workflow facilitate its use for clinical applications. Thus, the biosensing platform offers a promising solution to achieve ultrasensitive, selective, multiplexed, and quantitative point-of-care detection of COVID-19.


title: Traveling wave enantio-selective electron paramagnetic resonance

abstract: We propose a novel method for enantio-selective electron paramagnetic resonance spectroscopy based on magneto-chiral anisotropy. We calculate the strength of this effect and propose a dedicated interferometer setup for its observation.


title: Techniques for Measuring Parallax and Proper Motion with VLBI

abstract: Astrometry at centimeter wavelengths using Very Long Baseline Interferometry is approaching accuracies of ~1 uas for the angle between a target and a calibrator source separated by <1 degree on the sky. The BeSSeL Survey and the Japanese VERA project are using this to map the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with massive, young stars. This paper outlines how micro-arcsecond astrometry is done, including details regarding the scheduling of observations, calibration of data, and measuring positions.


title: Inferring Needless Write Memory Accesses on Ethereum Bytecode (Extended Version)

abstract: Efficiency is a fundamental property of any type of program, but it is even more so in the context of the programs executing on the blockchain (known as smart contracts). This is because optimizing smart contracts has direct consequences on reducing the costs of deploying and executing the contracts, as there are fees to pay related to their bytes-size and to their resource consumption (called gas). Optimizing memory usage is considered a challenging problem that, among other things, requires a precise inference of the memory locations being accessed. This is also the case for the Ethereum Virtual Machine (EVM) bytecode generated by the most-widely used compiler, \texttt{solc}, whose rather unconventional and low-level memory usage challenges automated reasoning. This paper presents a static analysis, developed at the level of the EVM bytecode generated by \texttt{solc}, that infers write memory accesses that are needless and thus can be safely removed. The application of our implementation on more than 19,000 real smart contracts has detected about 6,200 needless write accesses in less than 4 hours. Interestingly, many of these writes were involved in memory usage patterns generated by \texttt{solc} that can be greatly optimized by removing entire blocks of bytecodes. To the best of our knowledge, existing optimization tools cannot infer such needless write accesses, and hence cannot detect these inefficiencies that affect both the deployment and the execution costs of Ethereum smart contracts.


title: High-Order Mixed Finite Element Variable Eddington Factor Methods

abstract: We apply high-order mixed finite element discretization techniques and their associated preconditioned iterative solvers to the Variable Eddington Factor (VEF) equations in two spatial dimensions. The mixed finite element VEF discretizations are coupled to a high-order Discontinuous Galerkin (DG) discretization of the Discrete Ordinates transport equation to form effective linear transport algorithms that are compatible with high-order (curved) meshes. This combination of VEF and transport discretizations is motivated by the use of high-order mixed finite element methods in hydrodynamics calculations at the Lawrence Livermore National Laboratory. Due to the mathematical structure of the VEF equations, the standard Raviart Thomas (RT) mixed finite elements cannot be used to approximate the vector variable in the VEF equations. Instead, we investigate three alternatives based on the use of continuous finite elements for each vector component, a non-conforming RT approach where DG-like techniques are used, and a hybridized RT method. We present numerical results that demonstrate high-order accuracy, compatibility with curved meshes, and robust and efficient convergence in iteratively solving the coupled transport-VEF system and in the preconditioned linear solvers used to invert the discretized VEF equations.


title: Difference equations and Omega functions

abstract: We introduce Omega functions that generalize Euler Gamma functions and study the functional difference equation they satisfy. Under a natural exponential growth condition, the vector space of meromorphic solutions of the functional equation is finite dimensional. We construct a basis of the space of solutions composed by Omega functions. Omega functions are defined as exponential periods. They have a meromorphic extension to the complex plane of order $1$ with simple poles at negative integers. They are characterized by their growth property on vertical strips and their functional equation. This generalizes Wielandt's characterization of Euler Gamma function. We also introduce Incomplete Omega functions that play an important role in the proofs.


title: Pragmatic Estimation of Sample Size for Number of Interviews for PRO development in the 2009 FDA PRO guidance

abstract: PROs developed de novo, using the FDA guidance may involve structured patient interviews or focus groups. Qualitative Research is a methodology for eliciting and coding interviews and produces concepts or themes. These concepts are used to develop items in a PRO for use as an endpoint in Clinical trials. A convention in the field is that interviews and code/concept elicitation are considered complete when subsequent interviews produces "no new concepts" -termed "saturation". FDA reviewers frequently challenge PRO developers whether there are sufficient patient interviews to confirm that saturation is achieved after occurrence of zero new concepts. Several authors have reported that concrete criteria are need for confirming that saturation is achieved (Francis 2010, Mason 2010, Marshall 2013). I provide statistical methodology for confirming saturation, suitable for review by a regulatory authority. Type I error for saturation, may occur if further interviews elicited more concepts after first occurrence of saturation. I use published data set on code elicitation (Guest, 2006) to demonstrate that saturation may occur more than once in a sequence of interviews. I provide a statistical definition for saturation in qualitative research, that addresses regulatory concerns for PRO's developed for use as a clinical trial endpoint in a regulatory submission.


title: NarrowBERT: Accelerating Masked Language Model Pretraining and Inference

abstract: Large-scale language model pretraining is a very successful form of self-supervised learning in natural language processing, but it is increasingly expensive to perform as the models and pretraining corpora have become larger over time. We propose NarrowBERT, a modified transformer encoder that increases the throughput for masked language model pretraining by more than $2\times$. NarrowBERT sparsifies the transformer model such that the self-attention queries and feedforward layers only operate on the masked tokens of each sentence during pretraining, rather than all of the tokens as with the usual transformer encoder. We also show that NarrowBERT increases the throughput at inference time by as much as $3.5\times$ with minimal (or no) performance degradation on sentence encoding tasks like MNLI. Finally, we examine the performance of NarrowBERT on the IMDB and Amazon reviews classification and CoNLL NER tasks and show that it is also comparable to standard BERT performance.


title: Adjoint-Based Estimation of Sensitivity of Clinical Measures to Boundary Conditions for Arteries

abstract: The use of adjoint solvers is considered in order to obtain the sensitivity of clinical measures in aneurysms to incomplete (or unknown) boundary conditions and/or geometry. It is shown that these techniques offer interesting theoretical insights and viable computational tools to obtain these sensitivities.

categories: physics.med-ph math.OC physics.bio-ph physics.flu-dyn


title: Behaviors of pairs of dimensions and depths of edge ideals

abstract: Edge ideals of finite simple graphs $G$ on $n$ vertices are the ideals $I(G)$ of the polynomial ring $S$ in $n$ variables generated by the quadratic monomials associated with the edges of $G$. In this paper, we consider the possible pairs of dimensions and depths of $S/I(G)$ for connected graphs with a fixed number of vertices. We discuss such pairs in the case where dimension is relatively large. As a corollary, we completely determine the pairs for connected graphs with small number of vertices. We also study the possible pairs for connected chordal graphs.


title: Analyzing Inexact Hypergradients for Bilevel Learning

abstract: Estimating hyperparameters has been a long-standing problem in machine learning. We consider the case where the task at hand is modeled as the solution to an optimization problem. Here the exact gradient with respect to the hyperparameters cannot be feasibly computed and approximate strategies are required. We introduce a unified framework for computing hypergradients that generalizes existing methods based on the implicit function theorem and automatic differentiation/backpropagation, showing that these two seemingly disparate approaches are actually tightly connected. Our framework is extremely flexible, allowing its subproblems to be solved with any suitable method, to any degree of accuracy. We derive a priori and computable a posteriori error bounds for all our methods, and numerically show that our a posteriori bounds are usually more accurate. Our numerical results also show that, surprisingly, for efficient bilevel optimization, the choice of hypergradient algorithm is at least as important as the choice of lower-level solver.


title: The role of hosting capacity study in power system advancements: a review

abstract: The fast depletion of conventional energy sources due to increased energy demands and environmental concern has motivated power utilities to integrate more renewable energy sources into their power systems. Due to the intermittent nature and low or non-existent inertial response of these sources, a high penetration of RESs can lead to various issues in the operation of power systems such as oscillations in power system voltage and frequency, increased harmonic distortion, failure of protective equipment, overloading of transformers and feeders, and increased line losses. Such problems arise when the hosting capacity, defined as the maximum RES capacity that can be installed without having any technical and operational problems, of the network exceeds its limit due to the increased integration of RESs to the existing network. This paper reviews the progress made in HC assessment and enhancement of electrical networks research and development since its inception. Attempts are also made to highlight the current and future issues involved in HC technology for the development of an affordable, in-exhaustive, clean and reliable power supply for longer term benefits


title: Thermal Uhlmann phase in a locally driven two-spin system

abstract: We study the geometric Uhlmann phase of mixed states at finite temperature in a system of two coupled spin-$\frac 1 2$ particles driven by a magnetic field applied to one of the spins. In the parameter space of temperature and coupling, we show the emergence of two topological Uhlmann phase transitions when the magnetic field evolves around the equator, where a winding number can characterize each temperature range. For small couplings, the width of the temperature gap of the non-trivial phase is roughly the critical temperature $T_c$ of one-dimensional fermion systems with two-band Hamiltonians. The first phase transition in the low-temperature regime and small values of the coupling corresponds to the peak of the \textit{Schottky anomaly} of the heat capacity, typical of a two-level system in solid-state physics involving the ground and first excited states. The second phase transition occurs at temperatures very close to the second maximum of the heat capacity associated with a multilevel system. We also derive analytical expressions for the thermal Uhlmann phase for both subsystems, showing that they exhibit phase transitions. In the driven subsystem, for minimal $g$, a topological phase transition phase appears at $T_c$ again. However, for larger values of $g$, the transitions occur at lower temperature values, and they disappear when the coupling reaches a critical value $g_c$. The latter is not the case for the undriven subsystem, where at low temperatures, a single phase transition occurs at $g_c$. Nevertheless, as the temperature rises, we demonstrate the emergence of two phase transitions defining a coupling gap, where the phase is non-trivial and vanishes as the temperature reaches a critical value.


title: RAD-Sim: Rapid Architecture Exploration for Novel Reconfigurable Acceleration Devices

abstract: With the continued growth in field-programmable gate array (FPGA) capacity and their incorporation into new environments such as datacenters, we have witnessed the introduction of a new class of reconfigurable acceleration devices (RADs) that go beyond conventional FPGA architectures. These devices combine a reconfigurable fabric with coarse-grained domain-specialized accelerator blocks all connected via a high-performance packet-switched network-on-chip (NoC) for efficient system-wide communication. However, we lack the tools necessary to efficiently explore the huge design space for RADs, study the complex interactions between their different components and evaluate various combinations of design choices. In this work, we develop RAD-Sim, a cycle-level architecture simulator that allows rapid application-driven exploration of the design space of novel RADs. To showcase the capabilities of RADSim, we map and simulate a state-of-the-art deep learning (DL) inference overlay on a RAD instance incorporating an FPGA fabric and a complex of hard matrix-vector multiplication engines, communicating over a system-wide NoC. Through this example, we show how RAD-Sim can help architects quantify the effect of changing specific architecture parameters on end-to-end application performance.


title: Closely estimating the entropy of sparse graph models

abstract: We introduce an algorithm for estimating the entropy of pairwise, probabilistic graph models by leveraging bridges between social communities and an accurate entropy estimator on sparse samples. We propose using a measure of investment from the sociological literature, Burt's structural constraint, as a heuristic for identifying bridges that partition a graph into conditionally independent components. We combine this heuristic with the Nemenman-Shafee-Bialek entropy estimator to obtain a faster and more accurate estimator. We demonstrate it on the pairwise maximum entropy, or Ising, models of judicial voting, to improve na\"ive entropy estimates. We use our algorithm to estimate the partition function closely, which we then apply to the problem of model selection, where estimating the likelihood is difficult. This serves as an improvement over existing methods that rely on point correlation functions to test fit can be extended to other graph models with a straightforward modification of the open-source implementation.