title: Mise en abyme with artificial intelligence: how to predict the accuracy of NN, applied to hyper-parameter tuning

abstract: In the context of deep learning, the costliest phase from a computational point of view is the full training of the learning algorithm. However, this process is to be used a significant number of times during the design of a new artificial neural network, leading therefore to extremely expensive operations. Here, we propose a low-cost strategy to predict the accuracy of the algorithm, based only on its initial behaviour. To do so, we train the network of interest up to convergence several times, modifying its characteristics at each training. The initial and final accuracies observed during this beforehand process are stored in a database. We then make use of both curve fitting and Support Vector Machines techniques, the latter being trained on the created database, to predict the accuracy of the network, given its accuracy on the primary iterations of its learning. This approach can be of particular interest when the space of the characteristics of the network is notably large or when its full training is highly time-consuming. The results we obtained are promising and encouraged us to apply this strategy to a topical issue: hyper-parameter optimisation (HO). In particular, we focused on the HO of a convolutional neural network for the classification of the databases MNIST and CIFAR-10. By using our method of prediction, and an algorithm implemented by us for a probabilistic exploration of the hyper-parameter space, we were able to find the hyper-parameter settings corresponding to the optimal accuracies already known in literature, at a quite low-cost.


title: Applications of Zalcman's lemma in $C^n$

abstract: The aim of this paper is to give some applications of Zalcman's Rescalling Lemma.


title: Blowing Up Solutions to the Zakharov System for Langmuir Waves

abstract: Langmuir waves take place in a quasi-neutral plasma and are modeled by the Zakharov system. The phenomenon of collapse, described by blowing up solutions plays a central role in their dynamics. We present in this article a review of the main mathematical properties of blowing up solutions. They include conditions for blowup in finite or infinite time, description of self-similar singular solutions and lower bounds for the rate of blowup of certain norms associated to the solutions.


title: A Unified Approach to Robust Mean Estimation

abstract: In this paper, we develop connections between two seemingly disparate, but central, models in robust statistics: Huber's epsilon-contamination model and the heavy-tailed noise model. We provide conditions under which this connection provides near-statistically-optimal estimators. Building on this connection, we provide a simple variant of recent computationally-efficient algorithms for mean estimation in Huber's model, which given our connection entails that the same efficient sample-pruning based estimators is simultaneously robust to heavy-tailed noise and Huber contamination. Furthermore, we complement our efficient algorithms with statistically-optimal albeit computationally intractable estimators, which are simultaneously optimally robust in both models. We study the empirical performance of our proposed estimators on synthetic datasets, and find that our methods convincingly outperform a variety of practical baselines.


title: Trading Off Computation with Transmission in Status Update Systems

abstract: This paper is motivated by emerging edge computing applications in which generated data are pre-processed at the source and then transmitted to an edge server. In such a scenario, there is typically a tradeoff between the amount of pre-processing and the amount of data to be transmitted. We model such a system by considering two non-preemptive queues in tandem whose service times are independent over time but the transmission service time is dependent on the computation service time in mean value. The first queue is in M/GI/1/1 form with a single server, memoryless exponential arrivals, general independent service and no extra buffer to save incoming status update packets. The second queue is in GI/M/1/2* form with a single server receiving packets from the first queue, memoryless service and a single data buffer to save incoming packets. Additionally, mean service times of the first and second queues are dependent through a deterministic monotonic function. We perform stationary distribution analysis in this system and obtain closed form expressions for average age of information (AoI) and average peak AoI. Our numerical results illustrate the analytical findings and highlight the tradeoff between average AoI and average peak AoI generated by the tandem nature of the queueing system with dependent service times.


title: Trimming Graphs Using Clausal Proof Optimization

abstract: We present a method to gradually compute a smaller and smaller unsatisfiable core of a propositional formula by minimizing proofs of unsatisfiability. The goal is to compute a minimal unsatisfiable core that is relatively small compared to other minimal unsatisfiable cores of the same formula. We try to achieve this goal by postponing deletion of arbitrary clauses from the formula as long as possible---in contrast to existing minimal unsatisfiable core algorithms. We applied this method to reduce the smallest known unit-distance graph with chromatic number 5 from 553 vertices and 2720 edges to 529 vertices and 2670 edges.


title: A Joint Optimization Approach of LiDAR-Camera Fusion for Accurate Dense 3D Reconstructions

abstract: Fusing data from LiDAR and camera is conceptually attractive because of their complementary properties. For instance, camera images are higher resolution and have colors, while LiDAR data provide more accurate range measurements and have a wider Field Of View (FOV). However, the sensor fusion problem remains challenging since it is difficult to find reliable correlations between data of very different characteristics (geometry vs. texture, sparse vs. dense). This paper proposes an offline LiDAR-camera fusion method to build dense, accurate 3D models. Specifically, our method jointly solves a bundle adjustment (BA) problem and a cloud registration problem to compute camera poses and the sensor extrinsic calibration. In experiments, we show that our method can achieve an averaged accuracy of 2.7mm and resolution of 70 points per square cm by comparing to the ground truth data from a survey scanner. Furthermore, the extrinsic calibration result is discussed and shown to outperform the state-of-the-art method.


title: Chemical potential of an antiferromagnetic magnon gas

abstract: Understanding the statistics of quasi-particle excitations in magnetic systems is essential for exploring new magnetic phases and collective quantum phenomena. While the chemical potential of a ferromagnetic gas has been extensively investigated both theoretically and experimentally, its antiferromagnetic counterpart remains uncharted. Here, we derive the statistics of a two-component U(1)-symmetric Bose gas and apply our results to an axially-symmetric antiferromagnetic insulator. We find that the two magnon eigenmodes of the system are described by an equal and opposite chemical potential, in analogy with a particle-antiparticle pair. Furthermore, we derive the thermomagnonic torques describing the interaction between the coherent and incoherent antiferromagnetic spin dynamics. Our results show that the magnitude and sign of the chemical potential can be tuned via an AC magnetic field driving resonantly one of the magnon modes. Finally, we propose NV-center relaxometry as a method to experimentally test our predictions.


title: A Framework For Identifying Group Behavior Of Wild Animals

abstract: Activity recognition and, more generally, behavior inference tasks are gaining a lot of interest. Much of it is work in the context of human behavior. New available tracking technologies for wild animals are generating datasets that indirectly may provide information about animal behavior. In this work, we propose a method for classifying these data into behavioral annotation, particularly collective behavior of a social group. Our method is based on sequence analysis with a direct encoding of the interactions of a group of wild animals. We evaluate our approach on a real world dataset, showing significant accuracy improvements over baseline methods.


title: Biased permutative equivariant categories

abstract: For a finite group G, we introduce the complete suboperad $Q_G$ of the categorical G-Barratt-Eccles operad $P_G$. We prove that $P_G$ is not finitely generated, but $Q_G$ is finitely generated and is a genuine $E_\infty$ G-operad (i.e., it is $N_\infty$ and includes all norms). For G cyclic of order 2 or 3, we determine presentations of the object operad of $Q_G$ and conclude with a discussion of algebras over $Q_G$, which we call biased permutative equivariant categories.


title: Ich wei{\ss}, was du n\"achsten Sommer getan haben wirst: Predictive Policing in \"Osterreich

abstract: Predictive policing is a data-based, predictive analytical technique used in law enforcement. In this paper, we give an overview of the current situation in Austria and discuss technical, sociopolitical and legal questions raised by the use of PP, such as the lack of awareness of discriminatory structures in society, the biases in data underlying PP and the lack of reflection on the basic premises and feedback mechanisms of PP. Violations of fundamental rights without cause are not allowed by the Austrian Code of Criminal Procedure (Strafproze{\ss}ordnung, StPO), the Security Police Act (Sicherheitspolizeigesetz, SPG) or the Act concerning Police Protection of the State (Polizeiliches Staatsschutzgesetz, PStSG); the principle of allowing police intervention only on the basis of concrete threats or suspicion must remain absolute. Considering the numerous problems (not least from the point of view of legal policy), we conclude that the use of PP should be eschewed and that resources and planning should instead be focussed on solving the social problems which actually cause crime. ----- Predictive Policing ist ein datenbasiertes und prognosegetriebenes Modell f\"ur Polizeiarbeit. Wir geben in diesem Artikel einen \"Uberblick \"uber den aktuellen Stand in \"Osterreich und diskutieren technische, politisch-gesellschaftliche und rechtliche Probleme, die sich daraus ergeben -- etwa das mangelhafte Bewusstsein f\"ur Prozesse gesellschaftlicher Diskriminierung, die verzerrte Datenbasis, die PP zugrundeliegt, und fehlende Reflexion \"uber zugrundeliegende Annahmen und R\"uckkopplungseffekte. Anlasslose Grundrechtseingriffe sind weder durch die StPO noch das SPG oder das PStSG gedeckt; dem Grundgedanken, dass Polizei erst bei konkreter Gefahrenlage oder Tatverdacht t\"atig werden darf, muss weiterhin Rechnung getragen werden. Aus unserer Sicht sollte angesichts der zahlreichen Probleme (und auch aus rechtspolitischen Erw\"agungen) auf PP verzichtet werden und stattdessen Ressourcen und \"Uberlegung in die L\"osung jener gesellschaftlicher Probleme investiert werden, die zu Kriminalit\"at f\"uhren.


title: One-Time Programs made Practical

abstract: A one-time program (OTP) works as follows: Alice provides Bob with the implementation of some function. Bob can have the function evaluated exclusively on a single input of his choosing. Once executed, the program will fail to evaluate on any other input. State-of-the-art one-time programs have remained theoretical, requiring custom hardware that is cost-ineffective/unavailable, or confined to adhoc/unrealistic assumptions. To bridge this gap, we explore how the Trusted Execution Environment (TEE) of modern CPUs can realize the OTP functionality. Specifically, we build two flavours of such a system: in the first, the TEE directly enforces the one-timeness of the program; in the second, the program is represented with a garbled circuit and the TEE ensures Bob's input can only be wired into the circuit once, equivalent to a smaller cryptographic primitive called one-time memory. These have different performance profiles: the first is best when Alice's input is small and Bob's is large, and the second for the converse.


