title: Low-Cost GPS-Aided LiDAR State Estimation and Map Building

abstract: Using different sensors in an autonomous vehicle (AV) can provide multiple constraints to optimize AV location estimation. In this paper, we present a low-cost GPS-assisted LiDAR state estimation system for AVs. Firstly, we utilize LiDAR to obtain highly precise 3D geometry data. Next, we use an inertial measurement unit (IMU) to correct point cloud misalignment caused by incorrect place recognition. The estimated LiDAR odometry and IMU measurement are then jointly optimized. We use a lost-cost GPS instead of a real-time kinematic (RTK) module to refine the estimated LiDAR-inertial odometry. Our low-cost GPS and LiDAR complement each other, and can provide highly accurate vehicle location information. Moreover, a low-cost GPS is much cheaper than an RTK module, which reduces the overall AV sensor cost. Our experimental results demonstrate that our proposed GPS-aided LiDAR-inertial odometry system performs very accurately. The accuracy achieved when processing a dataset collected in an industrial zone is approximately 0.14 m.


title: Lithium in red giant stars: Constraining non-standard mixing with large surveys in the Gaia era

abstract: Lithium is extensively known to be a good tracer of non-standard mixing processes occurring in stellar interiors. We present the results of a new large Lithium survey in red giant stars and combine it with surveys from the literature to probe the impact of rotation-induced mixing and thermohaline double-diffusive instability along stellar evolution. We determined the surface Li abundance for a sample of 829 giant stars with accurate Gaia parallaxes for a large sub-sample (810 stars) complemented with accurate Hipparcos parallaxes (19 stars). The spectra of our sample of northern and southern giant stars were obtained in three ground-based observatories (OHP, ESO-La Silla, Mc Donald). We determined the atmospheric parameters (Teff, log(g), [Fe/H]), and the Li abundance. We used Gaia parallaxes and photometry to determine the luminosity of our objects and we estimated the mass and evolution status of each sample star with a maximum-likelihood technique using stellar evolution models computed with the STAREVOL code. We compared the observed Li behaviour with predictions from stellar models, including rotation and thermohaline mixing. The same approach was used for stars from selected Li surveys from the literature. Rotation-induced mixing accounts nicely for the lithium behaviour in stars warmer than about 4200K, independently of the mass domain. For stars with masses lower than 2Msun thermohaline mixing leads to further Li depletion below the Teff of the RGB bump (about 4000K), and on the early AGB, as observed. Depending on the definition we adopt, we find between 0.8 and 2.2% of Li-rich giants in our new sample. Gaia puts a new spin on the understanding of mixing processes in stars, and our study confirms the importance of rotation-induced processes and of thermohaline mixing. However asteroseismology is required to definitively pinpoint the actual evolution status of Li-rich giants.


title: Constrained minimizers of the von Neumann entropy and their characterization

abstract: We consider in this work the problem of minimizing the von Neumann entropy under the constraints that the density of particles, the current, and the kinetic energy of the system is fixed at each point of space. The unique minimizer is a self-adjoint positive trace class operator, and our objective is to characterize its form. We will show that this minimizer is solution to a self-consistent nonlinear eigenvalue problem. One of the main difficulties in the proof is to parametrize the feasible set in order to derive the Euler-Lagrange equation, and we will proceed by constructing an appropriate form of perturbations of the minimizer. The question of deriving quantum statistical equilibria is at the heart of the quantum hydrody-namical models introduced by Degond and Ringhofer in [5]. An original feature of the problem is the local nature of constraints, i.e. they depend on position, while more classical models consider the total number of particles, the total current and the total energy in the system to be fixed.


title: A Semi-Automated Approach for Information Extraction, Classification and Analysis of Unstructured Data

abstract: In this paper, we show how Quantitative Narrative Analysis and simple Natural Language Processing techniques apply to the extraction and categorization of data in a sample case study of the Diary of the former President of the Italian Republic (PoR), Giorgio Napolitano. The Diary contains a record of all his institutional meetings. This information, if properly handled, allows for an analysis of how the PoR used his so-called soft-powers to influence the Italian political system during his first mandate. In this paper, we propose a way to use simple, yet very effective, Natural Language Processing techniques - such as Regular Expressions and Named Entity Recognition - to extract information from the Diary. Then, we propose an innovative way to organize the extracted data relying on the methodological framework of Quantitative Narrative Analysis. Finally, we show how to analyze the structured data under different levels of detail using PC-ACE (Program for Computer-Assisted Coding of Events), a software developed specifically for this task and for data visualization.


title: Learning to Recommend from Sparse Data via Generative User Feedback

abstract: Traditional collaborative filtering (CF) based recommender systems tend to perform poorly when the user-item interactions/ratings are highly scarce. To address this, we propose a learning framework that improves collaborative filtering with a synthetic feedback loop (CF-SFL) to simulate the user feedback. The proposed framework consists of a "recommender" and a "virtual user". The "recommender" is formulated as a CF model, recommending items according to observed user preference. The "virtual user" estimates rewards from the recommended items and generates a \emph{feedback} in addition to the observed user preference. The "recommender" connected with the "virtual user" constructs a closed loop, that recommends users with items and imitates the \emph{unobserved} feedback of the users to the recommended items. The synthetic feedback is used to augment the observed user preference and improve recommendation results. Theoretically, such model design can be interpreted as inverse reinforcement learning, which can be learned effectively via rollout (simulation). Experimental results show that the proposed framework is able to enrich the learning of user preference and boost the performance of existing collaborative filtering methods on multiple datasets.


title: A Framework for Data-Driven Computational Mechanics Based on Nonlinear Optimization

abstract: Data-Driven Computational Mechanics is a novel computing paradigm that enables the transition from standard data-starved approaches to modern data-rich approaches. At this early stage of development, one can distinguish two mainstream directions. The first one relies on a discrete-continuous optimization problem and seeks to assign to each material point a point in the phase space that satisfies compatibility and equilibrium, while being closest to the data set provided. The second one is a data driven inverse approach that seeks to reconstruct a constitutive manifold from data sets by manifold learning techniques, relying on a well-defined functional structure of the underlying constitutive law. In this work, we propose a third route that combines the strengths of the two existing directions and mitigates some of their weaknesses. This is achieved by the formulation of an approximate nonlinear optimization problem, which can be robustly solved, is computationally efficient, and does not rely on any special functional structure of the reconstructed constitutive manifold. Additional benefits include the natural incorporation of kinematic constraints and the possibility to operate with implicitly defined stress-strain relations. We discuss important mathematical aspects of our approach for a data-driven truss element and investigate its key numerical behavior for a data-driven beam element that makes use of all components of our methodology.


title: Quantum effects due to a moving Dirichlet point

abstract: We study quantum effects induced by a point-like object that imposes Dirichlet boundary conditions along its world-line, on a real scalar field $\varphi$ in 1, 2 and 3 spatial dimensions. The boundary conditions result from the strong coupling limit of a term quadratic in the field and localized on the particle's trajectory. We discuss the renormalization issues that appear and evaluate the effective action. Special attention is paid to the case of 2 spatial dimensions where the coupling constant is adimensional.


title: Distinguishing numbers and distinguishing indices of oriented graphs

abstract: A distinguishing r-vertex-labelling (resp. r-edge-labelling) of an undirected graph G is a mapping $\lambda$ from the set of vertices (resp. the set of edges) of G to the set of labels {1,. .. , r} such that no non-trivial automorphism of G preserves all the vertex (resp. edge) labels. The distinguishing number D(G) and the distinguishing index D (G) of G are then the smallest r for which G admits a distinguishing r-vertex-labelling or r-edge-labelling, respectively. The distinguishing chromatic number D $\chi$ (G) and the distinguishing chromatic index D $\chi$ (G) are defined similarly, with the additional requirement that the corresponding labelling must be a proper colouring. These notions readily extend to oriented graphs, by considering arcs instead of edges. In this paper, we study the four corresponding parameters for oriented graphs whose underlying graph is a path, a cycle, a complete graph or a bipartite complete graph. In each case, we determine their minimum and maximum value, taken over all possible orientations of the corresponding underlying graph, except for the minimum values for unbalanced complete bipartite graphs K m,n with m = 2, 3 or 4 and n > 3, 6 or 13, respectively, or m $\ge$ 5 and n > 2 m -- m 2 , for which we only provide upper bounds.


title: Qubits as edge state detectors: illustration using the SSH model

abstract: As is well known, qubits are the fundamental building blocks of quantum computers, and more generally, of quantum information. A major challenge in the development of quantum devices arises because the information content in any quantum state is rather fragile, as no system is completely isolated from its environment. Generally, such interactions degrade the quantum state, resulting in a loss of information. Topological edge states are promising in this regard because they are in ways more robust against noise and decoherence. But creating and detecting edge states can be challenging. We describe a composite system consisting of a two-level system (the qubit) interacting with a finite Su-Schrieffer-Heeger chain (a hopping model with alternating hopping parameters) attached to an infinite chain. In this model, the dynamics of the qubit changes dramatically depending on whether or not an edge state exists. Thus, the qubit can be used to determine whether or not an edge state exists in this model.


title: Sequence-to-sequence Automatic Speech Recognition with Word Embedding Regularization and Fused Decoding

abstract: In this paper, we investigate the benefit that off-the-shelf word embedding can bring to the sequence-to-sequence (seq-to-seq) automatic speech recognition (ASR). We first introduced the word embedding regularization by maximizing the cosine similarity between a transformed decoder feature and the target word embedding. Based on the regularized decoder, we further proposed the fused decoding mechanism. This allows the decoder to consider the semantic consistency during decoding by absorbing the information carried by the transformed decoder feature, which is learned to be close to the target word embedding. Initial results on LibriSpeech demonstrated that pre-trained word embedding can significantly lower ASR recognition error with a negligible cost, and the choice of word embedding algorithms among Skip-gram, CBOW and BERT is important.


title: Machine-Learning-Assisted Metasurface Design for High-Efficiency Thermal Emitter Optimization

