title: Unidirectional acoustic metamaterials based on nonadiabatic holonomic quantum transformations

abstract: Nonadiabatic holonomic quantum transformations (NHQTs) have attracted wide attention and have been applied in many aspects of quantum computation, whereas related research is usually limited to the field of quantum physics. Here we bring NHQTs into constructing a unidirectional acoustic metamaterial (UDAM) for shaping classical beams. The UDAM is made up of an array of three-waveguide couplers, where the propagation of acoustic waves mimics the evolution of NHQTs. The excellent agreement among analytical predictions, numerical simulations, and experimental measurements confirms the great applicability of NHQTs in acoustic metamaterial engineering. The present work extends research on NHQTs from quantum physics to the field of classical waves for designing metamaterials with simple structures and may pave a new way to design UDAMs that would be of potential applications in acoustic isolation, communication, and stealth.


title: Towards Live Video Analytics with On-Drone Deeper-yet-Compatible Compression

abstract: In this work, we present DCC(Deeper-yet-Compatible Compression), one enabling technique for real-time drone-sourced edge-assisted video analytics built on top of the existing codec. DCC tackles an important technical problem to compress streamed video from the drone to the edge without scarifying accuracy and timeliness of video analytical tasks performed at the edge. DCC is inspired by the fact that not every bit in streamed video is equally valuable to video analytics, which opens new compression room over the conventional analytics-oblivious video codec technology. We exploit drone-specific context and intermediate hints from object detection to pursue adaptive fidelity needed to retain analytical quality. We have prototyped DCC in one showcase application of vehicle detection and validated its efficiency in representative scenarios. DCC has reduced transmission volume by 9.5-fold over the baseline approach and 19-683% over the state-of-the-art with comparable detection accuracy.


title: Electrostatic forces above graphene nanoribbons and edges interpreted as partly hydrogen-free

abstract: Graphene nanoribbons' electronic transport properties strongly depend on the type of edge, armchair, zigzag or other, and on edge functionalization that can be used for band-gap engineering. For only partly hydrogenated edges interesting magnetic properties are predicted. Electric charge accumulates at edges and corners. Scanning force microscopy has so far shown the centre of graphene nanoribbons with atomic resolution using a quartz crystal tuning fork sensor of high stiffness. Weak long-range electrostatic forces related to the charge accumulation on the edges of graphene nanoribbons could not be imaged so far. Here, we show the electrostatic forces at the corners and edges of graphene nanoribbons are amenable to measurement. We use soft cantilevers and a bimodal imaging technique to combine enhanced sensitivity to weak long-range electrostatic forces with the high resolution of the second-frequency shift. Additionally, in our work the edges of the nanoribbons are mainly hydrogen-free, opening to the route to investigations of partly hydrogenated magnetic nanoribbons.


title: Dense Unsupervised Learning for Video Segmentation

abstract: We present a novel approach to unsupervised learning for video object segmentation (VOS). Unlike previous work, our formulation allows to learn dense feature representations directly in a fully convolutional regime. We rely on uniform grid sampling to extract a set of anchors and train our model to disambiguate between them on both inter- and intra-video levels. However, a naive scheme to train such a model results in a degenerate solution. We propose to prevent this with a simple regularisation scheme, accommodating the equivariance property of the segmentation task to similarity transformations. Our training objective admits efficient implementation and exhibits fast training convergence. On established VOS benchmarks, our approach exceeds the segmentation accuracy of previous work despite using significantly less training data and compute power.


title: AlphaDDA: Strategies for Adjusting the Playing Strength of a Fully Trained AlphaZero System to a Suitable Human Training Partner

abstract: Artificial intelligence (AI) has achieved superhuman performance in board games such as Go, chess, and Othello (Reversi). In other words, the AI system surpasses the level of a strong human expert player in such games. In this context, it is difficult for a human player to enjoy playing the games with the AI. To keep human players entertained and immersed in a game, the AI is required to dynamically balance its skill with that of the human player. To address this issue, we propose AlphaDDA, an AlphaZero-based AI with dynamic difficulty adjustment (DDA). AlphaDDA consists of a deep neural network (DNN) and a Monte Carlo tree search, as in AlphaZero. AlphaDDA learns and plays a game the same way as AlphaZero, but can change its skills. AlphaDDA estimates the value of the game state from only the board state using the DNN. AlphaDDA changes a parameter dominantly controlling its skills according to the estimated value. Consequently, AlphaDDA adjusts its skills according to a game state. AlphaDDA can adjust its skill using only the state of a game without any prior knowledge regarding an opponent. In this study, AlphaDDA plays Connect4, Othello, and 6x6 Othello with other AI agents. Other AI agents are AlphaZero, Monte Carlo tree search, the minimax algorithm, and a random player. This study shows that AlphaDDA can balance its skill with that of the other AI agents, except for a random player. The DDA ability of AlphaDDA is based on an accurate estimation of the value from the state of a game. We believe that the AlphaDDA approach for DDA can be used for any game AI system if the DNN can accurately estimate the value of the game state and we know a parameter controlling the skills of the AI system.


title: The Harmless Set Problem

abstract: Given a graph $G = (V,E)$, a threshold function $t~ :~ V \rightarrow \mathbb{N}$ and an integer $k$, we study the Harmless Set problem, where the goal is to find a subset of vertices $S \subseteq V$ of size at least $k$ such that every vertex $v\in V$ has less than $t(v)$ neighbors in $S$. We enhance our understanding of the problem from the viewpoint of parameterized complexity. Our focus lies on parameters that measure the structural properties of the input instance. We show that the problem is W[1]-hard parameterized by a wide range of fairly restrictive structural parameters such as the feedback vertex set number, pathwidth, treedepth, and even the size of a minimum vertex deletion set into graphs of pathwidth and treedepth at most three. On dense graphs, we show that the problem is W[1]-hard parameterized by cluster vertex deletion number. We also show that the Harmless Set problem with majority thresholds is W[1]-hard when parameterized by the treewidth of the input graph. We prove that the Harmless Set problem can be solved in polynomial time on graph with bounded cliquewidth. On the positive side, we obtain fixed-parameter algorithms for the problem with respect to neighbourhood diversity, twin cover and vertex integrity of the input graph. We show that the problem parameterized by the solution size is fixed parameter tractable on planar graphs. We thereby resolve two open questions stated in C. Bazgan and M. Chopin (2014) concerning the complexity of {\sc Harmless Set} parameterized by the treewidth of the input graph and on planar graphs with respect to the solution size.


title: Raman spectroscopy in open world learning settings using the Objectosphere approach

abstract: Raman spectroscopy in combination with machine learning has significant promise for applications in clinical settings as a rapid, sensitive, and label-free identification method. These approaches perform well in classifying data that contains classes that occur during the training phase. However, in practice, there are always substances whose spectra have not yet been taken or are not yet known and when the input data are far from the training set and include new classes that were not seen at the training stage, a significant number of false positives are recorded which limits the clinical relevance of these algorithms. Here we show that these obstacles can be overcome by implementing recently introduced Entropic Open Set and Objectosphere loss functions. To demonstrate the efficiency of this approach, we compiled a database of Raman spectra of 40 chemical classes separating them into 20 biologically relevant classes comprised of amino acids, 10 irrelevant classes comprised of bio-related chemicals, and 10 classes that the Neural Network has not seen before, comprised of a variety of other chemicals. We show that this approach enables the network to effectively identify the unknown classes while preserving high accuracy on the known ones, dramatically reducing the number of false positives while preserving high accuracy on the known classes, which will allow this technique to bridge the gap between laboratory experiments and clinical applications.

categories: cs.LG cs.AI physics.app-ph physics.bio-ph physics.comp-ph physics.med-ph


title: Photo-acoustic inversion using plasmonic contrast agents: The full Maxwell model

abstract: We analyze the inversion of the photo-acoustic imaging modality using electromagnetic plasmonic nano-particles as contrast agents. We prove that the generated pressure, before and after injecting the plasmonic nano-particles, measured at a single point, located away from the targeted inhomogeneity to image, and at a given band of incident frequencies is enough to reconstruct the (eventually complex valued) permittivity. Indeed, from these measurements, we define an indicator function which depends on the used incident frequency and the time of measurement. This indicator function has differentiating behaviors in terms of both time and frequency. First, from the behavior in terms of time, we can estimate the arrival time of the pressure from which we can localize the injected nano-particle. Second, we show that this indicator function has maximum picks at incident frequencies close to the plasmonic resonances. This allows us to estimate these resonances from which we construct the permittivity. To justify these properties, we derive the dominant electric field created by the injected nano-particle when the incident frequency is close to plasmonic resonances. This is done for the full Maxwell system. To this end, we use a natural spectral decomposition of the vector space (L^2(D))^3 based on the spectra of the Newtonian and the Magnetization operators. As another key argument, we provide the singularity analysis of the Green's tensor of the Maxwell problem with varying permittivity. Such singularity is unusual if compared to the ones of the elliptic case (as the acoustic or elastic models). In addition, we show how the derived approximation of the electric fields propagates, as a source, in the induced pressure with the time.


title: Quantum Correlations in the Minimal Scenario

abstract: In the minimal scenario of quantum correlations, two parties can choose from two observables with two possible outcomes each. Probabilities are specified by four marginals and four correlations. The resulting four-dimensional convex body of correlations, denoted $\mathcal{Q}$, is fundamental for quantum information theory. We review and systematize what is known about $\Qm$, and add many details, visualizations, and complete proofs. In particular, we provide a detailed description of the boundary, which consists of three-dimensional faces isomorphic to elliptopes and sextic algebraic manifolds of exposed extreme points. These patches are separated by cubic surfaces of non-exposed extreme points. We provide a trigonometric parametrization of all extreme points, along with their exposing Tsirelson inequalities and quantum models. All non-classical extreme points (exposed or not) are self-testing, i.e., realized by an essentially unique quantum model. Two principles, which are specific to the minimal scenario, allow a quick and complete overview: The first is the pushout transformation, i.e., the application of the sine function to each coordinate. This transforms the classical correlation polytope exactly into the correlation body $\mathcal{Q}$, also identifying the boundary structures. The second principle, self-duality, is an isomorphism between $\Qm$ and its polar dual, i.e., the set of affine inequalities satisfied by all quantum correlations (``Tsirelson inequalities''). The same isomorphism links the polytope of classical correlations contained in $\Qm$ to the polytope of no-signalling correlations, which contains $\Qm$. We also discuss the sets of correlations achieved with fixed Hilbert space dimension, fixed state or fixed observables, and establish a new non-linear inequality for $\Qm$ involving the determinant of the correlation matrix.


