title: Effects of simultaneous real-time fMRI and EEG neurofeedback in major depressive disorder evaluated with brain electromagnetic tomography

abstract: Recently, we reported an emotion self-regulation study (Zotev et al., 2020), in which patients with major depressive disorder (MDD) used simultaneous real-time fMRI and EEG neurofeedback (rtfMRI-EEG-nf) to upregulate two fMRI and two EEG activity measures, relevant to MDD. The target measures included fMRI activities of the left amygdala and left rostral anterior cingulate cortex, and frontal EEG asymmetries in the alpha band (FAA) and high-beta band (FBA). Here we apply the exact low resolution brain electromagnetic tomography (eLORETA) to investigate EEG source activities during the rtfMRI-EEG-nf procedure. The exploratory analyses reveal significant changes in hemispheric lateralities of upper alpha and high-beta current source densities in the prefrontal regions, consistent with upregulation of the FAA and FBA during the rtfMRI-EEG-nf task. Similar laterality changes are observed for current source densities in the amygdala. Prefrontal upper alpha current density changes show significant negative correlations with anhedonia severity. Changes in prefrontal high-beta current density are consistent with reduction in comorbid anxiety. Comparisons with results of previous LORETA studies suggest that the rtfMRI-EEG-nf training is beneficial to MDD patients, and may have the ability to correct functional deficiencies associated with anhedonia and comorbid anxiety in MDD.


abstract: We present a novel technique for sparse principal component analysis. This method, named Eigenvectors from Eigenvalues Sparse Principal Component Analysis (EESPCA), is based on the formula for computing squared eigenvector loadings of a Hermitian matrix from the eigenvalues of the full matrix and associated sub-matrices. We explore two versions of the EESPCA method: a version that uses a fixed threshold for inducing sparsity and a version that selects the threshold via cross-validation. Relative to the state-of-the-art sparse PCA methods of Witten et al., Yuan & Zhang and Tan et al., the fixed threshold EESPCA technique offers an order-of-magnitude improvement in computational speed, does not require estimation of tuning parameters via cross-validation, and can more accurately identify true zero principal component loadings across a range of data matrix sizes and covariance structures. Importantly, the EESPCA method achieves these benefits while maintaining out-of-sample reconstruction error and PC estimation error close to the lowest error generated by all evaluated approaches. EESPCA is a practical and effective technique for sparse PCA with particular relevance to computationally demanding statistical problems such as the analysis of high-dimensional data sets or application of statistical techniques like resampling that involve the repeated calculation of sparse PCs.


title: Asynchronous Distributed Averaging: A Switched System Framework for Average Error Analysis

abstract: This paper investigates an expected average error for distributed averaging problems under asynchronous updates. The asynchronism in this context implies no existence of a global clock as well as random characteristics in communication uncertainty such as communication delays and packet drops. Although some previous works contributed to the design of average consensus protocols to guarantee the convergence to an exact average, these methods may increase computational burdens due to extra works. Sometimes it is thus beneficial to make each agent exchange information asynchronously without modifying the algorithm, which causes randomness in the average value as a trade-off. In this study, an expected average error is analyzed based on the switched system framework, to estimate an upper bound of the asynchronous average compared to the exact one in the expectation sense. Numerical examples are provided to validate the proposed results.


title: Would You Like to Hear the News? Investigating Voice-BasedSuggestions for Conversational News Recommendation

abstract: One of the key benefits of voice-based personal assistants is the potential to proactively recommend relevant and interesting information. One of the most valuable sources of such information is the News. However, in order for the user to hear the news that is useful and relevant to them, it must be recommended in an interesting and informative way. However, to the best of our knowledge, how to present a news item for a voice-based recommendation remains an open question. In this paper, we empirically compare different ways of recommending news, or specific news items, in a voice-based conversational setting. Specifically, we study the user engagement and satisfaction with five different variants of presenting news recommendations: (1) a generic news briefing; (2) news about a specific entity relevant to the current conversation; (3) news about an entity from a past conversation; (4) news on a trending news topic; and (5) the default - a suggestion to talk about news in general. Our results show that entity-based news recommendations exhibit 29% higher acceptance compared to briefing recommendations, and almost 100% higher acceptance compared to recommending generic or trending news. Our investigation into the presentation of news recommendations and the resulting insights could make voice assistants more informative and engaging.


title: Incorporating Physical Knowledge into Machine Learning for Planetary Space Physics

abstract: Recent improvements in data collection volume from planetary and space physics missions have allowed the application of novel data science techniques. The Cassini mission for example collected over 600 gigabytes of scientific data from 2004 to 2017. This represents a surge of data on the Saturn system. Machine learning can help scientists work with data on this larger scale. Unlike many applications of machine learning, a primary use in planetary space physics applications is to infer behavior about the system itself. This raises three concerns: first, the performance of the machine learning model, second, the need for interpretable applications to answer scientific questions, and third, how characteristics of spacecraft data change these applications. In comparison to these concerns, uses of black box or un-interpretable machine learning methods tend toward evaluations of performance only either ignoring the underlying physical process or, less often, providing misleading explanations for it. We build off a previous effort applying a semi-supervised physics-based classification of plasma instabilities in Saturn's magnetosphere. We then use this previous effort in comparison to other machine learning classifiers with varying data size access, and physical information access. We show that incorporating knowledge of these orbiting spacecraft data characteristics improves the performance and interpretability of machine learning methods, which is essential for deriving scientific meaning. Building on these findings, we present a framework on incorporating physics knowledge into machine learning problems targeting semi-supervised classification for space physics data in planetary environments. These findings present a path forward for incorporating physical knowledge into space physics and planetary mission data analyses for scientific discovery.


title: Water Assisted Proton Transport in Confined Nanochannels

abstract: Hydrated excess protons under hydrophobic confinement are a critical component of charge transport behavior and reactivity in nanoporous materials and biomolecular systems. Herein excess proton confinement effects are computationally investigated for sub-2 nm hydrophobic nanopores by varying the diameters (d = 0.81, 0.95, 1.09, 1.22, 1.36, 1.63, and 1.90 nm), lengths (l ~3 and ~5 nm), curvature, and chirality of cylindrical carbon nanotube (CNT) nanopores. CNTs with a combination of different diameter segments are also explored. The spatial distribution of water molecules under confinement are diameter-dependent; however, proton solvation and transport is consistently found to occur in the water layer adjacent to the pore wall, showing an "amphiphilic" character of the hydrated excess proton hydronium-like structure. The proton transport free energy barrier also decreases significantly as the nanopore diameter increases and proton transport becomes almost barrierless in the d > 1 nm nanopores. Among the nanopores studied, the Zundel cation (${H_{5}O_{2}}^{+}$) is populated only in the d = 0.95 nm CNT (7,7) nanopore. The presence of the hydrated excess proton and $K^{+}$ inside the CNT (7,7) nanopore induces a water density increase by 40 and 20%, respectively. The $K^{+}$ transport through CNT nanopores is also consistently higher in free energy barrier than proton transport. Interestingly, the evolution of excess protonic charge defect distribution reveals a "frozen" single water wire configuration in the d = 0.81 nm CNT (6,6) nanopore (or segment), through which hydrated excess protons can only shuttle via the Grotthuss mechanism. Vehicular diffusion becomes relevant to proton transport in the "flat" free energy regions and in the wider nanopores, where protons do not primarily shuttle in the axial direction.


title: A Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

abstract: Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation-based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation, will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. It is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.


title: Effect of axially varying sandpaper roughness on bubbly drag reduction in Taylor-Couette turbulence

abstract: We experimentally investigate the influence of alternating rough and smooth walls on bubbly drag reduction (DR). We apply rough sandpaper bands of width $s$ between $48.4\,mm$ and $148.5\,mm$, and roughness height $k = 695\,{\mu}m$, around the smooth inner cylinder (IC) of the Twente Turbulent Taylor-Couette facility. Between sandpaper bands, the IC is left uncovered over similar width $s$, resulting in alternating rough and smooth bands, a constant pattern in axial direction. We measure the DR in water that originates from introducing air bubbles to the fluid at (shear) Reynolds numbers $\textit{Re}_s$ ranging from $0.5 \times 10^6$ to $1.8 \times 10^6$. Results are compared to bubbly DR measurements with a completely smooth IC and an IC that is completely covered with sandpaper of the same roughness $k$. The outer cylinder is left smooth for all variations. Results are also compared to bubbly DR measurements where a smooth outer cylinder is rotating in opposite direction to the smooth IC. This counter rotation induces secondary flow structures that are very similar to those observed when the IC is composed of alternating rough and smooth bands. For the measurements with roughness, the bubbly DR is found to initially increase more strongly with $\textit{Re}_s$, before levelling off to reach a value that no longer depends on $\textit{Re}_s$. This is attributed to a more even axial distribution of the air bubbles, resulting from the increased turbulence intensity of the flow compared to flow over a completely smooth wall at the same $\textit{Re}_s$. The air bubbles are seen to accumulate at the rough wall sections in the flow. Here, locally, the drag is largest and so the drag reducing effect of the bubbles is felt strongest. Therefore, a larger maximum value of bubbly DR is found for the alternating rough and smooth walls compared to the completely rough wall.


title: Boshernitzan's condition, factor complexity, and an application

