title: On the Power of False Negative Awareness in Indicator-based Caching Systems

abstract: Distributed caching systems such as content distribution networks often advertise their content via lightweight approximate indicators (e.g., Bloom filters) to efficiently inform clients where each datum is likely cached. While false-positive indications are necessary and well understood, most existing works assume no false-negative indications. Our work illustrates practical scenarios where false-negatives are unavoidable and ignoring them has a significant impact on system performance. Specifically, we focus on false-negatives induced by indicator staleness, which arises whenever the system advertises the indicator only periodically, rather than immediately reporting every change in the cache. Such scenarios naturally occur, e.g., in bandwidth-constraint environments or when latency impedes the ability of each client to obtain an updated indicator. Our work introduces novel false-negative aware access policies that continuously estimate the false-negative ratio and sometimes access caches despite negative indications. We present optimal policies for homogeneous settings and provide approximation guarantees for our algorithms in heterogeneous environments. We further perform an extensive simulation study with multiple real system traces. We show that our false-negative aware algorithms incur a significantly lower access cost than existing approaches or match the cost of these approaches while requiring an order of magnitude fewer resources (e.g., caching capacity or bandwidth).


title: A survey of optimal polynomial approximants, applications to digital filter design, and related open problems

abstract: In the last few years, the notion of optimal polynomial approximant has appeared in the mathematics literature in connection with Hilbert spaces of analytic functions of one or more variables. In the 70s, researchers in engineering and applied mathematics introduced least-squares inverses in the context of digital filters in signal processing. It turns out that in the Hardy space $H^2$ these objects are identical. This paper is a survey of known results about optimal polynomial approximants. In particular, we will examine their connections with orthogonal polynomials and reproducing kernels in weighted spaces and digital filter design. We will also describe what is known about the zeros of optimal polynomial approximants, their rates of decay, and convergence results. Throughout the paper, we state many open questions that may be of interest.


title: Invariant hyperbolic curves: determinantal representations and applications to the numerical range

abstract: Here we study the space of real hyperbolic plane curves that are invariant under actions of the cyclic and dihedral groups and show they have determinantal representations that certify this invariance. We show an analogue of Nuij's theorem for the set of invariant hyperbolic polynomials of a given degree. The main theorem is that every invariant hyperbolic plane curve has a determinantal representation using a block cyclic weighted shift matrix. This generalizes previous work by Lentzos and the first author, as well as by Chien and Nakazato. One consequence is that if the numerical range of a matrix is invariant under rotation, then it is the numerical range of a block cyclic weighted shift matrix.


title: Pecan: An Automated Theorem Prover for Automatic Sequences using B\"uchi Automata

abstract: Pecan is an automated theorem prover for reasoning about properties of Sturmian words, an important object in the field of combinatorics on words. It is capable of efficiently proving non-trivial mathematical theorems about all Sturmian words.


title: A Novel Transfer Learning-Based Approach for Screening Pre-existing Heart Diseases Using Synchronized ECG Signals and Heart Sounds

abstract: Diagnosing pre-existing heart diseases early in life is important as it helps prevent complications such as pulmonary hypertension, heart rhythm problems, blood clots, heart failure and sudden cardiac arrest. To identify such diseases, phonocardiogram (PCG) and electrocardiogram (ECG) waveforms convey important information. Therefore, effectively using these two modalities of data has the potential to improve the disease screening process. We evaluate this hypothesis on a subset of the PhysioNet Challenge 2016 Dataset which contains simultaneously acquired PCG and ECG recordings. Our novel Dual-Convolutional Neural Network based approach uses transfer learning to tackle the problem of having limited amounts of simultaneous PCG and ECG data that is publicly available, while having the potential to adapt to larger datasets. In addition, we introduce two main evaluation frameworks named record-wise and sample-wise evaluation which leads to a rich performance evaluation for the transfer learning approach. Comparisons with methods which used single or dual modality data show that our method can lead to better performance. Furthermore, our results show that individually collected ECG or PCG waveforms are able to provide transferable features which could effectively help to make use of a limited number of synchronized PCG and ECG waveforms and still achieve significant classification performance.


title: Majorizing Measures, Sequential Complexities, and Online Learning

abstract: We introduce the technique of generic chaining and majorizing measures for controlling sequential Rademacher complexity. We relate majorizing measures to the notion of fractional covering numbers, which we show to be dominated in terms of sequential scale-sensitive dimensions in a horizon-independent way, and, under additional complexity assumptions establish a tight control on worst-case sequential Rademacher complexity in terms of the integral of sequential scale-sensitive dimension. Finally, we establish a tight contraction inequality for worst-case sequential Rademacher complexity. The above constitutes the resolution of a number of outstanding open problems in extending the classical theory of empirical processes to the sequential case, and, in turn, establishes sharp results for online learning.


title: On Greedy Approaches to Hierarchical Aggregation

abstract: We analyze greedy algorithms for the Hierarchical Aggregation (HAG) problem, a strategy introduced in [Jia et al., KDD 2020] for speeding up learning on Graph Neural Networks (GNNs). The idea of HAG is to identify and remove redundancies in computations performed when training GNNs. The associated optimization problem is to identify and remove the most redundancies. Previous work introduced a greedy approach for the HAG problem and claimed a 1-1/e approximation factor. We show by example that this is not correct, and one cannot hope for better than a 1/2 approximation factor. We prove that this greedy algorithm does satisfy some (weaker) approximation guarantee, by showing a new connection between the HAG problem and maximum matching problems in hypergraphs. We also introduce a second greedy algorithm which can out-perform the first one, and we show how to implement it efficiently in some parameter regimes. Finally, we introduce some greedy heuristics that are much faster than the above greedy algorithms, and we demonstrate that they perform well on real-world graphs.


title: Study on the stability of thermoelastic Bresse and Timoshenko type systems with Gurtin-Pipkin's law via the vertical displacements

abstract: The objective of this paper is to study the stability of a linear one-dimensional thermoelastic Bresse system in a bounded domain, where the coupling is given through the first component of the Bresse model with the heat conduction of Gurtin-Pipkin type. Two kinds of coupling are considered; the first coupling is of order one with respect to space variable, and the second one is of order zero. We state the well-posedness and show the polynomial stability of the systems, where the decay rates depend on the smoothness of initial data. Moreover, in case of coupling of order one, we prove the equivalence between the exponential stability and some new conditions on the parameters of the system. However, when the coupling is of order zero, we prove the non-exponential stability independently of the parameters of the system. Applications to the corresponding particular Timoshenko models are also given, where we prove that both couplings lead to the exponential stability if and only if some conditions on the parameters of the systems are satisfied, and both couplings guarantee the polynomial stability independently of the parameters of the systems. The proof is based on the semigroup theory and a combination of the energy method and the frequency domain approach.


title: Truly Sparse Neural Networks at Scale

abstract: Recently, sparse training methods have started to be established as a de facto approach for training and inference efficiency in artificial neural networks. Yet, this efficiency is just in theory. In practice, everyone uses a binary mask to simulate sparsity since the typical deep learning software and hardware are optimized for dense matrix operations. In this paper, we take an orthogonal approach, and we show that we can train truly sparse neural networks to harvest their full potential. To achieve this goal, we introduce three novel contributions, specially designed for sparse neural networks: (1) a parallel training algorithm and its corresponding sparse implementation from scratch, (2) an activation function with non-trainable parameters to favour the gradient flow, and (3) a hidden neurons importance metric to eliminate redundancies. All in one, we are able to break the record and to train the largest neural network ever trained in terms of representational power -- reaching the bat brain size. The results show that our approach has state-of-the-art performance while opening the path for an environmentally friendly artificial intelligence era.


title: FedProf: Selective Federated Learning with Representation Profiling

abstract: Federated Learning (FL) has shown great potential as a privacy-preserving solution to learning from decentralized data that are only accessible to end devices (i.e., clients). In many scenarios, however, a large proportion of the clients are probably in possession of low-quality data that are biased, noisy or even irrelevant. As a result, they could significantly slow down the convergence of the global model we aim to build and also compromise its quality. In light of this, we propose FedProf, a novel algorithm for optimizing FL under such circumstances without breaching data privacy. The key of our approach is a distributional representation profiling and matching scheme that uses the global model to dynamically profile data representations and allows for low-cost, lightweight representation matching. Based on the scheme we adaptively score each client and adjust its participation probability so as to mitigate the impact of low-value clients on the training process. We have conducted extensive experiments on public datasets using various FL settings. The results show that the selective behaviour of our algorithm leads to a significant reduction in the number of communication rounds and the amount of time (up to 2.4x speedup) for the global model to converge and also provides accuracy gain.


title: A 3D radiative transfer framework: XII. Many-core, vector and GPU methods

abstract: 3D detailed radiative transfer is computationally taxing, since the solution of the radiative transfer equation involves traversing the six dimensional phase space of the 3D domain. With modern supercomputers the hardware available for wallclock speedup is rapidly changing, mostly in response to requirements to minimize the cost of electrical power. Given the variety of modern computing architectures, we aim to develop and adapt algorithms for different computing architectures to improve performance on a wide variety of platforms. We implemented the main time consuming kernels for solving 3D radiative transfer problems for vastly different computing architectures using MPI, OpenMP, OpenACC and vector algorithms. Adapted algorithms lead to massively improved speed for all architectures, making extremely large model calculations easily feasible. These calculations would have previously been considered impossible or prohibitively expensive. Efficient use of modern computing devices is entirely feasible, but unfortunately requires the implementation of specialized algorithms for them.


title: On the decay rates for two Cauchy thermoelastic laminated Timoshenko problems of type III with interfacial slip