title: Multi-Resolution Elevation Mapping and Safe Landing Site Detection with Applications to Planetary Rotorcraft

abstract: In this paper, we propose a resource-efficient approach to provide an autonomous UAV with an on-board perception method to detect safe, hazard-free landing sites during flights over complex 3D terrain. We aggregate 3D measurements acquired from a sequence of monocular images by a Structure-from-Motion approach into a local, robot-centric, multi-resolution elevation map of the overflown terrain, which fuses depth measurements according to their lateral surface resolution (pixel-footprint) in a probabilistic framework based on the concept of dynamic Level of Detail. Map aggregation only requires depth maps and the associated poses, which are obtained from an onboard Visual Odometry algorithm. An efficient landing site detection method then exploits the features of the underlying multi-resolution map to detect safe landing sites based on slope, roughness, and quality of the reconstructed terrain surface. The evaluation of the performance of the mapping and landing site detection modules are analyzed independently and jointly in simulated and real-world experiments in order to establish the efficacy of the proposed approach.


title: Machine Learning-Based Optimization of Chiral Photonic Nanostructures: Evolution- and Neural Network-Based Design

abstract: Chiral photonics opens new pathways to manipulate light-matter interactions and tailor the optical response of meta-surfaces and -materials by nanostructuring nontrivial patterns. Chirality of matter, such as that of molecules, and light, which in the simplest case is given by the handedness of circular polarization, have attracted much attention for applications in chemistry, nanophotonics and optical information processing. We report the design of chiral photonic structures using two machine learning methods, the evolutionary algorithm and neural network approach, for rapid and efficient optimization of optical properties for dielectric metasurfaces. The design recipes obtained for visible light in the range of transition-metal dichalcogenide exciton resonances show a frequency-dependent modification in the reflected light's degree of circular polarization, that is represented by the difference between left- and right-circularly polarized intensity. Our results suggest the facile fabrication and characterization of optical nanopatterned reflectors for chirality-sensitive light-matter coupling scenarios employing tungsten disulfide as possible active material with features such as valley Hall effect and optical valley coherence.


title: Quasinormal modes of charged black holes with corrections from nonlinear electrodynamics

abstract: We study quasinormal modes related to gravitational and electromagnetic perturbations of spherically symmetric charged black holes in nonlinear electrodynamics. Beyond the linear Maxwell electrodynamics, we consider a class of Lagrangian with higher-order corrections written by the electromagnetic field strength and its Hodge dual with arbitrary coefficients, and we parametrize the corrections for quasinormal frequencies in terms of the coefficients. It is confirmed that the isospectrality of quasinormal modes under parity is generally violated due to nonlinear electrodynamics. As applications, the corrections for quasinormal frequencies in Euler-Heisenberg and Born-Infeld electrodynamics are calculated, then it is clarified that the nonlinear effects act to lengthen the oscillation period and enhance the damping rate of the quasinormal modes.


title: Energy spreading, equipartition and chaos in lattices with non-central forces

abstract: We numerically study a one dimensional, nonlinear lattice model which in the linear limit is relevant to the study of bending (flexural) waves. In contrast with the classic one dimensional mass-spring system, the linear dispersion relation of the considered model has different characteristics in the low frequency limit. By introducing disorder in the masses of the lattice particles, we investigate how different nonlinearities (cubic, quartic and their combination) lead to energy delocalization, equipartition and chaotic dynamics. We excite the lattice using single site initial momentum excitations corresponding to a strongly localized linear mode and increase the initial energy of excitation. Beyond a certain energy threshold, when the cubic nonlinearity is present, the system is found to reach energy equipartition and total delocalization. On the other hand, when only the quartic nonlinearity is activated, the system remains localized and away from equipartition at least for the energies and evolution times considered here. However, for large enough energies for all types of nonlinearities we observe chaos. This chaotic behavior is combined with energy delocalization when cubic nonlinearities are present, while the appearance of only quartic nonlinearity leads to energy localization. Our results reveal a rich dynamical behavior and show differences with the relevant Fermi-Pasta-Ulam-Tsingou model. Our findings pave the way for the study of models relevant to bending (flexural) waves in the presence of nonlinearity and disorder, anticipating different energy transport behaviors.


title: Exceptional points of PT-symmetric reflectionless states in complex scattering systems

abstract: We investigate experimentally and analytically the coalescence of reflectionless (RL) states in symmetric complex wave-scattering systems. We observe RL-exceptional points (EPs), first, with a conventional Fabry-Perot system for which the scattering strength within the system is tuned symmetrically, and then with single- and multi-channel symmetric disordered systems. We identify that an EP of the parity-time (PT)-symmetric RL-operator is obtained when the spacing between central frequencies of two natural resonances of the system is equal to the decay rate into incoming and outgoing channels. Finally, we leverage the transfer functions associated with RL and RL-EPs states to implement first- and second-order analog differentiation.


title: 6D Pose Estimation with Combined Deep Learning and 3D Vision Techniques for a Fast and Accurate Object Grasping

abstract: Real-time robotic grasping, supporting a subsequent precise object-in-hand operation task, is a priority target towards highly advanced autonomous systems. However, such an algorithm which can perform sufficiently-accurate grasping with time efficiency is yet to be found. This paper proposes a novel method with a 2-stage approach that combines a fast 2D object recognition using a deep neural network and a subsequent accurate and fast 6D pose estimation based on Point Pair Feature framework to form a real-time 3D object recognition and grasping solution capable of multi-object class scenes. The proposed solution has a potential to perform robustly on real-time applications, requiring both efficiency and accuracy. In order to validate our method, we conducted extensive and thorough experiments involving laborious preparation of our own dataset. The experiment results show that the proposed method scores 97.37% accuracy in 5cm5deg metric and 99.37% in Average Distance metric. Experiment results have shown an overall 62% relative improvement (5cm5deg metric) and 52.48% (Average Distance metric) by using the proposed method. Moreover, the pose estimation execution also showed an average improvement of 47.6% in running time. Finally, to illustrate the overall efficiency of the system in real-time operations, a pick-and-place robotic experiment is conducted and has shown a convincing success rate with 90% of accuracy. This experiment video is available at https://sites.google.com/view/dl-ppf6dpose/.


title: A weak to strong type T1 theorem for general smooth Calder\'on-Zygmund operators with doubling weights, II

abstract: We consider the weak to strong type problem for two weight norm inequalities for Calder\'on-Zygmund operators with doubling weights. We show that if a Calder\'on-Zygmund operator T is weak type (2,2) with doubling weights, then it is strong type (2,2) if and only if the dual cube testing condition for T^{*} holds, alternatively if and only if the dual cancellation condition of Stein holds. The testing condition can be taken with respect to either cubes or balls, and more generally, this is extended to a weak form of Tb theorem. Finally, we show that for all pairs of locally finite positive Borel measures, and all Stein elliptic Calder\'on-Zygmund operators T, the weak type (2,2) inequalities for T and and its associated maximal truncations operator T_{*} are equivalent. Thus the characterization of weak type for T_{*} in [LaSaUr1] applies to T as well.


title: On Nash-solvability of finite $n$-person deterministic graphical games; Catch 22

abstract: We consider finite $n$-person deterministic graphical (DG) games. These games are modelled by finite directed graphs (digraphs) $G$ which may have directed cycles and, hence, infinite plays. Yet, it is assumed that all these plays are equivalent and form a single outcome $c$, while the terminal vertices $V_T = \{a_1, \ldots, a_p\}$ form $p$ remaining outcomes. We study the existence of Nash equilibria (NE) in pure stationary strategies. It is known that NE exist when $n=2$ and may fail to exist when $n > 2$. Yet, the question becomes open for $n > 2$ under the following extra condition: (C) For each of $n$ players, $c$ is worse than each of $p$ terminal outcomes. In other words, all players are interested in terminating the play, which is a natural assumption. Moreover, Nash-solvability remains open even if we replace (C) by a weaker condition: (C22) There exist no two players for whom $c$ is better than (at least) two terminal outcomes. We conjecture that such two players exist in each NE-free DG game, or in other words, that (C22) implies Nash-solvability, for all $n$. Recently, the DG games were extended to a wider class of the DG multi-stage (DGMS) games, whose outcomes are the strongly connected components (SCC) of digraph $G$. Merging all outcomes of a DGMS game that correspond to its non-terminal SCCs we obtain a DG game. Clearly, this operation respects Nash-solvability (NS). Basic conditions and conjectures related to NS can be extended from the DG to DGMS games: in both cases NE exist if $n=2$ and may fail to exist when $n > 2$; furthermore, we modify conditions (C) and (C22) to adapt them for the DGMS games. Keywords: $n$-person deterministic graphical (multi-stage) games, Nash equilibrium, Nash-solvability, pure stationary strategy, digraph, directed cycle, strongly connected component.


title: An ensemble approach to the structure-function problem in microbial communities

