title: On the Benefit of Adversarial Training for Monocular Depth Estimation

abstract: In this paper we address the benefit of adding adversarial training to the task of monocular depth estimation. A model can be trained in a self-supervised setting on stereo pairs of images, where depth (disparities) are an intermediate result in a right-to-left image reconstruction pipeline. For the quality of the image reconstruction and disparity prediction, a combination of different losses is used, including L1 image reconstruction losses and left-right disparity smoothness. These are local pixel-wise losses, while depth prediction requires global consistency. Therefore, we extend the self-supervised network to become a Generative Adversarial Network (GAN), by including a discriminator which should tell apart reconstructed (fake) images from real images. We evaluate Vanilla GANs, LSGANs and Wasserstein GANs in combination with different pixel-wise reconstruction losses. Based on extensive experimental evaluation, we conclude that adversarial training is beneficial if and only if the reconstruction loss is not too constrained. Even though adversarial training seems promising because it promotes global consistency, non-adversarial training outperforms (or is on par with) any method trained with a GAN when a constrained reconstruction loss is used in combination with batch normalisation. Based on the insights of our experimental evaluation we obtain state-of-the art monocular depth estimation results by using batch normalisation and different output scales.


title: A semiclassical Thomas-Fermi model to tune the metallicity of electrodes in molecular simulations

abstract: Spurred by the increasing needs in electrochemical energy storage devices, the electrode/electrolyte interface has received a lot of interest in recent years. Molecular dynamics simulations play a proeminent role in this field since they provide a microscopic picture of the mechanisms involved. The current state-of-the-art consists in treating the electrode as a perfect conductor, precluding the possibility to analyze the effect of its metallicity on the interfacial properties. Here we show that the Thomas-Fermi model provides a very convenient framework to account for the screening of the electric field at the interface and differenciating good metals such as gold from imperfect conductors such as graphite. All the interfacial properties are modified by screening within the metal: the capacitance decreases significantly and both the structure and dynamics of the adsorbed electrolyte are affected. The proposed model opens the door for quantitative predictions of the capacitive properties of materials for energy storage.


title: Search for footprints of quantum spacetime in black hole QNM spectrum

abstract: Black hole (BH) perturbation is followed by a ringdown phase which is dominated by quasinormal modes (QNM). These modes may provide key signature of the gravitational waves. The presence of a deformed spacetime structure may distort this signal. In order to account for such effects, we consider a toy model consisting of a noncommutative charged scalar field propagating in a realistic black hole background. We then analyse the corresponding field dynamics by applying the methods of the Hopf algebra deformation by Drinfeld twist. The latter framework is well suited for incorporating deformed symmetries into a study of this kind. As a result, we obtain the BH QNM spectrum that, besides containing the intrinsic information about a black hole that is being analysed, also carry the information about the underlying structure of spacetime.


title: Solubility of nitrogen in methane, ethane, and mixtures of methane and ethane at Titan-like conditions: A molecular dynamics study

abstract: We have studied the temperature dependence of the solubility of nitrogen in methane, ethane, and mixtures of methane and ethane using vapor-liquid equilibrium simulations of binary and ternary mixtures of nitrogen, methane and ethane for a range of temperatures between 90K and 110K at a pressure of 1.5 atm, thermodynamic conditions that may exist on the Saturn's giant moon, Titan. We find that --(i) the solubility of nitrogen in both methane and ethane decreases with increasing temperature; (ii) the solubility of nitrogen in methane is much larger compared to that in ethane at low temperatures, (iii) solubility of nitrogen in a ternary mixture of methane, ethane, and nitrogen increases upon increasing mole-fraction of methane. Our results are in quantitative agreement with the recent experimental measurement of the solubility of nitrogen in methane, ethane, and a mixture of methane and ethane. Furthermore, we find a strong temperature-dependent surface adsorption of nitrogen at the nitrogen-hydrocarbon interface, previously unknown. We have also investigated surface tension of the gas-liquid interface and find that it decreases upon decreasing temperature. Moreover, we find that the interfacial layer of adsorbed nitrogen and ethane show a preferential orientational ordering at the interface.


title: Influence of doping level on Brillouin oscillations in GaAs

abstract: Time-domain Brillouin scattering has proved to be an unique tool for determining depth dependent material properties. Here, we show the influence of doping level in GaAs on Brillouin oscillations. Measurements were performed on intrinsic, n-type and p-type GaAs samples. The results show high sensitivity of the amplitude of Brillouin oscillations to the doping concentration. The theoretical calculations are in a good agreement with the experimental data. This work provides an insight into the specific dopant profiling as a function of depth.


title: Replay Spoofing Countermeasure Using Autoencoder and Siamese Network on ASVspoof 2019 Challenge

abstract: Automatic Speaker Verification (ASV) is the process of identifying a person based on the voice presented to a system. Different synthetic approaches allow spoofing to deceive ASV systems (ASVs), whether using techniques to imitate a voice or recunstruct the features. Attackers try to beat up the ASVs using four general techniques; impersonation, speech synthesis, voice conversion, and replay. The last technique is considered as a common and high potential tool for spoofing purposes since replay attacks are more accessible and require no technical knowledge from adversaries. In this study, we introduce a novel replay spoofing countermeasure for ASVs. Accordingly, we used the Constant Q Cepstral Coefficient (CQCC) features fed into an autoencoder to attain more informative features and to consider the noise information of spoofed utterances for discrimination purpose. Finally, different configurations of the Siamese network were used for the first time in this context for classification. The experiments performed on ASVspoof challenge 2019 dataset using Equal Error Rate (EER) and Tandem Detection Cost Function (t-DCF) as evaluation metrics show that the proposed system improved the results over the baseline by 10.73% and 0.2344 in terms of EER and t-DCF, respectively.


title: Intelligent-Unrolling: Exploiting Regular Patterns in Irregular Applications

abstract: Modern optimizing compilers are able to exploit memory access or computation patterns to generate vectorization codes. However, such patterns in irregular applications are unknown until runtime due to the input dependence. Thus, either compiler's static optimization or profile-guided optimization based on specific inputs cannot predict the patterns for any common input, which leads to suboptimal code generation. To address this challenge, we develop Intelligent-Unroll, a framework to automatically optimize irregular applications with vectorization. Intelligent-Unroll allows the users to depict the computation task using \textit{code seed} with the memory access and computation patterns represented in \textit{feature table} and \textit{information-code tree}, and generates highly efficient codes. Furthermore, Intelligent-Unroll employs several novel optimization techniques to optimize reduction operations and gather/scatter instructions. We evaluate Intelligent-Unroll with sparse matrix-vector multiplication (SpMV) and graph applications. Experimental results show that Intelligent-Unroll is able to generate more efficient vectorization codes compared to the state-of-the-art implementations.


title: Jump balls, rating falls, and elite status: A sensitivity analysis of three quarterback rating statistics

abstract: Quarterback performance can be difficult to rank, and much effort has been spent in creating new rating systems. However, the input statistics for such ratings are subject to randomness and factors outside the quarterback's control. To investigate this variance, we perform a sensitivity analysis of three quarterback rating statistics: the Traditional 1971 rating by Smith, the Burke, and the Wages of Wins ratings. The comparisons are made at the team level for the 32 NFL teams from 2002-2015, thus giving each case an even 16 games. We compute quarterback ratings for each offense with 1-5 additional touchdowns, 1-5 fewer interceptions, 1-5 additional sacks, and a 1-5 percent increase in the passing completion rate. Our sensitivity analysis provides insight into whether an elite passing team could seem mediocre or vice versa based on random outcomes. The results indicate that the Traditional rating is the most sensitive statistic with respect to touchdowns, interceptions, and completions, whereas the Burke rating is most sensitive to sacks. The analysis suggests that team passing offense rankings are highly sensitive to aspects of football that are out of the quarterback's hands (e.g., deflected passes that lead to interceptions). Thus, on the margins, we show arguments about whether a specific quarterback has entered the elite or remains mediocre are irrelevant.


title: Sequential image processing methods for improving semantic video segmentation algorithms

abstract: Recently, semantic video segmentation gained high attention especially for supporting autonomous driving systems. Deep learning methods made it possible to implement real time segmentation and object identification algorithms on videos. However, most of the available approaches process each video frame independently disregarding their sequential relation in time. Therefore their results suddenly miss some of the object segments in some of the frames even if they were detected properly in the earlier frames. Herein we propose two sequential probabilistic video frame analysis approaches to improve the segmentation performance of the existing algorithms. Our experiments show that using the information of the past frames we increase the performance and consistency of the state of the art algorithms.


title: E2-Train: Training State-of-the-art CNNs with Over 80% Energy Savings

abstract: Convolutional neural networks (CNNs) have been increasingly deployed to edge devices. Hence, many efforts have been made towards efficient CNN inference in resource-constrained platforms. This paper attempts to explore an orthogonal direction: how to conduct more energy-efficient training of CNNs, so as to enable on-device training. We strive to reduce the energy cost during training, by dropping unnecessary computations from three complementary levels: stochastic mini-batch dropping on the data level; selective layer update on the model level; and sign prediction for low-cost, low-precision back-propagation, on the algorithm level. Extensive simulations and ablation studies, with real energy measurements from an FPGA board, confirm the superiority of our proposed strategies and demonstrate remarkable energy savings for training. For example, when training ResNet-74 on CIFAR-10, we achieve aggressive energy savings of >90% and >60%, while incurring a top-1 accuracy loss of only about 2% and 1.2%, respectively. When training ResNet-110 on CIFAR-100, an over 84% training energy saving is achieved without degrading inference accuracy.