abstract: The subject of this paper is to study the decay of solutions for two systems of laminated Timoshenko beams with interfacial slip in the whole space R subject to a thermal effect of type III acting only on one component. When the thermal effect is acting via the second or third component of the laminated Timoshenko beam (rotation angle displacement or dynamic of the slip), we prove that both systems are polynomially stable and obtain stability estimates in the L2 -norm of solutions and their higher order derivatives with respect of the space variable. The decay rates, as well as the absence and presence of the regularity-loss type property, depend on the regularity of the initial data and the speeds of wave propagations. However, when the thermal effect is acting via the first comoponent (transversal displacement), we introduce a new stability number \c{hi} and prove that the stability of systems is equivalent to \c{hi} 6= 0. An application to a case of lower order coupling terms will be also given. To prove our results, we use the energy method in Fourier space combined with well chosen weight functions to build appropriate Lyapunov functionals.


title: Open Quantum Systems for Quarkonia

abstract: I review recent applications of the open quantum system framework in the understanding of quarkonium suppression in heavy-ion collisions, which has been used as a probe of the quark-gluon plasma for decades. The derivation of the Lindblad equations for quarkonium in both the quantum Brownian motion and the quantum optical limits and their semiclassical counterparts is explained. The hierarchy of time scales assumed in the derivation is justified from the separation of energy scales in nonrelativistic effective field theories of QCD. Physical implications of the open quantum system approach are also discussed. Finally, I list some open questions for future studies.


title: Canonical Form of Lyapunov Second Method in Mathematical Modelling and Control Design

abstract: The objective of the paper is to put canonical Lyapunov function(CLF), canonizing diffeomorphism (CD) and canonical form of dynamical systems (CFDS), which have led to the generalization of the Lyapunov second method, in perspective of their high efficiency for Mathematical Modelling and Control Design. We show how the symbiosis of the ideas of Henri Poincare and Nikolay Chetaev leads us to CD, CFDS and CLF. Our approach successfully translates into mathematical modelling and control design for special two-angles synchronized longitudinal maneuvering of a thrust-vectored aircraft. The essentially nonlinear five-dimensional mathematical model of the longitudinal flight dynamics of a thrust-vectored aircraft in a wing-body coordinate system with two controls, namely the angular deflections of a movable horizontal stabilizer and a turbojet engine nozzle, is investigated. The wide-sense robust and stable in the large tracking control law is designed. Its core is the hierarchical cascade of two controlling attractor-mediators and two controlling terminal attractors embedded in the extended phase space of the mathematical model of the aircraft longitudinal motion. The detailed demonstration of the elaborated technique of designing wide-sense robust tracking control for the nonlinear multidimensional mathematical model constitutes the quintessence of the paper.


title: Noise and the frontier of quantum supremacy

abstract: Noise is the defining feature of the NISQ era, but it remains unclear if noisy quantum devices are capable of quantum speedups. Quantum supremacy experiments have been a major step forward, but gaps remain between the theory behind these experiments and their actual implementations. In this work we initiate the study of the complexity of quantum random circuit sampling experiments with realistic amounts of noise. Actual quantum supremacy experiments have high levels of uncorrected noise and exponentially decaying fidelities. It is natural to ask if there is any signal of exponential complexity in these highly noisy devices. Surprisingly, we show that it remains hard to compute the output probabilities of noisy random quantum circuits without error correction. More formally, so long as the noise rate of the device is below the error detection threshold, we show it is #P-hard to compute the output probabilities of random circuits with a constant rate of noise per gate. This hardness persists even though these probabilities are exponentially close to uniform. Interestingly these hardness results also have implications for the complexity of experiments in a low-noise setting. The issue here is that prior hardness results for computing output probabilities of random circuits are not robust enough to imprecision to connect with the Stockmeyer argument for hardness of sampling from circuits with constant fidelity. We exponentially improve the robustness of prior results to imprecision, both in the cases of Random Circuit Sampling and BosonSampling. In the latter case we bring the proven hardness within a constant factor in the exponent of the robustness required for hardness of sampling for the first time. We then show that our results are in tension with one another -- the high-noise result implies the low-noise result is essentially optimal, even with generalizations of our techniques.


title: An enhanced parametric nonlinear reduced order model for imperfect structures using Neumann expansion

abstract: We present an enhanced version of the parametric nonlinear reduced order model for shape imperfections in structural dynamics we studied in a previous work [1]. The model is computed intrusively and with no training using information about the nominal geometry of the structure and some user-defined displacement fields representing shape defects, i.e. small deviations from the nominal geometry parametrized by their respective amplitudes. The linear superposition of these artificial displacements describe the defected geometry and can be embedded in the strain formulation in such a way that, in the end, nonlinear internal elastic forces can be expressed as a polynomial function of both these defect fields and the actual displacement field. This way, a tensorial representation of the internal forces can be obtained and, owning the reduction in size of the model given by a Galerkin projection, high simulation speed-ups can be achieved. We show that by adopting a rigorous deformation framework we are able to achieve better accuracy as compared to the previous work. In particular, exploiting Neumann expansion in the definition of the Green-Lagrange strain tensor, we show that our previous model is a lower order approximation with respect to the one we present now. Two numerical examples of a clamped beam and a MEMS gyroscope finally demonstrate the benefits of the method in terms of speed and increased accuracy.


title: Reliability Analysis of Artificial Intelligence Systems Using Recurrent Events Data from Autonomous Vehicles

abstract: Artificial intelligence (AI) systems have become increasingly common and the trend will continue. Examples of AI systems include autonomous vehicles (AV), computer vision, natural language processing, and AI medical experts. To allow for safe and effective deployment of AI systems, the reliability of such systems needs to be assessed. Traditionally, reliability assessment is based on reliability test data and the subsequent statistical modeling and analysis. The availability of reliability data for AI systems, however, is limited because such data are typically sensitive and proprietary. The California Department of Motor Vehicles (DMV) oversees and regulates an AV testing program, in which many AV manufacturers are conducting AV road tests. Manufacturers participating in the program are required to report recurrent disengagement events to California DMV. This information is being made available to the public. In this paper, we use recurrent disengagement events as a representation of the reliability of the AI system in AV, and propose a statistical framework for modeling and analyzing the recurrent events data from AV driving tests. We use traditional parametric models in software reliability and propose a new nonparametric model based on monotonic splines to describe the event process. We develop inference procedures for selecting the best models, quantifying uncertainty, and testing heterogeneity in the event process. We then analyze the recurrent events data from four AV manufacturers, and make inferences on the reliability of the AI systems in AV. We also describe how the proposed analysis can be applied to assess the reliability of other AI systems.


title: Lifshitz scaling, ringing black holes, and superradiance

abstract: We investigate the ringdown waveform and reflectivity of a Lifshitz scalar field around a fixed Schwarzschild black hole. The radial wave equation is modified due to the Lorentz breaking terms, which leads to a diversity of ringdown waveforms. Also, it turns out that Lifshitz waves scattered by the Schwarzschild black hole exhibits superradiance. The Lorentz breaking terms lead to superluminal propagation and high-frequency modes can enter and leave the interior of the Killing horizon where negativity of energy is not prohibited. This allows the Lifshitz waves to carry out additional positive energy to infinity while leaving negative energy inside the Killing horizon, similar to the Penrose process in the ergosphere of a Kerr spacetime. Another interesting phenomenon is emergence of long-lived quasinormal modes, associated with roton-type dispersion relations. These effects drastically modify the greybody factor of a microscopic black hole, whose Hawking temperature is comparable with or higher than the Lifshitz energy scale.


title: cissa(): A MATLAB Function for Signal Extraction

abstract: cissa() is a MATLAB function for signal extraction by Circulant Singular Spectrum Analysis, a procedure proposed in Bogalo et al (2021). cissa() extracts the underlying signals in a time series identifying their frequency of oscillation in an automated way, by just introducing the data and the window length. This solution can be applied to stationary as well as to non-stationary and non-linear time series. Additionally, in this paper, we solve some technical issues regarding the beginning and end of sample data points. We also introduce novel criteria in order to reconstruct the underlying signals grouping some of the extracted components. The output of cissa() is the input of the function group() to reconstruct the desired signals by further grouping the extracted components. group() allows a novel user to create standard signals by automated grouping options while an expert user can decide on the number of groups and their composition. To illustrate its versatility and performance in several fields we include 3 examples: an AM-FM synthetic signal, an example of the physical world given by a voiced speech signal and an economic time series. Possible applications include de-noising, de-seasonalizing, de-trending and extracting business cycles, among others.


title: On singular values of Hankel operators on Bergman spaces

abstract: In this paper, we study the behavior of the singular values of Hankel operators on weighted Bergman spaces $A^2_{\omega _\varphi}$, where $\omega _\varphi= e^{-\varphi}$ and $\varphi$ is a subharmonic function. We consider compact Hankel operators $H_{\overline {\phi}}$, with anti-analytic symbols ${\overline {\phi}}$, and give estimates of the trace of $h(|H_{\overline \phi}|)$ for any convex function $h$. This allows us to give asymptotic estimates of the singular values $(s_n(H_{\overline {\phi}}))_n$ in terms of decreasing rearrangement of $|\phi '|/\sqrt{\Delta \varphi}$. For the radial weights, we first prove that the critical decay of $(s_n(H_{\overline {\phi}}))_n$ is achieved by $(s_n (H_{\overline{z}}))_n$. Namely, we establish that if $s_n(H_{\overline {\phi}})= o (s_n(H_{\overline {z}}))$, then $H_{\overline {\phi}} = 0$. Then, we show that if $\Delta \varphi (z) \asymp \frac{1}{(1-|z|^2)^{2+\beta}}$ with $\beta \geq 0$, then $s_n(H_{\overline {\phi}}) = O(s_n(H_{\overline {z}}))$ if and only if $\phi '$ belongs to the Hardy space $H^p$, where $p= \frac{2(1+\beta)}{2+\beta}$. Finally, we compute the asymptotics of $s_n(H_{\overline {\phi}})$ whenever $ \phi ' \in H^{p }$.


