title: Learning the Step-size Policy for the Limited-Memory Broyden-Fletcher-Goldfarb-Shanno Algorithm

abstract: We consider the problem of how to learn a step-size policy for the Limited-Memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm. This is a limited computational memory quasi-Newton method widely used for deterministic unconstrained optimization but currently avoided in large-scale problems for requiring step sizes to be provided at each iteration. Existing methodologies for the step size selection for L-BFGS use heuristic tuning of design parameters and massive re-evaluations of the objective function and gradient to find appropriate step-lengths. We propose a neural network architecture with local information of the current iterate as the input. The step-length policy is learned from data of similar optimization problems, avoids additional evaluations of the objective function, and guarantees that the output step remains inside a pre-defined interval. The corresponding training procedure is formulated as a stochastic optimization problem using the backpropagation through time algorithm. The performance of the proposed method is evaluated on the training of classifiers for the MNIST database for handwritten digits and for CIFAR-10. The results show that the proposed algorithm outperforms heuristically tuned optimizers such as ADAM, RMSprop, L-BFGS with a backtracking line search, and L-BFGS with a constant step size. The numerical results also show that a learned policy can be used as a warm-start to train new policies for different problems after a few additional training steps, highlighting its potential use in multiple large-scale optimization problems.


title: Use Cases of Quantum Optimization for Finance

abstract: In this paper we briefly review two recent use-cases of quantum optimization algorithms applied to hard problems in finance and economy. Specifically, we discuss the prediction of financial crashes as well as dynamic portfolio optimization. We comment on the different types of quantum strategies to carry on these optimizations, such as those based on quantum annealers, universal gate-based quantum processors, and quantum-inspired Tensor Networks.


title: The Fourier signatures of memristive hysteresis

abstract: While resistors with memory, sometimes called memristive elements (such as ReRAM cells), are often studied under conditions of periodic driving, little attention has been paid to the Fourier features of their memory response (hysteresis). Here we demonstrate experimentally that the hysteresis of memristive systems can be unambiguously distinguished from the linear or non-linear response of systems without hysteresis by the values of certain Fourier series coefficients. We also show that the Fourier series convergence depends on driving conditions, and introduce a measure of hysteresis. These results may be used to quantify the memory content of resistive memories, and tune their Fourier spectrum according to the excitation signal.


title: On almost quasi-negative holomorphic sectional curvature

abstract: A recent celebrated theorem of Diverio-Trapani and Wu-Yau states that a compact K\"ahler manifold admitting a K\"ahler metric of quasi-negative holomorphic sectional curvature has an ample canonical line bundle, confirming a conjecture of Yau. In this paper we shall consider a natural notion of almost quasi-negative holomorphic sectional curvature and extend this theorem to compact K\"ahler manifolds of almost quasi-negative holomorphic sectional curvature. We also obtain a gap-type theorem for the inequality $\int_Xc_1(K_X)^n>0$ in terms of the holomorphic sectional curvature. In the discussions, we introduce a capacity notion for the negative part of holomorphic sectional curvature, which plays a key role in studying the relation between the almost quasi-negative holomorphic sectional curvature and ampleness of the canonical line bundle.


title: A simulation environment for drone cinematography

abstract: In this paper, we present a workflow for the simulation of drone operations exploiting realistic background environments constructed within Unreal Engine 4 (UE4). Methods for environmental image capture, 3D reconstruction (photogrammetry) and the creation of foreground assets are presented along with a flexible and user-friendly simulation interface. Given the geographical location of the selected area and the camera parameters employed, the scanning strategy and its associated flight parameters are first determined for image capture. Source imagery can be extracted from virtual globe software or obtained through aerial photography of the scene (e.g. using drones). The latter case is clearly more time consuming but can provide enhanced detail, particularly where coverage of virtual globe software is limited. The captured images are then used to generate 3D background environment models employing photogrammetry software. The reconstructed 3D models are then imported into the simulation interface as background environment assets together with appropriate foreground object models as a basis for shot planning and rehearsal. The tool supports both free-flight and parameterisable standard shot types along with programmable scenarios associated with foreground assets and event dynamics. It also supports the exporting of flight plans. Camera shots can also be designed to provide suitable coverage of any landmarks which need to appear in-shot. This simulation tool will contribute to enhanced productivity, improved safety (awareness and mitigations for crowds and buildings), improved confidence of operators and directors and ultimately enhanced quality of viewer experience.


title: Accounts, Accountability and Agency for Safe and Ethical AI

abstract: We examine the problem of explainable AI (xAI) and explore what delivering xAI means in practice, particularly in contexts that involve formal or informal and ad-hoc collaboration where agency and accountability in decision-making are achieved and sustained interactionally. We use an example from an earlier study of collaborative decision-making in screening mammography and the difficulties users faced when trying to interpret the behavior of an AI tool to illustrate the challenges of delivering usable and effective xAI. We conclude by setting out a study programme for future research to help advance our understanding of xAI requirements for safe and ethical AI.


title: Differential interferometry of close binary of supermassive black holes in an elliptical configuration

abstract: The Very Large Telescope Interferometer (VLTI), and the Extremely Large Telescope (ELT) will be a robust astrophysics suite offering the opportunity of probing the structure and dynamics of CB-SMBH at high spectral and angular resolution. Here, we explore and illustrate the application of differential interferometry on unresolved the CB-SMBH systems in elliptical orbital configurations and a single SMBH with clouds in elliptical orbital motion. Photocenter displacements between each SMBH and regions in their disc-like broad line regions (BLR) appear as small interferometric differential phase variability. To investigate the application of interferometric phases for the detection of CB-SMBH systems, we simulate a series of differential interferometry signatures, based on our model comprising ensembles of clouds surrounding each of supermassive black hole in a CB-SMBH. Setting model to the parameters of a single SMBH with elliptical cloud motion, we also calculated a series of differential interferometry observables for this case. We found various deviations from the canonical S-shaped of CB-SMBH phases profile for elliptically configured CB-SMBH systems. The amplitude and specific shape of the interferometry observables depend on orbital configurations of the CB-SMBH system. We get distinctive results when considering antialigned angular momenta of cloud orbits regarding total CB-SMBH angular momentum. Some simulated spectral lines from our model closely resemble observations of Pa \alpha line got from near-infrared AGN surveys. We found differences between differential phases zoo of a single SMBH and CB-SMBH systems. The differential phases zoo for a single SMBH comprises deformed S shape. We also showed how their differential phase shape, amplitude, and slope evolve with various sets of cloud orbital parameters and observer position.


title: Gaussian Vector: An Efficient Solution for Facial Landmark Detection

abstract: Significant progress has been made in facial landmark detection with the development of Convolutional Neural Networks. The widely-used algorithms can be classified into coordinate regression methods and heatmap based methods. However, the former loses spatial information, resulting in poor performance while the latter suffers from large output size or high post-processing complexity. This paper proposes a new solution, Gaussian Vector, to preserve the spatial information as well as reduce the output size and simplify the post-processing. Our method provides novel vector supervision and introduces Band Pooling Module to convert heatmap into a pair of vectors for each landmark. This is a plug-and-play component which is simple and effective. Moreover, Beyond Box Strategy is proposed to handle the landmarks out of the face bounding box. We evaluate our method on 300W, COFW, WFLW and JD-landmark. That the results significantly surpass previous works demonstrates the effectiveness of our approach.


title: Deep learning algorithms for solving high dimensional nonlinear backward stochastic differential equations

abstract: In this work, we propose a new deep learning-based scheme for solving high dimensional nonlinear backward stochastic differential equations (BSDEs). The idea is to reformulate the problem as a global optimization, where the local loss functions are included. Essentially, we approximate the unknown solution of a BSDE using a deep neural network and its gradient with automatic differentiation. The approximations are performed by globally minimizing the quadratic local loss function defined at each time step, which always includes the terminal condition. This kind of loss functions are obtained by iterating the Euler discretization of the time integrals with the terminal condition. Our formulation can prompt the stochastic gradient descent algorithm not only to take the accuracy at each time layer into account, but also converge to a good local minima. In order to demonstrate performances of our algorithm, several high-dimensional nonlinear BSDEs including pricing problems in finance are provided.


title: New Revival Phenomena for Linear Integro-Differential Equations

abstract: We present and analyse a novel manifestation of the revival phenomenon for linear spatially periodic evolution equations, in the concrete case of three nonlocal equations that arise in water wave theory and are defined by convolution kernels. Revival in these cases is manifested in the form of dispersively quantised cusped solutions at rational times. We give an analytic description of this phenomenon, and present illustrative numerical simulations.


title: Temporal Logic of Minkowski Spacetime

abstract: We present the proof that the temporal logic of two-dimensional Minkowski spacetime is decidable, PSPACE-complete. The proof is based on a type of two-dimensional mosaic. Then we present the modification of the proof so as to work for slower-than-light signals. Finally, a subframe of the slower-than-light Minkowski frame is used to prove the new result that the temporal logic of real intervals with during as the accessibility relation is also PSPACE-complete.


title: Barrier methods for minimal submanifolds in the Gibbons-Hawking ansatz

abstract: We describe a barrier argument for compact minimal submanifolds in the multi-Eguchi-Hanson and in the multi-Taub-NUT spaces, which are hyperkaehler 4-manifolds given by the Gibbons-Hawking ansatz. This approach is used to obtain results towards a classification of compact minimal submanifolds in this setting. We also prove a converse of Tsai and Wang's result that relates the strong stability condition to the convexity of the distance function.


title: The Design of Tangible Digital Musical Instruments

