title: The Quantum Random Energy Model as a Limit of p-Spin Interactions

abstract: We consider the free energy of a mean-field quantum spin glass described by a $ p $-spin interaction and a transversal magnetic field. Recent rigorous results for the case $ p= \infty $, i.e. the quantum random energy model (QREM), are reviewed. We show that the free energy of the $ p $-spin model converges in a joint thermodynamic and $ p \to \infty $ limit to the free energy of the QREM.


title: Realisation of a quantum current standard at liquid helium temperature with sub-ppm reproducibility

abstract: A silicon electron pump operating at the temperature of liquid helium has demonstrated repeatable operation with sub-ppm accuracy. The pump current, approximately 168 pA, is measured by three laboratories, and the measurements agree with the expected current ef within the uncertainties which range from 0.2 ppm to 1.3 ppm. All the measurements are carried out in zero applied magnetic field, and the pump drive signal is a sine wave. The combination of simple operating conditions with high accuracy demonstrates the possibility that an electron pump can operate as a current standard in a National Measurement Institute. We also discuss other practical aspects of using the electron pump as a current standard, such as testing its robustness to changes in the control parameters, and using a rapid tuning procedure to locate the optimal operation point..


title: Modelling equilibration of local many-body quantum systems by random graph ensembles

abstract: We introduce structured random matrix ensembles, constructed to model many-body quantum systems with local interactions. These ensembles are employed to study equilibration of isolated many-body quantum systems, showing that rather complex matrix structures, well beyond Wigner's full or banded random matrices, are required to faithfully model equilibration times. Viewing the random matrices as connectivities of graphs, we analyse the resulting network of classical oscillators in Hilbert space with tools from network theory. One of these tools, called the maximum flow value, is found to be an excellent proxy for equilibration times. Since maximum flow values are less expensive to compute, they give access to approximate equilibration times for system sizes beyond those accessible by exact diagonalisation.


title: Sequences of consecutive factoradic happy numbers

abstract: Given a positive integer $n$, the factorial base representation of $n$ is given by $n=\sum_{i=1}^ka_i\cdot i!$, where $a_k\neq 0$ and $0\leq a_i\leq i$ for all $1\leq i\leq k$. For $e\geq 1$, we define $S_{e,!}:\mathbb{Z}_{\geq0}\to\mathbb{Z}_{\geq0}$ by $S_{e,!}(0) = 0$ and $S_{e,!}(n)=\sum_{i=0}^{n}a_i^e$, for $n \neq 0$. For $\ell\geq 0$, we let $S_{e,!}^\ell(n)$ denote the $\ell$-th iteration of $S_{e,!}$, while $S_{e,!}^0(n)=n$. If $p\in\mathbb{Z}^+$ satisfies $S_{e,!}(p)=p$, then we say that $p$ is an $e$-power factoradic fixed point of $S_{e,!}$. Moreover, given $x\in \mathbb{Z}^+$, if $p$ is an $e$-power factoradic fixed point and if there exists $\ell\in \mathbb{Z}_{\geq 0}$ such that $S_{e,!}^\ell(x)=p$, then we say that $x$ is an $e$-power factoradic $p$-happy number. Note an integer $n$ is said to be an $e$-power factoradic happy number if it is an $e$-power factoradic $1$-happy number. In this paper, we prove that all positive integers are $1$-power factoradic happy and, for $2\leq e\leq 4$, we prove the existence of arbitrarily long sequences of $e$-power factoradic $p$-happy numbers. A curious result establishes that for any $e\geq 2$ the $e$-power factoradic fixed points of $S_{e,!}$ that are greater than $1$, always appear in sets of consecutive pairs. Our last contribution, provides the smallest sequences of $m$ consecutive $e$-power factoradic happy numbers for $2\leq e\leq 5$, for some values of $m$.


title: Privacy-Preserving Search for a Similar Genomic Makeup in the Cloud

abstract: In this paper, we attempt to provide a privacy-preserving and efficient solution for the "similar patient search" problem among several parties (e.g., hospitals) by addressing the shortcomings of previous attempts. We consider a scenario in which each hospital has its own genomic dataset and the goal of a physician (or researcher) is to search for a patient similar to a given one (based on a genomic makeup) among all the hospitals in the system. To enable this search, we let each hospital encrypt its dataset with its own key and outsource the storage of its dataset to a public cloud. The physician can get authorization from multiple hospitals and send a query to the cloud, which efficiently performs the search across authorized hospitals using a privacy-preserving index structure. We propose a hierarchical index structure to index each hospital's dataset with low memory requirements. Furthermore, we develop a novel privacy-preserving index merging mechanism that generates a common search index from individual indices of each hospital to significantly improve the search efficiency. We also consider the storage of medical information associated with genomic data of a patient (e.g., diagnosis and treatment). We allow access to this information via a fine-grained access control policy that we develop through the combination of standard symmetric encryption and ciphertext policy attribute-based encryption. Using this mechanism, a physician can search for similar patients and obtain medical information about the matching records if the access policy holds. We conduct experiments on large-scale genomic data and show the efficiency of the proposed scheme. Notably, we show that under our experimental settings, the proposed scheme is more than $60$ times faster than Wang et al.'s protocol and $95$ times faster than Asharov et al.'s solution.


title: Elementary laser-less quantum logic operations with (anti-)protons in Penning traps

abstract: Static magnetic field gradients superimposed on the electromagnetic trapping potential of a Penning trap can be used to implement laser-less spin-motion couplings that allow the realization of elementary quantum logic operations in the radio-frequency regime. An important scenario of practical interest is the application to $g$-factor measurements with single (anti-)protons to test the fundamental charge, parity, time reversal (CPT) invariance as pursued in the BASE collaboration [Smorra et al., Eur. Phys. J. Spec. Top. 224, 3055-3108 (2015), Smorra et al., Nature 550, 371-374 (2017), Schneider et al., Science 358, 1081-1084 (2017)]. We discuss the classical and quantum behavior of a charged particle in a Penning trap with a superimposed magnetic field gradient. Using analytic and numerical calculations, we find that it is possible to carry out a SWAP gate between the spin and the motional qubit of a single (anti-)proton with high fidelity, provided the particle has been initialized in the motional ground state. We discuss the implications of our findings for the realization of quantum logic spectroscopy in this system.


title: Neural Machine Translation: A Review and Survey

abstract: The field of machine translation (MT), the automatic translation of written text from one natural language into another, has experienced a major paradigm shift in recent years. Statistical MT, which mainly relies on various count-based models and which used to dominate MT research for decades, has largely been superseded by neural machine translation (NMT), which tackles translation with a single neural network. In this work we will trace back the origins of modern NMT architectures to word and sentence embeddings and earlier examples of the encoder-decoder network family. We will conclude with a survey of recent trends in the field.


title: Siamese Natural Language Tracker: Tracking by Natural Language Descriptions with Siamese Trackers

abstract: We propose a novel Siamese Natural Language Tracker (SNLT), which brings the advancements in visual tracking to the tracking by natural language (NL) descriptions task. The proposed SNLT is applicable to a wide range of Siamese trackers, providing a new class of baselines for the tracking by NL task and promising future improvements from the advancements of Siamese trackers. The carefully designed architecture of the Siamese Natural Language Region Proposal Network (SNL-RPN), together with the Dynamic Aggregation of vision and language modalities, is introduced to perform the tracking by NL task. Empirical results over tracking benchmarks with NL annotations show that the proposed SNLT improves Siamese trackers by 3 to 7 percentage points with a slight tradeoff of speed. The proposed SNLT outperforms all NL trackers to-date and is competitive among state-of-the-art real-time trackers on LaSOT benchmarks while running at 50 frames per second on a single GPU.


title: On Even Rainbow or Nontriangular Directed Cycles

abstract: Let $G = (V, E)$ be an $n$-vertex edge-colored graph. In 2013, H. Li proved that if every vertex $v \in V$ is incident to at least $(n+1)/2$ distinctly colored edges, then $G$ admits a rainbow triangle. We establish a corresponding result for fixed even rainbow $\ell$-cycles $C_{\ell}$: if every vertex $v \in V$ is incident to at least $(n+5)/3$ distinctly colored edges, where $n \geq n_0(\ell)$ is sufficiently large, then $G$ admits an even rainbow $\ell$-cycle $C_{\ell}$. This result is best possible whenever $\ell \not\equiv 0$ (mod 3). Correspondingly, we also show that for a fixed (even or odd) integer $\ell \geq 4$, every large $n$-vertex oriented graph $\vec{G} = (V, \vec{E})$ with minimum outdegree at least $(n+1)/3$ admits a (consistently) directed $\ell$-cycle $\vec{C}_{\ell}$. Our latter result relates to one of Kelly, K\"uhn, and Osthus, who proved a similar statement for oriented graphs with large semi-degree. Our proofs are based on the stability method.


title: SimAS: A Simulation-assisted Approach for the Scheduling Algorithm Selection under Perturbations

abstract: Many scientific applications consist of large and computationally-intensive loops. Dynamic loop self-scheduling (DLS) techniques are used to parallelize and to balance the load during the execution of such applications. Load imbalance arises from variations in the loop iteration (or tasks) execution times, caused by problem, algorithmic, or systemic characteristics. The variations in systemic characteristics are referred to as perturbations, and can be caused by other applications or processes that share the same resources, or a temporary system fault or malfunction. Therefore, the selection of the most efficient DLS technique is critical to achieve the best application performance. The following question motivates this work: Given an application, an HPC system, and their characteristics and interplay, which DLS technique will achieve improved performance under unpredictable perturbations? Existing studies focus on variations in the delivered computational speed only as the source of perturbations in the system. However, perturbations in available network bandwidth or latency are inevitable on production HPC systems. A Simulator-assisted scheduling (SimAS) is introduced as a new control-theoretic-inspired approach to dynamically select DLS techniques that improve the performance of applications executing on heterogeneous HPC systems under perturbations. The present work examines the performance of seven applications on a heterogeneous system under all the above system perturbations. SimAS is evaluated as a proof of concept using native and simulative experiments. The performance results confirm the original hypothesis that no single DLS technique can deliver the absolute best performance in all scenarios, whereas the SimAS-based DLS selection resulted in improved application performance in most experiments.


