title: Unraveling the Excitation Spectrum of Many-Body Systems from Quantum Quenches

abstract: Quenches are now routinely used in synthetic quantum systems to study a variety of fundamental effects, including ergodicity breaking, light-cone-like spreading of information, and dynamical phase transitions. It was shown recently that the dynamics of equal-time correlators may be related to ground-state phase transitions and some properties of the system excitations. Here, we show that the full low-lying excitation spectrum of a generic many-body quantum system can be extracted from the after-quench dynamics of equal-time correlators. We demonstrate it for a variety of one-dimensional lattice models amenable to exact numerical calculations, including Bose and spin models, with short- or long-range interactions. The approach also applies to higher dimensions, correlated fermions, and continuous models. We argue that it provides an alternative approach to standard pump-probe spectroscopic methods and discuss its advantages.


title: Farkas layers: don't shift the data, fix the geometry

abstract: Successfully training deep neural networks often requires either batch normalization, appropriate weight initialization, both of which come with their own challenges. We propose an alternative, geometrically motivated method for training. Using elementary results from linear programming, we introduce Farkas layers: a method that ensures at least one neuron is active at a given layer. Focusing on residual networks with ReLU activation, we empirically demonstrate a significant improvement in training capacity in the absence of batch normalization or methods of initialization across a broad range of network sizes on benchmark datasets.


title: Magnonics and Supermagnonics

abstract: The magnetic community continues to discuss the possibility to observe the magnetic superfluidity, despite the fact that it has been discovered long time ago. It was observed in antiferromagnetic states of superfluid $^3$He in 1984. In this article we reminds the main principles of spin superfluidity and related Bose-Einstain magnon condensation. We discuss applications of this phenomenon in supermagnonic devises.


title: On Rational Invariant Summation of p-Adic Power Series with Binomial Coefficient

abstract: In the work we have considered p-adic functional series with binomial coefficients and discussed its p-adic convergence. Then we have derived a recurrence relation following with a summation formula which is invariant for rational argument. More precisely, we have investigated a sufficient condition under which the p-adic power series converges to a rational invariant sum 0. Finally we have shown application of the invariant summation formula to get some intersting relations involving Bernoulli numbers and Bernoulli polynomials.


title: Frame Soft Shrinkage as Proximity Operator

abstract: Let $\mathcal H$ and $\mathcal K$ be real Hilbert spaces and $T \in \mathcal{B} (\mathcal H,\mathcal K)$ an injective operator with closed range and Moore-Penrose inverse $T^\dagger$. Based on the well-known characterization of proximity operators by Moreau, we prove that for any proximity operator $\text{Prox} \colon \mathcal K \to \mathcal K$ the operator $T^\dagger \, \text{Prox} \, T$ is a proximity operator on the linear space $\mathcal H$ equipped with a suitable norm. In particular, it follows for the frequently applied soft shrinkage operator $\text{Prox} = S_{\lambda}\colon \ell_2 \rightarrow \ell_2$ and any frame analysis operator $T\colon \mathcal H \to \ell_2$, that the frame shrinkage operator $T^\dagger\, S_\lambda\, T$ is a proximity operator in a suitable Hilbert space.


title: DeshadowGAN: A Deep Learning Approach to Remove Shadows from Optical Coherence Tomography Images

abstract: Purpose: To remove retinal shadows from optical coherence tomography (OCT) images of the optic nerve head(ONH). Methods:2328 OCT images acquired through the center of the ONH using a Spectralis OCT machine for both eyes of 13 subjects were used to train a generative adversarial network (GAN) using a custom loss function. Image quality was assessed qualitatively (for artifacts) and quantitatively using the intralayer contrast: a measure of shadow visibility ranging from 0 (shadow-free) to 1 (strong shadow) and compared to compensated images. This was computed in the Retinal Nerve Fiber Layer (RNFL), the Inner Plexiform Layer (IPL), the Photoreceptor layer (PR) and the Retinal Pigment Epithelium (RPE) layers. Results: Output images had improved intralayer contrast in all ONH tissue layers. On average the intralayer contrast decreased by 33.7$\pm$6.81%, 28.8$\pm$10.4%, 35.9$\pm$13.0%, and43.0$\pm$19.5%for the RNFL, IPL, PR, and RPE layers respectively, indicating successful shadow removal across all depths. This compared to 70.3$\pm$22.7%, 33.9$\pm$11.5%, 47.0$\pm$11.2%, 26.7$\pm$19.0%for compensation. Output images were also free from artifacts commonly observed with compensation. Conclusions: DeshadowGAN significantly corrected blood vessel shadows in OCT images of the ONH. Our algorithm may be considered as a pre-processing step to improve the performance of a wide range of algorithms including those currently being used for OCT image segmentation, denoising, and classification. Translational Relevance: DeshadowGAN could be integrated to existing OCT devices to improve the diagnosis and prognosis of ocular pathologies.


title: Combining docking pose rank and structure with deep learning improves protein-ligand binding mode prediction

abstract: We present a simple, modular graph-based convolutional neural network that takes structural information from protein-ligand complexes as input to generate models for activity and binding mode prediction. Complex structures are generated by a standard docking procedure and fed into a dual-graph architecture that includes separate sub-networks for the ligand bonded topology and the ligand-protein contact map. This network division allows contributions from ligand identity to be distinguished from effects of protein-ligand interactions on classification. We show, in agreement with recent literature, that dataset bias drives many of the promising results on virtual screening that have previously been reported. However, we also show that our neural network is capable of learning from protein structural information when, as in the case of binding mode prediction, an unbiased dataset is constructed. We develop a deep learning model for binding mode prediction that uses docking ranking as input in combination with docking structures. This strategy mirrors past consensus models and outperforms the baseline docking program in a variety of tests, including on cross-docking datasets that mimic real-world docking use cases. Furthermore, the magnitudes of network predictions serve as reliable measures of model confidence


title: What is transiting HD 139139 ?

abstract: The NASA mission Kepler has detected 28 transits with depths and durations hours in the light curve of HD139139 during a 87 days campaign . Their arrival times are erratic. Rappaport et al. (2019) discard ten explanations. It is not clear if the transits are for HD 139139 or for a star B at 3.3 arcsec. New radial velocity variation data give RV for HD139139 in 4 days (F. Bouchy and S. Udry, private communication) excluding a close-in planet orbiting HD 139139, as proposed earlier (Schneider 2019). (However, this explanation is still valid for the star B for which radial velocity data are very poor). Here I thus explore new tentative explanations and their likelihood: 1/ An eccentric transiters belt around HD 139139 2/ Interstellar transiters 3/ Solar System objects.


title: CAN Radar: Sensing Physical Devices in CAN Networks based on Time Domain Reflectometry

abstract: The presence of security vulnerabilities in automotive networks has already been shown by various publications in recent years. Due to the specification of the Controller Area Network (CAN) as a broadcast medium without security mechanisms, attackers are able to read transmitted messages without being noticed and to inject malicious messages. In order to detect potential attackers within a network or software system as early as possible, Intrusion Detection Systems (IDSs) are prevalent. Many approaches for vehicles are based on techniques which are able to detect deviations from specified CAN network behaviour regarding protocol or payload properties. However, it is challenging to detect attackers who secretly connect to CAN networks and do not actively participate in bus traffic. In this paper, we present an approach that is capable of successfully detecting unknown CAN devices and determining the distance (cable length) between the attacker device and our sensing unit based on Time Domain Reflectometry (TDR) technique. We evaluated our approach on a real vehicle network.


title: Three-Dimensional Stability of Current Sheets Supported by Electron Pressure Anisotropy

abstract: The stability of electron current sheets embedded within the reconnection exhaust is studied with a 3D fully kinetic particle-in-cell simulation. The electron current layers studied here form self-consistently in a reconnection regime with a moderate guide field, are supported by electron pressure anisotropy with the pressure component parallel to the magnetic field direction larger than the perpendicular components, and extend well beyond electron kinetic scales. In 3D, in addition to drift instabilities common to nearly all reconnection exhausts, the regime considered also exhibits an electromagnetic instability driven by the electron pressure anisotropy. While the fluctuations modulate the current density on small scales, they do not break apart the general structure of the extended electron current layers. The elongated current sheets should therefore persist long enough to be observed both in space observations and in laboratory experiments.


title: Space-efficient quantum multiplication of polynomials for binary finite fields with sub-quadratic Toffoli gate count

abstract: Multiplication is an essential step in a lot of calculations. In this paper we look at multiplication of 2 binary polynomials of degree at most $n-1$, modulo an irreducible polynomial of degree $n$ with $2n$ input and $n$ output qubits, without ancillary qubits, assuming no errors. With straightforward schoolbook methods this would result in a quadratic number of Toffoli gates and a linear number of CNOT gates. This paper introduces a new algorithm that uses the same space, but by utilizing space-efficient variants of Karatsuba multiplication methods it requires only $O(n^{\log_2(3)})$ Toffoli gates at the cost of a higher CNOT gate count: theoretically up to $O(n^2)$ but in examples the CNOT gate count looks a lot better.


title: A Study of Power Suppressed Contributions in $J/\psi\rightarrow p\bar{p}$ Decay

abstract: The power suppressed amplitude which describes the Pauli ($\sigma^{\mu\nu}$) coupling in the $J/\psi\to p\bar p$ decay is calculated within the effective field theory framework. It is shown that at the leading-order approximation this contribution is factorisable and the overlap with the hadronic final state can be described by collinear matrix elements. The obtained contribution depends on the nucleon light-cone distribution amplitudes of twist-3 and twist-4. This result is used for a qualitative phenomenological analysis of existing data for $J/\psi\to p\bar p$ decay: branching ratio and the angular distribution in the cross section $e^+e^-\to J/\psi\to p\bar p$.


title: ER-index: a referential index for encrypted genomic databases