title: Tuning in to the radio environment of HD189733b

abstract: The hot Jupiter HD189733b is expected to be a source of strong radio emission, due to its close proximity to its magnetically active host star. Here, we model the stellar wind of its host star, based on reconstructed surface stellar magnetic field maps. We use the local stellar wind properties at the planetary orbit obtained from our models to compute the expected radio emission from the planet. Our findings show that the planet emits with a peak flux density within the detection capabilities of LOFAR. However, due to absorption by the stellar wind itself, this emission may be attenuated significantly. We show that the best time to observe the system is when the planet is near primary transit of the host star, as the attenuation from the stellar wind is lowest in this region.


title: Admiring the Great Mountain: A Celebration Special Issue in Honor of Stephen Grossbergs 80th Birthday

abstract: This editorial summarizes selected key contributions of Prof. Stephen Grossberg and describes the papers in this 80th birthday special issue in his honor. His productivity, creativity, and vision would each be enough to mark a scientist of the first caliber. In combination, they have resulted in contributions that have changed the entire discipline of neural networks. Grossberg has been tremendously influential in engineering, dynamical systems, and artificial intelligence as well. Indeed, he has been one of the most important mentors and role models in my career, and has done so with extraordinary generosity and encouragement. All authors in this special issue have taken great pleasure in hereby commemorating his extraordinary career and contributions.


title: Equipartitions with Wedges and Cones

abstract: A famous result about mass partitions is the so called \emph{Ham-Sandwich theorem}. It states that any $d$ mass distributions in $\mathbb{R}^d$ can be simultaneously bisected by a single hyperplane. In this work, we study two related questions. The first one is how many masses we can simultaneously partition with a $k$-fan, that is, $k$ half-hyperplanes in $\mathbb{R}^d$, emanating from a common $(d-2)$-dimensional apex. This question was extensively studied in the plane, but in higher dimensions the only known results are for the case where $k$ is an odd prime. We extend these results to a larger family of values of $k$. We further present a new result for $k=2$, which generalizes to cones. The second question considers bisections with double wedges or, equivalently, Ham-Sandwich cuts after projective transformations. Here we prove that given $d$ families of $d+1$ point sets each, there is always a projective transformation such that after the transformation, each family has a Ham-Sandwich cut. We further prove a result on partitions with parallel hyperplanes after a projective transformation.


title: Bondi-Hoyle-Lyttleton accretion by binary stars

abstract: Binary stars often move through an ambient medium from which they accrete material and angular momentum, as in triple-star systems, star-forming clouds, young globular clusters and in the centres of galaxies. A binary form of Bondi-Hoyle-Lyttleton accretion results whereby the accretion rate depends on the binary properties: the stellar masses and separation, and the relative wind speed. We present the results of simulations performed with the hydrodynamic code GANDALF, to determine the mass accretion rates over a range of binary separations, inclinations and mass ratios. When the binary separation is short, the binary system accretes like a single star, while accretion onto stars in wide binaries is barely affected by their companion. We investigate intermediate-separation systems in some detail, finding that as the binary separation is increased, accretion rates smoothly decrease from the rate equal to that of a single star to the rate expected from two isolated stars. The form of this decrease depends on the relative centre-of-mass velocity of the binary and the gas, with faster-moving binaries showing a shallower decrease. Accretion rates vary little with orbital inclination, except when the orbit is side-on and the stars pass through each others' wakes. The specific angular momentum accretion rate also depends on the inclination but is never sufficient to prevent the binary orbit from contracting. Our results may be applied to accretion onto protostars, pollution of stars in globular and nuclear clusters, and wind mass-transfer in multiple stellar systems.


title: Wettability patterning in microfluidic devices using thermally-enhanced hydrophobic recovery of PDMS

abstract: Spatial control of wettability is key to many applications of microfluidic devices, ranging from double emulsion generation to localized cell adhesion. A number of techniques, often based on masking, have been developed to produce spatially-resolved wettability patterns at the surface of poly(dimethylsiloxane) (PDMS) elastomers. A major impediment they face is the natural hydrophobic recovery of PDMS: hydrophilized PDMS surfaces tend to return to hydrophobicity with time, mainly because of diffusion of low molecular weight silicone species to the surface. Instead of trying to avoid this phenomenon, we propose in this work to take advantage of hydrophobic recovery to modulate spatially the surface wettability of PDMS. Because temperature speeds up the rate of hydrophobic recovery, we show that space-resolved hydrophobic patterns can be produced by locally heating a plasma-hydrophilized PDMS surface with microresistors. Importantly, local wettability is quantified in microchannels using a fluorescent probe. This "thermo-patterning" technique provides a simple route to in situ wettability patterning in closed PDMS chips, without requiring further surface chemistry.


title: Simultaneous Interference-Data Transmission for Secret Key Generation in Distributed IoT Sensor Networks

abstract: Internet of Things (IoT) networks for smart sensor nodes in the next generation of smart wireless sensing systems require a distributed security scheme to prevent the passive (eavesdropping) or active (jamming and interference) attacks from untrusted sensor nodes. This paper concerns advancing the security of the IoT system to address their vulnerability to being attacked or compromised by the advancement of future supercomputers. In this work, a novel embedded architecture has been designed and implemented for a distributed IoT network that utilizes a master-slave full-duplex communication to exchange the random and continuous modulated phase shift as the secret key to be used in higher-layer encryptions.


title: Real-time Timepix3 data clustering, visualization and classification with a new Clusterer framework

abstract: With the next-generation Timepix3 hybrid pixel detector, new possibilities and challenges have arisen. The Timepix3 segments active sensor area of 1.98 $cm^2$ into a square matrix of 256 x 256 pixels. In each pixel, the Time of Arrival (ToA, with a time binning of 1.56 $ns$) and Time over Threshold (ToT, energy) are measured simultaneously in a data-driven, i.e. self-triggered, read-out scheme. This contribution presents a framework for data acquisition, real-time clustering, visualization, classification and data saving. All of these tasks can be performed online, directly from multiple readouts through UDP protocol. Clusters are reconstructed on a pixel-by-pixel decision from the stream of not-necessarily chronologically sorted pixel data. To achieve quick spatial pixel-to-cluster matching, non-trivial data structures (quadtree) are utilized. Furthermore, parallelism (i.e multi-threaded architecture) is used to further improve the performance of the framework. Such real-time clustering offers the advantages of online filtering and classification of events. Versatility of the software is ensured by supporting all major operating systems (macOS, Windows and Linux) with both graphical and command-line interfaces. The performance of the real-time clustering and applied filtration methods are demonstrated using data from the Timepix3 network installed in the ATLAS and MoEDAL experiments at CERN.


title: Roles of repertoire diversity in robustness of humoral immune response

abstract: The adaptive immune system relies on diversity of its repertoire of receptors to protect the organism from a great variety of pathogens. Since the initial repertoire is the result of random gene rearrangement, binding of receptors is not limited to pathogen-associated antigens but also includes self antigens. There is a fine balance between having a diverse repertoire, protecting from many different pathogens, and yet reducing its self-reactivity as far as possible to avoid damage to self. In the ageing immune system this balance is altered, manifesting in reduced specificity of response to pathogens or vaccination on a background of higher self-reactivity. To answer the question whether age-related changes of repertoire in the diversity and self/non-self affinity balance of antibodies could explain the reduced efficacy of the humoral response in older people, we construct a minimal mathematical model of the humoral immune response. The principle of least damage allows us, for a given repertoire of antibodies, to resolve a tension between the necessity to neutralise target antigens as quickly as possible and the requirement to limit the damage to self antigens leading to an optimal dynamics of immune response. The model predicts slowing down of immune response for repertoires with reduced diversity and increased self-reactivity.


title: On distance covariance in metric and Hilbert spaces

abstract: Distance covariance is a measure of dependence between two random variables that take values in two, in general different, metric spaces, see Sz\'ekely, Rizzo and Bakirov (2007) and Lyons (2013). It is known that the distance covariance, and its generalization $\alpha$-distance covariance, can be defined in several different ways that are equivalent under some moment conditions. The present paper considers four such definitions and find minimal moment conditions for each of them, together with some partial results when these conditions are not satisfied. The paper also studies the special case when the variables are Hilbert space valued, and shows under weak moment conditions that two such variables are independent if and only if their ($\alpha$-)distance covariance is 0; this extends results by Lyons (2013) and Dehling et al. (2018+). The proof uses a new definition of distance covariance in the Hilbert space case, generalizing the definition for Euclidean spaces using characteristic functions by Sz\'ekely, Rizzo and Bakirov (2007).


title: Crystal-based intensive gamma-ray light sources

abstract: We discuss design and practical realization of novel gamma-ray Crystal-based Light Sources (CLS) that can be constructed through exposure of oriented crystals (linear, bent, periodically bent) to beams of ultrarelativistic charged particles. In an exemplary case study, we estimate brilliance of radiation emitted in a Crystalline Undulator (CU) LS by available positron beams. Intensity of CU radiation in the photon energy range $10^0-10^1$ MeV, which is inaccessible to conventional synchrotrons, undulators and XFELs, greatly exceeds that of laser-Compton scattering LSs and can be higher than predicted in the Gamma Factory proposal to CERN. Brilliance of CU-LSs can be boosted by up to 8 orders of magnitude through the process of superradiance by a pre-bunched beam. Construction of novel CLSs is a challenging task which constitutes a highly interdisciplinary field entangling a broad range of correlated activities. CLSs provide a low-cost alternative to conventional LSs and have enormous number of applications.