title: Asymptotics for Strassen's Optimal Transport Problem

abstract: In this paper, we consider Strassen's version of optimal transport (OT) problem, which concerns minimizing the excess-cost probability (i.e., the probability that the cost is larger than a given value) over all couplings of two given distributions. We derive large deviation, moderate deviation, and central limit theorems for this problem. Our proof is based on Strassen's dual formulation of the OT problem, Sanov's theorem on the large deviation principle (LDP) of empirical measures, as well as the moderate deviation principle (MDP) and central limit theorems (CLT) of empirical measures. In order to apply the LDP, MDP, and CLT to Strassen's OT problem, nested formulas for Strassen's OT problem are derived. Based on these nested formulas and using a splitting technique, we construct asymptotically optimal solutions to Strassen's OT problem and its dual formulation.


title: Predictions for $B_s \to \bar{K}^* \ell \,\ell$ in non-universal $Z'$ models

abstract: The lepton flavor universality violating (LFUV) measurements $R_K$ and $R_{K^*}$ in $B$ meson decays can be accounted for in non-universal $Z'$ models. We constrain the couplings of these $Z'$ models by performing a global fit to correlated $b \to s \ell \ell$ and $b \to d \ell \ell $ processes, and calculate their possible implications for $B_s \to \bar{K}^*\ell \ell$ observables. For real new physics (NP) couplings, the 1-$\sigma$ favored parameters allow the corresponding LFUV ratio $R_{K^*}^{(s)}$ in $B_s \to \bar{K}^*\ell \ell$ to range between 0.8 -- 1.2 at low $q^2$. Complex NP couplings improve the best fit only marginally, however they allow a significant enhancement of the branching ratio, while increasing the range of $R_{K^*}^{(s)}$ at low $q^2$ to 0.8 -- 1.8. We find that NP could cause zero-crossing in the forward-backward asymmetry $A_{FB}$ to shift towards lower $q^2$ values, and enhancement in the magnitude of integrated $A_{FB}$. The $CP$ asymmetry $A_{CP}$ may be suppressed and even change sign. The simultaneous measurements of integrated $R_{K^*}^{(s)}$ and $A_{CP}$ values to 0.1 and 1% respectively, would help in constraining the effective NP Wilson coefficient $C_9$ in $ b \to d \mu \mu$ interactions.


title: Tensor-network method to simulate strongly interacting quantum thermal machines

abstract: We present a methodology to simulate the quantum thermodynamics of thermal machines which are built from an interacting working medium in contact with fermionic reservoirs at fixed temperature and chemical potential. Our method works at finite temperature, beyond linear response and weak system-reservoir coupling, and allows for non-quadratic interactions in the working medium. The method uses mesoscopic reservoirs, continuously damped towards thermal equilibrium, in order to represent continuum baths and a novel tensor network algorithm to simulate the steady-state thermodynamics. Using the example of a quantum-dot heat engine, we demonstrate that our technique replicates the well known Landauer-B\"uttiker theory for efficiency and power. We then go beyond the quadratic limit to demonstrate the capability of our method by simulating a three-site machine with non-quadratic interactions. Remarkably, we find that such interactions lead to power enhancement, without being detrimental to the efficiency. Furthermore, we demonstrate the capability of our method to tackle complex many-body systems by extracting the super-diffusive exponent for high-temperature transport in the isotropic Heisenberg model. Finally, we discuss transport in the gapless phase of the anisotropic Heisenberg model at finite temperature and its connection to charge conjugation-parity, going beyond the predictions of single-site boundary driving configurations.


title: New physics searches with top quarks

abstract: This is a brief pedagogic introduction to searches for new physics with top quarks. It covers indirect searches for heavy new particles based on standard model effective theory and direct searches for new signatures of a light hidden sector. LHC and flavor observables complement and strengthen each other in this endeavor.


title: Ahlfors regular spaces have regular subspaces of any dimension

abstract: We characterize $Q$-dimensional Ahlfors regular spaces among trees' boundaries and show how to construct, for each $0 < \alpha < Q$, an $\alpha$-regular subspace. As an application, we give an alternative simple proof of the existence of $\alpha$-regular subspaces of a $Q$-dimensional complete Ahlfors regular metric space $(X,\rho)$, which was proved in \cite{JJKRRS}.


title: Exploring the latitude and depth dependence of solar Rossby waves using ring-diagram analysis

abstract: Global-scale Rossby waves have recently been unambiguously identified on the Sun. Here we study the latitude and depth dependence of the Rossby wave eigenfunctions. By applying helioseismic ring-diagram analysis and granulation tracking to SDO/HMI observations, we compute maps of the radial vorticity of flows in the upper solar convection zone (down to depths of more than $16$ Mm). We use a Fourier transform in longitude to separate the different azimuthal orders m in the range $3 \le m \le 15$. At each $m$ we obtain the phase and amplitude of the Rossby waves as a function of depth using the helioseismic data. At each $m$ we also measure the latitude dependence of the eigenfunctions by calculating the covariance between the equator and other latitudes. We then study the horizontal and radial dependences of the radial vorticity eigenfunctions. The horizontal eigenfunctions are complex. As observed previously, the real part peaks at the equator and switches sign near $\pm 30^\circ$, thus the eigenfunctions show significant non-sectoral contributions. The imaginary part is smaller than the real part. The phase of the radial eigenfunctions varies by only roughly $\pm 5^\circ$ over the top $15$ Mm. The amplitude of the radial eigenfunctions decreases by about $10\%$ from the surface down to $8$ Mm (the region where ring-diagram analysis is most reliable, as seen by comparing with the rotation rate measured by global-mode seismology). The radial dependence of the radial vorticity eigenfunctions deduced from ring-diagram analysis is consistent with a power-law down to $8$ Mm and is unreliable at larger depths. However, the observations provide only weak constraints on the power-law exponents. For the real part, the latitude dependence of the eigenfunctions is consistent with previous work (using granulation tracking). The imaginary part is smaller than the real part but significantly nonzero.


title: RTN: Reparameterized Ternary Network

abstract: To deploy deep neural networks on resource-limited devices, quantization has been widely explored. In this work, we study the extremely low-bit networks which have tremendous speed-up, memory saving with quantized activation and weights. We first bring up three omitted issues in extremely low-bit networks: the squashing range of quantized values; the gradient vanishing during backpropagation and the unexploited hardware acceleration of ternary networks. By reparameterizing quantized activation and weights vector with full precision scale and offset for fixed ternary vector, we decouple the range and magnitude from the direction to extenuate the three issues. Learnable scale and offset can automatically adjust the range of quantized values and sparsity without gradient vanishing. A novel encoding and computation pat-tern are designed to support efficient computing for our reparameterized ternary network (RTN). Experiments on ResNet-18 for ImageNet demonstrate that the proposed RTN finds a much better efficiency between bitwidth and accuracy, and achieves up to 26.76% relative accuracy improvement compared with state-of-the-art methods. Moreover, we validate the proposed computation pattern on Field Programmable Gate Arrays (FPGA), and it brings 46.46x and 89.17x savings on power and area respectively compared with the full precision convolution.


title: Logarithmic corrections to $\mathbf{a^2}$ scaling in lattice Yang Mills theory

abstract: We analyse the leading logarithmic corrections to the $a^2$ scaling of lattice artefacts in QCD, following the seminal work of Balog, Niedermayer and Weisz in the O(n) non-linear sigma model. Restricting our attention to contributions from the action, the leading logarithmic corrections can be determined by the anomalous dimensions of a minimal on-shell basis of mass-dimension 6 operators. We present results for the SU(N) pure gauge theory. In this theory the logarithmic corrections reduce the cutoff effects. These computations are the first step towards a study of full lattice QCD at O($a^2$), which is in progress.


title: Learning to Dynamically Coordinate Multi-Robot Teams in Graph Attention Networks

abstract: Increasing interest in integrating advanced robotics within manufacturing has spurred a renewed concentration in developing real-time scheduling solutions to coordinate human-robot collaboration in this environment. Traditionally, the problem of scheduling agents to complete tasks with temporal and spatial constraints has been approached either with exact algorithms, which are computationally intractable for large-scale, dynamic coordination, or approximate methods that require domain experts to craft heuristics for each application. We seek to overcome the limitations of these conventional methods by developing a novel graph attention network formulation to automatically learn features of scheduling problems to allow their deployment. To learn effective policies for combinatorial optimization problems via machine learning, we combine imitation learning on smaller problems with deep Q-learning on larger problems, in a non-parametric framework, to allow for fast, near-optimal scheduling of robot teams. We show that our network-based policy finds at least twice as many solutions over prior state-of-the-art methods in all testing scenarios.


title: Limit polygons of convex domains in the projective plane

abstract: For any sequence of properly convex domains in the real projective plane such that the zeros of Pick differentials have bounded multiplicity and get further and further apart, we determined all Hausdorff limit domains that one can obtain after normalizing each member of the sequence by a projective transformation. We then show that the result can be applied to convex domains generated by projective triangular reflection groups.


title: The CMS Outer Tracker for the High Luminosity LHC upgrade

abstract: The High Luminosity Large Hadron Collider (HL-LHC) at CERN is expected to collide protons at a centre-of-mass energy of 14\,TeV and to reach the unprecedented peak instantaneous luminosity of 5\,$-$\,7.5\,x\,$10^{34}$\,cm$^{-2}$s$^{-1}$ with an average number of pileup events of 140\,$-$\,200. This will allow the ATLAS and CMS experiments to collect integrated luminosities up to 3000\,$-$\,4000\,fb$^{-1}$ during the project lifetime. To cope with this extreme scenario the CMS detector will be substantially upgraded before starting the HL-LHC, with a plan known as CMS Phase-2 upgrade. The CMS Tracker detector will have to be replaced in order to fully exploit the delivered luminosity and cope with the demanding operating conditions. The new detector will provide robust tracking as well as input for the first level trigger. This paper is focused on the replacement of the CMS Outer Tracker system, describing the new layout and the technological choices together with some highlights of module assembly and quality assurance aspects.


title: Unraveling the degradation mechanism of FIrpic based blue OLEDs: I. A theoretical investigation