abstract: Huge DBMSs storing genomic information are being created and engineerized for doing large-scale, comprehensive and in-depth analysis of human beings and their diseases. However, recent regulations like the GDPR require that sensitive data are stored and elaborated thanks to privacy-by-design methods and software. We designed and implemented ER-index, a new full-text index in minute space which was optimized for compressing and encrypting collections of genomic sequences, and for performing on them fast pattern-search queries. Our new index complements the E2FM-index, which was introduced to compress and encrypt collections of nucleotide sequences without relying on a reference sequence. When used on collections of highly similar sequences, the ER-index allows to obtain compression ratios which are an order of magnitude smaller than those achieved with the E2FM-index, but maintaining its very good search performance. Moreover, thanks to the ER-index multi-user and multiple-keys encryption model, a single index can store the sequences related to a population of individuals so that users may perform search operations only on the sequences to which they were granted access. The ER-index C++ source code plus scripts and data to assess the tool performance are available at: https://github.com/EncryptedIndexes/erindex.


title: Geometrical Nonlinearity of Circular Plates and Membranes: an Alternative Method

abstract: We apply the well-established theoretical method developed for geometrical nonlinearities of micro/nano-mechanical clamped beams to circular drums. The calculation is performed under the same hypotheses, the extra difficulty being to analytically describe the (coordinate-dependent) additional stress generated in the structure by the motion. Specifically, the model applies to non-axisymmetric mode shapes. An analytic expression is produced for the Duffing (hardening) nonlinear coefficient, which requires only the knowledge of the mode shape functions to be evaluated. This formulation is simple to handle, and does not rely on complex numerical methods. Moreover, no hypotheses are made on the drive scheme and the nature of the in-plane stress: it is not required to be of electrostatic origin. We confront our predictions with both typical experimental devices and relevant theoretical results from the literature. Generalization of the presented method to Duffing-type mode-coupling should be a straightforward extension of this work. We believe that the presented modeling will contribute to the development of nonlinear physics implemented in 2D micro/nano-mechanical structures.


title: Probing Unification With Precision Higgs Physics

abstract: We propose a novel approach of probing grand unification through precise measurements on the Higgs Yukawa couplings at the LHC. This idea is well motivated by the appearance of effective operators not suppressed by the mass scale of unification $M_{\rm{U}}$ in realistic models of unification with the minimal structure of Yukawa sector. Such operators modify the Higgs Yukawa couplings in correlated patterns at scale $M_{\rm{U}}$ that hold up to higher-order corrections. The coherences reveal a feature that, the deviation of tau Yukawa coupling relative to its standard model value at the weak scale is the largest one among the third-generation Yukawa couplings. This feature, if verified by the future LHC, can serve as a hint of unification.


title: Mode stability for the Teukolsky equation on extremal and subextremal Kerr spacetimes

abstract: We prove that there are no exponentially growing modes nor modes on the real axis for the Teukolsky equation on Kerr black hole spacetimes, both in the extremal and subextremal case. We also give a quantitative refinement of mode stability. As an immediate application, we show that the transmission and reflection coefficients of the scattering problem are bounded, independently of the specific angular momentum of the black hole, in any compact set of real frequencies excluding zero frequency and the superradiant threshold. While in the subextremal setting these results were known previously, the extremal case is more involved and has remained an open problem. Ours are the first results outside axisymmetry and could serve as a preliminary step towards understanding boundedness, scattering and decay properties of general solutions to the Teukolsky equation on extremal Kerr black holes.


title: NuSTAR observations of the unidentified INTEGRAL sources: constraints on the Galactic population of HMXBs

abstract: The NuSTAR Legacy program 'Unidentified INTEGRAL sources' targeted faint hard X-ray sources revealed by INTEGRAL in the Galactic plane in order to provide conclusive identification of their nature and insights on the population of faint hard X-ray sources. The NuSTAR and Swift/XRT observations obtained in 2015-2017 contributed to the successful identification of five persistent sources. Here, we report on the spectral and variability analyses which helped to consolidate the classifications of IGR J10447-6027, IGR J16181-5407 and IGR J20569+4940 as active galactic nuclei, and IGR J17402-3656 as an intermediate polar. An optical spectrum of the blazar IGR J20569+4940 is also presented. Combining these results with successful identifications of other such faint and persistent INTEGRAL sources reported in the literature, we investigate possible implications for the population of persistent high mass X-ray binaries (HMXBs) below the identification completion limit of the INTEGRAL survey. The current trend hints at a deficit of persistent HMXBs below $F_{\rm 17-60keV} = 10^{-11}$ erg cm$^{-2}$ s$^{-1}$, but additional efforts dedicated to classifying faint hard X-ray sources are needed before we can draw solid conclusions.


title: Combinatorial considerations on the invariant measure of a stochastic matrix

abstract: The invariant measure is a fundamental object in the theory of Markov processes. In finite dimensions a Markov process is defined by transition rates of the corresponding stochastic matrix. The Markov tree theorem provides an explicit representation of the invariant measure of a stochastic matrix. In this note, we given a simple and purely combinatorial proof of the Markov tree theorem. In the symmetric case of detailed balance, the statement and the proof simplifies even more.


title: All-at-once versus reduced iterative methods for time dependent inverse problems

abstract: In this paper we investigate all-at-once versus reduced regularization of dynamic inverse problems on finite time intervals $(0,T)$. In doing so, we concentrate on iterative methods and nonlinear problems, since they have already been shown to exhibit considerable differences in their reduced and all-at-once versions, whereas Tikhonov regularization is basically the same in both settings. More precisely, we consider Landweber iteration, the iteratively regularized Gauss-Newton method, and the Landweber-Kaczmarz method, the latter relying on cyclic iteration over a subdivision of the problem into subsequent subintervals of $(0,T)$. Part of the paper is devoted to providing an appropriate function space setting as well as establishing the required differentiability results needed for well-definedness and convergence of the methods under consideration. Based on this, we formulate and compare the above mentioned iterative methods in their all-at-once and their reduced version. Finally, we provide some convergence results in the framework of Hilbert space regularization theory and illustrate the convergence conditions by an example of an inverse source problem for a nonlinear diffusion equation.


title: FLEXI: A high order discontinuous Galerkin framework for hyperbolic-parabolic conservation laws

abstract: High order (HO) schemes are attractive candidates for the numerical solution of multiscale problems occurring in fluid dynamics and related disciplines. Among the HO discretization variants, discontinuous Galerkin schemes offer a collection of advantageous features which have lead to a strong increase in interest in them and related formulations in the last decade. The methods have matured sufficiently to be of practical use for a range of problems, for example in direct numerical and large eddy simulation of turbulence. However, in order to take full advantage of the potential benefits of these methods, all steps in the simulation chain must be designed and executed with HO in mind. Especially in this area, many commercially available closed-source solutions fall short. In this work, we therefor present the FLEXI framework, a HO consistent, open-source simulation tool chain for solving the compressible Navier-Stokes equations in a high performance computing setting. We describe the numerical algorithms and implementation details and give an overview of the features and capabilities of all parts of the framework. Beyond these technical details, we also discuss the important, but often overlooked issues of code stability, reproducibility and user-friendliness. The benefits gained by developing an open-source framework are discussed, with a particular focus on usability for the open-source community. We close with sample applications that demonstrate the wide range of use cases and the expandability of FLEXI and an overview of current and future developments.


title: Assessing and Visualizing Matrix Variate Normality

abstract: A framework for assessing the matrix variate normality of three-way data is developed. The framework comprises a visual method and a goodness of fit test based on the Mahalanobis squared distance (MSD). The MSD of multivariate and matrix variate normal estimators, respectively, are used as an assessment tool for matrix variate normality. Specifically, these are used in the form of a distance-distance (DD) plot as a graphical method for visualizing matrix variate normality. In addition, we employ the popular Kolmogorov-Smirnov goodness of fit test in the context of assessing matrix variate normality for three-way data. Finally, an appropriate simulation study spanning a large range of dimensions and data sizes shows that for various settings, the test proves itself highly robust.


title: Cable Manipulation with a Tactile-Reactive Gripper

abstract: Cables are complex, high dimensional, and dynamic objects. Standard approaches to manipulate them often rely on conservative strategies that involve long series of very slow and incremental deformations, or various mechanical fixtures such as clamps, pins or rings. We are interested in manipulating freely moving cables, in real time, with a pair of robotic grippers, and with no added mechanical constraints. The main contribution of this paper is a perception and control framework that moves in that direction, and uses real-time tactile feedback to accomplish the task of following a dangling cable. The approach relies on a vision-based tactile sensor, GelSight, that estimates the pose of the cable in the grip, and the friction forces during cable sliding. We achieve the behavior by combining two tactile-based controllers: 1) Cable grip controller, where a PD controller combined with a leaky integrator regulates the gripping force to maintain the frictional sliding forces close to a suitable value; and 2) Cable pose controller, where an LQR controller based on a learned linear model of the cable sliding dynamics keeps the cable centered and aligned on the fingertips to prevent the cable from falling from the grip. This behavior is possible by a reactive gripper fitted with GelSight-based high-resolution tactile sensors. The robot can follow one meter of cable in random configurations within 2-3 hand regrasps, adapting to cables of different materials and thicknesses. We demonstrate a robot grasping a headphone cable, sliding the fingers to the jack connector, and inserting it. To the best of our knowledge, this is the first implementation of real-time cable following without the aid of mechanical fixtures.


title: Spectral sparsification of matrix inputs as a preprocessing step for quantum algorithms