title: On the supersymmetric XXX spin chains associated to $\mathfrak{gl}_{1|1}$

abstract: We study the $\mathfrak{gl}_{1|1}$ supersymmetric XXX spin chains. We give an explicit description of the algebra of Hamiltonians acting on any cyclic tensor products of polynomial evaluation $\mathfrak{gl}_{1|1}$ Yangian modules. It follows that there exists a bijection between common eigenvectors (up to proportionality) of the algebra of Hamiltonians and monic divisors of an explicit polynomial written in terms of the Drinfeld polynomials. In particular our result implies that each common eigenspace of the algebra of Hamiltonians has dimension one. We also give dimensions of the generalized eigenspaces. We show that when the tensor product is irreducible, then all eigenvectors can be constructed using Bethe ansatz. We express the transfer matrices associated to symmetrizers and anti-symmetrizers of vector representations in terms of the first transfer matrix and the center of the Yangian.


title: A 2000-year record of arsenic variability from a South Pole ice core

abstract: Using a South Pole ice core covering the last ~2050 years we present a high-resolution (4-27 samples/year), continuous record of natural and anthropogenic arsenic (As) deposition. Our results show that volcanic emissions (notably from Mt. Erebus) are a significant source of As for South Pole, and human activities potentially contributed to As deposition as early as 225 C.E. The most significant anthropogenic source As enrichment in the record, starting in 1975 C.E., is linked to increased copper production in Chile and at least partially to coal combustion from throughout the Southern Hemisphere. East and West Antarctic ice core delivered deposition differences are a result of differences in the atmospheric circulation patterns that transport As to these regions.


title: Invariant Galerkin Ansatz Spaces and Davison-Maki Methods for the Numerical Solution of Differential Riccati Equations

abstract: The differential Riccati equation appears in different fields of applied mathematics like control and system theory. Recently Galerkin methods based on Krylov subspaces were developed for the autonomous differential Riccati equation. These methods overcome the prohibitively large storage requirements and computational costs of the numerical solution. In view of memory efficient approximation, we review and extend known solution formulas and identify invariant subspaces for a possibly low-dimensional solution representation. Based on these theoretical findings, we propose a Galerkin projection onto a space related to a low-rank approximation of the algebraic Riccati equation. For the numerical implementation, we provide an alternative interpretation of the modified \emph{Davison-Maki method} via the transformed flow of the differential Riccati equation, which enables us to rule out known stability issues of the method in combination with the proposed projection scheme. We present numerical experiments for large-scale autonomous differential Riccati equations and compare our approach with high-order splitting schemes.


title: Floquet topological insulator laser

abstract: We introduce a class of topological lasers based on the photonic Floquet topological insulator concept. The proposed system is realized as a truncated array of the lasing helical waveguides, where the pseudo-magnetic field arises due to twisting of the waveguides along the propagation direction that breaks the time-reversal symmetry and opens up a topological gap. When sufficient gain is provided in the edge channels of the array then the system lases into the topological edge states. Topological lasing is stable only in certain intervals of the Bloch momenta, that ensure a dynamic, but stable balance between the linear amplification and nonlinear absorption leading to the formation of the breathing edge states. We also illustrate topological robustness of the edge currents by simulating lattice defects and triangular arrangements of the waveguides.


title: Regular polygonal equilibrium configurations on S^1 and stability of the associated relative equilibria

abstract: For the curved n-body problem in S^3, we show that a regular polygonal configuration for n masses on a geodesic is an equilibrium configuration if and only if n is odd and the masses are equal. The equilibrium configuration is associated with a one-parameter family (depending on the angular velocity) of relative equilibria, which take place on S^1 embedded in S^2. We then study the stability of the associated relative equilibria on two invariant manifolds, T^*((\S^1)^n\D) and T^*((\S^2)^n\D). We show that they are Lyapunov stable on S^1, they are Lyapunov stable on S^2 if the absolute value of angular velocity is larger than a certain value, and that they are linearly unstable on S^2 if the absolute value of angular velocity is smaller than that certain value.


title: ANAIS-112 status: two years results on annual modulation

abstract: ANAIS (Annual modulation with NaI Scintillators) is a dark matter direct detection experiment located at the Canfranc Underground Laboratory (LSC), in Spain. The goal is to confirm or refute in a model independent way the DAMA/LIBRA positive result: an annual modulation in the low-energy detection rate compatible with the expected signal induced by dark matter particles in the galactic halo. This signal, observed for about 20 years, is in strong tension with the negative results of other very sensitive experiments, but a direct comparison using the same target material, NaI(Tl), was still lacking. ANAIS-112, consisting of 112.5 kg of NaI(Tl) scintillators, is taking data at the LSC since August 2017. Here we present the preliminary annual modulation analysis corresponding to two years of data (exposure of 220.69 kg x y) and the ANAIS-112 projected sensitivity for the scheduled 5 y of operation.


title: Voronoi tesselation analysis of sets of randomly placed finite-size spheres

abstract: The purpose of this note is to clarify the effect of the finite size of spherical particles upon the characteristics of their spatial distribution through a random Poisson process (RPP). This information is of special interest when using RPP data as a reference for the analysis of the spatial structure of a given (non-RPP) particulate system, in which case ignoring finite-size effects upon the former may yield misleading conclusions. We perform Monte Carlo simulations in triply-periodic spatial domains, and then analyze the particle-centered Voronoi tesselations for solid volume fractions ranging from 10^(-5) to 0.3. We show that the standard-deviation of these volumes decreases with the solid volume fraction, the deviation from the value of point sets being reasonably approximated by an exponential function. As can be expected, the domain size for which the random assemblies of finite-size particles are generated has a constraining effect if the number of particles per realization is chosen too small. This effect is quantified, and recommendations are given. We have also revisited the case of random point sets (i.e. the limit of vanishing particle diameter), for which we have confirmed the accuracy of the earlier data by Tanemura [Forma, 18(4):221-247, 2003].


title: Derivation and Analysis of Fast Bilinear Algorithms for Convolution

abstract: The prevalence of convolution in applications within signal processing, deep neural networks, and numerical solvers has motivated the development of numerous fast convolution algorithms. In many of these problems, convolution is performed on terabytes or petabytes of data, so even constant factors of improvement can significantly reduce the computation time. We leverage the formalism of bilinear algorithms to describe and analyze all of the most popular approaches. This unified lens permits us to study the relationship between different variants of convolution as well as to derive error bounds and analyze the cost of the various algorithms. We provide new derivations, which predominantly leverage matrix and tensor algebra, to describe the Winograd family of convolution algorithms as well as reductions between 1D and multidimensional convolution. We provide cost and error bounds as well as experimental numerical studies. Our experiments for two of these algorithms, the overlap-add approach and Winograd convolution algorithm with polynomials of degree greater than one, show that fast convolution algorithms can rival the accuracy of the fast Fourier transform (FFT) without using complex arithmetic. These algorithms can be used for convolution problems with multidimensional inputs or for filters larger than size of four, extending the state-of-the-art in Winograd-based convolution algorithms.


title: Viscosity, non-conformal equation of state and sound velocity in Landau hydrodynamics

abstract: We find an analytical solution to relativistic viscous hydrodynamics for a 1+1 dimensional Landau flow profile. We consider relativistic Navier-Stokes form of the dissipative hydrodynamic equation, for a non-conformal system with a constant speed of sound, and employ the obtained solution to fit rapidity spectrum of observed pions in $\sqrt{s_{NN}}=$ 200, 17.3, 12.3, 8.76, 7.62, 6.27, 4.29, 3.83, 3.28 and 2.63 GeV collision energies. We find that at the freeze-out hypersurface with improved Landau's freeze-out prescription, the viscous corrections do not affect the rapidity spectra. We demonstrate that the solution of the non-conformal Landau flow lead to a better agreement with the experimental data compared to the conformal ideal solution. We also extract speed of sound from fit to the rapidity spectra for various collision energies and find a monotonous decrease with decreasing collision energies. Appealing to the fact that viscosity has negligible effect on rapidity spectra for Landau's freeze-out scenario, we argue that our calculations provides a framework for extracting the average value of speed of sound in relativistic heavy-ion collisions.


title: A Hamilton-Jacobi Reachability-Based Framework for Predicting and Analyzing Human Motion for Safe Planning

abstract: Real-world autonomous systems often employ probabilistic predictive models of human behavior during planning to reason about their future motion. Since accurately modeling human behavior a priori is challenging, such models are often parameterized, enabling the robot to adapt predictions based on observations by maintaining a distribution over the model parameters. Although this enables data and priors to improve the human model, observation models are difficult to specify and priors may be incorrect, leading to erroneous state predictions that can degrade the safety of the robot motion plan. In this work, we seek to design a predictor which is more robust to misspecified models and priors, but can still leverage human behavioral data online to reduce conservatism in a safe way. To do this, we cast human motion prediction as a Hamilton-Jacobi reachability problem in the joint state space of the human and the belief over the model parameters. We construct a new continuous-time dynamical system, where the inputs are the observations of human behavior, and the dynamics include how the belief over the model parameters change. The results of this reachability computation enable us to both analyze the effect of incorrect priors on future predictions in continuous state and time, as well as to make predictions of the human state in the future. We compare our approach to the worst-case forward reachable set and a stochastic predictor which uses Bayesian inference and produces full future state distributions. Our comparisons in simulation and in hardware demonstrate how our framework can enable robust planning while not being overly conservative, even when the human model is inaccurate.