abstract: The metabolic activity of microbes has played an essential role in the evolution and persistence of life on Earth. Microbial metabolism plays a primary role in the flow of carbon, nitrogen and other elements through the biosphere on a global scale. Microbes perform these metabolic activities in the context of complex communities comprised of many species that interact in dynamic and spatially-structured environments. Molecular genetics has revealed many of the metabolic pathways microbes utilize to generate energy and biomass. However, most of this knowledge is derived from model organisms, so we have a limited view of role of the massive genomic diversity in the wild on metabolic phenotypes. Further, we are only beginning to glimpse the principles governing how the metabolism of a community emerges from the collective action of its constituent members. As a result, one of the biggest challenges in the field is to understand how the metabolic activity of a community emerges from the genomic structure of the constituents. Here we propose an approach to this problem that rests on the quantitative analysis of metabolic activity in ensembles of microbial communities. We propose quantifying metabolic fluxes in diverse communities, either in the laboratory or the wild. We suggest that using sequencing data to quantify the genomic, taxonomic or transcriptional variation across an ensemble of communities can reveal low-dimensional descriptions of community structure that can explain or predict their emergent metabolic activity. We survey the types of communities for which this approach might be best suited and then review the analytical techniques available for quantifying metabolite dynamics in communities. Finally, we discuss what types of data analysis approaches might be lucrative for learning the structure-function mapping in communities from these data.


title: Scale dependent anisotropy of electric field fluctuations in solar wind turbulence

abstract: We study the variation of average powers and spectral indices of electric field fluctuations with respect to the angle between average flow direction and the mean magnetic field in solar wind turbulence. Cluster spacecraft data from the years 2002 and 2007 are used for the present analysis. We perform a scale dependent study with respect to the local mean magnetic field using wavelet analysis technique. Prominent anisotropies are found for both the spectral index and power levels of the electric power spectra. Similar to the magnetic field fluctuations, the parallel (or antiparallel) electric fluctuation spectrum is found to be steeper than the perpendicular spectrum. However the parallel (or antiparallel) electric power is found to be greater than the perpendicular one. Below 0.1 Hz, the slope of the parallel electric power spectra deviates substantially from that of the total magnetic power spectra, supporting the existence of Alfv\'enic turbulence.


title: Nonconvex flexible sparsity regularization: theory and monotone numerical schemes

abstract: Flexible sparsity regularization means stably approximating sparse solutions of operator equations by using coefficient-dependent penalizations. We propose and analyse a general nonconvex approach in this respect, from both theoretical and numerical perspectives. Namely, we show convergence of the regularization method and establish convergence properties of a couple of majorization approaches for the associated nonconvex problems. We also test a monotone algorithm for an academic example where the operator is an $M$ matrix, and on a time-dependent optimal control problem, pointing out the advantages of employing variable penalties over a fixed penalty.


title: Compositions and decompositions of binary relations

abstract: It is well-known that to every binary relation on a non-void set I there can be assigned its incidence matrix, also in the case when I is infinite. We show that a certain kind of "multiplication" of such incidence matrices corresponds to the composition of the corresponding relations. Using this fact we investigate the solvability of the equation R o X = S for given binary relations R and S on I and derive an algorithm for solving this equation by using the connections between the corresponding incidence matrices. Moreover, we describe how one can obtain the incidence matrix of a product of binary relations from the incidence matrices of its factors.


title: DropGNN: Random Dropouts Increase the Expressiveness of Graph Neural Networks

abstract: This paper studies Dropout Graph Neural Networks (DropGNNs), a new approach that aims to overcome the limitations of standard GNN frameworks. In DropGNNs, we execute multiple runs of a GNN on the input graph, with some of the nodes randomly and independently dropped in each of these runs. Then, we combine the results of these runs to obtain the final result. We prove that DropGNNs can distinguish various graph neighborhoods that cannot be separated by message passing GNNs. We derive theoretical bounds for the number of runs required to ensure a reliable distribution of dropouts, and we prove several properties regarding the expressive capabilities and limits of DropGNNs. We experimentally validate our theoretical findings on expressiveness. Furthermore, we show that DropGNNs perform competitively on established GNN benchmarks.


title: Optimal Physical Sorting of Mobile Agents

abstract: Given a collection of red and blue mobile agents located on two grid rows, we seek to move all the blue agents to the far left side and all the red agents to the far right side, thus \textit{physically sorting} them according to color. The agents all start on the bottom row. They move simultaneously at discrete time steps and must not collide. Our goal is to design a centralized algorithm that controls the agents so as to sort them in the least number of time steps. We derive an \textbf{exact} lower bound on the amount of time any algorithm requires to sort a given initial configuration of agents. We find an instance optimal algorithm that provably matches this lower bound, attaining the best possible sorting time for any initial configuration. Surprisingly, we find that whenever the leftmost agent is red and the rightmost agent is blue, a straightforward decentralized and local sensing-based algorithm is at most $1$ time step slower than the centralized instance-optimal algorithm.


title: Stable solutions to the fractional Allen-Cahn equation in the nonlocal perimeter regime

abstract: We study stable solutions to the fractional Allen-Cahn equation \linebreak $(-\Delta)^{s/2} u = u-u^3$, $|u|<1$ in $\mathbb{R}^n$. For every $s\in (0,1)$ and dimension $n\geq 2$, we establish sharp energy estimates, density estimates, and the convergence of blow-downs to stable nonlocal $s$-minimal cones. As a consequence, we obtain a new classification result: if for some pair $(n,s)$, with $n\ge 3$, hyperplanes are the only stable nonlocal $s$-minimal cones in $\mathbb{R}^n\setminus\{0\}$, then every stable solution to the fractional Allen-Cahn equation in $\mathbb{R}^n$ is 1D, namely, its level sets are parallel hyperplanes. Combining this result with the classification of stable $s$-minimal cones in $\mathbb{R}^3\setminus\{0\}$ for $s\sim 1$ obtained by the authors in a recent paper, we give positive answers to the "stability conjecture" in $\mathbb{R}^3$ and to the "De Giorgi conjecture" in $\mathbb{R}^4$ for the fractional Allen-Cahn equation when the order $s\in (0,1)$ of the operator is sufficiently close to $1$.


title: Comparison of the semiclassical and quantum optical field dynamics in a pulse-excited optical cavity with a finite number of quantum emitters

abstract: The spectral and temporal response of a set of $N$ quantum emitters embedded in a photonic cavity is studied. Quantum mechanically, such systems can be described by the Tavis-Cummings (TC) model of $N$ two-level systems coupled to a single light mode. Here we compare the full quantum solution of the TC model for different numbers of quantum emitters with its semiclassical limit after a pulsed excitation of the cavity mode. Considering different pulse amplitudes, we find that the spectra obtained from the TC model approach the semiclassical one for an increasing number of emitters $N$. Furthermore they match very well for small pulse amplitudes. While we observe a very good agreement in the temporal dynamics for photon numbers much smaller than $N$, considerable deviations occur in the regime of photon numbers similar to or larger than $N$, which are linked to collapse and revival phenomena. Wigner functions of the light mode are calculated for different scenarios to analyze the quantum state of the light field. We find strong deviations from a coherent state even if the dynamics of the expectation values are still well described by the semiclassical limit. For higher pulse amplitudes Wigner functions similar to those of Schr\"odinger cat states between two or more quasi-coherent contributions build up.


title: Relaxation for an optimal design problem in $BD(\Omega)$

abstract: We obtain a measure representation for a functional arising in the context of optimal design problems under linear growth conditions. The functional in question corresponds to the relaxation with respect to a pair $(\chi,u)$, where $\chi$ is the characteristic function of a set of finite perimeter and $u$ is a function of bounded deformation, of an energy with a bulk term depending on the symmetrised gradient as well as a perimeter term.


title: MaTIC: a Mathematical Theory of Inferential Communication

abstract: The arrival of Immersive Virtual and Augmented Reality hardware to the consumer market suggests seamless multi-modal communication between human participants and autonomous interactive characters is an achievable goal in the near future. This possibility is further reinforced by the rapid improvements in the automated analysis of speech, facial expressions and body language, as well as improvements in character animation and speech synthesis techniques. However, we do not have a formal theory that allows us to compare, on one side, interactive social scenarios among human users and autonomous virtual characters and, on the other side, pragmatic inference mechanisms as they occur in non-mediated communication. Grices' and Sperbers' model of inferential communication does explain the nature of everyday communication through cognitive mechanisms that support spontaneous inferences performed in pragmatic communication. However, such a theory is not based on a mathematical framework with a precision comparable to classical information theory. To address this gap, in this article we introduce a Mathematical Theory of Inferential Communication (MaTIC). MaTIC formalises some assumptions of inferential communication, it explores its theoretical consequences and outlines the practical steps needed to use it in different application scenarios.


title: Controlling transition metal atomic ordering in two-dimensional Mo$_{1-x}$W$_{x}$S$_{2}$ alloys

abstract: The unique optical and electronic properties of two-dimensional transition metal dichalcogenides (2D TMDs) make them promising materials for applications in (opto-)electronics, catalysis and more. Specifically, alloys of 2D TMDs have broad potential applications owing to their composition-controlled properties. Several important challenges remain regarding controllable and scalable fabrication of these alloys, such as achieving control over their atomic ordering (i.e. clustering or random mixing of the transition metal atoms within the 2D layers). In this work, atomic layer deposition (ALD) is used to synthesize the TMD alloy Mo$_{1-x}$W$_{x}$S$_{2}$ with excellent composition control along the complete composition range 0 $\leq$ x $\leq$ 1. Importantly, this composition control allows us to control the atomic ordering of the alloy from well-mixed to clustered while keeping the alloy composition fixed, as is confirmed directly through atomic-resolution HAADF-STEM imaging. The control over atomic ordering leads to tuning of the bandgap, as is demonstrated using optical transmission spectroscopy. The relation between this tuning of the electronic structure and the atomic ordering of the alloy was further confirmed through ab-initio calculations. Furthermore, as the atomic ordering modulates from clustered to well-mixed, the typical MoS$_{2}$ and WS$_{2}$ A$_{1g}$ vibrational modes converge. Our results demonstrate that atomic ordering is an important parameter that can be tuned experimentally to finely tune the fundamental properties of 2D TMD alloys for specific applications.


title: Fairness, Integrity, and Privacy in a Scalable Blockchain-based Federated Learning System