abstract: We report a detailed ab-initio study of the of the microscopic degradation mechanism of FIrpic, a popular blue emitter in OLED devices. We simulate the \emph{operando} conditions of FIrpic by adding an electron-hole pair (exciton) to the system. We perform both static calculations with the TDDFT framework and we also simulate the evolution of the system at finite temperature via Car-Parrinello molecular dynamics. We found triplet excitons are very effective in reducing the Ir-N bond breaking barrier of the picolinate moiety. After the first bond breaking, the two oxygen of picolinate swap their position and FIrpic can either remain stable in an "open" configuration, or loose a picolinate fragment, which at a later stage can evolve a CO$_2$ molecule. Our method can be applied to other light emitting Ir-complexes in order to quickly estimate their stability in OLED devices. In Paper~II we complement our theoretical study with a parallel experimental investigation of the key degradation steps of FIrpic in an aged device.


title: 3-Dimensional TQFTs From Non-Semisimple Modular Categories

abstract: We use modified traces to renormalize Lyubashenko's closed 3-manifold invariants coming from twist non-degenerate finite unimodular ribbon categories. Our construction produces new topological invariants which we upgrade to 2+1-TQFTs under the additional assumption of factorizability. The resulting functors provide monoidal extensions of Lyubashenko's mapping class group representations, as discussed in arXiv:2010.14852. This general framework encompasses important examples of non-semisimple modular categories coming from the representation theory of quasi-Hopf algebras, which were left out of previous non-semisimple TQFT constructions.


title: Volterra Equations Driven by Rough Signals

abstract: This article is devoted to the extension of the theory of rough paths in the context of Volterra equations with possibly singular kernels. We begin to describe a class of two parameter functions defined on the simplex called Volterra paths. These paths are used to construct a so-called Volterra-signature, analogously to the signature used in Lyon's theory of rough paths. We provide a detailed algebraic and analytic description of this object. Interestingly, the Volterra signature does not have a multiplicative property similar to the classical signature, and we introduce an integral product behaving like a convolution extending the classical tensor product. We show that this convolution product is well defined for a large class of Volterra paths, and we provide an analogue of the extension theorem from the theory of rough paths (which guarantees in particular the existence of a Volterra signature). Moreover the concept of convolution product is essential in the construction of Volterra controlled paths, which is the natural class of processes to be integrated with respect to the driving noise in our situation. This leads to a rough integral given as a functional of the Volterra signature and the Volterra controlled paths, combined through the convolution product. The rough integral is then used in the construction of solutions to Volterra equations driven by H\"older noises with singular kernels. An example concerning Brownian noises and a singular kernel is treated


title: Safety and Robustness in Decision Making: Deep Bayesian Recurrent Neural Networks for Somatic Variant Calling in Cancer

abstract: The genomic profile underlying an individual tumor can be highly informative in the creation of a personalized cancer treatment strategy for a given patient; a practice known as precision oncology. This involves next generation sequencing of a tumor sample and the subsequent identification of genomic aberrations, such as somatic mutations, to provide potential candidates of targeted therapy. The identification of these aberrations from sequencing noise and germline variant background poses a classic classification-style problem. This has been previously broached with many different supervised machine learning methods, including deep-learning neural networks. However, these neural networks have thus far not been tailored to give any indication of confidence in the mutation call, meaning an oncologist could be targeting a mutation with a low probability of being true. To address this, we present here a deep bayesian recurrent neural network for cancer variant calling, which shows no degradation in performance compared to standard neural networks. This approach enables greater flexibility through different priors to avoid overfitting to a single dataset. We will be incorporating this approach into software for oncologists to obtain safe, robust, and statistically confident somatic mutation calls for precision oncology treatment choices.


title: Gaussian trajectory approach to dissipative phase transitions: the case of quadratically driven photonic lattices

abstract: We apply the Gaussian trajectories approach to the study of the critical behavior of two-dimensional dissipative arrays of nonlinear photonic cavities, in presence of two-photon driving and in regimes of sizable loss rates. In spite of the highly mixed character of the density matrix of this system, the numerical approach is able to provide precise estimations of the steady-state expectation values, even for large lattices made of more than 100 sites. By performing a finite-size scaling of the relevant properties of the steady state, we extrapolate the behavior of the system in the thermodynamic limit and we show the emergence of a second-order dissipative phase transition, belonging to the universality class of thermal Ising model. This result indicates the occurrence of a crossover when the loss rate is increased from the weak-loss limit, in which the phase transition belongs to the universality class of the quantum Ising model


title: Improved Bhabha cross section at LEP and the number of light neutrino species

abstract: In $\rm e^+ \rm e^-$ collisions, the integrated luminosity is generally measured from the rate of low-angle Bhabha interactions $\rm e^+ e^- \to e^+ e^- $. In the published LEP results, the inferred theoretical uncertainty of $\pm 0.061\%$ on the predicted rate is significantly larger than the reported experimental uncertainties. We present an updated and more accurate prediction of the Bhabha cross section in this letter, which is found to reduce the Bhabha cross section by about $0.048\%$, and its uncertainty to $\pm 0.037\%$. When accounted for, these changes modify the number of light neutrino species (and its accuracy), as determined from the LEP measurement of the hadronic cross section at the Z peak, to $N_\nu = 2.9963 \pm 0.0074$. The 20-years-old $2 \sigma$ tension with the Standard Model is gone.


title: Embedding graphs into larger graphs: results, methods, and problems

abstract: Extremal Graph Theory is a very deep and wide area of modern combinatorics. It is very fast developing, and in this long but relatively short survey we select some of those results which either we feel very important in this field or which are new breakthrough results, or which --- for some other reasons --- are very close to us. Some results discussed here got stronger emphasis, since they are connected to Lov\'asz (and sometimes to us).


title: Graph signal interpolation with Positive Definite Graph Basis Functions

abstract: For the interpolation of graph signals with generalized shifts of a graph basis function (GBF), we introduce the concept of positive definite functions on graphs. This concept merges kernel-based interpolation with spectral theory on graphs and can be regarded as a graph analog of radial basis function interpolation in euclidean spaces or spherical basis functions. We provide several descriptions of positive definite functions on graphs, the most relevant one is a Bochner-type characterization in terms of positive Fourier coefficients. These descriptions allow us to design GBF's and to study GBF interpolation in more detail: we are able to characterize the native spaces of the interpolants, we provide explicit estimates for the interpolation error and obtain bounds for the numerical stability. As a final application, we show how GBF interpolation can be used to get quadrature formulas on graphs.


title: On the emergence of gravitational dynamics from tensor networks

abstract: Tensor networks are powerful techniques that widely used in condensed matter physics. In this language, the wave function of a quantum manybody system is described by a network of tensors with specific entanglement structures. Recently, it is shown that tensor network can generate the anti-de Sitter (AdS) geometry by using the entanglement renormalization approach. However, whether the dynamical connections can be found between the tensor network and the gravity is an important unsolved problem. In this paper, we give a novel proposal to integrate ideas from tensor networks, entanglement entropy, canonical quantization of quantum gravity and the holographic principle and argue that the gravitational dynamics can be generated from a tensor network if the wave function of the latter satisfies the Wheeler-DeWitt equation.


title: Optimization in Software Engineering -- A Pragmatic Approach

abstract: Empirical software engineering is concerned with the design and analysis of empirical studies that include software products, processes, and resources. Optimization is a form of data analytics in support of human decision-making. Optimization methods are aimed to find the best decision alternatives. Empirical studies serve both as a model and as data input for optimization. In addition, the complexity of the models used for optimization trigger further studies on explaining and validating the results in real-world scenarios. The goal of this chapter is to give an overview of the as-is and of the to-be usage of optimization in software engineering. The emphasis is on pragmatic use of optimization, and not so much on describing the most recent algorithmic innovations and tool developments. The usage of optimization covers a wide range of questions from different types of software engineering problems along the whole life-cycle. To facilitate its more comprehensive and more effective usage, a checklist for a guided process is described. The chapter uses a running example Asymmetric Release Planning to illustrate the whole process. A Return-on-Investment analysis is proposed as part of the problem scoping. This helps to decide on the depth and breadth of analysis in relation to the effort needed to run the analysis and the projected value of the solution.


title: Finding entries of maximum absolute value in low-rank tensors

abstract: We present an iterative method for the search of extreme entries in low-rank tensors which is based on a power iteration combined with a binary search. In this work we use the HT-format for low-rank tensors but other low-rank formats can be used verbatim. We present two different approaches to accelerate the basic power iteration: an orthogonal projection of Rayleigh-Ritz type, as well as an acceleration of the power iteration itself which can be achieved due to the diagonal structure of the underlying eigenvalue problem. Finally the maximizing index is determined by a binary search based on the proposed iterative method for the approximation of the maximum norm. The iterative method for the maximum norm estimation inherits the linear complexity of the involved tensor arithmetic in the HT-format w.r.t. the tensor order, which is also verified by numerical tests.


title: Temperature-dependent fluorescence emission spectra of acrylic (PMMA) and tetraphenyl butadiene (TPB) excited with UV light

abstract: Acrylic (poly(methyl methacrylate) or PMMA) is commonly used as a vessel to hold scintillating liquids in rare-event searches. Certain types of PMMA can fluoresce with a low efficiency at room temperature. We have investigated the fluorescence spectra under 280nm ultraviolet (UV) excitation of the PMMA used in the DEAP dark matter search at various low temperatures. Fluorescence of this PMMA is observed, with an increasing intensity as the sample is cooled from a temperature of 300K down to 4K. The common wavelength-shifter tetraphenyl butadiene (TPB) is also measured for use as a reference.


title: AlgaeDICE: Policy Gradient from Arbitrary Experience

abstract: In many real-world applications of reinforcement learning (RL), interactions with the environment are limited due to cost or feasibility. This presents a challenge to traditional RL algorithms since the max-return objective involves an expectation over on-policy samples. We introduce a new formulation of max-return optimization that allows the problem to be re-expressed by an expectation over an arbitrary behavior-agnostic and off-policy data distribution. We first derive this result by considering a regularized version of the dual max-return objective before extending our findings to unregularized objectives through the use of a Lagrangian formulation of the linear programming characterization of Q-values. We show that, if auxiliary dual variables of the objective are optimized, then the gradient of the off-policy objective is exactly the on-policy policy gradient, without any use of importance weighting. In addition to revealing the appealing theoretical properties of this approach, we also show that it delivers good practical performance.