title: PCN calculations for Z=111 to Z=118

abstract: In previous publications, we presented evidence for the importance of spin in determining capture and evaporation residue cross sections in the synthesis of heavy nuclei. We extend the previous calculations which dealt with nuclei where ZCN is less than 110 to the region of ZCN = 111-118. We deduce a new systematics of the fusion probability for these reactions


title: Semantic Product Search

abstract: We study the problem of semantic matching in product search, that is, given a customer query, retrieve all semantically related products from the catalog. Pure lexical matching via an inverted index falls short in this respect due to several factors: a) lack of understanding of hypernyms, synonyms, and antonyms, b) fragility to morphological variants (e.g. "woman" vs. "women"), and c) sensitivity to spelling errors. To address these issues, we train a deep learning model for semantic matching using customer behavior data. Much of the recent work on large-scale semantic search using deep learning focuses on ranking for web search. In contrast, semantic matching for product search presents several novel challenges, which we elucidate in this paper. We address these challenges by a) developing a new loss function that has an inbuilt threshold to differentiate between random negative examples, impressed but not purchased examples, and positive examples (purchased items), b) using average pooling in conjunction with n-grams to capture short-range linguistic patterns, c) using hashing to handle out of vocabulary tokens, and d) using a model parallel training architecture to scale across 8 GPUs. We present compelling offline results that demonstrate at least 4.7% improvement in Recall@100 and 14.5% improvement in mean average precision (MAP) over baseline state-of-the-art semantic search methods using the same tokenization method. Moreover, we present results and discuss learnings from online A/B tests which demonstrate the efficacy of our method.


title: Dynamical confirmation of a black hole in MAXI J1820+070

abstract: We present time-resolved 10.4-m GTC and 4.2-m WHT intermediate resolution spectroscopy of the X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained during its decline to the quiescent state. Cross-correlation of the 21 individual spectra against late-type templates reveals a sinusoidal velocity modulation with a period of 0.68549 +/- 0.00001 d and semi-amplitude of 417.7 +/- 3.9 km/s. We derive a mass function f(M) = 5.18 +/- 0.15 Msun, dynamically confirming the black hole nature of the compact object. Our analysis of the stellar absorption features supports a K3-5 spectral classification for the donor star, which contributes 20% of the total flux at 5200-6800 Angs. The photometric 0.703 +/- 0.003 d periodicity observed during outburst is 2.6% longer than the orbital period supporting the presence of a superhump modulation in the outburst lightcurves. In line with this interpretation, we constrain the binary mass ratio to be q=0.12. In addition, we observe a sharp increase in the Halpha emission line equivalent width during inferior conjunction of the donor star that we interpret as a grazing eclipse of the accretion disc and allows us to constrain the binary inclination to > 69 deg. On the other hand, the absence of X-ray eclipses during outburst imply i < 77 deg. These inclination limits, together with our dynamical solution, lead to a black hole mass in the range 7-8 Msun. We also measure a systemic velocity = -21.6 +/- 2.3 km/s which, combined with the Gaia DR2 proper motion and parallax, implies a large peculiar velocity of 100 km/s.


title: Pano Popups: Indoor 3D Reconstruction with a Plane-Aware Network

abstract: In this work we present a method to train a plane-aware convolutional neural network for dense depth and surface normal estimation as well as plane boundaries from a single indoor $360^\circ$ image. Using our proposed loss function, our network outperforms existing methods for single-view, indoor, omnidirectional depth estimation and provides an initial benchmark for surface normal prediction from $360^\circ$ images. Our improvements are due to the use of a novel plane-aware loss that leverages principal curvature as an indicator of planar boundaries. We also show that including geodesic coordinate maps as network priors provides a significant boost in surface normal prediction accuracy. Finally, we demonstrate how we can combine our network's outputs to generate high quality 3D "pop-up" models of indoor scenes.


title: Charge Smoothening and Band Flattening due to Hartree corrections in Twisted Bilayer Graphene