abstract: Here we present guidelines that highlight the impact of haptic feedback upon the experiences of computer musicians using Digital Musical Instruments (DMIs). In this context, haptic feedback offers a tangible, bi-directional exchange between a musician and a DMI. We propose that by adhering to and exploring these guidelines the application of haptic feedback can enhance and augment the physical and affective experiences of a musician in interactions with these devices. It has been previously indicated that in the design of haptic DMIs, the experiences and expectations of a musician must be considered for the creation of tangible DMIs and that haptic feedback can be used to address the physical-digital divide that currently exists between users of such instruments.


title: Superdeterminism: A Guide for the Perplexed

abstract: Superdeterminism is presently the only known consistent description of nature that is local, deterministic, and can give rise to the observed correlations of quantum mechanics. I here want to explain what makes this approach promising and offer the reader some advice for how to avoid common pitfalls. In particular, I explain why superdeterminism is not a threat to science, is not necessarily finetuned, what the relevance of future input is, and what the open problems are.


title: Eisenstein series, p-adic modular functions, and overconvergence

abstract: Let $p$ be a prime $\ge 5$. We establish explicit rates of overconvergence for members of the "Eisenstein family", notably for the $p$-adic modular function $V(E_{(1,0)}^{\ast})/E_{(1,0)}^{\ast}$ ($V$ the $p$-adic Frobenius operator) that plays a pi\-votal role in Coleman's theory of $p$-adic families of modular forms. The proof goes via an in-depth analysis of rates of overconvergence of $p$-adic modular functions of form $V(E_k)/E_k$ where $E_k$ is the classical Eisenstein series of level $1$ and weight $k$ divisible by $p-1$. Under certain conditions, we extend the latter result to a vast generalization of a theorem of Coleman--Wan regarding the rate of overconvergence of $V(E_{p-1})/E_{p-1}$. We also comment on previous results in the literature. These include applications of our results for the primes $5$ and $7$.


title: Digital Musical Instrument Analysis: The Haptic Bowl

abstract: This experiment is a case study that applies a HCI-informed DMI Evaluation Framework. This framework applies existing HCI evaluation methods to the assessment of prototype Digital Musical Instruments (DMIs). The overall study will involve a three-part analysis - a description and categorisation of the device, a functionality evaluation that included an examination of usability and user experience, and finally an exploration of the device's effectiveness as a digital instrument. Here we present the findings of the first two parts of the framework, outlining the constituent components of the interface and testing the functionality of the device. The final stage of analysis will involve a longitudinal study and will be carried out in order to assess the musical affordances of the device.


title: A Toy Model for Local and Deterministic Wave-function Collapse

abstract: A local, deterministic toy model for quantum mechanics is introduced and discussed. It is demonstrated that, when averaged over the hidden variables, the model produces the same predictions as quantum mechanics. In the model considered here, the dynamics depends only on the settings of the measurement device at the detection time, not how those settings were chosen. As a consequence, the model is locally causal but violates statistical independence. We show that it is neither fine-tuned nor allows for superluminal signalling.


title: HCI Models for Digital Musical Instruments: Methodologies for Rigorous Testing of Digital Musical Instruments

abstract: Here we present an analysis of literature relating to the evaluation methodologies of Digital Musical Instruments (DMIs) derived from the field of Human-Computer Interaction (HCI). We then apply choice aspects from these existing evaluation models and apply them to an optimized evaluation for assessing new DMIs.


title: Multi-Step Adversarial Perturbations on Recommender Systems Embeddings

abstract: Recommender systems (RSs) have attained exceptional performance in learning users' preferences and helping them in finding the most suitable products. Recent advances in adversarial machine learning (AML) in the computer vision domain have raised interests in the security of state-of-the-art model-based recommenders. Recently, worrying deterioration of recommendation accuracy has been acknowledged on several state-of-the-art model-based recommenders (e.g., BPR-MF) when machine-learned adversarial perturbations contaminate model parameters. However, while the single-step fast gradient sign method (FGSM) is the most explored perturbation strategy, multi-step (iterative) perturbation strategies, that demonstrated higher efficacy in the computer vision domain, have been highly under-researched in recommendation tasks. In this work, inspired by the basic iterative method (BIM) and the projected gradient descent (PGD) strategies proposed in the CV domain, we adapt the multi-step strategies for the item recommendation task to study the possible weaknesses of embedding-based recommender models under minimal adversarial perturbations. Letting the magnitude of the perturbation be fixed, we illustrate the highest efficacy of the multi-step perturbation compared to the single-step one with extensive empirical evaluation on two widely adopted recommender datasets. Furthermore, we study the impact of structural dataset characteristics, i.e., sparsity, density, and size, on the performance degradation issued by presented perturbations to support RS designer in interpreting recommendation performance variation due to minimal variations of model parameters. Our implementation and datasets are available at https://anonymous.4open.science/r/9f27f909-93d5-4016-b01c-8976b8c14bc5/.


title: On the group of holomorphic automorphisms of model surfaces

abstract: It is proved that the group of holomorphic automorphisms of holomorphically homogeneous nondegenerate (finite Bloom-Graham type + holomorphic nondegenaracy) model surface Q is a subgroup of the group of birational automorphisms of the ambient space (Cremona group) with uniformly bounded degree. An estimate of the degree of the automorphisms in terms of the dimension of the ambient space is given (Theorem 4). It is proved that none of the conditions of the theorem can be omitted. We consider also the question about connectedness of the automorphism group of Q (Theorem 7).


title: Neighborhood Matters: Influence Maximization in Social Networks with Limited Access

abstract: Influence maximization (IM) aims at maximizing the spread of influence by offering discounts to influential users (called seeding). In many applications, due to user's privacy concern, overwhelming network scale etc., it is hard to target any user in the network as one wishes. Instead, only a small subset of users is initially accessible. Such access limitation would significantly impair the influence spread, since IM often relies on seeding high degree users, which are particularly rare in such a small subset due to the power-law structure of social networks. In this paper, we attempt to solve the limited IM in real-world scenarios by the adaptive approach with seeding and diffusion uncertainty considered. Specifically, we consider fine-grained discounts and assume users accept the discount probabilistically. The diffusion process is depicted by the independent cascade model. To overcome the access limitation, we prove the set-wise friendship paradox (FP) phenomenon that neighbors have higher degree in expectation, and propose a two-stage seeding model with the FP embedded, where neighbors are seeded. On this basis, for comparison we formulate the non-adaptive case and adaptive case, both proven to be NP-hard. In the non-adaptive case, discounts are allocated to users all at once. We show the monotonicity of influence spread w.r.t. discount allocation and design a two-stage coordinate descent framework to decide the discount allocation. In the adaptive case, users are sequentially seeded based on observations of existing seeding and diffusion results. We prove the adaptive submodularity and submodularity of the influence spread function in two stages. Then, a series of adaptive greedy algorithms are proposed with constant approximation ratio.


title: Inter- and intra-band Coulomb interactions between holes in silicon nanostructures

abstract: We present a full derivation of the interaction Hamiltonian for holes in silicon within the six-band envelope-function scheme, which appropriately describes the valence band close to the $\boldsymbol{\Gamma}$ point. The full structure of the single-hole eigenstates is taken into account, including the Bloch part. The scattering processes caused by the Coulomb interaction are shown to be both intraband and interband, the latter being mostly short-ranged. In the asymptotic long-range limit, the effective potential tends to the screened Coulomb potential, and becomes purely intraband, as assumed in previous models. We apply our model to compute the excitation spectra of two interacting holes in prototypical silicon quantum dots, taking into account different dielectric environments. It is shown that, in the highly screened regime, short-range interactions (both intra- and inter-band) can be very relevant, while they lose importance when there is no screening other than the one proper of the bulk silicon crystal. In the latter case, we predict the formation of hole Wigner molecules.


title: EGMM: an Evidential Version of the Gaussian Mixture Model for Clustering

abstract: The Gaussian mixture model (GMM) provides a simple yet principled framework for clustering, with properties suitable for statistical inference. In this paper, we propose a new model-based clustering algorithm, called EGMM (evidential GMM), in the theoretical framework of belief functions to better characterize cluster-membership uncertainty. With a mass function representing the cluster membership of each object, the evidential Gaussian mixture distribution composed of the components over the powerset of the desired clusters is proposed to model the entire dataset. The parameters in EGMM are estimated by a specially designed Expectation-Maximization (EM) algorithm. A validity index allowing automatic determination of the proper number of clusters is also provided. The proposed EGMM is as simple as the classical GMM, but can generate a more informative evidential partition for the considered dataset. The synthetic and real dataset experiments show that the proposed EGMM performs better than other representative clustering algorithms. Besides, its superiority is also demonstrated by an application to multi-modal brain image segmentation.


title: Space-Charge Limited Current from a Finite Emitter in Nano- and Microdiodes

abstract: We simulate numerically the classical charge dynamics in a microscopic, planar, vacuum diode with a finite emitter area and a finite number of electrons in the gap. We assume electrons are emitted under space-charge limited conditions with a fixed potential applied to the diode. The Coulomb interaction between all electrons is included using the method of molecular dynamics. We compare our results to the conventional two-dimensional Child-Langmuir and explain how it is limited in applicability for sub-micron diameter emitters. Finally, we offer some simple relations for understanding space-charge limited flow from very small emitters.


title: Wealth and Poverty: The Effect of Poverty on Communities