title: Charge density wave sliding driven by an interplay of conventional and Hall voltages in NbSe$_3$ microbridges

abstract: Collective charge-density wave (CDW) transport was measured under a high magnetic field in NbSe$_3$ microbridges which have been cut transversely and at an angle to the chains' direction. We give evidences that the CDW sliding is driven by the Hall voltage generated by the inter-chain current of normal carriers. We have discovered a re-entrance effect of the Hall-driven sliding above a crossover temperature at which the Hall constant has been known to change sign. For the narrow channel, cut at 45$^\circ$ relative to the chain axis, we observed an evolution from the Hall-driven sliding at low temperatures, to the conventional sliding at higher temperatures, which corroborates with falling of the Hall constant. In this course, the nonlinear contribution to the conductivity coming from the collective sliding changes sign. The quantization of Shapiro-steps, generated presumably by a coherent sequence of phase slips, indicates that their governing changes from the applied voltage to the current.


title: UEFA against the champions? An evaluation of the recent reform of the Champions League qualification

abstract: The UEFA Champions League is the major European club football competition, however, most national champions have to play qualification matches to receive a slot in the group stage of the tournament. The paper evaluates the impact of the only reform in the Champions Path of the qualifying system since 2009, effective from the 2018/19 season. While it is anticipated that the reduction in the number of berths decreases the probability of advancing to the group stage, the distribution of the effects among the national associations is revealed here via Monte-Carlo simulations. In contrast to previous studies, our methodology considers five seasons instead of only one to filter out any possible season-specific attributes. Almost all countries are found to gain less prize money on average. Several champions, including the Bulgarian, the Scottish, and, especially, the Swiss, might lose million Euros. As the consequences depend to a large extent on the arbitrary cutoffs in the access list, we propose to introduce some randomness into the determination of entry stages in the qualification.


title: PDC -- a probabilistic distributional clustering algorithm: a case study on suicide articles in PubMed

abstract: The need to organize a large collection in a manner that facilitates human comprehension is crucial given the ever-increasing volumes of information. In this work, we present PDC (probabilistic distributional clustering), a novel algorithm that, given a document collection, computes disjoint term sets representing topics in the collection. The algorithm relies on probabilities of word co-occurrences to partition the set of terms appearing in the collection of documents into disjoint groups of related terms. In this work, we also present an environment to visualize the computed topics in the term space and retrieve the most related PubMed articles for each group of terms. We illustrate the algorithm by applying it to PubMed documents on the topic of suicide. Suicide is a major public health problem identified as the tenth leading cause of death in the US. In this application, our goal is to provide a global view of the mental health literature pertaining to the subject of suicide, and through this, to help create a rich environment of multifaceted data to guide health care researchers in their endeavor to better understand the breadth, depth and scope of the problem. We demonstrate the usefulness of the proposed algorithm by providing a web portal that allows mental health researchers to peruse the suicide-related literature in PubMed.


title: Enumerating the class of minimally path connected simplicial complexes

abstract: In 1983 Kalai proved an incredible generalisation of Cayley's formula for the number of trees on a labelled vertex set to a formula for a class of $r$-dimensional simplicial complexes. These simplicial complexes generalise trees by means of being homologically $\Bbb Q$-acyclic. In this text we consider a different generalisation of trees to the class of pure dimensional simplicial complexes that \textit{minimally connect a vertex set} (in the sense that the removal of any top dimensional face disconnects the complex). Our main result provides an upper and lower bound for the number of these minimally connected complexes on a labelled vertex set. We also prove that they are potentially vastly more topologically complex than the generalisation of Kalai. As an application of our bounds we compute the threshold probability for the connectivity of a random simplicial complex.


title: FocusNet++: Attentive Aggregated Transformations for Efficient and Accurate Medical Image Segmentation

abstract: We propose a new residual block for convolutional neural networks and demonstrate its state-of-the-art performance in medical image segmentation. We combine attention mechanisms with group convolutions to create our group attention mechanism, which forms the fundamental building block of our network, FocusNet++. We employ a hybrid loss based on balanced cross entropy, Tversky loss and the adaptive logarithmic loss to enhance the performance along with fast convergence. Our results show that FocusNet++ achieves state-of-the-art results across various benchmark metrics for the ISIC 2018 melanoma segmentation and the cell nuclei segmentation datasets with fewer parameters and FLOPs.


title: Steiner symmetrization for anisotropic quasilinear equations via partial discretization