abstract: We study the potential utility of classical techniques of spectral sparsification of graphs as a preprocessing step for digital quantum algorithms, in particular, for Hamiltonian simulation. Our results indicate that spectral sparsification of a graph with $n$ nodes through a sampling method, e.g.\ as in \cite{Spielman2011resistances} using effective resistances, gives, with high probability, a locally computable matrix $\tilde H$ with row sparsity at most $\mathcal{O}(\text{poly}\log n)$. For a symmetric matrix $H$ of size $n$ with $m$ non-zero entries, a one-time classical runtime overhead of $\mathcal{O}(m||H||t\log n/\epsilon)$ expended in spectral sparsification is then found to be useful as a way to obtain a sparse matrix $\tilde H$ that can be used to approximate time evolution $e^{itH}$ under the Hamiltonian $H$ to precision $\epsilon$. Once such a sparsifier is obtained, it could be used with a variety of quantum algorithms in the query model that make crucial use of row sparsity. We focus on the case of efficient quantum algorithms for sparse Hamiltonian simulation, since Hamiltonian simulation underlies, as a key subroutine, several quantum algorithms, including quantum phase estimation and recent ones for linear algebra. Finally, we also give two simple quantum algorithms to estimate the row sparsity of an input matrix, which achieve a query complexity of $\mathcal{O}(n^{3/2})$ as opposed to $\mathcal{O}(n^2)$ that would be required by any classical algorithm for the task.


title: On The Local Sum Conjecture in Two Dimensions

abstract: In this paper we give an elementary proof of the local sum conjecture in two dimensions. In a remarkable paper [CMN, arXiv:1810.11340], this conjecture has been established in all dimensions using sophisticated, powerful techniques from a research area blending algebraic geometry with ideas from logic. The purpose of this paper is to give an elementary proof of this conjecture which will be accessbile to a broad readership.


title: Provenance tracking in the LHCb software

abstract: Even though computational reproducibility is widely accepted as necessary for research validation and reuse, it is often not considered during the research process. This is because reproducibility tools are typically stand-alone and require additional training to be employed. In this article, we present a solution to foster reproducibility, which is integrated within existing scientific software that is actively used in the LHCb collaboration. Our provenance tracking service captures metadata of a dataset, which is then saved inside the output data file on the disk. The captured information allows a complete understanding of how the file was produced and enables a user to reproduce the dataset, even when the original input code (that was used to initially produce the dataset) is altered or lost. This article describes the implementation of the service and gives examples of its application.


title: Possible lightest $\Xi$ Hypernucleus with Modern $\Xi N$ Interactions

abstract: Experimental evidence exists that the $\Xi$-nucleus interaction is attractive. We search for $NN\Xi$ and $NNN\Xi$ bound systems on the basis of the AV8 $NN$ potential combined with either a phenomenological Nijmegen $\Xi N$ potential or a first principles HAL QCD $\Xi N$ potential. The binding energies of the three-body and four-body systems (below the $d+\Xi$ and $^3{\rm H}$/$^3{\rm He}+\Xi$ thresholds, respectively) are calculated by a high precision variational approach, the Gaussian Expansion Method. Although the two $\Xi N$ potentials have significantly different isospin ($T$) and spin ($S$) dependence, the $NNN\Xi$ system with quantum numbers $(T=0, J^{\pi}=1^+$) appears to be bound (one deep for Nijmegen and one shallow for HAL QCD) below the $^3{\rm H}$/$^3{\rm He}+\Xi$ threshold. Experimental implications for such a state are discussed.


title: Nematic alignment of self-propelled particles in the macroscopic regime

abstract: Starting from a particle model describing self-propelled particles interacting through nematic alignment, we derive a macroscopic model for the particle density and mean direction of motion. We first propose a mean-field kinetic model of the particle dynamics. After diffusive rescaling of the kinetic equation, we formally show that the distribution function converges to an equilibrium distribution in particle direction, whose local density and mean direction satisfies a cross-diffusion system. We show that the system is consistent with symmetries typical of a nematic material. The derivation is carried over by means of a Hilbert expansion. It requires the inversion of the linearized collision operator for which we show that the generalized collision invariants, a concept introduced to overcome the lack of momentum conservation of the system, plays a central role. This cross diffusion system poses many new challenging questions.


title: Nonparametric principal subspace regression

abstract: In scientific applications, multivariate observations often come in tandem with temporal or spatial covariates, with which the underlying signals vary smoothly. The standard approaches such as principal component analysis and factor analysis neglect the smoothness of the data, while multivariate linear or nonparametric regression fail to leverage the correlation information among multivariate response variables. We propose a novel approach named nonparametric principal subspace regression to overcome these issues. By decoupling the model discrepancy, a simple and general two-step framework is introduced, which leaves much flexibility in choice of model fitting. We establish theoretical property of the general framework, and offer implementation procedures that fulfill requirements and enjoy the theoretical guarantee. We demonstrate the favorable finite-sample performance of the proposed method through simulations and a real data application from an electroencephalogram study.


title: Characterization of the equality of weak efficiency and efficiency on convex free disposal hulls

abstract: In solving a multi-objective optimization problem by scalarization techniques, solutions to a scalarized problem are, in general, weakly efficient rather than efficient to the original problem. Thus, it is crucial to understand what problem ensures that all weakly efficient solutions are efficient. In this paper, we give a characterization of the equality of the weakly efficient set and the efficient set, provided that the free disposal hull of the domain is convex. By using this characterization, we obtain various mathematical applications. As a practical application, we show that all weakly efficient solutions to a multi-objective LASSO with mild modification are efficient.


title: A Differential Approximation Model For Passive Scalar Turbulence

abstract: Two dimensional passive scalar turbulence is studied by means of a k-space diffusion model based on a third order differential approximation. This simple description of local nonlinear interactions in Fourier space is shown to present a general expression, in line with previous seminal works, and appears to be suitable for various 2D turbulence problems. Steady state solutions for the spectral energy density of the flow is shown to recover the Kraichnan-Kolmogorov phenomenology of the dual cascade, while various passive scalar spectra, such as Batchelor or Obukhov-Corssin spectra are recovered as steady state solutions of the spectral energy density of the passive scalar. These analytical results are then corroborated by numerical solutions of the time evolving problem with energy and passive scalar injection and dissipation on a logarithmic wavenumber space grid over a large range of scales. The particular power law spectra that is obtained is found to depend mainly on the location of the kinetic and passive scalar energy injections.


title: Anomalous Interlayer Exciton Diffusion in Twist-Angle-Dependent Moir\'{e} Potentials of WS$_2$-WSe$_2$ Heterobilayers

abstract: The nanoscale periodic potentials introduced by moir\'{e} patterns in semiconducting van der Waals (vdW) heterostructures provide a new platform for designing exciton superlattices. To realize these applications, a thorough understanding of the localization and delocalization of interlayer excitons in the moir\'{e} potentials is necessary. Here, we investigated interlayer exciton dynamics and transport modulated by the moir\'{e} potentials in WS$_2$-WSe$_2$ heterobilayers in time, space, and momentum domains using transient absorption microscopy combined with first-principles calculations. Experimental results verified the theoretical prediction of energetically favorable K-Q interlayer excitons and unraveled exciton-population dynamics that was controlled by the twist-angle-dependent energy difference between the K-Q and K-K excitons. Spatially- and temporally-resolved exciton-population imaging directly visualizes exciton localization by twist-angle-dependent moir\'{e} potentials of ~100 meV. Exciton transport deviates significantly from normal diffusion due to the interplay between the moir\'{e} potentials and strong many-body interactions, leading to exciton-density- and twist-angle-dependent diffusion length. These results have important implications for designing vdW heterostructures for exciton and spin transport as well as for quantum communication applications.


title: Electric double layers with surface charge modulations: Novel exact Poisson-Boltzmann solutions

abstract: Poisson-Boltzmann theory is the cornerstone for soft matter electrostatics. We provide novel exact analytical solutions to this non-linear mean-field approach, for the diffuse layer of ions in the vicinity of a planar or a cylindrical macroion. While previously known solution are for homogeneously charged objects, the cases worked out exhibit a modulated surface charge --or equivalently surface potential-- on the macroion (wall) surface. In addition to asymptotic features at large distances from the wall, attention is paid to the fate of the contact theorem, relating the contact density of ions to the local wall charge density. For salt-free systems (counterions only), we make use of results pertaining to the two-dimensional Liouville equation, supplemented by an inverse approach. When salt is present, we invoke the exact two-soliton solution to the 2D sinh-Gordon equation. This leads to inhomogeneous charge patterns, that are either localized or periodic in space. Without salt, the electrostatic signature of a charge pattern on the macroion fades exponentially with distance for a planar macroion, while it decays as an inverse power-law for a cylindrical macroion. With salt, our study is limited to the planar geometry, and reveals that pattern screening is exponential.


title: Stationary distributions via decomposition of stochastic reaction networks

abstract: We examine reaction networks (CRNs) through their associated continuous-time Markov processes. Studying the dynamics of such networks is in general hard, both analytically and by simulation. In particular, stationary distributions of stochastic reaction networks are only known in some cases. We analyze class properties of the underlying continuous-time Markov chain of CRNs under the operation of join and examine conditions such that the form of the stationary distributions of a CRN is derived from the parts of the decomposed CRNs. The conditions can be easily checked in examples and allow recursive application. The theory developed enables sequential decomposition of the Markov processes and calculations of stationary distributions. Since the class of processes expressible through such networks is big and only few assumptions are made, the principle also applies to other stochastic models. We give examples of interest from CRN theory to highlight the decomposition.


title: Taylor spectrum approach to Brownian-type operators with quasinormal entry