abstract: With the emergence of new photonic and plasmonic materials with optimized properties as well as advanced nanofabrication techniques, nanophotonic devices are now capable of providing solutions to global challenges in energy conversion, information technologies, chemical/biological sensing, space exploration, quantum computing, and secure communication. Addressing grand challenges poses inherently complex, multi-disciplinary problems with a manifold of stringent constraints in conjunction with the required system's performance. Conventional optimization techniques have long been utilized as powerful tools to address multi-constrained design tasks. One example is so-called topology optimization that has emerged as a highly successful architect for the advanced design of non-intuitive photonic structures. Despite many advantages, this technique requires substantial computational resources and thus has very limited applicability to highly constrained optimization problems within high-dimensions parametric space. In our approach, we merge the topology optimization method with machine learning algorithms such as adversarial autoencoders and show substantial improvement of the optimization process by providing unparalleled control of the compact design space representations. By enabling efficient, global optimization searches within complex landscapes, the proposed compact hyperparametric representations could become crucial for multi-constrained problems. The proposed approach could enable a much broader scope of the optimal designs and data-driven materials synthesis that goes beyond photonic and optoelectronic applications.


title: A Gaussian process related to the mass spectrum of the near-critical Ising model

abstract: Let $\Phi^h(x)$ with $x=(t,y)$ denote the near-critical scaling limit of the planar Ising magnetization field. We take the limit of $\Phi^h$ as the spatial coordinate $y$ scales to infinity with $t$ fixed and prove that it is a stationary Gaussian process $X(t)$ whose covariance function is the Laplace transform of a mass spectral measure $\rho$ of the relativistic quantum field theory associated to the Euclidean field $\Phi^h$. Our analysis of the small distance/time behavior of the covariance functions of $\Phi^h$ and $X(t)$ shows that $\rho$ is finite but has infinite first moment.


title: The analytic theory of vectorial Drinfeld modular forms