title: Particle acceleration in kink-unstable jets

abstract: Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instability first generates large scale current sheets, which at later times break up into small-scale turbulence. Reconnection in these sheets proceeds in the strong guide field regime, which results in a formation of steep power laws in the particle spectra. Later evolution shows heating of the plasma, which is driven by small-amplitude turbulence induced by the kink instability. These two processes energize particles due to a combination of ideal and non-ideal electric fields.


title: Magneto-optical Kerr effect in spin split two-dimensional massive Dirac materials

abstract: Two-dimensional (2D) massive Dirac electrons possess a finite Berry curvature, with Chern number $\pm 1/2$, that entails both a quantized dc Hall response and a subgap full-quarter Kerr rotation. The observation of these effects in 2D massive Dirac materials such as gapped graphene, hexagonal boron nitride or transition metal dichalcogenides (TMDs) is obscured by the fact that Dirac cones come in pairs with opposite sign Berry curvatures, leading to a vanishing Chern number. Here, we show that the presence of spin-orbit interactions, combined with an exchange spin splitting induced either by diluted magnetic impurities or by proximity to a ferromagnetic insulator, gives origin to a net magneto-optical Kerr effect in such systems. We focus on the case of TMD monolayers and study the dependence of Kerr rotation on frequency and exchange spin splitting. The role of the substrate is included in the theory and found to critically affect the results. Our calculations indicate that state-of-the-art magneto-optical Kerr spectroscopy can detect a single magnetic impurity in diluted magnetic TMDs.


title: Predicting gait events from tibial acceleration in rearfoot running: a structured machine learning approach

abstract: Gait event detection of the initial contact and toe off is essential for running gait analysis, allowing the derivation of parameters such as stance time. Heuristic-based methods exist to estimate these key gait events from tibial accelerometry. However, these methods are tailored to very specific acceleration profiles, which may offer complications when dealing with larger data sets and inherent biological variability. Therefore, this paper investigates whether a structured machine learning approach can achieve a more accurate prediction of running gait event timings from tibial accelerometry. Force-based event detection acted as the criterion measure in order to assess the accuracy, repeatability and sensitivity of the predicted gait events. A heuristic method and two structured machine learning methods were employed to derive initial contact, toe off and stance time from tibial acceleration signals. Both a structured perceptron model (median absolute error of stance time estimation: 10.00 $\pm$ 8.73 ms) and a structured recurrent neural network model (median absolute error of stance time estimation: 6.50 $\pm$ 5.74 ms) significantly outperformed the existing heuristic approach (median absolute error of stance time estimation: 11.25 $\pm$ 9.52 ms) on data from 93 rearfoot runners. Thus, results indicate that a structured recurrent neural network machine learning model offers the most accurate and consistent estimation of the gait events and its derived stance time during level overground running. The machine learning methods seem less affected by intra- and inter-subject variation within the data, allowing for accurate and efficient automated data output during rearfoot overground running. Furthermore offering possibilities for real-time monitoring and biofeedback during prolonged measurements, even outside the laboratory.


title: Decomposing Collectives for Exploiting Multi-lane Communication

abstract: Many modern, high-performance systems increase the cumulated node-bandwidth by offering more than a single communication network and/or by having multiple connections to the network. Efficient algorithms and implementations for collective operations as found in, e.g., MPI must be explicitly designed for such multi-lane capabilities. We discuss a model for the design of multi-lane algorithms, and in particular give a recipe for converting any standard, one-ported, (pipelined) communication tree algorithm into a multi-lane algorithm that can effectively use $k$ lanes simultaneously. We first examine the problem from the perspective of \emph{self-consistent performance guidelines}, and give simple, \emph{full-lane, mock-up implementations} of the MPI broadcast, reduction, scan, gather, scatter, allgather, and alltoall operations using only similar operations of the given MPI library itself in such a way that multi-lane capabilities can be exploited. These implementations which rely on a decomposition of the communication domain into communicators for nodes and lanes are full-fledged and readily usable implementations of the MPI collectives. The mock-up implementations, contrary to expectation, in many cases show surprising performance improvements with different MPI libraries on a small 36-node dual-socket, dual-lane Intel OmniPath cluster, indicating severe problems with the native MPI library implementations. Our full-lane implementations are in many cases considerably more than a factor of two faster than the corresponding MPI collectives. We see similar results on the larger Vienna Scientific Cluster, VSC-3. These experiments indicate considerable room for improvement of the MPI collectives in current libraries including more efficient use of multi-lane communication.


title: Instantons and L-space surgeries

abstract: We prove that instanton L-space knots are fibered and strongly quasipositive. Our proof differs conceptually from proofs of the analogous result in Heegaard Floer homology, and includes a new decomposition theorem for cobordism maps in framed instanton Floer homology akin to the $\textrm{Spin}^c$ decompositions of cobordism maps in other Floer homology theories. As our main application, we prove (modulo a mild nondegeneracy condition) that for $r$ a positive rational number and $K$ a nontrivial knot in the $3$-sphere, there exists an irreducible homomorphism \[\pi_1(S^3_r(K)) \to SU(2)\] unless $r \geq 2g(K)-1$ and $K$ is both fibered and strongly quasipositive, broadly generalizing results of Kronheimer and Mrowka. We also answer a question of theirs from 2004, proving that there is always an irreducible homomorphism from the fundamental group of 4-surgery on a nontrivial knot to $SU(2)$. In another application, we show that a slight enhancement of the A-polynomial detects infinitely many torus knots, including the trefoil.


title: Origin of the Magnetic Spin Hall Effect: Spin Current Vorticity in the Fermi Sea

abstract: The interplay of spin-orbit coupling (SOC) and magnetism gives rise to a plethora of charge-to-spin conversion phenomena that harbor great potential for spintronics applications. In addition to the spin Hall effect, magnets may exhibit a magnetic spin Hall effect (MSHE), as was recently discovered [Kimata \textit{et al.}, Nature \textbf{565}, 627-630 (2019)]. To date, the MSHE is still awaiting its intuitive explanation. Here we relate the MSHE to the vorticity of spin currents in the Fermi sea, which explains pictorially the origin of the MSHE. For all magnetic Laue groups that allow for nonzero spin current vorticities the related tensor elements of the MSH conductivity are given. Minimal requirements for the occurrence of a MSHE are compatibility with either a magnetization or a magnetic toroidal quadrupole. This finding implies in particular that the MSHE is expected in all ferromagnets with sufficiently large SOC. To substantiate our symmetry analysis, we present various models, in particular a two-dimensional magnetized Rashba electron gas, that corroborate an interpretation by means of spin current vortices. Considering thermally induced spin transport and the magnetic spin Nernst effect in magnetic insulators, which are brought about by magnons, our findings for electron transport can be carried over to the realm of spincaloritronics, heat-to-spin conversion, and energy harvesting.


title: How Much Can We See? A Note on Quantifying Explainability of Machine Learning Models

abstract: One of the most popular approaches to understanding feature effects of modern black box machine learning models are partial dependence plots (PDP). These plots are easy to understand but only able to visualize low order dependencies. The paper is about the question 'How much can we see?': A framework is developed to quantify the explainability of arbitrary machine learning models, i.e. up to what degree the visualization as given by a PDP is able to explain the predictions of the model. The result allows for a judgement whether an attempt to explain a black box model is sufficient or not.


title: Radiation reaction in electron-beam interactions with high-intensity lasers

abstract: Charged particles accelerated by electromagnetic fields emit radiation, which must, by the conservation of momentum, exert a recoil on the emitting particle. The force of this recoil, known as radiation reaction, strongly affects the dynamics of ultrarelativistic electrons in intense electromagnetic fields. Such environments are found astrophysically, e.g. in neutron star magnetospheres, and will be created in laser-matter experiments in the next generation of high-intensity laser facilities. In many of these scenarios, the energy of an individual photon of the radiation can be comparable to the energy of the emitting particle, which necessitates modelling not only of radiation reaction, but quantum radiation reaction. The worldwide development of multi-petawatt laser systems in large-scale facilities, and the expectation that they will create focussed electromagnetic fields with unprecedented intensities $> 10^{23}~\mathrm{W}\text{cm}^{-2}$, has motivated renewed interest in these effects. In this paper I review theoretical and experimental progress towards understanding radiation reaction, and quantum effects on the same, in high-intensity laser fields that are probed with ultrarelativistic electron beams. In particular, we will discuss how analytical and numerical methods give insight into new kinds of radiation-reaction-induced dynamics, as well as how the same physics can be explored in experiments at currently existing laser facilities.


title: Random plane partitions and corner distributions

abstract: We explore some probabilistic applications arising in connections with $K$-theoretic symmetric functions. For instance, we determine certain corner distributions of random lozenge tilings and plane partitions. We also introduce some distributions that are naturally related to the corner growth model. Our main tools are dual symmetric Grothendieck polynomials and normalized Schur functions.