abstract: This paper analyzes the differences in poverty in high wealth communities and low wealth communities. We first discuss methods of measuring poverty and analyze the causes of individual poverty and poverty in the Bay Area. Three cases are considered regarding relative poverty. The first two cases involve neighborhoods in the Bay Area while the third case evaluates two neighborhoods within the city of San Jose, CA. We find that low wealth communities have more crime, more teen births, and more cost-burdened renters because of high concentrations of temporary and seasonal workers, extensive regulations on greenhouse gas emissions, minimum wage laws, and limited housing supply. In the conclusion, we review past attempts to alleviate the effects of poverty and give suggestions on how future policy can be influenced to eventually create a future free of poverty.


title: Laser-processing of grinded and mechanically abraded Nb-surfaces

abstract: The effect of pulsed laser polishing on rough niobium surfaces was investigated. We created different well-defined roughness profiles with standard emery papers and subsequently remelted random surface areas with a size of about 2x2 mm^2 with nanosecond laser pulses (wavelength of 1064 nm, pulse length of 10 ns). Pristine as well as laser-treated surfaces were investigated using optical profilometry and atomic force microscopy, and the surface topography was described by means of correlation functions. Uniformly rough and highly smooth surface geometries were achieved for fractals above and below 7 {\mu}m, respectively. Moreover, the behavior of foreign particles during the laser processing was investigated in detail. The polishing procedure was also monitored point by point by detecting electrical signals, i.e., sample charging, which resulted from the intense laser illumination. The measured electrical charges were found to be correlated with the local surface texture. Thus, regions with initially high roughness profiles and regions with extensive laser-induced defects could be directly identified from the detected electrical signals.


title: Bitcoin and its impact on the economy

abstract: The purpose of this paper is to review the concept of cryptocurrencies in our economy. First, Bitcoin and alternative cryptocurrencies' histories are analyzed. We then study the implementation of Bitcoin in the airline and real estate industries. Our study finds that many Bitcoin companies partner with airlines in order to decrease processing times, to provide ease of access for spending in international airports, and to reduce fees on foreign exchanges for fuel expenses, maintenance, and flight operations. Bitcoin transactions have occurred in the real estate industry, but many businesses are concerned with Bitcoin's potential interference with the U.S. government and its high volatility. As Bitcoin's price has been growing rapidly, we assessed Bitcoin's real value; Bitcoin derives value from its scarcity, utility, and public trust. In the conclusion, we discuss Bitcoin's future and conclude that Bitcoin may change from a short-term profit investment to a more steady industry as we identify Bitcoin with the "greater fool theory", and as the number of available Bitcoins to be mined dwindles and technology becomes more expensive.


title: A Possible kilonova powered by magnetic wind from a newborn black hole

abstract: The merger of binary neutron stars (NS-NS) as the progenitor of short Gamma-ray bursts (GRBs) has been confirmed by the discovery of the association of the gravitational wave (GW) event GW170817 with GRB 170817A. However, the merger product of binary NS remains an open question. An X-ray plateau followed by a steep decay ("internal plateau") has been found in some short GRBs, implying that a supra-massive magnetar operates as the merger remnant and then collapses into a newborn black hole (BH) at the end of the plateau. X-ray bump or second-plateau following the "internal plateau" is considered as the expected signature from the fallback accretion onto this newborn BH through Blandford-Znajek mechanism (BZ). At the same time, a nearly isotropic wind driven by Blandford-Paynemechanism (BP) from the newborn BH's disk can produce a bright kilonova. Therefore, the bright kilonova observation for a short GRB with "internal plateau" (and followed by X-ray bump or second-plateau) provides further evidence for this scenario. In this paper, we find that GRB 160821B is a candidate of such a case, and the kilonova emission of GRB 160821B is possibly powered by the BP wind from a newborn BH. Future GW detection of GRB 160821B-like events may provide further support to this scenario, enable us to investigate the properties of the magnetar and the newborn BH, and constrain the equation of state of neutron stars.


title: Weighted Sum-Rate Maximization for Multi-IRS-assisted Full-Duplex Systems with Hardware Impairments

abstract: Smart and reconfigurable wireless communication environments can be established by exploiting well-designed intelligent reflecting surfaces (IRSs) to shape the communication channels. In this paper, we investigate how multiple IRSs affect the performance of multi-user full-duplex communication systems under hardware impairment at each node, wherein the base station (BS) and the uplink users are subject to maximum transmission power constraints. Firstly, the uplink-downlink system weighted sum-rate (SWSR) is derived which serves as a system performance metric. Then, we formulate the resource allocation design for the maximization of SWSR as an optimization problem which jointly optimizes the beamforming and the combining vectors at the BS, the transmit powers of the uplink users, and the phase shifts of multiple IRSs. Since the SWSR optimization problem is non-convex, an efficient iterative alternating approach is proposed to obtain a suboptimal solution for the design problem considered and its complexity is also discussed. In particular, we firstly reformulate the main problem into an equivalent weighted minimum mean-square-error form and then transform it into several convex sub-problems which can be analytically solved for given phase shifts. Then, the IRSs phases are optimized via a gradient ascent-based algorithm. Finally, numerical results are presented to clarify how multiple IRSs enhance the performance metric under hardware impairment.


title: Electrostatic and hydrophobic patches on Amyloid-\b{eta} oligomers govern their fractal self-assembly: Implication to proteins

abstract: We present a generic model to describe the fractal self-assembly of proteins in terms of electrostatic and hydrophobic interactions. The predictions of the model were correlated with the simulated fractals obtained using patchy diffusion-limited aggregation, and the experimentally observed Amyloid-\b{eta} (A\b{eta}) fractal self-assembly using confocal microscopy. The molecular docking was used to determine the properties of the patches on A\b{eta} oligomers. Similar patch properties were used to design the particles for simulation. In agreement with the model predictions, the free energy of formation for the simulated fractal self-assembly was proportional to its fractal dimension; moreover, their morphologies were similar to the fractal morphologies of A\b{eta} observed experimentally at different pH.


title: The Hardy-Littlewood property and maximal operators associated with the inverse Gauss measure

abstract: In this paper we characterize the Banach lattices with the Hardy-Littlewood property by using maximal operators defined by semigroups of operators associated with the inverse Gauss measure.


title: End-to-End Training of CNN Ensembles for Person Re-Identification

abstract: We propose an end-to-end ensemble method for person re-identification (ReID) to address the problem of overfitting in discriminative models. These models are known to converge easily, but they are biased to the training data in general and may produce a high model variance, which is known as overfitting. The ReID task is more prone to this problem due to the large discrepancy between training and test distributions. To address this problem, our proposed ensemble learning framework produces several diverse and accurate base learners in a single DenseNet. Since most of the costly dense blocks are shared, our method is computationally efficient, which makes it favorable compared to the conventional ensemble models. Experiments on several benchmark datasets demonstrate that our method achieves state-of-the-art results. Noticeable performance improvements, especially on relatively small datasets, indicate that the proposed method deals with the overfitting problem effectively.


title: A Variational Information Bottleneck Based Method to Compress Sequential Networks for Human Action Recognition

abstract: In the last few years, compression of deep neural networks has become an important strand of machine learning and computer vision research. Deep models require sizeable computational complexity and storage, when used for instance for Human Action Recognition (HAR) from videos, making them unsuitable to be deployed on edge devices. In this paper, we address this issue and propose a method to effectively compress Recurrent Neural Networks (RNNs) such as Gated Recurrent Units (GRUs) and Long-Short-Term-Memory Units (LSTMs) that are used for HAR. We use a Variational Information Bottleneck (VIB) theory-based pruning approach to limit the information flow through the sequential cells of RNNs to a small subset. Further, we combine our pruning method with a specific group-lasso regularization technique that significantly improves compression. The proposed techniques reduce model parameters and memory footprint from latent representations, with little or no reduction in the validation accuracy while increasing the inference speed several-fold. We perform experiments on the three widely used Action Recognition datasets, viz. UCF11, HMDB51, and UCF101, to validate our approach. It is shown that our method achieves over 70 times greater compression than the nearest competitor with comparable accuracy for the task of action recognition on UCF11.


title: Codes for Distributed Storage

abstract: This chapter deals with the topic of designing reliable and efficient codes for the storage and retrieval of large quantities of data over storage devices that are prone to failure. For long, the traditional objective has been one of ensuring reliability against data loss while minimizing storage overhead. More recently, a third concern has surfaced, namely of the need to efficiently recover from the failure of a single storage unit, corresponding to recovery from the erasure of a single code symbol. We explain here, how coding theory has evolved to tackle this fresh challenge.


title: A Geometry-Inspired Attack for Generating Natural Language Adversarial Examples

abstract: Generating adversarial examples for natural language is hard, as natural language consists of discrete symbols, and examples are often of variable lengths. In this paper, we propose a geometry-inspired attack for generating natural language adversarial examples. Our attack generates adversarial examples by iteratively approximating the decision boundary of Deep Neural Networks (DNNs). Experiments on two datasets with two different models show that our attack fools natural language models with high success rates, while only replacing a few words. Human evaluation shows that adversarial examples generated by our attack are hard for humans to recognize. Further experiments show that adversarial training can improve model robustness against our attack.


title: Period change of binary pulsars due to the accretion of dark matter particles

abstract: As binary systems move inside galaxies, they interact with the dark matter halo. This interaction leads to the accretion of dark matter particles inside binary components. The accretion of dark matter particles increases the mass of the binary components and then, the total mass of the system. Increased mass by this way can affect other physical parameters of the systems, like orbital periods of the systems. In this work, we estimated the period change of some known compact binary systems due to the accretion of dark matter particles inside them. Our conclustion is that, for WIMP particles with masses in the range $\simeq 10 \: GeV.c^{-2}$ and dark matter density as high as the dark matter density around the sun, period change due to accretion of dark matter partices is negligible compared to the period change due to the emission of gravitaional waves from the systems.