abstract: Boshernitzan found a decay condition on the measure of cylinder sets that implies unique ergodicity for minimal subshifts. Interest in the properties of subshifts satisfying this condition has grown recently, due to a connection with the study of discrete Schr\"odinger operators. Of particular interest is the question of how restrictive Boshernitzan's condition is. While it implies zero topological entropy, our main theorem shows how to construct minimal subshifts satisfying the condition whose factor complexity grows faster than any pre-assigned subexponential rate. As an application, via a theorem of Damanik and Lenz, we show that there is no subexponentially growing sequence for which the spectra of all discrete Schr\"odinger operators associated with subshifts whose complexity grows faster than the given sequence, have only finitely many gaps.


title: Nonlinear nonlocal Douglas identity

abstract: We give Hardy-Stein and Douglas identities for nonlinear nonlocal Sobolev-Bregman integral forms with unimodal L\'evy measures. We prove that the corresponding Poisson integral defines an extension operator for the Sobolev-Bregman spaces. We also show that the Poisson integrals are quasiminimizers of the Sobolev-Bregman forms.


title: Hierarchical Clustering: a 0.585 Revenue Approximation

abstract: Hierarchical Clustering trees have been widely accepted as a useful form of clustering data, resulting in a prevalence of adopting fields including phylogenetics, image analysis, bioinformatics and more. Recently, Dasgupta (STOC 16') initiated the analysis of these types of algorithms through the lenses of approximation. Later, the dual problem was considered by Moseley and Wang (NIPS 17') dubbing it the Revenue goal function. In this problem, given a nonnegative weight $w_{ij}$ for each pair $i,j \in [n]=\{1,2, \ldots ,n\}$, the objective is to find a tree $T$ whose set of leaves is $[n]$ that maximizes the function $\sum_{i<j \in [n]} w_{ij} (n -|T_{ij}|)$, where $|T_{ij}|$ is the number of leaves in the subtree rooted at the least common ancestor of $i$ and $j$. In our work we consider the revenue goal function and prove the following results. First, we prove the existence of a bisection (i.e., a tree of depth 2 in which the root has two children, each being a parent of $n/2$ leaves) which approximates the general optimal tree solution up to a factor of $\frac{1}{2}$ (which is tight). Second, we apply this result in order to prove a $\frac{2}{3}p$ approximation for the general revenue problem, where $p$ is defined as the approximation ratio of the Max-Uncut Bisection problem. Since $p$ is known to be at least 0.8776 (Wu et al., 2015, Austrin et al., 2016), we get a 0.585 approximation algorithm for the revenue problem. This improves a sequence of earlier results which culminated in an 0.4246-approximation guarantee (Ahmadian et al., 2019).


title: Electron and muon $g-2$ anomalies in general flavour conserving two Higgs doublets models

abstract: In general two Higgs doublet models (2HDMs) without scalar flavour changing neutral couplings (SFCNC) in the lepton sector, the electron, muon and tau interactions can be decoupled in a robust framework, stable under renormalization group evolution. In this framework, the breaking of lepton flavour universality (LFU) goes beyond the mass proportionality, opening the possibility to accommodate a different behaviour among charged leptons. We analyze the electron and muon $(g-2)$ anomalies in the context of these general flavour conserving models in the leptonic sector (g$\ell$FC). We consider two different models, I-g$\ell$FC and II-g$\ell$FC, in which the quark Yukawa couplings coincide, respectively, with the ones in type I and in type II 2HDMs. We find two types of solutions that fully reproduce both $(g-2)$ anomalies, and which are compatible with experimental constraints from LEP and LHC, from LFU, from flavour and electroweak physics, and with theoretical constraints in the scalar sector. In the first type of solution, all the new scalars have masses in the 1--2.5 TeV range, the vacuum expectation values (vevs) of both doublets are quite similar in magnitude, and both anomalies are dominated by two loop Barr-Zee contributions. This solution appears in both models. In a second type of solution, one loop contributions are dominant in the muon anomaly, all new scalars have masses below 1 TeV, and the ratio of vevs is in the range 10--100. The second neutral scalar $H$ is the lighter among the new scalars, with a mass in the 210--390 GeV range while the pseudoscalar $A$ is the heavier, with a mass in the range 400--900 GeV. The new charged scalar $H^\pm$ is almost degenerate either with the scalar or with the pseudoscalar. This second type of solution only appears in the I-g$\ell$FC model. Both solutions require the soft breaking of the $\mathbb{Z}_{2}$ symmetry of the Higgs potential.


title: Analysis of the Schwarz domain decomposition method for the conductor-like screening continuum model

abstract: We study the Schwarz overlapping domain decomposition method applied to the Poisson problem on a special family of domains, which by construction consist of a union of a large number of fixed-size subdomains. These domains are motivated by applications in computational chemistry where the subdomains consist of van der Waals balls. As is usual in the theory of domain decomposition methods, the rate of convergence of the Schwarz method is related to a stable subspace decomposition. We derive such a stable decomposition for this family of domains and analyze how the stability "constant" depends on relevant geometric properties of the domain. For this, we introduce new descriptors that are used to formalize the geometry for the family of domains. We show how, for an increasing number of subdomains, the rate of convergence of the Schwarz method depends on specific local geometry descriptors and on one global geometry descriptor. The analysis also naturally provides lower bounds in terms of the descriptors for the smallest eigenvalue of the Laplace eigenvalue problem for this family of domains.


title: An Alternative Metric for Detecting Anomalous Ship Behavior Using a Variation of the DBSCAN Clustering Algorithm

abstract: There is a growing need to quickly and accurately identify anomalous behavior in ships. This paper applies a variation of the Density Based Spatial Clustering Among Noise (DBSCAN) algorithm to identify such anomalous behavior given a ship's Automatic Identification System (AIS) data. This variation of the DBSCAN algorithm has been previously introduced in the literature, and in this study, we elucidate and explore the mathematical details of this algorithm and introduce an alternative anomaly metric which is more statistically informative than the one previously suggested.


title: An Extremal Problem Motivated by Triangle-Free Strongly Regular Graphs

abstract: We introduce the following combinatorial problem. Let $G$ be a triangle-free regular graph with edge density $\rho$. What is the minimum value $a(\rho)$ for which there always exist two non-adjacent vertices such that the density of their common neighborhood is $\leq a(\rho)$? We prove a variety of upper bounds on the function $a(\rho)$ that are tight for the values $\rho=2/5,\ 5/16,\ 3/10,\ 11/50$, with $C_5$, Clebsch, Petersen and Higman-Sims being respective extremal configurations. Our proofs are entirely combinatorial and are largely based on counting densities in the style of flag algebras. For small values of $\rho$, our bound attaches a combinatorial meaning to Krein conditions that might be interesting in its own right. We also prove that for any $\epsilon>0$ there are only finitely many values of $\rho$ with $a(\rho)\geq\epsilon$ but this finiteness result is somewhat purely existential (the bound is double exponential in $1/\epsilon$).


title: Nurse is Closer to Woman than Surgeon? Mitigating Gender-Biased Proximities in Word Embeddings

abstract: Word embeddings are the standard model for semantic and syntactic representations of words. Unfortunately, these models have been shown to exhibit undesirable word associations resulting from gender, racial, and religious biases. Existing post-processing methods for debiasing word embeddings are unable to mitigate gender bias hidden in the spatial arrangement of word vectors. In this paper, we propose RAN-Debias, a novel gender debiasing methodology which not only eliminates the bias present in a word vector but also alters the spatial distribution of its neighbouring vectors, achieving a bias-free setting while maintaining minimal semantic offset. We also propose a new bias evaluation metric - Gender-based Illicit Proximity Estimate (GIPE), which measures the extent of undue proximity in word vectors resulting from the presence of gender-based predilections. Experiments based on a suite of evaluation metrics show that RAN-Debias significantly outperforms the state-of-the-art in reducing proximity bias (GIPE) by at least 42.02%. It also reduces direct bias, adding minimal semantic disturbance, and achieves the best performance in a downstream application task (coreference resolution).


title: A New Chaos and Permutation Based Algorithm for Image and Video Encryption

abstract: Images and video sequences carry large volumes of highly correlated and redundant data. Applications like military and telecommunication require encryption methods to protect the data from unwanted access. This requirement in most cases needs to be realized in real-time. In this paper, we propose a fast new Fiestal-structured approach for image and video encryption based on a chaotic random sequence generator and a Permutation-Inverse Permutation (PIP) pixel transform. This approach utilizes mathematical functions and transforms with low complexity. The algorithm at the same time, ensures no drastic pay off in terms of encryption quality. This renders the algorithm with promising scope for real time applications and easy hardware implementation. MATLAB simulation of the algorithm establishes its high quality of encryption in terms of elevated entropy values and negligible correlation of the encrypted data with the original. Simulation results also show high sensitivity to slight variation in keys ensuring high security.


title: The Art of Modeling Stellar Mergers and the Case of the B[e] Supergiant R4 in the Small Magellanic Cloud

abstract: Most massive stars exchange mass with a companion, leading to evolution which is altered drastically from that expected of stars in isolation. Such systems are the result of unusual binary evolution pathways and, as such, may be used to place stringent constraints on the physics of these interactions. We use the R4 system's B[e] supergiant, which has been postulated to be the product of a binary stellar merger, to guide our understanding of such outcomes by comparing observations of R4 to the results of simulations of mergers performed with the 3d hydrodynamics code FLASH. Our approach tailors the simulation initial conditions to the observed properties of R4 and implements realistic stellar profiles generated by the 1d stellar evolution code MESA onto the 3d grid, resolving the merger inspiral to within $0.02\, R_{\odot}$. We then map the merger remnant into MESA to track its evolution on the HR diagram over a period of $10^4$ years. This generates models for a B[e] supergiant with stellar properties, age, and nebula structure in qualitative agreement with that of the R4 system. Our calculations provide concrete evidence to support the idea that R4 was originally a member of a triple system in which the inner binary merged after its most massive member evolved off the main sequence, producing a new object that is of similar mass yet significantly more luminous than the A supergiant companion. The potential applications of the code framework presented in this paper are wide ranging and can be used to generate models of a variety of merger stellar remnants.