title: A 3D short-characteristics method for continuum and line scattering problems in the winds of hot stars

abstract: Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the UV resonance line formation in the winds of rapidly rotating O stars. Methods: We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic B\'ezier interpolations, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated $\Lambda$-iteration scheme. Results: For the spherical test models, the mean relative error of the source function is on the $5-20\,\%$ level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion. Conclusions: Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.


title: Inhomogeneous cosmology in an anisotropic Universe

abstract: With the era of precision cosmology upon us, and upcoming surveys expected to further improve the precision of our observations below the percent level, ensuring the accuracy of our theoretical cosmological model is of the utmost importance. Current tensions between our observations and predictions from the standard cosmological model have sparked curiosity in extending the model to include new physics. Although, some suggestions include simply accounting for aspects of our Universe that are ignored in the standard model. One example acknowledges the fact that our Universe contains significant density contrasts on small scales; in the form of galaxies, galaxy clusters, filaments, and voids. This small-scale structure is smoothed out in the standard model, by assuming large-scale homogeneity of the matter distribution, which could have a measurable effect due to the nonlinearity of Einstein's equations. This backreaction of small-scale structures on the large-scale dynamics has been suggested to explain the measured accelerating expansion rate of the Universe. Current standard cosmological simulations ignore the effects of General Relativity by assuming purely Newtonian dynamics. In this thesis, we take the first steps towards quantifying the backreaction of small-scale structures by performing cosmological simulations that solve Einstein's equations directly. Simulations like these will allow us to quantify potentially important effects on our observations that could become measurable as the precision of these observations increases into the future.


title: Local Wiener's Theorem and Coherent Sets of Frequencies

abstract: Using a local analog of the Wiener-Levi theorem, we investigate the class of measures on Euclidean space with discrete support and spectrum. Also, we find a new sufficient conditions for a discrete set in Euclidean space to be a coherent set of frequencies.


title: New Quantum Structure of the Space-Time

abstract: We start from quantum theory (instead of general relativity) to approach quantum gravity within a minimal setting and promote the space-time coordinates to quantum non-commuting operators. Comparison to the harmonic oscillator global phase space is enlighting. The phase space instanton (X, P = iT) describes the hyperbolic quantum space-time structure and generates the quantum light cone. The classical Minkowski space-time null generators X = T dissapear at the quantum level replaced by four Planck scale hyperbolae and a new quantum Planck scale vacuum region emerges. We describe the quantum Rindler and quantum Schwarzshild-Kruskal space-time structures. The horizons and the r=0 singularity are quantum mechanically erased. The four Kruskal regions merge inside a single quantum Planck scale "world". The quantum space-time structure consists of hyperbolic discrete levels of odd numbers (2n + 1) (in Planck units ), n = 0, 1, 2.... . The mass levels are the square root of (2n + 1). Large n are semiclassical tending towards a classical continuum space-time. Low n are quantum, the lowest mode (n = 0) being the Planck scale. Two dual branches are present covering the whole mass spectrum: from the largest astrophysical objects to the quantum elementary particles passing by the Planck mass. Black holes belong to both branches. (Abridged)


title: XFe4Ge2 (X = Y, Lu) and Mn3Pt: Filling-enforced magnetic topological metals

abstract: Magnetism, coupled with nontrivial band topology, can bring about many interesting and exotic phenomena, so that magnetic topological materials have attracted persistent research interest. However, compared with non-magnetic topological materials (TMs), the magnetic TMs are less studied, since their magnetic structures and topological phase transitions are usually complex and the first-principles predictions are usually sensitive on the effect of Coulomb interaction. In this work, we present a comprehensive investigation of XFe4Ge2 (X = Y, Lu) and Mn3Pt, and find these materials to be filling-enforced magnetic topological metals. Our first-principles calculations show that XFe4Ge2 (X = Y, Lu) host Dirac points near the Fermi level at high symmetry point S. These Dirac points are protected by PT symmetry (P and T are inversion and time-reversal transformations, respectively) and a 2-fold screw rotation symmetry. Moreover, through breaking PT symmetry, the Dirac points would split into Weyl nodes. Mn3Pt is found to host 4-fold degenerate band crossings in the whole high symmetry path of A-Z. We also utilize the GGA+U scheme to take into account the effect of Coulomb repulsion and find that the filling-enforced topological properties are naturally insensitive on U.


title: Sums of two squares in short intervals

abstract: We show that there are short intervals $[x,x+y]$ containing $\gg y^{1/10}$ numbers expressible as the sum of two squares, which is many more than the average when $y=o( (\log{x})^{5/9})$. We obtain similar results for sums of two squares in short arithmetic progressions.


title: Hipsters and the Cool: A Game Theoretic Analysis of Social Identity, Trends and Fads

abstract: Cultural trends and popularity cycles can be observed all around us, yet our theories of social influence and identity expression do not explain what perpetuates these complex, often unpredictable social dynamics. We propose a theory of social identity expression based on the opposing, but not mutually exclusive, motives to conform and to be unique among one's neighbors in a social network. We then model the social dynamics that arise from these motives. We find that the dynamics typically enter random walks or stochastic limit cycles rather than converging to a static equilibrium. We also prove that without social network structure or, alternatively, without the uniqueness motive, reasonable adaptive dynamics would necessarily converge to equilibrium. Thus, we show that nuanced psychological assumptions (recognizing preferences for uniqueness along with conformity) and realistic social network structure are both necessary for explaining how complex, unpredictable cultural trends emerge.


title: NC Algorithms for Popular Matchings in One-Sided Preference Systems and Related Problems

abstract: The popular matching problem is of matching a set of applicants to a set of posts, where each applicant has a preference list, ranking a non-empty subset of posts in the order of preference, possibly with ties. A matching M is popular if there is no other matching M' such that more applicants prefer M' to M. We give the first NC algorithm to solve the popular matching problem without ties. We also give an NC algorithm that solves the maximum-cardinality popular matching problem. No NC or RNC algorithms were known for the matching problem in preference systems prior to this work. Moreover, we give an NC algorithm for a weaker version of the stable matching problem, that is, the problem of finding the "next" stable matching given a stable matching.


title: Direct Optical-Structure Correlation in Atomically Thin Dichalcogenides and Heterostructures

abstract: Atomically thin transition metal dichalcogenides (TMDs) have distinct opto-electronic properties including enhanced luminescence and high on-off current ratios, which can be further modulated by making more complex TMD heterostructures. However, resolution limits of conventional optical methods do not allow for direct optical-structural correlation measurements in these materials, particularly of buried interfaces in TMD heterostructures. Here we use, for the first time, electron beam induced cathodoluminescence in a scanning transmission electron microscope (CL-STEM) to measure optical properties of monolayer TMDs (WS2, MoS2 and WSSe alloy) encapsulated between layers of hBN. We observe dark areas resulting from localized (~ 100 nm) imperfect interfaces and monolayer folding, which shows that the intimate contact between layers in this application-relevant heterostructure is required for proper inter layer coupling. We also realize a suitable imaging method that minimizes electron-beam induced changes and provides measurement of intrinsic properties. To overcome the limitation of small electron interaction volume in TMD monolayer (and hence low photon yield), we find that encapsulation of TMD monolayers with hBN and subsequent annealing is important. CL-STEM offers to be a powerful method to directly measure structure-optical correspondence in lateral or vertical heterostructures and alloys.


title: Neutrino clustering in the Milky Way and beyond

abstract: The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent $N$-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.


title: Learning Sparse Distributions using Iterative Hard Thresholding

abstract: Iterative hard thresholding (IHT) is a projected gradient descent algorithm, known to achieve state of the art performance for a wide range of structured estimation problems, such as sparse inference. In this work, we consider IHT as a solution to the problem of learning sparse discrete distributions. We study the hardness of using IHT on the space of measures. As a practical alternative, we propose a greedy approximate projection which simultaneously captures appropriate notions of sparsity in distributions, while satisfying the simplex constraint, and investigate the convergence behavior of the resulting procedure in various settings. Our results show, both in theory and practice, that IHT can achieve state of the art results for learning sparse distributions.


title: Ferromagnetic Kitaev interaction and the origin of large magnetic anisotropy in $\alpha$-RuCl$_3$

abstract: $\alpha$-RuCl$_3$ is drawing much attention as a promising candidate Kitaev quantum spin liquid. However, despite intensive research efforts, controversy remains about the form of the basic interactions governing the physics of this material. Even the sign of the Kitaev interaction (the bond-dependent anisotropic interaction responsible for Kitaev physics) is still under debate, with conflicting results from theoretical and experimental studies. The significance of the symmetric off-diagonal exchange interaction (referred to as the $\Gamma$ term) is another contentious question. Here, we present resonant elastic x-ray scattering data that provides unambiguous experimental constraints to the two leading terms in the magnetic interaction Hamiltonian. We show that the Kitaev interaction ($K$) is ferromagnetic, and that the $\Gamma$ term is antiferromagnetic and comparable in size to the Kitaev interaction. Our findings also provide a natural explanation for the large anisotropy of the magnetic susceptibility in $\alpha$-RuCl$_3$ as arising from the large $\Gamma$ term. We therefore provide a crucial foundation for understanding the interactions underpinning the exotic magnetic behaviours observed in $\alpha$-RuCl$_3$.


title: Bessel $F$-isocrystals for reductive groups