title: Leave No User Behind: Towards Improving the Utility of Recommender Systems for Non-mainstream Users

abstract: In a collaborative-filtering recommendation scenario, biases in the data will likely propagate in the learned recommendations. In this paper we focus on the so-called mainstream bias: the tendency of a recommender system to provide better recommendations to users who have a mainstream taste, as opposed to non-mainstream users. We propose NAECF, a conceptually simple but effective idea to address this bias. The idea consists of adding an autoencoder (AE) layer when learning user and item representations with text-based Convolutional Neural Networks. The AEs, one for the users and one for the items, serve as adversaries to the process of minimizing the rating prediction error when learning how to recommend. They enforce that the specific unique properties of all users and items are sufficiently well incorporated and preserved in the learned representations. These representations, extracted as the bottlenecks of the corresponding AEs, are expected to be less biased towards mainstream users, and to provide more balanced recommendation utility across all users. Our experimental results confirm these expectations, significantly improving the recommendations for non-mainstream users while maintaining the recommendation quality for mainstream users. Our results emphasize the importance of deploying extensive content-based features, such as online reviews, in order to better represent users and items to maximize the de-biasing effect.


title: The climate system and the second law of thermodynamics

abstract: The second law of thermodynamics implies a relationship between the net entropy export by the Earth and its internal irreversible entropy production. The application of this constraint for the purpose of understanding Earth's climate is reviewed. Both radiative processes and material processes are responsible for irreversible entropy production in the climate system. Focusing on material processes, an entropy budget for the climate system is derived which accounts for the multi-phase nature of the hydrological cycle. The entropy budget facilitates a heat-engine perspective of atmospheric circulations that has been used to propose theories for convective updraft velocities, tropical cyclone intensity, and the atmospheric meridional heat transport. Such theories can only be successful, however, if they properly account for the irreversible entropy production associated with water in all its phases in the atmosphere. Irreversibility associated with such moist processes is particularly important in the context of global climate change, for which the concentration of water vapor in the atmosphere is expected to increase, and recent developments toward understanding the response of the atmospheric heat engine to climate change are discussed. Finally, the application of variational approaches to the climate and geophysical flows is briefly reviewed, including the use of equilibrium statistical mechanics to predict behavior of long-lived coherent structures, and the controversial maximum entropy production principle.


title: Inference of the Selective Auditory Attention using Sequential LMMSE Estimation

abstract: Attentive listening in a multispeaker environment such as a cocktail party requires suppression of the interfering speakers and the noise around. People with normal hearing perform remarkably well in such situations. Analysis of the cortical signals using electroencephalography (EEG) has revealed that the EEG signals track the envelope of the attended speech stronger than that of the interfering speech. This has enabled the development of algorithms that can decode the selective attention of a listener in controlled experimental settings. However, often these algorithms require longer trial duration and computationally expensive calibration to obtain a reliable inference of attention. In this paper, we present a novel framework to decode the attention of a listener within trial durations of the order of two seconds. It comprises of three modules: 1) Dynamic estimation of the temporal response functions (TRF) in every trial using a sequential linear minimum mean squared error (LMMSE) estimator, 2) Extract the N1-P2 peak of the estimated TRF that serves as a marker related to the attentional state and 3) Obtain a probabilistic measure of the attentional state using a support vector machine followed by a logistic regression. The efficacy of the proposed decoding framework was evaluated using EEG data collected from 27 subjects. The total number of electrodes required to infer the attention was four: One for the signal estimation, one for the noise estimation and the other two being the reference and the ground electrodes. Our results make further progress towards the realization of neuro-steered hearing aids.


title: Real-time rendering of complex fractals

abstract: This chapter describes how to use intersection and closest-hit shaders to implement real-time visualizations of complex fractals using distance functions. The Mandelbulb and Julia Sets are used as examples.


title: Near-Optimal Offline Reinforcement Learning via Double Variance Reduction

abstract: We consider the problem of offline reinforcement learning (RL) -- a well-motivated setting of RL that aims at policy optimization using only historical data. Despite its wide applicability, theoretical understandings of offline RL, such as its optimal sample complexity, remain largely open even in basic settings such as \emph{tabular} Markov Decision Processes (MDPs). In this paper, we propose Off-Policy Double Variance Reduction (OPDVR), a new variance reduction based algorithm for offline RL. Our main result shows that OPDVR provably identifies an $\epsilon$-optimal policy with $\widetilde{O}(H^2/d_m\epsilon^2)$ episodes of offline data in the finite-horizon stationary transition setting, where $H$ is the horizon length and $d_m$ is the minimal marginal state-action distribution induced by the behavior policy. This improves over the best known upper bound by a factor of $H$. Moreover, we establish an information-theoretic lower bound of $\Omega(H^2/d_m\epsilon^2)$ which certifies that OPDVR is optimal up to logarithmic factors. Lastly, we show that OPDVR also achieves rate-optimal sample complexity under alternative settings such as the finite-horizon MDPs with non-stationary transitions and the infinite horizon MDPs with discounted rewards.


title: Vehicle trajectory prediction in top-view image sequences based on deep learning method

abstract: Annually, a large number of injuries and deaths around the world are related to motor vehicle accidents. This value has recently been reduced to some extent, via the use of driver-assistance systems. Developing driver-assistance systems (i.e., automated driving systems) can play a crucial role in reducing this number. Estimating and predicting surrounding vehicles' movement is essential for an automated vehicle and advanced safety systems. Moreover, predicting the trajectory is influenced by numerous factors, such as drivers' behavior during accidents, history of the vehicle's movement and the surrounding vehicles, and their position on the traffic scene. The vehicle must move over a safe path in traffic and react to other drivers' unpredictable behaviors in the shortest time. Herein, to predict automated vehicles' path, a model with low computational complexity is proposed, which is trained by images taken from the road's aerial image. Our method is based on an encoder-decoder model that utilizes a social tensor to model the effect of the surrounding vehicles' movement on the target vehicle. The proposed model can predict the vehicle's future path in any freeway only by viewing the images related to the history of the target vehicle's movement and its neighbors. Deep learning was used as a tool for extracting the features of these images. Using the HighD database, an image dataset of the road's aerial image was created, and the model's performance was evaluated on this new database. We achieved the RMSE of 1.91 for the next 5 seconds and found that the proposed method had less error than the best path-prediction methods in previous studies.


title: Manifold Repairing, Reconstruction and Denoising from Scattered Data in High-Dimension

abstract: We consider a problem of great practical interest: the repairing and recovery of a low-dimensional manifold embedded in high-dimensional space from noisy scattered data. Suppose that we observe a point cloud sampled from the low-dimensional manifold, with noise, and let us assume that there are holes in the data. Can we recover missing information inside the holes? While in low-dimension the problem was extensively studied, manifold repairing in high dimension is still an open problem. We introduce a new approach, called Repairing Manifold Locally Optimal Projection (R-MLOP), that expands the MLOP method introduced by Faigenbaum-Golovin et al. in 2020, to cope with manifold repairing in low and high-dimensional cases. The proposed method can deal with multiple holes in a manifold. We prove the validity of the proposed method, and demonstrate the effectiveness of our approach by considering different manifold topologies, for single and multiple holes repairing, in low and high dimensions.


title: Distributed Conditional Generative Adversarial Networks (GANs) for Data-Driven Millimeter Wave Communications in UAV Networks

abstract: In this paper, a novel framework is proposed to perform data-driven air-to-ground (A2G) channel estimation for millimeter wave (mmWave) communications in an unmanned aerial vehicle (UAV) wireless network. First, an effective channel estimation approach is developed to collect mmWave channel information, allowing each UAV to train a stand-alone channel model via a conditional generative adversarial network (CGAN) along each beamforming direction. Next, in order to expand the application scenarios of the trained channel model into a broader spatial-temporal domain, a cooperative framework, based on a distributed CGAN architecture, is developed, allowing each UAV to collaboratively learn the mmWave channel distribution in a fully-distributed manner. To guarantee an efficient learning process, necessary and sufficient conditions for the optimal UAV network topology that maximizes the learning rate for cooperative channel modeling are derived, and the optimal CGAN learning solution per UAV is subsequently characterized, based on the distributed network structure. Simulation results show that the proposed distributed CGAN approach is robust to the local training error at each UAV. Meanwhile, a larger airborne network size requires more communication resources per UAV to guarantee an efficient learning rate. The results also show that, compared with a stand-alone CGAN without information sharing and two other distributed schemes, namely: A multi-discriminator CGAN and a federated CGAN method, the proposed distributed CGAN approach yields a higher modeling accuracy while learning the environment, and it achieves a larger average data rate in the online performance of UAV downlink mmWave communications.


title: Continuous Wasserstein-2 Barycenter Estimation without Minimax Optimization

abstract: Wasserstein barycenters provide a geometric notion of the weighted average of probability measures based on optimal transport. In this paper, we present a scalable algorithm to compute Wasserstein-2 barycenters given sample access to the input measures, which are not restricted to being discrete. While past approaches rely on entropic or quadratic regularization, we employ input convex neural networks and cycle-consistency regularization to avoid introducing bias. As a result, our approach does not resort to minimax optimization. We provide theoretical analysis on error bounds as well as empirical evidence of the effectiveness of the proposed approach in low-dimensional qualitative scenarios and high-dimensional quantitative experiments.


title: Debiased Inference on Heterogeneous Quantile Treatment Effects with Regression Rank-Scores