title: Vanishing Flats: A Combinatorial Viewpoint on the Planarity of Functions and Their Application

abstract: For a function $f$ from $\mathbb{F}_2^n$ to $\mathbb{F}_2^n$, the planarity of $f$ is usually measured by its differential uniformity and differential spectrum. In this paper, we propose the concept of vanishing flats, which supplies a combinatorial viewpoint on the planarity. First, the number of vanishing flats of $f$ can be regarded as a measure of the distance between $f$ and the set of almost perfect nonlinear functions. In some cases, the number of vanishing flats serves as an "intermediate" concept between differential uniformity and differential spectrum, which contains more information than differential uniformity, however less than the differential spectrum. Secondly, the set of vanishing flats forms a combinatorial configuration called partial quadruple system, since it convey detailed structural information about $f$. We initiate this study by considering the number of vanishing flats and the partial quadruple systems associated with monomials and Dembowski-Ostrom polynomials. In addition, we present an application of vanishing flats to the partition of a vector space into disjoint equidimensional affine spaces. We conclude the paper with several further questions and challenges.


title: On alternative approximating distributions in the multivariate version of Kolmogorov's second uniform limit theorem

abstract: The aim of the present work is to show that recent results of the authors on the approximation of distributions of sums of independent summands by the infinitely divisible laws on convex polyhedra can be shown via an alternative class of approximating infinitely divisible distributions. We will also generalize the results to the infinite-dimensional case.


title: Ear2Face: Deep Biometric Modality Mapping

abstract: In this paper, we explore the correlation between different visual biometric modalities. For this purpose, we present an end-to-end deep neural network model that learns a mapping between the biometric modalities. Namely, our goal is to generate a frontal face image of a subject given his/her ear image as the input. We formulated the problem as a paired image-to-image translation task and collected datasets of ear and face image pairs from the Multi-PIE and FERET datasets to train our GAN-based models. We employed feature reconstruction and style reconstruction losses in addition to adversarial and pixel losses. We evaluated the proposed method both in terms of reconstruction quality and in terms of person identification accuracy. To assess the generalization capability of the learned mapping models, we also run cross-dataset experiments. That is, we trained the model on the FERET dataset and tested it on the Multi-PIE dataset and vice versa. We have achieved very promising results, especially on the FERET dataset, generating visually appealing face images from ear image inputs. Moreover, we attained a very high cross-modality person identification performance, for example, reaching 90.9% Rank-10 identification accuracy on the FERET dataset.


title: Designing Differentially Private Estimators in High Dimensions

abstract: We study differentially private mean estimation in a high-dimensional setting. Existing differential privacy techniques applied to large dimensions lead to computationally intractable problems or estimators with excessive privacy loss. Recent work in high-dimensional robust statistics has identified computationally tractable mean estimation algorithms with asymptotic dimension-independent error guarantees. We incorporate these results to develop a strict bound on the global sensitivity of the robust mean estimator. This yields a computationally tractable algorithm for differentially private mean estimation in high dimensions with dimension-independent privacy loss. Finally, we show on synthetic data that our algorithm significantly outperforms classic differential privacy methods, overcoming barriers to high-dimensional differential privacy.


title: Continual Learning of Predictive Models in Video Sequences via Variational Autoencoders

abstract: This paper proposes a method for performing continual learning of predictive models that facilitate the inference of future frames in video sequences. For a first given experience, an initial Variational Autoencoder, together with a set of fully connected neural networks are utilized to respectively learn the appearance of video frames and their dynamics at the latent space level. By employing an adapted Markov Jump Particle Filter, the proposed method recognizes new situations and integrates them as predictive models avoiding catastrophic forgetting of previously learned tasks. For evaluating the proposed method, this article uses video sequences from a vehicle that performs different tasks in a controlled environment.


title: Broadband Terahertz Wave Emission from Liquid Metal

abstract: Metals have been studied as terahertz sources for decades. Recent research has shown the potential of metals in generating extremely high THz pulse energy excited by intense laser pulses. To avoid the metal surface debris caused by laser pulses, here, we report the results of the broadband terahertz wave emission from a flowing liquid metal line excited by sub-picosecond laser pulses. The THz signal emitted from the liquid gallium line shows stronger field with broader bandwidth comparing with the signal from water under the identical optical excitation conditions. Our preliminary study suggests that the liquid metals have the potential to serve as efficient and powerful THz sources for the intense lasers with a high repetition rate.


title: Study of muonium emission from laser-ablated silica aerogel

abstract: The emission of muonium ($\mu^+e^-$) atoms into vacuum from silica aerogel with laser ablation on its surface was studied with various ablation structures at room temperature using the subsurface muon beams at TRIUMF and Japan Proton Accelerator Research Complex (J-PARC). Laser ablation was applied to produce holes or grooves with typical dimensions of a few hundred $\mu$m to a few mm, except for some extreme conditions. The measured emission rate tends to be higher for larger fractions of ablation opening and for shallower depths. More than a few ablation structures reach the emission rates similar to the highest achieved in the past measurements. The emission rate is found to be stable at least for a couple of days. Measurements of spin precession amplitudes for the produced muonium atoms and remaining muons in a magnetic field determine a muonium formation fraction of $(65.5 \pm 1.8)$%. The precession of the polarized muonium atoms is also observed clearly in vacuum. A projection of the emission rates measured at TRIUMF to the corresponding rates at J-PARC is demonstrated taking the different beam condition into account reasonably.


title: Understanding the MiniBooNE and the muon and electron $g-2$ anomalies with a light $Z'$ and a second Higgs doublet

abstract: Two of the most widely studied extensions of the Standard Model (SM) are $a)$ the addition of a new $U(1)$ symmetry to its existing gauge groups, and $b)$ the expansion of its scalar sector to incorporate a second Higgs doublet. We show that when combined, they allow us to understand the electron-like event excess seen in the MiniBooNE (MB) experiment as well as account for the observed anomalous values of the muon magnetic moment. A light $Z'$ associated with an additional $U(1)$ coupled to baryons and to the dark sector, with flavor non-universal couplings to leptons, in conjunction with a second Higgs doublet is capable of explaining the MB excess. The $Z'$ obtains its mass from a dark singlet scalar, which mixes with the two Higgs doublets. Choosing benchmark parameter values, we show that $U(1)_{B-3L_\tau}$, which is anomaly-free, and $U(1)_B$, both provide (phenomenologically) equally good solutions to the excess. We also point out the other (anomaly-free) $U(1)$ choices that may be possible upon fuller exploration of the parameter space. We obtain very good matches to the energy and angular distributions for neutrinos and anti-neutrinos in MB. The extended Higgs sector has two light CP-even scalars, $h'$ and $H$, and their masses and couplings are such that in principle, both contribute to help explain the MB excess as well as the present observed values of the muon and electron $g-2$. We discuss the constraints on our model as well as future tests. Our work underlines the role that light scalars may play in understanding present-day low-energy anomalies. It also points to the possible existence of portals to the dark sector, i.e., a light gauge boson field $(Z')$ and a dark neutrino which mixes with the active neutrinos, as well as a dark sector light scalar which mixes with the extended Higgs sector.


title: Asymptotically Scale-invariant Multi-resolution Quantization

abstract: A multi-resolution quantizer is a sequence of quantizers where the output of a coarser quantizer can be deduced from the output of a finer quantizer. In this paper, we propose an asymptotically scale-invariant multi-resolution quantizer, which performs uniformly across any choice of average quantization step, when the length of the range of input numbers is large. Scale invariance is especially useful in worst case or adversarial settings, ensuring that the performance of the quantizer would not be affected greatly by small changes of storage or error requirements. We also show that the proposed quantizer achieves a tradeoff between rate and error that is arbitrarily close to the optimum.


title: Transient Work Function Gating: A New Photoemission Regime

abstract: We present the theoretical basis for a new photoemission regime, transient work function gating (TWFG), that temporally and energetically gates photoemission and produces near-threshold photoelectrons with thermally limited emittance, percent-level quantum efficiency, and control over temporal coherence. The technique consists of actively gating the work function of a generalized photocathode using non-ionizing long-wavelength optical field to produce an adiabatic modulation of the carrier density at their surface. We examine TWFG as a means to circumvent the long-standing trade-off between low emittance and high quantum efficiency, untethered to particle source or photocathode specifics. TWFG promises new opportunities in photoemission physics for next generation electron and accelerator-based x-ray photon sources.

categories: physics.app-ph physics.acc-ph physics.optics quant-ph


title: Guided matter wave inertial sensing in a miniature physics package

abstract: We describe an ultra-compact ($\sim 10$ cm$^3$ physics package) inertial sensor based on atomic matter waves that are guided within an optical lattice during almost the entire interferometer cycle. We demonstrate large momentum transfer (LMT) of up to 8$\hbar k$ photon momentum with a combination of Bragg pulses and Bloch oscillations with scalability to larger numbers of photons. Between momentum transfer steps, we maintain the atoms in a co-moving optical lattice waveguide so that the atoms are in free space only during the Bragg pulses. Our guided matter wave approach paves the way for atomic inertial sensing in dynamic environments in which untrapped atoms would otherwise quickly collide with the walls of a miniature chamber.


title: Learning Active Task-Oriented Exploration Policies for Bridging the Sim-to-Real Gap