abstract: In this paper, we introduce operators that are represented by upper triangular $2\times 2$ block matrices whose entries satisfy some algebraic constraints. We call them Brownian-type operators of class $\mathcal Q,$ briefly operators of class $\mathcal Q.$ These operators emerged from the study of Brownian isometries performed by Agler and Stankus via detailed analysis of the time shift operator of the modified Brownian motion process. It turns out that the class $\mathcal Q$ is closely related to the Cauchy dual subnormality problem which asks whether the Cauchy dual of a completely hyperexpansive operator is subnormal. Since the class $\mathcal Q$ is closed under the operation of taking the Cauchy dual, the problem itself becomes a part of a more general question of investigating subnormality in this class. This issue, along with the analysis of nonstandard moment problems, covers a large part of the paper. Using the Taylor spectrum technique culminates in a full characterization of subnormal operators of class $\mathcal Q.$ As a consequence, we solve the Cauchy dual subnormality problem for expansive operators of class $\mathcal Q$ in the affirmative, showing that the original problem can surprisingly be extended to a class of operators that are far from being completely hyperexpansive. The Taylor spectrum approach turns out to be fruitful enough to allow us to characterize other classes of operators including $m$-isometries. We also study linear operator pencils associated with operators of class $\mathcal Q$ proving that the corresponding regions of subnormality are closed intervals with explicitly described endpoints.


title: Two-Dimensional Optomechanical Crystal Cavity with High Quantum Cooperativity

abstract: Optomechanical systems offer new opportunities in quantum information processing and quantum sensing. Many solid-state quantum devices operate at millikelvin temperatures -- however, it has proven challenging to operate nanoscale optomechanical devices at these ultralow temperatures due to their limited thermal conductance and parasitic optical absorption. Here, we demonstrate a two-dimensional optomechanical crystal resonator capable of achieving large cooperativity $C$ and small effective bath occupancy $n_b$, resulting in a quantum cooperativity $C_{\text{eff}}\equiv C/n_b \approx 1.3 > 1$ under continuous-wave optical driving. This is realized using a two-dimensional phononic bandgap structure to host the optomechanical cavity, simultaneously isolating the acoustic mode of interest in the bandgap while allowing heat to be removed by phonon modes outside of the bandgap. This achievement paves the way for a variety of applications requiring quantum-coherent optomechanical interactions, such as transducers capable of bi-directional conversion of quantum states between microwave frequency superconducting quantum circuits and optical photons in a fiber optic network.


title: Field-based Coordination with the Share Operator

abstract: Field-based coordination has been proposed as a model for coordinating collective adaptive systems, promoting a view of distributed computations as functions manipulating data structures spread over space and evolving over time, called computational fields. The field calculus is a formal foundation for field computations, providing specific constructs for evolution (time) and neighbor interaction (space), which are handled by separate operators (called rep and nbr, respectively). This approach, however, intrinsically limits the speed of information propagation that can be achieved by their combined use. In this paper, we propose a new field-based coordination operator called share, which captures the space-time nature of field computations in a single operator that declaratively achieves: (i) observation of neighbors' values; (ii) reduction to a single local value; and (iii) update and converse sharing to neighbors of a local variable. We show that for an important class of self-stabilising computations, share can replace all occurrences of rep and nbr constructs. In addition to conceptual economy, use of the share operator also allows many prior field calculus algorithms to be greatly accelerated, which we validate empirically with simulations of frequently used network propagation and collection algorithms.


title: The asymptotic spectrum of the Hessian of DNN throughout training

abstract: The dynamics of DNNs during gradient descent is described by the so-called Neural Tangent Kernel (NTK). In this article, we show that the NTK allows one to gain precise insight into the Hessian of the cost of DNNs. When the NTK is fixed during training, we obtain a full characterization of the asymptotics of the spectrum of the Hessian, at initialization and during training. In the so-called mean-field limit, where the NTK is not fixed during training, we describe the first two moments of the Hessian at initialization.


title: Reinforcement Learning with Structured Hierarchical Grammar Representations of Actions

abstract: From a young age humans learn to use grammatical principles to hierarchically combine words into sentences. Action grammars is the parallel idea, that there is an underlying set of rules (a "grammar") that govern how we hierarchically combine actions to form new, more complex actions. We introduce the Action Grammar Reinforcement Learning (AG-RL) framework which leverages the concept of action grammars to consistently improve the sample efficiency of Reinforcement Learning agents. AG-RL works by using a grammar inference algorithm to infer the "action grammar" of an agent midway through training. The agent's action space is then augmented with macro-actions identified by the grammar. We apply this framework to Double Deep Q-Learning (AG-DDQN) and a discrete action version of Soft Actor-Critic (AG-SAC) and find that it improves performance in 8 out of 8 tested Atari games (median +31%, max +668%) and 19 out of 20 tested Atari games (median +96%, maximum +3,756%) respectively without substantive hyperparameter tuning. We also show that AG-SAC beats the model-free state-of-the-art for sample efficiency in 17 out of the 20 tested Atari games (median +62%, maximum +13,140%), again without substantive hyperparameter tuning.


title: Heap and Ternary Self-Distributive Cohomology

abstract: Heaps are para-associative ternary operations bijectively exemplified by groups via the operation $(x,y,z) \mapsto x y^{-1} z$. They are also ternary self-distributive, and have a diagrammatic interpretation in terms of framed links. Motivated by these properties, we define para-associative and heap cohomology theories and also a ternary self-distributive cohomology theory with abelian heap coefficients. We show that one of the heap cohomologies is related to group cohomology via a long exact sequence. Moreover we construct maps between second cohomology groups of normalized group cohomology and heap cohomology, and show that the latter injects into the ternary self-distributive second cohomology group. We proceed to study heap objects in symmetric monoidal categories providing a characterization of pointed heaps as involutory Hopf monoids in the given category. Finally we prove that heap objects are also "categorically" self-distributive in an appropriate sense.


title: Topological Origin of Non-Hermitian Skin Effects

abstract: A unique feature of non-Hermitian systems is the skin effect, which is the extreme sensitivity to the boundary conditions. Here, we reveal that the skin effect originates from intrinsic non-Hermitian topology. Such a topological origin not merely explains the universal feature of the known skin effect, but also leads to new types of the skin effects---symmetry-protected skin effects. In particular, we discover the $\mathbb{Z}_{2}$ skin effect protected by time-reversal symmetry. On the basis of topological classification, we also discuss possible other skin effects in arbitrary dimensions. Our work provides a unified understanding about the bulk-boundary correspondence and the skin effects in non-Hermitian systems.

categories: cond-mat.mes-hall cond-mat.str-el math-ph math.MP physics.optics quant-ph


title: Direct Visualization of Electric Field induced Structural Dynamics in Monolayer Transition Metal Dichalcogenides

abstract: Layered transition metal dichalcogenides (TMDs) offer many attractive features for next-generation low-dimensional device geometries. Due to the practical and fabrication challenges related to in situ methods, the atomistic dynamics that give rise to realizable macroscopic device properties are often unclear. In this study, in situ transmission electron microscopy techniques are utilized in order to understand the structural dynamics at play, especially at interfaces and defects, in the prototypical film of monolayer MoS2 under electrical bias. Through our sample fabrication process, we clearly identify the presence of mass transport in the presence of a lateral electric field. In particular, we observe that the voids present at grain boundaries combine to induce structural deformation. The electric field mediates a net vacancy flux from the grain boundary interior to the exposed surface edge sites that leaves molybdenum clusters in its wake. Following the initial biasing cycles, however, the mass flow is largely diminished, and the resultant structure remains stable over repeated biasing. We believe insights from this work can help explain observations of non-uniform heating and preferential oxidation at grain boundary sites in these materials.


title: Heat diffusion related damping process in highly precise coarse-grained model for SWCNT's nonlinearity

abstract: Second sound and heat diffusion in single-walled carbon nanotubes (SWCNT) is well-known phenomena which is related to the high thermal conductivity of this material. In this paper, we have shown that the heat diffusion along the tube axis affects the macroscopic motion of SWCNT and adapting this phenomena to coarse-grained model can improve the precision of the coarse-grained molecular dynamics (CGMD) exceptionally. The nonlinear macroscopic motion of SWCNT in the free thermal vibration condition in adiabatic environment is demonstrated in the most simplified version of CG modeling as maintaining finite temperature and total energy with suggested dissipation process derived from internal heat diffusion. The internal heat diffusion related to the cross correlated momentum from different potential energy functions is considered, and it can reproduce the nonlinear dynamic nature of SWCNTs without external thermostatting in CG model. Memory effect and thermostat with random noise distribution are not included, and the effect of heat diffusion on memory effect is reviewed through Mori-Zwanzig formalism. This diffusion shows perfect syncronization of the motion between that of CGMD and MD simulation, which is started with initial conditions from the molecular dynamics (MD) simulation. The heat diffusion related to this process has shown the same dispersive characteristics to second wave in SWCNT. This replication with good precision indicates that the internal heat diffusion process is the essential cause of the nonlinearity of the tube. The nonlinear dynamic characteristics from the various scale of simple beads systems are examined with expanding its time step and node length simultaneously.

categories: cond-mat.soft cond-mat.mes-hall nlin.CD physics.comp-ph


title: Axial-vector exchange contribution to the hadronic light-by-light piece of the muon anomalous magnetic moment

abstract: In this work we study the axial contributions to the hadronic light-by-light piece of the muon anomalous magnetic moment using the framework of resonance chiral theory. As a result, we obtain $a_{\mu}^{\textrm{HLbL};A} = \left(0.8^{+3.5}_{-0.1}\right)\cdot 10^{-11}$, that might suggest a smaller value than most recent calculations, underlining the need of future work along this direction. In particular, we find that our results depend critically on the asymptotic behavior of the form factors, and as such, emphasizes the relevance of future experiments for large photon virtualities. In addition, we present general results regarding the involved axial form factors description, comprehensively examining (and relating) the current approaches, that shall be of general interest.


title: Anisotropic particles in two-dimensional convective turbulence