abstract: Understanding treatment effect heterogeneity is vital to many scientific fields because the same treatment may affect different individuals differently. Quantile regression provides a natural framework for modeling such heterogeneity. We propose a new method for inference on heterogeneous quantile treatment effects in the presence of high-dimensional covariates. Our estimator combines an $\ell_1$-penalized regression adjustment with a quantile-specific bias correction scheme based on rank scores. We study the theoretical properties of this estimator, including weak convergence and semiparametric efficiency of the estimated heterogeneous quantile treatment effect process. We illustrate the finite-sample performance of our approach through simulations and an empirical example, dealing with the differential effect of statin usage for lowering low-density lipoprotein cholesterol levels for the Alzheimer's disease patients who participated in the UK Biobank study.


title: LSSED: a large-scale dataset and benchmark for speech emotion recognition

abstract: Speech emotion recognition is a vital contributor to the next generation of human-computer interaction (HCI). However, current existing small-scale databases have limited the development of related research. In this paper, we present LSSED, a challenging large-scale english speech emotion dataset, which has data collected from 820 subjects to simulate real-world distribution. In addition, we release some pre-trained models based on LSSED, which can not only promote the development of speech emotion recognition, but can also be transferred to related downstream tasks such as mental health analysis where data is extremely difficult to collect. Finally, our experiments show the necessity of large-scale datasets and the effectiveness of pre-trained models. The dateset will be released on https://github.com/tobefans/LSSED.


title: Measuring solar neutrinos over Gigayear timescales with Paleo Detectors

abstract: Measuring the solar neutrino flux over gigayear timescales could provide a new window to inform the Solar Standard Model as well as studies of the Earth's long-term climate. We demonstrate the feasibility of measuring the time-evolution of the $^8$B solar neutrino flux over gigayear timescales using paleo detectors, naturally occurring minerals which record neutrino-induced recoil tracks over geological times. We explore suitable minerals and identify track lengths of 15--30 nm to be a practical window to detect the $^8$B solar neutrino flux. A collection of ultra-radiopure minerals of different ages, each some 0.1 kg by mass, can be used to probe the rise of the $^8$B solar neutrino flux over the recent gigayear of the Sun's evolution. We also show that models of the solar abundance problem can be distinguished based on the time-integrated tracks induced by the $^8$B solar neutrino flux.

categories: hep-ph astro-ph.HE astro-ph.IM astro-ph.SR hep-ex


title: A Novel Magnetic Respiratory Sensor for Human Healthcare

abstract: Breathing is vital to life. Therefore, the real-time monitoring of breathing pattern of a patient is crucial to respiratory rehabilitation therapies such as magnetic resonance exams for respiratory-triggered imaging, chronic pulmonary disease treatment, and synchronized functional electrical stimulation. While numerous respiratory devices have been developed, they are often in direct contact with a patient, which can yield inaccurate or limited data. In this study, we developed a novel, non-invasive, and contactless magnetic sensing platform that can precisely monitor breathing, movement, or sleep patterns of a patient, thus providing efficient monitoring at a clinic or home. A magneto-LC resonance (MLCR) sensor converts the magnetic oscillations generated by breathing of the patient into an impedance spectrum, which allows for a deep analysis of breath variation to identify respiratory-related diseases like COVID-19. Owing to its ultrahigh sensitivity, the MLCR sensor yields a distinct breathing pattern for each patient tested. The sensor also provides an accurate measure of the strength of breath at multiple stages as well as anomalous variations in respiratory rate and amplitude. This suggests that the MLCR sensor can detect symptoms of COVID-19 in a patient due to shortness of breath or difficulty breathing as well as track the progress of the disease in real time.


title: The Multilingual TEDx Corpus for Speech Recognition and Translation

abstract: We present the Multilingual TEDx corpus, built to support speech recognition (ASR) and speech translation (ST) research across many non-English source languages. The corpus is a collection of audio recordings from TEDx talks in 8 source languages. We segment transcripts into sentences and align them to the source-language audio and target-language translations. The corpus is released along with open-sourced code enabling extension to new talks and languages as they become available. Our corpus creation methodology can be applied to more languages than previous work, and creates multi-way parallel evaluation sets. We provide baselines in multiple ASR and ST settings, including multilingual models to improve translation performance for low-resource language pairs.


title: Effective out-of-plane g-factor in strained-Ge/SiGe quantum dots

abstract: Recently, lithographic quantum dots in strained-Ge/SiGe have become a promising candidate for quantum computation, with a remarkably quick progression from demonstration of a quantum dot to qubit logic demonstrations. Here we present a measurement of the out-of-plane $g$-factor for single-hole quantum dots in this material. As this is a single-hole measurement, this is the first experimental result that avoids the strong orbital effects present in the out-of-plane configuration. In addition to verifying the expected $g$-factor anisotropy between in-plane and out-of-plane magnetic ($B$)-fields, variations in the $g$-factor dependent on the occupation of the quantum dot are observed. These results are in good agreement with calculations of the $g$-factor using the heavy- and light-hole spaces of the Luttinger Hamiltonian, especially the first two holes, showing a strong spin-orbit coupling and suggesting dramatic $g$-factor tunability through both the $B$-field and the charge state.


title: Ansatz-Independent Variational Quantum Classifier

abstract: The paradigm of variational quantum classifiers (VQCs) encodes \textit{classical information} as quantum states, followed by quantum processing and then measurements to generate classical predictions. VQCs are promising candidates for efficient utilization of a near-term quantum device: classifiers involving $M$-dimensional datasets can be implemented with only $\lceil \log_2 M \rceil$ qubits by using an amplitude encoding. A general framework for designing and training VQCs, however, has not been proposed, and a fundamental understanding of its power and analytical relationships with classical classifiers are not well understood. An encouraging specific embodiment of VQCs, quantum circuit learning (QCL), utilizes an ansatz: it expresses the quantum evolution operator as a circuit with a predetermined topology and parametrized gates; training involves learning the gate parameters through optimization. In this letter, we first address the open questions about VQCs and then show that they, including QCL, fit inside the well-known kernel method. Based on such correspondence, we devise a design framework of efficient ansatz-independent VQCs, which we call the unitary kernel method (UKM): it directly optimizes the unitary evolution operator in a VQC. Thus, we show that the performance of QCL is bounded from above by the UKM. Next, we propose a variational circuit realization (VCR) for designing efficient quantum circuits for a given unitary operator. By combining the UKM with the VCR, we establish an efficient framework for constructing high-performing circuits. We finally benchmark the relatively superior performance of the UKM and the VCR via extensive numerical simulations on multiple datasets.


title: A Speaker Verification Backend with Robust Performance across Conditions

abstract: In this paper, we address the problem of speaker verification in conditions unseen or unknown during development. A standard method for speaker verification consists of extracting speaker embeddings with a deep neural network and processing them through a backend composed of probabilistic linear discriminant analysis (PLDA) and global logistic regression score calibration. This method is known to result in systems that work poorly on conditions different from those used to train the calibration model. We propose to modify the standard backend, introducing an adaptive calibrator that uses duration and other automatically extracted side-information to adapt to the conditions of the inputs. The backend is trained discriminatively to optimize binary cross-entropy. When trained on a number of diverse datasets that are labeled only with respect to speaker, the proposed backend consistently and, in some cases, dramatically improves calibration, compared to the standard PLDA approach, on a number of held-out datasets, some of which are markedly different from the training data. Discrimination performance is also consistently improved. We show that joint training of the PLDA and the adaptive calibrator is essential -- the same benefits cannot be achieved when freezing PLDA and fine-tuning the calibrator. To our knowledge, the results in this paper are the first evidence in the literature that it is possible to develop a speaker verification system with robust out-of-the-box performance on a large variety of conditions.


title: Deep Convolutional Neural Networks to Predict Mutual Coupling Effects in Metasurfaces

abstract: Metasurfaces have provided a novel and promising platform for the realization of compact and large-scale optical devices. The conventional metasurface design approach assumes periodic boundary conditions for each element, which is inaccurate in most cases since the near-field coupling effects between elements will change when surrounded by non-identical structures. In this paper, we propose a deep learning approach to predict the actual electromagnetic (EM) responses of each target meta-atom placed in a large array with near-field coupling effects taken into account. The predicting neural network takes the physical specifications of the target meta-atom and its neighbors as input, and calculates its phase and amplitude in milliseconds. This approach can be applied to explain metasurfaces' performance deterioration caused by mutual coupling and further used to optimize their efficiencies once combined with optimization algorithms. To demonstrate the efficacy of this methodology, we obtain large improvements in efficiency for a beam deflector and a metalens over the conventional design approach. Moreover, we show the correlations between a metasurface's performance and its design errors caused by mutual coupling are not bound to certain specifications (materials, shapes, etc.). As such, we envision that this approach can be readily applied to explore the mutual coupling effects and improve the performance of various metasurface designs.


title: Profinite genus of fundamental groups of compact flat manifolds with the cyclic holonomy group of square-free order

abstract: In this article we study the extent to which an $n$-dimensional compact flat manifold with the cyclic holonomy group of square-free order may be distinguished by the finite quotients of its fundamental group. In particular, we display a formula for the cardinality of profinite genus of the fundamental group of an $n$-dimensional compact flat manifold with the cyclic holonomy group of square-free order.


title: Optimization of the high-frequency magnetoimpedance response in melt-extracted Co-rich microwires through novel multiple-step Joule heating