categories: astro-ph.HE astro-ph.GA astro-ph.SR hep-th


title: Computationally sound Bitcoin tokens

abstract: We propose a secure and efficient implementation of fungible tokens on Bitcoin. Our technique is based on a small extension of the Bitcoin script language, which allows the spending conditions in a transaction to depend on the neighbour transactions. We show that our implementation is computationally sound: that is, adversaries can make tokens diverge from their ideal functionality only with negligible probability.


title: Sustained oscillations in Interstellar chemistry models

abstract: Non-linear behavior in interstellar chemical models has been recognized for 25 years now. Different mechanisms account for the possibility of multiple fixed-points at steady state, characterized by the ionization degree of the gas. Chemical oscillations are also a natural behaviour of non-linear chemical models. We study under which conditions spontaneous sustained chemical oscillations are possible, and what kind of biffurcations lead to, or quench, the occurrence of such oscillations. Methods. The well known Ordinary Differential Equations (ODE) integrator VODE is used to explore initial conditions and parameter space in a gas phase chemical model of a dark interstellar cloud. We recall that the time evolution of the various chemical abundances under fixed temperature conditions depends on the density over cosmic ionization rate nH/{\zeta} ratio. We also report the occurrence of naturally sustained oscillations for a limited but well defined range of control parameters. The period of oscillations is within the range of characteristic time scales of interstellar processes and could lead to spectacular resonances in time dependent models. Reservoir species (C, CO, NH3, ...) oscillation amplitudes are generally less than a factor two. However, these amplitudes reach a factor ten to thousand for low abundance species, e.g. HCN, ND3, that may play a key role for diagnostic purposes.The mechanism responsible for oscillations is tightly linked to the chemistry of nitrogen, and requires long chains of reactions such as found in multi-deuteration processes.


title: Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $\gamma$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $\gamma$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $\gamma$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $\gamma$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $\gamma$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $\gamma$-ray cosmology.


title: Duality theory for generalized summing linear operators

abstract: Generalizing classical results of the theory of absolutely summing operators, in this paper we characterize the duals of a quite large class of Banach operator ideals defined or characterized by the transformation of vector-valued sequences.


title: Disentangling causal effects for hierarchical reinforcement learning

abstract: Exploration and credit assignment under sparse rewards are still challenging problems. We argue that these challenges arise in part due to the intrinsic rigidity of operating at the level of actions. Actions can precisely define how to perform an activity but are ill-suited to describe what activity to perform. Instead, causal effects are inherently composable and temporally abstract, making them ideal for descriptive tasks. By leveraging a hierarchy of causal effects, this study aims to expedite the learning of task-specific behavior and aid exploration. Borrowing counterfactual and normality measures from causal literature, we disentangle controllable effects from effects caused by other dynamics of the environment. We propose CEHRL, a hierarchical method that models the distribution of controllable effects using a Variational Autoencoder. This distribution is used by a high-level policy to 1) explore the environment via random effect exploration so that novel effects are continuously discovered and learned, and to 2) learn task-specific behavior by prioritizing the effects that maximize a given reward function. In comparison to exploring with random actions, experimental results show that random effect exploration is a more efficient mechanism and that by assigning credit to few effects rather than many actions, CEHRL learns tasks more rapidly.


title: Control of the exciton valley dynamics in van der Waals heterostructures

abstract: The exciton valley dynamics in van der Waals heterostructures with transition metal dichalcogenides monolayers is driven by the long-range exchange interaction between the electron and the hole in the exciton. It couples the states active in the opposite circular polarizations resulting in the longitudinal-transverse splitting of excitons propagating in the monolayer plane. Here we study theoretically the effect of the dielectric environment on the long-range exchange interaction and demonstrate how the encapsulation in the hexagonal boron nitride modifies the exciton longitudinal-transverse splitting. We calculate the exciton spin/valley polarization relaxation due to the long-range exchange interaction and demonstrate that the variation of the monolayer environment results in significant, up to five-fold, enhancement of the exciton valley polarization lifetime.


title: CPTM symmetry, closed time paths and cosmological constant problem in the formalism of extended manifold

abstract: The problem of the cosmological constant is considered in the formalism of an extended space-time consisting of the extended classical solution of Einstein equations. The different regions of the extended manifold are proposed to be related by the charge, parity, time and mass (CPTM) reversal symmetry applied with respect to the metric fields of the manifolds. There are interactions between the points of the extended manifold provided by scalar fields present separately in the different patches of the extended solution. The value of the constant is obtained equal to zero at the classical level due the mutual contribution of the fields in the vacuum energy, it's non-zero value is due the quantum interactions between the fields. There are few possible scenario for the actions of the fields are discussed. Each from the obtained variants is similar to the closed time path approach of non-equilibrium condensed matter physics and among these possibilities for the closed paths, there is a variant of the action equivalent to the formalism of Keldysh. Accordingly, we consider and shortly discuss the application of the proposed formalism to the problem of smallness of the cosmological constant and singularities problem.


title: Interstitial Carbon in bcc HfNbTiVZr high entropy alloy from first principles

abstract: The remarkable mechanical properties of high entropy alloys can be further improved by interstitial alloying. In this work we employ density functional theory calculations to study the solution energies of dilute carbon interstitial atoms in tetrahedral and octahedral sites in bcc HfNbTiVZr. Our results indicate that carbon interstitials in tetrahedral sites are unstable, and the preferred octahedral sites present a large spread in the energy of solution. The inclusion of carbon interstitials induces large structural relaxations with long-range effects. The effect of local chemical environment on the energy of solution is investigated by performing a local cluster expansion including studies of its correlation with the carbon atomic Voronoi volume. However, the spread in solution energetics can not be explained with a local environment analysis only pointing towards a complex, long-range influence of interstitial carbon in this alloy.


title: Direct band gap and strong Rashba effect in van der Waals heterostructures of InSe and Sb single layers

abstract: Van der Waals heterostructures formed by stacking different types of 2D materials are attracting increasing attention due to new emergent physical properties such as interlayer excitons. Recently synthesized atomically thin indium selenide (InSe) and antimony (Sb) individually exhibit interesting electronic properties such as high electron mobility in the former and high hole mobility in the latter. In this work, we present a first-principles investigation on the stability and electronic properties of ultrathin bilayer heterostructures composed of InSe and Sb single layers. The calculated electronic band structures reveal a direct band gap semiconducting nature of the InSe/Sb heterostructures independent of stacking pattern. Taking spin-orbit coupling into account, we find a large Rashba spin splitting at the bottom of conduction band, which originates from the atomic spin-orbit coupling with the symmetry breaking in the heterostructure. The strength of the Rashba spin splitting can be tuned by applying in-plane biaxial strain or an out-of-plane external electric field. The presence of large Rashba spin splitting together with a suitable band gap in InSe/Sb bilayer heterostructures make them promising candidates for spin field-effect transistor and optoelectronic device applications.


title: Expectigrad: Fast Stochastic Optimization with Robust Convergence Properties

abstract: Many popular adaptive gradient methods such as Adam and RMSProp rely on an exponential moving average (EMA) to normalize their stepsizes. While the EMA makes these methods highly responsive to new gradient information, recent research has shown that it also causes divergence on at least one convex optimization problem. We propose a novel method called Expectigrad, which adjusts stepsizes according to a per-component unweighted mean of all historical gradients and computes a bias-corrected momentum term jointly between the numerator and denominator. We prove that Expectigrad cannot diverge on every instance of the optimization problem known to cause Adam to diverge. We also establish a regret bound in the general stochastic nonconvex setting that suggests Expectigrad is less susceptible to gradient variance than existing methods are. Testing Expectigrad on several high-dimensional machine learning tasks, we find it often performs favorably to state-of-the-art methods with little hyperparameter tuning.


title: Actionet: An Interactive End-To-End Platform For Task-Based Data Collection And Augmentation In 3D Environment

abstract: The problem of task planning for artificial agents remains largely unsolved. While there has been increasing interest in data-driven approaches for the study of task planning for artificial agents, a significant remaining bottleneck is the dearth of large-scale comprehensive task-based datasets. In this paper, we present ActioNet, an interactive end-to-end platform for data collection and augmentation of task-based dataset in 3D environment. Using ActioNet, we collected a large-scale comprehensive task-based dataset, comprising over 3000 hierarchical task structures and videos. Using the hierarchical task structures, the videos are further augmented across 50 different scenes to give over 150,000 video. To our knowledge, ActioNet is the first interactive end-to-end platform for such task-based dataset generation and the accompanying dataset is the largest task-based dataset of such comprehensive nature. The ActioNet platform and dataset will be made available to facilitate research in hierarchical task planning.


title: Solving Quasiparticle Band Spectra of Real Solids using Neural-Network Quantum States

abstract: Establishing a predictive ab initio method for solid systems is one of the fundamental goals in condensed matter physics and computational materials science. The central challenge is how to encode a highly-complex quantum-many-body wave function compactly. Here, we demonstrate that artificial neural networks, known for their overwhelming expressibility in the context of machine learning, are excellent tool for first-principles calculations of extended periodic materials. We show that the ground-state energies in real solids in one-, two-, and three-dimensional systems are simulated precisely, reaching their chemical accuracy. The highlight of our work is that the quasiparticle band spectra, which are both essential and peculiar to solid-state systems, can be efficiently extracted with a computational technique designed to exploit the low-lying energy structure from neural networks. This work opens up a path to elucidate the intriguing and complex many-body phenomena in solid-state systems.


title: Perplexity-free Parametric t-SNE