abstract: In this paper we obtain comparison results for the quasilinear equation $-\Delta_{p,x} u - u_{yy} = f$ with homogeneous Dirichlet boundary conditions by Steiner rearrangement in variable $x$, thus solving a long open problem. In fact, we study a broader class of anisotropic problems. Our approach is based on a finite-differences discretization in $y$, and the proof of a comparison principle for the discrete version of the auxiliary problem $A U - U_{yy} \le \int_0^s f^*$, where $AU = (n\omega^{1/n}s^{1/n'} )^p (- U_{ss})^{p-1}$. We show that this operator is T-accretive in $L^\infty$. We extend our results for $-\Delta_{p,x}$ to general operators of the form $-\mathrm{div} (a(|\nabla_x u|) \nabla_x u)$ where $a$ is non-decreasing and behaves like $| \cdot |^{p-2}$ at infinity.


title: Gorenstein-projective modules over short local algebras

abstract: Following the well-established terminology in commutative algebra, any (not necessarily commutative) finite-dimensional local algebra $A$ with radical $J$ will be said to be short provided $J^3 = 0$. As in the commutative case, we show: if a short local algebra $A$ has an indecomposable non-projective Gorenstein-projective module $M$, then either $A$ is self-injective (so that all modules are Gorenstein-projective) and then $|J^2| \le 1$, or else $|J^2| = |J/J^2| - 1$ and $|JM| = |J^2||M/JM|$. More generally, we focus the attention to semi-Gorenstein-projective and $\infty$-torsionfree modules, even to $\mho$-paths of length 2, 3 and 4. In particular, we show that the existence of a non-projective reflexive module implies that $|J^2| < |J/J^2|$ and further restrictions. In addition, we consider exact complexes of projective modules with a non-projective image. Again, as in the commutative case, we see that if such a complex exists, then $A$ is self-injective or satisfies the condition $|J^2| = |J/J^2| - 1.$ Also, we show that any non-projective semi-Gorenstein-projective module $M$ satisfies $Ext^1(M,M) \neq 0$. In this way, we prove the Auslander-Reiten conjecture (one of the classical homological conjectures) for arbitrary short local algebras. Many arguments used in the commutative case actually work in general, but there are interesting differences and some of our results may be new also in the commutative case.


title: Periodic Homogenization of a L\'evy-Type Process with Small Jumps

abstract: In this article, we consider the problem of periodic homogenization of a Feller process generated by a pseudo-differential operator, the so-called L\'evy-type process. Under the assumptions that the generator has rapidly periodically oscillating coefficients, and that it admits "small jumps" only (that is, the jump kernel has finite second moment), we prove that the appropriately centered and scaled process converges weakly to a Brownian motion with covariance matrix given in terms of the coefficients of the generator. The presented results generalize the classical and well-known results related to periodic homogenization of a diffusion process.


title: Evaluating the Data Quality of Eye Tracking Signals from a Virtual Reality System: Case Study using SMI's Eye-Tracking HTC Vive

abstract: We evaluated the data quality of SMI's tethered eye-tracking head-mounted display based on the HTC Vive (ET-HMD) during a random saccade task. We measured spatial accuracy, spatial precision, temporal precision, linearity, and crosstalk. We proposed the use of a non-parametric spatial precision measure based on the median absolute deviation (MAD). Our linearity analysis considered both the slope and adjusted R-squared of a best-fitting line. We were the first to test for a quadratic component to crosstalk. We prepended a calibration task to the random saccade task and evaluated 2 methods to employ this user-supplied calibration. For this, we used a unique binning approach to choose samples to be included in the recalibration analyses. We compared our quality measures between the ET-HMD and our EyeLink 1000 (SR-Research, Ottawa, Ontario, CA). We found that the ET-HMD had significantly better spatial accuracy and linearity fit than our EyeLink, but both devices had similar spatial precision and linearity slope. We also found that, while the EyeLink had no significant crosstalk, the ET-HMD generally exhibited quadratic crosstalk. Fourier analysis revealed that the binocular signal was a low-pass filtered version of the monocular signal. Such filtering resulted in the binocular signal being useless for the study of high-frequency components such as saccade dynamics.


title: Mining Domain Knowledge: Improved Framework towards Automatically Standardizing Anatomical Structure Nomenclature in Radiotherapy

abstract: The automatic standardization of nomenclature for anatomical structures in radiotherapy (RT) clinical data is a critical prerequisite for data curation and data-driven research in the era of big data and artificial intelligence, but it is currently an unmet need. Existing methods either cannot handle cross-institutional datasets or suffer from heavy imbalance and poor-quality delineation in clinical RT datasets. To solve these problems, we propose an automated structure nomenclature standardization framework, 3D Non-local Network with Voting (3DNNV). This framework consists of an improved data processing strategy, namely, adaptive sampling and adaptive cropping (ASAC) with voting, and an optimized feature extraction module. The framework simulates clinicians' domain knowledge and recognition mechanisms to identify small-volume organs at risk (OARs) with heavily imbalanced data better than other methods. We used partial data from an open-source head-and-neck cancer dataset to train the model, then tested the model on three cross-institutional datasets to demonstrate its generalizability. 3DNNV outperformed the baseline model, achieving higher average true positive rates (TPR) overall categories on the three test datasets (+8.27%, +2.39%, and +5.53%, respectively). More importantly, the 3DNNV outperformed the baseline on the test dataset, 28.63% to 91.17%, in terms of F1 score for a small-volume OAR with only 9 training samples. The results show that 3DNNV can be applied to identify OARs, even error-prone ones. Furthermore, we discussed the limitations and applicability of the framework in practical scenarios. The framework we developed can assist in standardizing structure nomenclature to facilitate data-driven clinical research in cancer radiotherapy.


title: Protecting Geolocation Privacy of Photo Collections

abstract: People increasingly share personal information, including their photos and photo collections, on social media. This information, however, can compromise individual privacy, particularly as social media platforms use it to infer detailed models of user behavior, including tracking their location. We consider the specific issue of location privacy as potentially revealed by posting photo collections, which facilitate accurate geolocation with the help of deep learning methods even in the absence of geotags. One means to limit associated inadvertent geolocation privacy disclosure is by carefully pruning select photos from photo collections before these are posted publicly. We study this problem formally as a combinatorial optimization problem in the context of geolocation prediction facilitated by deep learning. We first demonstrate the complexity both by showing that a natural greedy algorithm can be arbitrarily bad and by proving that the problem is NP-Hard. We then exhibit an important tractable special case, as well as a more general approach based on mixed-integer linear programming. Through extensive experiments on real photo collections, we demonstrate that our approaches are indeed highly effective at preserving geolocation privacy.


title: Constraining the Milky Way Mass Profile with Phase-Space Distribution of Satellite Galaxies

abstract: We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations.


title: Third components with elliptical orbits in the eclipsing binaries

abstract: In our research, we studied 9 eclipsing binary stars: AR Lac, U CrB, S Equ, SU Boo, VV UMa, WW Gem, YY Eri, V0404 Lyr, HP Aur. We collected large sets of moments of minima from BRNO and observational data from AAVSO databases. Then, we obtained moments of minima for AAVSO observations (totally - 397 minima) using method of approximation with symmetric polynomial, realized in software MAVKA. This software was provided by Kateryna D. Andrych and Ivan L. Andronov. After we combined obtained moments of minima with data taken from BRNO and plotted O-C diagrams. For all stars these diagrams represented sinusoidal-like oscillations with superposition of linear (for SU Boo, VV UMa, WW Gem, V0404 Lyr and HP Aur) or parabolic trend (for AR Lac, U CrB, S Equ and YY Eri). The oscillations could be described as presence of the third component, which does not take part in eclipses, but causes the well-known light-time effect (LTE). That effect could be easily detected using long data series of observations. However, the oscillations have clear asymmetry, which we interpreted as elliptical shape of the third components orbit. Parabolic trend we explained as mass transfer between components of binary system. For all stars we computed minimal possible mass of the third component. In addition, we developed our own code in computed language Python and using it we computed orbital elements of the third component. Moreover, for stars with parabolic trend we calculated rate of the mass transfer. Finally, for all computed values we estimated errors.


title: Heavy-Light Mesons on the Light Front

abstract: We study the heavy-light mesons within basis light-front quantization. The resulting mass spectra of $D$, $D_s$, $B$, and $B_s$ agree reasonably well with experiments. We also predict states which could be measured in the near future. In the light-front formalism, we calculate the light-front wave functions and additional experimental observables, such as parton distribution functions, distribution amplitudes, and decay constants by means of integrations over light-front wave functions. We also provide ratios of decay constants for selected pseudoscalar meson decays ($D_s$ to $D$ and $B_s$ to $B$) as they may prove to be theoretically more robust and more reliably determined in experiments. We find that our ratios are systematically smaller than existing experiment and other approaches by $5-18\%$.


title: On the Security of A Remote Cloud Storage Integrity Checking Protocol

abstract: Data security and privacy is an important but challenging problem in cloud computing. One of the security concerns from cloud users is how to efficiently verify the integrity of their data stored on the cloud server. Third Party Auditing (TPA) is a new technique proposed in recent years to achieve this goal. In a recent paper (IEEE Transactions on Computers 62(2): 362-375 (2013)), Wang et al. proposed a highly efficient and scalable TPA protocol and also a Zero Knowledge Public Auditing protocol which can prevent offline guessing attacks. However, in this paper, we point out several security weaknesses in Wang et al's protocols: first, we show that an attacker can arbitrarily modify the cloud data without being detected by the auditor in the integrity checking process, and the attacker can achieve this goal even without knowing the content of the cloud data or any verification metadata maintained by the cloud server; secondly, we show that the Zero Knowledge Public Auditing protocol cannot achieve its design goal, that is to prevent offline guessing attacks.


title: Diffeological statistical models, the Fisher metric and probabilistic mappings

abstract: In this note we introduce the notion of a $C^k$-diffeological statistical model, which allows us to apply the theory of diffeological spaces to (possibly singular) statistical models. In particular, we introduce a class of almost 2-integrable $C^k$-diffeological statistical models that encompasses all known statistical models for which the Fisher metric is defined. This class contains a statistical model which does not appear in the Ay-Jost-L\^e-Schwachh\"ofer theory of parametrized measure models. Then we show that for any positive integer $k$ the class of almost 2-integrable $C^k$-diffeological statistical models is preserved under probabilistic mappings. Furthermore, the monotonicity theorem for the Fisher metric also holds for this class. As a consequence, the Fisher metric on an almost 2-integrable $C^k$-diffeological statistical model $P \subset {\cal P}({\cal X})$ is preserved under any probabilistic mapping $T: {\cal X}\leadsto {\cal Y}$ that is sufficient w.r.t. $P$. Finally we extend the Cram\'er-Rao inequality to the class of 2-integrable $C^k$-diffeological statistical models.


title: Applications of resonance theory without analyticity assumption

abstract: We prove that the results in scattering theory that involve resonances are still valid for non-analytic potentials, even if the notion of resonance is not defined in this setting. More precisely, we show that if the potential of a semiclassical Schr\"odinger operator is supposed to be smooth and to decrease at infinity, the usual formulas relating scattering quantities and resonances still hold. The main ingredient for the proofs is a resolvent estimate of a new type, relating the resolvent of an operator with the resolvent of its cut-off counterpart.


title: Investigating the signatures of long-range persistence in seismic sequences along Circum-Pacific subduction zones

abstract: In the present paper, we analyze the signatures of long-range persistence in seismic sequences along Circum-Pacific subduction zones, from Chile to Kermadec, extracted from the National Earthquake Information Center (NEIC) catalog. This region, known as the Pacific Ring of Fire, is the world's most active fault line, containing about 90$\%$ of the world's earthquakes. We used the classical rescaled range ($R/S$) analysis to estimate the long-term persistence signals derived from a scaling parameter called the Hurst exponent, $H$. We measured the referred exponent and obtained values of $H>0.5$, indicating that a long-term memory effect exists. We found a possible fractal relationship between $H$ and the $b_{s}(q)$-index, which emerges from the non-extensive Gutenberg-Richter law as a function of the asperity. Therefore, $H$ can be associated with a mechanism that controls the level of seismic activity. Finally, we concluded that the dynamics associated with fragment-asperity interactions can be classified as a self-affine fractal phenomenon.


title: Implementing a smooth exact penalty function for general constrained nonlinear optimization

abstract: We build upon Estrin et al. (2019) to develop a general constrained nonlinear optimization algorithm based on a smooth penalty function proposed by Fletcher (1970, 1973b). Although Fletcher's approach has historically been considered impractical, we show that the computational kernels required are no more expensive than those in other widely accepted methods for nonlinear optimization. The main kernel for evaluating the penalty function and its derivatives solves structured linear systems. When the matrices are available explicitly, we store a single factorization each iteration. Otherwise, we obtain a factorization-free optimization algorithm by solving each linear system iteratively. The penalty function shows promise in cases where the linear systems can be solved efficiently, e.g., PDE-constrained optimization problems when efficient preconditioners exist. We demonstrate the merits of the approach, and give numerical results on several PDE-constrained and standard test problems.


title: A Step Towards Exposing Bias in Trained Deep Convolutional Neural Network Models

abstract: We present Smooth Grad-CAM++, a technique which combines two recent techniques: SMOOTHGRAD and Grad-CAM++. Smooth Grad-CAM++ has the capability of either visualizing a layer, subset of feature maps, or subset of neurons within a feature map at each instance. We experimented with few images, and we discovered that Smooth Grad-CAM++ produced more visually sharp maps with larger number of salient pixels highlighted in the given input images when compared with other methods. Smooth Grad-CAM++ will give insight into what our deep CNN models (including models trained on medical scan or imagery) learn. Hence informing decisions on creating a representative training set.


title: Open source QGIS toolkit for the Advanced Research WRF modelling system

abstract: The Advanced Research WRF (Weather Research and Forecasting) model is a popular atmospheric model used for research and Numerical Weather Prediction (NWP). However, despite its popularity, its set-up and configuration often demand several interdisciplinary skills that go beyond the understanding of physical processes. Pre-processing tasks, such as importing custom high-resolution datasets in the WRF Pre-processing System (WPS), still require a considerable effort from the user. We present GIS4WRF, a free, open-source, and cross-platform QGIS Python plug-in to help scientists and practitioners with their Advanced Research WRF modelling workflows. GIS4WRF incorporates new and existing tools for data-processing, configuration, simulation, and visualization into a single graphical environment, and offers WRF-CMake binary distributions for Windows, macOS, and Linux. We highlight its main features and provide useful insights into several key approaches and techniques used in its development. We end with two example applications highlighting the contributions of GIS4WRF in simplifying several WRF-related tasks.


title: Learning Multi-Object Tracking and Segmentation from Automatic Annotations

abstract: In this work we contribute a novel pipeline to automatically generate training data, and to improve over state-of-the-art multi-object tracking and segmentation (MOTS) methods. Our proposed track mining algorithm turns raw street-level videos into high-fidelity MOTS training data, is scalable and overcomes the need of expensive and time-consuming manual annotation approaches. We leverage state-of-the-art instance segmentation results in combination with optical flow predictions, also trained on automatically harvested training data. Our second major contribution is MOTSNet - a deep learning, tracking-by-detection architecture for MOTS - deploying a novel mask-pooling layer for improved object association over time. Training MOTSNet with our automatically extracted data leads to significantly improved sMOTSA scores on the novel KITTI MOTS dataset (+1.9%/+7.5% on cars/pedestrians), and MOTSNet improves by +4.1% over previously best methods on the MOTSChallenge dataset. Our most impressive finding is that we can improve over previous best-performing works, even in complete absence of manually annotated MOTS training data.


title: Optimal Green Hybrid Attacks in Secure IoT

abstract: In pursuance to understand the behavior of a potential green hybrid attacker in secure internet-of-things (IoT), this letter investigates optimal energy utilization from the attacker's viewpoint. Specifically, we propose a novel framework for optimizing the efficacy of a hybrid attacker, possessing capability to both eavesdrop and jam, while considering the underlying energy consumption. In particular, we maximize the attacker energy efficiency (AEE) in secure IoT by deriving the analytical solutions for jointly global-optimal attacking mode, eavesdropping rate, and jamming power. Numerical results, validating analytical claims, reveal that the proposed green design can offer about 45% improvement in average AEE over the relevant benchmarks.


title: Expressiveness and Learning of Hidden Quantum Markov Models

abstract: Extending classical probabilistic reasoning using the quantum mechanical view of probability has been of recent interest, particularly in the development of hidden quantum Markov models (HQMMs) to model stochastic processes. However, there has been little progress in characterizing the expressiveness of such models and learning them from data. We tackle these problems by showing that HQMMs are a special subclass of the general class of observable operator models (OOMs) that do not suffer from the \emph{negative probability problem} by design. We also provide a feasible retraction-based learning algorithm for HQMMs using constrained gradient descent on the Stiefel manifold of model parameters. We demonstrate that this approach is faster and scales to larger models than previous learning algorithms.


title: Energy upper bound for structurally stable N-passive states

abstract: Passive states are special configurations of a quantum system which exhibit no energy decrement at the end of an arbitrary cyclic driving of the model Hamiltonian. When applied to an increasing number of copies of the initial density matrix, the requirement of passivity induces a hierarchical ordering which, in the asymptotic limit of infinitely many elements, pinpoints ground states and thermal Gibbs states. In particular, for large values of N the energy content of a N-passive state which is also structurally stable (i.e. capable to maintain its passivity status under small perturbations of the model Hamiltonian), is expected to be close to the corresponding value of the thermal Gibbs state which has the same entropy. In the present paper we provide a quantitative assessment of this fact, by producing an upper bound for the energy of an arbitrary N-passive, structurally stable state which only depends on the spectral properties of the Hamiltonian of the system. We also show the condition under which our inequality can be saturated. A generalization of the bound is finally presented that, for sufficiently large N, applies to states which are N-passive, but not necessarily structurally stable.

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


title: Current-Generating Mechanism for El Ni\~no & La Ni\~na, an Data Evidence from Integrated ESMD Method

abstract: Different to classical theories which emphasize east--west feedback along the equator, data analyses via the extreme-point symmetric mode decomposition (ESMD) method support a new viewpoint that the El Ni\~{n}o--Southern Oscillation (ENSO) cycle is horizontally dependent, rather than vertically dependent. The consistency between the ocean current changes and sea-surface temperature anomalies results in a current-generating mechanism for El Ni\~{n}o and La Ni\~{n}a events. The corresponding new findings are as follows: (1) The appearance of the Central Equatorial Countercurrent (CECC) and the enhancement of the North Equatorial Countercurrent (NECC) are primary causes for the occurrences of positive anomalies in the central and eastern Pacific, respectively. (2) In addition to the contribution from the CECC, initiations of the eastern and central Pacific types of El Ni\~{n}o events are mainly caused by the southward transport of warm water from the high-temperature sea-area off Central America, and from the eastern low-speed gap of the NECC, respectively. (3) The northern and southern branches of the South Equatorial Current contribute to the negative anomalies in the eastern and central Pacific, respectively. Particularly, the later transfers cold water from the southeast to the central equatorial Pacific in a direct way, and this implies the crack of the ENSO-asymmetry problem.


title: The Zeeman, Spin-Orbit, and Quantum Spin-Hall Interactions in Anisotropic and Low-Dimensional Conductors

abstract: When an electron or hole is in a conduction band of a crystal, it can be very different from 2, depending upon the crystalline anisotropy and the direction of the applied magnetic induction ${\bf B}$. In fact, it can even be 0! To demonstrate this quantitatively, the Dirac equation is extended for a relativistic electron or hole in an orthorhombically-anisotropic conduction band with effective masses $m_j$ for $j=1,2,3$ with geometric mean $m_g=(m_1m_2m_3)^{1/3}$. The appropriate Foldy-Wouthuysen transformations are extended to evaluate the non-relativistic Hamiltonian to $O({\rm m}c^2)^{-4}$, where ${\rm m}c^2$ is the particle's Einstein rest energy. For ${\bf B}||\hat{\bf e}_{\mu}$, the Zeeman $g_{\mu}$ factor is $2{\rm m}\sqrt{m_{\mu}}/m_g^{3/2} + O({\rm m}c^2)^{-2}$. While propagating in a two-dimensional (2D) conduction band with $m_3\gg m_1,m_2$, $g_{||}<<2$, consistent with recent measurements of the temperature $T$ dependence of the parallel upper critical induction $B_{c2,||}(T)$ in superconducting monolayer NbSe$_2$ and in twisted bilayer graphene. While a particle is in its conduction band of an atomically thin one-dimensional metallic chain along $\hat{\bf e}_{\mu}$, $g<<2$ for all ${\bf B}={\bf\nabla}\times{\bf A}$ directions and vanishingly small for ${\bf B}||\hat{\bf e}_{\mu}$. The quantum spin Hall Hamiltonian for 2D metals with $m_1=m_2=m_{||}$ is $K[{\bf E}\times({\bf p}-q{\bf A})]_{\perp}\sigma_{\perp}+O({\rm m}c^2)^{-4}$, where ${\bf E}$ and ${\bf p}-q{\bf A}$ are the planar electric field and gauge-invariant momentum, $q=\mp|e|$ is the particle's charge, $\sigma_{\perp}$ is the Pauli matrix normal to the layer, $K=\pm\mu_B/(2m_{||}c^2)$, and $\mu_B$ is the Bohr magneton.


title: Artificial Intelligence for Low-Resource Communities: Influence Maximization in an Uncertain World

abstract: The potential of Artificial Intelligence (AI) to tackle challenging problems that afflict society is enormous, particularly in the areas of healthcare, conservation and public safety and security. Many problems in these domains involve harnessing social networks of under-served communities to enable positive change, e.g., using social networks of homeless youth to raise awareness about Human Immunodeficiency Virus (HIV) and other STDs. Unfortunately, most of these real-world problems are characterized by uncertainties about social network structure and influence models, and previous research in AI fails to sufficiently address these uncertainties. This thesis addresses these shortcomings by advancing the state-of-the-art to a new generation of algorithms for interventions in social networks. In particular, this thesis describes the design and development of new influence maximization algorithms which can handle various uncertainties that commonly exist in real-world social networks. These algorithms utilize techniques from sequential planning problems and social network theory to develop new kinds of AI algorithms. Further, this thesis also demonstrates the real-world impact of these algorithms by describing their deployment in three pilot studies to spread awareness about HIV among actual homeless youth in Los Angeles. This represents one of the first-ever deployments of computer science based influence maximization algorithms in this domain. Our results show that our AI algorithms improved upon the state-of-the-art by 160% in the real-world. We discuss research and implementation challenges faced in deploying these algorithms, and lessons that can be gleaned for future deployment of such algorithms. The positive results from these deployments illustrate the enormous potential of AI in addressing societally relevant problems.


title: Examples of metric spaces with asymptotic property $C$

abstract: We construct a class of metric spaces whose transfinite asymptotic dimension and complementary-finite asymptotic dimension are both $\omega+k$ for any $k\in\mathbb{N}$, where $\omega$ is the smallest infinite ordinal number and a metric space whose transfinite asymptotic dimension and complementary-finite asymptotic dimension are both $2\omega$. Moreover, we study the relationship between asymptotic dimension growth, transfinite asymptotic dimension and finite decomposition complexity.


title: $(\phi,\tau)$-modules diff\'erentiels et repr\'esentations potentiellement semi-stables

abstract: Soit $K$ un corps $p$-adique et soit $V$ une repr\'esentation $p$-adique de $\mathcal{G}_K = \mathrm{Gal}(\bar{K}/K)$. La surconvergence des $(\phi,\tau)$-modules nous permet d'attacher \`a $V$ un $\phi$-module diff\'erentiel \`a connexion $D_{\tau,\mathrm{rig}}^\dagger(V)$ sur l'anneau de Robba $\mathbf{B}_{\tau,\mathrm{rig},K}^\dagger$. On montre dans cet article comment retrouver les invariants $D_{\mathrm{cris}}(V)$ et $D_{\mathrm{st}}(V)$ \`a partir de $D_{\tau,\mathrm{rig}}^\dagger(V)$, et comment caract\'eriser les repr\'esentations potentiellement semi-stables, ainsi que celles de $E$-hauteur finie, \`a partir de la connexion. Let $K$ be a $p$-adic field and let $V$ be a $p$-adic representation of $\mathcal{G}_K=\mathrm{Gal}(\bar{K}/K)$. The overconvergence of $(\phi,\tau)$-modules allows us to attach to $V$ a differential $\phi$-module $D_{\tau,\mathrm{rig}}^\dagger(V)$ on the Robba ring $\mathbf{B}_{\tau,\mathrm{rig},K}^\dagger$ that comes equipped with a connection. We show in this paper how to recover the invariants $D_{\mathrm{cris}}(V)$ and $D_{\mathrm{st}}(V)$ from $D_{\tau,\mathrm{rig}}^\dagger(V)$, and give a characterization of both potentially semi-stable representations of $\mathcal{G}_K$ and finite $E$-height representations in terms of the connection operator.


title: Influence Maximization for Social Good: Use of Social Networks in Low Resource Communities

abstract: This thesis proposal makes the following technical contributions: (i) we provide a definition of the Dynamic Influence Maximization Under Uncertainty (or DIME) problem, which models the problem faced by homeless shelters accurately; (ii) we propose a novel Partially Observable Markov Decision Process (POMDP) model for solving the DIME problem; (iii) we design two scalable POMDP algorithms (PSINET and HEALER) for solving the DIME problem, since conventional POMDP solvers fail to scale up to sizes of interest; and (iv) we test our algorithms effectiveness in the real world by conducting a pilot study with actual homeless youth in Los Angeles. The success of this pilot (as explained later) shows the promise of using influence maximization for social good on a larger scale.


title: High-precision limits on $W$-$W'$ and $Z$-$Z'$ mixing from diboson production using the full LHC Run 2 ATLAS data set

abstract: The full ATLAS Run 2 data set with time-integrated luminosity of 139 fb$^{-1}$ in the diboson channels in hadronic final states is used to probe a simple model with an extended gauge sector (EGM), proposed by Altarelli et al., and often taken as a convenient benchmark by experimentalists. This model accommodates new charged $W'$ and neutral $Z'$ vector bosons with modified trilinear Standard Model gauge couplings, decaying into electroweak gauge boson pairs $WZ$ or $WW$, where $W$/$Z$ decay hadronically. Exclusion limits at the 95% C.L. on the $Z'$ and $W'$ resonance production cross section times branching ratio to electroweak gauge boson pairs in the mass range of $\sim$ 1 - 5 TeV are here converted to constraints on $W$-$W'$ and $Z$-$Z'$ mixing parameters and masses for the EGM. We present exclusion regions on the parameter space of the $W'$ and $Z'$ by using the full Run 2 data set comprised of $pp$ collisions at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector at the CERN LHC. The obtained exclusion regions are significantly extended compared to those obtained from the previous analysis performed with Tevatron data as well as with LHC data collected at 7 and 8 TeV in Run 1 and are the most stringent bounds to date.


title: Exact solution of a anisotropic $J_1-J_2$ spin chain with antiperiodic boundary condition

abstract: The exact solution of an integrable anisotropic Heisenberg spin chain with nearest-neighbour, next-nearest-neighbour and scalar chirality couplings is studied, where the boundary condition is the antiperiodic one. The detailed construction of Hamiltonian and the proof of integrability are given. The antiperiodic boundary condition breaks the $U(1)$-symmetry of the system and we use the off-diagonal Bethe Ansatz to solve it. The energy spectrum is characterized by the inhomogeneous $T-Q$ relations and the contribution of the inhomogeneous term is studied. The ground state energy and the twisted boundary energy in different regions are obtained. We also find that the Bethe roots at the ground state form the string structure if the coupling constant $J=-1$ although the Bethe Ansatz equations are the inhomogeneous ones.


title: Multivariate scale-mixed stable distributions and related limit theorems

abstract: In the paper, multivariate probability distributions are considered that are representable as scale mixtures of multivariate elliptically contoured stable distributions. It is demonstrated that these distributions form a special subclass of scale mixtures of multivariate elliptically contoured normal distributions. Some properties of these distributions are discussed. Main attention is paid to the representations of the corresponding random vectors as products of independent random variables. In these products, relations are traced of the distributions of the involved terms with popular probability distributions. As examples of distributions of the class of scale mixtures of multivariate elliptically contoured stable distributions, multivariate generalized Linnik distributions are considered in detail. Their relations with multivariate `ordinary' Linnik distributions, multivariate normal, stable and Laplace laws as well as with univariate Mittag-Leffler and generalized Mittag-Leffler distributions are discussed. Limit theorems are proved presenting necessary and sufficient conditions for the convergence of the distributions of random sequences with independent random indices (including sums of a random number of random vectors and multivariate statistics constructed from samples with random sizes) to scale mixtures of multivariate elliptically contoured stable distributions. The property of scale-mixed multivariate stable distributions to be both scale mixtures of a non-trivial multivariate stable distribution and a normal scale mixture is used to obtain necessary and sufficient conditions for the convergence of the distributions of random sums of random vectors with both infinite or finite covariance matrices to the multivariate generalized Linnik distribution.


title: Quantifying Urban Canopy Cover with Deep Convolutional Neural Networks

abstract: Urban canopy cover is important to mitigate the impact of climate change. Yet, existing quantification of urban greenery is either manual and not scalable, or use traditional computer vision methods that are inaccurate. We train deep convolutional neural networks (DCNNs) on datasets used for self-driving cars to estimate urban greenery instead, and find that our semantic segmentation and direct end-to-end estimation method are more accurate and scalable, reducing mean absolute error of estimating the Green View Index (GVI) metric from 10.1% to 4.67%. With the revised DCNN methods, the Treepedia project was able to scale and analyze canopy cover in 22 cities internationally, sparking interest and action in public policy and research fields.


title: Isospin amplitudes in $\Lambda_b^0\to J/\psi \Lambda(\Sigma^0)$ and $\Xi_b^0\to J/\psi \Xi^0(\Lambda)$ decays

abstract: Ratios of isospin amplitudes in hadron decays are a useful probe of the interplay between weak and strong interactions, and allow searches for physics beyond the Standard Model. We present the first results on isospin amplitudes in $b$-baryon decays, using data corresponding to an integrated luminosity of 8.5 fb$^{-1}$, collected with the LHCb detector in $pp$ collisions at center of mass energies of 7, 8 and 13 TeV. The isospin amplitude ratio $|A_1(\Lambda_b^0\to J/\psi \Sigma^0)/A_0(\Lambda_b^0\to J/\psi\Lambda)|$, where the subscript on $A$ indicates the final-state isospin, is measured to be less than 1/21.8 at 95\% confidence level. The Cabibbo suppressed $\Xi_b^0\to J/\psi\Lambda$ decay is observed for the first time, allowing for the measurement $|A_0(\Xi_b^0\to J/\psi\Lambda)/A_{1/2}(\Xi_b^0\to J/\psi\Xi^0)| =0.37 \pm 0.06\pm 0.02$, where the uncertainties are statistical and systematic, respectively.


title: Recollements induced by good silting objects

abstract: Let $U$ be a silting object in a derived category over a dg-algebra $A$, and let $B$ be the endomorphism dg-algebra of $U$. Under some appropriate hypotheses, we show that if $U$ is good, then there exist a dg-algebra $C$, a homological epimorphism $B\rightarrow C$ and a recollement among the (unbounded) derived categories $\mathbf{D}(C,d)$ of $C$, $\mathbf{D}(B,d)$ of $B$ and $\mathbf{D}(A,d)$ of $A$. In particular, the kernel of the left derived functor $-\otimes^{\mathbb{L}}_{B}U$ is triangle equivalent to the derived category $\mathbf{D}(C,d)$. Conversely, if $-\otimes^{\mathbb{L}}_{B}U$ admits a fully faithful left adjoint functor, then $U$ is good. Moreover, we establish a criterion for the existence of a recollement of the derived category of a dg-algebra relative to two derived categories of weak non-positive dg-algebras. Finally, some applications are given related to good cosilting objects, good 2-term silting complexes, good tilting complexes and modules, which recovers a recent result by Chen and Xi.


title: On value sets of fractional ideals

abstract: The aim of this work is to study duality of fractional ideals with respect to a fixed ideal and to investigate the relationship between value sets of pairs of dual ideals in admissible rings, a class of rings that contains the local rings of algebraic curves at singular points. We characterize canonical ideals by means of a symmetry relation between lengths of certain quotients of associated ideals to a pair of dual ideals. In particular, we extend the symmetry among absolute and relative maximals in the sets of values of pairs of dual fractional ideals to other kinds of maximal points. Our results generalize and complement previous ones by other authors.


title: Diagnosing pulsar winds in black-widow, redback, and other binary millisecond pulsar systems

abstract: Binary systems composed of a recycled millisecond pulsar and a stellar companion in close orbit could be excellent sites to diagnose pulsar winds. In such systems, the pulsar outflow irradiates and heats up the companion atmosphere, which can lead to the observation of strong day/night modulations in temperature. We demonstrate with particle shower simulations that the particle energy of the wind affects the heating depth in the atmosphere: the wind heat can be deposited above or below the photosphere, leading to different signatures in the observed spectra. We apply our method to four specific systems: We find that systems with cool night side companions showing strong temperature variations can give interesting lower limits on the particle energy in the winds. In particular, if the companion night side of PSR B1957+20 were to be suddenly irradiated, deep heating would only take place if particles with energy > 100 TeV were present. Observational evidence of deep heating in this system thus suggests that i) such particles exist in the pulsar wind and/or ii) binary evolution non-trivially takes the companion to the observed temperature asymmetry. Besides, the observed temperature difference can be maintained only with particle energies of the order of 100 MeV.