abstract: Training robotic policies in simulation suffers from the sim-to-real gap, as simulated dynamics can be different from real-world dynamics. Past works tackled this problem through domain randomization and online system-identification. The former is sensitive to the manually-specified training distribution of dynamics parameters and can result in behaviors that are overly conservative. The latter requires learning policies that concurrently perform the task and generate useful trajectories for system identification. In this work, we propose and analyze a framework for learning exploration policies that explicitly perform task-oriented exploration actions to identify task-relevant system parameters. These parameters are then used by model-based trajectory optimization algorithms to perform the task in the real world. We instantiate the framework in simulation with the Linear Quadratic Regulator as well as in the real world with pouring and object dragging tasks. Experiments show that task-oriented exploration helps model-based policies adapt to systems with initially unknown parameters, and it leads to better task performance than task-agnostic exploration.


title: Actions of groups of foliated homeomorphisms on spaces of leaves

abstract: Let $\Delta$ be a foliation on a topological manifold $X$, $Y$ be the space of leaves, and $p: X \to Y$ be the natural projection. Endow $Y$ with the factor topology with respect to $p$. Then the group $\mathcal{H}(X, \Delta)$ of foliated (i.e. mapping leaves onto leaves) homeomorphisms of $X$ naturally acts on the space of leaves $Y$, which gives a homomorphism $\psi: \mathcal{H}(X, \Delta) \to \mathcal{H}(Y)$. We present sufficient conditions when $\psi$ is continuous with respect to the corresponding compact open topologies. In fact similar results hold not only for foliations but for a more general class of partitions $\Delta$ of locally compact Hausdorff spaces $X$.


title: Unidimensional and Multidimensional Methods for Recurrence Quantification Analysis with crqa

abstract: Recurrence quantification analysis is a widely used method for characterizing patterns in time series. This article presents a comprehensive survey for conducting a wide range of recurrence-based analyses to quantify the dynamical structure of single and multivariate time series, and to capture coupling properties underlying leader-follower relationships. The basics of recurrence quantification analysis (RQA) and all its variants are formally introduced step-by-step from the simplest auto-recurrence to the most advanced multivariate case. Importantly, we show how such RQA methods can be deployed under a single computational framework in R using a substantially renewed version our crqa 2.0 package. This package includes implementations of several recent advances in recurrence-based analysis, among them applications to multivariate data, and improved entropy calculations for categorical data. We show concrete applications of our package to example data, together with a detailed description of its functions and some guidelines on their usage.


title: Tightness of radially-symmetric solutions to 2D aggregation-diffusion equations with weak interaction forces

abstract: We prove the tightness of radially-symmetric solutions to 2D aggregation-diffusion equations, where the pairwise attraction force is possibly degenerate at large distance. We first reduce the problem into the finiteness of a time integral in the density on a bounded region, by introducing a new assumption on the interaction potential called the essentially radially contractive (ERC) property. Then we prove this finiteness by using the 2-Wasserstein gradient flow structure, combining with the continuous Steiner symmetrization curves and the local clustering curve. This is the first tightness result on the 2D aggregation-diffusion equations for a general class of weakly confining potentials, i.e., those $W(r)$ with $\lim_{r\rightarrow\infty}W(r)<\infty$, and serves as an important step towards the study of equilibration.


title: Interaction between kinks and antikinks with double long-range tails

abstract: We explore a class of $\phi^{4n}$ models with kink and antikink solutions that have long-range tails on both sides, specializing to the cases with $n=2$ and $n=3$. A recently developed method of an accelerating kink ansatz is used to estimate the force between the kink and the antikink. We use state-of-the-art numerical methods to initialize the system in a kink-antikink configuration with a finite initial velocity and to evolve the system according to the equations of motion. Among these methods, we propose a computationally efficient way to initialize the velocity field of the system. Interestingly, we discover that, for this class of models, $\phi^{4n}$ with $n>1$, the kink-antikink annihilation behaves differently from the archetypal $\phi^4$ model or even the kinks with one long-range tail because there is neither long-lived bion formation nor resonance windows and the critical velocity is ultrarelativistic.


title: Workflow-as-a-Service Cloud Platform and Deployment of Bioinformatics Workflow Applications

abstract: Workflow management systems (WMS) support the composition and deployment of workflow-oriented applications in distributed computing environments. They hide the complexity of managing large-scale applications, which includes the controlling data pipelining between tasks, ensuring the application's execution, and orchestrating the distributed computational resources to get a reasonable processing time. With the increasing trends of scientific workflow adoption, the demand to deploy them using a third-party service begins to increase. Workflow-as-a-service (WaaS) is a term representing the platform that serves the users who require to deploy their workflow applications on third-party cloud-managed services. This concept drives the existing WMS technology to evolve towards the development of the WaaS cloud platform. Based on this requirement, we extend CloudBus WMS functionality to handle the workload of multiple workflows and develop the WaaS cloud platform prototype. We implemented the Elastic Budget-constrained resource Provisioning and Scheduling algorithm for Multiple workflows (EBPSM) algorithm that is capable of scheduling multiple workflows and evaluated the platform using two bioinformatics workflows. Our experimental results show that the platform is capable of efficiently handling multiple workflows execution and gaining its purpose to minimize the makespan while meeting the budget.


title: Nucleus Decomposition in Probabilistic Graphs: Hardness and Algorithms

abstract: Finding dense components in graphs is of great importance in analyzing the structure of networks. Popular and computationally feasible frameworks for discovering dense subgraphs are core and truss decompositions. Recently, Sariyuce et al. introduced nucleus decomposition, a generalization which uses higher-order structures and can reveal interesting subgraphs that can be missed by core and truss decompositions. In this paper, we present nucleus decomposition in probabilistic graphs. We study the most interesting case of nucleus decomposition, k-(3,4)-nucleus, which asks for maximal subgraphs where each triangle is contained in k 4-cliques. The major questions we address are: How to define meaningfully nucleus decomposition in probabilistic graphs? How hard is computing nucleus decomposition in probabilistic graphs? Can we devise efficient algorithms for exact or approximate nucleus decomposition in large graphs? We present three natural definitions of nucleus decomposition in probabilistic graphs: local, global, and weakly-global. We show that the local version is in PTIME, whereas global and weakly-global are #P-hard and NP-hard, respectively. We present an efficient and exact dynamic programming approach for the local case and furthermore, present statistical approximations that can scale to large datasets without much loss of accuracy. For global and weakly-global decompositions, we complement our intractability results by proposing efficient algorithms that give approximate solutions based on search space pruning and Monte-Carlo sampling. Our extensive experimental results show the scalability and efficiency of our algorithms. Compared to probabilistic core and truss decompositions, nucleus decomposition significantly outperforms in terms of density and clustering metrics.


title: NewtonianVAE: Proportional Control and Goal Identification from Pixels via Physical Latent Spaces

abstract: Learning low-dimensional latent state space dynamics models has been a powerful paradigm for enabling vision-based planning and learning for control. We introduce a latent dynamics learning framework that is uniquely designed to induce proportional controlability in the latent space, thus enabling the use of much simpler controllers than prior work. We show that our learned dynamics model enables proportional control from pixels, dramatically simplifies and accelerates behavioural cloning of vision-based controllers, and provides interpretable goal discovery when applied to imitation learning of switching controllers from demonstration.


title: Mass loss from massive globular clusters in tidal fields

abstract: Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for tens of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles.


title: Automatic Counting of Restricted Dyck Paths via (Numeric and Symbolic) Dynamic Programming

abstract: Dyck paths are one of the most important objects in enumerative combinatorics, and there are many papers devoted to counting selected families of Dyck paths. Here we present two approaches for the automatic counting of many such families, using both a "dumb" approach (driven by numeric dynamic programming) that often works in practice, and a "clever" approach, needed for larger problems, driven by "symbolic" dynamic programming. Both approaches are fully automated and implemented in Maple.


title: Characterizing an Analogical Concept Memory for Architectures Implementing the Common Model of Cognition

abstract: Architectures that implement the Common Model of Cognition - Soar, ACT-R, and Sigma - have a prominent place in research on cognitive modeling as well as on designing complex intelligent agents. In this paper, we explore how computational models of analogical processing can be brought into these architectures to enable concept acquisition from examples obtained interactively. We propose a new analogical concept memory for Soar that augments its current system of declarative long-term memories. We frame the problem of concept learning as embedded within the larger context of interactive task learning (ITL) and embodied language processing (ELP). We demonstrate that the analogical learning methods implemented in the proposed memory can quickly learn a diverse types of novel concepts that are useful not only in recognition of a concept in the environment but also in action selection. Our approach has been instantiated in an implemented cognitive system \textsc{Aileen} and evaluated on a simulated robotic domain.


title: Multi-level Graph Convolutional Networks for Cross-platform Anchor Link Prediction

abstract: Cross-platform account matching plays a significant role in social network analytics, and is beneficial for a wide range of applications. However, existing methods either heavily rely on high-quality user generated content (including user profiles) or suffer from data insufficiency problem if only focusing on network topology, which brings researchers into an insoluble dilemma of model selection. In this paper, to address this problem, we propose a novel framework that considers multi-level graph convolutions on both local network structure and hypergraph structure in a unified manner. The proposed method overcomes data insufficiency problem of existing work and does not necessarily rely on user demographic information. Moreover, to adapt the proposed method to be capable of handling large-scale social networks, we propose a two-phase space reconciliation mechanism to align the embedding spaces in both network partitioning based parallel training and account matching across different social networks. Extensive experiments have been conducted on two large-scale real-life social networks. The experimental results demonstrate that the proposed method outperforms the state-of-the-art models with a big margin.