abstract: The optimization of high frequency giant magnetoimpedance (GMI) effect and its magnetic field sensitivity in melt-extracted Co69.25Fe4.25Si13B12.5Nb1 amorphous microwires, through a multi-step Joule annealing (MSA) technique, was systematically studied. The surface morphology, microstructure, surface magnetic property, and high frequency GMI response of the Co-rich microwires were explored using scanning electron microscopy (SEM), magneto-optical Kerr effect (MOKE) magnetometry, transmission electron microscopy (TEM), and impedance analyzer, respectively. An initial dc current (idc) of 20 mA, which was then increased by 20 mA at every time-step (10 min) up to 300 mA, was applied to the microwires. The MSA of 20 mA to 100 mA remarkably improved the GMI ratio and its field sensitivity up to 760% (1.75 time of that of the as-prepared), and 925%/Oe (more than 17.92 times of that of the as-prepared) at an operating frequency of 20 MHz, respectively. Our study indicates that the MSA technique can enhance the microstructures and the surface magnetic domain structures of the Co-rich magnetic microwires, giving rise to the GMI enhancement. This technique is suitable for improving the GMI sensitivity at small magnetic fields, which is highly promising for biomedical sensing and healthcare monitoring.


title: MANA: Microarchitecting an Instruction Prefetcher

abstract: L1 instruction (L1-I) cache misses are a source of performance bottleneck. Sequential prefetchers are simple solutions to mitigate this problem; however, prior work has shown that these prefetchers leave considerable potentials uncovered. This observation has motivated many researchers to come up with more advanced instruction prefetchers. In 2011, Proactive Instruction Fetch (PIF) showed that a hardware prefetcher could effectively eliminate all of the instruction-cache misses. However, its enormous storage cost makes it an impractical solution. Consequently, reducing the storage cost was the main research focus in the instruction prefetching in the past decade. Several instruction prefetchers, including RDIP and Shotgun, were proposed to offer PIF-level performance with significantly lower storage overhead. However, our findings show that there is a considerable performance gap between these proposals and PIF. While these proposals use different mechanisms for instruction prefetching, the performance gap is largely not because of the mechanism, and instead, is due to not having sufficient storage. Prior proposals suffer from one or both of the following shortcomings: (1) a large number of metadata records to cover the potential, and (2) a high storage cost of each record. The first problem causes metadata miss, and the second problem prohibits the prefetcher from storing enough records within reasonably-sized storage.


title: Insights into the Emergence of Mass from Studies of Pion and Kaon Structure

abstract: There are two mass generating mechanisms in the standard model of particle physics (SM). One is related to the Higgs boson and fairly well understood. The other is embedded in quantum chromodynamics (QCD), the SM's strong interaction piece; and although responsible for emergence of the roughly 1 GeV mass scale that characterises the proton and hence all observable matter, the source and impacts of this emergent hadronic mass (EHM) remain puzzling. As bound states seeded by a valence-quark and -antiquark, pseudoscalar mesons present a simpler problem in quantum field theory than that associated with the nucleon. Consequently, there is a large array of robust predictions for pion and kaon properties whose empirical validation will provide a clear window onto many effects of both mass generating mechanisms and the constructive interference between them. This has now become significant because new-era experimental facilities, in operation, construction, or planning, are capable of conducting such tests and thereby contributing greatly to resolving the puzzles of EHM. These aspects of experiment, phenomenology, and theory, along with contemporary successes and challenges, are sketched herein, simultaneously highlighting the potential gains that can accrue from a coherent effort aimed at finally reaching an understanding of the character and structure of Nature's Nambu-Goldstone modes.


title: Novel one-shot inner bounds for unassisted fully quantum channels via rate splitting

abstract: We prove the first non-trivial one-shot inner bounds for sending quantum information over an entanglement unassisted two-sender quantum multiple access channel (QMAC) and an unassisted two-sender two-receiver quantum interference channel (QIC). Previous works only studied the unassisted QMAC in the limit of many independent and identical uses of the channel also known as the asymptotic iid limit, and did not study the unassisted QIC at all. We employ two techniques, rate splitting and successive cancellation}, in order to obtain our inner bound. Rate splitting was earlier used to obtain inner bounds, avoiding time sharing, for classical channels in the asymptotic iid setting. Our main technical contribution is to extend rate splitting from the classical asymptotic iid setting to the quantum one-shot setting. In the asymptotic iid limit our one-shot inner bound for QMAC approaches the rate region of Yard, Devetak and Hayden. For the QIC we get novel non-trivial rate regions in the asymptotic iid setting. All our results also extend to the case where limited entanglement assistance is provided, in both one-shot and asymptotic iid settings. The limited entanglement results for one-setting for both QMAC and QIC are new. For the QIC the limited entanglement results are new even in the asymptotic iid setting.


title: Automatic analysis of artistic paintings using information-based measures

abstract: The artistic community is increasingly relying on automatic computational analysis for authentication and classification of artistic paintings. In this paper, we identify hidden patterns and relationships present in artistic paintings by analysing their complexity, a measure that quantifies the sum of characteristics of an object. Specifically, we apply Normalized Compression (NC) and the Block Decomposition Method (BDM) to a dataset of 4,266 paintings from 91 authors and examine the potential of these information-based measures as descriptors of artistic paintings. Both measures consistently described the equivalent types of paintings, authors, and artistic movements. Moreover, combining the NC with a measure of the roughness of the paintings creates an efficient stylistic descriptor. Furthermore, by quantifying the local information of each painting, we define a fingerprint that describes critical information regarding the artists' style, their artistic influences, and shared techniques. More fundamentally, this information describes how each author typically composes and distributes the elements across the canvas and, therefore, how their work is perceived. Finally, we demonstrate that regional complexity and two-point height difference correlation function are useful auxiliary features that improve current methodologies in style and author classification of artistic paintings. The whole study is supported by an extensive website (http://panther.web.ua.pt) for fast author characterization and authentication.


title: CHAMPION: Chalmers Hierarchical Atomic, Molecular, Polymeric & Ionic Analysis Toolkit

abstract: We present CHAMPION: a software developed to automatically detect time-dependent bonds between atoms based on their dynamics, classify the local graph topology around them, and analyze the physicochemical properties of these topologies by statistical physics. In stark contrast to methodologies where bonds are detected based on static conditions such as cut-off distances, CHAMPION considers pairs of atoms to be bound only if they move together and act as a bound pair over time. Furthermore, the time-dependent global bond graph is possible to split into dynamically shifting connected components or subgraphs around a certain chemical motif and thereby allow the physicochemical properties of each such topology to be analyzed by statistical physics. Applicable to condensed matter and liquids in general, and electrolytes in particular, this allows both quantitative and qualitative descriptions of local structure, as well as dynamical processes such as speciation and diffusion. We present here a detailed overview of CHAMPION, including its underlying methodology, implementation and capabilities.


title: A Basis Approach to Surface Clustering

abstract: This paper presents a novel method for clustering surfaces. The proposal involves first using basis functions in a tensor product to smooth the data and thus reduce the dimension to a finite number of coefficients, and then using these estimated coefficients to cluster the surfaces via the k-means algorithm. An extension of the algorithm to clustering tensors is also discussed. We show that the proposed algorithm exhibits the property of strong consistency, with or without measurement errors, in correctly clustering the data as the sample size increases. Simulation studies suggest that the proposed method outperforms the benchmark k-means algorithm which uses the original vectorized data. In addition, an EGG real data example is considered to illustrate the practical application of the proposal.


title: A privacy-preserving approach to streaming eye-tracking data

abstract: Eye-tracking technology is being increasingly integrated into mixed reality devices. Although critical applications are being enabled, there are significant possibilities for violating user privacy expectations. We show that there is an appreciable risk of unique user identification even under natural viewing conditions in virtual reality. This identification would allow an app to connect a user's personal ID with their work ID without needing their consent, for example. To mitigate such risks we propose a framework that incorporates gatekeeping via the design of the application programming interface and via software-implemented privacy mechanisms. Our results indicate that these mechanisms can reduce the rate of identification from as much as 85% to as low as 30%. The impact of introducing these mechanisms is less than 1.5$^\circ$ error in gaze position for gaze prediction. Gaze data streams can thus be made private while still allowing for gaze prediction, for example, during foveated rendering. Our approach is the first to support privacy-by-design in the flow of eye-tracking data within mixed reality use cases.


title: Secret Key Agreement and Secure Omniscience of Tree-PIN Source with Linear Wiretapper

abstract: While the wiretap secret key capacity remains unknown for general source models even in the two-user case, we obtained a single-letter characterization for a large class of multi-user source models with a linear wiretapper who can observe any linear combinations of the source. We introduced the idea of irreducible sources to show existence of an optimal communication scheme that achieves perfect omniscience with minimum leakage of information to the wiretapper. This implies a duality between the problems of wiretap secret key agreement and secure omniscience, and such duality potentially holds for more general sources.


title: Gyro-Kinematic Ages for around 30,000 Kepler Stars

abstract: Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, and crucially K and M dwarfs, using their rotation and kinematic properties. We apply the simple assumption that the velocity dispersion of stars increases over time and adopt an age--velocity--dispersion relation (AVR) to estimate average stellar ages for groupings of coeval stars. We calculate the vertical velocity dispersion of stars in bins of absolute magnitude, temperature, rotation period, and Rossby number and then convert velocity dispersion to kinematic age via an AVR. Using this method, we estimate gyro-kinematic ages for 29,949 Kepler stars with measured rotation periods. We are able to estimate ages for clusters and asteroseismic stars with an RMS of 1.22 Gyr and 0.26 Gyr respectively. With our Astraea machine learning algorithm, which predicts rotation periods, we suggest a new selection criterion (a weight of 0.15) to increase the size of the McQuillian et al. (2014) catalog of Kepler rotation periods by up to 25%. Using predicted rotation periods, we estimated gyro-kinematic ages for stars without measured rotation periods and found promising results by comparing 12 detailed age--element abundance trends with literature values.


title: Reducing number fluctuations in an ultracold atomic sample using Faraday rotation and iterative feedback