title: Keyword Aware Influential Community Search in Large Attributed Graphs

abstract: We introduce a novel keyword-aware influential community query KICQ that finds the most influential communities from an attributed graph, where an influential community is defined as a closely connected group of vertices having some dominance over other groups of vertices with the expertise (a set of keywords) matching with the query terms (words or phrases). We first design the KICQ that facilitates users to issue an influential CS query intuitively by using a set of query terms, and predicates (AND or OR). In this context, we propose a novel word-embedding based similarity model that enables semantic community search, which substantially alleviates the limitations of exact keyword based community search. Next, we propose a new influence measure for a community that considers both the cohesiveness and influence of the community and eliminates the need for specifying values of internal parameters of a network. Finally, we propose two efficient algorithms for searching influential communities in large attributed graphs. We present detailed experiments and a case study to demonstrate the effectiveness and efficiency of the proposed approaches.


title: Heavy quark flow as better probes of QGP properties

abstract: In earlier studies we have proposed that most parton $v_2$ comes from the anisotropic escape of partons, not from the hydrodynamic flow, even for semi-central Au+Au collisions at $\sqrt {s_{NN}}=200$ GeV. Here we study the flavor dependence of this escape mechanism with a multi-phase transport model. In contrast to naive expectations, we find that the charm $v_2$ is much more sensitive to the hydrodynamic flow than the lighter quark $v_2$, and the fraction of $v_2$ from the escape mechanism decreases strongly with the quark mass for large collision systems. We also find that the light quark collective flow is essential for the charm quark $v_2$. Our finding thus suggests that heavy quark flows are better probes of the quark-gluon-plasma properties than light quark flows.