title: From two rolling shutters to one global shutter

abstract: Most consumer cameras are equipped with electronic rolling shutter, leading to image distortions when the camera moves during image capture. We explore a surprisingly simple camera configuration that makes it possible to undo the rolling shutter distortion: two cameras mounted to have different rolling shutter directions. Such a setup is easy and cheap to build and it possesses the geometric constraints needed to correct rolling shutter distortion using only a sparse set of point correspondences between the two images. We derive equations that describe the underlying geometry for general and special motions and present an efficient method for finding their solutions. Our synthetic and real experiments demonstrate that our approach is able to remove large rolling shutter distortions of all types without relying on any specific scene structure.


title: Disease and information spreading at different speeds in multiplex networks

abstract: Nowadays, one of the challenges we face when carrying out modeling of epidemic spreading is to develop methods to control disease transmission. In this article we study how the spreading of knowledge of a disease affects the propagation of that disease in a population of interacting individuals. For that, we analyze the interaction between two different processes on multiplex networks: the propagation of an epidemic using the susceptible-infected-susceptible dynamics and the dissemination of information about the disease --and its prevention methods-- using the unaware-aware-unaware dynamics, so that informed individuals are less likely to be infected. Unlike previous related models where disease and information spread at the same time scale, we introduce here a parameter that controls the relative speed between the propagation of the two processes. We study the behavior of this model using a mean-field approach that gives results in good agreement with Monte Carlo simulations on homogeneous complex networks. We find that increasing the rate of information dissemination reduces the disease prevalence, as one may expect. However, increasing the speed of the information process as compared to that of the epidemic process has the counter intuitive effect of increasing the disease prevalence. This result opens an interesting discussion about the effects of information spreading on disease propagation.


title: The Typology of Polysemy: A Multilingual Distributional Framework

abstract: Lexical semantic typology has identified important cross-linguistic generalizations about the variation and commonalities in polysemy patterns---how languages package up meanings into words. Recent computational research has enabled investigation of lexical semantics at a much larger scale, but little work has explored lexical typology across semantic domains, nor the factors that influence cross-linguistic similarities. We present a novel computational framework that quantifies semantic affinity, the cross-linguistic similarity of lexical semantics for a concept. Our approach defines a common multilingual semantic space that enables a direct comparison of the lexical expression of concepts across languages. We validate our framework against empirical findings on lexical semantic typology at both the concept and domain levels. Our results reveal an intricate interaction between semantic domains and extra-linguistic factors, beyond language phylogeny, that co-shape the typology of polysemy across languages.


title: Grafted network for person re-identification

abstract: Convolutional neural networks have shown outstanding effectiveness in person re-identification (re-ID). However, the models always have large number of parameters and much computation for mobile application. In order to relieve this problem, we propose a novel grafted network (GraftedNet), which is designed by grafting a high-accuracy rootstock and a light-weighted scion. The rootstock is based on the former parts of ResNet-50 to provide a strong baseline, while the scion is a new designed module, composed of the latter parts of SqueezeNet, to compress the parameters. To extract more discriminative feature representation, a joint multi-level and part-based feature is proposed. In addition, to train GraftedNet efficiently, we propose an accompanying learning method, by adding an accompanying branch to train the model in training and removing it in testing for saving parameters and computation. On three public person re-ID benchmarks (Market1501, DukeMTMC-reID and CUHK03), the effectiveness of GraftedNet are evaluated and its components are analyzed. Experimental results show that the proposed GraftedNet achieves 93.02%, 85.3% and 76.2% in Rank-1 and 81.6%, 74.7% and 71.6% in mAP, with only 4.6M parameters.


title: Interpretable Network Propagation with Application to Expanding the Repertoire of Human Proteins that Interact with SARS-CoV-2

abstract: Background: Network propagation has been widely used for nearly 20 years to predict gene functions and phenotypes. Despite the popularity of this approach, little attention has been paid to the question of provenance tracing in this context, e.g., determining how much any experimental observation in the input contributes to the score of every prediction. Results: We design a network propagation framework with two novel components and apply it to predict human proteins that directly or indirectly interact with SARS-CoV-2 proteins. First, we trace the provenance of each prediction to its experimentally validated sources, which in our case are human proteins experimentally determined to interact with viral proteins. Second, we design a technique that helps to reduce the manual adjustment of parameters by users. We find that for every top-ranking prediction, the highest contribution to its score arises from a direct neighbor in a human protein-protein interaction network. We further analyze these results to develop functional insights on SARS-CoV-2 that expand on known biology such as the connection between endoplasmic reticulum stress, HSPA5, and anti-clotting agents. Conclusions: We examine how our provenance tracing method can be generalized to a broad class of network-based algorithms. We provide a useful resource for the SARS-CoV-2 community that implicates many previously undocumented proteins with putative functional relationships to viral infection. This resource includes potential drugs that can be opportunistically repositioned to target these proteins. We also discuss how our overall framework can be extended to other, newly-emerging viruses.


title: REL: An Entity Linker Standing on the Shoulders of Giants

abstract: Entity linking is a standard component in modern retrieval system that is often performed by third-party toolkits. Despite the plethora of open source options, it is difficult to find a single system that has a modular architecture where certain components may be replaced, does not depend on external sources, can easily be updated to newer Wikipedia versions, and, most important of all, has state-of-the-art performance. The REL system presented in this paper aims to fill that gap. Building on state-of-the-art neural components from natural language processing research, it is provided as a Python package as well as a web API. We also report on an experimental comparison against both well-established systems and the current state-of-the-art on standard entity linking benchmarks.


title: Impact of indirect transitions on valley polarization in WS$_2$ and WSe$_2$

abstract: Controlling the momentum of carriers in semiconductors, known as valley polarization, is a new resource for optoelectronics and information technologies. Materials exhibiting high polarization are needed for valley-based devices. Few-layer WS$_2$ shows a remarkable spin-valley polarization above 90%, even at room temperature. In stark contrast, polarization is absent for few-layer WSe$_2$ despite the expected material similarities. Here, we explain the origin of valley polarization in both materials based on the interplay between two indirect optical transitions. We show that the relative energy minima at the $\Lambda$- and K-valleys in the conduction band determine the spin-valley polarization of the direct K-K transition. Polarization appears as the energy of the K-valley rises above the $\Lambda$-valley as a function of temperature and number of layers. Our results advance the understanding of the high spin-valley polarization in WS$_2$. This insight will impact the design of both passive and tunable valleytronic devices operating at room temperature.


title: Cold HI ejected into the Magellanic Stream

abstract: We report the direct detection of cold HI gas in a cloud ejected from the Small Magellanic Cloud (SMC) towards the Magellanic Stream. The cloud is part of a fragmented shell of HI gas on the outskirts of the SMC. This is the second direct detection of cold HI associated with the Magellanic Stream using absorption. The cold gas was detected using 21-cm HI absorption-line observations with the Australia Telescope Compact Array (ATCA) towards the extra-galactic source PMN J0029$-$7228. We find a spin (excitation) temperature for the gas of $68 \pm 20$ K. We suggest that breaking super shells from the Magellanic Clouds may be a source of cold gas to supply the rest of the Magellanic Stream.


title: Theory of cavity QED with 2D atomic arrays

abstract: We develop a quantum optical formalism to treat a two-dimensional array of atoms placed in an optical cavity. Importantly, and in contrast to typical treatments, we account for cooperative dipole-dipole effects mediated by the interaction of the atoms with the outside, non-cavity-confined modes. Based on the observation that scattering to these modes is largely inhibited due to these cooperative effects, we construct a generic formalism, independent of the specific cavity structure, and apply it to an array of non-saturated atoms. By further considering the atomic motion, we show that the inhibited damping can lead to a favorable scaling of the optomechanical parameters of an atom-array membrane placed within a cavity. The developed formalism lays the basis for further investigation of many-body QED with atom arrays in transversely confined geometries.


title: Cavity quantum optomechanics with an atom-array membrane

abstract: We consider a quantum optomechanical scheme wherein an ordered two-dimensional array of laser-trapped atoms is used as a movable membrane. The extremely light mass of the atoms yields very strong optomechanical coupling, while their spatial order largely eliminates scattering losses. We show that this combination opens the way for quantum optomechanical nonlinearities, well within the ultimate single-photon strong-coupling regime. As an example, we analyze the possibility to observe optomechanically induced quantum effects such as photon blockade and time-delayed non-classical correlations. We discuss novel opportunities opened by the optomechanical backaction on the internal states of the array atoms.


title: Countering hate on social media: Large scale classification of hate and counter speech

abstract: Hateful rhetoric is plaguing online discourse, fostering extreme societal movements and possibly giving rise to real-world violence. A potential solution to this growing global problem is citizen-generated counter speech where citizens actively engage in hate-filled conversations to attempt to restore civil non-polarized discourse. However, its actual effectiveness in curbing the spread of hatred is unknown and hard to quantify. One major obstacle to researching this question is a lack of large labeled data sets for training automated classifiers to identify counter speech. Here we made use of a unique situation in Germany where self-labeling groups engaged in organized online hate and counter speech. We used an ensemble learning algorithm which pairs a variety of paragraph embeddings with regularized logistic regression functions to classify both hate and counter speech in a corpus of millions of relevant tweets from these two groups. Our pipeline achieved macro F1 scores on out of sample balanced test sets ranging from 0.76 to 0.97---accuracy in line and even exceeding the state of the art. On thousands of tweets, we used crowdsourcing to verify that the judgments made by the classifier are in close alignment with human judgment. We then used the classifier to discover hate and counter speech in more than 135,000 fully-resolved Twitter conversations occurring from 2013 to 2018 and study their frequency and interaction. Altogether, our results highlight the potential of automated methods to evaluate the impact of coordinated counter speech in stabilizing conversations on social media.