abstract: Federated machine learning (FL) allows to collectively train models on sensitive data as only the clients' models and not their training data need to be shared. However, despite the attention that research on FL has drawn, the concept still lacks broad adoption in practice. One of the key reasons is the great challenge to implement FL systems that simultaneously achieve fairness, integrity, and privacy preservation for all participating clients. To contribute to solving this issue, our paper suggests a FL system that incorporates blockchain technology, local differential privacy, and zero-knowledge proofs. Our implementation of a proof-of-concept with multiple linear regression illustrates that these state-of-the-art technologies can be combined to a FL system that aligns economic incentives, trust, and confidentiality requirements in a scalable and transparent system.


title: Related Work on Image Quality Assessment

abstract: Due to the existence of quality degradations introduced in various stages of visual signal acquisition, compression, transmission and display, image quality assessment (IQA) plays a vital role in image-based applications. According to whether the reference image is complete and available, image quality evaluation can be divided into three categories: Full-Reference(FR), Reduced- Reference(RR), and Non- Reference(NR). This article will review the state-of-the-art image quality assessment algorithms.


title: Alternating twisted mutilayer graphene: generic partition rules, double flat bands, and orbital magnetoelectric effect

abstract: Twisted graphene systems have draw significant attention due to the discoveries of various correlated and topological phases. In particular, recently the alternating twisted trilayer graphene is discovered to exhibit unconventional superconductivity, which motivates us to study the electronic structures and possible interesting correlation effects of this class of alternating twisted graphene systems. In this work we consider generic alternating twisted multilayer graphene (ATMG) systems with $M$-$L$-$N$ stacking configurations, in which the $M$ ($L$) graphene layers and the $L$ ($N$) layers are twisted by an angle $\theta$ (-$\theta$). Based on analysis from a simplified $\textbf{k}\!\cdot\!\textbf{p}$ model approach, we analytically derive generic partition rules for the low-energy electronic structures, which exhibit various intriguing band dispersions including one pair of flat bands, two pairs of flat bands, as well as flat bands co-existing with with Dirac cones, quadratic bands, or more generally $E(\mathbf{k})\!\sim\!k^J$ dispersions ($J$ is positive integer) for each spin and valley. Such unusual non-interacting electronic structures may have unconventional correlation effects. Especially for a mirror symmetric ATMG system with two pairs of flat bands (per spin per valley), we find that Coulomb interactions may drive the system into a state breaking both time-reversal and mirror symmetries, which can exhibit a novel type of orbital magnetoelectric effect due to the interwining of electric polarization and orbital magnetization orders in the symmetry-breaking state.


title: Towards understanding the magnetic properties of the breathing pyrochlore compound Ba3Yb2Zn5O11: A single crystal study

abstract: Ba3Yb2Zn5O11 is unique among breathing pyrochlore compounds for being in the nearly decoupled limit where inter-tetrahedron interactions are weak, hosting isolated clusters or "molecular magnet" like tetrahedra of magnetic ytterbium (Yb3+) ions. In this work, we present the first study carried out on single-crystal samples of the breathing pyrochlore Ba3Yb2Zn5O11, using a variety of magnetometry and neutron scattering techniques along with theoretical modeling. We employ inelastic neutron scattering to investigate the magnetic dynamics as a function of applied field (with respect to both magnitude and direction) down to a temperature of 70 mK, where inelastic scattering reveals dispersionless bands of excitations as found in earlier powder sample studies, in good agreement with a single-tetrahedron model. However, diffuse neutron scattering at zero field and dc-susceptibility at finite field exhibit features suggesting the presence of excitations at low-energy that are not captured by the single tetrahedron model. Analysis of the local structure down to 2 K via pair distribution function analysis finds no evidence of structural disorder. We conclude that effects beyond the single tetrahedron model are important in describing the low-energy, low temperature physics of Ba3Yb2Zn5O11, but their nature remains undetermined.


title: On the Maximal Number of Columns of a $\Delta$-modular Integer Matrix: Bounds and Computations

abstract: We study the maximal number of pairwise distinct columns in a $\Delta$-modular integer matrix with $m$ rows. Recent results by Lee et al. provide an asymptotically tight upper bound of $O(m^2)$ for fixed $\Delta$. We complement this and obtain an upper bound of the form $O(\Delta)$ for fixed $m$, and with the implied constant depending polynomially on $m$.


title: On constraint preservation and strong hyperbolicity

abstract: We use partial differential equations (PDEs) to describe physical systems. In general, these equations include evolution and constraint equations. One method used to find solutions to these equations is the Free-evolution approach, which consists in obtaining the solutions of the entire system by solving only the evolution equations. Certainly, this is valid only when the chosen initial data satisfies the constraints and the constraints are preserved in the evolution. In this paper, we establish the sufficient conditions required for the PDEs of the system to guarantee the constraint preservation. This is achieved by considering quasi-linear first-order PDEs, assuming the sufficient condition and deriving strongly hyperbolic first-order partial differential evolution equations for the constraints. We show that, in general, these constraint evolution equations correspond to a family of equations parametrized by a set of free parameters. We also explain how these parameters fix the propagation velocities of the constraints. As application examples of this framework, we study the constraint conservation of the Maxwell electrodynamics and the wave equations in arbitrary space-times. We conclude that the constraint evolution equations are unique in the Maxwell case and a family in the wave equation case.


title: Search for supersymmetry in final states with two or three soft leptons and missing transverse momentum in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search for supersymmetry in events with two or three low-momentum leptons and missing transverse momentum is performed. The search uses proton-proton collisions at $\sqrt{s} =$ 13 TeV collected in the three-year period 2016-2018 by the CMS experiment at the LHC and corresponding to an integrated luminosity of up to 137 fb$^{-1}$. The data are found to be in agreement with expectations from standard model processes. The results are interpreted in terms of electroweakino and top squark pair production with a small mass difference between the produced supersymmetric particles and the lightest neutralino. For the electroweakino interpretation, two simplified models are used, a wino-bino model and a higgsino model. Exclusion limits at 95% confidence level are set on $\widetilde{\chi}^0_2 / \widetilde{\chi}^\pm_1$ masses up to 275 GeV for a mass difference of 10 GeV in the wino-bino case, and up to 205 (150) GeV for a mass difference of 7.5 (3) GeV in the higgsino case. The results for the higgsino are further interpreted using a phenomenological minimal supersymmetric standard model, excluding the higgsino mass parameter $\mu$ up to 180 GeV with the bino mass parameter $M_1$ at 800 GeV. In the top squark interpretation, exclusion limits are set at top squark masses up to 540 GeV for four-body top squark decays and up to 480 GeV for chargino-mediated decays with a mass difference of 30 GeV.


title: Bourgain-Brezis-Mironescu-Maz'ya-Shaposhnikova limit formulae for fractional Sobolev spaces via interpolation and extrapolation

abstract: The real interpolation spaces between $L^{p}({\mathbb{R}}^{n})$ and $\dot {H}^{t,p}({\mathbb{R}}^{n})$ (resp. $H^{t,p}({\mathbb{R}}^{n})$), $t>0,$ are characterized in terms of fractional moduli of smoothness, and the underlying seminorms are shown to be " the correct" fractional generalization of the classical Gagliardo seminorms. This is confirmed by the fact that, using the new spaces combined with interpolation and extrapolation methods, we are able to extend the Bourgain-Brezis-Mironescu-Maz'ya-Shaposhnikova limit formulae, as well as the Bourgain-Brezis-Mironescu convergence theorem, to fractional Sobolev spaces. On the other hand, we disprove a conjecture of \cite{Braz} suggesting fractional convergence results given in terms of classical Gagliardo seminorms. We also solve a problem proposed in \cite{Braz} concerning sharp forms of the fractional Sobolev embedding.


title: The Entanglement Wedge of Unknown Couplings

abstract: The black hole interior is a mysterious region of spacetime where non-perturbative effects are sometimes important. These non-perturbative effects are believed to be highly theory-dependent. We sharpen these statements by considering a setup where the state of the black hole is in a superposition of states corresponding to boundary theories with different couplings, entangled with a reference which keeps track of those couplings. The entanglement wedge of the reference can then be interpreted as the bulk region most sensitive to the values of the couplings. In simple bulk models, e.g., JT gravity + a matter BCFT, the QES formula implies that the reference contains the black hole interior at late times. We also analyze the Renyi-2 entropy of the reference, which can be viewed as a diagnostic of chaos via the Loschmidt echo. We find explicitly the replica wormhole that diagnoses the island and restores unitarity. Numerical and analytical evidence of these statements in the SYK model is presented. Similar considerations are expected to apply in higher dimensional AdS/CFT, for marginal and even irrelevant couplings.


title: Approximating Sparsest Cut in Low-Treewidth Graphs via Combinatorial Diameter

abstract: The fundamental sparsest cut problem takes as input a graph $G$ together with the edge costs and demands, and seeks a cut that minimizes the ratio between the costs and demands across the cuts. For $n$-node graphs~$G$ of treewidth~$k$, \chlamtac, Krauthgamer, and Raghavendra (APPROX 2010) presented an algorithm that yields a factor-$2^{2^k}$ approximation in time $2^{O(k)} \cdot \operatorname{poly}(n)$. Later, Gupta, Talwar and Witmer (STOC 2013) showed how to obtain a $2$-approximation algorithm with a blown-up run time of $n^{O(k)}$. An intriguing open question is whether one can simultaneously achieve the best out of the aforementioned results, that is, a factor-$2$ approximation in time $2^{O(k)} \cdot \operatorname{poly}(n)$. In this paper, we make significant progress towards this goal, via the following results: (i) A factor-$O(k^2)$ approximation that runs in time $2^{O(k)} \cdot \operatorname{poly}(n)$, directly improving the work of Chlamt\'a\v{c} et al. while keeping the run time single-exponential in $k$. (ii) For any $\varepsilon>0$, a factor-$O(1/\varepsilon^2)$ approximation whose run time is $2^{O(k^{1+\varepsilon}/\varepsilon)} \cdot \operatorname{poly}(n)$, implying a constant-factor approximation whose run time is nearly single-exponential in $k$ and a factor-$O(\log^2 k)$ approximation in time $k^{O(k)} \cdot \operatorname{poly}(n)$. Key to these results is a new measure of a tree decomposition that we call combinatorial diameter, which may be of independent interest.