abstract: We demonstrate a method to reduce number fluctuations in an ultracold atomic sample using real-time feedback. By measuring the Faraday rotation of an off-resonant probe laser beam with a pair of avalanche photodetectors in a polarimetric setup we produce a proxy for the number of atoms in the sample. We iteratively remove a fraction of the excess atoms from the sample to converge on a target proxy value in a way that is insensitive to environmental perturbations and robust to errors in light polarization. Using absorption imaging for out-of-loop verification, we demonstrate a reduction in the number fluctuations from $3\%$ to $0.45\%$ for samples at a temperature of 16.4 $\mu$K over the time-scale of several hours which is limited by temperature fluctuations, beam pointing noise, and photon shot noise.


title: starry_process: Interpretable Gaussian processes for stellar light curves

abstract: In this note we present the starry_process code, which implements an interpretable Gaussian process (GP) for modeling variability in stellar light curves. As dark starspots rotate in and out of view, the total flux received from a distant star will change over time. Unresolved flux time series therefore encode information about the spatial structure of features on the stellar surface. The starry_process software package allows one to easily model the flux variability due to starspots, whether one is interested in understanding the properties of these spots or marginalizing over the stellar variability when it is treated as a nuisance signal. The main difference between the GP implemented here and typical GPs used to model stellar variability is the explicit dependence of our GP on physical properties of the star, such as its period, inclination, and limb darkening coefficients, and on properties of the spots, such as their radius and latitude distributions. This code is the Python implementation of the interpretable GP algorithm developed in Luger, Foreman-Mackey, and Hedges (2021).


title: Safe Search for Stackelberg Equilibria in Extensive-Form Games

abstract: Stackelberg equilibrium is a solution concept in two-player games where the leader has commitment rights over the follower. In recent years, it has become a cornerstone of many security applications, including airport patrolling and wildlife poaching prevention. Even though many of these settings are sequential in nature, existing techniques pre-compute the entire solution ahead of time. In this paper, we present a theoretically sound and empirically effective way to apply search, which leverages extra online computation to improve a solution, to the computation of Stackelberg equilibria in general-sum games. Instead of the leader attempting to solve the full game upfront, an approximate "blueprint" solution is first computed offline and is then improved online for the particular subgames encountered in actual play. We prove that our search technique is guaranteed to perform no worse than the pre-computed blueprint strategy, and empirically demonstrate that it enables approximately solving significantly larger games compared to purely offline methods. We also show that our search operation may be cast as a smaller Stackelberg problem, making our method complementary to existing algorithms based on strategy generation.


title: Improved Treatment of Host-Galaxy Correlations in Cosmological Analyses With Type Ia Supernovae

abstract: Improving the use of Type Ia supernovae (SNIa) as standard candles requires a better approach to incorporate the relationship between SNIa and the properties of their host galaxies. Using a spectroscopically-confirmed sample of $\sim$1600 SNIa, we develop the first empirical model of underlying populations for SNIa light-curve properties that includes their dependence on host-galaxy stellar mass. These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNIa-host correlations, and we validate this technique on simulated data samples. We recover the input relationship between SNIa luminosity and host-galaxy stellar mass (the mass step, $\gamma$) to within 0.004 mags, which is a factor of 5 improvement over the previous method that results in a $\gamma$-bias of ${\sim}0.02$. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data ($\gamma = 0.017 \pm 0.008$), and for simulations ($\gamma =0.006 \pm 0.007$). Analysing simulated SNIa, the biases on the dark energy equation-of-state, $w$, vary from $\Delta w = 0.006(5)$ to $0.010(5)$ with our new BBC method; these biases are significantly smaller than the $0.02(5)$ $w$-bias using previous BBC methods that ignore SNIa-host correlations.


title: Kinetics of Intrinsic Stress in Nanocrystalline Films

abstract: Conventional polycrystalline materials acquire high levels of intrinsic mechanical stress (ranging from MPa to a few GPa) during preparation and use, but this stress decays quickly (~minutes) to small residual values (~kPa) under standard resting conditions. Nanocrystalline materials reach similar or even higher levels of intrinsic stress, but surprisingly retain a significant portion of this stress over much longer time scales (~hours). This behavior directly contradicts current theoretical models that predict stress relaxation through diffusive currents. Diffusive currents, which flow mainly on the surfaces of grains, are expected to produce faster relaxation kinetics when the stress to be released is higher. In this work, we study the kinetics of intrinsic stress relaxation in nanocrystalline films and identify the limitations of this process as a preliminary step towards designing a strategy for high-stress stabilization in nanostructured systems.


title: Averaging with the Divisor Function: $\ell^p$-improving and Sparse Bounds