title: Sparsification of Directed Graphs via Cut Balance

abstract: In this paper, we consider the problem of designing cut sparsifiers and sketches for directed graphs. To bypass known lower bounds, we allow the sparsifier/sketch to depend on the balance of the input graph, which smoothly interpolates between undirected and directed graphs. We give nearly matching upper and lower bounds for both for-all (cf. Bencz\'ur and Karger, STOC 1996) and for-each (Andoni et al., ITCS 2016) cut sparsifiers/sketches as a function of cut balance, defined the maximum ratio of the cut value in the two directions of a directed graph (Ene et al., STOC 2016). We also show an interesting application of digraph sparsification via cut balance by using it to give a very short proof of a celebrated maximum flow result of Karger and Levine (STOC 2002).


title: Noise robustness and experimental demonstration of a quantum generative adversarial network for continuous distributions

abstract: The potential advantage of machine learning in quantum computers is a topic of intense discussion in the literature. Theoretical, numerical and experimental explorations will most likely be required to understand its power. There has been different algorithms proposed to exploit the probabilistic nature of variational quantum circuits for generative modelling. In this paper, we employ a hybrid architecture for quantum generative adversarial networks (QGANs) and study their robustness in the presence of noise. We devise a simple way of adding different types of noise to the quantum generator circuit, and numerically simulate the noisy hybrid quantum generative adversarial networks (HQGANs) to learn continuous probability distributions, and show that the performance of HQGANs remain unaffected. We also investigate the effect of different parameters on the training time to reduce the computational scaling of the algorithm and simplify its deployment on a quantum computer. We then perform the training on Rigetti's Aspen-4-2Q-A quantum processing unit, and present the results from the training. Our results pave the way for experimental exploration of different quantum machine learning algorithms on noisy intermediate scale quantum devices.


title: Preventing Denial of Service Attacks in IoT Networks through Verifiable Delay Functions

abstract: Permissionless distributed ledgers provide a promising approach to deal with the Internet of Things (IoT) paradigm. Since IoT devices mostly generate data transactions and micropayments, distributed ledgers that use fees to regulate the network access are not an optimal choice. In this paper, we study a feeless architecture developed by IOTA and designed specifically for the IoT. Due to the lack of fees, malicious nodes can exploit this feature to generate an unbounded number of transactions and perform a denial of service attacks. We propose to mitigate these attacks through verifiable delay functions. These functions, which are non-parallelizable, hard to compute, and easy to verify, have been formulated only recently. In our work, we design a denial of service prevention mechanism which addresses network heterogeneity, limited node computational capabilities, and hardware-specific implementation optimizations. Verifiable delay functions have mostly been studied from a theoretical point of view, but little has been done in tangible applications. Hence, this paper can be considered as a pioneer work in the field, since it builds a bridge between this theoretical mathematical framework and a real-world problem.


title: Detection of spark discharges in an agitated Mars dust simulant isolated from foreign surfaces

abstract: Numerous laboratory experiments, starting in the Viking Lander era, have reported that frictional interactions between Martian analog dust grains can catalyze electrostatic processes (i.e. triboelectrification). Such findings have been cited to suggest that Martian dust devils and dust storms may sustain lightning storms, glow discharges, and other complex electrostatic phenomena. However, in many cases (if not most), these experiments allowed Martian dust simulant grains to contact foreign surfaces (for instance, the wall of an environmental chamber or other chemically dissimilar particles). A number of authors have noted that such interactions could produce charging that is not representative of processes occurring near the surface of Mars. In this work, we experimentally characterize the triboelectrification of a Martian dust simulant resulting from both isolated particle-particle collisions and particle-wall collisions under a simulated Martian environment. For the first time, we report the direct detection of spark discharges under Martian surface conditions in a flow composed solely of natural basalt and isolated from man-made surfaces. The charge densities acquired by the fluidized grains are found to be of order 10^-6 Cm^-2 (in excess of the theoretical maximum charge density for the near-surface Martian environment). Additionally, we demonstrate that the interaction of simulant particles with experimental walls can modulate the polarity of spark discharges. Our work supports the idea that small-scale spark discharges may indeed be present in Martian granular flows and may be qualitatively similar to small-scale discharges in terrestrial volcanic vents.


title: Consistent Investment of Sophisticated Rank-Dependent Utility Agents in Continuous Time

abstract: We study portfolio selection in a complete continuous-time market where the preference is dictated by the rank-dependent utility. As such a model is inherently time inconsistent due to the underlying probability weighting, we study the investment behavior of sophisticated consistent planners who seek (subgame perfect) intra-personal equilibrium strategies. We provide sufficient conditions under which an equilibrium strategy is a replicating portfolio of a final wealth. We derive this final wealth profile explicitly, which turns out to be in the same form as in the classical Merton model with the market price of risk process properly scaled by a deterministic function in time. We present this scaling function explicitly through the solution to a highly nonlinear and singular ordinary differential equation, whose existence of solutions is established. Finally, we give a necessary and sufficient condition for the scaling function to be smaller than 1 corresponding to an effective reduction in risk premium due to probability weighting.


title: On the Equivalence between Online and Private Learnability beyond Binary Classification

abstract: Alon et al. [2019] and Bun et al. [2020] recently showed that online learnability and private PAC learnability are equivalent in binary classification. We investigate whether this equivalence extends to multi-class classification and regression. First, we show that private learnability implies online learnability in both settings. Our extension involves studying a novel variant of the Littlestone dimension that depends on a tolerance parameter and on an appropriate generalization of the concept of threshold functions beyond binary classification. Second, we show that while online learnability continues to imply private learnability in multi-class classification, current proof techniques encounter significant hurdles in the regression setting. While the equivalence for regression remains open, we provide non-trivial sufficient conditions for an online learnable class to also be privately learnable.


title: Training End-to-End Analog Neural Networks with Equilibrium Propagation

abstract: We introduce a principled method to train end-to-end analog neural networks by stochastic gradient descent. In these analog neural networks, the weights to be adjusted are implemented by the conductances of programmable resistive devices such as memristors [Chua, 1971], and the nonlinear transfer functions (or `activation functions') are implemented by nonlinear components such as diodes. We show mathematically that a class of analog neural networks (called nonlinear resistive networks) are energy-based models: they possess an energy function as a consequence of Kirchhoff's laws governing electrical circuits. This property enables us to train them using the Equilibrium Propagation framework [Scellier and Bengio, 2017]. Our update rule for each conductance, which is local and relies solely on the voltage drop across the corresponding resistor, is shown to compute the gradient of the loss function. Our numerical simulations, which use the SPICE-based Spectre simulation framework to simulate the dynamics of electrical circuits, demonstrate training on the MNIST classification task, performing comparably or better than equivalent-size software-based neural networks. Our work can guide the development of a new generation of ultra-fast, compact and low-power neural networks supporting on-chip learning.


title: The magnetic helicity density patterns from non-axisymmetric solar dynamo

abstract: In the paper we study the helicity density patterns which can result from the emerging bipolar regions. Using the relevant dynamo model and the magnetic helicity conservation law we find that the helicity density pattern around the bipolar regions depends on the configuration of the ambient large-scale magnetic field, and in general they show the quadrupole distribution. The position of this pattern relative to the equator can depend on the tilt of the bipolar region. We compute the time-latitude diagrams of the helicity density evolution. The longitudinally averaged effect of the bipolar regions show two bands of sign for the density distribution in each hemisphere. Similar helicity density patterns are provided by the helicity density flux from the emerging bipolar regions subjected to the surface differential rotation. Examining effect of helicity fluxes from the bipolar regions on the large-scale dynamo we find that its effect to the dynamo saturation is negligible.


title: Quantifying the Uncertainty in Model Parameters Using Gaussian Process-Based Markov Chain Monte Carlo: An Application to Cardiac Electrophysiological Models

abstract: Estimation of patient-specific model parameters is important for personalized modeling, although sparse and noisy clinical data can introduce significant uncertainty in the estimated parameter values. This importance source of uncertainty, if left unquantified, will lead to unknown variability in model outputs that hinder their reliable adoptions. Probabilistic estimation model parameters, however, remains an unresolved challenge because standard Markov Chain Monte Carlo sampling requires repeated model simulations that are computationally infeasible. A common solution is to replace the simulation model with a computationally-efficient surrogate for a faster sampling. However, by sampling from an approximation of the exact posterior probability density function (pdf) of the parameters, the efficiency is gained at the expense of sampling accuracy. In this paper, we address this issue by integrating surrogate modeling into Metropolis Hasting (MH) sampling of the exact posterior pdfs to improve its acceptance rate. It is done by first quickly constructing a Gaussian process (GP) surrogate of the exact posterior pdfs using deterministic optimization. This efficient surrogate is then used to modify commonly-used proposal distributions in MH sampling such that only proposals accepted by the surrogate will be tested by the exact posterior pdf for acceptance/rejection, reducing unnecessary model simulations at unlikely candidates. Synthetic and real-data experiments using the presented method show a significant gain in computational efficiency without compromising the accuracy. In addition, insights into the non-identifiability and heterogeneity of tissue properties can be gained from the obtained posterior distributions.


title: The Power Graph of a Torsion-Free Group Determines the Directed Power Graph

abstract: The directed power graph $\vec{\mathcal G}(\mathbf G)$ of a group $\mathbf G$ is the simple digraph with vertex set $G$ such that $x\rightarrow y$ if $y$ is a power of $x$. The power graph of $\mathbf G$, denoted with $\mathcal G(\mathbf G)$, is the underlying simple graph. In this paper, for groups $\mathbf G$ and $\mathbf H$, the following is proved. If $\mathbf G$ has no quasicyclic subgroup $\mathbf C_{p^\infty}$ which has trivial intersection with every cyclic subgroup $\mathbf K$ of $\mathbf G$ such that $\mathbf K\not\leq\mathbf C_{p^\infty}$, then $\mathcal G(\mathbf G)\cong \mathcal G(\mathbf H)$ implies $\vec{\mathcal G}(\mathbf G)\cong \vec{\mathcal G}(\mathbf H)$. Consequently, any two torsion-free groups having isomorphic power graphs have isomorphic directed power graphs.


title: Finite Horizon Robust Synthesis Using Integral Quadratic Constraints

abstract: We present a robust synthesis algorithm for uncertain linear time-varying (LTV) systems on finite horizons. The uncertain system is described as an interconnection of a known LTV system and a perturbation. The input-output behavior of the perturbation is specified by time-domain Integral Quadratic Constraints (IQCs). The objective is to synthesize a controller to minimize the worst-case performance. This leads to a non-convex optimization. The proposed approach alternates between an LTV synthesis step and an IQC analysis step. Both induced $\mathcal{L}_2$ and terminal Euclidean norm penalties on output are considered for finite horizon performance. The proposed algorithm ensures that the robust performance is non-increasing at each iteration step. The effectiveness of this method is demonstrated using numerical examples.


title: Robust and efficient mean estimation: an approach based on the properties of self-normalized sums

abstract: Let $X$ be a random variable with unknown mean and finite variance. We present a new estimator of the mean of $X$ that is robust with respect to the possible presence of outliers in the sample, provides tight sub-Gaussian deviation guarantees without any additional assumptions on the shape or tails of the distribution, and moreover is asymptotically efficient. This is the first estimator that provably combines all these qualities in one package. Our construction is inspired by robustness properties possessed by the self-normalized sums. Theoretical findings are supplemented by numerical simulations highlighting strong performance of the proposed estimator in comparison with previously known techniques.


title: Impact-Aware Task-Space Quadratic-Programming Control

abstract: Robots usually establish contacts at rigid surfaces with near-zero relative velocities. Otherwise, impact-induced energy propagates in the robot's linkage and may cause irreversible damage to the hardware. Moreover, abrupt changes in task-space contact velocity and peak impact forces also result in abrupt changes in robot joint velocities and torques; which can compromise controllers' stability, especially for those based on smooth models. In reality, several tasks would require establishing contact with moderately high velocity. We propose to enhance task-space multi-objective controllers formulated as a quadratic program to be resilient to frictional impacts in three dimensions. We devise new constraints and reformulate the usual ones to be robust to the abrupt joint state changes mentioned earlier. The impact event becomes a controlled process once the optimal control search space is aware of: (1) the hardware-affordable impact bounds and (2) analytically-computed feasible set (polyhedra) that constrain post-impact critical states. Prior to and nearby the targeted contact spot, we assume, at each control cycle, that the impact will occur at the next iteration. This somewhat one-step preview makes our controller robust to impact time and location. To assess our approach, we experimented its resilience to moderate impacts with the Panda manipulator and achieved swift grabbing tasks with the HRP-4 humanoid robot.


title: Electron in bilayer graphene with magnetic fields leading to shape invariant potentials

abstract: The quantum behavior of electrons in bilayer graphene with applied magnetic fields is addressed. By using second-order supersymmetric quantum mechanics the problem is transformed into two intertwined one dimensional stationary Schr\"odinger equations whose potentials are required to be shape invariant. Analytical solutions for the energy bound states are obtained for several magnetic fields. The associated spectrum is analyzed, and the probability and current densities are determined.

categories: math-ph cond-mat.mes-hall hep-th math.MP quant-ph


title: Suppressed self-discharge of an aqueous supercapacitor using Earth-abundant materials

abstract: Supercapacitors (aka electrostatic double-layer capacitors -- EDLCs) offer excellent power storage capacity and kinetics, but suffer under rapid self-discharge. We introduced a zeolite framework into the active capacitor electrode, with the goal to tailor the free desorption energy and thus the self-discharge rate of supercapacitors. Low-cost carbonaceous materials and benchtop production methods were used to create supercapacitor electrodes with a measured specific capacitance of 17.25 F g$^{-1}$, a coulombic efficiency of 100%, and charge retention of over 25% over 24 hours determined by galvanostatic charge/discharge curve measurements. This charge retention was an enhancement of $\sim$350% compared with electrodes without zeolite coating.


title: Dynamic Inchworm Crawling: Performance Analysis and Optimization of a Three-link Robot

abstract: Inchworm crawling allows for both quasistatic and dynamic gaits at a wide range of actuation frequencies. This locomotion mechanism is common in nonskeletal animals and exploited extensively in the bio-inspired field of soft robotics. In this work we develop and simulate the hybrid dynamic crawling of a three-link robot, with passive frictional contacts. We fabricate and experimentally test such robot under periodic inputs of joints' angles, with good agreement to the theoretical predictions. This allows to comprehend and exploit the effects of inertia in order to find optimal performance in inputs' parameters. A simple criterion of robustness to uncertainties in friction is proposed. Tuning the inputs according to this criterion improves the robustness of low-frequency actuation, while increasing the frequency allows for gaits with both high advancement velocity and robustness. Finally, the advantages of uneven mass distribution are studied. Time-scaling technique is introduced to shape inputs that achieve similar effect without reassembling the robot. A machine-learning based optimization is applied to these inputs to further improve the robot's performance in traveling distance.


title: Detecting and Understanding Real-World Differential Performance Bugs in Machine Learning Libraries

abstract: Programming errors that degrade the performance of systems are widespread, yet there is little tool support for analyzing these bugs. We present a method based on differential performance analysis---we find inputs for which the performance varies widely, despite having the same size. To ensure that the differences in the performance are robust (i.e. hold also for large inputs), we compare the performance of not only single inputs, but of classes of inputs, where each class has similar inputs parameterized by their size. Thus, each class is represented by a performance function from the input size to performance. Importantly, we also provide an explanation for why the performance differs in a form that can be readily used to fix a performance bug. The two main phases in our method are discovery with fuzzing and explanation with decision tree classifiers, each of which is supported by clustering. First, we propose an evolutionary fuzzing algorithm to generate inputs. For this fuzzing task, the unique challenge is that we not only need the input class with the worst performance, but rather a set of classes exhibiting differential performance. We use clustering to merge similar input classes which significantly improves the efficiency of our fuzzer. Second, we explain the differential performance in terms of program inputs and internals. We adapt discriminant learning approaches with clustering and decision trees to localize suspicious code regions. We applied our techniques to a set of applications. On a set of micro-benchmarks, we show that our approach outperforms state-of-the-art fuzzers in finding inputs to characterize the differential performance. On a set of case-studies, we discover and explain multiple performance bugs in popular machine learning frameworks. Four of these bugs, reported first in this paper, have since been fixed by the developers.


title: Giant impacts stochastically change the internal pressures of terrestrial planets

abstract: Pressure is a key parameter in the physics and chemistry of planet formation and evolution. Previous studies have erroneously assumed that internal pressures monotonically increase with the mass of a body. Using smoothed particle hydrodynamics and potential field method calculations, we demonstrate that the hot, rapidly-rotating bodies produced by giant impacts can have much lower internal pressures than cool, slowly-rotating planets of the same mass. Pressures subsequently increase due to thermal and rotational evolution of the body. Using the Moon-forming impact as an example, we show that the internal pressures after the collision could have been less than half that in present-day Earth. The current pressure profile was not established until Earth cooled and the Moon receded, a process that may take up to 10s Myr after the last giant impact. Our work defines a new paradigm for pressure evolution during accretion of terrestrial planets: stochastic changes driven by impacts.


title: PILArNet: Public Dataset for Particle Imaging Liquid Argon Detectors in High Energy Physics

abstract: Rapid advancement of machine learning solutions has often coincided with the production of a test public data set. Such datasets reduce the largest barrier to entry for tackling a problem -- procuring data -- while also providing a benchmark to compare different solutions. Furthermore, large datasets have been used to train high-performing feature finders which are then used in new approaches to problems beyond that initially defined. In order to encourage the rapid development in the analysis of data collected using liquid argon time projection chambers, a class of particle detectors used in high energy physics experiments, we have produced the PILArNet, first 2D and 3D open dataset to be used for a couple of key analysis tasks. The initial dataset presented in this paper contains 300,000 samples simulated and recorded in three different volume sizes. The dataset is stored efficiently in sparse 2D and 3D matrix format with auxiliary information about simulated particles in the volume, and is made available for public research use. In this paper we describe the dataset, tasks, and the method used to procure the sample.


title: Dynamic Merkle B-tree with Efficient Proofs

abstract: We propose and define a recursive Merkle structure with q-mercurial commitments, in order to create a concise B-Merkle tree. This Merkle B-Tree builds on previous work of q-ary Merkle trees which use concise, constant size, q-mercurial commitments for intranode proofs. Although these q-ary trees reduce the branching factor and height, they still have their heights based on the key length, and are forced to fixed heights. Instead of basing nodes on q-ary prefix trees, the B Merkle Tree incorporates concise intranode commitments within a self-balancing tree. The tree is based on the ordering of elements, which requires extra information to determine element placement, but it enables significantly smaller proof sizes. This allows for much lower tree heights (directed by the order of elements, not the size of the key), and therefore creates smaller and more efficient proofs and operations. Additionally, the B Merkle Tree is defined with subset queries that feature similar communication costs to non-membership proofs. Our scheme has the potential to benefit outsourced database models, like blockchain, which use authenticated data structures and database indices to ensure immutability and integrity of data. We present potential applications in key-value stores, relational databases, and, in part, Merkle forests.


title: Chemical composition in the IRAS 16562--3959 high-mass star-forming region

abstract: We have analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) cycle 2 data of band 6 toward the G345.4938+01.4677 massive young protostellar object (G345.5+1.47 MYSO) in the IRAS 16562--3959 high-mass star-forming region with an angular resolution of $\sim 0.3"$, corresponding to $\sim 760$ au. We spatially resolve the central region which consists of three prominent molecular emission cores. A hypercompact (HC) H$_{\rm {II}}$ region (Core A) and two molecule-rich cores (Core B and Core C) are identified using the moment zero images of the H30$\alpha$ line and a CH$_{3}$OH line, respectively. Various oxygen-bearing complex organic molecules (COMs), such as (CH$_{3}$)$_{2}$CO and CH$_{3}$OCHO, have been detected toward the positions of Core B and Core C, while nitrogen-bearing species, CH$_{3}$CN, HC$_{3}$N and its $^{13}$C isotopologues, have been detected toward all of the cores. We discuss the formation mechanisms of H$_{2}$CO by comparing the spatial distribution of C$^{18}$O with that of H$_{2}$CO. The $^{33}$SO emission, on the other hand, shows a ring-like structure surrounding Core A, and it peaks on the outer edge of the H30$\alpha$ emission region. These results imply that SO is enhanced in a shock produced by the expanding motion of the ionized region.


title: Suppression of Gate Screening on Edge Magnetoplasmons by Highly Resistive ZnO Gate

abstract: We investigate a way to suppress high-frequency coupling between a gate and low-dimensional electron systems in the gigahertz range by measuring the velocity of edge magnetoplasmons (EMPs) in InAs quantum Hall systems.We compare the EMPvelocity in three samples with different electromagnetic environments-one has a highly resistive zinc oxide (ZnO) top gate, another has a normal metal (Ti/Au) top gate, and the other does not have a gate. The measured EMP velocity in the ZnO gate sample is one order of magnitude larger than that in the Ti/Au gate sample and almost the same as that in the ungated sample. As is well known, the smaller velocity in the Ti/Au gate sample is due to the screening of the electric field in EMPs. The suppression of the gate screening effect in the ZnO gate sample allows us to measure the velocity of unscreened EMPs while changing the electron density. It also offers a way to avoid unwanted high-frequency coupling between quantum Hall edge channels and gate electrodes.


title: Automatic Text Summarization of COVID-19 Medical Research Articles using BERT and GPT-2

abstract: With the COVID-19 pandemic, there is a growing urgency for medical community to keep up with the accelerating growth in the new coronavirus-related literature. As a result, the COVID-19 Open Research Dataset Challenge has released a corpus of scholarly articles and is calling for machine learning approaches to help bridging the gap between the researchers and the rapidly growing publications. Here, we take advantage of the recent advances in pre-trained NLP models, BERT and OpenAI GPT-2, to solve this challenge by performing text summarization on this dataset. We evaluate the results using ROUGE scores and visual inspection. Our model provides abstractive and comprehensive information based on keywords extracted from the original articles. Our work can help the the medical community, by providing succinct summaries of articles for which the abstract are not already available.


title: Finiteness of $\mathbb Q$-Fano compactifications of semisimple group with K\"ahler-Einstein metrics

abstract: In this note, we give a way to classify $\mathbb Q$-Fano compactifications of a semisimple group $G$. We will prove that there are only finitely many such $\mathbb Q$-Fano $G$-compactifications, which admits (singular) K\"ahler-Einstein metrics. As an application, this improves a former result in arXiv:2001.11320.


title: Measurements of an AC-LGAD strip sensor with a 120 GeV proton beam

abstract: The development of detectors that provide high resolution in four dimensions has attracted wide-spread interest in the scientific community for several applications in high-energy physics, nuclear physics, medical imaging, mass spectroscopy as well as quantum information. In addition to high time resolution and thanks to the AC-coupling of the electrodes, LGAD silicon sensors can provide high resolution in the measurement of spatial coordinates of an incident minimum ionizing particle. Such AC-coupled LGADs, also known as AC-LGADs, are therefore considered as candidates for future detectors to provide 4-dimensional measurements in a single sensing device with 100$\%$ fill factor. This article presents the first characterization of an AC-LGAD sensor with a proton beam of 120 GeV momentum at Fermilab. The sensor consists of strips with 80 $\mu$m width, fabricated at Brookhaven National Laboratory. The signal properties, efficiency, spatial, and time resolution are presented. The experimental results show that the time resolution of such an AC-LGAD is compatible to standard LGADs with similar gain, and that AC-LGADs can be segmented with fine pitches as standard strip or pixel detectors.


title: MultiXNet: Multiclass Multistage Multimodal Motion Prediction

abstract: One of the critical pieces of the self-driving puzzle is understanding the surroundings of a self-driving vehicle (SDV) and predicting how these surroundings will change in the near future. To address this task we propose MultiXNet, an end-to-end approach for detection and motion prediction based directly on lidar sensor data. This approach builds on prior work by handling multiple classes of traffic actors, adding a jointly trained second-stage trajectory refinement step, and producing a multimodal probability distribution over future actor motion that includes both multiple discrete traffic behaviors and calibrated continuous position uncertainties. The method was evaluated on large-scale, real-world data collected by a fleet of SDVs in several cities, with the results indicating that it outperforms existing state-of-the-art approaches.


title: Learning with CVaR-based feedback under potentially heavy tails

abstract: We study learning algorithms that seek to minimize the conditional value-at-risk (CVaR), when all the learner knows is that the losses incurred may be heavy-tailed. We begin by studying a general-purpose estimator of CVaR for potentially heavy-tailed random variables, which is easy to implement in practice, and requires nothing more than finite variance and a distribution function that does not change too fast or slow around just the quantile of interest. With this estimator in hand, we then derive a new learning algorithm which robustly chooses among candidates produced by stochastic gradient-driven sub-processes. For this procedure we provide high-probability excess CVaR bounds, and to complement the theory we conduct empirical tests of the underlying CVaR estimator and the learning algorithm derived from it.


title: Characteristics and performance of the Multigap Resistive Plate Chambers of the EEE experiment

abstract: The Extreme Energy Events (EEE) experiment, dedicated to the study of secondary cosmic rays, is arguably the largest detector system in the world implemented by Multigap Resistive Plate Chambers. The EEE network consists of 60 telescopes distributed over all the Italian territory; each telescope is made of three MRPCs and allows to reconstruct the trajectory of cosmic muons with high efficiency and optimal angular resolution. A distinctive feature of the EEE network is that almost all telescopes are housed in High Schools and managed by groups of students and teachers, who previously took care of their construction at CERN. This peculiarity is a big plus for the experiment, which combines the scientific relevance of its objectives with effective outreach activities. The unconventional location of the detectors, mainly in standard classrooms of school buildings, with heterogeneous maintenance conditions and without controlled temperature and dedicated power lines, is a unique test field to verify the robustness, the low aging characteristics and the long-lasting performance of MRPC technology for particle monitoring and timing. Finally, it is reported how the spatial resolution, efficiency, tracking capability and stability of these chambers behave in time.


title: Open-Set Recognition with Gaussian Mixture Variational Autoencoders

abstract: In inference, open-set classification is to either classify a sample into a known class from training or reject it as an unknown class. Existing deep open-set classifiers train explicit closed-set classifiers, in some cases disjointly utilizing reconstruction, which we find dilutes the latent representation's ability to distinguish unknown classes. In contrast, we train our model to cooperatively learn reconstruction and perform class-based clustering in the latent space. With this, our Gaussian mixture variational autoencoder (GMVAE) achieves more accurate and robust open-set classification results, with an average F1 improvement of 29.5%, through extensive experiments aided by analytical results.


title: Topological Hall Effect and Skyrmion-like Bubbles at a Charge-transfer Interface

abstract: Exploring exotic interface magnetism due to charge transfer and strong spin-orbit coupling has profound application in future development of spintronic memory. Here, the emergence, tuning and interpretation of hump-shape Hall Effect from a CaMnO3/CaIrO3/CaMnO3 trilayer structure are studied in detail. The hump signal can be recognized as Topological Hall Effect suggesting the presence of Skyrmion-like magnetic bubbles; but the debated alternative interpretation where the signal being an artefact between two cancelling Anomalous Hall Effect loops is also discussed. Firstly, by tilting the magnetic field direction, the evolution of Hall signal suggests transformation of the bubbles topology into a more trivial kind. Secondly, by varying the thickness of CaMnO3, the optimal thicknesses for the hump signal emergence are found, suggesting a tuning of charge transfer fraction. Using high-resolution transmission electron microscopy, a stacking fault is also identified, which distinguishes the top and bottom CaMnO3/CaIrO3 interfaces in terms of charge transfer fraction and possible interfacial Dzyaloshinskii-Moriya Interaction. Finally, a spin-transfer torque experiment revealed a low threshold current density of ~10^9 A/m^2 for initiating the motion of bubbles. This discovery opens a possible route for integrating Skyrmions with antiferromagnetic spintronics.