abstract: The t-distributed Stochastic Neighbor Embedding (t-SNE) algorithm is a ubiquitously employed dimensionality reduction (DR) method. Its non-parametric nature and impressive efficacy motivated its parametric extension. It is however bounded to a user-defined perplexity parameter, restricting its DR quality compared to recently developed multi-scale perplexity-free approaches. This paper hence proposes a multi-scale parametric t-SNE scheme, relieved from the perplexity tuning and with a deep neural network implementing the mapping. It produces reliable embeddings with out-of-sample extensions, competitive with the best perplexity adjustments in terms of neighborhood preservation on multiple data sets.


title: Practical Precoding via Asynchronous Stochastic Successive Convex Approximation

abstract: We consider stochastic optimization of a smooth non-convex loss function with a convex non-smooth regularizer. In the online setting, where a single sample of the stochastic gradient of the loss is available at every iteration, the problem can be solved using the proximal stochastic gradient descent (SGD) algorithm and its variants. However in many problems, especially those arising in communications and signal processing, information beyond the stochastic gradient may be available thanks to the structure of the loss function. Such extra-gradient information is not used by SGD, but has been shown to be useful, for instance in the context of stochastic expectation-maximization, stochastic majorization-minimization, and stochastic successive convex approximation (SCA) approaches. By constructing a stochastic strongly convex surrogates of the loss function at every iteration, the stochastic SCA algorithms can exploit the structural properties of the loss function and achieve superior empirical performance as compared to the SGD. In this work, we take a closer look at the stochastic SCA algorithm and develop its asynchronous variant which can be used for resource allocation in wireless networks. While the stochastic SCA algorithm is known to converge asymptotically, its iteration complexity has not been well-studied, and is the focus of the current work. The insights obtained from the non-asymptotic analysis allow us to develop a more practical asynchronous variant of the stochastic SCA algorithm which allows the use of surrogates calculated in earlier iterations. We characterize precise bound on the maximum delay the algorithm can tolerate, while still achieving the same convergence rate. We apply the algorithm to the problem of linear precoding in wireless sensor networks, where it can be implemented at low complexity but is shown to perform well in practice.


title: Dynamics and Allocation of Transaction Cost in Multiagent Industrial Symbiosis

abstract: This paper discusses the dynamics of Transaction Cost (TC) in Industrial Symbiosis Institutions (ISI) and provides a fair and stable mechanism for TC allocation among the involved firms in a given ISI. In principle, industrial symbiosis, as an implementation of the circular economy paradigm in the context of industrial relation, is a practice aiming at reducing the material/energy footprint of the firm. The well-engineered form of this practice is proved to decrease the transaction costs at a collective level. This can be achieved using information systems for: identifying potential synergies, evaluating mutually beneficial ones, implementing the contracts, and governing the behavior of the established relations. Then the question is "how to distribute the costs for maintaining such an information system in a fair and stable manner?" We see such a cost as a collective transaction cost and employ an integrated method rooted in cooperative game theory and multiagent systems research to develop a fair and stable allocation mechanism for it. The novelty is twofold: in developing analytical multiagent methods for capturing the dynamics of transaction costs in industrial symbiosis and in presenting a novel game-theoretic mechanism for its allocation in industrial symbiosis institutions. While the former contributes to the theories of industrial symbiosis (methodological contribution), the latter supports decision makers aiming to specify fair and stable industrial symbiosis contracts (practical contribution).


title: COVID-19 Classification of X-ray Images Using Deep Neural Networks

abstract: In the midst of the coronavirus disease 2019 (COVID-19) outbreak, chest X-ray (CXR) imaging is playing an important role in the diagnosis and monitoring of patients with COVID-19. Machine learning solutions have been shown to be useful for X-ray analysis and classification in a range of medical contexts. The purpose of this study is to create and evaluate a machine learning model for diagnosis of COVID-19, and to provide a tool for searching for similar patients according to their X-ray scans. In this retrospective study, a classifier was built using a pre-trained deep learning model (ReNet50) and enhanced by data augmentation and lung segmentation to detect COVID-19 in frontal CXR images collected between January 2018 and July 2020 in four hospitals in Israel. A nearest-neighbors algorithm was implemented based on the network results that identifies the images most similar to a given image. The model was evaluated using accuracy, sensitivity, area under the curve (AUC) of receiver operating characteristic (ROC) curve and of the precision-recall (P-R) curve. The dataset sourced for this study includes 2362 CXRs, balanced for positive and negative COVID-19, from 1384 patients (63 +/- 18 years, 552 men). Our model achieved 89.7% (314/350) accuracy and 87.1% (156/179) sensitivity in classification of COVID-19 on a test dataset comprising 15% (350 of 2326) of the original data, with AUC of ROC 0.95 and AUC of the P-R curve 0.94. For each image we retrieve images with the most similar DNN-based image embeddings; these can be used to compare with previous cases.


title: Very ample linear series on quadrigonal curves

abstract: In this paper we describe all complete and very ample linear series on 4-gonal smooth curves.


title: Determination of the Local Hubble Constant from Virgo Infall Using TRGB Distances

abstract: An independent determination of $H_0$ is crucial given the growing tension between the Hubble constant, $H_0$, derived locally and that determined from the modeling of the cosmic microwave background (CMB) originating in the early universe. In this work, we present a new determination of $H_0$ using velocities and tip of the red giant branch (TRGB) distances to 33 galaxies located between the Local Group and the Virgo cluster. We use a model of the infall pattern of the local Hubble flow modified by the Virgo mass, which is given as a function of the cosmological constants ($H_0$, $\Omega_\Lambda$), the radius of the zero-velocity surface $R_0$, and the intrinsic velocity dispersion, $\sigma_v$. Fitting velocities and TRGB distances of 33 galaxies to the model, we obtain $H_0 =$ 65.8 $\pm$ 3.5($stat$) $\pm$ 2.4($sys$) km s$^{-1}$ Mpc$^{-1}$ and $R_0 =$ 6.76 $\pm$ 0.35 Mpc. Our local $H_0$ is consistent with the global $H_0$ determined from CMB radiation, showing no tension. In addition, we present new TRGB distances to NGC 4437 and NGC 4592 which are located near the zero-velocity surface: D = 9.28 $\pm$ 0.39 Mpc and D = 9.07 $\pm$ 0.27 Mpc, respectively. Their spatial separation is 0.29$^{+0.30}_{-0.03}$ Mpc, suggesting that they form a physical pair.


title: On the monophonic rank of a graph

abstract: A set of vertices $S$ of a graph $G$ is $monophonically \ convex$ if every induced path joining two vertices of $S$ is contained in $S$. The $monophonic \ convex \ hull$ of $S$, $\langle S \rangle$, is the smallest monophonically convex set containing $S$. A set $S$ is $monophonic \ convexly \ independent$ if $v \not\in \langle S - \{v\} \rangle$ for every $v \in S$. The $monophonic \ rank$ of $G$ is the size of the largest monophonic convexly independent set of $G$. We present a characterization of the monophonic convexly independent sets. Using this result, we show how to determine the monophonic rank of graph classes like bipartite, cactus, triangle-free and line graphs in polynomial time. Furthermore, we show that this parameter can be computed in polynomial time for $1$-starlike graphs, $i.e.$, for split graphs, and that its determination is $NP$-complete for $k$-starlike graphs for any fixed $k \ge 2$, a subclass of chordal graphs. We also consider this problem on the graphs whose intersection graph of the maximal prime subgraphs is a tree.


title: Properties of Rotational Symmetric multiple valued functions and their Reed-Muller-Fourier spectra

abstract: The concept of rotation symmetric functions from the Boolean domain is extended to the multiple-valued (MV) domain. It is shown that symmetric functions are a subset of the rotation symmetric functions. Functions exhibiting these kinds of symmetry may be given a compact value vector representation. It is shown that the Reed-Muller-Fourier spectrum of a function preserves the kind of symmetry and therefore it may be given a compact vector representation of the same length as the compact value vector of the corresponding function. A method is presented for calculating the RMF spectrum of symmetric and rotation symmetric functions from their compact representations. Examples are given for 3-valued and 4-valued functions.


title: EECBS: A Bounded-Suboptimal Search for Multi-Agent Path Finding

abstract: Multi-Agent Path Finding (MAPF), i.e., finding collision-free paths for multiple robots, is important for many applications where small runtimes are necessary, including the kind of automated warehouses operated by Amazon. CBS is a leading two-level search algorithm for solving MAPF optimally. ECBS is a bounded-suboptimal variant of CBS that uses focal search to speed up CBS by sacrificing optimality and instead guaranteeing that the costs of its solutions are within a given factor of optimal. In this paper, we study how to decrease its runtime even further using inadmissible heuristics. Motivated by Explicit Estimation Search (EES), we propose Explicit Estimation CBS (EECBS), a new bounded-suboptimal variant of CBS, that uses online learning to obtain inadmissible estimates of the cost of the solution of each high-level node and uses EES to choose which high-level node to expand next. We also investigate recent improvements of CBS and adapt them to EECBS. We find that EECBS with the improvements runs significantly faster than the state-of-the-art bounded-suboptimal MAPF algorithms ECBS, BCP-7, and eMDD-SAT on a variety of MAPF instances. We hope that the scalability of EECBS enables additional applications for bounded-suboptimal MAPF algorithms.


title: Modified gravity or imperfect dark matter: a model-independent discrimination