title: Ground-state properties of the narrowest zigzag graphene nanoribbon from quantum Monte Carlo and comparison with density functional theory

abstract: By means of quantum Monte Carlo (QMC) calculations from first principles, we study the ground-state properties of the narrowest zigzag graphene nanoribbon, with an infinite linear acene structure. We show that this quasi-one-dimensional system is correlated and its ground state is made of localized $\pi$ electrons whose spins are antiferromagnetically (AFM) ordered. The AFM stablization energy (36(3) meV per carbon atom) and the absolute magnetization (1.13(1) $\mu_\textrm{B}$ per unit cell) predicted by QMC are sizable, and they suggest the survival of antiferromagnetic correlations above room temperature. These values can be reproduced to some extent by density functional theory (DFT) only by assuming strong interactions, either within the DFT+U framework or using hybrid functionals. Based on our QMC results, we then provide the strength of Hubbard repulsion in DFT+U suitable for this class of systems.


title: Developing and Publishing Code for Trusted Research Environments: Best Practices and Ways of Working

abstract: This report discusses 3 distinct, but overlapping topics. Firstly, it recommends the tools and best practices for research software engineering and data science that are most relevant to the researchers working on the Wales Multimorbidity Machine Learning (WMML) project. Secondly, it expands upon these recommendations for the specific use case of Trusted Research Environments (TREs), with development workflows for computational research in TREs offered that respect and complement existing best practices. Finally, it discusses the considerations around publishing research code that is developed to run within a TRE on sensitive data, offering practical advice that researchers using TREs can follow.


title: On Recovering the Best Rank-r Approximation from Few Entries

abstract: In this note, we investigate how well we can reconstruct the best rank-$r$ approximation of a large matrix from a small number of its entries. We show that even if a data matrix is of full rank and cannot be approximated well by a low-rank matrix, its best low-rank approximations may still be reliably computed or estimated from a small number of its entries. This is especially relevant from a statistical viewpoint: the best low-rank approximations to a data matrix are often of more interest than itself because they capture the more stable and oftentimes more reproducible properties of an otherwise complicated data-generating model. In particular, we investigate two agnostic approaches: the first is based on spectral truncation; and the second is a projected gradient descent based optimization procedure. We argue that, while the first approach is intuitive and reasonably effective, the latter has far superior performance in general. We show that the error depends on how close the matrix is to being of low rank. Both theoretical and numerical evidence is presented to demonstrate the effectiveness of the proposed approaches.


title: Comment on "Shear-free barotropic perfect fluids cannot rotate and expand simultaneously" by R. Goswami and G.F.R. Ellis

abstract: We point out an error in a recent paper by Goswami and Ellis. As a consequence the question of whether shear-free barotropic perfect fluids (with $p+\mu\neq 0$) can or cannot rotate and expand simultaneously, is still wide open


title: Joint Radar-Communications Processing from a Dual-Blind Deconvolution Perspective

abstract: We consider a general spectral coexistence scenario, wherein the channels and transmit signals of both radar and communications systems are unknown at the receiver. In this \textit{dual-blind deconvolution} (DBD) problem, a common receiver admits the multi-carrier wireless communications signal that is overlaid with the radar signal reflected-off multiple targets. When the radar receiver is not collocated with the transmitter, such as in passive or multistatic radars, the transmitted signal is also unknown apart from the target parameters. Similarly, apart from the transmitted messages, the communications channel may also be unknown in dynamic environments such as vehicular networks. As a result, the estimation of unknown target and communications parameters in a DBD scenario is highly challenging. In this work, we exploit the sparsity of the channel to solve DBD by casting it as an atomic norm minimization problem. Our theoretical analyses and numerical experiments demonstrate perfect recovery of continuous-valued range-time and Doppler velocities of multiple targets as well as delay-Doppler communications channel parameters using uniformly-spaced time samples in the dual-blind receiver.


title: Two regimes of confinement in photonic nanocavities: bulk confinement versus lightning rods

abstract: We present a theoretical study of dielectric bowtie cavities and show that they are governed by two essentially different confinement regimes. The first is confinement inside the bulk dielectric and the second is a local lightning-rod regime where the field is locally enhanced at sharp corners and may yield a vanishing mode volume without necessarily enhancing the mode inside the bulk dielectric. We show that while the bulk regime is reminiscent of the confinement in conventional nanocavities, the most commonly used definition of the mode volume gauges in fact the lightning-rod effect when applied to ultra-compact cavities, such as bowties. Distinguishing between these two regimes will be crucial for future research on nanocavities, and our insights show how to obtain strongly enhanced light-matter interaction over large bandwidths.


title: Automatically identifying a mobile phone user's position within a vehicle

abstract: Traffic-related injuries and fatalities are major health risks in the United States. Mobile phone use while driving quadruples the risk for a motor vehicle crash. This work demonstrates the feasibility of using the mobile phone camera to passively detect the location of the phone's user within a vehicle. In a large, varied dataset we were able correctly identify if the user was in the driver's seat or one of the passenger seats with 94.9% accuracy. This model could be used by application developers to selectively change or lock functionality while a user is driving, but not if the user is a passenger in a moving vehicle.


title: Logical limit laws for layered permutations and related structures

abstract: We show that several classes of ordered structures (namely, convex linear orders, layered permutations, and compositions) admit first-order logical limit laws.


title: Online Discrepancy with Recourse for Vectors and Graphs

abstract: The vector-balancing problem is a fundamental problem in discrepancy theory: given T vectors in $[-1,1]^n$, find a signing $\sigma(a) \in \{\pm 1\}$ of each vector $a$ to minimize the discrepancy $\| \sum_{a} \sigma(a) \cdot a \|_{\infty}$. This problem has been extensively studied in the static/offline setting. In this paper we initiate its study in the fully-dynamic setting with recourse: the algorithm sees a stream of T insertions and deletions of vectors, and at each time must maintain a low-discrepancy signing, while also minimizing the amortized recourse (the number of times any vector changes its sign) per update. For general vectors, we show algorithms which almost match Spencer's $O(\sqrt{n})$ offline discrepancy bound, with ${O}(n\cdot poly\!\log T)$ amortized recourse per update. The crucial idea is to compute a basic feasible solution to the linear relaxation in a distributed and recursive manner, which helps find a low-discrepancy signing. To bound recourse we argue that only a small part of the instance needs to be re-computed at each update. Since vector balancing has also been greatly studied for sparse vectors, we then give algorithms for low-discrepancy edge orientation, where we dynamically maintain signings for 2-sparse vectors. Alternatively, this can be seen as orienting a dynamic set of edges of an n-vertex graph to minimize the absolute difference between in- and out-degrees at any vertex. We present a deterministic algorithm with $O(poly\!\log n)$ discrepancy and $O(poly\!\log n)$ amortized recourse. The core ideas are to dynamically maintain an expander-decomposition with low recourse and then to show that, as the expanders change over time, a natural local-search algorithm converges quickly (i.e., with low recourse) to a low-discrepancy solution. We also give strong lower bounds for local-search discrepancy minimization algorithms.


title: Failure of standard density functional theory to describe the phase behavior of a fluid of hard right isosceles triangles

abstract: A fluid of hard right isosceles triangles was studied using an extension of Scaled-Particle Density-Functional Theory which includes the exact third virial coefficient. We show that the only orientationally ordered stable liquid-crystal phase predicted by the theory is the uniaxial nematic phase, in agreement with the second-order virial theory. By contrast, Monte Carlo simulations predict exotic liquid-crystal phases exhibiting tetratic and octatic correlations, with orientational distribution functions having four and eight equivalent peaks, respectively. This demonstrates the failure of the standard Density Functional Theory based on two and three-body correlations to describe high-symmetry orientational phases in two-dimensional hard right-triangle fluids, and points to the necessity to reformulate the theory to take into account high-order body correlations and ultimately particle self-assembling and clustering effects. This avenue may represent a great challenge for future research, and we discuss some fundamental ideas to construct a modified version of Density-Functional Theory to account for these clustering effects.


title: Self-Normalized Importance Sampling for Neural Language Modeling

abstract: To mitigate the problem of having to traverse over the full vocabulary in the softmax normalization of a neural language model, sampling-based training criteria are proposed and investigated in the context of large vocabulary word-based neural language models. These training criteria typically enjoy the benefit of faster training and testing, at a cost of slightly degraded performance in terms of perplexity and almost no visible drop in word error rate. While noise contrastive estimation is one of the most popular choices, recently we show that other sampling-based criteria can also perform well, as long as an extra correction step is done, where the intended class posterior probability is recovered from the raw model outputs. In this work, we propose self-normalized importance sampling. Compared to our previous work, the criteria considered in this work are self-normalized and there is no need to further conduct a correction step. Through self-normalized language model training as well as lattice rescoring experiments, we show that our proposed self-normalized importance sampling is competitive in both research-oriented and production-oriented automatic speech recognition tasks.


title: On the number of cycles in commutators of random permutations

abstract: We present general links between statistics of non-Hermitian random matrices and the distribution of the number of cycles of some specific random permutations. In particular, we derive explicit formulas for the generating functions of the number of cycles in the commutator $[\sigma,\tau] = \sigma \tau \sigma^{-1} \tau^{-1}$ where $\sigma$ is uniformly distributed, and $\tau$ is either one cycle, the product of two cycles of same size, or the product of many transpositions.


title: Implicit SVD for Graph Representation Learning