abstract: The orientational dynamics of inertialess anisotropic particles transported by two-dimensional convective turbulent flows display a coexistence of regular and chaotic features. We numerically demonstrate that very elongated particles (rods) align preferentially with the direction of the fluid flow, i.e., horizontally close to the isothermal walls and dominantly vertically in the bulk. This behaviour is due to the the presence of a persistent large scale circulation flow structure, which induces strong shear at wall boundaries and in up/down-welling regions. The near-wall horizontal alignment of rods persists at increasing the Rayleigh number, while the vertical orientation in the bulk is progressively weakened by the corresponding increase of turbulence intensity. Furthermore, we show that very elongated particles are nearly orthogonal to the orientation of the temperature gradient, an alignment independent of the system dimensionality and which becomes exact only in the limit of infinite Prandtl number. Tumbling rates are extremely vigorous adjacent to the walls, where particles roughly perform Jeffery orbits. This implies that the root-mean-square near-wall tumbling rates for spheres are much stronger than for rods, up to $\mathcal{O}(10)$ times at $Ra\simeq 10^9$. In the turbulent bulk the situation reverses and rods tumble slightly faster than isotropic particles, in agreement with earlier observations in two-dimensional turbulence.


title: The transverse nucleon single-spin asymmetry for the semi-inclusive production of photons in lepton-nucleon scattering

abstract: We study the semi-inclusive production of real, high-$p_T$, isolated photons in unpolarized and polarized lepton-proton collisions, $\ell p\to \ell\gamma X$. In particular we analyze the transverse nucleon single-spin asymmetry within the collinear twist-3 formalism in perturbative QCD to leading order (LO) accuracy. We find that this spin asymmetry is generated by twist-3 dynamical quark-gluon-quark ($qgq$) correlations in the nucleon through the so-called soft-fermion and hard pole contributions. Hence, this process unprecedentedly allows for a point-by-point scan of the support of the dynamical $qgq$ twist-3 matrix elements $F_{FT}(x,x^\prime)$ and $G_{FT}(x,x^\prime)$ in lepton-nucleon scattering experiments.


title: Introduction to Concentration Inequalities

abstract: In this report, we aim to exemplify concentration inequalities and provide easy to understand proofs for it. Our focus is on the inequalities which are helpful in the design and analysis of machine learning algorithms.


abstract: In his paper Almost-Primes Represented by Quadratic Polynomials, Iwaniec proved that the polynomial n^2 + 1 takes on values with at most two prime factors (counted with multiplicity) infinitely often. He states that "in order to avoid technical complications, we shall restrict our proof to the polynomial n^2 + 1.". In this exposition, we follow Iwaniec's proof and show that for any irreducible quadratic polynomial G(n) (satisfying some obviously necessary hypotheses), G(n) has at most two prime factors for infinitely many values of n.


title: Measuring relative abundances in the solar corona with optimised linear combinations of spectral lines

abstract: Context. Elemental abundances in some coronal structures differ significantly from photospheric abundances, with a dependence on the first ionization potential (FIP) of the element. Measuring these FIP-dependent abundance biases is important for coronal and heliospheric physics. Aims. We aim to build a method for optimal determination of FIP biases in the corona from spectroscopic observations in a way that is in practice independent from Differential Emission Measure (DEM) inversions. Methods. We optimised linear combinations of spectroscopic lines of low-FIP and high-FIP elements so that the ratio of the corresponding radiances yields the relative FIP bias with good accuracy for any DEM in a small set of typical DEMs. Results. These optimised linear combinations of lines allow retrieval of a test FIP bias map with good accuracy for all DEMs in the map. The results also compare well with a FIP bias map obtained from observations using a DEM-dependent method. Conclusions. The method provides a convenient, fast, and accurate way of computing relative FIP bias maps. It can be used to optimise the use of existing observations and the design of new observations and instruments.


title: Asymptotic Analysis Of Determinant Of Discrete Laplacian

abstract: In this paper, we study the relation between the partition function of the free scalar field theory on hypercubes with boundary conditions and asymptotics of discrete partition functions on a sequence of "lattices" which approximate the hypercube as the mesh approaches to zero. More precisely, we show that the logarithm of the zeta regularized determinant of Laplacian on the hypercube with Dirichlet boundary condition appears as the constant term in the asymptotic expansion of the log-determinant of the discrete Laplacian up to an explicitly computable constant. We also investigate similar problems for the massive Laplacian on tori.


title: A note on $\mu$-stabilizers in ACVF

abstract: We study $\mu$-stabilizers for groups definable in ACVF in the valued field sort. We prove that $\mathrm{Stab}^\mu(p)$ is an infinite unbounded definable subgroup of $G$ when $p$ is standard and unbounded. In the particular case when $G$ is linear algebraic, we show that $\mathrm{Stab}^\mu(p)$ is a solvable algebraic subgroup of $G$, with $\mathrm{dim}(\mathrm{Stab}^\mu(p))=\mathrm{dim}(p)$ when $p$ is $\mu$-reduced and unbounded.


title: Low-redshift 21cm Cosmology in Canada

abstract: Line-intensity mapping of the 21cm line is a powerful probe of large scale structure at z<6, tracing large-scale structure via neutral hydrogen content that is found within galaxies. In principle, it enables cost-efficient surveys of the matter distribution up to z~6, unlocking orders of magnitude more modes for observational cosmology. Canada has been a traditional leader in this field, having led the first detections of the cosmological 21cm signal via cross-correlations with optical galaxy surveys and having constructed the Canadian Hydrogen Intensity Mapping Experiment (CHIME). The field is now entering a new era where data is abundant, allowing studies in how to overcome systematics to be tackled in an empirical, head-on fashion. In the next few years, this will produce the first detection of the 21cm auto power spectrum, which will pave the way towards a large suite of scientific possibilities. These potentially include precision measurements on the dark energy equation of state and other LCDM parameters, constraints on how HI mass traces dark matter, a detection of neutrino effects on large-scale structure, and the use of 21cm lensing to further constrain cosmology. To turn these promising directions into reality, we recommend a sustained program of investment in 21cm cosmology, starting with funding for the Canadian Hydrogen Observatory and Radio transient Detector (CHORD), followed by small-scale development efforts targeting next-generation hardware and sustained support for theory and technical staff support. Additionally, Canada should invest in complementary line-intensity mapping efforts (such as with CO or [CII] lines) and maintain participation in next-generation international efforts such as the Packed Ultra-wideband Mapping Array (PUMA) and the Square Kilometre Array (SKA).


title: MADHAT: Model-Agnostic Dark Halo Analysis Tool

abstract: We present the Model-Agnostic Dark Halo Analysis Tool (MADHAT), a numerical tool which implements a Fermi-LAT data-driven, model-independent analysis of gamma-ray emission from dwarf satellite galaxies and dwarf galaxy candidates due to dark matter annihilation, dark matter decay, or other nonstandard or unknown astrophysics. This tool efficiently provides statistical upper bounds on the number of observed photons in excess of the number expected, based on empirical determinations of foregrounds and backgrounds, using a stacked analysis of any selected set of dwarf targets. It also calculates the resulting bounds on the properties of dark matter under any assumptions the user makes regarding dark sector particle physics or astrophysics. As an application, we determine new bounds on Sommerfeld-enhanced dark matter annihilation in a set of eight dwarfs. MADHAT v1.0 includes 58 dwarfs and dwarf candidate targets, and we discuss future planned developments. MADHAT is available and will be maintained at https://github.com/MADHATdm


title: Phase and amplitude asymmetry in the quasi-biennial oscillation of solar Halpha flare activity

abstract: Quasi-biennial oscillation (QBO) of solar magnetic activities is intrinsic to dynamo mechanism, but still far from fully understood. In this work, the phase and amplitude asymmetry of solar QBO of Halpha flare activity in the northern and southern hemispheres is studied by the ensemble empirical mode decomposition, the cross-correlation analysis, and the wavelet transform technique. The following results are found: (1) solar QBO of Halpha flare index in the two hemispheres has a complicated phase relationship, but does not show any systematic regularity; (2) the solar cycle mode of solar Halpha flare index in the northern hemisphere generally leads that in the southern one by 9 months for the time interval from 1966 to 2014. The possible origin of these results is discussed.


title: Chemical abundances along the 1G sequence of the chromosome maps: The Globular Cluster NGC 3201

abstract: The Hubble Space Telescope (HST) UV Legacy Survey of Galactic Globular Clusters (GCs) has investigated multiple stellar populations by means of the "chromosome map" (ChM) diagnostic tool that maximises the separation between stars with different chemical composition. One of the most challenging features revealed by ChMs analysis is the apparent inhomogeneity among stars belonging to the first population, a phenomenon largely attributed to He variations. However, this explanation is not supported by the uniformity in p-capture elements of these stars. The HST survey has revealed that the GC NGC 3201 shows an exceptionally wide coverage in the Delta(F275W,F814W) parameter of the ChM. We present a chemical abundance analysis of 24 elements in 18 giants belonging to the first population of this GC, and having a wide range in Delta(F275W,F814W). As far as the p-capture elements are concerned, the chemical abundances are typical of 1G stars, as expected from the location of our targets in the ChM. Based on radial velocities and chemical abundances arguments, we find that the three stars with the lowest Delta(F275W,F814W) values are binary candidates. This suggests that, at least those stars could be explained with binarity. These results are consistent with evidence inferred from multi-band photometry that evolved blue stragglers populate the bluest part of the 1G sequence in the ChM. The remaining 15 spectroscopic targets show a small range in the overall metallicity by ~0.10 dex, with stars at higher Delta(F275W,F814W) values having higher absolute abundances. We suggest that a small variation in metals and binarity govern the color spread of the 1G in the ChM, and that evolved blue stragglers contribute to the bluest tail of the 1G sequence.


title: Parallel Iterative Edit Models for Local Sequence Transduction