abstract: In this paper we generalise the notion of Drinfeld modular form for the group $\Gamma$ := GL2(Fq[$\theta$]) to a vector-valued setting, where the target spaces are certain modules over positive characteristic Banach algebras over which are defined what we call the 'representations of the first kind'. Under quite reasonable restrictions, we show that the spaces of such modular forms are finite-dimensional, are endowed with certain generalisations of Hecke operators, with differential operators{\`a} la Serre etc. The crucial point of this work is the introduction of a 'field of uniformisers', a valued field in which we can study the expansions at the cusp infinity of our modular forms and which is wildly ramified. Examples of such modular forms are given through the construction of Poincar{\'e} and Eisenstein series. After this the paper continues with a more detailed analysis of the special case of modular forms associated to a restricted class of representations $\rho$ * $\Sigma$ of $\Gamma$ which has more importance in arithmetical applications. More structure results are given in this case, and a harmonic product formula is obtained which allows, with the help of three conjectures on the structure of an Fp-algebra of A-periodic multiple sums, to produce conjectural formulas for Eisenstein series. Some of these formulas can be proved.


title: On approximating $\nabla f$ with neural networks

abstract: Consider a feedforward neural network $\psi: \mathbb{R}^d\rightarrow \mathbb{R}^d$ such that $\psi\approx \nabla f$, where $f:\mathbb{R}^d \rightarrow \mathbb{R}$ is a smooth function, therefore $\psi$ must satisfy $\partial_j \psi_i = \partial_i \psi_j$ pointwise. We prove a theorem that a $\psi$ network with more than one hidden layer can only represent one feature in its first hidden layer; this is a dramatic departure from the well-known results for one hidden layer. The proof of the theorem is straightforward, where two backward paths and a weight-tying matrix play the key roles. We then present the alternative, the implicit parametrization, where the neural network is $\phi: \mathbb{R}^d \rightarrow \mathbb{R}$ and $\nabla \phi \approx \nabla f$; in addition, a "soft analysis" of $\nabla \phi$ gives a dual perspective on the theorem. Throughout, we come back to recent probabilistic models that are formulated as $\nabla \phi \approx \nabla f$, and conclude with a critique of denoising autoencoders.


title: Machine learning based data retrieval for inverse scattering problems with incomplete data

abstract: We are concerned with the inverse scattering problems associated with incomplete measurement data. It is a challenging topic of increasing importance in many practical applications. Based on a prototypical working model, we propose a machine learning based inverse scattering scheme, which integrates a CNN (convolution neural network) for the data retrieval. The proposed method can effectively cope with the reconstruction under limited-aperture and/or phaseless far-field data. Numerical experiments verify the promising features of our new scheme.


title: A Lyapunov-based small-gain theorem for infinite networks

abstract: This paper presents a small-gain theorem for networks composed of a countably infinite number of finite-dimensional subsystems. Assuming that each subsystem is exponentially input-to-state stable, we show that if the gain operator, collecting all the information about the internal Lyapunov gains, has a spectral radius less than one, the overall infinite network is exponentially input-to-state stable. The effectiveness of our result is illustrated through several examples including nonlinear spatially invariant systems with sector nonlinearities and a road traffic network.


title: Introduction to local certification

abstract: A distributed graph algorithm is basically an algorithm where every node of a graph can look at its neighborhood at some distance in the graph and chose its output. As distributed environment are subject to faults, an important issue is to be able to check that the output is correct, or in general that the network is in proper configuration with respect to some predicate. One would like this checking to be very local, to avoid using too much resources. Unfortunately most predicates cannot be checked this way, and that is where certification comes into play. Local certification (also known as proof-labeling schemes, locally checkable proofs or distributed verification) consists in assigning labels to the nodes, that certify that the configuration is correct. There are several point of view on this topic: it can be seen as a part of self-stabilizing algorithms, as labeling problem, or as a non-deterministic distributed decision. This paper is an introduction to the domain of local certification, giving an overview of the history, the techniques and the current research directions.


title: A Study of Machine Learning Models in Predicting the Intention of Adolescents to Smoke Cigarettes

abstract: The use of electronic cigarette (e-cigarette) is increasing among adolescents. This is problematic since consuming nicotine at an early age can cause harmful effects in developing teenager's brain and health. Additionally, the use of e-cigarette has a possibility of leading to the use of cigarettes, which is more severe. There were many researches about e-cigarette and cigarette that mostly focused on finding and analyzing causes of smoking using conventional statistics. However, there is a lack of research on developing prediction models, which is more applicable to anti-smoking campaign, about e-cigarette and cigarette. In this paper, we research the prediction models that can be used to predict an individual e-cigarette user's (including non-e-cigarette users) intention to smoke cigarettes, so that one can be early informed about the risk of going down the path of smoking cigarettes. To construct the prediction models, five machine learning (ML) algorithms are exploited and tested for their accuracy in predicting the intention to smoke cigarettes among never smokers using data from the 2018 National Youth Tobacco Survey (NYTS). In our investigation, the Gradient Boosting Classifier, one of the prediction models, shows the highest accuracy out of all the other models. Also, with the best prediction model, we made a public website that enables users to input information to predict their intentions of smoking cigarettes.


title: HIDRA: Head Initialization across Dynamic targets for Robust Architectures

abstract: The performance of gradient-based optimization strategies depends heavily on the initial weights of the parametric model. Recent works show that there exist weight initializations from which optimization procedures can find the task-specific parameters faster than from uniformly random initializations and that such a weight initialization can be learned by optimizing a specific model architecture across similar tasks via MAML (Model-Agnostic Meta-Learning). Current methods are limited to populations of classification tasks that share the same number of classes due to the static model architectures used during meta-learning. In this paper, we present HIDRA, a meta-learning approach that enables training and evaluating across tasks with any number of target variables. We show that Model-Agnostic Meta-Learning trains a distribution for all the neurons in the output layer and a specific weight initialization for the ones in the hidden layers. HIDRA explores this by learning one master neuron, which is used to initialize any number of output neurons for a new task. Extensive experiments on the Miniimagenet and Omniglot data sets demonstrate that HIDRA improves over standard approaches while generalizing to tasks with any number of target variables. Moreover, our approach is shown to robustify low-capacity models in learning across complex tasks with a high number of classes for which regular MAML fails to learn any feasible initialization.


title: Deep-Learning-Based Image Segmentation Integrated with Optical Microscopy for Automatically Searching for Two-Dimensional Materials

abstract: Deep-learning algorithms enable precise image recognition based on high-dimensional hierarchical image features. Here, we report the development and implementation of a deep-learning-based image segmentation algorithm in an autonomous robotic system to search for two-dimensional (2D) materials. We trained the neural network based on Mask-RCNN on annotated optical microscope images of 2D materials (graphene, hBN, MoS2, and WTe2). The inference algorithm is run on a 1024 x 1024 px2 optical microscope images for 200 ms, enabling the real-time detection of 2D materials. The detection process is robust against changes in the microscopy conditions, such as illumination and color balance, which obviates the parameter-tuning process required for conventional rule-based detection algorithms. Integrating the algorithm with a motorized optical microscope enables the automated searching and cataloging of 2D materials. This development will allow researchers to utilize unlimited amounts of 2D materials simply by exfoliating and running the automated searching process.


title: Dissipative solutions to a system for the flow of magnetoviscoelastic materials

abstract: We address the question of global in time existence of solutions to a magnetoviscoelastic system with general initial data. We show that the notion of dissipative solutions allows to prove such an existence in two and three dimensions. This extends an earlier result for the viscoelastic subsystem to the setting which includes the magnetization vector and its evolution in terms of a Landau-Lifshitz-Gilbert equation.


title: Analytical classical density functionals from an equation learning network

abstract: We explore the feasibility of using machine learning methods to obtain an analytic form of the classical free energy functional for two model fluids, hard rods and Lennard--Jones, in one dimension . The Equation Learning Network proposed in Ref. 1 is suitably modified to construct free energy densities which are functions of a set of weighted densities and which are built from a small number of basis functions with flexible combination rules. This setup considerably enlarges the functional space used in the machine learning optimization as compared to previous work 2 where the functional is limited to a simple polynomial form. As a result, we find a good approximation for the exact hard rod functional and its direct correlation function. For the Lennard--Jones fluid, we let the network learn (i) the full excess free energy functional and (ii) the excess free energy functional related to interparticle attractions. Both functionals show a good agreement with simulated density profiles for thermodynamic parameters inside and outside the training region.


title: Patient-specific fine-tuning of CNNs for follow-up lesion quantification

abstract: Convolutional neural network (CNN) methods have been proposed to quantify lesions in medical imaging. Commonly more than one imaging examination is available for a patient, but the serial information in these images often remains unused. CNN-based methods have the potential to extract valuable information from previously acquired imaging to better quantify current imaging of the same patient. A pre-trained CNN can be updated with a patient's previously acquired imaging: patient-specific fine-tuning. In this work, we studied the improvement in performance of lesion quantification methods on MR images after fine-tuning compared to a base CNN. We applied the method to two different approaches: the detection of liver metastases and the segmentation of brain white matter hyperintensities (WMH). The patient-specific fine-tuned CNN has a better performance than the base CNN. For the liver metastases, the median true positive rate increases from 0.67 to 0.85. For the WMH segmentation, the mean Dice similarity coefficient increases from 0.82 to 0.87. In this study we showed that patient-specific fine-tuning has potential to improve the lesion quantification performance of general CNNs by exploiting the patient's previously acquired imaging.


title: Deep XMM-Newton observations of the northern disk of M31 II: Tracing the hot interstellar medium

abstract: Aims: We use new deep XMM-Newton observations of the northern disk of M 31 to trace the hot interstellar medium (ISM) in unprecedented detail and to characterise the physical properties of the X-ray emitting plasmas. Methods: We used all XMM-Newton data up to and including our new observations to produce the most detailed image yet of the hot ISM plasma in a grand design spiral galaxy such as our own. We compared the X-ray morphology to multi-wavelength studies in the literature to set it in the context of the multi-phase ISM. We performed spectral analyses on the extended emission using our new observations as they offer sufficient depth and count statistics to constrain the plasma properties. Data from the Panchromatic Hubble Andromeda Treasury were used to estimate the energy injected by massive stars and their supernovae. We compared these results to the hot gas properties. Results: The brightest emission regions were found to be correlated with populations of massive stars, notably in the 10 kpc star-forming ring. The plasma temperatures in the ring regions are ~0.2 keV up to ~0.6 keV. We suggest this emission is hot ISM heated in massive stellar clusters and superbubbles. We derived X-ray luminosities, densities, and pressures for the gas in each region. We also found large extended emission filling low density gaps in the dust morphology of the northern disk, notably between the 5 kpc and 10 kpc star-forming rings. We propose that the hot gas was heated and expelled into the gaps by the populations of massive stars in the rings. Conclusions. It is clear that the massive stellar populations are responsible for heating the ISM to X-ray emitting temperatures, filling their surroundings, and possibly driving the hot gas into the low density regions. Overall, the morphology and spectra of the hot gas in the northern disk of M 31 is similar to other galaxy disks.


title: Two recent p-adic approaches towards the (effective) Mordell conjecture

abstract: We give an introductory account of two recent approaches towards an effective proof of the Mordell conjecture, due to Lawrence--Venkatesh and Kim. The latter method, which is usually called the method of Chabauty--Kim or non-abelian Chabauty in the literature, has the advantage that in some cases it has been turned into an effective method to determine the set of rational points on a curve, and we illustrate this by presenting three new examples of modular curves where this set can be determined.


title: Fast classification rates without standard margin assumptions

abstract: We consider the classical problem of learning rates for classes with finite VC dimension. It is well known that fast learning rates up to $O\left(\frac{d}{n}\right)$ are achievable by the empirical risk minimization algorithm (ERM) if low noise or margin assumptions are satisfied. These usually require the optimal Bayes classifier to be in the class, and it has been shown that when this is not the case, the fast rates cannot be achieved even in the noise free case. In this paper, we further investigate the question of the fast rates under the misspecification, when the Bayes classifier is not in the class (also called the agnostic setting). First, we consider classification with a reject option, namely Chow's reject option model, and show that by slightly lowering the impact of hard instances, a learning rate of order $O\left(\frac{d}{n}\log \frac{n}{d}\right)$ is always achievable in the agnostic setting by a specific learning algorithm. Similar results were only known under special versions of margin assumptions. We also show that the performance of the proposed algorithm is never worse than the performance of ERM. Based on those results, we derive the necessary and sufficient conditions for classification (without a reject option) with fast rates in the agnostic setting achievable by improper learners. This simultaneously extends the work of Massart and N\'{e}d\'{e}lec (Ann. of Statistics, 2006), which studied this question in the case where the Bayesian optimal rule belongs to the class, and the work of Ben-David and Urner (COLT, 2014), which allows the misspecification but is limited to the no noise setting. Our result also provides the first general setup in statistical learning theory in which an improper learning algorithm may significantly improve the learning rate for non-convex losses.


title: A Large-Scale Deep Architecture for Personalized Grocery Basket Recommendations

abstract: With growing consumer adoption of online grocery shopping through platforms such as Amazon Fresh, Instacart, and Walmart Grocery, there is a pressing business need to provide relevant recommendations throughout the customer journey. In this paper, we introduce a production within-basket grocery recommendation system, RTT2Vec, which generates real-time personalized product recommendations to supplement the user's current grocery basket. We conduct extensive offline evaluation of our system and demonstrate a 9.4% uplift in prediction metrics over baseline state-of-the-art within-basket recommendation models. We also propose an approximate inference technique 11.6x times faster than exact inference approaches. In production, our system has resulted in an increase in average basket size, improved product discovery, and enabled faster user check-out


title: A randomly determined unpredictable function

abstract: In this paper, we construct a new unpredictable function. Our approach is based on adapting the concept of symbolic dynamics to introduce a map on the space of infinite sequences generated by the discrete distribution. We show that there exists an unpredictable sequence on the space and then use the sequence to construct an unpredictable function.


title: Is Hyperon Polarization in Relativistic Heavy Ion Collisions Connected to Axial U(1) Symmetry Breaking at High Temperature?

abstract: Experiments at the Relativistic Heavy Ion Collider (RHIC) have measured the net polarization of $\Lambda$ and $\bar{\Lambda}$ hyperons and attributed it to a coupling between their spin and the vorticity of the fluid created in heavy ion collisions, but how the spin comes to equilibrium with vorticity is an open problem. Recently we found that vorticity fluctuations and helicity flip of strange quarks in quark-gluon plasma through perturbative QCD processes resulted in equilibration times far too long to be relevant. Here we consider the Nambu--Jona-Lasinio model with the inclusion of the six-quark Kobayashi--Maskawa--'t Hooft interaction which breaks axial U(1). Using instanton inspired models for the temperature dependence of the axial symmetry breaking, we find that constituent strange quarks can reach spin equilibrium at temperatures below about 170 MeV, just before they hadronize to form hyperons.


title: Quantum enhancements for deep reinforcement learning in large spaces

abstract: In the past decade, the field of quantum machine learning has drawn significant attention due to the prospect of bringing genuine computational advantages to now widespread algorithmic methods. However, not all domains of machine learning have benefited equally from quantum enhancements. Notably, deep learning and reinforcement learning, despite their tremendous success in the classical domain, both individually and combined, remain relatively unaddressed by the quantum community. Arguably, one reason behind this is the systematic use in these domains of models and methods without prominent computational bottlenecks, leaving little room for quantum improvements. In this work, we study the state-of-the-art neural-network approaches for reinforcement learning with quantum enhancements in mind. We demonstrate the substantial learning advantage that models with a sampling bottleneck can provide over conventional neural network architectures in complex learning environments. These so-called energy-based models, like deep energy-based reinforcement learning, and deep projective simulation that we also introduce in this work, effectively allow to trade off learning performance for efficiency of computation. To alleviate the additional computational costs, we propose to leverage future and near-term quantum algorithms, resulting in overall more advantageous learning algorithms. This is achieved using cutting-edge and new quantum computing machinery to speed-up classical sampling methods and by employing generalized models to gain an additional quantum advantage.


title: On the compression of the fullerene shell by an extra positive charge at its center

abstract: In this Letter, we investigate the variation of endohedral A@CN potential due to addition at the center of it a positive charge, for example, in the process of atom A photoionization. Using a reasonable model to describe the fullerenes shell, we managed to calculate the variation that is a consequence of the monopole polarization of CN shell. We analyze model potentials with flat and non-flat bottoms and demonstrate that the phenomenological potentials that properly simulates the C60 shell potential should belong to a family of potentials with a non-flat bottom. As concrete example, we use the Lorentz-bubble model potential. By varying the thickness of this potential, we describe the various degrees of the monopole polarization of the C60 shell by positive electric charge in the center of the shell. We calculated the photoionization cross-sections of He, Ar and Xe atoms located at the center of C60 shell with and without taking into account accompanying this process monopole polarization of the fullerenes shell. Unexpectedly, we found that the monopole polarization do not affect the photoionization cross sections of these endohedral atoms.


title: Dark matter to baryon ratio from scalar triplets decay in type-II seesaw

abstract: We propose a minimal model for the cosmic coincidence problem $\Omega_{\rm DM}/\Omega_B \sim 5$ and neutrino mass in a type-II seesaw scenario. We extend the standard model of particle physics with a $\rm SU(2)$ singlet leptonic Dirac fermion $\chi$, which represents the candidate of dark matter (DM), and two triplet scalars $\Delta_{1,2}$ with hierarchical masses. In the early Universe, the CP violating out-of-equilibrium decay of lightest $\Delta$ generates a net $B-L$ asymmetry in the visible sector (comprising of SM fields), where $B$ and $L$ represents the total baryon and lepton number respectively. A part of this asymmetry gets transferred to the dark sector (comprising of DM $\chi$) through a dimension eight operator which conserves $B-L$. Above the electroweak phase transition, the $B-L$ asymmetry of the visible sector gets converted to a net $B$-asymmetry by the $B+L$ violating sphalerons, while the $B-L$ asymmetry of the dark sector remains untouched which we see today as relics of DM. We show that the observed DM abundance can be explained for a DM mass about 8 GeV. We then introduce an additional singlet scalar field $\phi$ which mixes with the SM-Higgs to annihilate the symmetric component of the DM resonantly which requires the singlet scalar mass to be twice the DM mass, {\it i.e.} around 16 GeV, which can be searched at collider experiments. In our model, the active neutrinos also get small masses by the induced vacuum expectation value (vev) of the triplet scalars $\Delta_{1,2}$. In the later part of the paper we discuss all the constraints on model parameters coming from invisible Higgs decay, Higgs signal strength, DM direct detection and relic density of DM.


title: On positionality of trigger strategies Nash Equilibria in SCAR

abstract: We study the positionality of \emph{trigger strategies} Nash equilibria $\overline{\sigma}$ for the $N$-player SCAR games $\Gamma_{N}(G|s_{0},\gamma,\varepsilon)$ (with $N\geq3$). Our study is exhaustive with respect to types of graphs $G$, initial states $s_{0}$ and values of $N,\gamma,\varepsilon$. We conclude that in the majority of cases, profiles $\overline{\sigma}$ are nonpositional. Whenever $\overline{\sigma}$ are positional a key role is played by paths and the $\varepsilon$, $\gamma$ values (especially whether $\varepsilon>0$ or not). A crucial concept in our analysis is the \emph{state cop number}, which is first introduced in the current paper.


title: Bounded integral and motivic Milnor fiber

abstract: We construct a new motivic integration morphism, the so-call bounded integral, that interpolates both the integration morphisms with and without volume forms of Hrushovski and Kazhdan. This is done within the framework of model theory of algebraically closed valued fields of equicharacteristic zero. As an application, we recover and extend some results of Hrushovski and Loeser about the motivic Milnor fiber.


title: Symmetry-breaking vortex-lattice of a binary superfluid in a rotating bucket

abstract: We study spontaneous-symmetry-broken phase-separated vortex lattice in a weakly interacting uniform rapidly rotating binary Bose superfluid contained in a quasi-two-dimensional circular or square bucket. For the inter-species repulsion above a critical value, the two superfluid components separate and form a demixed phase with practically no overlap in the vortex lattices of the two components, which will permit an efficient experimental observation of such vortices and study their properties. In case of a circular bucket with equal intra-species energies of the two components, the two components separate into two non-overlapping semicircular domains for all frequencies of rotation $\Omega$ generating distinct demixed vortex lattices. In case of a binary Bose superfluid in both circular and square buckets, (a) the number of vortices increases linearly with $\Omega$ in agreement with a suggestion by Feynman, and (b) the rotational energy in the rotating frame decreases quadratically with $\Omega$ in agreement with a suggestion by Fetter.


title: Implementation of annihilation and junction reactions in vector density-based continuum dislocation dynamics

abstract: In a continuum dislocation dynamics formulation by Xia and El-Azab, dislocations are represented by a set of vector density fields, one per crystallographic slip systems. The space-time evolution of these densities is obtained by solving a set of dislocation transport equations coupled with crystal mechanics. Here, we present an approach for incorporating dislocation annihilation and junction reactions into the dislocation transport equations. These reactions consume dislocations and result in nothing as in the annihilation reactions, or produce new dislocations of different types as in the case of junction reactions. Collinear annihilation, glissile junctions, and sessile junctions are particularly emphasized here. A generalized energy-based criterion for junction reactions is established in terms of the dislocation density and Burgers vectors of the reacting species, and the reaction rate terms for junction reactions are formulated in terms of the dislocation densities. In order to illustrate how the dislocation network changes as a result of junction formation and annihilation in a continuum dislocation dynamics setting, we present some numerical examples focusing on the reactions processes themselves. The results show that our modeling approach is able to capture the respective dislocation network changes associated with dislocation reactions in FCC crystals: dislocations of opposite line directions encountering each other on collinear slip systems annihilate to connect the dislocations on the two slip systems, glissile junctions form on new slip system behave like Frank-Read sources, and sessile junctions form and expand along the intersection of the slip planes of the reacting dislocation species. A collective-dynamics test showing the frequency of occurrence of junctions of different types relative to each other is also presented.


title: Latency and timeliness in multi-hop satellite networks

abstract: The classical definition of network delay has been recently augmented by the concept of information timeliness, or Age of Information (AoI). We analyze the network delay and the AoI in a multi-hop satellite network that relays status updates from satellite 1, receiving uplink traffic from ground devices, to satellite K, using K-2 intermediate satellite nodes. The last node, K, is the closest satellite with connectivity to a ground station. The satellite formation is modeled as a queue network of M/M/1 systems connected in series. The scenario is then generalized for the case in which all satellites receive uplink traffic from ground, and work at the same time as relays of the packets from the previous nodes. The results show that the minimum average AoI is experienced at a decreasing system utilization when the number of nodes is increased. Furthermore, unloading the first nodes of the chain reduces the queueing time and therefore the average AoI. These findings provide insights for designing multi-hop satellite networks for latency-sensitive applications.


title: Highlights of EPS HEP 2019

abstract: An opinionated and informal recap of highlights from the EPS HEP 2019 conference in Ghent, including some aspects of flavour physics, neutrinos, high-density QCD, astrophysics and energy frontier collider physics, and some thoughts about the future.


title: H-COUP Version 2: a program for one-loop corrected Higgs boson decays in non-minimal Higgs sectors

abstract: We present the concept of H-COUP_ver 2, which evaluates the decay rates (including higher order corrections) for the Higgs boson with a mass of 125 GeV in various extended Higgs models. In the previous version (H-COUP_1.0), only a full set of the Higgs boson vertices are evaluated at one-loop level in a gauge invariant manner in these models. H-COUP_ver 2 contains all the functions of H-COUP_1.0. After shortly introducing these extended Higgs models and discussing their theoretical and experimental constraints, we summarize formulae for the renormalized vertices and the decay rates. We then explain how to install and run H-COUP_ver 2 with some numerical examples.


title: Skip-Clip: Self-Supervised Spatiotemporal Representation Learning by Future Clip Order Ranking

abstract: Deep neural networks require collecting and annotating large amounts of data to train successfully. In order to alleviate the annotation bottleneck, we propose a novel self-supervised representation learning approach for spatiotemporal features extracted from videos. We introduce Skip-Clip, a method that utilizes temporal coherence in videos, by training a deep model for future clip order ranking conditioned on a context clip as a surrogate objective for video future prediction. We show that features learned using our method are generalizable and transfer strongly to downstream tasks. For action recognition on the UCF101 dataset, we obtain 51.8% improvement over random initialization and outperform models initialized using inflated ImageNet parameters. Skip-Clip also achieves results competitive with state-of-the-art self-supervision methods.


title: Look globally, age locally: Face aging with an attention mechanism

abstract: Face aging is of great importance for cross-age recognition and entertainment-related applications. Recently, conditional generative adversarial networks (cGANs) have achieved impressive results for face aging. Existing cGANs-based methods usually require a pixel-wise loss to keep the identity and background consistent. However, minimizing the pixel-wise loss between the input and synthesized images likely resulting in a ghosted or blurry face. To address this deficiency, this paper introduces an Attention Conditional GANs (AcGANs) approach for face aging, which utilizes attention mechanism to only alert the regions relevant to face aging. In doing so, the synthesized face can well preserve the background information and personal identity without using the pixel-wise loss, and the ghost artifacts and blurriness can be significantly reduced. Based on the benchmarked dataset Morph, both qualitative and quantitative experiment results demonstrate superior performance over existing algorithms in terms of image quality, personal identity, and age accuracy.


title: Extension Technique for Functions of Diffusion Operators: a stochastic approach

abstract: It has recently been shown that complete Bernstein functions of the Laplace operator map the Dirichlet boundary condition of a related elliptic PDE to the Neumann boundary condition. The importance of this mapping consists in being able to convert problems involving non-local operators, like fractional Laplacians, into ones that only involve differential operators. We generalise this result to diffusion operators associated with stochastic differential equations, using a method which is entirely based on stochastic analysis.


title: Diverse properties of Ly-alpha emission in low-redshift compact star-forming galaxies with extremely high [OIII]/[OII] ratios

abstract: We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of eight compact star-forming galaxies at redshifts z=0.02811-0.06540, with low oxygen abundances 12+log(O/H)=7.43-7.82 and extremely high emission-line flux ratios O32=[OIII]5007/[OII]3727~22-39, aiming to study the properties of Ly-alpha emission in such conditions. We find a diversity in Ly-alpha properties. In five galaxies Ly-alpha emission line is strong, with equivalent width (EW) in the range 45-190A. In the remaining galaxies, weak Ly-alpha emission with EW(Ly-alpha)~2-7A is superposed on a broad Ly-alpha absorption line, indicating a high neutral hydrogen column density N(HI)~(1-3)x10^21 cm^-2. We examine the relation between the Ly-alpha escape fraction fesc(Ly-alpha) and the Lyman continuum escape fraction fesc(LyC), using direct measures of the latter in eleven low-redshift LyC leakers, to verify whether fesc(Ly-alpha) can be an indirect measure of escaping LyC radiation. The usefulness of O32, of the Ly-alpha equivalent width EW(Ly-alpha) and of the Ly-alpha peak separation Vsep as indirect indicators of Ly-alpha leakage is also discussed. It is shown that there is no correlation between O32 and fesc(Ly-alpha). We find an increase of fesc(Ly-alpha) with increasing EW(Ly-alpha) for EW(Ly-alpha)<100A, but for higher EW(Ly-alpha)>150A the fesc(Ly-alpha) is nearly constant attaining the value of ~0.25. We find an anticorrelation between fesc(Ly-alpha) and Vsep, though not as tight as the one found earlier between fesc(LyC) and Vsep. This finding makes Vsep a promising indirect indicator of both the Ly-alpha and ionizing radiation leakage.


title: Missing Not at Random in Matrix Completion: The Effectiveness of Estimating Missingness Probabilities Under a Low Nuclear Norm Assumption

abstract: Matrix completion is often applied to data with entries missing not at random (MNAR). For example, consider a recommendation system where users tend to only reveal ratings for items they like. In this case, a matrix completion method that relies on entries being revealed at uniformly sampled row and column indices can yield overly optimistic predictions of unseen user ratings. Recently, various papers have shown that we can reduce this bias in MNAR matrix completion if we know the probabilities of different matrix entries being missing. These probabilities are typically modeled using logistic regression or naive Bayes, which make strong assumptions and lack guarantees on the accuracy of the estimated probabilities. In this paper, we suggest a simple approach to estimating these probabilities that avoids these shortcomings. Our approach follows from the observation that missingness patterns in real data often exhibit low nuclear norm structure. We can then estimate the missingness probabilities by feeding the (always fully-observed) binary matrix specifying which entries are revealed or missing to an existing nuclear-norm-constrained matrix completion algorithm by Davenport et al. [2014]. Thus, we tackle MNAR matrix completion by solving a different matrix completion problem first that recovers missingness probabilities. We establish finite-sample error bounds for how accurate these probability estimates are and how well these estimates debias standard matrix completion losses for the original matrix to be completed. Our experiments show that the proposed debiasing strategy can improve a variety of existing matrix completion algorithms, and achieves downstream matrix completion accuracy at least as good as logistic regression and naive Bayes debiasing baselines that require additional auxiliary information.


title: The conditional censored graphical lasso estimator

abstract: In many applied fields, such as genomics, different types of data are collected on the same system, and it is not uncommon that some of these datasets are subject to censoring as a result of the measurement technologies used, such as data generated by polymerase chain reactions and flow cytometer. When the overall objective is that of network inference, at possibly different levels of a system, information coming from different sources and/or different steps of the analysis can be integrated into one model with the use of conditional graphical models. In this paper, we develop a doubly penalized inferential procedure for a conditional Gaussian graphical model when data can be subject to censoring. The computational challenges of handling censored data in high dimensionality are met with the development of an efficient Expectation-Maximization algorithm, based on approximate calculations of the moments of truncated Gaussian distributions and on a suitably derived two-step procedure alternating graphical lasso with a novel block-coordinate multivariate lasso approach. We evaluate the performance of this approach on an extensive simulation study and on gene expression data generated by RT-qPCR technologies, where we are able to integrate network inference, differential expression detection and data normalization into one model.


title: Normal-Superfluid Phase Separation in Spin-Half Bosons at Finite Temperature

abstract: For pseudospin-half bosons with inter-spin attraction and intra-spin repulsion, normal phase and Bose condensed phase can coexist at finite temperature. The homogeneous system is unstable against the spinodal decomposition within a medium density interval, and consequently, a normal-superfluid phase separation takes place. The isothermal equation-of-state shows a characteristic plateau in the P-V (pressure-volume) diagram, which is reminiscent of a classical gas-liquid transition, although, unlike the latter, the coexistence lines never terminate at a critical point as temperature increases. In a harmonic trap, the phase separation can be revealed by the density profile of the atomic cloud, which exhibits a sudden jump across the phase boundary.


title: Noise in stimulated Raman scattering measurement: From basics to practice

abstract: We revisit laser intensity noise in the context of stimulated Raman scattering (SRS), which has recently proved to be a key technique to provide label free images of chemical bonds in biological and medical samples. Contrary to most microscopy techniques, which detect a weak photon flux resulting from light matter interactions, SRS is a pump-probe scheme that works in the high flux regime and happens as a weak modulation ($10^{-4}-10^{-6}$) in a strong laser field. As a result, laser noise is a key issue in SRS detection. This practical tutorial provides the experimentalists with the tools required to assess the amount of noise and the ultimate SRS detection limit in a conventional lock-in-based SRS system. We first define the quantities that are relevant when discussing intensity noise, and illustrate them through a conventional model of light detection by a photodiode. Stimulated Raman Scattering is then introduced in its lock-in-based implementation, and the model presented is adapted in this particular case. The power spectral density (PSD), relative intensity noise (RIN), signal to noise ratio (SNR), and sensitivity of the system are derived and discussed. Two complementary methods are presented that allow measurement of the RIN and assessment of the performance of a SRS system. Such measurements are illustrated on two commercial laser systems. Finally, the consequences of noise in SRS are discussed, and future developments are suggested. The presentation is made simple enough for under-graduated, graduated students, and newcomers in the field of stimulated Raman, and more generally in pump-probe based schemes.


title: A First-Order Algorithmic Framework for Wasserstein Distributionally Robust Logistic Regression

abstract: Wasserstein distance-based distributionally robust optimization (DRO) has received much attention lately due to its ability to provide a robustness interpretation of various learning models. Moreover, many of the DRO problems that arise in the learning context admits exact convex reformulations and hence can be tackled by off-the-shelf solvers. Nevertheless, the use of such solvers severely limits the applicability of DRO in large-scale learning problems, as they often rely on general purpose interior-point algorithms. On the other hand, there are very few works that attempt to develop fast iterative methods to solve these DRO problems, which typically possess complicated structures. In this paper, we take a first step towards resolving the above difficulty by developing a first-order algorithmic framework for tackling a class of Wasserstein distance-based distributionally robust logistic regression (DRLR) problem. Specifically, we propose a novel linearized proximal ADMM to solve the DRLR problem, whose objective is convex but consists of a smooth term plus two non-separable non-smooth terms. We prove that our method enjoys a sublinear convergence rate. Furthermore, we conduct three different experiments to show its superb performance on both synthetic and real-world datasets. In particular, our method can achieve the same accuracy up to 800+ times faster than the standard off-the-shelf solver.


title: Using braids to quantify interface growth and coherence in a rotor-oscillator flow

abstract: The growth rate of material interfaces is an important proxy for mixing and reaction rates in fluid dynamics, and can also be used to identify regions of coherence. Estimating such growth rates can be difficult, since they depend on detailed properties of the velocity field, such as its derivatives, that are hard to measure directly. When an experiment gives only sparse trajectory data, it is natural to encode planar trajectories as mathematical braids, which are topological objects that contain information on the mixing characteristics of the flow, in particular through their action on topological loops. We test such braid methods on an experimental system, the rotor-oscillator flow, which is well-described by a theoretical model. We conduct a series of laboratory experiments to collect particle tracking and particle image velocimetry data, and use the particle tracks to identify regions of coherence within the flow that match the results obtained from the model velocity field. We then use the data to estimate growth rates of material interface, using both the braid approach and numerical simulations. The interface growth rates follow similar qualitative trends in both the experiment and model, but have significant quantitative differences, suggesting that the two are not as similar as first seems. Our results shows that there are challenges in using the braid approach to analyze data, in particular the need for long trajectories, but that these are not insurmountable.


title: Non-Centralized Navigation for Source Localization by Cooperative UAVs

abstract: In this paper, we propose a distributed solution to the navigation of a population of unmanned aerial vehicles (UAVs) to best localize a static source. The network is considered heterogeneous with UAVs equipped with received signal strength (RSS) sensors from which it is possible to estimate the distance from the source and/or the direction of arrival through ad-hoc rotations. This diversity in gathering and processing RSS measurements mitigates the loss of localization accuracy due to the adoption of low-complexity sensors. The UAVs plan their trajectories on-the-fly and in a distributed fashion. The collected data are disseminated through the network via multi-hops, therefore being subject to latency. Since not all the paths are equal in terms of information gathering rewards, the motion planning is formulated as a minimization of the uncertainty of the source position under UAV kinematic and anti-collision constraints and performed by 3D non-linear programming. The proposed analysis takes into account non-line-of-sight (NLOS) channel conditions as well as measurement age caused by the latency constraints in communication.


title: Empirical Analysis of Session-Based Recommendation Algorithms

abstract: Recommender systems are tools that support online users by pointing them to potential items of interest in situations of information overload. In recent years, the class of session-based recommendation algorithms received more attention in the research literature. These algorithms base their recommendations solely on the observed interactions with the user in an ongoing session and do not require the existence of long-term preference profiles. Most recently, a number of deep learning based ("neural") approaches to session-based recommendations were proposed. However, previous research indicates that today's complex neural recommendation methods are not always better than comparably simple algorithms in terms of prediction accuracy. With this work, our goal is to shed light on the state-of-the-art in the area of session-based recommendation and on the progress that is made with neural approaches. For this purpose, we compare twelve algorithmic approaches, among them six recent neural methods, under identical conditions on various datasets. We find that the progress in terms of prediction accuracy that is achieved with neural methods is still limited. In most cases, our experiments show that simple heuristic methods based on nearest-neighbors schemes are preferable over conceptually and computationally more complex methods. Observations from a user study furthermore indicate that recommendations based on heuristic methods were also well accepted by the study participants. To support future progress and reproducibility in this area, we publicly share the session-rec evaluation framework that was used in our research.


title: Zeta functions with respect to general coined quantum walk of periodic graphs

abstract: We define a zeta function of a graph by using the time evolution matrix of a general coined quantum walk on it, and give a determinant expression for the zeta function of a finite graph. Furthermore, we present a determinant expression for the zeta function of an (infinite) periodic graph.


title: Distributed Networked Learning with Correlated Data

abstract: We consider a distributed estimation method in a setting with heterogeneous streams of correlated data distributed across nodes in a network. In the considered approach, linear models are estimated locally (i.e., with only local data) subject to a network regularization term that penalizes a local model that differs from neighboring models. We analyze computation dynamics (associated with stochastic gradient updates) and information exchange (associated with exchanging current models with neighboring nodes). We provide a finite-time characterization of convergence of the weighted ensemble average estimate and compare this result to federated learning, an alternative approach to estimation wherein a single model is updated by locally generated gradient updates. This comparison highlights the trade-off between speed vs precision: while model updates take place at a faster rate in federated learning, the proposed networked approach to estimation enables the identification of models with higher precision. We illustrate the method's general applicability in two examples: estimating a Markov random field using wireless sensor networks and modeling prey escape behavior of flocking birds based on a publicly available dataset.


title: Complex Spacing Ratios: A Signature of Dissipative Quantum Chaos

abstract: We introduce a complex-plane generalization of the consecutive level-spacing distribution, used to distinguish regular from chaotic quantum spectra. Our approach features the distribution of complex-valued ratios between nearest- and next-to-nearest neighbor spacings. We show that this quantity can successfully detect the chaotic or regular nature of complex-valued spectra. This is done in two steps. First, we show that, if eigenvalues are uncorrelated, the distribution of complex spacing ratios is flat within the unit circle, whereas random matrices show a strong angular dependence in addition to the usual level repulsion. The universal fluctuations of Gaussian Unitary and Ginibre Unitary universality classes in the large-matrix-size limit are shown to be well described by Wigner-like surmises for small-size matrices with eigenvalues on the circle and on the two-torus, respectively. To study the latter case, we introduce the Toric Unitary Ensemble, characterized by a flat joint eigenvalue distribution on the two-torus. Second, we study different physical situations where nonhermitian matrices arise: dissipative quantum systems described by a Lindbladian, non-unitary quantum dynamics described by nonhermitian Hamiltonians, and classical stochastic processes. We show that known integrable models have a flat distribution of complex spacing ratios whereas generic cases, expected to be chaotic, conform to Random Matrix Theory predictions. Specifically, we were able to clearly distinguish chaotic from integrable dynamics in boundary-driven dissipative spin-chain Liouvillians and in the classical asymmetric simple exclusion process and to differentiate localized from delocalized phases in a nonhermitian disordered many-body system.

categories: cond-mat.stat-mech math-ph math.MP nlin.CD quant-ph


title: Detection of a Water Tracer in Interstellar Comet 2I/Borisov

abstract: We present high spectral resolution optical spectra obtained with the ARCES instrument at Apache Point Observatory showing detection of the [OI]6300 A line in interstellar comet 2I/Borisov. We employ the observed flux in this line to derive an H$_2$O production rate of (6.3$\pm$1.5)$\times$10$^{26}$ mol s$^{-1}$. Comparing to previously reported observations of CN, this implies a CN/H$_2$O ratio of $\sim$0.3-0.6%. The lower end of this range is consistent with the average value in comets, while the upper end is higher than the average value for Solar System comets, but still within the range of observed values. C$_2$/H$_2$O is depleted, with a value likely less than 0.1%. The dust-to-gas ratio is consistent with the normal value for Solar System comets. Using a simple sublimation model we estimate an H$_2$O active area of 1.7 km$^2$, which for current estimates for the size of Borisov suggests active fractions between 1-150%, consistent with values measured in Solar System comets. More detailed characterization of 2I/Borisov, including compositional information and properties of the nucleus, is needed to fully interpret the observed H$_2$O production rate.


title: Infrared spectroscopy of the recent outburst in V1047 Cen (Nova Centauri 2005)

abstract: Fourteen years after its eruption as a classical nova (CN), V1047 Cen (Nova Cen 2005) began an unusual re-brightening in 2019 April. The amplitude of the brightening suggests that this is a dwarf nova (DN) eruption in a CN system. Very few CNe have had DN eruptions within decades of the main CN outburst. The 14 years separating the CN and DN eruptions of V1047 Cen is the shortest of all instances recorded thus far. Explaining this rapid succession of CN and DN outbursts in V1047 Cen may be challenging within the framework of standard theories for DN outbursts. Following a CN eruption, the mass accretion rate is believed to remain high $(\dot{M}\sim10^{-8}$M$_\odot$yr$^{-1})$ for a few centuries, due to the irradiation of the secondary star by the still-hot surface of the white dwarf. Thus a DN eruption is not expected to occur during this high mass accretion phase as DN outbursts, which result from thermal instabilities in the accretion disk, and arise during a regime of low mass accretion rate $(\dot{M}\sim10^{-10}$M$_\odot$yr$^{-1})$. Here we present near-infrared spectroscopy to show that the present outburst is most likely a DN eruption, and discuss the possible reasons for its early occurrence. Even if the present re-brightening is later shown to be due to a cause other than a DN outburst, the present study provides invaluable documentation of this unusual event.


title: Weak equivalence of higher-dimensional automata

abstract: This paper introduces a notion of equivalence for higher-dimensional automata, called weak equivalence. Weak equivalence focuses mainly on a traditional trace language and a new homology language, which captures the overall independence structure of an HDA. It is shown that weak equivalence is compatible with both the tensor product and the coproduct of HDAs and that, under certain conditions, HDAs may be reduced to weakly equivalent smaller ones by merging and collapsing cubes.


title: The Zariski closure of integral points on varieties parametrizing periodic continued fractions

abstract: Let $R$ be the ring of $S$-integers in a number field $K$. Let $\mathcal{B}=\{\beta, \beta^{\ast}\}$ be the multi-set of roots of a nonzero quadratic polynomial over $R$. There are varieties $V(\mathcal{B})_{N,k}$ defined over $R$ parametrizing periodic continued fractions $[b_1,\ldots , b_N,\overline{a_1,\ldots ,a_k}]$ for $\beta$ or $\beta^{\ast}$. We study the $R$-points on these varieties, finding contrasting behavior according to whether groups of units are infinite or not. If $R$ is the rational integers or the ring of integers in an imaginary quadratic field, we prove that the $R$-points of $V(\mathcal{B})_{N,k}$ are not Zariski dense. On the other hand, suppose that $\beta\not\in K\cup\{\infty\}$, $R^\times$ is infinite, and that there are infinitely many units in the (left) order $R_\beta$ of $\beta R+R\subseteq K(\beta)$ with norm to $K$ equal to $(-1)^k$. Then we prove that the $R$-points on $V(\mathcal{B})_{1,k}$ are Zariski dense for $k\geq 8$ and the $R$-points on $V(\mathcal{B})_{0,k}$ are Zariski dense for $k\geq 9$. We also prove that $V(\mathcal{B})_{1,k}$ and $V(\mathcal{B})_{0,k}$ are $K$-rational irreducible varieties for $k$ sufficiently large.


title: Constructions of k-uniform and absolutely maximally entangled states beyond maximum distance codes

abstract: Pure multipartite quantum states of n parties and local dimension q are called k-uniform if all reductions to k parties are maximally mixed. These states are relevant for our understanding of multipartite entanglement, quantum information protocols, and the construction of quantum error correction codes. To our knowledge, the only known systematic construction of these quantum states is based on classical error correction codes. We present a systematic method to construct other examples of k-uniform states and show that the states derived through our construction are not equivalent to any k-uniform state constructed from the so-called maximum distance separable error correction codes. Furthermore, we use our method to construct several examples of absolutely maximally entangled states whose existence was open so far.


title: Permutations encoding the local shape of level curves of real polynomials via generic projections

abstract: The non-convexity of a smooth and compact connected component of a real algebraic plane curve can be measured by a combinatorial object called the Poincare-Reeb tree associated to the curve and to a direction of projection. In this paper we show that if the chosen projection avoids the bitangents and the inflectional tangencies to the small enough level curves of a real bivariate polynomial function near a strict local minimum at the origin, then the asymptotic Poincare-Reeb tree becomes a complete binary tree and its vertices become endowed with a total order relation. Such a projection direction is called generic. We prove that for any such asymptotic family of level curves, there are finitely many intervals on the real projective line outside of which all the directions are generic with respect to all the curves in the family. If the choice of the direction of projection is generic, then the local shape of the curves can be encoded in terms of alternating permutations, that we call snakes. The snakes offer an effective description of the local geometry and topology, well-suited for further computations.


title: Off-Shell CHY Amplitudes and Feynman Graphs

abstract: A polynomial form is established for the off-shell CHY scattering equations proposed by Lam and Yao. Re-expressing this in terms of independent Mandelstam invariants provides a new expression for the polynomial scattering equations, immediately valid off shell, which makes it evident that they yield the off-shell amplitudes given by massless $\phi^3$ Feynman graphs. A CHY expression for individual Feynman graphs, valid even off shell, is established through a recurrence relation.


title: A vector valued almost sure invariance principle for time dependent non-uniformly expanding dynamical systems

abstract: We prove a vector-valued almost sure invariance principle for some classes of time dependent non-uniformly distance expanding dynamical systems. The models we have in mind are certain sequential versions of the smooth non-uniformly distance expanding maps considered in \cite{castro} and \cite{Vara}, as well as certain types of sequences of covering maps. Our results rely on the theory of complex projective metrics which was developed in \cite{Rug}, together with the spectral methods of Gou\"ezel \cite{GO}. A big advantage in applying the theory of complex cones here is that it also yields additional probabilistic limit theorems for random dynamical systems, as described at the last section of this paper.


title: Bounded degree complexes of forests

abstract: Given an arbitrary sequence of non-negative integers $\vec{\lambda}=(\lambda_1,\dots,\lambda_n)$ and a graph $G$ with vertex set $\{v_1,\dots,v_n\}$, the bounded degree complex, denoted $\text{BD}^{\vec{\lambda}}(G)$, is a simplicial complex whose faces are the subsets $H\subseteq E(G)$ such that for each $i \in \{1,\dots,n\}$, the degree of vertex $v_i$ in the induced subgraph $G[H]$ is at most $\lambda_i$. When $\lambda_i=k$ for all $i$, the bounded degree complex $\text{BD}^{\vec{\lambda}}(G)$ is called the $k$-matching complex, denoted $M_k(G)$. In this article, we determine the homotopy type of bounded degree complexes of forests. In particular, we show that, for all $k\geq 1$, the $k$-matching complexes of caterpillar graphs are either contractible or homotopy equivalent to a wedge of spheres, thereby proving a conjecture of Julianne Vega \cite[Conjecture 7.3]{Vega19}. We also give a closed form formula for the homotopy type of the bounded degree complexes of those caterpillar graphs in which every non-leaf vertex is adjacent to at least one leaf vertex.


title: Stein Variational Gradient Descent With Matrix-Valued Kernels

abstract: Stein variational gradient descent (SVGD) is a particle-based inference algorithm that leverages gradient information for efficient approximate inference. In this work, we enhance SVGD by leveraging preconditioning matrices, such as the Hessian and Fisher information matrix, to incorporate geometric information into SVGD updates. We achieve this by presenting a generalization of SVGD that replaces the scalar-valued kernels in vanilla SVGD with more general matrix-valued kernels. This yields a significant extension of SVGD, and more importantly, allows us to flexibly incorporate various preconditioning matrices to accelerate the exploration in the probability landscape. Empirical results show that our method outperforms vanilla SVGD and a variety of baseline approaches over a range of real-world Bayesian inference tasks.


title: Learning Data Manipulation for Augmentation and Weighting

abstract: Manipulating data, such as weighting data examples or augmenting with new instances, has been increasingly used to improve model training. Previous work has studied various rule- or learning-based approaches designed for specific types of data manipulation. In this work, we propose a new method that supports learning different manipulation schemes with the same gradient-based algorithm. Our approach builds upon a recent connection of supervised learning and reinforcement learning (RL), and adapts an off-the-shelf reward learning algorithm from RL for joint data manipulation learning and model training. Different parameterization of the "data reward" function instantiates different manipulation schemes. We showcase data augmentation that learns a text transformation network, and data weighting that dynamically adapts the data sample importance. Experiments show the resulting algorithms significantly improve the image and text classification performance in low data regime and class-imbalance problems.


title: Recurrence of a Weighted Random Walk on a Circle Packing with Parabolic Carrier

abstract: In this paper we show that given a circle packing of an infinite planar triangulation such that its carrier is parabolic, placing weights on the edges according to a certain natural way introduced by Dubejko, makes the random walk recurrent. We also propose a higher-dimensional analogue of the Dubejko weights.


title: Testing Equivalence of Clustering

abstract: In this paper, we test whether two datasets share a common clustering structure. As a leading example, we focus on comparing clustering structures in two independent random samples from two mixtures of multivariate normal distributions. Mean parameters of these normal distributions are treated as potentially unknown nuisance parameters and are allowed to differ. Assuming knowledge of mean parameters, we first determine the phase diagram of the testing problem over the entire range of signal-to-noise ratios by providing both lower bounds and tests that achieve them. When nuisance parameters are unknown, we propose tests that achieve the detection boundary adaptively as long as ambient dimensions of the datasets grow at a sub-linear rate with the sample size.


title: Arithmetic conjectures suggested by the statistical behavior of modular symbols

abstract: Suppose $E$ is an elliptic curve over $\mathbb{Q}$ and $\chi$ is a Dirichlet character. We use statistical properties of modular symbols to estimate heuristically the probability that $L(E,\chi,1) = 0$. Via the Birch and Swinnerton-Dyer conjecture, this gives a heuristic estimate of the probability that the Mordell-Weil rank grows in abelian extensions of $\mathbb{Q}$. Using this heuristic we find a large class of infinite abelian extensions $F$ where we expect $E(F)$ to be finitely generated. Our work was inspired by earlier conjectures (based on random matrix heuristics) due to David, Fearnley, and Kisilevsky. Where our predictions and theirs overlap, the predictions are consistent.


title: Deep learning is adaptive to intrinsic dimensionality of model smoothness in anisotropic Besov space

abstract: Deep learning has exhibited superior performance for various tasks, especially for high-dimensional datasets, such as images. To understand this property, we investigate the approximation and estimation ability of deep learning on anisotropic Besov spaces. The anisotropic Besov space is characterized by direction-dependent smoothness and includes several function classes that have been investigated thus far. We demonstrate that the approximation error and estimation error of deep learning only depend on the average value of the smoothness parameters in all directions. Consequently, the curse of dimensionality can be avoided if the smoothness of the target function is highly anisotropic. Unlike existing studies, our analysis does not require a low-dimensional structure of the input data. We also investigate the minimax optimality of deep learning and compare its performance with that of the kernel method (more generally, linear estimators). The results show that deep learning has better dependence on the input dimensionality if the target function possesses anisotropic smoothness, and it achieves an adaptive rate for functions with spatially inhomogeneous smoothness.


title: Attenuating Random Noise in Seismic Data by a Deep Learning Approach

abstract: In the geophysical field, seismic noise attenuation has been considered as a critical and long-standing problem, especially for the pre-stack data processing. Here, we propose a model to leverage the deep-learning model for this task. Rather than directly applying an existing de-noising model from ordinary images to the seismic data, we have designed a particular deep-learning model, based on residual neural networks. It is named as N2N-Seismic, which has a strong ability to recover the seismic signals back to intact condition with the preservation of primary signals. The proposed model, achieving with great success in attenuating noise, has been tested on two different seismic datasets. Several metrics show that our method outperforms conventional approaches in terms of Signal-to-Noise-Ratio, Mean-Squared-Error, Phase Spectrum, etc. Moreover, robust tests in terms of effectively removing random noise from any dataset with strong and weak noises have been extensively scrutinized in making sure that the proposed model is able to maintain a good level of adaptation while dealing with large variations of noise characteristics and intensities.


title: Charge Density Wave and Superconductivity in Transition Metal Dichalcogenides

abstract: Competing orders in condensed matter give rise to the emergence of fascinating, new phenomena. Here, we investigate the competition between superconductivity and charge density wave in the context of layered-metallic compounds, transition metal dichalcogenides, in which the superconducting state is usually suppressed by the charge density wave. We show, using real-space self-consistent Bogoliubov-de Gennes calculations and momentum-space calculations involving density-functional theory and dynamical mean-field theory, that there is a surprising reappearance of superconductivity in the presence of non-magnetic disorder fluctuations, as observed in recent experiments.


title: Model-Free Mean-Field Reinforcement Learning: Mean-Field MDP and Mean-Field Q-Learning

abstract: We study infinite horizon discounted Mean Field Control (MFC) problems with common noise through the lens of Mean Field Markov Decision Processes (MFMDP). We allow the agents to use actions that are randomized not only at the individual level but also at the level of the population. This common randomization allows us to establish connections between both closed-loop and open-loop policies for MFC and Markov policies for the MFMDP. In particular, we show that there exists an optimal closed-loop policy for the original MFC. Building on this framework and the notion of state-action value function, we then propose reinforcement learning (RL) methods for such problems, by adapting existing tabular and deep RL methods to the mean-field setting. The main difficulty is the treatment of the population state, which is an input of the policy and the value function. We provide convergence guarantees for tabular algorithms based on discretizations of the simplex. Neural network based algorithms are more suitable for continuous spaces and allow us to avoid discretizing the mean field state space. Numerical examples are provided.


title: Homotopy Types of Random Cubical Complexes

abstract: We study the topology of a random cubical complex associated to Bernoulli site percolation on a cubical grid. We begin by establishing a limit law for homotopy types. More precisely, looking within an expanding window, we define a sequence of normalized counting measures (counting connected components according to homotopy type), and we show that this sequence of random probability measures converges in probability to a deterministic probability measure. We then investigate the dependence of the limiting homotopy measure on the coloring probability $p$, and our results show a qualitative change in the homotopy measure as $p$ crosses the percolation threshold $p=p_c$. Specializing to the case of $d=2$ dimensions, we also present empirical results that raise further questions on the $p$-dependence of the limiting homotopy measure.


title: Joint Indoor Localization and Navigation of UAVs for Network Formation Control

abstract: In this paper, we propose a joint indoor localization and navigation algorithm to enable a swarm of unmanned aerial vehicles (UAVs) to deploy in a specific spatial formation in indoor environments. In the envisioned scenario, we consider a static user acting as a central unit whose main task is to acquire all the UAV measurements carrying position-dependent information and to estimate the UAV positions when there is no existing infrastructure for positioning. Subsequently, the user exploits the estimated positions as inputs for the navigation control with the aim of deploying the UAVs in a desired formation in space (formation shaping). The user plans the trajectory of each UAV in real time, guaranteeing a safe navigation in the presence of obstacles. The proposed algorithm guides the UAVs to their desired final locations with good accuracy.


title: Relaxation Effects in Twisted Bilayer Graphene: a Multi-Scale Approach

abstract: We present a multi-scale density functional theory (DFT) informed molecular dynamics and tight-binding (TB) approach to capture the interdependent atomic and electronic structures of twisted bilayer graphene. We calibrate the flat band magic angle to be at $\theta_{\rm M} = 1.08^{\circ}$ by rescaling the interlayer tunneling for different atomic structure relaxation models as a way to resolve the indeterminacy of existing atomic and electronic structure models whose predicted magic angles vary widely between $0.9^\circ \sim 1.3^\circ$. The interatomic force fields are built using input from various stacking and interlayer distance dependent DFT total energies including the exact exchange and random phase approximation (EXX+RPA). We use a Fermi velocity of $\upsilon_{\rm F} \simeq 10^{6}$~m/s for graphene that is enhanced by about $\sim 15\%$ over the local density approximation (LDA) values. Based on this atomic and electronic structure model we obtain high-resolution spectral functions comparable with experimental angle-resolved photoemission spectra (ARPES). Our analysis of the interdependence between the atomic and electronic structures indicates that the intralayer elastic parameters compatible with the DFT-LDA, which are stiffer by $\sim$30\% than widely used reactive empirical bond order force fields, can combine with EXX+RPA interlayer potentials to yield the magic angle at $\sim 1.08^{\circ}$ without further rescaling of the interlayer tunneling.


title: An Ensemble Approach toward Automated Variable Selection for Network Anomaly Detection

abstract: While variable selection is essential to optimize the learning complexity by prioritizing features, automating the selection process is preferred since it requires laborious efforts with intensive analysis otherwise. However, it is not an easy task to enable the automation due to several reasons. First, selection techniques often need a condition to terminate the reduction process, for example, by using a threshold or the number of features to stop, and searching an adequate stopping condition is highly challenging. Second, it is uncertain that the reduced variable set would work well; our preliminary experimental result shows that well-known selection techniques produce different sets of variables as a result of reduction (even with the same termination condition), and it is hard to estimate which of them would work the best in future testing. In this paper, we demonstrate the potential power of our approach to the automation of selection process that incorporates well-known selection methods identifying important variables. Our experimental results with two public network traffic data (UNSW-NB15 and IDS2017) show that our proposed method identifies a small number of core variables, with which it is possible to approximate the performance to the one with the entire variables.


title: Better Exploration with Optimistic Actor-Critic

abstract: Actor-critic methods, a type of model-free Reinforcement Learning, have been successfully applied to challenging tasks in continuous control, often achieving state-of-the art performance. However, wide-scale adoption of these methods in real-world domains is made difficult by their poor sample efficiency. We address this problem both theoretically and empirically. On the theoretical side, we identify two phenomena preventing efficient exploration in existing state-of-the-art algorithms such as Soft Actor Critic. First, combining a greedy actor update with a pessimistic estimate of the critic leads to the avoidance of actions that the agent does not know about, a phenomenon we call pessimistic underexploration. Second, current algorithms are directionally uninformed, sampling actions with equal probability in opposite directions from the current mean. This is wasteful, since we typically need actions taken along certain directions much more than others. To address both of these phenomena, we introduce a new algorithm, Optimistic Actor Critic, which approximates a lower and upper confidence bound on the state-action value function. This allows us to apply the principle of optimism in the face of uncertainty to perform directed exploration using the upper bound while still using the lower bound to avoid overestimation. We evaluate OAC in several challenging continuous control tasks, achieving state-of the art sample efficiency.


title: Modelling X-ray RMS spectra I: intrinsically variable AGN

abstract: We present simple XSPEC models for fitting excess variance spectra of AGN. Using a simple Monte-Carlo approach, we simulate a range of spectra corresponding to physical parameters varying, then calculate the resulting variance spectra. Starting from a variable power-law, we build up a set of models corresponding to the different physical processes that can affect the final excess variance spectrum. We show that the complex excess variance spectrum of IRAS 13224-3809 can be well described by such an intrinsic variability model, where the power-law variability is damped by relativistic reflection and enhanced by an ultra fast outflow. The reflection flux is correlated with that of the power-law, but not perfectly. We argue that this correlation is stronger at high frequencies, where reverberation lags are detected, while excess variance spectra are typically dominated by low frequency variability.


title: Minimax Weight and Q-Function Learning for Off-Policy Evaluation

abstract: We provide theoretical investigations into off-policy evaluation in reinforcement learning using function approximators for (marginalized) importance weights and value functions. Our contributions include: (1) A new estimator, MWL, that directly estimates importance ratios over the state-action distributions, removing the reliance on knowledge of the behavior policy as in prior work (Liu et al., 2018). (2) Another new estimator, MQL, obtained by swapping the roles of importance weights and value-functions in MWL. MQL has an intuitive interpretation of minimizing average Bellman errors and can be combined with MWL in a doubly robust manner. (3) Several additional results that offer further insights into these methods, including the sample complexity analyses of MWL and MQL, their asymptotic optimality in the tabular setting, how the learned importance weights depend the choice of the discriminator class, and how our methods provide a unified view of some old and new algorithms in RL.


title: Automatic Mapping of the Indoor World with Personal Radars

abstract: Digital maps will revolutionize our experience of perceiving and navigating indoor environments. While today we rely only on the representation of the outdoors, the mapping of indoors is mainly a part of the traditional SLAM problem where robots discover the surrounding and perform self-localization. Nonetheless, robot deployment prevents from a large diffusion and fast mapping of indoors and, further, they are usually equipped with laser and vision technology that fail in scarce visibility conditions. To this end, a possible solution is to turn future personal devices into personal radars as a milestone towards the automatic generation of indoor maps using massive array technology at millimeter-waves, already in place for communications. In this application-oriented paper, we will describe the main achievements attained so far to develop the personal radar concept, using ad-hoc collected experimental data, and by discussing possible future directions of investigation.


title: Classical tests on a charged Weyl black hole: bending of light, Shapiro delay and Sagnac effect

abstract: In this paper, we apply the classical test of general relativity on a charged Weyl black hole, whose exterior geometry is defined by altering the spherically symmetric solutions of Weyl conformal theory of gravity. The tests are basically founded on scrutinizing the angular geodesics of light rays propagating in the gravitating system caused by the black hole. In this investigation, we bring detailed discussions about the bending of light, together with two other relativistic effects, known as the Shapiro and the Sagnac effects. We show that the results are in good conformity with the general relativistic effects, besides giving long-distance corrections caused by the cosmological nature of the background geometry under study.


title: Pauls rectifiable and purely Pauls unrectifiable smooth hypersurfaces

abstract: This paper is related to the problem of finding a good notion of rectifiability in sub-Riemannian geometry. In particular, we study which kind of results can be expected for smooth hypersurfaces in Carnot groups. Our main contribution will be a consequence of the following result: there exists a $C^{\infty}$ hypersurface $S$ without characteristic points that has uncountably many pairwise non-isomorphic tangent groups on every positive-measure subset. The example is found in a Carnot group of topological dimension 8, it has Hausdorff dimension 12 and so we use on it the Hausdorff measure $\mathcal{H}^{12}$. As a consequence, we show that for every Carnot group of Hausdorff dimension 12, any Lipschitz map defined on a subset of it with values in $S$ has $\mathcal{H}^{12}$-null image. In particular, we deduce that this smooth hypersurface cannot be Lipschitz parametrizable by countably many maps each defined on some subset of some Carnot group of Hausdorff dimension $12$. As main consequence we have that a notion of rectifiability proposed by S.Pauls is not equivalent to one proposed by B.Franchi, R.Serapioni and F.Serra Cassano, at least for arbitrary Carnot groups. In addition, we show that, given a subset $U$ of a homogeneous subgroup of Hausdorff dimension $12$ of a Carnot group, every bi-Lipschitz map $f:U\to S$ satisfies $\mathcal{H}^{12}(f(U))=0$. Finally, we prove that such an example does not exist in Heisenberg groups: we prove that all $C^{\infty}$-hypersurfaces in $\mathbb H^n$ with $n\geq 2$ are countably $\mathbb{H}^{n-1}\times\mathbb R$-rectifiabile according to Pauls' definition, even with bi-Lipschitz maps.


title: NGC6240: A triple nucleus system in the advanced or final state of merging

abstract: NGC6240 is a well-studied nearby galaxy system in the process of merging. Based on optical, X-ray, and radio observations, it is thought to harbor two active nuclei We carried out a detailed optical 3D spectroscopic study to investigate the inner region of this system in connection with existing MERLIN and VLBA data. We observed NGC6240 with very high spatial resolution using the MUSE instrument in the Narrow-Field Mode with the four-laser GALACSI adaptive optics system on the ESO VLT under seeing conditions of 0.49 arcsec. Our 3D spectra cover the wavelength range from 4725 to 9350 AA at a spatial resolution of ~75 mas. We report the discovery of three nuclei in the final state of merging within a region of only 1 kpc in the NGC6240 system.Thanks to MUSE we are able to show that the formerly unresolved southern component actually consists of two distinct nuclei separated by only 198 pc. In combination with Gaia data we reach an absolute positional accuracy of only 30 mas that is essential to compare optical spectra with MERLIN and VLBA radio positions. The verification and detailed study of a system with three nuclei, two of which are active and each with a mass in excess of $9\times10^{7} M_{\odot}$, is of great importance for the understanding of hierarchical galaxy formation via merging processes since multiple mergers lead to a faster evolution of massive galaxies in comparison to binary mergers. So far it has been suggested that the formation of galactic nuclei with multiple supermassive black holes (SMBHs) is expected to be rare in the local universe.Triple massive black hole systems might be of fundamental importance for the coalescence of massive black hole binaries in less than a Hubble time leading to the loudest sources of gravitational waves in the millihertz regime.


title: Dynamic Radar Networks of UAVs

abstract: In the incoming years, the low-aerial space will be crowded by unmanned aerial vehicles (UAVs), which will be providing different services. In this expected context, an emerging problem is to detect and track unauthorized or malicious mini/micro UAVs. In contrast to current solutions mainly based on fixed terrestrial radars, this article will put forth the idea of a dynamic radar network (DRN) composed of UAVs able to smartly adapt their formation-navigation control to best track malicious UAVs in real-time, with high accuracy, and in a distributed fashion. To this end, some technological solutions and the main methods for target detection and tracking will be described. Further, an optimized navigation scheme will be developed according to an information-seeking approach. Some examples of simulation results and future directions of work will be finally presented highlighting the advantages of dynamic and reconfigurable networks over static ones.


title: Approximate Bayesian Computation with the Sliced-Wasserstein Distance

abstract: Approximate Bayesian Computation (ABC) is a popular method for approximate inference in generative models with intractable but easy-to-sample likelihood. It constructs an approximate posterior distribution by finding parameters for which the simulated data are close to the observations in terms of summary statistics. These statistics are defined beforehand and might induce a loss of information, which has been shown to deteriorate the quality of the approximation. To overcome this problem, Wasserstein-ABC has been recently proposed, and compares the datasets via the Wasserstein distance between their empirical distributions, but does not scale well to the dimension or the number of samples. We propose a new ABC technique, called Sliced-Wasserstein ABC and based on the Sliced-Wasserstein distance, which has better computational and statistical properties. We derive two theoretical results showing the asymptotical consistency of our approach, and we illustrate its advantages on synthetic data and an image denoising task.