abstract: Doping twisted bilayer graphene away from charge neutrality leads to an enormous buildup of charge inhomogeneities within each Moir\'{e} unit cell. Here we show, using unbiased real-space self-consistent Hartree calculations on a relaxed lattice, that Coulomb interactions smoothen this charge imbalance by changing the occupation of earlier identified `ring' orbitals in the AB/BA region and `center' orbitals at the AA region. For hole doping, this implies an increase of the energy of the states at the ${\bf \Gamma}$ point, leading to a further flattening of the flat bands and a pinning of the Van Hove singularity at the Fermi level. The charge smoothening will affect the subtle competition between different possible correlated phases.


title: Global Pixel Transformers for Virtual Staining of Microscopy Images

abstract: Visualizing the details of different cellular structures is of great importance to elucidate cellular functions. However, it is challenging to obtain high quality images of different structures directly due to complex cellular environments. Fluorescence staining is a popular technique to label different structures but has several drawbacks. In particular, label staining is time consuming and may affect cell morphology, and simultaneous labels are inherently limited. This raises the need of building computational models to learn relationships between unlabeled microscopy images and labeled fluorescence images, and to infer fluorescence labels of other microscopy images excluding the physical staining process. We propose to develop a novel deep model for virtual staining of unlabeled microscopy images. We first propose a novel network layer, known as the global pixel transformer layer, that fuses global information from inputs effectively. The proposed global pixel transformer layer can generate outputs with arbitrary dimensions, and can be employed for all the regular, down-sampling, and up-sampling operators. We then incorporate our proposed global pixel transformer layers and dense blocks to build an U-Net like network. We believe such a design can promote feature reusing between layers. In addition, we propose a multi-scale input strategy to encourage networks to capture features at different scales. We conduct evaluations across various fluorescence image prediction tasks to demonstrate the effectiveness of our approach. Both quantitative and qualitative results show that our method outperforms the state-of-the-art approach significantly. It is also shown that our proposed global pixel transformer layer is useful to improve the fluorescence image prediction results.


title: Cryogenic Memory Architecture Integrating Spin Hall Effect based Magnetic Memory and Superconductive Cryotron Devices

abstract: One of the most challenging obstacles to realizing exascale computing is minimizing the energy consumption of L2 cache, main memory, and interconnects to that memory. For promising cryogenic computing schemes utilizing Josephson junction superconducting logic, this obstacle is exacerbated by the cryogenic system requirements that expose the technology's lack of high-density, high-speed and power-efficient memory. Here we demonstrate an array of cryogenic memory cells consisting of a non-volatile three-terminal magnetic tunnel junction element driven by the spin Hall effect, combined with a superconducting heater-cryotron bit-select element. The write energy of these memory elements is roughly 8 pJ with a bit-select element, designed to achieve a minimum overhead power consumption of about 30%. Individual magnetic memory cells measured at 4 K show reliable switching with write error rates below $10^{-6}$, and a 4x4 array can be fully addressed with bit select error rates of $10^{-6}$. This demonstration is a first step towards a full cryogenic memory architecture targeting energy and performance specifications appropriate for applications in superconducting high performance and quantum computing control systems, which require significant memory resources operating at 4 K.


title: Estimating brain age based on a healthy population with deep learning and structural MRI

abstract: Numerous studies have established that estimated brain age, as derived from statistical models trained on healthy populations, constitutes a valuable biomarker that is predictive of cognitive decline and various neurological diseases. In this work, we curate a large-scale heterogeneous dataset (N = 10,158, age range 18 - 97) of structural brain MRIs in a healthy population from multiple publicly-available sources, upon which we train a deep learning model for brain age estimation. The availability of the large-scale dataset enables a more uniform age distribution across adult life-span for effective age estimation with no bias toward certain age groups. We demonstrate that the age estimation accuracy, evaluated with mean absolute error (MAE) and correlation coefficient (r), outperforms previously reported methods in both a hold-out test set reflective of the custom population (MAE = 4.06 years, r = 0.970) and an independent life-span evaluation dataset (MAE = 4.21 years, r = 0.960) on which a previous study has evaluated. We further demonstrate the utility of the estimated age in life-span aging analysis of cognitive functions. Furthermore, we conduct extensive ablation tests and employ feature-attribution techniques to analyze which regions contribute the most predictive value, demonstrating the prominence of the frontal lobe as well as pattern shift across life-span. In summary, we achieve superior age estimation performance confirming the efficacy of deep learning and the added utility of training with data both in larger number and more uniformly distributed than in previous studies. We demonstrate the regional contribution to our brain age predictions through multiple routes and confirm the association of divergence between estimated and chronological brain age with neuropsychological measures.


title: Extremes and Records

abstract: These are lecture notes from a course offered at the Bangalore School on Statistical Physics - X, during 17-28 June 2019, [ https://www.icts.res.in/program/bssp2019 ] at International centre of theoretical physics (ICTS), Bangalore. These pedagogical lectures are at the introductory level, intended mainly for master/Ph.D. students or researchers from outside the field. In these lectures, we discuss about the limit laws for the sample mean and the maximum of a set of independent and identically distributed (i.i.d.) random variables as well as random walks / Brownian motion. The density of near-extreme events is also discussed. Finally, we discuss the statistics of records for an i.i.d. random sequence as well as random walks in discrete and continuous time. Some exercises are provided for the students to work out. The video recording of the lectures are available at https://www.icts.res.in/program/bssp2019/talks


title: ICDAR2019 Robust Reading Challenge on Multi-lingual Scene Text Detection and Recognition -- RRC-MLT-2019

abstract: With the growing cosmopolitan culture of modern cities, the need of robust Multi-Lingual scene Text (MLT) detection and recognition systems has never been more immense. With the goal to systematically benchmark and push the state-of-the-art forward, the proposed competition builds on top of the RRC-MLT-2017 with an additional end-to-end task, an additional language in the real images dataset, a large scale multi-lingual synthetic dataset to assist the training, and a baseline End-to-End recognition method. The real dataset consists of 20,000 images containing text from 10 languages. The challenge has 4 tasks covering various aspects of multi-lingual scene text: (a) text detection, (b) cropped word script classification, (c) joint text detection and script classification and (d) end-to-end detection and recognition. In total, the competition received 60 submissions from the research and industrial communities. This paper presents the dataset, the tasks and the findings of the presented RRC-MLT-2019 challenge.


title: Uniserial Noetherian Centrally Essential Rings

abstract: It is proved that a ring $R$ is a right uniserial, right Noetherian centrally essential ring if and only if $R$ is a commutative discrete valuation domain or a left and right Artinian, left and right uniserial ring. It is also proved that there exist non-commutative uniserial Artinian centrally essential rings. Victor Markov is supported by the Russian Foundation for Basic Research, project 17-01-00895-A. Askar Tuganbaev is supported by Russian Scientific Foundation, project 16-11-10013.


title: Automorphism-Liftable Modules

abstract: In this paper, we describe all automorphism-liftable torsion modules over non-primitive hereditary Noetherian prime rings. We also study automorphism-liftable non-torsion modules over not necessarily commutative Dedekind prime rings


title: Spin Diffusion Equations for Magnetized or Orbital Polarized Systems

abstract: Charge and spin transport in spintronics devices can be described by a spin diffusion equation suitable for modelling scales much larger than the scattering and atomic scales. This work concerns the coarse graining procedure used to compute the coefficients of the diffusion equation, which are sensitive the details of individual atoms and impurities. We show with two simple examples that in spintronics devices which have both a spin-orbit interaction and magnetization, standard coarse graining can easily obtain diffusion equations which fail to conserve electronic charge. The same failure can occur in systems with both a spin-orbit interaction and orbital polarization. We show that linear response theory, coupled with the self-consistent Born approximation and ladder diagrams, offers an improved way of calculating diffusion equations. We show that the resulting equations satisfy a Ward-Takahashi identity that guarantees charge conservation.


title: Optimization on flag manifolds

abstract: A flag is a sequence of nested subspaces. Flags are ubiquitous in numerical analysis, arising in finite elements, multigrid, spectral, and pseudospectral methods for numerical PDE; they arise in the form of Krylov subspaces in matrix computations, and as multiresolution analysis in wavelets constructions. They are common in statistics too --- principal component, canonical correlation, and correspondence analyses may all be viewed as methods for extracting flags from a data set. The main goal of this article is to develop the tools needed for optimizing over a set of flags, which is a smooth manifold called the flag manifold, and it contains the Grassmannian as the simplest special case. We will derive closed-form analytic expressions for various differential geometric objects required for Riemannian optimization algorithms on the flag manifold; introducing various systems of extrinsic coordinates that allow us to parameterize points, metrics, tangent spaces, geodesics, distance, parallel transport, gradients, Hessians in terms of matrices and matrix operations; and thereby permitting us to formulate steepest descent, conjugate gradient, and Newton algorithms on the flag manifold using only standard numerical linear algebra.


title: Explaining the Human Visual Brain Challenge 2019 -- receptive fields and surrogate features

abstract: In this paper I review the submission to the Explaining the Human Visual Brain Challenge 2019 in both the fMRI and MEG tracks. The goal was to construct neural network features which generate the so-called representational dissimilarity matrix (RDM) which is most similar to the one extracted from fMRI and MEG data upon viewing a set of images. I review exploring the optimal granularity of the receptive field, a construction of intermediate surrogate features using Multidimensional Scaling and modelling them using neural network features. I also point out some peculiarities of the RDM construction which have to be taken into account.


title: Colength, multiplicity, and ideal closure operations

abstract: In a formally unmixed Noetherian local ring, if the colength and multiplicity of an integrally closed ideal agree, then $R$ is regular. We deduce this using the relationship between multiplicity and various ideal closure operations.


title: Coexistence of ultraheavy and ultrarelativistic Dirac quasiparticles in sandwiched trilayer graphene

abstract: Electrons in quantum materials exhibiting coexistence of dispersionless (flat) bands piercing dispersive (steep) bands can give rise to strongly correlated phenomena, and are associated with unconventional superconductivity. It is known that in twisted trilayer graphene steep Dirac cones can coexist with band flattening, but the phenomenon is not stable under layer misalignments. Here we show that such a twisted sandwiched graphene (TSWG) -- a three-layer van der Waals heterostructure with a twisted middle layer -- can have very stable flat bands coexisting with Dirac cones near the Fermi energy when twisted to 1.5 degrees. These flat bands require a specific high-symmetry stacking order, and our atomistic calculations predict that TSWG always relaxes to it. Additionally, with external fields, we can control the relative energy offset between the Dirac cone vertex and the flat bands. Our work establishes twisted sandwiched graphene as a new platform for research into strongly interacting phases, and topological transport beyond Dirac and Weyl semimetals.


title: Stochastic Latent Actor-Critic: Deep Reinforcement Learning with a Latent Variable Model

abstract: Deep reinforcement learning (RL) algorithms can use high-capacity deep networks to learn directly from image observations. However, these high-dimensional observation spaces present a number of challenges in practice, since the policy must now solve two problems: representation learning and task learning. In this work, we tackle these two problems separately, by explicitly learning latent representations that can accelerate reinforcement learning from images. We propose the stochastic latent actor-critic (SLAC) algorithm: a sample-efficient and high-performing RL algorithm for learning policies for complex continuous control tasks directly from high-dimensional image inputs. SLAC provides a novel and principled approach for unifying stochastic sequential models and RL into a single method, by learning a compact latent representation and then performing RL in the model's learned latent space. Our experimental evaluation demonstrates that our method outperforms both model-free and model-based alternatives in terms of final performance and sample efficiency, on a range of difficult image-based control tasks. Our code and videos of our results are available at our website.


title: Observation of the $\Lambda_b^0\rightarrow \chi_{c1}(3872)pK^-$ decay

abstract: Using proton-proton collision data, collected with the LHCb detector and corresponding to 1.0, 2.0 and 1.9fb$^{-1}$ of integrated luminosity at the centre-of-mass energies of 7, 8, and 13 TeV, respectively, the decay $\Lambda_b^0\to \chi_{c1}(3872)pK^-$ with $\chi_{c1}\to J/\psi\pi^+\pi^-$ is observed for the first time. The significance of the observed signal is in excess of seven standard deviations. It is found that $(58\pm15)\%$ of the decays proceed via the two-body intermediate state $\chi_{c1}(3872)\Lambda(1520)$. The~branching fraction with respect to that of the $\Lambda_b\rightarrow\psi(2S)p K^{-}$ decay mode, where the $\psi(2S)$~meson is reconstructed in the $J/\psi \pi^+\pi^-$ final state, is measured to be: \begin{equation*} \frac{\Lambda_b^0\to\chi_{c1}(3872)pK^-}{\Lambda_b\to\psi(2S)p K^-} \times \frac{\mathcal{B}(\chi_{c1} \to J/\psi \pi^+\pi^-)}{\mathcal{B}(\psi(2S)\to J/\psi \pi^+\pi^-)} = \left(5.4 \pm 1.1 \pm 0.2\right)\times 10^{-2}\,, \end{equation*} where the first uncertainty is statistical and the second is systematic.


title: Extremal horizons stationary to the second order: new constraints

abstract: We consider non-expanding shear free (NE-SF) null surface geometries embeddable as extremal Killing horizons to the second order in Einstein vacuum spacetimes. A NE-SF null surface geometry consists of a degenerate metric tensor and a consistent torsion free covariant derivative. We derive the constraints implied by the existence of an embedding. The first constraint is well known as the near horizon geometry equation. The second constraint we find is new. The constraints lead to a complete characterization of those NE-SF null geometries that are embeddable in the extremal Kerr spacetime. Our results are also valid for spacetimes with a cosmological constant.


title: Fast Uniform Dispersion of a Crash-prone Swarm

abstract: We consider the problem of completely covering an unknown discrete environment with a swarm of asynchronous, frequently-crashing autonomous mobile robots. We represent the environment by a discrete graph, and task the robots with occupying every vertex and with constructing an implicit distributed spanning tree of the graph. The robotic agents activate independently at random exponential waiting times of mean $1$ and enter the graph environment over time from a source location. They grow the environment's coverage by 'settling' at empty locations and aiding other robots' navigation from these locations. The robots are identical and make decisions driven by the same simple and local rule of behaviour. The local rule is based only on the presence of neighbouring robots, and on whether a settled robot points to the current location. Whenever a robot moves, it may crash and disappear from the environment. Each vertex in the environment has limited physical space, so robots frequently obstruct each other. Our goal is to show that even under conditions of asynchronicity, frequent crashing, and limited physical space, the simple mobile robots complete their mission in linear time asymptotically almost surely, and time to completion degrades gracefully with the frequency of the crashes. Our model and analysis are based on the well-studied "totally asymmetric simple exclusion process" in statistical mechanics.


title: Magnetic skyrmion artificial synapse for neuromorphic computing

abstract: Since the experimental discovery of magnetic skyrmions achieved one decade ago, there have been significant efforts to bring the virtual particles into all-electrical fully functional devices, inspired by their fascinating physical and topological properties suitable for future low-power electronics. Here, we experimentally demonstrate such a device: electrically-operating skyrmion-based artificial synaptic device designed for neuromorphic computing. We present that controlled current-induced creation, motion, detection and deletion of skyrmions in ferrimagnetic multilayers can be harnessed in a single device at room temperature to imitate the behaviors of biological synapses. Using simulations, we demonstrate that such skyrmion-based synapses could be used to perform neuromorphic pattern-recognition computing using handwritten recognition data set, reaching to the accuracy of ~89 percents, comparable to the software-based training accuracy of ~94 percents. Chip-level simulation then highlights the potential of skyrmion synapse compared to existing technologies. Our findings experimentally illustrate the basic concepts of skyrmion-based fully functional electronic devices while providing a new building block in the emerging field of spintronics-based bio-inspired computing.


title: Baryon as a Quantum Hall Droplet and the Cheshire Cat Principle

abstract: We show that the recent proposal to describe the $N_f=1$ baryon in the large number of color limit as a quantum Hall droplet, can be understood as a chiral bag in a 1+2 dimensional strip using the Cheshire cat principle. For a small bag radius, the bag reduces to a vortex line which is the smile of the cat with flowing gapless quarks all spinning in the same direction. The disc enclosed by the smile is described by a topological field theory due to the Callan-Harvey anomaly out-flow. The chiral bag carries naturally unit baryon number and spin $\frac 12 N_c$. The generalization to arbitrary $N_f$ is discussed.


title: Single-Path Mobile AutoML: Efficient ConvNet Design and NAS Hyperparameter Optimization

abstract: Can we reduce the search cost of Neural Architecture Search (NAS) from days down to only few hours? NAS methods automate the design of Convolutional Networks (ConvNets) under hardware constraints and they have emerged as key components of AutoML frameworks. However, the NAS problem remains challenging due to the combinatorially large design space and the significant search time (at least 200 GPU-hours). In this work, we alleviate the NAS search cost down to less than 3 hours, while achieving state-of-the-art image classification results under mobile latency constraints. We propose a novel differentiable NAS formulation, namely Single-Path NAS, that uses one single-path over-parameterized ConvNet to encode all architectural decisions based on shared convolutional kernel parameters, hence drastically decreasing the search overhead. Single-Path NAS achieves state-of-the-art top-1 ImageNet accuracy (75.62%), hence outperforming existing mobile NAS methods in similar latency settings (~80ms). In particular, we enhance the accuracy-runtime trade-off in differentiable NAS by treating the Squeeze-and-Excitation path as a fully searchable operation with our novel single-path encoding. Our method has an overall cost of only 8 epochs (24 TPU-hours), which is up to 5,000x faster compared to prior work. Moreover, we study how different NAS formulation choices affect the performance of the designed ConvNets. Furthermore, we exploit the efficiency of our method to answer an interesting question: instead of empirically tuning the hyperparameters of the NAS solver (as in prior work), can we automatically find the hyperparameter values that yield the desired accuracy-runtime trade-off? We open-source our entire codebase at: https://github.com/dstamoulis/single-path-nas.


title: Going Deeper with Lean Point Networks

abstract: In this work we introduce Lean Point Networks (LPNs) to train deeper and more accurate point processing networks by relying on three novel point processing blocks that improve memory consumption, inference time, and accuracy: a convolution-type block for point sets that blends neighborhood information in a memory-efficient manner; a crosslink block that efficiently shares information across low- and high-resolution processing branches; and a multiresolution point cloud processing block for faster diffusion of information. By combining these blocks, we design wider and deeper point-based architectures. We report systematic accuracy and memory consumption improvements on multiple publicly available segmentation tasks by using our generic modules as drop-in replacements for the blocks of multiple architectures (PointNet++, DGCNN, SpiderNet, PointCNN).


title: On the development of symmetry-preserving finite element schemes for ordinary differential equations

abstract: In this paper we introduce a procedure, based on the method of equivariant moving frames, for formulating continuous Galerkin finite element schemes that preserve the Lie point symmetries of initial value problems for ordinary differential equations. Our methodology applies to projectable and non-projectable symmetry group actions, to ordinary differential equations of arbitrary order, and finite element approximations of arbitrary polynomial degree. Several examples are included to illustrate various features of the symmetry-preserving process. We summarise extensive numerical experiments showing that symmetry-preserving finite element schemes may provide better long term accuracy than their non-invariant counterparts and can be implemented on larger elements.


title: Claim Extraction in Biomedical Publications using Deep Discourse Model and Transfer Learning

abstract: Claims are a fundamental unit of scientific discourse. The exponential growth in the number of scientific publications makes automatic claim extraction an important problem for researchers who are overwhelmed by this information overload. Such an automated claim extraction system is useful for both manual and programmatic exploration of scientific knowledge. In this paper, we introduce a new dataset of 1,500 scientific abstracts from the biomedical domain with expert annotations for each sentence indicating whether the sentence presents a scientific claim. We introduce a new model for claim extraction and compare it to several baseline models including rule-based and deep learning techniques. Moreover, we show that using a transfer learning approach with a fine-tuning step allows us to improve performance from a large discourse-annotated dataset. Our final model increases F1-score by over 14 percent points compared to a baseline model without transfer learning. We release a publicly accessible tool for discourse and claims prediction along with an annotation tool. We discuss further applications beyond biomedical literature.


title: Limit theorems for additive functionals of continuous time random walks

abstract: For a continuous-time random walk $X=\{X_t,t\ge 0\}$ (in general non-Markov), we study the asymptotic behavior, as $t\rightarrow \infty$, of the normalized additive functional $c_t\int_0^{t} f(X_s)ds$, $t\ge 0$. Similarly to the Markov situation, assuming that the distribution of jumps of $X$ belongs to the domain of attraction to $\alpha$-stable law with $\alpha>1$, we establish the convergence to the local time at zero of an $\alpha$-stable L\'evy motion. We further study a situation where $X$ is delayed by a random environment given by the Poisson shot-noise potential: $\Lambda(x,\gamma)= e^{-\sum_{y\in \gamma} \phi(x-y)},$ where $\phi\colon\mathbb R\to [0,\infty)$ is a bounded function decaying sufficiently fast, and $\gamma$ is a homogeneous Poisson point process, independent of $X$. We find that in this case the weak limit has both "quenched" component depending on $\Lambda$, and a component, where $\Lambda$ is "averaged".


title: Tur\'an theorems for unavoidable patterns

abstract: We prove Tur\'an-type theorems for two related Ramsey problems raised by Bollob\'as and by Fox and Sudakov. First, for $t \ge 3$, we show that any two-colouring of the complete graph on $n$ vertices that is $\delta$-far from being monochromatic contains an \emph{unavoidable $t$-colouring} when $\delta \gg n^{-1/t}$, where an unavoidable $t$-colouring is any two-colouring of a clique of order $2t$ in which one colour forms either a clique of order $t$ or two disjoint cliques of order $t$. Next, for $ t\ge 3$, we show that any tournament on $n$ vertices that is $\delta$-far from being transitive contains an \emph{unavoidable $t$-tournament} when $\delta \gg n^{-1/\lceil t/2 \rceil}$, where an unavoidable $t$-tournament is the blow-up of a cyclic triangle obtained by replacing each vertex of the triangle by a transitive tournament of order $t$. Conditional on a well-known conjecture about bipartite Tur\'an numbers, both results are sharp up to implied constants and hence determine the order of magnitude of the corresponding off-diagonal Ramsey numbers.


title: On-demand thermoelectric generation of equal-spin Cooper pairs

abstract: Superconducting spintronics is based on the creation of spin-triplet Cooper pairs in ferromagnet-superconductor (F-S) hybrid junctions. Previous proposals to manipulate spin-polarized supercurrents on-demand typically require the ability to carefully control magnetic materials. We, instead, propose a quantum heat engine that generates equal-spin Cooper pairs and drives supercurrents on-demand without manipulating magnetic components. We consider a S-F-S junction, connecting two leads at different temperatures, on top of the helical edge of a two-dimensional topological insulator. Heat and charge currents generated by the thermal bias are caused by different transport processes, where electron cotunneling is responsible for the heat flow to the cold lead and, strikingly, only crossed Andreev reflections contribute to the charge current. Such a purely nonlocal Andreev thermoelectric effect injects spin-polarized Cooper pairs at the superconductors, generating a supercurrent that can be switched on/off by tuning their relative phase. We further demonstrate that signatures of spin-triplet pairing are facilitated by rather low fluctuations of the thermoelectric current for temperature gradients smaller than the superconducting gap.


title: Shock acceleration efficiency in radio relics

abstract: Radio relics in galaxy clusters are giant diffuse synchrotron sources powered in cluster outskirts by merger shocks. Although the relic-shock connection has been consolidated in recent years by a number of observations, the details of the mechanisms leading to the formation of relativistic particles in this environment are still not well understood. The diffusive shock acceleration (DSA) theory is a commonly adopted scenario to explain the origin of cosmic rays at astrophysical shocks, including those in radio relics in galaxy clusters. However, in a few specific cases it has been shown that the energy dissipated by cluster shocks is not enough to reproduce the luminosity of the relics via DSA of thermal particles. Studies based on samples of radio relics are required to further address this limitation of the mechanism. In this paper, we focus on ten well-studied radio relics with underlying shocks observed in the X-rays and calculate the electron acceleration efficiency of these shocks that is necessary to reproduce the observed radio luminosity of the relics. We find that in general the standard DSA cannot explain the origin of the relics if electrons are accelerated from the thermal pool with an efficiency significantly smaller than 10%. Our results show that other mechanisms, such as shock re-acceleration of supra-thermal seed electrons or a modification of standard DSA, are required to explain the formation of radio relics.


title: Autonomous CPS mobility securely designed

abstract: In the last years the interconnection and ongoing development of physical systems combined with cyber resources has led to increasing automation. Through this progress in technology, autonomous vehicles, especially autonomous trains are getting more attention from industry and are already under test. The use of autonomous trains is known for increasing operation efficiency and reduction of personnel and infrastructure costs, which is mostly considered for main tracks. However, for less-used secondary lines, autonomous trains and their underlying sensor infrastructure are not yet considered. Thus, a system needs to be developed, which is less expensive for installation and operation of these trains and underlying infrastructure for secondary lines. Therefore, this position paper describes the process of how to derive an approach to help develop a digital interlocking system at design time for the use with secondary railway lines. In this work, we motivate the necessary research by investigating gaps in existing work as well as presenting a possible solution for this problem, a meta-model. The model considers safety, security as well as interoperability like 5G and socio-technical aspects to provide a holistic modeling approach for the development of the interlocking system for industrial secondary line use cases.


title: Energy-efficient Wireless Analog Sensing for Persistent Underwater Environmental Monitoring

abstract: The design of sensors or "things" as part of the new Internet of Underwater Things (IoUTs) paradigm comes with multiple challenges including limited battery capacity, not polluting the water body, and the ability to track continuously phenomena with high temporal/spatial variability. We claim that traditional digital sensors are incapable to meet these demands because of their high power consumption, high complexity (cost), and the use of non-biodegradable materials. To address the above challenges, we propose a novel architecture consisting of a sensing substrate of dense analog biodegradable sensors over which lies the traditional Wireless Sensor Network (WSN). The substrate analog biodegradable sensors perform Shannon mapping (a data-compression technique) using just a single Field Effect Transistor (FET) without the need for power-hungry Analog-to-Digital Converters (ADCs) resulting in much lower power consumption, complexity, and the ability to be powered using only sustainable energy-harvesting techniques. A novel and efficient decoding technique is also presented. Both encoding/decoding techniques have been verified via Spice and MATLAB simulations accounting for underwater acoustic channel variations.


title: Phonon Confinement and Transport in Ultrathin Films

abstract: Thermal transport by phonons in films with thicknesses of less than 10 nm is investigated in a soft system (Lennard-Jones argon) and a stiff system (Tersoff silicon) using two-dimensional lattice dynamics calculations and the Boltzmann transport equation. This approach uses a unit cell that spans the film thickness, which removes approximations related to the finite cross-plane dimension required in typical three-dimensional-based approaches. Molecular dynamics simulations are performed to obtain finite-temperature structures for the lattice dynamics calculations and to predict thermal conductivity benchmarks. Thermal conductivity decreases with decreasing film thickness until the thickness reaches four unit cells (2.1 nm) for argon and three unit cells (1.6 nm) for silicon. With a further decrease in film thickness, thermal conductivity plateaus in argon while it increases in silicon. This unexpected behavior, which we identify as a signature of phonon confinement, is a result of an increased contribution from low-frequency phonons, whose density of states increases as the film thickness decreases. Phonon mode-level analysis suggests that confinement effects emerge below thicknesses of ten unit cells (5.3 nm) for argon and six unit cells (3.2 nm) for silicon. Thermal conductivity predictions based on the bulk phonon properties combined with a boundary scattering model do not capture the low thickness behavior. To match the two-dimensional lattice dynamics and molecular dynamics predictions for larger thicknesses, the three-dimensional lattice dynamics calculations require a finite specularity parameter that in some cases approaches unity. These findings point to the challenges associated with interpreting experimental thermal conductivity measurements of ultrathin silicon films, where surface roughness and a native oxide layer impact phonon transport.


title: Probing superfluid $^4\mathrm{He}$ with high-frequency nanomechanical resonators down to $\mathrm{mK}$ temperatures

abstract: Superfluids, such as superfluid $^3\mathrm{He}$ and $^4\mathrm{He}$, exhibit a broad range of quantum phenomena and excitations which are unique to these systems. Nanoscale mechanical resonators are sensitive and versatile force detectors with the ability to operate over many orders of magnitude in damping. Using nanomechanical-doubly clamped beams of extremely high quality factors ($Q>10^6$), we probe superfluid $^4\mathrm{He}$ from the superfluid transition temperature down to $\mathrm{mK}$ temperatures at frequencies up to $11.6 \, \mathrm{MHz}$. Our studies show that nanobeam damping is dominated by hydrodynamic viscosity of the normal component of $^4\mathrm{He}$ above $1\,\mathrm{K}$. In the temperature range $0.3-0.8\,\mathrm{K}$, the ballistic quasiparticles (phonons and rotons) determine the beams' behavior. At lower temperatures, damping saturates and is determined either by magnetomotive losses or acoustic emission into helium. It is remarkable that all these distinct regimes can be extracted with just a single device, despite damping changing over six orders of magnitude.


title: Universal audio synthesizer control with normalizing flows

abstract: The ubiquity of sound synthesizers has reshaped music production and even entirely defined new music genres. However, the increasing complexity and number of parameters in modern synthesizers make them harder to master. Hence, the development of methods allowing to easily create and explore with synthesizers is a crucial need. Here, we introduce a novel formulation of audio synthesizer control. We formalize it as finding an organized latent audio space that represents the capabilities of a synthesizer, while constructing an invertible mapping to the space of its parameters. By using this formulation, we show that we can address simultaneously automatic parameter inference, macro-control learning and audio-based preset exploration within a single model. To solve this new formulation, we rely on Variational Auto-Encoders (VAE) and Normalizing Flows (NF) to organize and map the respective auditory and parameter spaces. We introduce the disentangling flows, which allow to perform the invertible mapping between separate latent spaces, while steering the organization of some latent dimensions to match target variation factors by splitting the objective as partial density evaluation. We evaluate our proposal against a large set of baseline models and show its superiority in both parameter inference and audio reconstruction. We also show that the model disentangles the major factors of audio variations as latent dimensions, that can be directly used as macro-parameters. We also show that our model is able to learn semantic controls of a synthesizer by smoothly mapping to its parameters. Finally, we discuss the use of our model in creative applications and its real-time implementation in Ableton Live


title: Efficient computation of collisional $\ell$-mixing rate coefficients in astrophysical plasmas

abstract: We present analytical expressions for direct evaluation of $\ell$-mixing rate coefficients in proton-excited hydrogen atom collisions and describe a software package for efficient numerical evaluation of the collisional rate coefficients. Comparisons between rate coefficients calculated with various levels of approximation are discussed, highlighting their range of validity. These rate coefficients are benchmarked via radio recombination lines for hydrogen, evaluating the corresponding departure coefficients from local thermal equilibrium.


title: Neural parameters estimation for brain tumor growth modeling

abstract: Understanding the dynamics of brain tumor progression is essential for optimal treatment planning. Cast in a mathematical formulation, it is typically viewed as evaluation of a system of partial differential equations, wherein the physiological processes that govern the growth of the tumor are considered. To personalize the model, i.e. find a relevant set of parameters, with respect to the tumor dynamics of a particular patient, the model is informed from empirical data, e.g., medical images obtained from diagnostic modalities, such as magnetic-resonance imaging. Existing model-observation coupling schemes require a large number of forward integrations of the biophysical model and rely on simplifying assumption on the functional form, linking the output of the model with the image information. In this work, we propose a learning-based technique for the estimation of tumor growth model parameters from medical scans. The technique allows for explicit evaluation of the posterior distribution of the parameters by sequentially training a mixture-density network, relaxing the constraint on the functional form and reducing the number of samples necessary to propagate through the forward model for the estimation. We test the method on synthetic and real scans of rats injected with brain tumors to calibrate the model and to predict tumor progression.


title: Spin dynamics of moving bodies in rotating black hole spacetimes

abstract: The dynamics of spinning test bodies, moving in rotating black hole (Kerr, Bardeen-like and Hayward-like) spacetimes, are investigated. In Kerr spacetime, all the spherical, zoom-whirl and unbound orbits are considered numerically. Along spherical orbits and for high spin, an amplitude modulation is found in the harmonic evolution of the spin precessional angular velocity, caused by the spin-curvature coupling. Along the discussed zoom-whirl and unbound orbits, the test body approaches the center so much that it passes through the ergosphere. Near and inside the ergosphere, the variation of the spin direction can be very rapid. The effects of the spin-curvature coupling is also investigated. The initial values are chosen such a way, that the body and its spin move in the equatorial plane of the coordinate space and of the comoving frame, respectively. Hence, a clear effect of the spin-curvature coupling is observed as the orbit and the spin vector leave the equatorial plane. Additional effects in the spin precessional angular velocity and in the evolution of the Boyer-Lindquist coordinate components of the spin vector is also considered. Finally, in case of different regular black holes, the spin-curvature coupling influences differently the orbit and the spin evolutions.


title: Stationary quantum correlations in a system with mean-field PT symmetry

abstract: A pair of coupled quantum harmonic oscillators, one subject to a gain one to a loss, is a paradigmatic setup to implement PT-symmetric, non-Hermitian Hamiltonians in that one such Hamiltonian governs the mean-field dynamics for equal gain and loss strengths. Through a full quantum description (so as to account for quantum noise) here is shown that when the system starts in any two-mode coherent state, including vacuum, there appear quantum correlations (QCs) without entanglement, as measured by the Gaussian discord. When the loss rate is above a threshold, once generated QCs no more decay. This occurs in a wide region of parameters, significantly larger than that where the full quantum dynamics is stable. For equal gain and loss rates, in particular, QCs decay in the exact phase (including the exceptional point) and are stable in the broken phase.

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


title: GASP. XX. From the loose spatially-resolved to the tight global SFR-Mass relation in local spiral galaxies

abstract: Exploiting the sample of 30 local star-forming, undisturbed late-type galaxies in different environments drawn from the GAs Stripping Phenomena in galaxies with MUSE (GASP), we investigate the spatially resolved Star Formation Rate-Mass ({\Sigma}SFR-{\Sigma}_star) relation. Our analysis includes also the galaxy outskirts (up to >4 effective radii, re), a regime poorly explored by other Integral Field Spectrograph surveys. Our observational strategy allows us to detect H{\alpha} out to more than 2.7re for 75% of the sample. Considering all galaxies together, the correlation between the {\Sigma}SFR and {\Sigma}_star is quite broad, with a scatter of 0.3 dex. It gets steeper and shifts to higher {\Sigma}_star values when external spaxels are excluded and moving from less to more massive galaxies. The broadness of the overall relation suggests galaxy-by-galaxy variations. Indeed, each object is characterized by a distinct {\Sigma}SFR-{\Sigma}_star relation and in some cases the correlation is very loose. The scatter of the relation mainly arises from the existence of bright off-center star-forming knots whose {\Sigma}SFR-{\Sigma}_star relation is systematically broader than that of the diffuse component. The {\Sigma}SFR-{\Sigma}tot gas (total gas surface density) relation is as broad as the {\Sigma}SFR-{\Sigma}_star relation, indicating that the surface gas density is not a primary driver of the relation. Even though a large galaxy-by-galaxy variation exists, mean {\Sigma}SFR and {\Sigma}_star values vary of at most 0.7 dex across galaxies. We investigate the relationship between the local and global SFR-M_star relation, finding that the latter is driven by the existence of the size-mass relation.


title: Rapid early coeval star formation and assembly of the most massive galaxies in the universe

abstract: The current consensus on the formation and evolution of the brightest cluster galaxies is that their stellar mass forms early ($z \gtrsim 4$) in separate galaxies that then eventually assemble the main structure at late times ($z \lesssim 1$). However, advances in observational techniques have led to the discovery of protoclusters out to $z \sim 7$, suggesting that the late-assembly picture may not be fully complete. If these protoclusters assemble rapidly in the early universe, they should form the brightest cluster galaxies much earlier than suspected by the late-assembly picture. Using a combination of observationally constrained hydrodynamical and dark-matter-only simulations, we show that the stellar assembly time of a sub-set of the brightest cluster galaxies occurs at high redshifts ($z > 3$) rather than at low redshifts ($z < 1$), as is commonly thought. We find, using isolated non-cosmological hydrodynamical simulations, that highly overdense protoclusters assemble their stellar mass into brightest cluster galaxies within $\sim 1$ $\mathrm{Gyr}$ of evolution -- producing massive blue elliptical galaxies at high redshifts ($z \gtrsim 1.5$). We argue that there is a downsizing effect on the cluster scale wherein some of the brightest cluster galaxies in the cores of the most-massive clusters assemble earlier than those in lower-mass clusters. In those clusters with $z = 0$ virial mass $\geqslant 5\times10^{14}$ M$_\mathrm{\odot}$, we find that $9.8$% have their cores assembly early, and a higher fraction of $16.4$% in those clusters above $10^{15}$ M$_\mathrm{\odot}$. The James Webb Space Telescope will be able to detect and confirm our prediction in the near future, and we discuss the implications to constraining the value of $\sigma_\mathrm{8}$.


title: Galactic cirri in deep optical imaging

abstract: The ubiquitous presence of Galactic cirri in deep optical images represents a major obstacle to study the low surface brightness features of extragalactic sources. To address this issue, we have explored the optical properties of cirri using g, r, i and z bands in the Sloan Digital Sky Survey (SDSS) Stripe82 region. Using state-of-the-art, custom made, image processing techniques, including the modeling and removal of the instrumental scattered light produced by the stars, we managed to isolate the optical diffuse emission by the cirri, allowing their photometric characterization. We find that their optical colors are driven by the dust column density: The cirri become redder as their 100 $\mu$m emission increases. Remarkably, the optical colors of the Galactic cirri differ significantly from those of extragalactic sources, with a characteristic bluer r-i color for a given g-r, allowing one to detect these by using a simple color relation. Our results show the high potential of deep multi-band optical photometry, on its own, identifying the presence of cirri at a higher spatial resolution than those provided by far-infrared observations. The combination of very deep data and multi-band photometry (as the one produced by LSST and Euclid) would make it possible to build dust maps of unprecedented quality.


title: Constraints on dark matter annihilation in dwarf spheroidal galaxies from low frequency radio observations

abstract: We present the first observational limits on the predicted synchrotron signals from particle Dark Matter annihilation models in dwarf spheroidal galaxies at radio frequencies below 1 GHz. We use a combination of survey data from the Murchison Widefield Array (MWA) and the Giant Metre-wave Radio Telescope (GMRT) to search for diffuse radio emission from 14 dwarf spheroidal galaxies. For in-situ magnetic fields of 1 $\mu G$ and any plausible value for the diffusion coefficient, our limits do not constrain any Dark Matter models. However, for stronger magnetic fields our data might provide constraints comparable to existing limits from gamma-ray and cosmic ray observations. Predictions for the sensitivity of the upgraded MWA show that models with Dark Matter particle mass up to $\sim$ 1.6 TeV (1 TeV) may be constrained for magnetic field of 2 $\mu G$ (1 $\mu G$). While much deeper limits from the future low frequency Square Kilometre Array (SKA) will challenge the LHC in searches for Dark Matter particles, the MWA provides a valuable first step toward the SKA at low frequencies.


title: Quasi-Dirac neutrino oscillations at DUNE and JUNO

abstract: Quasi-Dirac neutrinos are obtained when the Lagrangian density of a neutrino mass model contains both Dirac and Majorana mass terms, and the Majorana terms are sufficiently small. This type of neutrinos introduces new mixing angles and mass splittings into the Hamiltonian, which will modify the standard neutrino oscillation probabilities. In this paper, we focus on the case where the new mass splittings are too small to be measured, but new angles and phases are present. We perform a sensitivity study for this scenario for the upcoming experiments DUNE and JUNO, finding that they will improve current bounds on the relevant parameters. Finally, we also explore the discovery potential of both experiments, assuming that neutrinos are indeed quasi-Dirac particles.


title: Evolution of black hole and galaxy growth in a semi-numerical galaxy formation model

abstract: We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts.


title: Positive and negative feedback of AGN outflows in NGC 5728

abstract: We present a spatially-resolved analysis of ionized and molecular gas in a nearby Seyfert 2 galaxy NGC 5728, using the VLT/MUSE and ALMA data. We find ionized gas outflows out to ~kpc scales, which encounter the star formation ring at 1 kpc radius. The star formation rate of the encountering region is significantly high (~1.8 M$_{\rm sol}/yr/kpc^2$) compared to other regions in the ring. In contrast, the CO (2-1) emission is significantly weaker by a factor of ~3.5, indicating very high star formation efficiency. These results support the positive feedback scenario that the AGN-driven outflows compress the ISM in the ring, enhancing the star formation activity. In addition, we detect outflow regions outside of spiral arms, in which gas is likely to be removed from the spiral arms and no clear sign of star formation is detected. The overall impact of AGN outflows on the global star formation in NGC 5728 is limited, suggesting the feedback of the low-luminosity AGN is insignificant.


title: Quasi-Conformal Models and the Early Universe

abstract: Extensions of the Standard Model and general relativity featuring a UV fixed point can leave observable implications at accessible energies. Although mass parameters such as the Planck scale can appear through dimensional transmutation, all fundamental dimension-4 operators can (at least approximately) respect Weyl invariance at finite energy. An example is the Weyl-squared term, whose consistency and observational consequences are studied. This quasi-conformal scenario emerges from the UV complete quadratic gravity and is a possible framework for inflation. We find two realizations. In the first one the inflaton is a fundamental scalar with a quasi-conformal non-minimal coupling to the Ricci scalar. In this case the field excursion must not exceed the Planck mass by far. An example discussed in detail is hilltop inflation. In the second realization the inflaton is a pseudo-Goldstone boson (natural inflation). In this case we show how to obtain an elegant UV completion within an asymptotically free QCD-like theory, in which the inflaton is a composite scalar due to new strong dynamics. We also show how efficient reheating can occur. Unlike the natural inflation based on Einstein gravity, the tensor-to-scalar ratio is well below the current bound set by Planck. In both realizations mentioned above, the basic inflationary formulae are computed analytically and, therefore, these possibilities can be used as simple benchmark models.


title: Inflaxion Dark Matter

abstract: A new mechanism for producing axion dark matter is proposed. By invoking low-scale inflation and a kinetic mixing between the axion and the inflaton, it is shown that the axion is driven to a field point slightly displaced from the potential minimum, which can give rise to the observed dark matter abundance. In this framework, different combinations of the axion and inflaton fields play various cosmological roles, including generating the cosmological perturbations, reheating the universe, and serving as dark matter. The kinetic mixing also relates the dark matter lifetime with the reheating temperature. The mechanism tames axions that would otherwise overdominate the universe, and thus opens up new windows in the axion parameter space, including decay constants at the GUT scale and higher.


title: Photometric Supernovae Redshift Systematics Requirements

abstract: Imaging surveys will find many tens to hundreds of thousands of Type Ia supernovae in the next decade, and measure their light curves. In addition to a need for characterizing their types and subtypes, a redshift is required to place them on a Hubble diagram to map the cosmological expansion. We investigate the requirements on redshift systematics control in order not to bias cosmological results, in particular dark energy parameter estimation. We find that additive and multiplicative systematics must be constrained at the few$\,\times 10^{-3}$ level, effectively requiring spectroscopic followup for robust use of photometric supernovae. Catastrophic outliers need control at the subpercent level. We also investigate sculpting the spectroscopic sample.


title: Estimating black hole masses: Accretion disk fitting versus reverberation mapping and single epoch

abstract: We selected a sample of 28 Type 1 AGNs for which a black hole mass has been inferred using the reverberation mapping technique and single epoch scaling relations. All 28 sources show clear evidence of the "Big Blue Bump" in the optical-UV band whose emission is produced by an accretion disk (AD) around a supermassive black hole. We fitted the spectrum of these sources with the relativistic thin AD model KERRBB in order to infer the black hole masses and compared them with those from Reverberation mapping and Single epoch methods, discussing the possible uncertainties linked to such a model by quantifying their weight on our results. We find that for the majority of the sources, KERRBB is a good description of the AD emission for a wide wavelength range. The overall uncertainty on the black hole mass estimated through the disk fitting procedure is $\sim 0.45$ dex (which includes the uncertainty on fitting parameters such as e.g. spin and viewing angle), comparable to the systematic uncertainty of reverberation mapping and single epoch methods; however, such an uncertainty can be $\leq 0.3$ dex if one of the parameters of the fit is well constrained. Although all of the estimates are affected by large uncertainties, the masses inferred using the three methods are compatible if the dimensionless scale factor $f$ (linked to the unknown kinematics and geometry of the Broad Line Region) is assumed to be larger than one. For the majority of the sources, the comparison between the results coming from the three methods favors small spin values. To check the goodness of the KERRBB results, we compared them with those inferred with other models, such as AGNSED, a model that also accounts for the emission originating from an X-ray corona: using two sources with a good data coverage in the X band, we find that the masses estimated with the two models differ at most by a factor of $\sim 0.2$ dex.


title: Applying Liouville's Theorem to Gaia Data

abstract: The Milky Way is filled with the tidally-disrupted remnants of globular clusters and dwarf galaxies. Determining the properties of these objects -- in particular, initial masses and density profiles -- is relevant to both astronomy and dark matter physics. However, most direct measures of mass cannot be applied to tidal debris, as the systems of interest are no longer in equilibrium. Since phase-space density is conserved during adiabatic phase mixing, Liouville's theorem provides a connection between stellar kinematics as measured by observatories such as Gaia and the original mass of the disrupted system. Accurately recovering the phase-space density is complicated by uncertainties resulting from measurement errors and orbital integration, which both effectively inject entropy into the system, preferentially decreasing the measured density. In this paper, we demonstrate that these two issues can be overcome. First, we measure the phase-space density of the globular cluster M4 in Gaia data, and use Liouville's theorem to derive its mass. We then show that, for tidally disrupted systems, the orbital parameters and thus phase-space density can be inferred by minimizing the phase-space entropy of cold stellar streams. This work is therefore a proof of principle that true phase-space density can be measured and the original properties of the star cluster reconstructed in systems of astrophysical interest.


title: PDFSense: Mapping the PDF sensitivity of future facilities (HL-LHC, LHeC, and EIC)

abstract: Particle and nuclear physics are moving toward a new generation of experiments to stress-test the Standard Model (SM), search for novel degrees of freedom, and comprehensively map the internal structure of hadrons. Due to the complex nature of QCD and wide array of past, present, and possible future experiments, measurements taken at these next-generation facilities will inhabit an expansive space of high-energy data. Maximizing the impact of each future collider program will depend on identifying its place within this sprawling landscape. As an initial exploration, we use the recently-developed PDFSense framework to assess the PDF sensitivity of two future high-energy facilities --- the high-luminosity upgrade to the LHC (HL-LHC) and the Large Hadron-electron Collider (LHeC) proposal --- as well as the electron-ion collider (EIC) proposed to map the few-GeV quark-hadron transition region. We report that each of these experimental facilities occupies a unique place in the kinematical parameter space with specialized pulls on particular collinear quantities. As such, there is a clear complementarity among these programs, with an opportunity for each to mutually reinforce and inform the others.


title: Beyond two-point statistics: using the minimum spanning tree as a tool for cosmology

abstract: Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, $N$-body particles or haloes are used as nodes to construct a unique graph, the MST, that traces skeletal structure. We study the dependence of the MST on cosmological parameters using haloes from a suite of COLA simulations with a box size of $250\ h^{-1}{\rm Mpc}$, varying the amplitude of scalar fluctuations $\left(A_{\rm s}\right)$, matter density $\left(\Omega_{\rm m}\right)$, and neutrino mass $\left(\sum m_{\nu}\right)$. The power spectrum $P$ and bispectrum $B$ are measured for wavenumbers between $0.125$ and $0.5$ $h{\rm Mpc}^{-1}$, while a corresponding lower cut of $\sim12.6$ $h^{-1}{\rm Mpc}$ is applied to the MST. The constraints from the individual methods are fairly similar but when combined we see improved $1\sigma$ constraints of $\sim 17\%$ ($\sim 12\%$) on $\Omega_{\rm m}$ and $\sim 12\%$ ($\sim 10\%$) on $A_{\rm s}$ with respect to $P$ ($P+B$) thus showing the MST is providing additional information. The MST can be applied to current and future spectroscopic surveys (BOSS, DESI, Euclid, PSF, WFIRST, and 4MOST) in 3D and photometric surveys (DES and LSST) in tomographic shells to constrain parameters and/or test systematics.


title: Properties of Dual Supertranslation Charges in Asymptotically Flat Spacetimes

abstract: We study several properties of some new charges of asymptotically flat spacetimes. These dual supertranslation charges are akin to the magnetic large $U(1)$ charges in QED. In this paper we find the symmetries associated with these charges and show that the global dual supertranslation charge is topological because it is invariant under globally defined, smooth variations of the asymptotic metric. We also exhibit spacetimes where the charge does not vanish and we find dynamical processes that interpolate between regions with different values of these charges.


title: Neutrino Non-Standard Interactions: A Status Report

abstract: This report summarizes the present status of neutrino non-standard interactions (NSI). After a brief overview, several aspects of NSIs are discussed, including connection to neutrino mass models, model-building and phenomenology of large NSI with both light and heavy mediators, NSI phenomenology in both short- and long-baseline neutrino oscillation experiments, neutrino cross-sections, complementarity of NSI with other low- and high-energy experiments, fits with neutrino oscillation and scattering data, DUNE sensitivity to NSI, effective field theory of NSI, as well as the relevance of NSI to dark matter and cosmology. We also discuss the open questions and interesting future directions that can be pursued by the community at large. This report is based on talks and discussions during the Neutrino Theory Network NSI workshop held at Washington University in St. Louis from May 29-31, 2019 (https://indico.cern.ch/event/812851/)

categories: hep-ph astro-ph.CO astro-ph.HE hep-ex nucl-th


title: String corrections to AdS amplitudes and the double-trace spectrum of N=4 SYM

abstract: We consider $\alpha'$ corrections to four-point correlators of half-BPS operators in $\mathcal{N}=4$ super Yang-Mills theory in the supergravity limit. By demanding the correct behaviour in the flat space limit, we find that the leading $(\alpha')^3$ correction to the Mellin amplitude is fixed for arbitrary charges of the external operators. By considering the mixing of double-trace operators we can find the $(\alpha')^3$ corrections to the double-trace spectrum which we give explicitly for $su(4)$-singlet operators. We observe striking patterns in the corrections to the spectra which hint at their common ten-dimensional origin. By extending the observed patterns and imposing them at order $(\alpha')^5$ we are able to reproduce the recently found result for the correction to the Mellin amplitude for $\langle \mathcal{O}_2 \mathcal{O}_2 \mathcal{O}_p \mathcal{O}_p \rangle$ correlators. By applying a similar logic to the $[0,1,0]$ channel of $su(4)$ we are able to deduce new results for the correlators of the form $\langle \mathcal{O}_2 \mathcal{O}_3 \mathcal{O}_{p-1} \mathcal{O}_p \rangle$.


title: Dynamical Emergence of Scalaron in Higgs Inflation

abstract: We point out that a light scalaron dynamically emerges if scalar fields have a sizable non-minimal coupling to the Ricci scalar as in the Higgs inflation model. We support this claim in two ways. One is based on the renormalization group equation; the non-minimal coupling inevitably induces a Ricci scalar quadratic term due to the renormalization group running. The other is based on scattering amplitudes; a scalar four-point amplitude develops a pole after summing over a certain class of diagrams, which we identify as the scalaron. Our result implies that the Higgs inflation is actually a two-field inflationary model. Another implication is that the Higgs inflation does not suffer from the unitarity issue since the scalaron pushes up the cut-off scale to the Planck scale.


title: Compact object binary mergers driven by cluster tides: a new channel for LIGO/Virgo gravitational wave events

abstract: The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive black hole) until gravitational wave emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH-BH, NS-BH and NS-NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several pc. For NS-NS binaries we find merger rates in the range $0.01-0.07$ Gpc$^{-3}$ yr$^{-1}$ from globular clusters and $0.1-0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters. For NS-BH and BH-BH binaries we find small merger rates from globular clusters, but a rate of $0.1 - 0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several per cent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future gravitational wave detectors.


title: Uniformity and Stability of the LSST Focal Plane

abstract: The LSST focal plane consists of 21 autonomous modules ("Raft Tower Modules", RTM), each of which contains nine thick, fully-depleted 4K x 4K CCDs with associated control and readout electronics. To enable LSST's repetitive short-exposure cadence while maintaining high duty factor and low read noise, the readout is highly parallelized into 3024 independent video channels (16 per CCD, 144 per RTM). Two vendors supplied the LSST sensors; the devices have compatible mechanical and electrical interfaces and meet the same electro-optic specifications, but each RTM is constructed with sensors from a single supplier. The full complement of rafts were assembled at Brookhaven National Laboratory during January 2017 - March 2019. Each unit underwent extensive electro-optic and metrology characterization at operating temperature, the results of which are presented here along with a discussion of uniformity and stability.


title: A Bayesian Approach to Matching Thermonuclear X-ray Burst Observations with Models

abstract: We present a new method of matching observations of Type I (thermonuclear) X-ray bursts with models, comparing the predictions of a semi-analytic ignition model with X-ray observations of the accretion-powered millisecond pulsar SAX J1808.4--3658 in outburst. We used a Bayesian analysis approach to marginalise over the parameters of interest and determine parameters such as fuel composition, distance/anisotropy factors, neutron star mass and neutron star radius. Our study includes a treatment of the system inclination effects, inferring that the rotation axis of the system is inclined $\left(69^{+4}_{-2}\right)^\circ$ from the observers line of sight, assuming a flat disc model. This method can be applied to any accreting source that exhibits Type I X-ray bursts. We find a hydrogen mass fraction of $0.57^{+0.13}_{-0.14}$ and CNO metallicity of $0.013^{+0.006}_{-0.004}$ for the accreted fuel is required by the model to match the observed burst energies, for a distance to the source of $3.3^{+0.3}_{-0.2}\,\mathrm{kpc}$. We infer a neutron star mass of $1.5^{+0.6}_{-0.3}\,\mathrm{M}_{\odot}$ and radius of $11.8^{+1.3}_{-0.9}\,\mathrm{km}$ for a surface gravity of $1.9^{+0.7}_{-0.4}\times10^{14}\,\mathrm{cm}\,\mathrm{s}^{-2}$ for SAX J1808.4--3658.


title: Top quark decay at next-to-leading order in the Standard Model Effective Field Theory

abstract: We consider top quark decay in the Standard Model Effective Field Theory (SMEFT). We present a calculation of the total decay width and the $W$-boson helicity fractions at next-to-leading order (NLO) in SMEFT. Our result includes the complete set of contributing four-fermion operators in addition to QCD dipole operators and bottom-mass suppressed effects. We show that operators that first appear at NLO in the SMEFT can be bounded by the current data as well as future data from both a high-luminosity LHC and a potential $e^+e^-$ collider, demonstrating the importance of going beyond leading order when studying the SMEFT. We discuss technical aspects of our calculation that we believe will be useful in future higher-order studies of the SMEFT, in particular the treatment of $\gamma_5$ in loop diagrams.


title: Geographical Security Questions for Fallback Authentication

abstract: Fallback authentication is the backup authentication method used when the primary authentication method (e.g., passwords, fingerprints, etc.) fails. Currently, widely-deployed fallback authentication methods (e.g., security questions, email resets, and SMS resets) suffer from documented security and usability flaws that threaten the security of accounts. These flaws motivate us to design and study Geographical Security Questions (GeoSQ), a system for fallback authentication. GeoSQ is an Android application that utilizes autobiographical location data for fallback authentication. We performed security and usability analyses of GeoSQ through an in-person two-session lab study (n=36,18 pairs). Our results indicate that GeoSQ exceeds the security of its counterparts, while its usability (specifically login time) has room for improvement.


title: Evidence against a general positive eddy feedback in atmospheric blocking

abstract: The eddy straining mechanism of Shutts (1983; S83) has long been considered a main process for explaining the maintenance of atmospheric blocking. As hypothesized in S83, incoming synoptic eddies experience a meridional straining effect when approaching a split jetstream, and as a result, enhanced PV fluxes reinforce the block. A two-layer QG model is adopted here as a minimal model to conduct mechanism-denial experiments. While transient eddies' forcing is clearly critical to the formation and maintenance of a block, using a large ensemble, the authors demonstrate that the straining of generic eddies does not maintain blocks, thus challenge the idea of eddy straining serving as a positive feedback for the blocks. These results indicate that specific configurations of the eddy field are required for the maintenance stage. The authors also remark on the main supporting evidence in S83: the second-order induced flow is sensitive to the location of the wavemaker.


title: Twisted Spin in Quantum Mechanics

abstract: In quantum mechanics, it is often thought that the spin of an object points in a fixed direction at any point in time. For example, after selecting the z-direction as the axis of quantization, a spin-1/2 object (such as an electron) may either point up or down. The spin can also be a linear combination of these two states, in which case, there is an axis in another direction in which the spin points in that direction. In this article, we focus on the spin-spin-1/2 case and point out that spin may not necessarily point in a fixed direction, a phenomenon that we call twisted spin. We argue that twisted spin occurs in nature, that at least some degree of twisting is generic, and propose an experiment to verify its existence.


title: Quantum Fishchain in $AdS_5$

abstract: In our previous paper we derived the holographic dual of the planar fishnet CFT in four dimensions. The dual model becomes classical in the strongly coupled regime of the CFT and takes the form of an integrable chain of particles in five dimensions. Here we study the theory at the quantum level. By applying the canonical quantization procedure with constraints, we show that the model describes a quantum chain of particles propagating in $AdS_5$. We prove the duality at the full quantum level in the ${\mathfrak u}(1)$ sector and reproduce exactly the spectrum for the cases when it is known analytically.


title: Detectability of Collective Neutrino Oscillation Signatures in the Supernova Explosion of a 8.8 $M_\odot$ star

abstract: In order to investigate the impact of collective neutrino oscillations (CNO) on the neutrino signal from a nearby supernova, we perform 3-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 \Msun progenitor star. We find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\bar{\nu}_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to $\sim 10$ kpc. On the other hand, for the normal mass hierarchy, the spectrum of $\nu_e$ becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to $\sim 10$ kpc. More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\bar{\nu}_e$ in HK becomes softer due to CNO, the spectrum of $\nu_e$ in DUNE becomes harder, and vice versa. This synergistic observations in $\bar{\nu}_e$ and $\nu_e$, by HK and DUNE respectively, will be an intriguing opportunity to test the occurrence of CNO.


title: Accurate, reliable and fast robustness evaluation

abstract: Throughout the past five years, the susceptibility of neural networks to minimal adversarial perturbations has moved from a peculiar phenomenon to a core issue in Deep Learning. Despite much attention, however, progress towards more robust models is significantly impaired by the difficulty of evaluating the robustness of neural network models. Today's methods are either fast but brittle (gradient-based attacks), or they are fairly reliable but slow (score- and decision-based attacks). We here develop a new set of gradient-based adversarial attacks which (a) are more reliable in the face of gradient-masking than other gradient-based attacks, (b) perform better and are more query efficient than current state-of-the-art gradient-based attacks, (c) can be flexibly adapted to a wide range of adversarial criteria and (d) require virtually no hyperparameter tuning. These findings are carefully validated across a diverse set of six different models and hold for L0, L1, L2 and Linf in both targeted as well as untargeted scenarios. Implementations will soon be available in all major toolboxes (Foolbox, CleverHans and ART). We hope that this class of attacks will make robustness evaluations easier and more reliable, thus contributing to more signal in the search for more robust machine learning models.


title: Symmetry Detection and Classification in Drawings of Graphs

abstract: Symmetry is a key feature observed in nature (from flowers and leaves, to butterflies and birds) and in human-made objects (from paintings and sculptures, to manufactured objects and architectural design). Rotational, translational, and especially reflectional symmetries, are also important in drawings of graphs. Detecting and classifying symmetries can be very useful in algorithms that aim to create symmetric graph drawings and in this paper we present a machine learning approach for these tasks. Specifically, we show that deep neural networks can be used to detect reflectional symmetries with 92% accuracy. We also build a multi-class classifier to distinguish between reflectional horizontal, reflectional vertical, rotational, and translational symmetries. Finally, we make available a collection of images of graph drawings with specific symmetric features that can be used in machine learning systems for training, testing and validation purposes. Our datasets, best trained ML models, source code are available online.


title: Algorithms and data structures for matrix-free finite element operators with MPI-parallel sparse multi-vectors

abstract: Traditional solution approaches for problems in quantum mechanics scale as $\mathcal O(M^3)$, where $M$ is the number of electrons. Various methods have been proposed to address this issue and obtain linear scaling $\mathcal O(M)$. One promising formulation is the direct minimization of energy. Such methods take advantage of physical localization of the solution, namely that the solution can be sought in terms of non-orthogonal orbitals with local support. In this work a numerically efficient implementation of sparse parallel vectors within the open-source finite element library deal.II is proposed. The main algorithmic ingredient is the matrix-free evaluation of the Hamiltonian operator by cell-wise quadrature. Based on an a-priori chosen support for each vector we develop algorithms and data structures to perform (i) matrix-free sparse matrix multivector products (SpMM), (ii) the projection of an operator onto a sparse sub-space (inner products), and (iii) post-multiplication of a sparse multivector with a square matrix. The node-level performance is analyzed using a roofline model. Our matrix-free implementation of finite element operators with sparse multivectors achieves the performance of 157 GFlop/s on Intel Cascade Lake architecture. Strong and weak scaling results are reported for a typical benchmark problem using quadratic and quartic finite element bases.

categories: cs.MS cs.DS cs.NA math.NA physics.comp-ph


title: Enumeration of Preferred Extensions in Almost Oriented Digraphs

abstract: In this paper, we present enumeration algorithms to list all preferred extensions of an argumentation framework. This task is equivalent to enumerating all maximal semikernels of a directed graph. For directed graphs on $n$ vertices, all preferred extensions can be enumerated in $O^*(3^{n/3})$ time and there are directed graphs with $\Omega(3^{n/3})$ preferred extensions. We give faster enumeration algorithms for directed graphs with at most $0.8004\cdot n$ vertices occurring in $2$-cycles. In particular, for oriented graphs (digraphs with no 2-cycles) one of our algorithms runs in time $O(1.2321^n)$, and we show that there are oriented graphs with $\Omega(3^{n/6}) > \Omega(1.2009^n)$ preferred extensions. A combination of three algorithms leads to the fastest enumeration times for various proportions of the number of vertices in $2$-cycles. The most innovative one is a new 2-stage sampling algorithm, combined with a new parameterized enumeration algorithm, analyzed with a combination of the recent monotone local search technique (STOC 2016) and an extension thereof (ICALP 2017).