abstract: We present a Parallel Iterative Edit (PIE) model for the problem of local sequence transduction arising in tasks like Grammatical error correction (GEC). Recent approaches are based on the popular encoder-decoder (ED) model for sequence to sequence learning. The ED model auto-regressively captures full dependency among output tokens but is slow due to sequential decoding. The PIE model does parallel decoding, giving up the advantage of modelling full dependency in the output, yet it achieves accuracy competitive with the ED model for four reasons: 1.~predicting edits instead of tokens, 2.~labeling sequences instead of generating sequences, 3.~iteratively refining predictions to capture dependencies, and 4.~factorizing logits over edits and their token argument to harness pre-trained language models like BERT. Experiments on tasks spanning GEC, OCR correction and spell correction demonstrate that the PIE model is an accurate and significantly faster alternative for local sequence transduction.


title: Challenges and Opportunities of Near-Term Quantum Computing Systems

abstract: The concept of quantum computing has inspired a whole new generation of scientists, including physicists, engineers, and computer scientists, to fundamentally change the landscape of information technology. With experimental demonstrations stretching back more than two decades, the quantum computing community has achieved a major milestone over the past few years: the ability to build systems that are stretching the limits of what can be classically simulated, and which enable cloud-based research for a wide range of scientists, thus increasing the pool of talent exploring early quantum systems. While such noisy near-term quantum computing systems fall far short of the requirements for fault-tolerant systems, they provide unique testbeds for exploring the opportunities for quantum applications. Here we highlight the facets associated with these systems, including quantum software, cloud access, benchmarking quantum systems, error correction and mitigation in such systems, and understanding the complexity of quantum circuits and how early quantum applications can run on near term quantum computers.


title: Science with the Murchison Widefield Array: Phase I Results and Phase II Opportunities

abstract: The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.


title: Complications in the interpretation of the charge asymmetry dependent $\pi$ flow for the chiral magnetic wave

abstract: The charge asymmetry ($A_{\rm ch}$) dependence of the $\pi^{-}$ and $\pi^{+}$ elliptic flow difference, $\Delta v_{2}(A_{\rm ch})$, has been regarded as a sensitive observable for the possible chiral magnetic wave (CMW) in relativistic heavy ion collisions. In this work, we first demonstrate that, due to non-flow backgrounds, the flow measurements by the Q-cumulant method using all charged particles as reference introduce a trivial linear term to $\Delta v_{2}(A_{\rm ch})$. The trivial slope can be negative in the triangle flow difference $\Delta v_{3}(A_{\rm ch})$ if the non-flow is dominated by back-to-back pairs. After eliminating the trivial term, we find that the non-flow between like-sign pairs gives rise to an additional positive slope to $\Delta v_{2}(A_{\rm ch})$ because of the larger dilution effect to $\pi^{+}$ ($\pi^{-}$) at positive (negative) $A_{\rm ch}$. We further find that the competition between different $\pi$ sources can introduce another non-trivial linear-$A_{\rm ch}$ term due to their different multiplicity fluctuations and anisotropic flows. We then study the effect of neutral cluster (resonance) decays as a mechanism for local charge conservation on the slope parameter of $\Delta v_{2}(A_{\rm ch})$. We find that the slope parameter is sensitive to the kinematics of those neutral clusters. Light resonances give positive slopes while heavy resonances give negative slopes. Local charge conservation from continuum cluster mass distribution can give a positive slope parameter comparable to experimental data. Our studies indicate that many non-CMW physics mechanisms can give rise to a $A_{\rm ch}$-dependent $\Delta v_{2}(A_{\rm ch})$ and the interpretation of $\Delta v_{2}(A_{\rm ch})$ in terms of the CMW is delicate.


title: Global well-posedness of the 4-d energy-critical stochastic nonlinear Schr\"{o}dinger equations with non-vanishing boundary condition

abstract: We consider the energy-critical stochastic cubic nonlinear Schr\"odinger equation on $\mathbb R^4$ with additive noise, and with the non-vanishing boundary conditions at spatial infinity. By viewing this equation as a perturbation to the energy-critical cubic nonlinear Schr\"odinger equation on $\mathbb R^4$, we prove global well-posedness in the energy space. Moreover, we establish unconditional uniqueness of solutions in the energy space.


title: Fast Fitting of Reflectivity Data of Growing Thin Films Using Neural Networks

abstract: X-ray reflectivity (XRR) is a powerful and popular scattering technique that can give valuable insight into the growth behavior of thin films. In this study, we show how a simple artificial neural network model can be used to predict the thickness, roughness and density of thin films of different organic semiconductors (diindenoperylene, copper(II) phthalocyanine and $\alpha$-sexithiophene) on silica from their XRR data with millisecond computation time and with minimal user input or a priori knowledge. For a large experimental dataset of 372 XRR curves, we show that a simple fully connected model can already provide good predictions with a mean absolute percentage error of 8-18 % when compared to the results obtained by a genetic least mean squares fit using the classical Parratt formalism. Furthermore, current drawbacks and prospects for improvement are discussed.


title: Investigating saturation effects in ultraperipheral collisions at the LHC with the color dipole model

abstract: We investigate saturation effects in $ep$ scattering as well as in ultraperipheral $p$A and AA collisions at small $x$ with four variants of the impact parameter dependent color dipole model: with and without gluon saturation and with and without a novel mechanism that suppresses unphysical dipole radii above the confinement scale, a problem not addressed by most implementations. We show that $ep$ scattering at HERA can be very well described by any of the four variants. When going from $ep$ to $e$A scattering, saturation effects are expected to increase as $\sim$A$^{1/3}$. In lieu of an electron-ion collider, we confront the different versions of the dipole model with data recorded in ultraperipheral collisions at the LHC in order to estimate the sensitivity of the data to gluon saturation in the target nuclei. We find that ultraperipheral PbPb collisions indicate strong saturation effects while $p$Pb collisions turn out to not have any discriminating power to distinguish saturation from non-saturation scenarios.


title: Deep Learning for mmWave Beam and Blockage Prediction Using Sub-6GHz Channels

abstract: Predicting the millimeter wave (mmWave) beams and blockages using sub-6GHz channels has the potential of enabling mobility and reliability in scalable mmWave systems. These gains attracted increasing interest in the last few years. Prior work, however, has focused on extracting spatial channel characteristics at the sub-6GHz band first and then use them to reduce the mmWave beam training overhead. This approach has a number of limitations: (i) It still requires a beam search at mmWave, (ii) its performance is sensitive to the error associated with extracting the sub-6GHz channel characteristics, and (iii) it does not normally account for the different dielectric properties at the different bands. In this paper, we first prove that under certain conditions, there exist mapping functions that can predict the optimal mmWave beam and correct blockage status directly from the sub-6GHz channel, which overcome the limitations in prior work. These mapping functions, however, are hard to characterize analytically which motivates exploiting deep neural network models to learn them. For that, we prove that a large enough neural network can use the sub-6GHz channel to directly predict the optimal mmWave beam and the correct blockage status with success probabilities that can be made arbitrarily close to one. Then, we develop an efficient deep learning model and empirically evaluate its beam/blockage prediction performance using the publicly available dataset DeepMIMO. The results show that the proposed solution can predict the mmWave blockages with more than 90$\%$ success probability. Further, these results confirm the capability of the proposed deep learning model in predicting the optimal mmWave beams and approaching the optimal data rates, that assume perfect channel knowledge, while requiring no beam training overhead...


title: Combined force-torque spectroscopy of proteins by means of multiscale molecular simulation

abstract: Assessing the structural properties of large proteins is important to gain an understanding of their function in, e.g., biological systems or biomedical applications. We propose a method to examine the mechanical properties of proteins subject to applied forces by means of multiscale simulation. We consider both stretching and torsional forces, which can be applied independently of each other. We apply torsional forces to a coarse-grained continuum model of the antibody protein immunoglobulin G (IgG) using Fluctuating Finite Element Analysis and identify the area of strongest deformation. This region is essential to the torsional properties of the molecule as a whole, as it represents the softest, most deformable domain. We subject this part of the molecule to torques and stretching forces on an atomistic level, using molecular dynamics simulations, in order to investigate its torsional properties. We calculate the torsional resistance as a function of the rotation of the domain, while subjecting it to various stretching forces. We learn how these obtained torsion profiles evolve with increasing stretching force and show that they exhibit torsion stiffening, which is in qualitative agreement with experimental findings. We argue that combining the torsion profiles for various stretching forces effectively results in a combined force-torque spectroscopy analysis, which may serve as a mechanical signature for the examined molecule.


title: Correlations between Word Vector Sets

abstract: Similarity measures based purely on word embeddings are comfortably competing with much more sophisticated deep learning and expert-engineered systems on unsupervised semantic textual similarity (STS) tasks. In contrast to commonly used geometric approaches, we treat a single word embedding as e.g. 300 observations from a scalar random variable. Using this paradigm, we first illustrate that similarities derived from elementary pooling operations and classic correlation coefficients yield excellent results on standard STS benchmarks, outperforming many recently proposed methods while being much faster and trivial to implement. Next, we demonstrate how to avoid pooling operations altogether and compare sets of word embeddings directly via correlation operators between reproducing kernel Hilbert spaces. Just like cosine similarity is used to compare individual word vectors, we introduce a novel application of the centered kernel alignment (CKA) as a natural generalisation of squared cosine similarity for sets of word vectors. Likewise, CKA is very easy to implement and enjoys very strong empirical results.


title: Families of Geometries, Real Algebras, and Transitions

abstract: This thesis details the results of four interrelated projects. The first of these presents a new proof of the theorem of Cooper, Danciger and Wienhard classifying the limits under conjugacy of the orthogonal groups in GL(n; R). The second provides a detailed investigation into Heisenberg geometry, which is the maximally degenerate such limit in dimension two. The remaining two projects concern understanding geometric transitions which do not occur naturally as limits under conjugacy in some ambient geometry. The third project describes a new degeneration of complex hyperbolic space, formed by degenerating the complex numbers as a real algebra, into the algebra R+R. Inspired by this example, the final project attempts to build the beginnings of a framework for studying transitions between geometries abstractly. As a first application of this, we generalize the previous result and describe a collection of new geometric transitions, defined by constructing analogs of familiar geometries (projective geometry, hyperbolic geometry, etc) over real algebras, and then allowing this algebra to vary.