abstract: We construct the Frobenius structure on a rigid connection $\mathrm{Be}_{\check{G}}$ on $\mathbb{G}_m$ for a split reductive group $\check{G}$ introduced by Frenkel-Gross. These data form a $\check{G}$-valued overconvergent $F$-isocrystal $\mathrm{Be}_{\check{G}}^{\dagger}$ on $\mathbb{G}_{m,\mathbb{F}_p}$, which is the $p$-adic companion of the Kloosterman $\check{G}$-local system $\mathrm{Kl}_{\check{G}}$ constructed by Heinloth-Ng\^o-Yun. By exploring the structure of the underlying differential equation, we calculate the monodromy group of $\mathrm{Be}_{\check{G}}^{\dagger}$ when $\check{G}$ is almost simple (which recovers the calculation of monodromy group of $\mathrm{Kl}_{\check{G}}$ due to Katz and Heinloth-Ng\^o-Yun), and establish functoriality between different Kloosterman $\check{G}$-local systems as conjectured by Heinloth-Ng\^o-Yun. We show that the Frobenius Newton polygons of $\mathrm{Kl}_{\check{G}}$ are generically ordinary for every $\check{G}$ and are everywhere ordinary on $|\mathbb{G}_{m,\mathbb{F}_p}|$ when $\check{G}$ is classical or $G_2$.


title: Two body interactions induces the axion-phason field in Weyl semimetals