title: Rapid technique for preparing conductive graphene-borosilicone putty

abstract: When borosilicone putty is filled with graphene sheets, the result can be a "g-putty"--an electrically conducting putty that is exquisitely sensitive to deformation. In this work, we present a rapid method for producing g-putty from commercially available ingredients. The process can be performed at room temperature and has methanol as its only waste product.


title: Dust masses and grain size distributions of a sample of Galactic pulsar wind nebulae

abstract: We calculate dust spectral energy distributions (SEDs) for a range of grain sizes and compositions, using physical properties appropriate for five pulsar wind nebulae (PWNe) from which dust emission associated with the ejecta has been detected. By fitting the observed dust SED with our models, with the number of grains of different sizes as the free parameters, we are able to determine the grain size distribution and total dust mass in each PWN. We find that all five PWNe require large ($\ge 0.1 \, {\rm \mu m}$) grains to make up the majority of the dust mass, with strong evidence for the presence of micron-sized or larger grains. Only two PWNe contain non-negligible quantities of small ($<0.01 \, {\rm \mu m}$) grains. The size distributions are generally well-represented by broken power laws, although our uncertainties are too large to rule out alternative shapes. We find a total dust mass of $0.02-0.28 \, {\rm M}_\odot$ for the Crab Nebula, depending on the composition and distance from the synchrotron source, in agreement with recent estimates. For three objects in our sample, the PWN synchrotron luminosity is insufficient to power the observed dust emission, and additional collisional heating is required, either from warm, dense gas as found in the Crab Nebula, or higher temperature shocked material. For G$54.1$+$0.3$, the dust is heated by nearby OB stars rather than the PWN. Inferred dust masses vary significantly depending on the details of the assumed heating mechanism, but in all cases large mass fractions of micron-sized grains are required.