abstract: Recent improvements in the performance of state-of-the-art (SOTA) methods for Graph Representational Learning (GRL) have come at the cost of significant computational resource requirements for training, e.g., for calculating gradients via backprop over many data epochs. Meanwhile, Singular Value Decomposition (SVD) can find closed-form solutions to convex problems, using merely a handful of epochs. In this paper, we make GRL more computationally tractable for those with modest hardware. We design a framework that computes SVD of \textit{implicitly} defined matrices, and apply this framework to several GRL tasks. For each task, we derive linear approximation of a SOTA model, where we design (expensive-to-store) matrix $\mathbf{M}$ and train the model, in closed-form, via SVD of $\mathbf{M}$, without calculating entries of $\mathbf{M}$. By converging to a unique point in one step, and without calculating gradients, our models show competitive empirical test performance over various graphs such as article citation and biological interaction networks. More importantly, SVD can initialize a deeper model, that is architected to be non-linear almost everywhere, though behaves linearly when its parameters reside on a hyperplane, onto which SVD initializes. The deeper model can then be fine-tuned within only a few epochs. Overall, our procedure trains hundreds of times faster than state-of-the-art methods, while competing on empirical test performance. We open-source our implementation at: https://github.com/samihaija/isvd


title: Modeling a Three-Stage SQUID System in Space with the First Micro-X Sounding Rocket Flight

abstract: The Micro-X sounding rocket is a NASA funded X-ray telescope payload that completed its first flight on July 22, 2018. This event marked the first operation of Transition Edge Sensors (TESs) and their SQUID-based multiplexing readout system in space. Unfortunately, due to an ACS pointing failure, the rocket was spinning during its five minute observation period and no scientific data was collected. However, data collected from the internal calibration source marked a partial success for the payload and offers a unique opportunity to study the response of TESs and SQUIDs in space. Of particular interest is the magnetic field response of the NIST MUX06a SQUID readout system to tumbling through Earth's magnetic field. We present a model to explain the baseline response of the SQUIDs, which lead to a subset of pixels failing to "lock" for the full observational period. Future flights of the Micro-X rocket will include the NIST MUX18b SQUID system with dramatically reduced magnetic susceptibility.


title: Proper Scoring Rules, Gradients, Divergences, and Entropies for Paths and Time Series

abstract: Many forecasts consist not of point predictions but concern the evolution of quantities. For example, a central bank might predict the interest rates during the next quarter, an epidemiologist might predict trajectories of infection rates, a clinician might predict the behaviour of medical markers over the next day, etc. The situation is further complicated since these forecasts sometimes only concern the approximate "shape of the future evolution" or "order of events". Formally, such forecasts can be seen as probability measures on spaces of equivalence classes of paths modulo time-parametrization. We leverage the statistical framework of proper scoring rules with classical mathematical results to derive a principled approach to decision making with such forecasts. In particular, we introduce notions of gradients, entropy, and divergence that are tailor-made to respect the underlying non-Euclidean structure.


title: Gradient-push algorithm for distributed optimization with event-triggered communications

abstract: Decentralized optimization problems consist of multiple agents connected by a network. The agents have each local cost function, and the goal is to minimize the sum of the functions cooperatively. It requires the agents communicate with each other, and reducing the cost for communication is desired for a communication-limited environment. In this work, we propose a gradient-push algorithm involving event-triggered communication on directed network. Each agent sends its state information to its neighbors only when the difference between the latest sent state and the current state is larger than a threshold. The convergence of the algorithm is established under a decay and a summability condition on a stepsize and a triggering threshold. Numerical experiments are presented to support the effectiveness and the convergence results of the algorithm.


title: Unsupervised Noise Adaptive Speech Enhancement by Discriminator-Constrained Optimal Transport

abstract: This paper presents a novel discriminator-constrained optimal transport network (DOTN) that performs unsupervised domain adaptation for speech enhancement (SE), which is an essential regression task in speech processing. The DOTN aims to estimate clean references of noisy speech in a target domain, by exploiting the knowledge available from the source domain. The domain shift between training and testing data has been reported to be an obstacle to learning problems in diverse fields. Although rich literature exists on unsupervised domain adaptation for classification, the methods proposed, especially in regressions, remain scarce and often depend on additional information regarding the input data. The proposed DOTN approach tactically fuses the optimal transport (OT) theory from mathematical analysis with generative adversarial frameworks, to help evaluate continuous labels in the target domain. The experimental results on two SE tasks demonstrate that by extending the classical OT formulation, our proposed DOTN outperforms previous adversarial domain adaptation frameworks in a purely unsupervised manner.


title: Faint intermediate luminosity optical transients (ILOTs) from engulfing exoplanets on the Hertzsprung gap

abstract: We follow the evolution of four observed exoplanets to the time when the respective parent star of each planet evolves off the main sequence and engulfs its planet to start a common envelope evolution (CEE), concluding that in each case this process powers an intermediate luminosity optical transient (ILOT; luminous red nova). We characterise the final thousands of days of the orbital decay towards a CEE and determine the properties of the star at the onset of the CEE. We scale the properties of the ILOT V1309 Scorpii to the properties of a planet that enters a CEE inside a star on and near the Hertzsprung gap to estimate the duration and luminosity of the expected ILOT. Based on these we estimate that for a planet of Jupiter mass the ILOT will last for several days and reach a luminosity of several thousand solar luminosity. This type of ILOTs are less luminous than classical novae. Because of the small amount of expected dust and the small amount of energy that an accretion process onto the planet can release, such ILOTs can teach us on the merger at the onset of CEE of stellar companions. Our study adds to the variety of ILOTs that planets can power as they interact with a more massive companion.


title: Multi-agent Reinforcement Learning for Cooperative Lane Changing of Connected and Autonomous Vehicles in Mixed Traffic

abstract: Autonomous driving has attracted significant research interests in the past two decades as it offers many potential benefits, including releasing drivers from exhausting driving and mitigating traffic congestion, among others. Despite promising progress, lane-changing remains a great challenge for autonomous vehicles (AV), especially in mixed and dynamic traffic scenarios. Recently, reinforcement learning (RL), a powerful data-driven control method, has been widely explored for lane-changing decision makings in AVs with encouraging results demonstrated. However, the majority of those studies are focused on a single-vehicle setting, and lane-changing in the context of multiple AVs coexisting with human-driven vehicles (HDVs) have received scarce attention. In this paper, we formulate the lane-changing decision making of multiple AVs in a mixed-traffic highway environment as a multi-agent reinforcement learning (MARL) problem, where each AV makes lane-changing decisions based on the motions of both neighboring AVs and HDVs. Specifically, a multi-agent advantage actor-critic network (MA2C) is developed with a novel local reward design and a parameter sharing scheme. In particular, a multi-objective reward function is proposed to incorporate fuel efficiency, driving comfort, and safety of autonomous driving. Comprehensive experimental results, conducted under three different traffic densities and various levels of human driver aggressiveness, show that our proposed MARL framework consistently outperforms several state-of-the-art benchmarks in terms of efficiency, safety and driver comfort.


title: Lower Bounds on Volumes of Hyperbolic 3-Manifolds via Decomposition

abstract: In a variety of settings we provide a method for decomposing a 3-manifold $M$ into pieces. When the pieces have the appropriate type of hyperbolicity, then the manifold $M$ is hyperbolic and its volume is bounded below by the sum of the appropriately defined hyperbolic volumes of the pieces. A variety of examples of appropriately hyperbolic pieces and volumes are provided, with many examples from link complements in the 3-sphere.


title: An Algebraic and Microlocal Approach to the Stochastic Non-linear Schr\"odinger Equation

abstract: In a recent work [DDRZ20], it has been developed a novel framework aimed at studying at a perturbative level a large class of non-linear, scalar, real, stochastic PDEs and inspired by the algebraic approach to quantum field theory. The main advantage is the possibility of computing the expectation value and the correlation functions of the underlying solutions accounting for renormalization intrinsically and without resorting to any specific regularization scheme. In this work we prove that it is possible to extend the range of applicability of this framework to cover also the stochastic non-linear Schroedinger equation in which randomness is codified by an additive, Gaussian, complex white noise.


title: Carbon nanotubes for the optical far-field readout of processes that are mediated by plasmonic near-fields

abstract: As science progresses at the nanoscopic level, it becomes more and more important to comprehend the interactions taking place at the nanoscale, where optical near-fields play a key role. Their phenomenology differs significantly from the propagative light we experience at the macroscopic level. This is particularly important in applications such as surface-enhanced spectroscopies for single-molecule detection, where often the optimization of the plasmonic structures and surfaces relies on far-field characterizations. The processes dominating in the far-field picture, though, are not the same dominating in the near-field. To highlight this, we resort to very simple metallic systems: Isolated gold nanorods in solution. We show how single-walled nanotubes can be exploited to read out processes occurring at the near-field level around metallic nanoparticles and make the information accessible in the far-field region. This is implemented by monitoring the spectral profile of the enhancement of the photoluminescence and Raman signal of the nanotubes for several excitation wavelengths. Through this excitation-resolved study, we show that the far-field optical read-out detects the transversal and longitudinal dipolar plasmonic oscillations of gold nanorods, whereas the near-field read-out through the nanotubes reveals other mechanisms to dominate. The spectral position of the maximum enhancement of the optical near-field mediated signals are located elsewhere than the far-field bands. This dichotomy between near-field and far-field response should be taken into account when optimizing plasmonic nanostructures for applications such as surface-enhanced spectroscopies.


title: Multisensor acoustic emission data set reflecting supervised untightening of bolts in a jointed vibrating structure