abstract: We analyze how to parametrize general modifications of the dark matter perturbations equations in a model-independent way. We prove that a general model with an imperfect and non-conserved dark matter fluid with bulk and shear viscosities and heat flux in a modified gravity scenario can be described with five general functions of time and scale. We focus on the sub-Hubble regime within the quasi-static approximation and calculate the observable power spectra of the galaxy distribution, galaxy velocities and weak lensing and find that these observables are only sensitive to three combinations of the initial five functions. Deviations of these three observable functions with respect to $\mathrm{\Lambda CDM}$ give us different characteristic signals which allow us to determine in which cases it is possible to discriminate a modification of gravity from an imperfect or non-conserved dark matter. Finally, we perform a Fisher forecast analysis for these three parameters and show an example for a particular model with shear viscosity.


title: Computational Separation Between Convolutional and Fully-Connected Networks

abstract: Convolutional neural networks (CNN) exhibit unmatched performance in a multitude of computer vision tasks. However, the advantage of using convolutional networks over fully-connected networks is not understood from a theoretical perspective. In this work, we show how convolutional networks can leverage locality in the data, and thus achieve a computational advantage over fully-connected networks. Specifically, we show a class of problems that can be efficiently solved using convolutional networks trained with gradient-descent, but at the same time is hard to learn using a polynomial-size fully-connected network.


title: Lyapunov-guided Deep Reinforcement Learning for Stable Online Computation Offloading in Mobile-Edge Computing Networks

abstract: Opportunistic computation offloading is an effective method to improve the computation performance of mobile-edge computing (MEC) networks under dynamic edge environment. In this paper, we consider a multi-user MEC network with time-varying wireless channels and stochastic user task data arrivals in sequential time frames. In particular, we aim to design an online computation offloading algorithm to maximize the network data processing capability subject to the long-term data queue stability and average power constraints. The online algorithm is practical in the sense that the decisions for each time frame are made without the assumption of knowing future channel conditions and data arrivals. We formulate the problem as a multi-stage stochastic mixed integer non-linear programming (MINLP) problem that jointly determines the binary offloading (each user computes the task either locally or at the edge server) and system resource allocation decisions in sequential time frames. To address the coupling in the decisions of different time frames, we propose a novel framework, named LyDROO, that combines the advantages of Lyapunov optimization and deep reinforcement learning (DRL). Specifically, LyDROO first applies Lyapunov optimization to decouple the multi-stage stochastic MINLP into deterministic per-frame MINLP subproblems. By doing so, it guarantees to satisfy all the long-term constraints by solving the per-frame subproblems that are much smaller in size. Then, LyDROO integrates model-based optimization and model-free DRL to solve the per-frame MINLP problems with low computational complexity. Simulation results show that under various network setups, the proposed LyDROO achieves optimal computation performance while stabilizing all queues in the system. Besides, it induces very low execution latency that is particularly suitable for real-time implementation in fast fading environments.


title: A Preconditioned Algorithm for Model-Based Iterative CT Reconstruction and Material Decomposition from Spectral CT Data

abstract: Model-based material decomposition is a statisticaliterative reconstruction framework where basis material densityimages are estimated directly from spectral CT data. This methoduses a physical model for polyenergetic x-ray transmission andattenuation and therefore it does not typically suffer frombeam-hardening artifacts. However, this estimation is a poorly-conditioned inverse problem due to the strong anticorrelationbetween basis materials. In this work we propose an precondi-tioned optimization algorithm for a nonlinear penalized weightedleast-squares objective function.


title: Vector Potential and Surface Magnetic Field in Magnetoelectric Antiferromagnetic Materials

abstract: A general formula for the average vector potential of bulk periodic systems is proposed and shown to set the boundary conditions at magnetic interfaces. For antiferromagnetic materials, the study reveals a unique relation between the macroscopic potential and the orientation-dependent magnetic quadrupole, as a result of the different crystalline and magnetic symmetries. In particular, at surfaces and interfaces of a truncated bulk without inversion and time-reversal symmetries, the average vector potential exhibits a discontinuity, which results in an interfacial magnetic field. In general, however, due to the surface and interface electronic and atomic relaxations, additional magnetization may result. For the experimentally-observed magnetoelectric antiferromagnets, in particular, our symmetry analysis suggest that the relaxation effects could well be a system response to the presence of such a potential discontinuity.


title: Ethna: Analyzing the Underlying Peer-to-Peer Network of the Ethereum Blockchain

abstract: The peer-to-peer (P2P) network of blockchain used to transport its transactions and blocks has a high impact on the efficiency and security of the system. The P2P network topologies of popular blockchains such as Bitcoin and Ethereum, therefore, deserve our highest attention. The current Ethereum blockchain explorers (e.g., Etherscan) focus on the tracking of block and transaction records but omit the characterization of the underlying P2P network. This work presents the Ethereum Network Analyzer (Ethna), a tool that probes and analyzes the P2P network of the Ethereum blockchain. Unlike Bitcoin that adopts an unstructured P2P network, Ethereum relies on the Kademlia DHT to manage its P2P network. Therefore, the existing analytical methods for Bitcoin-like P2P networks are not applicable to Ethereum. Ethna implements a novel method that accurately measures the degrees of Ethereum nodes. Furthermore, it incorporates an algorithm that derives the latency metrics of message propagation in the Ethereum P2P network. We ran Ethna on the Ethereum Mainnet and conducted extensive experiments to analyze the topological features of its P2P network. Our analysis shows that the Ethereum P2P network possesses a certain effect of small-world networks, and the degrees of nodes follow a power-law distribution that characterizes scale-free networks.


title: Exponential Lower Bounds for Planning in MDPs With Linearly-Realizable Optimal Action-Value Functions

abstract: We consider the problem of local planning in fixed-horizon and discounted Markov Decision Processes (MDPs) with linear function approximation and a generative model under the assumption that the optimal action-value function lies in the span of a feature map that is available to the planner. Previous work has left open the question of whether there exist sound planners that need only poly(H,d) queries regardless of the MDP, where H is the horizon and d is the dimensionality of the features. We answer this question in the negative: we show that any sound planner must query at least $\min(\exp({\Omega}(d)), {\Omega}(2^H))$ samples in the fized-horizon setting and $\exp({\Omega}(d))$ samples in the discounted setting. We also show that for any ${\delta}>0$, the least-squares value iteration algorithm with $O(H^5d^{H+1}/{\delta}^2)$ queries can compute a ${\delta}$-optimal policy in the fixed-horizon setting. We discuss implications and remaining open questions.


title: A new result on the existence and non-existence of positive solutions for two-parametric systems of quasilinear elliptic equations

abstract: This paper is devoted to the existence and non-existence of positive solutions for a $(p, q)$-Laplacian system with indefinite nonlinearity depending on two parameters $(\lambda,\mu)$. By using the sub-supersolution method together with adaptations of ideas found in \cite{AAS,CRMZ}, we extend some previous result.


title: Sparse Quantized Spectral Clustering

abstract: Given a large data matrix, sparsifying, quantizing, and/or performing other entry-wise nonlinear operations can have numerous benefits, ranging from speeding up iterative algorithms for core numerical linear algebra problems to providing nonlinear filters to design state-of-the-art neural network models. Here, we exploit tools from random matrix theory to make precise statements about how the eigenspectrum of a matrix changes under such nonlinear transformations. In particular, we show that very little change occurs in the informative eigenstructure even under drastic sparsification/quantization, and consequently that very little downstream performance loss occurs with very aggressively sparsified or quantized spectral clustering. We illustrate how these results depend on the nonlinearity, we characterize a phase transition beyond which spectral clustering becomes possible, and we show when such nonlinear transformations can introduce spurious non-informative eigenvectors.


title: A Sum-Product Estimate for Well Spaced Sets

abstract: We study the $\delta$-discretized sum-product estimates for well spaced sets. Our main result is: for a fixed $\alpha\in(1,\frac{3}{2}]$, we prove that for any $\sim|A|^{-1}$-separated set $A\subset[1,2]$ and $\delta=|A|^{-\alpha}$, we have: $\mathcal{N}(A+A, \delta)\cdot \mathcal{N}(AA, \delta) \gtrsim_{\epsilon}|A|\delta^{-1+\epsilon}$.


title: MFPC-Net: Multi-fidelity Physics-Constrained Neural Process

abstract: In this work, we propose a network which can utilize computational cheap low-fidelity data together with limited high-fidelity data to train surrogate models, where the multi-fidelity data are generated from multiple underlying models. The network takes a context set as input (physical observation points, low fidelity solution at observed points) and output (high fidelity solution at observed points) pairs. It uses the neural process to learn a distribution over functions conditioned on context sets and provide the mean and standard deviation at target sets. Moreover, the proposed framework also takes into account the available physical laws that govern the data and imposes them as constraints in the loss function. The multi-fidelity physical constraint network (MFPC-Net) (1) takes datasets obtained from multiple models at the same time in the training, (2) takes advantage of available physical information, (3) learns a stochastic process which can encode prior beliefs about the correlation between two fidelity with a few observations, and (4) produces predictions with uncertainty. The ability of representing a class of functions is ensured by the property of neural process and is achieved by the global latent variables in the neural network. Physical constraints are added to the loss using Lagrange multipliers. An algorithm to optimize the loss function is proposed to effectively train the parameters in the network on an ad hoc basis. Once trained, one can obtain fast evaluations at the entire domain of interest given a few observation points from a new low-and high-fidelity model pair. Particularly, one can further identify the unknown parameters such as permeability fields in elliptic PDEs with a simple modification of the network. Several numerical examples for both forward and inverse problems are presented to demonstrate the performance of the proposed method.


title: Symmetry breaking patterns, tricriticalities and quadruple points in quantum Rabi model with bias and nonlinear interaction