title: Thermodynamic Evidence of Proximity to a Kitaev Spin-Liquid in Ag$_{3}$LiIr$_{2}$O$_{6}$

abstract: Kitaev magnets are materials with bond-dependent Ising interactions between localized spins on a honeycomb lattice. Such interactions could lead to a quantum spin-liquid (QSL) ground state at zero temperature. Recent theoretical studies suggest two potential signatures of a QSL at finite temperatures, namely a scaling behavior of thermodynamic quantities in the presence of quenched disorder, and a two-step release of the magnetic entropy. Here, we present both signatures in Ag$_{3}$LiIr$_{2}$O$_{6}$ which is synthesized from $\alpha$-Li$_{2}$IrO$_{3}$ by replacing the inter-layer Li atoms with Ag atoms. In addition, the DC susceptibility data confirm absence of a long-range order, and the AC susceptibility data rule out a spin-glass transition. These observations suggest a closer proximity to the QSL in Ag$_{3}$LiIr$_{2}$O$_{6}$ compared to its parent compound $\alpha$-Li$_{2}$IrO$_{3}$ that orders at 15 K. We discuss an enhanced spin-orbit coupling due to a mixing between silver d and oxygen p orbitals as a potential underlying mechanism.


title: A Path Towards Quantum Advantage in Training Deep Generative Models with Quantum Annealers

abstract: The development of quantum-classical hybrid (QCH) algorithms is critical to achieve state-of-the-art computational models. A QCH variational autoencoder (QVAE) was introduced in Ref. [1] by some of the authors of this paper. QVAE consists of a classical auto-encoding structure realized by traditional deep neural networks to perform inference to, and generation from, a discrete latent space. The latent generative process is formalized as thermal sampling from either a quantum or classical Boltzmann machine (QBM or BM). This setup allows quantum-assisted training of deep generative models by physically simulating the generative process with quantum annealers. In this paper, we have successfully employed D-Wave quantum annealers as Boltzmann samplers to perform quantum-assisted, end-to-end training of QVAE. The hybrid structure of QVAE allows us to deploy current-generation quantum annealers in QCH generative models to achieve competitive performance on datasets such as MNIST. The results presented in this paper suggest that commercially available quantum annealers can be deployed, in conjunction with well-crafted classical deep neutral networks, to achieve competitive results in unsupervised and semisupervised tasks on large-scale datasets. We also provide evidence that our setup is able to exploit large latent-space (Q)BMs, which develop slowly mixing modes. This expressive latent space results in slow and inefficient classical sampling, and paves the way to achieve quantum advantage with quantum annealing in realistic sampling applications.


title: Improved Grid Optimization and Fitting in Least Squares Tensor Hypercontraction

abstract: A new method for generating fitting grids for least-squares tensor hypercontraction (LS-THC) is presented. This method draws inspiration from the related interpolative separable density fitting (ISDF) technique, but uses only a pivoted Cholesky decomposition of the metric matrix, S, already computed as a matter of course in LS-THC. The size and quality of the resulting grid is controlled by a user-defined cutoff parameter and the size of the starting grid. Additionally, the Cholesky-based method provides an alternative and possible more numerically stable method for performing the least-squares fit. The quality of the grids produced is evaluated for LS-DF-THC-MP2 calculations on retinal and benzene, the former with a large starting grid and small cc-pVDZ basis set, and the latter with a wide range of grids and basis sets. The error and grid size is found to be well-controlled by either the cutoff parameter (with a large starting grid) or the starting grid size (with a tight cutoff) and highly predictable. The Cholesky-based method is also able to generate unique grids tailored to different charge distributions, for example the (ab|, (ai|, and (ij| distributions that arise in the molecular orbital integrals. While only the (ai| grid directly affects the MP2 energy, the relative sizes of the other grids are examined.


title: Light-Driven Raman Coherence as a Non-Thermal Route to Ultrafast Topology Switching

abstract: A grand challenge underlies the entire field of topology-enabled quantum logic and information science: how to establish topological control principles driven by quantum coherence and understand the time-dependence of such periodic driving? Here we demonstrate a THz pulse-induced phase transition in Dirac materials that is periodically driven by vibrational coherence due to excitation of the lowest Raman-active mode. Above a critical field threshold, there emerges a long-lived metastable phase with unique Raman coherent phonon-assisted switching dynamics, absent for optical pumping. The switching also manifest itself by non-thermal spectral shape, relaxation slowing down near the Lifshitz transition where the critical Dirac point (DP) occurs, and diminishing signals at the same temperature that the Berry curvature induced Anomalous Hall Effect varnishes. These results, together with first-principles modeling, identify a mode-selective Raman coupling that drives the system from strong to weak topological insulators, STI to WTI, with a Dirac semimetal phase established at a critical atomic displacement controlled by the phonon pumping. Harnessing of vibrational coherence can be extended to steer symmetry-breaking transitions, i.e., Dirac to Weyl ones, with implications on THz topological quantum gate and error correction applications.


title: PRESAS: Block-Structured Preconditioning of Iterative Solvers within a Primal Active-Set Method for fast MPC

abstract: Model predictive control (MPC) for linear dynamical systems requires solving an optimal control structured quadratic program (QP) at each sampling instant. This paper proposes a primal active-set strategy (PRESAS) for the efficient solution of such block-sparse QPs, based on a preconditioned iterative solver to compute the search direction in each iteration. Rank-one factorization updates of the preconditioner result in a per-iteration computational complexity of $\mathcal{O}(N m^2)$, where $m$ denotes the number of state and control variables and $N$ the number of control intervals. Three different block-structured preconditioning techniques are presented and their numerical properties are studied further. In addition, an augmented Lagrangian based implementation is proposed to avoid a costly initialization procedure to find a primal feasible starting point. Based on a standalone C code implementation, we illustrate the computational performance of PRESAS against current state of the art QP solvers for multiple linear and nonlinear MPC case studies. We also show that the solver is real-time feasible on a dSPACE MicroAutoBox-II rapid prototyping unit for vehicle control applications, and numerical reliability is illustrated based on experimental results from a testbench of small-scale autonomous vehicles.


title: Methods in Higher Forcing Axioms (Workshop Notes)

abstract: Methods of Higher Forcing Axioms was a small workshop in Norwich, taking place between 10--12 of September, 2019. The goal was to encourage future collaborations, and create more focused threads of research on the topic of higher forcing axioms. This is an improved version of the notes taken during the meeting by Asaf Karagila.


title: Characterizing decoherence rates of a superconducting qubit by direct microwave scattering

abstract: We experimentally investigate a superconducting qubit coupled to the end of an open transmission line, in a regime where the qubit decay rates to the transmission line and to its own environment are comparable. We perform measurements of coherent and incoherent scattering, on- and off-resonant fluorescence, and time-resolved dynamics to determine the decay and decoherence rates of the qubit. In particular, these measurements let us discriminate between non-radiative decay and pure dephasing. We combine and contrast results across all methods and find consistent values for the extracted rates. The results show that the pure dephasing rate is one order of magnitude smaller than the non-radiative decay rate for our qubit. Our results indicate a pathway to benchmark decoherence rates of superconducting qubits in a resonator-free setting.


title: Towards Sustainable Architecture: 3D Convolutional Neural Networks for Computational Fluid Dynamics Simulation and Reverse DesignWorkflow

abstract: We present a general and flexible approximation model for near real-time prediction of steady turbulent flow in a 3D domain based on residual Convolutional Neural Networks (CNNs). This approach can provide immediate feedback for real-time iterations at the early stage of architectural design. This work-flow is then reversed and offers a designer a tool that generates building volumes based on target wind flow.


title: Qualitative methods in condensed matter physics

abstract: Expanding on our former hypothesis that, in the current information age, teaching physics should become more intuition-based and aiming at pattern recognition skills, we present multiple examples of qualitative methods in condensed matter physics. They include the subjects of phonons, thermal and electronic properties of matter, electron-phonon interactions and some properties of semiconductors.


title: Directly Mapping RDF Databases to Property Graph Databases

abstract: RDF triplestores and property graph databases are two approaches for data management which are based on modeling, storing, and querying graph-like data. In spite of such common principles, they present special features that complicate the task of database interoperability. While there exist some methods to transform RDF graphs into property graphs, and vice versa, they lack compatibility and a solid formal foundation. This paper presents three direct mappings (schema-dependent and schema-independent) for transforming an RDF database into a property graph database, including data and schema. We show that two of the proposed mappings satisfy the properties of semantics preservation and information preservation. The existence of both mappings allows us to conclude that the property graph data model subsumes the information capacity of the RDF data model.