abstract: Following our results that the two-body interaction can induce a space and time dependent topological axion term $ \frac{e^2}{2\pi \hbar}\theta(z,t)(\vec{E}\cdot\vec{B})$, we show that by applying this theory to a Weyl semimetal with two nodes in a magnetic field a one-dimensional sliding charge density wave (CDW) in agreement with the recent experimental finding of J. Gooth et al. [Nature, https://doi.org/10.1038/s41586-019-1630-4 (2019)]. We show that the theory is equivalent to a time and space dependent of the topological angle which represents the phason.


title: Constrained Reinforcement Learning Has Zero Duality Gap

abstract: Autonomous agents must often deal with conflicting requirements, such as completing tasks using the least amount of time/energy, learning multiple tasks, or dealing with multiple opponents. In the context of reinforcement learning~(RL), these problems are addressed by (i)~designing a reward function that simultaneously describes all requirements or (ii)~combining modular value functions that encode them individually. Though effective, these methods have critical downsides. Designing good reward functions that balance different objectives is challenging, especially as the number of objectives grows. Moreover, implicit interference between goals may lead to performance plateaus as they compete for resources, particularly when training on-policy. Similarly, selecting parameters to combine value functions is at least as hard as designing an all-encompassing reward, given that the effect of their values on the overall policy is not straightforward. The later is generally addressed by formulating the conflicting requirements as a constrained RL problem and solved using Primal-Dual methods. These algorithms are in general not guaranteed to converge to the optimal solution since the problem is not convex. This work provides theoretical support to these approaches by establishing that despite its non-convexity, this problem has zero duality gap, i.e., it can be solved exactly in the dual domain, where it becomes convex. Finally, we show this result basically holds if the policy is described by a good parametrization~(e.g., neural networks) and we connect this result with primal-dual algorithms present in the literature and we establish the convergence to the optimal solution.


title: Spectroscopy, lifetime and decay modes of the $T^-_{bb}$ tetraquark

abstract: We present the first full-fledged study of the flavor-exotic isoscalar $T_{bb}^-\equiv b b \bar u \bar d$ tetraquark with spin and parity $J^P=1^+$. We report accurate solutions of the four-body problem in a quark model, characterizing the structure of the state as a function of the ratio $M_Q/m_q$ of the heavy to light quark masses. For such a standard constituent model, $T_{bb}^-$ lies approximately 150 MeV below the strong decay threshold $B^-\bar {B^*}^{0}$ and 105 MeV below the electromagnetic decay threshold $B^- \bar B^0 \gamma$. We evaluate the lifetime of $T_{bb}^-$, identifying the promising decay modes where the tetraquark might be looked for in future experiments. Its total decay width is $\Gamma \approx 87 \times 10^{-15}$ GeV and therefore its lifetime $\tau \approx$ 7.6 ps. The promising final states are ${B^*}^{-}\, {D^*}^{+} \, \ell^- \, \bar \nu_\ell$ and $\bar {B^*}^{0} \, {D^*}^{0} \, \ell^- \, \bar \nu_\ell $ among the semileptonic decays, and ${B^*}^{-} \, {D^*}^{+} \, {D_s^*}^-$, $\bar {B^*}^{0} \, {D^*}^{0} \, {D_s^*}^- $, and ${B^*}^{-} \, {D^*}^{+} \, \rho^-$ among the nonleptonic ones. The semileptonic decay to the isoscalar $J^P=0^+$ tetraquark $T_{bc}^0$ is also relevant but it is not found to be dominant. There is a broad consensus about the existence of this tetraquark, and its detection will validate our understanding of the low-energy realizations of Quantum Chromodynamics (QCD) in the multiquark sector.


title: Dynamics Learning with Cascaded Variational Inference for Multi-Step Manipulation

abstract: The fundamental challenge of planning for multi-step manipulation is to find effective and plausible action sequences that lead to the task goal. We present Cascaded Variational Inference (CAVIN) Planner, a model-based method that hierarchically generates plans by sampling from latent spaces. To facilitate planning over long time horizons, our method learns latent representations that decouple the prediction of high-level effects from the generation of low-level motions through cascaded variational inference. This enables us to model dynamics at two different levels of temporal resolutions for hierarchical planning. We evaluate our approach in three multi-step robotic manipulation tasks in cluttered tabletop environments given high-dimensional observations. Empirical results demonstrate that the proposed method outperforms state-of-the-art model-based methods by strategically interacting with multiple objects.


title: $Z$ boson production in Pb+Pb collisions at $\sqrt{s_{\textrm{NN}}}$= 5.02 TeV measured by the ATLAS experiment

abstract: The production yield of $Z$ bosons is measured in the electron and muon decay channels in Pb+Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV with the ATLAS detector. Data from the 2015 LHC run corresponding to an integrated luminosity of 0.49 $\textrm{nb}^{-1}$ are used for the analysis. The $Z$ boson yield, normalised by the total total number of minimum-bias events and the mean nuclear thickness function, is measured as a function of dilepton rapidity and event centrality. The measurements in Pb+Pb collisions are compared with similar measurements made in proton$-$proton collisions at the same centre-of-mass energy. The nuclear modification factor is found to be consistent with unity for all centrality intervals. The results are compared with theoretical predictions obtained at next-to-leading order using nucleon and nuclear parton distribution functions. The normalised $Z$ boson yields in Pb+Pb collisions lie 1$-3$$\sigma$ above the predictions. The nuclear modification factor measured as a function of rapidity agrees with unity and is consistent with a next-to-leading-order QCD calculation including the isospin effect.


title: Reproducing Scientific Experiment with Cloud DevOps

abstract: The reproducibility of scientific experiment is vital for the advancement of disciplines based on previous work. To achieve this goal, many researchers focus on complex methodology and self-invented tools which have difficulty in practical usage. In this article, we introduce the Cloud DevOps infrastructure from software engineering community and shows how it can be used effectively for heterogeneous agents to reproduce experiments for computer science related disciplines. DevOps can be enabled using freely available cloud computing machines for medium-sized experiment and self-hosted computing engines for large-scale computing, thus powering researchers to share their experiment result with others in a more reliable way.


title: Stein's Lemma for the Reparameterization Trick with Exponential Family Mixtures

abstract: Stein's method (Stein, 1973; 1981) is a powerful tool for statistical applications, and has had a significant impact in machine learning. Stein's lemma plays an essential role in Stein's method. Previous applications of Stein's lemma either required strong technical assumptions or were limited to Gaussian distributions with restricted covariance structures. In this work, we extend Stein's lemma to exponential-family mixture distributions including Gaussian distributions with full covariance structures. Our generalization enables us to establish a connection between Stein's lemma and the reparamterization trick to derive gradients of expectations of a large class of functions under weak assumptions. Using this connection, we can derive many new reparameterizable gradient-identities that goes beyond the reach of existing works. For example, we give gradient identities when expectation is taken with respect to Student's t-distribution, skew Gaussian, exponentially modified Gaussian, and normal inverse Gaussian.


title: Robust Model-free Reinforcement Learning with Multi-objective Bayesian Optimization

abstract: In reinforcement learning (RL), an autonomous agent learns to perform complex tasks by maximizing an exogenous reward signal while interacting with its environment. In real-world applications, test conditions may differ substantially from the training scenario and, therefore, focusing on pure reward maximization during training may lead to poor results at test time. In these cases, it is important to trade-off between performance and robustness while learning a policy. While several results exist for robust, model-based RL, the model-free case has not been widely investigated. In this paper, we cast the robust, model-free RL problem as a multi-objective optimization problem. To quantify the robustness of a policy, we use delay margin and gain margin, two robustness indicators that are common in control theory. We show how these metrics can be estimated from data in the model-free setting. We use multi-objective Bayesian optimization (MOBO) to solve efficiently this expensive-to-evaluate, multi-objective optimization problem. We show the benefits of our robust formulation both in sim-to-real and pure hardware experiments to balance a Furuta pendulum.


title: Derivation of the HOMFLYPT knot polynomial via helicity and geometric quantization

abstract: In this Letter we propose a semiclassical interpretation of the HOMFLYPT polynomial building on the Liu-Ricca hydrodynamical approach to the latter and on the Besana-S. symplectic approach to framing via Brylinski's manifold of mildly singular links.


title: Model enhancement and personalization using weakly supervised learning for multi-modal mobile sensing

abstract: Always-on sensing of mobile device user's contextual information is critical to many intelligent use cases nowadays such as healthcare, drive assistance, voice UI. State-of-the-art approaches for predicting user context have proved the value to leverage multiple sensing modalities for better accuracy. However, those context inference algorithms that run on application processor nowadays tend to drain heavy amount of power, making them not suitable for an always-on implementation. We claim that not every sensing modality is suitable to be activated all the time and it remains challenging to build an inference engine using power friendly sensing modalities. Meanwhile, due to the diverse population, we find it challenging to learn a context inference model that generalizes well, with limited training data, especially when only using always-on low power sensors. In this work, we propose an approach to leverage the opportunistically-on counterparts in device to improve the always-on prediction model, leading to a personalized solution. We model this problem using a weakly supervised learning framework and provide both theoretical and experimental results to validate our design. The proposed framework achieves satisfying result in the IMU based activity recognition application we considered.


title: Digits of primes

abstract: We survey some different results on the digits of prime numbers, giving a simplified proof of weak forms of a result of Maynard and Mauduit-Rivat.


title: Generalized parabolic structures over smooth curves with many components and principal bundles over reducible nodal curves

abstract: Let $Y_{1},\dots,Y_{l}$ be smooth irreducible projective curves and let $Y$ be its disjoint union. Given a semisimple reductive algebraic group $G$ and a faithful representation $\rho:G\hookrightarrow \textrm{SL}(V)$ we construct a projective moduli space of $(\underline{\kappa},\delta)$-(semi)stable singular principal $G$-bundles with generalized parabolic structure of type $\underline{e}$. In case $Y$ is the normalization of a connected and reducible projective nodal curve $X$, there is a closed subscheme coarsely representing the subfunctor corresponding to descending bundles. We prove that the descent operation induces a birational, surjective and proper morphism onto the schematic closure of the space of $\delta$-stable singular principal $G$-bundles whose associated torsion free sheaf is of local type $\underline{e}$.


title: Measurement of the $B_c^-$ meson production fraction and asymmetry in 7 and 13 TeV $pp$ collisions

abstract: The production fraction of the $B_c^-$ meson with respect to the sum of $B^-$ and $\bar{B}^0$ mesons is measured in both 7 and 13 TeV center-of-mass energy $pp$ collisions produced by the Large Hadron Collider (LHC), using the LHCb detector. The rate, approximately 3.7 per mille, does not change with energy, but shows a transverse momentum dependence. The $B_c^- - B_c^+$ production asymmetry is also measured, and is consistent with zero within the determined statistical and systematic uncertainties of a few percent.


title: Experimental Comparison of Bohm-like Theories with Different Ontologies

abstract: The de Broglie-Bohm theory is a hidden variable interpretation of quantum mechanics which involves particles moving through space with definite trajectories. This theory singles out position as the primary ontological variable. Mathematically, it is possible to construct a similar theory where particles are moving through momentum space, and momentum is singled out as the primary ontological variable. In this paper we experimentally show how the two theories lead to different ontological descriptions. We construct the putative particle trajectories for a two-slit experiment in both the position and momentum space theories by simulating particle dynamics with coherent light. Using a method for constructing trajectories through the primary and derived (i.e. non-primary) spaces, we compare the ontological pictures offered by the two theories and show that they do not agree. This contradictory behaviour brings into question which ontology for Bohmian mechanics is to be preferred.


title: Generalization of Reinforcement Learners with Working and Episodic Memory

abstract: Memory is an important aspect of intelligence and plays a role in many deep reinforcement learning models. However, little progress has been made in understanding when specific memory systems help more than others and how well they generalize. The field also has yet to see a prevalent consistent and rigorous approach for evaluating agent performance on holdout data. In this paper, we aim to develop a comprehensive methodology to test different kinds of memory in an agent and assess how well the agent can apply what it learns in training to a holdout set that differs from the training set along dimensions that we suggest are relevant for evaluating memory-specific generalization. To that end, we first construct a diverse set of memory tasks that allow us to evaluate test-time generalization across multiple dimensions. Second, we develop and perform multiple ablations on an agent architecture that combines multiple memory systems, observe its baseline models, and investigate its performance against the task suite.


title: On-shell constructibility of Born amplitudes in spontaneously broken gauge theories

abstract: We perform a comprehensive study of on-shell recursion relations for Born amplitudes in spontaneously broken gauge theories and identify the minimal shifts required to construct amplitudes with a given particle content and spin quantum numbers. We show that two-line or three-line shifts are sufficient to construct all amplitudes with five or more particles, apart from amplitudes involving longitudinal vector bosons or scalars, which may require at most five-line shifts. As an application, we revisit selection rules for multi-boson amplitudes using on-shell recursion and little-group transformations.


title: A framework for deep learning emulation of numerical models with a case study in satellite remote sensing

abstract: Numerical models based on physics represent the state-of-the-art in earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest-generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning, across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in earth systems. A difficult test for deep learning-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A deep learning emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here we examine, with a case study in satellite-based remote sensing, the hypothesis that deep learning approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the deep learning emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of deep learning as well as the growing desire for higher-resolution simulations in the earth sciences.


title: GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

abstract: On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.


title: Mocking the u-plane integral

abstract: The u-plane integral is the contribution of the Coulomb branch to correlation functions of N=2 gauge theory on a compact four-manifold. We consider the u-plane integral for correlators of point and surface observables of topologically twisted theories with gauge group SU(2), for an arbitrary four-manifold with (b1,b2+)=(0,1). The u-plane contribution equals the full correlator in the absence of Seiberg-Witten contributions at strong coupling, and coincides with the mathematically defined Donaldson invariants in such cases. We demonstrate that the u-plane correlators are efficiently determined using mock modular forms for point observables, and Appell-Lerch sums for surface observables. We use these results to discuss the asymptotic behavior of correlators as function of the number of observables. Our findings suggest that the vev of exponentiated point and surface observables is an entire function of the fugacities.


title: Femtometer-resolved simultaneous measurement of multiple laser wavelengths in a speckle wavemeter

abstract: Many areas of optical science require an accurate measurement of optical spectra. Devices based on laser speckle promise compact wavelength measurement, with attometer-level sensitivity demonstrated for single wavelength laser fields. The measurement of multimode spectra using this approach would be attractive, yet this is currently limited to picometer resolution. Here, we present a method to improve the resolution and precision of speckle-based multi-wavelength measurements. We measure multiple wavelengths simultaneously, in a device comprising a single 1 m-long step-index multimode fiber and a fast camera. Independent wavelengths separated by as little as 1 fm are retrieved with 0.2 fm precision using Principal Component Analysis. The method offers a viable way to measure sparse spectra containing multiple individual lines and is likely to find application in the tracking of multiple lasers in fields such as portable quantum technologies and optical telecommunications.


title: Observation of spin-momentum locked surface states in amorphous Bi$_{2}$Se$_{3}$

abstract: Crystalline symmetries have played a central role in the identification of topological materials. The use of symmetry indicators and band representations have enabled a classification scheme for crystalline topological materials, leading to large scale topological materials discovery. In this work we address whether amorphous topological materials, which lie beyond this classification due to the lack of long-range structural order, exist in the solid state. We study amorphous Bi$_2$Se$_3$ thin films, which show a metallic behavior and an increased bulk resistance. The observed low field magnetoresistance due to weak antilocalization demonstrates a significant number of two dimensional surface conduction channels. Our angle-resolved photoemission spectroscopy data is consistent with a dispersive two-dimensional surface state that crosses the bulk gap. Spin resolved photoemission spectroscopy shows this state has an anti-symmetric spin texture resembling that of the surface state of crystalline Bi$_2$Se$_3$. These experimental results are consistent with theoretical photoemission spectra obtained with an amorphous tight-binding model that utilizes a realistic amorphous structure. This discovery of amorphous materials with topological properties uncovers an overlooked subset of topological matter outside the current classification scheme, enabling a new route to discover materials that can enhance the development of scalable topological devices.


title: Feature relevance quantification in explainable AI: A causal problem

abstract: We discuss promising recent contributions on quantifying feature relevance using Shapley values, where we observed some confusion on which probability distribution is the right one for dropped features. We argue that the confusion is based on not carefully distinguishing between observational and interventional conditional probabilities and try a clarification based on Pearl's seminal work on causality. We conclude that unconditional rather than conditional expectations provide the right notion of dropping features in contradiction to the theoretical justification of the software package SHAP. Parts of SHAP are unaffected because unconditional expectations (which we argue to be conceptually right) are used as approximation for the conditional ones, which encouraged others to `improve' SHAP in a way that we believe to be flawed.


title: Pronounced non-Markovian features in multiply-excited, multiple-emitter waveguide-QED: Retardation-induced anomalous population trappin

abstract: The Markovian approximation is widely applied in the field of quantum optics due to the weak frequency dependence of the vacuum field amplitude, and in consequence non-Markovian effects are typically regarded to play a minor role in the optical electron-photon interaction. Here, we give an example where non-Markovianity changes the qualitative behavior of a quantum optical system, rendering the Markovian approximation quantitatively and qualitatively insufficient. Namely, we study a multiple-emitter, multiple-excitation waveguide quantum-electrodynamic (waveguide-QED) system and include propagation time delay. In particular, we demonstrate anomalous population trapping as a result of the retardation in the excitation exchange between the waveguide and three initially excited emitters. Allowing for local phases in the emitter-waveguide coupling, this population trapping cannot be recovered using a Markovian treatment, proving the essential role of non-Markovian dynamics in the process. Furthermore, this time-delayed excitation exchange allows for a novel steady state, in which one emitter decays entirely to its ground state while the other two remain partially excited.


title: Analyzing Hack Subnetworks in the Bitcoin Transaction Graph

abstract: Hacks are one of the most damaging types of cryptocurrency related crime, accounting for billions of dollars in stolen funds since 2009. Professional investigators at Chainalysis have traced these stolen funds from the initial breach on an exchange to off-ramps, i.e. services where criminals are able to convert the stolen funds into fiat or other cryptocurrencies. We analyzed six hack subnetworks of bitcoin transactions known to belong to two prominent hacking groups. We analyze each hack according to eight network features, both static and temporal, and successfully classify each hack to its respective hacking group through our newly proposed method. We find that the static features, such as node balance, in degree, and out degree are not as useful in classifying the hacks into hacking groups as temporal features related to how quickly the criminals cash out. We validate our operating hypothesis that the key distinction between the two hacking groups is the acceleration with which the funds exit through terminal nodes in the subnetworks.


title: Spectroscopy of short-lived radioactive molecules: A sensitive laboratory for new physics

abstract: The study of molecular systems provides exceptional opportunities for the exploration of the fundamental laws of nature and for the search for physics beyond the Standard Model of particle physics. Measurements of molecules composed of naturally occurring nuclei have provided the most stringent upper bounds to the electron electric dipole moment to date, and offer a route to investigate the violation of fundamental symmetries with unprecedented sensitivity. Radioactive molecules - where one or more of their atoms possesses a radioactive nucleus - can contain heavy and deformed nuclei, offering superior sensitivity for EDM measurements as well as for other symmetry-violating effects. Radium monofluoride, RaF, is of particular interest as it is predicted to have an appropriate electronic structure for direct laser cooling. Furthermore, some Ra isotopes are known to be octupole deformed, thereby resulting in a large enhancement of their symmetry-violating nuclear moments. Until now,however, no experimental measurements of RaF have been performed, and their study is impeded by major experimental challenges, as no stable isotopes of radium exist. Here, we present a novel experimental approach to study short-lived radioactive molecules using the highly sensitive collinear resonance ionisation method. With this technique we have measured, for the first time, the energetically low-lying electronic states for each of the isotopically pure RaF molecules at the ISOLDE-CERN. Our results provide strong evidence of the existence of a suitable laser-cooling scheme for these molecules and constitute a pivotal step towards high-precision studies in these systems. Our findings open up new opportunities in the synthesis, manipulation and study of short-lived radioactive molecules, which will have a direct impact in many-body physics, astrophysics, nuclear structure, and fundamental physics research.


title: Lifting bicategories through the Grothendieck construction

abstract: We treat the problem of lifting bicategories into double categories through categories of vertical morphisms. We make use of a specific instance of the Grothendieck construction to provide, for every bicategory equipped with a possible vertical category, together with a suitable monoidal pre-cosheaf relating these two structures, a double category lifting the decorated bicategory along the category of vertical morphisms provided as set of initial conditions. We prove in particular that every decorated bicategory admits a lift to a double category. We study relations of instances of our construction to foldings, cofoldings, framed bicategories and globularily generated double categories.


title: Wasserstein $F$-tests and Confidence Bands for the Fr\`echet Regression of Density Response Curves

abstract: Data consisting of samples of probability density functions are increasingly prevalent, necessitating the development of methodologies for their analysis that respect the inherent nonlinearities associated with densities. In many applications, density curves appear as functional response objects in a regression model with vector predictors. For such models, inference is key to understand the importance of density-predictor relationships, and the uncertainty associated with the estimated conditional mean densities, defined as conditional Fr\'echet means under a suitable metric. Using the Wasserstein geometry of optimal transport, we consider the Fr\'echet regression of density curve responses and develop tests for global and partial effects, as well as simultaneous confidence bands for estimated conditional mean densities. The asymptotic behavior of these objects is based on underlying functional central limit theorems within Wasserstein space, and we demonstrate that they are asymptotically of the correct size and coverage, with uniformly strong consistency of the proposed tests under sequences of contiguous alternatives. The accuracy of these methods, including nominal size, power, and coverage, is assessed through simulations, and their utility is illustrated through a regression analysis of post-intracerebral hemorrhage hematoma densities and their associations with a set of clinical and radiological covariates.


title: A distributional approach to fractional Sobolev spaces and fractional variation: asymptotics I

abstract: We continue the study of the space $BV^\alpha(\mathbb{R}^n)$ of functions with bounded fractional variation in $\mathbb{R}^n$ of order $\alpha\in(0,1)$ introduced in arXiv:1809.08575, by dealing with the asymptotic behaviour of the fractional operators involved. After some technical improvements of certain results of our previous work, we prove that the fractional $\alpha$-variation converges to the standard De Giorgi's variation both pointwise and in the $\Gamma$-limit sense as $\alpha\to1^-$. We also prove that the fractional $\beta$-variation converges to the fractional $\alpha$-variation both pointwise and in the $\Gamma$-limit sense as $\beta\to\alpha^-$ for any given $\alpha\in(0,1)$.


title: Punctual Hilbert schemes for Kleinian singularities as quiver varieties

abstract: For a finite subgroup $\Gamma\subset \mathrm{SL}(2,\mathbb{C})$ and $n\geq 1$, we construct the (reduced scheme underlying the) Hilbert scheme of $n$ points on the Kleinian singularity $\mathbb{C}^2/\Gamma$ as a Nakajima quiver variety for the framed McKay quiver of $\Gamma$, taken at a specific non-generic stability parameter. We deduce that this Hilbert scheme is irreducible (a result previously due to Zheng), normal, and admits a unique symplectic resolution. More generally, we introduce a class of algebras obtained from the preprojective algebra of the framed McKay quiver by a process called cornering, and we show that fine moduli spaces of cyclic modules over these new algebras are isomorphic to quiver varieties for the framed McKay quiver and certain non-generic choices of stability parameter.


title: Linear random walks on the torus

abstract: We prove a quantitative equidistribution result for linear random walks on the torus, similar to a theorem of Bourgain, Furman, Lindenstrauss and Mozes, but without any proximality assumption.


title: On the reconstruction accuracy of multi-coil MRI with orthogonal projections

abstract: MRI signal acquisition with multiple coils in a phased array is nowadays commonplace. The use of multiple receiver coils increases the signal-to-noise ratio (SNR) and enables accelerated parallel imaging methods. Some of these methods, like GRAPPA or SPIRiT, yield individual coil images in the k-space domain which need to be combined to form a final image. Coil combination is often the last step of the image reconstruction, where the root sum of squares (rSOS) is frequently used. This straightforward method works well for coil images with high SNR, but can yield problems in images with artifacts or low SNR in all individual coils. We aim to analyze the final coil combination step in the framework of linear compression, including principal component analysis (PCA). With two data sets, a simulated and an in-vivo, we use random projections as a representation of the whole space of orthogonal projections. This allows us to study the impact of linear compression in the image space with diverse measures of reconstruction accuracy. In particular, the $L_2$ error, variance, SNR, and visual results serve as performance measures to describe the final image quality. We study their relationships and observe that the $L_2$ error and variance strongly correlate, but as expected minimal $L_2$ error does not necessarily correspond to the best visual results. In terms of visual evaluation and SNR, the compression with PCA outperforms all other methods, including rSOS on the uncompressed image space data.

categories: physics.med-ph cs.NA eess.IV eess.SP math.NA


title: Topological representations of motion groups and mapping class groups -- a unified functorial construction

abstract: For groups of a topological origin, such as braid groups and mapping class groups, an important source of interesting and highly non-trivial representations is given by their actions on the twisted homology of associated spaces; these are known as homological representations. Representations of this kind have proved themselves especially important for the question of linearity, a key example being the family of topologically-defined representations introduced by Lawrence and Bigelow, and used by Bigelow and Krammer to prove that braid groups are linear. In this paper, we give a unified foundation for the construction of homological representations using a functorial approach. Namely, we introduce homological representation functors encoding a large class of homological representations, defined on categories containing all mapping class groups and motion groups in a fixed dimension. These source categories are defined using a topological enrichment of the Quillen bracket construction applied to categories of decorated manifolds. This approach unifies many previously-known constructions, including those of Lawrence-Bigelow, and yields many new representations.


abstract: We study a critical case of Coagulation-Fragmentation equations with multiplicative coagulation kernel and constant fragmentation kernel. Our method is based on the study of viscosity solutions to a new singular Hamilton-Jacobi equation, which results from applying the Bernstein transform to the original Coagulation-Fragmentation equation. Our results include wellposedness, regularity and long-time behaviors of viscosity solutions to the Hamilton-Jacobi equation in certain regimes, which have implications to wellposedness and long-time behaviors of \emph{mass-conserving} solutions to the Coagulation-Fragmentation equation.