abstract: Quantum Rabi model (QRM) is fascinating not only because of its broad relevance and but also due to its few-body quantum phase transition. In practice both the bias and the nonlinear coupling in QRM are important controlling parameters in experimental setups. We study the interplay of the bias and the nonlinear interaction with the linear coupling in the ground state which exhibits various patterns of symmetry breaking and different orders of transitions. Several situations of tricriticality are unveiled in the low frequency limit and at finite frequencies. We find that the full quantum-mechanical effect leads to novel transitions, tricriticalities and quadruple points, which are much beyond the semiclassical picture. We clarify the underlying mechanisms by analyzing the energy competitions and the essential changeovers of the quantum states, which enables us to extract most analytic phase boundaries.


title: Predicting traffic overflows on private peering

abstract: Large content providers and content distribution network operators usually connect with large Internet service providers (eyeball networks) through dedicated private peering. The capacity of these private network interconnects is provisioned to match the volume of the real content demand by the users. Unfortunately, in case of a surge in traffic demand, for example due to a content trending in a certain country, the capacity of the private interconnect may deplete and the content provider/distributor would have to reroute the excess traffic through transit providers. Although, such overflow events are rare, they have significant negative impacts on content providers, Internet service providers, and end-users. These include unexpected delays and disruptions reducing the user experience quality, as well as direct costs paid by the Internet service provider to the transit providers. If the traffic overflow events could be predicted, the Internet service providers would be able to influence the routes chosen for the excess traffic to reduce the costs and increase user experience quality. In this article we propose a method based on an ensemble of deep learning models to predict overflow events over a short term horizon of 2-6 hours and predict the specific interconnections that will ingress the overflow traffic. The method was evaluated with 2.5 years' traffic measurement data from a large European Internet service provider resulting in a true-positive rate of 0.8 while maintaining a 0.05 false-positive rate. The lockdown imposed by the COVID-19 pandemic reduced the overflow prediction accuracy. Nevertheless, starting from the end of April 2020 with the gradual lockdown release, the old models trained before the pandemic perform equally well.


title: Cubic Spline Smoothing Compensation for Irregularly Sampled Sequences

abstract: The marriage of recurrent neural networks and neural ordinary differential networks (ODE-RNN) is effective in modeling irregularly-observed sequences. While ODE produces the smooth hidden states between observation intervals, the RNN will trigger a hidden state jump when a new observation arrives, thus cause the interpolation discontinuity problem. To address this issue, we propose the cubic spline smoothing compensation, which is a stand-alone module upon either the output or the hidden state of ODE-RNN and can be trained end-to-end. We derive its analytical solution and provide its theoretical interpolation error bound. Extensive experiments indicate its merits over both ODE-RNN and cubic spline interpolation.


title: A Taxonomy of Polytomous Item Response Models

abstract: A common framework is provided that comprises classical ordinal item response models as the cumulative, sequential and adjacent categories models as well as nominal response models and item response tree models. The taxonomy is based on the ways binary models can be seen as building blocks of the various models. In particular one can distinguish between conditional and unconditional model components. Conditional models are by far the larger class of models containing the adjacent categories model and the whole class of hierarchically structured models. The latter is introduced as a class of models that comprises binary trees and hierarchically structured models that use ordinal models conditionally. The study of the binary models contained in latent trait models clarifies the relation between models and the interpretation of item parameters. It is also used to distinguish between ordinal and nominal models by giving a conceptualization of ordinal models. The taxonomy differs from previous taxonomies by focusing on the structured use of dichotomizations instead of the role of parameterizations.


title: Spectral Fractional Laplacian with Inhomogeneous Dirichlet Data: Questions, Problems, Solutions

abstract: In this paper we discuss the topic of correct setting for the equation $(-\Delta )^s u=f$, with $0<s <1$. The definition of the fractional Laplacian on the whole space $\mathbb R^n$, $n=1,2,3$ is understood through the Fourier transform, see, e.g., Lischke et.al. (J. Comp. Phys., 2020). The real challenge however represents the case when this equation is posed in a bounded domain $\Omega$ and proper boundary conditions are needed for the correctness of the corresponding problem. Let us mention here that the case of inhomogeneous boundary data has been neglected up to the last years. The reason is that imposing nonzero boundary conditions in the nonlocal setting is highly nontrivial. There exist at least two different definitions of fractional Laplacian, and there is still ongoing research about the relations of them. They are not equivalent. The focus of our study is a new characterization of the spectral fractional Laplacian. One of the major contributions concerns the case when the right hand side $f$ is a Dirac $\delta$ function. For comparing the differences between the solutions in the spectral and Riesz formulations, we consider an inhomogeneous fractional Dirichlet problem. The provided theoretical analysis is supported by model numerical tests.


title: An adaptive Hessian approximated stochastic gradient MCMC method

abstract: Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their scalability to handle large datasets and the ability to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be inefficient when the random variables in the target posterior densities have scale differences or are highly correlated. In this work, we present an adaptive Hessian approximated stochastic gradient MCMC method to incorporate local geometric information while sampling from the posterior. The idea is to apply stochastic approximation to sequentially update a preconditioning matrix at each iteration. The preconditioner possesses second-order information and can guide the random walk of a sampler efficiently. Instead of computing and saving the full Hessian of the log posterior, we use limited memory of the sample and their stochastic gradients to approximate the inverse Hessian-vector multiplication in the updating formula. Moreover, by smoothly optimizing the preconditioning matrix, our proposed algorithm can asymptotically converge to the target distribution with a controllable bias under mild conditions. To reduce the training and testing computational burden, we adopt a magnitude-based weight pruning method to enforce the sparsity of the network. Our method is user-friendly and is scalable to standard SG-MCMC updating rules by implementing an additional preconditioner. The sparse approximation of inverse Hessian alleviates storage and computational complexities for large dimensional models. The bias introduced by stochastic approximation is controllable and can be analyzed theoretically. Numerical experiments are performed on several problems.


title: Limitations of Sums of Bounded-Read Formulas

abstract: Proving super polynomial size lower bounds for various classes of arithmetic circuits computing explicit polynomials is a very important and challenging task in algebraic complexity theory. We study representation of polynomials as sums of weaker models such as read once formulas (ROFs) and read once oblivious algebraic branching programs (ROABPs). We prove: (1) An exponential separation between sum of ROFs and read-$k$ formulas for some constant $k$. (2) A sub-exponential separation between sum of ROABPs and syntactic multilinear ABPs. Our results are based on analysis of the partial derivative matrix under different distributions. These results highlight richness of bounded read restrictions in arithmetic formulas and ABPs. Finally, we consider a generalization of multilinear ROABPs known as strict-interval ABPs defined in [Ramya-Rao, MFCS2019]. We show that strict-interval ABPs are equivalent to ROABPs upto a polynomial size blow up. In contrast, we show that interval formulas are different from ROFs and also admit depth reduction which is not known in the case of strict-interval ABPs.


title: Flood Evacuation During Pandemic: A multi-objective Framework to Handle Compound Hazard

abstract: The evacuation of the population from flood-affected regions is a non-structural measure to mitigate flood hazards. Shelters used for this purpose usually accommodate a large number of flood evacuees for a temporary period. Floods during pandemic result in a compound hazard. Evacuations under such situations are difficult to plan as social distancing is nearly impossible in the highly crowded shelters. This results in a multi-objective problem with conflicting objectives of maximizing the number of evacuees from flood-prone regions and minimizing the number of infections at the end of the shelter's stay. To the best of our knowledge, such a problem is yet to be explored in literature. Here we develop a simulation-optimization framework, where multiple objectives are handled with a max-min approach. The simulation model consists of an extended Susceptible Exposed Infectious Recovered Susceptible (SEIRS) model.We apply the proposed model to the flood-prone Jagatsinghpur district in the state of Odisha, India. We find that the proposed approach can provide an estimate of people required to be evacuated from individual flood-prone villages to reduce flood hazards during the pandemic. At the same time, this does not result in an uncontrolled number of new infections. The proposed approach can generalize to different regions and can provide a framework to stakeholders to manage conflicting objectives in disaster management planning and to handle compound hazards.


title: DuoBFT: Resilience vs. Performance Trade-off in Byzantine Fault Tolerance

abstract: This paper presents DuoBFT, a Byzantine fault-tolerant protocol that uses trusted components to provide commit decisions in the Hybrid fault model in addition to commit decisions in the BFT model. By doing so, it enables the clients to choose the response fault model for its commands. Internally, DuoBFT commits each client command under both the hybrid and Byzantine models, but since hybrid commits take fewer communication steps and use smaller quorums than BFT commits, clients can benefit from the low-latency commits in the hybrid model. DuoBFT uses a common view-change change protocol to handle both fault models. To achieve this, we enable a notion called Flexible Quorums in the hybrid fault model by revisiting the quorum intersection requirements in hybrid protocols. The flexible quorum technique enables having a hybrid view change quorum that is of the same size as a BFT view-change quorum. This paves a path for efficiently combining both the fault models within a single unified protocol. Our evaluation on a wide-area deployment reveal that DuoBFT can provide hybrid commits with 30% lower latency to existing protocols without sacrificing throughput. In absolute terms, DuoBFT provides sub-200-millisecond latency in a geographically replicated deployment.


title: Online Neural Networks for Change-Point Detection

abstract: Moments when a time series changes its behaviour are called change points. Detection of such points is a well-known problem, which can be found in many applications: quality monitoring of industrial processes, failure detection in complex systems, health monitoring, speech recognition and video analysis. Occurrence of change point implies that the state of the system is altered and its timely detection might help to prevent unwanted consequences. In this paper, we present two online change-point detection approaches based on neural networks. These algorithms demonstrate linear computational complexity and are suitable for change-point detection in large time series. We compare them with the best known algorithms on various synthetic and real world data sets. Experiments show that the proposed methods outperform known approaches.