abstract: The data set presented in this work, called ORION-AE, is made of raw AE data streams collected by three different AE sensors and a laser vibrometer during five campaigns of measurements by varying the tightening conditions of two bolted plates submitted to harmonic vibration tests. With seven different operating conditions, this data set was designed to challenge supervised and unsupervised machine/deep learning as well as signal processing methods which are developed for material characterization or Structural Health Monitoring (SHM). One motivation of this work was to create a common benchmark for comparing data-driven methods dedicated to AE data interpretation. The data set is made of time-series collected during an experiment designed to reproduce the loosening phenomenon observed in aeronautics, automotive or civil engineering structures where parts are assembled together by means of bolted joints. Monitoring loosening in jointed structures during operation remains challenging because contact and friction in bolted joints induce a nonlinear stochastic behavior. ORION-AE data set is available on a shared directory for Research Data, Harvard Dataverse (https://doi.org/10.7910/DVN/FBRDU0). A Matlab code is provided to extract the data stream from each sensor. This article is a preprint, the published version is available at https://doi.org/10.3390/data7030031.


title: An Online Multi-Index Approach to Human Ergonomics Assessment in the Workplace

abstract: Work-related musculoskeletal disorders (WMSDs) remain one of the major occupational safety and health problems in the European Union nowadays. Thus, continuous tracking of workers' exposure to the factors that may contribute to their development is paramount. This paper introduces an online approach to monitor kinematic and dynamic quantities on the workers, providing on the spot an estimate of the physical load required in their daily jobs. A set of ergonomic indexes is defined to account for multiple potential contributors to WMSDs, also giving importance to the subject-specific requirements of the workers. To evaluate the proposed framework, a thorough experimental analysis was conducted on twelve human subjects considering tasks that represent typical working activities in the manufacturing sector. For each task, the ergonomic indexes that better explain the underlying physical load were identified, following a statistical analysis, and supported by the outcome of a surface electromyography (sEMG) analysis. A comparison was also made with a well-recognised and standard tool to evaluate human ergonomics in the workplace, to highlight the benefits introduced by the proposed framework. Results demonstrate the high potential of the proposed framework in identifying the physical risk factors, and therefore to adopt preventive measures. Another equally important contribution of this study is the creation of a comprehensive database on human kinodynamic measurements, which hosts multiple sensory data of healthy subjects performing typical industrial tasks.


title: How player and opponent personalities influence cooperative gameplay

abstract: Research has shown that digital game players often feel engagement and rapport with a game hero or character when they can channel their own ambitions and goals through the hero's journey in the game world; in essence, they feel a sense of accomplishment and fulfilment whenever they put the game mechanics to use to help the hero reach a positive ending to the game quests. In the case of cooperative gameplay, rapport also has to do with their perception of their peers' skills, gameplay style and behaviour within the game. In this paper, we describe an experiment to identify whether matching players with different personalities, as characterized by the OCEAN or Big-5 personality model, can influence their player experience with a custom-made, cooperative game.


title: Measurement catastrophe and ballistic spread of charge density with vanishing current

abstract: One of the features of many-body quantum systems with Hilbert-space fragmentation is the existence of stationary states manifesting quantum jamming. It was recently shown that these are "states with memory", in which, e.g., measuring a localised observable has everlasting macroscopic effects. We study such a measurement catastrophe with an example that stands out for its clarity. We show in particular that at late times the expectation value of a charge density becomes a nontrivial function of the ratio between distance and time notwithstanding the corresponding current approaching zero.


title: Integrating psychotherapy practices and gamified elements in novel game mechanics for stress relief

abstract: We explore novel game mechanics and techniques in the domain of gamified and game-based mobile mental health applications. By combining modern game design elements with techniques applied by practitioners (e.g., therapists) and known mechanics used in relevant games, we developed an integrated mobile game. Playtesting with a group of individuals showed a positive response towards the study's claims and a promising direction for further research.


title: Simulating High-Dimensional Multivariate Data using the bigsimr R Package

abstract: It is critical to accurately simulate data when employing Monte Carlo techniques and evaluating statistical methodology. Measurements are often correlated and high dimensional in this era of big data, such as data obtained in high-throughput biomedical experiments. Due to the computational complexity and a lack of user-friendly software available to simulate these massive multivariate constructions, researchers resort to simulation designs that posit independence or perform arbitrary data transformations. To close this gap, we developed the Bigsimr Julia package with R and Python interfaces. This paper focuses on the R interface. These packages empower high-dimensional random vector simulation with arbitrary marginal distributions and dependency via a Pearson, Spearman, or Kendall correlation matrix. bigsimr contains high-performance features, including multi-core and graphical-processing-unit-accelerated algorithms to estimate correlation and compute the nearest correlation matrix. Monte Carlo studies quantify the accuracy and scalability of our approach, up to $d=10,000$. We describe example workflows and apply to a high-dimensional data set -- RNA-sequencing data obtained from breast cancer tumor samples.


title: Stationary Behavior of Constant Stepsize SGD Type Algorithms: An Asymptotic Characterization

abstract: Stochastic approximation (SA) and stochastic gradient descent (SGD) algorithms are work-horses for modern machine learning algorithms. Their constant stepsize variants are preferred in practice due to fast convergence behavior. However, constant step stochastic iterative algorithms do not converge asymptotically to the optimal solution, but instead have a stationary distribution, which in general cannot be analytically characterized. In this work, we study the asymptotic behavior of the appropriately scaled stationary distribution, in the limit when the constant stepsize goes to zero. Specifically, we consider the following three settings: (1) SGD algorithms with smooth and strongly convex objective, (2) linear SA algorithms involving a Hurwitz matrix, and (3) nonlinear SA algorithms involving a contractive operator. When the iterate is scaled by $1/\sqrt{\alpha}$, where $\alpha$ is the constant stepsize, we show that the limiting scaled stationary distribution is a solution of an integral equation. Under a uniqueness assumption (which can be removed in certain settings) on this equation, we further characterize the limiting distribution as a Gaussian distribution whose covariance matrix is the unique solution of a suitable Lyapunov equation. For SA algorithms beyond these cases, our numerical experiments suggest that unlike central limit theorem type results: (1) the scaling factor need not be $1/\sqrt{\alpha}$, and (2) the limiting distribution need not be Gaussian. Based on the numerical study, we come up with a formula to determine the right scaling factor, and make insightful connection to the Euler-Maruyama discretization scheme for approximating stochastic differential equations.


title: A Systematic Look at the Temperature Gradient Contribution to the Dayside Magnetopause Current

abstract: Magnetopause diamagnetic currents arise from density and temperature driven pressure gradients across the boundary layer. While theoretically recognized, the temperature contributions to the magnetopause current system have not yet been systematically studied. To bridge this gap, we used a database of Magnetospheric Multiscale (MMS) magnetopause crossings to analyze diamagnetic current densities and their contributions across the dayside and flank magnetopause. Our results indicate that the ion temperature gradient component makes up to 37% of the ion diamagnetic current density along the magnetopause and typically opposes the classical Chapman-Ferraro current direction, interfering destructively with the density gradient component, thus lowering the total diamagnetic current density. This effect is most pronounced on the flank magnetopause. The electron diamagnetic current was found to be 5 to 14 times weaker than the ion diamagnetic current on average.


title: Improving the Solar Wind Density Model Used in Processing of Spacecraft Ranging Observations

abstract: Solar plasma as a cause of radio signal delay has been playing an important role in solar and planetary science. Early experiments studying the distribution of electrons near the Sun from spacecraft ranging measurements were designed so that the radio signal was passing close to the Sun. At present, processing of spacecraft tracking observations serves a different goal: precise (at meter level) determination of orbits of planets, most importantly Mars. Solar plasma adds a time-varying delay to those observations, which is, in this case, unwanted and must be subtracted prior to putting the data into planetary solution. Present planetary ephemeris calculate the delay assuming symmetric stationary power-law model of solar plasma. The present work, based on a custom variant of the EPM lunar-planetary ephemeris, raises the question of accuracy and correctness of that assumption and examines alternative models based on in situ data provided by OMNI and on the ENLIL numerical model of solar wind.

categories: physics.space-ph astro-ph.EP astro-ph.IM astro-ph.SR


title: Towards an Efficient Voice Identification Using Wav2Vec2.0 and HuBERT Based on the Quran Reciters Dataset

abstract: Current authentication and trusted systems depend on classical and biometric methods to recognize or authorize users. Such methods include audio speech recognitions, eye, and finger signatures. Recent tools utilize deep learning and transformers to achieve better results. In this paper, we develop a deep learning constructed model for Arabic speakers identification by using Wav2Vec2.0 and HuBERT audio representation learning tools. The end-to-end Wav2Vec2.0 paradigm acquires contextualized speech representations learnings by randomly masking a set of feature vectors, and then applies a transformer neural network. We employ an MLP classifier that is able to differentiate between invariant labeled classes. We show several experimental results that safeguard the high accuracy of the proposed model. The experiments ensure that an arbitrary wave signal for a certain speaker can be identified with 98% and 97.1% accuracies in the cases of Wav2Vec2.0 and HuBERT, respectively.


title: Is the Hemispheric Asymmetry of Monthly Sunspot Area an Irregular Process with Long-Term Memory?

abstract: The hemispheric asymmetry of the sunspot cycle is a real feature of the Sun. However, its origin is still not well understood. Here we perform nonlinear time series analysis of the sunspot area (and number) asymmetry to explore its dynamics. By measuring the correlation dimension of the sunspot area asymmetry, we conclude that there is no strange attractor in the data. Further computing Higuchi's dimension, we conclude that the hemispheric asymmetry is largely governed by stochastic noise. However, the behaviour of Hurst exponent reveals that the time series is not completely determined by a memory-less stochastic noise, rather there is a long-term persistence, which can go beyond two solar cycles. Therefore, our study suggests that the hemispheric asymmetry of the sunspot cycle is predominantly originated due to some irregular process in the solar dynamo. The long-term persistence in the solar cycle asymmetry suggests that the solar magnetic field has some memory in the convection zone.


title: Nucleon Axial Form Factor from Domain Wall on HISQ

abstract: The Deep Underground Neutrino Experiment (DUNE) is an upcoming neutrino oscillation experiment that is poised to answer key questions about the nature of neutrinos. Lattice QCD has the ability to make significant impact upon DUNE, beginning with computations of nucleon-neutrino interactions with weak currents. Nucleon amplitudes involving the axial form factor are part of the primary signal measurement process for DUNE, and precise calculations from LQCD can significantly reduce the uncertainty for inputs into Monte Carlo generators. Recent calculations of the nucleon axial charge have demonstrated that sub-percent precision is possible on this vital quantity. In these proceedings, we discuss preliminary results for the CalLat collaboration's calculation of the axial form factor of the nucleon. These computations are performed with M\"obius domain wall valence quarks on HISQ sea quark ensembles generated by the MILC and CalLat collaborations. The results use a variety of ensembles including several at physical pion mass.


title: Identification of Fine-Grained Location Mentions in Crisis Tweets

abstract: Identification of fine-grained location mentions in crisis tweets is central in transforming situational awareness information extracted from social media into actionable information. Most prior works have focused on identifying generic locations, without considering their specific types. To facilitate progress on the fine-grained location identification task, we assemble two tweet crisis datasets and manually annotate them with specific location types. The first dataset contains tweets from a mixed set of crisis events, while the second dataset contains tweets from the global COVID-19 pandemic. We investigate the performance of state-of-the-art deep learning models for sequence tagging on these datasets, in both in-domain and cross-domain settings.


title: Sparse grid approximation in weighted Wiener spaces

abstract: We study approximation properties of multivariate periodic functions from weighted Wiener spaces by sparse grids methods constructed with the help of quasi-interpolation operators. The class of such operators includes classical interpolation and sampling operators, Kantorovich-type operators, scaling expansions associated with wavelet constructions, and others. We obtain the rate of convergence of the corresponding sparse grids methods in weighted Wiener norms as well as analogues of the Littlewood-Paley-type characterizations in terms of families of quasi-interpolation operators.


title: Character-level HyperNetworks for Hate Speech Detection

abstract: The massive spread of hate speech, hateful content targeted at specific subpopulations, is a problem of critical social importance. Automated methods of hate speech detection typically employ state-of-the-art deep learning (DL)-based text classifiers-large pretrained neural language models of over 100 million parameters, adapting these models to the task of hate speech detection using relevant labeled datasets. Unfortunately, there are only a few public labeled datasets of limited size that are available for this purpose. We make several contributions with high potential for advancing this state of affairs. We present HyperNetworks for hate speech detection, a special class of DL networks whose weights are regulated by a small-scale auxiliary network. These architectures operate at character-level, as opposed to word or subword-level, and are several orders of magnitude smaller compared to the popular DL classifiers. We further show that training hate detection classifiers using additional large amounts of automatically generated examples is beneficial in general, yet this practice especially boosts the performance of the proposed HyperNetworks. We report the results of extensive experiments, assessing the performance of multiple neural architectures on hate detection using five public datasets. The assessed methods include the pretrained language models of BERT, RoBERTa, ALBERT, MobileBERT and CharBERT, a variant of BERT that incorporates character alongside subword embeddings. In addition to the traditional setup of within-dataset evaluation, we perform cross-dataset evaluation experiments, testing the generalization of the various models in conditions of data shift. Our results show that the proposed HyperNetworks achieve performance that is competitive, and better in some cases, than these pretrained language models, while being smaller by orders of magnitude.


title: Optimal control of families of quantum gates

abstract: Quantum Optimal Control (QOC) enables the realization of accurate operations, such as quantum gates, and support the development of quantum technologies. To date, many QOC frameworks have been developed but those remain only naturally suited to optimize a single targeted operation at a time. We extend this concept to optimal control with a continuous family of targets, and demonstrate that an optimization based on neural networks can find families of time-dependent Hamiltonians that realize desired classes of quantum gates in minimal time.


title: On the parameterized complexity of Compact Set Packing

abstract: The Set Packing problem is, given a collection of sets $\mathcal{S}$ over a ground set $\mathcal{U}$, to find a maximum collection of sets that are pairwise disjoint. The problem is among the most fundamental NP-hard optimization problems that have been studied extensively in various computational regimes. The focus of this work is on parameterized complexity, Parameterized Set Packing (PSP): Given $r \in {\mathbb N}$, is there a collection $ \mathcal{S}' \subseteq \mathcal{S}: |\mathcal{S}'| = r$ such that the sets in $\mathcal{S}'$ are pairwise disjoint? Unfortunately, the problem is not fixed parameter tractable unless $\mathsf{W[1] = FPT}$, and, in fact, an "enumeration" running time of $|\mathcal{S}|^{\Omega(r)}$ is required unless the exponential time hypothesis (ETH) fails. This paper is a quest for tractable instances of Set Packing from parameterized complexity perspectives. We say that the input $(\mathcal{U},\mathcal{S})$ is "compact" if $|\mathcal{U}| = f(r)\cdot\Theta(\textsf{poly}( \log |\mathcal{S}|))$, for some $f(r) \ge r$. In the Compact Set Packing problem, we are given a compact instance of PSP. In this direction, we present a "dichotomy" result of PSP: When $|\mathcal{U}| = f(r)\cdot o(\log |\mathcal{S}|)$, PSP is in $\textsf{FPT}$, while for $|\mathcal{U}| = r\cdot\Theta(\log (|\mathcal{S}|))$, the problem is $W[1]$-hard; moreover, assuming ETH, Compact PSP does not even admit $|\mathcal{S}|^{o(r/\log r)}$ time algorithm. Although certain results in the literature imply hardness of compact versions of related problems such as Set $r$-Covering and Exact $r$-Covering, these constructions fail to extend to Compact PSP. A novel contribution of our work is the identification and construction of a gadget, which we call Compatible Intersecting Set System pair, that is crucial in obtaining the hardness result for Compact PSP.


title: Supporting and Controlling Complex Concurrency in Fault- Tolerant Distributed Systems

abstract: Distributed computing often gives rise to complex concurrent and interacting activities. In some cases several concurrent activities may be working together, i.e. cooperating, to solve a given problem; in other cases, the activities may be independent but needing to share common system resources for which they must compete. Many difficulties and limitations occur in the widely advocated objects and (trans)actions model when it is supposed to support cooperating activities. We have introduced previously the concept of coordinated atomic (CA) actions [Xu et al. 1995]; this paper analyzes and examines the derived objects and CA actions model for constructing fault-tolerant distributed systems and providing unified support for both cooperative and competitive concurrency. Our investigation reveals and clarifies several significant problems that have not previously been studied extensively, including the problem of ensuring consistent access to shared objects from a joint action as opposed to a set of independent actions. Conceptual and implementation-related solutions are proposed and illustrated.


title: ParaFold: Paralleling AlphaFold for Large-Scale Predictions

abstract: AlphaFold predicts protein structures from the amino acid sequence at or near experimental resolution, solving the 50-year-old protein folding challenge, leading to progress by transforming large-scale genomics data into protein structures. AlphaFold will also greatly change the scientific research model from low-throughput to high-throughput manner. The AlphaFold framework is a mixture of two types of workloads: MSA construction based on CPUs and model inference on GPUs. The first CPU stage dominates the overall runtime, taking hours for a single protein due to the large database sizes and I/O bottlenecks. However, GPUs in this CPU stage remain idle, resulting in low GPU utilization and restricting the capacity of large-scale structure predictions. Therefore, we proposed ParaFold, an open-source parallel version of AlphaFold for high throughput protein structure predictions. ParaFold separates the CPU and GPU parts to enable large-scale structure predictions. ParaFold also effectively reduces the CPU and GPU runtime with two optimizations without compromising the quality of prediction results: using multi-threaded parallelism on CPUs and using optimized JAX compilation on GPUs. We evaluated ParaFold with three datasets of different size and protein lengths. We evaluated the accuracy and efficiency of optimizations on CPUs and GPUs, and showed the large-scale prediction capability by running ParaFold inferences of 19,704 small proteins in five hours on one NVIDIA DGX-2. Using the JAX compile optimization, ParaFold attained a 13.8X average speedup over AlphaFold. ParaFold offers a rapid and effective approach for high-throughput structure predictions, leveraging the predictive power by running on supercomputers, with shorter time, and at a lower cost. The development of ParaFold will greatly speed up high-throughput studies and render the protein "structure-omics" feasible.


title: Stochastic diffusion of electrons interacting with whistler-mode waves in the solar wind

abstract: Effects of increasing whistler amplitude and propagation angle are studied through a variational test particle simulation and calculations of the resonance width. While high amplitude and oblique whistlers in typical 1 AU solar wind parameters are capable of forming an isotropic population without any additional processes, anomalous interactions with quasi-parallel whistlers may be essential to the process of halo formation near the Sun. High amplitude and quasi-parallel whistlers can scatter strahl electrons to low velocities (less than the wave phase velocity) to form a halo population, as long as their amplitude is sufficiently high. We also present in detail a careful treatment of the sensitivity to initial conditions based on calculations of the phase space volume, which is necessary for numerical calculations of highly stochastic motion due to resonant interactions with large amplitude waves. Our method ensures that the volume-preserving characteristic of the Boris algorithm is consistently applied for simulations of both stochastic and non-stochastic particle motion.

categories: astro-ph.SR nlin.CD physics.plasm-ph physics.space-ph


title: Driver-Specific Risk Recognition in Interactive Driving Scenarios using Graph Representation

abstract: This paper presents a driver-specific risk recognition framework for autonomous vehicles that can extract inter-vehicle interactions. This extraction is carried out for urban driving scenarios in a driver-cognitive manner to improve the recognition accuracy of risky scenes. First, clustering analysis is applied to the operation data of drivers for learning the subjective assessment of risky scenes of different drivers and generating the corresponding risk label for each scene. Second, the graph representation model (GRM) is adopted to unify and construct the features of dynamic vehicles, inter-vehicle interactions and static traffic markings in real driving scenes into graphs. The driver-specific risk label provides ground truth to capture the risk evaluation criteria of different drivers. In addition, the graph model represents multiple features of the driving scenes. Therefore, the proposed framework can learn the risk-evaluating pattern of driving scenes of different drivers and establish driver-specific risk identifiers. Last, the performance of the proposed framework is evaluated via experiments conducted using real-world urban driving datasets collected by multiple drivers. The results show that the risks and their levels in real driving environments can be accurately recognized by the proposed framework.