title: New Concept in Moisture detection with Unconventional pore morphology Design

abstract: A break in traditional pore morphology approach, is presented here to see its niche merit over the conventional sensors for water vapour detection. Tubular pores were replaced with normal cone for trace- and inverse cone for RH- level detection. The normal conical pore was fabricated by sheer manipulation of reaction rates of electrolytes, anodic polarization rate and time; and the procedure made reversed in case of inverse cone structure. Sensor with normal cone geometry exhibits response in ppm level with sensitivity of 13pF/ppm, lower detection limit(LOD)~120 ppm with excellent response/recovery time. Lowering LOD further requires alteration of conical geometric parameters in tandem with kinetic theory of water vapour molecules. In contrast, sensor developed from inverse conical structure shows response in RH level and LOD touches down to even less than 20 RH% unlike 45 RH% in conventional RH sensors. Linear response characteristics with sensitivity of 5.14 pF/RH%; surprisingly, the limitations such as nonlinear response, large response recovery time and high hysteresis as observed in conventional anodic alumina based humidity sensors have been removed. Sensing mechanism in both the structures have been suitably demonstrated and ratified with experimental data. Trace level detection is interpreted with the statistical probabilistic approach in the light of kinetic theory of gases and Brownian energy. A correlation between top surface pore diameter (through which water molecule enters) and the optimized mean free path of vapour molecule is established, and demonstrated its effectiveness for humidity detection in trace level. Results are encouraging and same concept may be tried for detection of other gaseous stimuli including organic vapours.


title: Quantales, persistence, and magnitude homology

abstract: We construct a nerve functor parametrized by a choice of quantale, exhibiting both the Vietoris-Rips complex and the magnitude nerve as instances of this nerve for different choices of monoidal structure on $\mathbb{R}$. Furthermore, the difference between how persistent homology processes the Vietoris-Rips complex and how magnitude homology processes the magnitude nerve is cast as a choice of whether or not to "localize" the corresponding nerves along $\mathbb{R}$ in a precise sense. Lastly, we mention some application-oriented observations naturally suggested by the perspective mentioned above.


title: A Framework for Measuring Weak-Lensing Magnification Using the Fundamental Plane

abstract: Galaxy-galaxy lensing is an essential tool for probing dark matter halos and constraining cosmological parameters. While galaxy-galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the fundamental plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been proposed for performing a magnification measurement. We present a formalism for using the FP residuals of elliptical galaxies to jointly estimate the foreground mass and background redshift errors for a stacked lens scenario. The FP residuals include information about weak-lensing magnification $\kappa$, and therefore foreground mass, since to first order, nonzero $\kappa$ affects galaxy size but not other FP properties. We also present a modular, extensible code that implements the formalism using emulated galaxy catalogs of a photometric galaxy survey. We find that combining FP information with observed number counts of the source galaxies constrains mass and photo-z error parameters significantly better than an estimator that includes number counts only. In particular, the constraint on the mass is 17.0\% if FP residuals are included, as opposed to 27.7\% when only number counts are included. The effective size noise for a foreground lens of mass $M_H=10^{14}M_\odot$, with a conservative selection function in size and surface brightness applied to the source population, is $\sigma_{\kappa,\mathrm{eff}}=0.250$. We discuss the improvements to our FP model necessary to make this formalism a practical companion to shear analyses in weak lensing surveys.


title: Compact binary coalescences: The subtle issue of angular momentum