title: Discrete shedding of secondary vortices along a modified Kaden spiral

abstract: When an object is accelerated in a fluid, a primary vortex is formed through the roll-up of a shear layer. This primary vortex does not grow indefinitely and will reach a limiting size and strength. Additional vorticity beyond the critical limit will end up in a trailing shear layer and accumulate into secondary vortices. The secondary vortices are typically considerably smaller than the primary vortex. Here, we focus on the formation, shedding, and trajectory of secondary vortices generated by a rotating rectangular plate in a quiescent fluid using time-resolved particle image velocimetry. The Reynolds number (Re) is varied from 840 to 11150. At low Re, the shear layer is a continuous uninterrupted layer of vorticity that rolls up into a single coherent primary vortex. At Re=1955, the shear layer becomes unstable. For Re>4000, secondary vortices are discretely released from the plate tip. First, we demonstrate that the roll-up of the shear layer, the trajectory of the primary vortex, and the path of secondary vortices can be predicted by a modified Kaden spiral. Second, the timing of the secondary vortex shedding is analysed. The time interval between the release of successive secondary vortices is not constant during the rotation but increases the more vortices have been shed. The shedding time interval also increases with decreasing Reynolds number. The increased time interval under both conditions is due to a reduced circulation feeding rate.


title: Broadcasting on Two-Dimensional Regular Grids

abstract: We study a specialization of the problem of broadcasting on directed acyclic graphs, namely, broadcasting on 2D regular grids. Consider a 2D regular grid with source vertex $X$ at layer $0$ and $k+1$ vertices at layer $k\geq 1$, which are at distance $k$ from $X$. Every vertex of the 2D regular grid has outdegree $2$, the vertices at the boundary have indegree $1$, and all other vertices have indegree $2$. At time $0$, $X$ is given a random bit. At time $k\geq 1$, each vertex in layer $k$ receives transmitted bits from its parents in layer $k-1$, where the bits pass through binary symmetric channels with noise level $\delta\in(0,1/2)$. Then, each vertex combines its received bits using a common Boolean processing function to produce an output bit. The objective is to recover $X$ with probability of error better than $1/2$ from all vertices at layer $k$ as $k \rightarrow \infty$. Besides their natural interpretation in communication networks, such broadcasting processes can be construed as 1D probabilistic cellular automata (PCA) with boundary conditions that limit the number of sites at each time $k$ to $k+1$. We conjecture that it is impossible to propagate information in a 2D regular grid regardless of the noise level and the choice of processing function. In this paper, we make progress towards establishing this conjecture, and prove using ideas from percolation and coding theory that recovery of $X$ is impossible for any $\delta$ provided that all vertices use either AND or XOR processing functions. Furthermore, we propose a martingale-based approach that establishes the impossibility of recovering $X$ for any $\delta$ when all NAND processing functions are used if certain supermartingales can be rigorously constructed. We also provide numerical evidence for the existence of these supermartingales by computing explicit examples for different values of $\delta$ via linear programming.


title: An Infinite, Converging, Sequence of Brocard Porisms

abstract: The Brocard porism is a known 1d family of triangles inscribed in a circle and circumscribed about an ellipse. Remarkably, the Brocard angle is invariant and the Brocard points are stationary at the foci of the ellipse. In this paper we show that a certain derived triangle spawns off a second, smaller, Brocard porism so that repeating this calculation produces an infinite, converging sequence of porisms. We also show that this sequence is embedded in a continuous family of porisms.


title: CardioXNet: A Novel Lightweight Deep Learning Framework for Cardiovascular Disease Classification Using Heart Sound Recordings

abstract: The alarmingly high mortality rate and increasing global prevalence of cardiovascular diseases signify the crucial need for early detection schemes. Phonocardiogram (PCG) signals have been historically applied in this domain owing to its simplicity and cost-effectiveness. In this paper, we propose CardioXNet, a novel lightweight end-to-end CRNN architecture for automatic detection of five classes of cardiac auscultation namely normal, aortic stenosis, mitral stenosis, mitral regurgitation and mitral valve prolapse using raw PCG signal. The process has been automated by the involvement of two learning phases. Three parallel CNN pathways have been implemented in the representation learning phase to learn the coarse and fine-grained features from the PCG and to explore the salient features from variable receptive fields involving 2D-CNN based squeeze-expansion. Thus, in the representation learning phase, the network extracts efficient time-invariant features and converges with great rapidity. In the sequential residual learning phase, with the bidirectional-LSTMs and the skip connection, the network can proficiently extract temporal features without performing any feature extraction on the signal. The obtained results demonstrate that the proposed end-to-end architecture yields outstanding performance in all the evaluation metrics compared to the previous state-of-the-art methods with up to 99.60% accuracy, 99.56% precision, 99.52% recall and 99.68% F1- score on an average while being computationally comparable. This model outperforms the previous works using the same dataset by a considerable margin. The high accuracy metrics on both primary and secondary dataset combined with a significantly low number of parameters and end-to-end prediction approach makes the proposed network suitable for point of care CVD screening in low resource setups using memory constraint mobile devices.


title: Consistent solution of Einstein-Cartan equations with torsion outside matter

abstract: The Einstein-Cartan equations in first-order action of torsion are considered. From Belinfante-Rosenfeld equation special consistence conditions are derived for the torsion parameters relating them to the metric. Inside matter the torsion is given by the spin which leads to an extended Oppenhaimer-Volkov equation. Outside matter a second solution is found besides the torsion-free Schwarzschild one with the torsion completely determined by the metric and vice-versa. This solution is shown to be of non-spherical origin and its uniqueness with respect to the consistence is demonstrated. Unusual properties are discussed in different coordinate systems where the cosmological constant assumes the role of the Friedman parameter in Friedman-Lama\^itre-Robertson-Walker cosmoses. Parameters are specified where wormholes are possible. Transformations are presented to explore and map regions of expanding and contracting universes to the form of static metrics. The autoparallel equations are solved exactly and compared with geodesic motion. The Weyl tensor reveals that the here found solution is of Petrov-D type.


title: A postprocessing technique for a discontinuous Galerkin discretization of time-dependent Maxwell's equations

abstract: We present a novel postprocessing technique for a discontinuous Galerkin (DG) discretization of time-dependent Maxwell's equations that we couple with an explicit Runge-Kutta time-marching scheme. The postprocessed electromagnetic field converges one order faster than the unprocessed solution in the H(curl)-norm. The proposed approach is local, in the sense that the enhanced solution is computed independently in each cell of the computational mesh, and at each time step of interest. As a result, it is inexpensive to compute, especially if the region of interest is localized, either in time or space. The key ideas behind this postprocessing technique stem from hybridizable discontinuous Galerkin (HDG) methods, which are equivalent to the analyzed DG scheme for specific choices of penalization parameters. We present several numerical experiments that highlight the superconvergence properties of the postprocessed electromagnetic field approximation.


title: Evidence for sulfur-bearing species on Callisto's leading hemisphere: Sourced from Jupiter's irregular satellites or Io?

abstract: We investigated whether sulfur-bearing species are present on the icy Galilean moon Callisto by analyzing eight near-infrared reflectance spectra collected over a wide range of sub-observer longitudes. We measured the band areas and depths of a 4-micron feature in these spectra, which has been attributed to sulfur dioxide (SO2), as well as carbonates, in previously collected datasets of this moon. All eight spectra we collected display the 4-micron band. The four spectra collected over Callisto's leading hemisphere display significantly stronger 4-micron bands compared to the four trailing hemisphere spectra (> 3-sigma difference). We compared the central wavelength position and shape of Callisto's 4-micron band to laboratory spectra of various sulfur-bearing species and carbonates. Our comparison demonstrates that Callisto's 4-micron band has a spectral signature similar to thermally-altered sulfur, as well as a 4.025 micron feature attributed to disulfanide (HS2). Our analysis therefore supports the presence of S-bearing species on Callisto but is not consistent with the presence of SO2. The significantly stronger 4-micron band detected on Callisto's leading hemisphere could result from collisions with H2S-rich dust grains that originate on Jupiter's retrograde irregular satellites or implantation of magnetospheric S ions that originate from volcanic activity on Io. Alternatively, S-bearing species could be native to Callisto and are exposed by dust collisions and larger impacts that drive regolith overturn, primarily on its leading side.


title: Regularized Bayesian calibration and scoring of the WD-FAB IRT model improves predictive performance over marginal maximum likelihood

abstract: Item response theory (IRT) is the statistical paradigm underlying a dominant family of generative probabilistic models for test responses, used to quantify traits in individuals relative to target populations. The graded response model (GRM) is a particular IRT model that is used for ordered polytomous test responses. Both the development and the application of the GRM and other IRT models require statistical decisions. For formulating these models (calibration), one needs to decide on methodologies for item selection, inference, and regularization. For applying these models (test scoring), one needs to make similar decisions, often prioritizing computational tractability and/or interpretability. In many applications, such as in the Work Disability Functional Assessment Battery (WD-FAB), tractability implies approximating an individual's score distribution using estimates of mean and variance, and obtaining that score conditional on only point estimates of the calibrated model. In this manuscript, we evaluate the calibration and scoring of models under this common use-case using Bayesian cross-validation. Applied to the WD-FAB responses collected for the National Institutes of Health, we assess the predictive power of implementations of the GRM based on their ability to yield, on validation sets of respondents, ability estimates that are most predictive of patterns of item responses. Our main finding indicates that regularized Bayesian calibration of the GRM outperforms the regularization-free empirical Bayesian procedure of marginal maximum likelihood. We also motivate the use of compactly supported priors in test scoring.