abstract: We study averages along the integers using the divisor function $d(n)$, and defined as $$K_N f (x) = \frac{1}{D(N)} \sum _{n \leq N} d(n) \,f(x+n) , $$ where $D(N) = \sum _{n=1} ^N d(n) $. We shall show that these averages satisfy a uniform, scale free $\ell^p$-improving estimate for $p \in (1,2)$, that is $$ \left( \frac{1}{N} \sum |K_Nf|^{p'} \right)^{1/p'} \lesssim \left(\frac{1}{N} \sum |f|^p \right)^{1/p} $$ as long as $f$ is supported on $[0,N]$. We also show that the associated maximal function $K^*f = \sup_N |K_N f|$ satisfies $(p,p)$ sparse founds for $p \in (1,2)$, which implies that $K^*$ is bounded on $\ell ^p (w)$ for $p \in (1, \infty )$, for all weights $w$ in the Muckenhoupt $A_p$ class.


title: An algebraic treatment of the Askey biorthogonal polynomials on the unit circle

abstract: A joint algebraic interpretation of the biorthogonal Askey polynomials on the unit circle and of the orthogonal Jacobi polynomials is offered. It ties their bispectral properties to an algebra called the meta-Jacobi algebra $m\mathfrak{J}$.


title: A metaheuristic for crew scheduling in a pickup-and-delivery problem with time windows

abstract: A vehicle routing and crew scheduling problem (VRCSP) consists of simultaneously planning the routes of a fleet of vehicles and scheduling the crews, where the vehicle-crew correspondence is not fixed through time. This allows a greater planning flexibility and a more efficient use of the fleet, but in counterpart, a high synchronisation is demanded. In this work, we present a VRCSP where pickup-and-delivery requests with time windows have to be fulfilled over a given planning horizon by using trucks and drivers. Crews can be composed of 1 or 2 drivers and any of them can be relieved in a given set of locations. Moreover, they are allowed to travel among locations with non-company shuttles, at an additional cost that is minimised. As our problem considers distinct routes for trucks and drivers, we have an additional flexibility not contemplated in other previous VRCSP given in the literature where a crew is handled as an indivisible unit. We tackle this problem with a two-stage sequential approach: a set of truck routes is computed in the first stage and a set of driver routes consistent with the truck routes is obtained in the second one. We design and evaluate the performance of a metaheuristic based algorithm for the latter stage. Our algorithm is mainly a GRASP with a perturbation procedure that allows reusing solutions already found in case the search for new solutions becomes difficult. This procedure together with other to repair infeasible solutions allow us to find high-quality solutions on instances of 100 requests spread across 15 cities with a fleet of 12-32 trucks (depending on the planning horizon) in less than an hour. We also conclude that the possibility of carrying an additional driver leads to a decrease of the cost of external shuttles by about 60% on average with respect to individual crews and, in some cases, to remove this cost completely.


title: Optimization of the Variational Quantum Eigensolver for Quantum Chemistry Applications

abstract: This work studies the variational quantum eigensolver algorithm, designed to determine the ground state of a quantum mechanical system by combining classical and quantum hardware. Methods of reducing the number of required qubit manipulations, prone to induce errors, for the variational quantum eigensolver are studied. We formally justify the qubit removal process as sketched by Bravyi, Gambetta, Mezzacapo and Temme [arXiv:1701.08213 (2017)]. Furthermore, different classical optimization and entangling methods, both gate based and native, are surveyed by computing ground state energies of H$_2$ and LiH. This paper aims to provide performance-based recommendations for entangling methods and classical optimization methods. Analyzing the VQE problem is complex, where the optimization algorithm, the method of entangling, and the dimensionality of the search space all interact. In specific cases however, concrete results can be shown, and an entangling method or optimization algorithm can be recommended over others. In particular we find that for high dimensionality (many qubits and/or entanglement depth) certain classical optimization algorithms outperform others in terms of energy error.


title: Dipolar spin relaxation of divacancy qubits in silicon carbide

abstract: Divacancy spins in silicon carbide implement qubits with outstanding characteristics and capabilities in an industrial semiconductor host. On the other hand, there are still numerous open questions about the physics of divacancy point defects, for instance, spin relaxation has not been thoroughly studied yet. Here, we carry out a theoretical study on environmental spin induced spin relaxation processes of divacancy qubits in 4H-SiC. We reveal all the relevant magnetic field values where the longitudinal spin relaxation time T$_1$ drops resonantly due to the coupling to either nuclear spins or electron spins. We quantitatively analyze the dependence of the T$_1$ time on the concentration of point defect spins and the applied magnetic field in the most relevant cases and provide an analytical expression. We demonstrate that dipolar spin relaxation plays a significant role both in as-grown and ion implanted samples and it often limits the coherence time in 4H-SiC.


title: Implementation of efficient quantum search algorithms on NISQ computers

abstract: Despite the advent of Grover's algorithm for the unstructured search, its successful implementation on near-term quantum devices is still limited. We apply three strategies to reduce the errors associated with implementing quantum search algorithms. Our improved search algorithms have been implemented on the IBM quantum processors. Using them, we demonstrate three- and four-qubit search algorithm with higher average success probabilities compared to previous works. We present the successful execution of the five-qubit search on the IBM quantum processor for the first time. The results have been benchmarked using degraded ratio, which is the ratio between the experimental and the theoretical success probabilities. The fast decay of the degraded ratio supports our divide-and-conquer strategy. Our proposed strategies are also useful for implementation of quantum search algorithms in the post-NISQ era.


title: Adaptive Frequency Band Analysis for Functional Time Series

abstract: The frequency-domain properties of nonstationary functional time series often contain valuable information. These properties are characterized through its time-varying power spectrum. Practitioners seeking low-dimensional summary measures of the power spectrum often partition frequencies into bands and create collapsed measures of power within bands. However, standard frequency bands have largely been developed through manual inspection of time series data and may not adequately summarize power spectra. In this article, we propose a framework for adaptive frequency band estimation of nonstationary functional time series that optimally summarizes the time-varying dynamics of the series. We develop a scan statistic and search algorithm to detect changes in the frequency domain. We establish theoretical properties of this framework and develop a computationally-efficient implementation. The validity of our method is also justified through numerous simulation studies and an application to analyzing electroencephalogram data in participants alternating between eyes open and eyes closed conditions.


title: Toward a more stringent test of gravity with redshift space power spectrum: simultaneous probe of growth and amplitude of large-scale structure

abstract: Redshift-space distortions (RSD) offers an exciting opportunity to test the gravity on cosmological scales. In the presence of galaxy bias, however, the RSD measurement at large scales, where the linear theory prediction is safely applied, is known to exhibit a degeneracy between the parameters of structure growth f and fluctuation amplitude sigma8, and one can only constrain the parameters in the form of fsigma8. In order to disentangle this degeneracy, in this paper, we go beyond the linear theory, and consider the model of RSD applicable to a weakly nonlinear regime. Based on the Fisher matrix analysis, we show explicitly that the degeneracy of the parameter fsigma8 can be broken, and sigma8 is separately estimated in the presence of galaxy bias. Performing further the Markov chain Monte Carlo analysis, we verify that our model correctly reproduces the fiducial values of fsigma8 and sigma8, with the statistical errors consistent with those estimated from the Fisher matrix analysis. We show that upcoming galaxy survey of the stage-IV class can unambiguously determine sigma8 at the precision down to 10% at higher redshifts even if we restrict the accessible scales to k<0.16h/Mpc


title: Impulse Response Function for Brownian Motion

abstract: Motivated from the central role of the mean-square displacement and its second time-derivative -- that is the velocity autocorrelation function $\left\langle v(0)v(t)\right\rangle=\frac{1}{2} \frac{\mathrm{d}^{2} \left\langle \Delta r^{2} (t)\right\rangle}{\mathrm{d}t^{2}} $ in the description of Brownian motion, we revisit the physical meaning of the first time-derivative of the mean-square displacement of Brownian particles. By employing a rheological analogue for Brownian motion, we show that the time-derivative of the mean-square displacement $\frac{\mathrm{d}\left\langle \Delta r^{2} (t) \right\rangle}{\mathrm{d}t}$ of Brownian microspheres with mass $m$ and radius $R$ immersed in any linear, isotropic viscoelastic material is identical to $\frac{N K_B T}{3 \pi R}h(t)$, where $h(t)$ is the impulse response function of a rheological network that is a parallel connection of the linear viscoelastic material with an inerter with distributed inertance $m_R=\frac{m}{6 \pi R}$. The impulse response function $h(t)=\frac{3\pi R}{N K_B T}\frac{\mathrm{d}\left\langle \Delta r^{2} (t) \right\rangle}{\mathrm{d}t}$ of the viscoelastic material-inerter parallel connection derived in this paper at the stress-strain level of the rheological analogue is essentially the response function $\chi(t)=\frac{h(t)}{6\pi R}$ of the Brownian particles expressed at the force-displacement level by Nishi \textit{et al.} (2018). By employing the viscoelastic material-inerter rheological analogue we derive the mean-square displacement and its time-derivatives of Brownian particles immersed in a viscoelastic material described with a Maxwell element connected in parallel with a dashpot which captures the high-frequency viscous behavior and we show that for Brownian motion in such fluid-like soft matter the impulse response function, $h(t)$ maintains a finite constant value in the long term.


title: A negative answer to Ulam's Problem 19 from the Scottish Book

abstract: We give a negative answer to Ulam's Problem 19 from the Scottish Book asking {\it is a solid of uniform density which will float in water in every position a sphere?} Assuming that the density of water is $1$, we show that there exists a strictly convex body of revolution $K\subset {\mathbb R^3}$ of uniform density $\frac{1}{2}$, which is not a Euclidean ball, yet floats in equilibrium in every orientation. We prove an analogous result in all dimensions $d\ge 3$.


title: Recurrent Neural Network for MoonBoard Climbing Route Classification and Generation

abstract: Classifying the difficulties of climbing routes and generating new routes are both challenging. Existing machine learning models not only fail to accurately predict a problem's difficulty, but they are also unable to generate reasonable problems. In this work, we introduced "BetaMove", a new move preprocessing pipeline we developed, in order to mimic a human climber's hand sequence. The preprocessed move sequences were then used to train both a route generator and a grade predictor. By preprocessing a MoonBoard problem into a proper move sequence, the accuracy of our grade predictor reaches near human-level performance, and our route generator produces new routes of much better quality compared to previous work. We demonstrated that with BetaMove, we are able to inject human insights into the machine learning problems, and this can be the foundations for future transfer learning on climbing style classification problems.


title: A new type of functional equations on semigroups with involutions

abstract: Let $S$ be a commutative semigroup, $K$ a quadratically closed commutative field of characteristic different from $2$, $G$ a $2$-cancellative abelian group and $H$ an abelian group uniquely divisible by $2$. The aim of this paper is to determine the general solution $f:S^2\to K$ of the d'Alembert type equation: $$ f(x+y,z+w)+f(x+\sigma(y),z+\tau(w)) =2f(x,z)f(y,w),\quad\quad (x,y,z,w\in S) $$ the general solution $f:S^2\to G$ of the Jensen type equation: $$ f(x+y,z+w)+f(x+\sigma(y),z+\tau(w)) =2f(x,z),\quad\quad (x,y,z,w\in S) $$ the general solution $f:S^2\to H$ of the quadratic type equation quation: $$ f(x+y,z+w)+f(x+\sigma(y),z+\tau(w)) =2f(x,z)+2f(y,w),\quad\quad (x,y,z,w\in S) $$ where $\sigma,\tau: S\to S$ are two involutions.


title: Couplings of the Random-Walk Metropolis algorithm

abstract: Couplings play a central role in contemporary Markov chain Monte Carlo methods and in the analysis of their convergence to stationarity. In most cases, a coupling must induce relatively fast meeting between chains to ensure good performance. In this paper we fix attention on the random walk Metropolis algorithm and examine a range of coupling design choices. We introduce proposal and acceptance step couplings based on geometric, optimal transport, and maximality considerations. We consider the theoretical properties of these choices and examine their implication for the meeting time of the chains. We conclude by extracting a few general principles and hypotheses on the design of effective couplings.


title: A single-layer based numerical method for the slender body boundary value problem

abstract: Fluid flows containing dilute or dense suspensions of thin fibers are widespread in biological and industrial processes. To describe the motion of a thin immersed fiber, or to describe the forces acting on it, it is convenient to work with one-dimensional fiber centerlines and force densities rather than two-dimensional surfaces and surface tractions. Slender body theories offer ways to model and simulate the motion of immersed fibers using only one-dimensional data. However, standard formulations can break down when the fiber surface comes close to intersecting itself or other fibers. In this paper we introduce a numerical method for a recently derived three-dimensional slender body boundary value problem that can be stated entirely in terms of a one-dimensional distribution of forces on the centerline. The method is based on a new completed single-layer potential formulation of fluid velocity which removes the nullspace associated with the unmodified single layer potential. We discretize the model and present numerical results demonstrating the good conditioning and improved performance of the method in the presence of near-intersections. To avoid the modeling and numerical choices involved with free ends, we consider closed fibers.


title: A Survey on Understanding, Visualizations, and Explanation of Deep Neural Networks

abstract: Recent advancements in machine learning and signal processing domains have resulted in an extensive surge of interest in Deep Neural Networks (DNNs) due to their unprecedented performance and high accuracy for different and challenging problems of significant engineering importance. However, when such deep learning architectures are utilized for making critical decisions such as the ones that involve human lives (e.g., in control systems and medical applications), it is of paramount importance to understand, trust, and in one word "explain" the argument behind deep models' decisions. In many applications, artificial neural networks (including DNNs) are considered as black-box systems, which do not provide sufficient clue on their internal processing actions. Although some recent efforts have been initiated to explain the behaviors and decisions of deep networks, explainable artificial intelligence (XAI) domain, which aims at reasoning about the behavior and decisions of DNNs, is still in its infancy. The aim of this paper is to provide a comprehensive overview on Understanding, Visualization, and Explanation of the internal and overall behavior of DNNs.


title: Complex Dirac-like Electronic Structure in Atomic Site Ordered Rh3In3.4Ge3.6

abstract: We report the synthesis via an indium flux method of a novel single-crystalline compound Rh3In3.4Ge3.6 that belongs to the cubic Ir3Ge7 structure type. In Rh3In3.4Ge3.6, the In and Ge atoms choose to preferentially occupy, respectively, the 12d and 16f sites of the Im-3m space group, thus creating a colored version of the Ir3Ge7 structure. Like the other compounds of the Ir3Ge7 family, Rh3In3.4Ge3.6 shows potential as a thermoelectric displaying a relatively large power factor, PF ~ 2 mW/cmK2, at a temperature T ~ 225 K albeit showing a modest figure of merit, ZT = 8 x 10-4, due to the lack of a finite band gap. These figures might improve through a use of chemical substitution strategies to achieve band gap opening. Remarkably, electronic band structure calculations reveal that this compound displays a complex Dirac-like electronic structure relatively close to the Fermi level. The electronic structure is composed of several Dirac type-I and type-II nodes, and even Dirac type-III nodes that result from the touching between a flat band and a linearly dispersing band. This rich Dirac-like electronic dispersion offers the possibility to observe Dirac type-III nodes and study their role in the physical properties of Rh3In3.4Ge3.6 and related Ir3Ge7-type materials.


title: On the Cantor and Hilbert Cube Frames and the Alexandroff-Hausdorff Theorem

abstract: The aim of this work is to give a pointfree description of the Cantor set. It can be shown that the Cantor set is homeomorphic to the $p$-adic integers $\mathbb{Z}_{p}:=\{x\in\mathbb{Q}_{p}: |x|_p\leq 1\}$ for every prime number $p$. To give a pointfree description of the Cantor set, we specify the frame of $\mathbb{Z}_{p}$ by generators and relations. We use the fact that the open balls centered at integers generate the open subsets of $\mathbb{Z}_{p}$ and thus we think of them as the basic generators; on this poset we impose some relations and then the resulting quotient is the frame of the Cantor set $\mathcal{L}(\mathbb{Z}_{p})$. We prove that $\mathcal{L}(\mathbb{Z}_{p})$ is a spatial frame whose space of points is homeomorphic to $\mathbb{Z}_{p}$. In particular, we show with pointfree arguments that $\mathcal{L}(\mathbb{Z}_{p})$ is $0$-dimensional, (completely) regular, compact, and metrizable (it admits a countably generated uniformity). Finally, we give a point-free counterpart of the Hausdorff-Alexandroff Theorem which states that \emph{every compact metric space is a continuous image of the Cantor space} (see, e.g. \cite{Alexandroff} and \cite{Hausdorff}). We prove the point-free analog: if $L$ is a compact metrizable frame, then there is an injective frame homomorphism from $L$ into $\mathcal{L}(\mathbb{Z}_{2})$.


title: Anomalies and parities for quintessential and ultra-light axions

abstract: We discuss the energy scales of the explicit breaking terms of the global symmetries \USW~ needed for the quinessential axion (QA) and the ultra-light axion (ULA). The appropriate scale of QA is about $10^{8}$ GeV.


title: Real-time exciton dynamics with time-dependent density-functional theory

abstract: Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as $-1/k^{2}$ to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate for various materials that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.


title: Data-Driven Secondary Control of Distributed Energy Resources

abstract: In this paper, we present a data-driven secondary controller for regulating to some desired values several variables of interest in a power system, namely, electrical frequency, voltage magnitudes at critical buses, and active power flows through critical lines. The power generation system is based on distributed energy resources (DERs) interfaced with either grid-forming (GFM) or grid-following (GFL) inverters. The secondary controller is based on online feedback optimization leveraging the learned sensitivities of the changes in the system frequency, voltage magnitudes at critical buses, and active power flows through critical lines to the changes in inverter active and reactive power setpoints. To learn the sensitivities accurately from data, the feedback optimization has a built-in mechanism for keeping the secondary control inputs persistently exciting without degrading its performance. The feedback optimization also utilizes the learned power-voltage characteristics of photovoltaic (PV) arrays to compute DC-link voltage setpoints so as to allow the PV arrays to track the power setpoints. To learn the power-voltage characteristics, we separately execute a data-driven approach that fits a concave polynomial to the collected power-voltage measurements by solving a sum-of-squares (SoS) optimization. We showcase the secondary controller using the modified IEEE-14 bus test system, in which conventional energy sources are replaced with inverter-interfaced DERs.


title: The role of grain-environment heterogeneity in normal grain growth: a stochastic approach

abstract: The size distribution of grains is a fundamental characteristic of polycrystalline solids. In the absence of deformation, the grain-size distribution is controlled by normal grain growth. The canonical model of normal grain growth, developed by Hillert, predicts a grain-size distribution that bears a systematic discrepancy with observed distributions. To address this, we propose a change to the Hillert model that accounts for the influence of heterogeneity in the local environment of grains. In our model, each grain evolves in response to its own local environment of neighbouring grains, rather than to the global population of grains. The local environment of each grain evolves according to an Ornstein-Uhlenbeck stochastic process. Our results are consistent with accepted grain-growth kinetics. Crucially, our model indicates that the size of relatively large grains evolves as a random walk due to the inherent variability in their local environments. This leads to a broader grain-size distribution than the Hillert model and indicates that heterogeneity has a critical influence on the evolution of microstructure.


title: Non-torsion Brauer groups in positive characteristic

abstract: Unlike the classical Brauer group of a field, the Brauer-Grothendieck group of a singular scheme need not be torsion. We show that there exist integral normal projective surfaces over a large field of positive characteristic with non-torsion Brauer group. In contrast, we demonstrate that such examples cannot exist over the algebraic closure of a finite field.


title: Optimal Intervention in Economic Networks using Influence Maximization Methods

abstract: We consider optimal intervention in the Elliott-Golub-Jackson network model \cite{jackson14} and we show that it can be transformed into an influence maximization-like form, interpreted as the reverse of a default cascade. Our analysis of the optimal intervention problem extends well-established targeting results to the economic network setting, which requires additional theoretical steps. We prove several results about optimal intervention: it is NP-hard and cannot be approximated to a constant factor in polynomial time. In turn, we show that randomizing failure thresholds leads to a version of the problem which is monotone submodular, for which existing powerful approximations in polynomial time can be applied. In addition to optimal intervention, we also show practical consequences of our analysis to other economic network problems: (1) it is computationally hard to calculate expected values in the economic network, and (2) influence maximization algorithms can enable efficient importance sampling and stress testing of large failure scenarios. We illustrate our results on a network of firms connected through input-output linkages inferred from the World Input Output Database.


title: A circular polarisation survey for radio stars with the Australian SKA Pathfinder

abstract: We present results from a circular polarisation survey for radio stars in the Rapid ASKAP Continuum Survey (RACS). RACS is a survey of the entire sky south of $\delta=+41^\circ$ being conducted with the Australian Square Kilometre Array Pathfinder telescope (ASKAP) over a 288 MHz wide band centred on 887.5 MHz. The data we analyse includes Stokes I and V polarisation products to an RMS sensitivity of 250 $\mu$Jy PSF$^{-1}$. We searched RACS for sources with fractional circular polarisation above 6 per cent, and after excluding imaging artefacts, polarisation leakage, and known pulsars we identified radio emission coincident with 33 known stars. These range from M-dwarfs through to magnetic, chemically peculiar A- and B-type stars. Some of these are well known radio stars such as YZ CMi and CU Vir, but 23 have no previous radio detections. We report the flux density and derived brightness temperature of these detections and discuss the nature of the radio emission. We also discuss the implications of our results for the population statistics of radio stars in the context of future ASKAP and Square Kilometre Array surveys.


title: Teams: Heterogeneity, Sorting, and Complementarity

abstract: How much do individuals contribute to team output? I propose an econometric framework to quantify individual contributions when only the output of their teams is observed. The identification strategy relies on following individuals who work in different teams over time. I consider two production technologies. For a production function that is additive in worker inputs, I propose a regression estimator and show how to obtain unbiased estimates of variance components that measure the contributions of heterogeneity and sorting. To estimate nonlinear models with complementarity, I propose a mixture approach under the assumption that individual types are discrete, and rely on a mean-field variational approximation for estimation. To illustrate the methods, I estimate the impact of economists on their research output, and the contributions of inventors to the quality of their patents.


title: Tracking fast and slow changes in synaptic weights from simultaneously observed pre- and postsynaptic spiking

abstract: Synapses change on multiple timescales, ranging from milliseconds to minutes, due to a combination of both short- and long-term plasticity. Here we develop an extension of the common Generalized Linear Model to infer both short- and long-term changes in the coupling between a pre- and post-synaptic neuron based on observed spiking activity. We model short-term synaptic plasticity using additive effects that depend on the presynaptic spike timing, and we model long-term changes in both synaptic weight and baseline firing rate using point process adaptive smoothing. Using simulations, we first show that this model can accurately recover time-varying synaptic weights 1) for both depressing and facilitating synapses, 2) with a variety of long-term changes (including realistic changes, such as due to STDP), 3) with a range of pre- and post-synaptic firing rates, and 4) for both excitatory and inhibitory synapses. We then apply our model to two experimentally recorded putative synaptic connections. We find that simultaneously tracking fast changes in synaptic weights, slow changes in synaptic weights, and unexplained variations in baseline firing is essential. Omitting any one of these factors can lead to spurious inferences for the others. Altogether, this model provides a flexible framework for tracking short- and long-term variation in spike transmission.


title: Proceedings 11th International Workshop on Computing with Terms and Graphs

abstract: Graphs, and graph transformation systems, are used in many areas within Computer Science: to represent data structures and algorithms, to define computation models, as a general modelling tool to study complex systems, etc. Research in term and graph rewriting ranges from theoretical questions to practical implementation issues. Relevant research areas include: the modelling of first- and higher-order term rewriting by graph rewriting, graphical frameworks such as interaction nets and sharing graphs (optimal reduction), rewrite calculi for the analysis of functional programs, graph reduction implementations of programming languages, graphical calculi modelling concurrent and mobile computations, object-oriented systems, graphs as a model of biological or chemical systems, and automated reasoning and symbolic computation systems working on shared structures. The aim of the TERMGRAPH workshop is to bring together researchers working in these different domains and to foster their interaction, to provide a forum for presenting new ideas and work in progress, and to enable newcomers to learn about current activities in this area.