abstract: In presence of gravitational radiation, the notion of angular momentum of an isolated system acquires an infinite dimensional supertranslation ambiguity. This fact has been emphasized in the mathematical general relativity literature over several decades. We analyze the issue in the restricted context of compact binary coalescence (CBC) where the initial total angular momentum of the binary and the final black hole spin generically refer to \emph{distinct} rotation subgroups of the Bondi-Metzner-Sachs group, related by \emph{supertranslations}. We show that this ambiguity can be quantified using gravitational memory and the `black hole kick'. Our results imply that, although the ambiguity is conceptually important, under assumptions normally made in the CBC literature, it can be ignored in practice for the current and foreseeable gravitational wave detectors.


title: Skeleton based simulation of turbidite channels system

abstract: A new approach to model turbidite channels using training images is presented, it is called skeleton based simulation. This is an object based model that uses some elements of multipoint geostatistics. The main idea is to simplify the representation of the training image by a one-dimensional object called training skeleton. From this new object, information about the direction and length of the channels is extracted and it is used to simulate others skeletons. These new skeletons are used to create a 3D model of channels inside a turbidite lobe.


title: Supersymmetric solutions of 7D maximal gauged supergravity

abstract: We study a number of supersymmetric solutions in the form of $Mkw_3\times S^3$- and $AdS_3\times S^3$-sliced domain walls in the maximal gauged supergravity in seven dimensions. These solutions require non-vanishing three-form fluxes to support the $AdS_3$ and $S^3$ subspaces. We consider solutions with $SO(4)$, $SO(3)$, $SO(2)\times SO(2)$ and $SO(2)$ symmetries in $CSO(p,q,5-p-q)$, $CSO(p,q,4-p-q)$ and $SO(2,1)\ltimes \mathbf{R}^4$ gauge groups. All of these solutions can be analytically obtained. For $SO(5)$ and $CSO(4,0,1)$ gauge groups, the complete truncation ansatze in terms of eleven-dimensional supergravity on $S^4$ and type IIA theory on $S^3$ are known. We give the full uplifted solutions to eleven and ten dimensions in this case. The solutions with an $AdS_3\times S^3$ slice are interpreted as two-dimensional surface defects in six-dimensional $N=(2,0)$ superconformal field theory in the case of $SO(5)$ gauge group or $N=(2,0)$ nonconformal field theories for other gauge groups. For $SO(4)$ symmetric solutions, it is possible to find solutions with both the three-form fluxes and $SO(3)$ gauge fields turned on. However, in this case, the solutions can be found only numerically. For $SO(3)$ symmetric solutions, the three-form fluxes and $SO(3)$ gauge fields cannot be non-vanishing simultaneously.


title: Scaled Autonomy: Enabling Human Operators to Control Robot Fleets

abstract: Autonomous robots often encounter challenging situations where their control policies fail and an expert human operator must briefly intervene, e.g., through teleoperation. In settings where multiple robots act in separate environments, a single human operator can manage a fleet of robots by identifying and teleoperating one robot at any given time. The key challenge is that users have limited attention: as the number of robots increases, users lose the ability to decide which robot requires teleoperation the most. Our goal is to automate this decision, thereby enabling users to supervise more robots than their attention would normally allow for. Our insight is that we can model the user's choice of which robot to control as an approximately optimal decision that maximizes the user's utility function. We learn a model of the user's preferences from observations of the user's choices in easy settings with a few robots, and use it in challenging settings with more robots to automatically identify which robot the user would most likely choose to control, if they were able to evaluate the states of all robots at all times. We run simulation experiments and a user study with twelve participants that show our method can be used to assist users in performing a simulated navigation task. We also run a hardware demonstration that illustrates how our method can be applied to a real-world mobile robot navigation task.


title: Time-dependent coupled harmonic oscillators: Comment

abstract: Macedo and Guedes showed recently how to solve a system of coupled harmonic oscillators with time dependent parameters [{ J. Math. Phys.} {\bf 53}, 052101 (2012)]. We show here that the way in which they get rid of the time dependent masses is incorrect and some terms lack in the transformed Hamiltonian. We also show a correct way of eliminating from the Hamiltonian the time dependent masses.


title: Increasing Expressivity of a Hyperspherical VAE

abstract: Learning suitable latent representations for observed, high-dimensional data is an important research topic underlying many recent advances in machine learning. While traditionally the Gaussian normal distribution has been the go-to latent parameterization, recently a variety of works have successfully proposed the use of manifold-valued latents. In one such work (Davidson et al., 2018), the authors empirically show the potential benefits of using a hyperspherical von Mises-Fisher (vMF) distribution in low dimensionality. However, due to the unique distributional form of the vMF, expressivity in higher dimensional space is limited as a result of its scalar concentration parameter leading to a 'hyperspherical bottleneck'. In this work we propose to extend the usability of hyperspherical parameterizations to higher dimensions using a product-space instead, showing improved results on a selection of image datasets.


title: Quantum mechanics at high redshift -- Modelling Damped Lyman-$\alpha$ absorption systems

abstract: For around 100 years, hydrogen spectral modelling has been based on Voigt profile fitting. The semi-classical Voigt profile is based on a 2-level atom approximation. Whilst the Voigt profile is excellent for many circumstances, the accuracy is insufficient for very high column density damped Lyman-$\alpha$ absorption systems. We have adapted the quantum-mechanical Kramers-Heisenberg model to include thermal broadening, producing a new profile, the KHT profile. Interactions involving multiple discrete atomic levels and continuum terms, not accounted for in the Voigt model, generate asymmetries in the Lyman line wings. If not modelled, this can lead to significant systematics in parameter estimation when modelling real data. There are important ramifications in particular for measurements of the primordial deuterium abundance. However, the KHT model is complicated. We therefore present a simplified formulation based on Taylor series expansions and look-up tables, quantifying the impact of the approximations made. The KHT profile has been implemented within the widely-used VPFIT code.


title: Galaxy Power Spectrum Multipoles Covariance in Perturbation Theory

abstract: We compute the covariance of the galaxy power spectrum multipoles in perturbation theory, including the effects of nonlinear evolution, nonlinear and nonlocal bias, radial redshift-space distortions, arbitrary survey window and shot noise. We rewrite the power spectrum FKP estimator in terms of the usual windowed galaxy fluctuations and the fluctuations in the number of galaxies inside the survey volume. We show that this leads to a stronger super-sample covariance than assumed in the literature and causes a substantial leakage of Gaussian information. We decompose the covariance matrix into several contributions that provide an insight into its behavior for different biased tracers. We show that for realistic surveys, the covariance of power spectrum multipoles is already dominated by shot noise and super survey mode coupling in the weakly non-linear regime. Both these effects can be accurately modeled analytically, making a perturbative treatment of the covariance very compelling. Our method allows for the covariance to be varied as a function of cosmology and bias parameters very efficiently, with survey geometry entering as fixed kernels that can be computed separately using fast fourier transforms (FFTs). We find excellent agreement between our analytic covariance and that estimated from BOSS DR12 Patchy mock catalogs in the whole range we tested, up to $k=0.6$ h/Mpc. This bodes well for application to future surveys such as DESI and Euclid. The CovaPT code that accompanies this paper is available at https://github.com/JayWadekar/CovaPT


title: Commonsense Knowledge Base Completion with Structural and Semantic Context

abstract: Automatic KB completion for commonsense knowledge graphs (e.g., ATOMIC and ConceptNet) poses unique challenges compared to the much studied conventional knowledge bases (e.g., Freebase). Commonsense knowledge graphs use free-form text to represent nodes, resulting in orders of magnitude more nodes compared to conventional KBs (18x more nodes in ATOMIC compared to Freebase (FB15K-237)). Importantly, this implies significantly sparser graph structures - a major challenge for existing KB completion methods that assume densely connected graphs over a relatively smaller set of nodes. In this paper, we present novel KB completion models that can address these challenges by exploiting the structural and semantic context of nodes. Specifically, we investigate two key ideas: (1) learning from local graph structure, using graph convolutional networks and automatic graph densification and (2) transfer learning from pre-trained language models to knowledge graphs for enhanced contextual representation of knowledge. We describe our method to incorporate information from both these sources in a joint model and provide the first empirical results for KB completion on ATOMIC and evaluation with ranking metrics on ConceptNet. Our results demonstrate the effectiveness of language model representations in boosting link prediction performance and the advantages of learning from local graph structure (+1.5 points in MRR for ConceptNet) when training on subgraphs for computational efficiency. Further analysis on model predictions shines light on the types of commonsense knowledge that language models capture well.


title: On the upper tail problem for random hypergraphs

abstract: The upper tail problem in a random graph asks to estimate the probability that the number of copies of some fixed subgraph in an Erd\H{o}s--R\'enyi random graph exceeds its expectation by some constant factor. There has been much exciting recent progress on this problem. We study the corresponding problem for hypergraphs, for which less is known about the large deviation rate. We present new phenomena in upper tail large deviations for sparse random hypergraphs that are not seen in random graphs. We conjecture a formula for the large deviation rate, i.e., the first order asymptotics of the log-probability that the number of copies of fixed subgraph $H$ in a sparse Erd\H{o}s--R\'enyi random $k$-uniform hypergraph exceeds its expectation by a constant factor. This conjecture turns out to be significantly more intricate compared to the case for graphs. We verify our conjecture when the fixed subgraph $H$ being counted is a clique, as well as when $H$ is the 3-uniform 6-vertex 4-edge hypergraph consisting of alternating faces of an octahedron, where new techniques are required.


title: Tight-Binding Studio: A Technical Software Package to Find the Parameters of Tight-Binding Hamiltonian

abstract: We present the Tight-Binding Studio (TBStudio) software package for calculating tight-binding Hamiltonian from a set of Bloch energy bands obtained from first principle theories such as density functional theory, Hartree-Fock calculations or Semi-empirical band structure theory. This will be helpful for scientists who are interested in studying the electronic properties of structures using Green's function theory in tight-binding approximation. TBStudio is a cross-platform application written in C++ with a graphical user interface design that is user-friendly and easy to work with. This software is powered by Linear Algebra Package C interface library for solving the eigenvalue problems and the standard high-performance OpenGL graphic library for real-time plotting. TBStudio and its examples together with the tutorials are available for download from tight-binding.com.


title: A Quadratic Programming Flux Correction Method for High-Order DG Discretizations of SN Transport

abstract: We present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin discretizations of the SN transport equation. This approach is sweep-compatible: as the transport sweep is performed, a local quadratic programming (QP) problem is solved in each spatial cell to ensure that the solution satisfies certain physical constraints, including local particle balance. The constraints can be chosen in two ways, leading to two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). The coefficients of the solution are constrained to be nonnegative in QPZ, and they are constrained by an approximate discrete maximum principle in QPMP. There are two primary takeaways in this paper. First, it is shown that the QPMP method, when used with the positive Bernstein basis, eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The second takeaway is that the Variable Eddington Factor (VEF) method can be used to accelerate the convergence of source iteration with fixup for problems with optically thick regions. Regardless of whether a fixup is applied, source iteration converges slowly when optically thick regions are present and acceleration is needed. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.


title: Multi-step Greedy Reinforcement Learning Algorithms

abstract: Multi-step greedy policies have been extensively used in model-based reinforcement learning (RL), both when a model of the environment is available (e.g.,~in the game of Go) and when it is learned. In this paper, we explore their benefits in model-free RL, when employed using multi-step dynamic programming algorithms: $\kappa$-Policy Iteration ($\kappa$-PI) and $\kappa$-Value Iteration ($\kappa$-VI). These methods iteratively compute the next policy ($\kappa$-PI) and value function ($\kappa$-VI) by solving a surrogate decision problem with a shaped reward and a smaller discount factor. We derive model-free RL algorithms based on $\kappa$-PI and $\kappa$-VI in which the surrogate problem can be solved by any discrete or continuous action RL method, such as DQN and TRPO. We identify the importance of a hyper-parameter that controls the extent to which the surrogate problem is solved and suggest a way to set this parameter. When evaluated on a range of Atari and MuJoCo benchmark tasks, our results indicate that for the right range of $\kappa$, our algorithms outperform DQN and TRPO. This shows that our multi-step greedy algorithms are general enough to be applied over any existing RL algorithm and can significantly improve its performance.


title: Topological Photonic States at a 1-D Binary-Quaternary Interface

abstract: The existence of topological interface states is investigated at the boundary between a binary photonic crystal and a quaternary photonic crystal, with each possessing inversion symmetric unit cells. Conditions are established that describe where the quaternary crystal can exist in parameter space subject to constraints. The closing of band gaps is discussed for different optical path ratios. When the binary and quaternary crystals share an interface, optical states appear at the interface when the two crystals have different signs for the surface impedance. The evolution of the states is displayed as the geometry of the quaternary crystal changes.


title: Renormalized energy between vortices in some Ginzburg-Landau models on 2-dimensional Riemannian manifolds

abstract: We study a variational Ginzburg-Landau type model depending on a small parameter $\varepsilon>0$ for (tangent) vector fields on a $2$-dimensional Riemannian manifold $S$. As $\varepsilon\to 0$, these vector fields tend to have unit length so they generate singular points, called vortices, of a (non-zero) index if the genus $\mathfrak{g}$ of $S$ is different than $1$. Our first main result concerns the characterization of canonical harmonic unit vector fields with prescribed singular points and indices. The novelty of this classification involves flux integrals constrained to a particular vorticity-dependent lattice in the $2\mathfrak{g}$-dimensional space of harmonic $1$-forms on $S$ if $\mathfrak{g}\geq 1$. Our second main result determines the interaction energy (called renormalized energy) between vortex points as a $\Gamma$-limit (at the second order) as $\varepsilon\to 0$. The renormalized energy governing the optimal location of vortices depends on the Gauss curvature of $S$ as well as on the quantized flux. The coupling between flux quantization constraints and vorticity, and its impact on the renormalized energy, are new phenomena in the theory of Ginzburg-Landau type models. We also extend this study to two other (extrinsic) models for embedded hypersurfaces $S\subset \mathbb{R}^3$, in particular, to a physical model for non-tangent maps to $S$ coming from micromagnetics.


title: Global sensitivity analysis of bulk properties of an atomic nucleus

abstract: We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus $^{16}$O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)]. With this method we compute the binding energy and charge radius of $^{16}$O at more than one million different values of the 16 low-energy constants in one hour on a standard laptop. For relatively small subspace projections, the root-mean-square error is about 1% compared to full space coupled-cluster results. We find that 58(1)% of the variance in the energy can be apportioned to a single contact-term in the $^3S_1$-wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation, and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up for an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segr\'e chart.


title: A Survey on Active Learning and Human-in-the-Loop Deep Learning for Medical Image Analysis

abstract: Fully automatic deep learning has become the state-of-the-art technique for many tasks including image acquisition, analysis and interpretation, and for the extraction of clinically useful information for computer-aided detection, diagnosis, treatment planning, intervention and therapy. However, the unique challenges posed by medical image analysis suggest that retaining a human end user in any deep learning enabled system will be beneficial. In this review we investigate the role that humans might play in the development and deployment of deep learning enabled diagnostic applications and focus on techniques that will retain a significant input from a human end user. Human-in-the-Loop computing is an area that we see as increasingly important in future research due to the safety-critical nature of working in the medical domain. We evaluate four key areas that we consider vital for deep learning in the clinical practice: (1) Active Learning to choose the best data to annotate for optimal model performance; (2) Interaction with model outputs - using iterative feedback to steer models to optima for a given prediction and offering meaningful ways to interpret and respond to predictions; (3) Practical considerations - developing full scale applications and the key considerations that need to be made before deployment; (4) Future Prospective and Unanswered Questions - knowledge gaps and related research fields that will benefit human-in-the-loop computing as they evolve. We offer our opinions on the most promising directions of research and how various aspects of each area might be unified towards common goals.