title: Realizing the braided Temperley-Lieb-Jones C*-tensor categories as Hilbert C*-modules

abstract: We associate to each Temperley-Lieb-Jones C*-tensor category $\mathcal{T}\!\mathcal{L}\mathcal{J}(\delta)$ with parameter $\delta$ in the discrete range $\{2\cos(\pi/(k+2))\,:\,k=1,2,\ldots\}\cup\{2\}$ a certain C*-algebra $\mathcal{B}$ of compact operators. We use the unitary braiding on $\mathcal{T}\!\mathcal{L}\mathcal{J}(\delta)$ to equip the category $\mathrm{Mod}_{\mathcal{B}}$ of (right) Hilbert $\mathcal{B}$-modules with the structure of a braided C*-tensor category. We show that $\mathcal{T}\!\mathcal{L}\mathcal{J}(\delta)$ is equivalent, as a braided C*-tensor category, to the full subcategory $\mathrm{Mod}_{\mathcal{B}}^f$ of $\mathrm{Mod}_{\mathcal{B}}$ whose objects are those modules which admit a finite orthonormal basis. Finally, we indicate how these considerations generalize to arbitrary finitely generated rigid braided C*-tensor categories.


title: A Generic Efficient Biased Optimizer for Consensus Protocols

abstract: Consensus is one of the most fundamental distributed computing problems. In particular, it serves as a building block in many replication based fault-tolerant systems and in particular in multiple recent blockchain solutions. Depending on its exact variant and other environmental assumptions, solving consensus requires multiple communication rounds. Yet, there are known optimistic protocols that guarantee termination in a single communication round under favorable conditions. In this paper we present a generic optimizer than can turn any consensus protocol into an optimized protocol that terminates in a single communication round whenever all nodes start with the same predetermined value and no Byzantine failures occur (although node crashes are allowed). This is regardless of the network timing assumptions and additional oracle capabilities assumed by the base consensus protocol being optimized. In the case of benign failures, our optimizer works whenever the number of faulty nodes $f<n/2$. For Byzantine behavior, our optimizer's resiliency depends on the validity variant sought. In the case of classical validity, it can accommodate $f<n/4$ Byzantine failures. With the more recent external validity function assumption, it works whenever $f<n/3$. Either way, our optimizer only relies on oral messages, thereby imposing very light-weight crypto requirements.


title: Shear Destruction of Frictional Aging and Memory

abstract: We simultaneously measure the static friction and the real area of contact between two solid bodies. Under static conditions both quantities increase logarithmically in time, a phenomenon coined aging. Indeed, frictional strength is traditionally considered equivalent to the real area of contact. Here we show that this equivalence breaks down when a static shear load is applied during aging. The addition of such a shear load accelerates frictional aging while the aging rate of the real area of contact is unaffected. Moreover, a negative static shear - pulling instead of pushing - slows frictional aging, but similarly does not affect the aging of contacts. The origin of this shear effect on aging is geometrical. When shear load is increased, minute relative tilts between the two blocks prematurely erase interfacial memory prior to sliding, negating the effect of aging. Modifying the loading point of the interface eliminates these tilts and as a result frictional aging rate becomes insensitive to shear. We also identify a secondary memory-erasure effect that remains even when all tilts are eliminated and show that this effect can be leveraged to accelerate aging by cycling between two static shear loads.


title: On Bianchi type VI$_0$ spacetimes with orthogonal perfect fluid matter

abstract: We study the asymptotic behaviour of Bianchi type VI$_0$ spacetimes with orthogonal perfect fluid matter satisfying Einstein's equations. In particular, we prove a conjecture due to Wainwright about the initial singularity of such spacetimes. Using the expansion-normalized variables of Wainwright-Hsu, we demonstrate that for a generic solution the initial singularity is vacuum dominated, anisotropic and silent. In addition, by employing known results on Bianchi backgrounds, we obtain convergence results on the asymptotics of solutions to the Klein-Gordon equation on all backgrounds of this type, except for one specific case.


title: Hybrid Precoding for Multi-Group Multicasting in mmWave Systems

abstract: Multicast beamforming is known to improve spectral efficiency. However, its benefits and challenges for hybrid precoders design in millimeter-wave (mmWave) systems remain understudied. To this end, this paper investigates the first joint design of hybrid transmit precoders (with an arbitrary number of finite-resolution phase shifts) and receive combiners for mmWave multi-group multicasting. Our proposed design leverages semidefinite relaxation (SDR), alternating optimization and Cholesky matrix factorization to sequentially optimize the digital/analog precoders at the transmitter and the combiners at each receiver. By considering receivers with multiple-antenna architecture, our design remarkably improves the overall system performance. Specifically, with only two receive antennas the average transmit power per received message improves by $ 16.8\% $ while the successful information reception is boosted by $ 60\% $. We demonstrate by means of extensive simulations that our hybrid precoder design performs very close to its fully-digital counterpart even under challenging scenarios (i.e., when co-located users belong to distinct multicast groups).


title: Variational principles and thermodynamics

abstract: Variational principles play a fundamental role in deriving evolution equations of physics. They are working well in case of nondissipative evolution but for dissipative systems they are not unique, not predictive and not constructive. With methods of modern nonequilibrium thermodynamics, one can derive evolution equations for dissipative phenomena and, surprisingly, can also reproduce the Euler-Lagrange form of the evolution equations for ideal processes. In this work, we examine some demonstrative examples and compare thermodynamic and variational techniques. Then, we argue that instead of searching for variational principles for dissipative systems, a different program can be more fruitful: the second law alone can be an effective tool to construct both dissipative and nondissipative evolution equations.


title: Interfacial contributions to spin-orbit torque and magnetoresistance in ferromagnet/heavy-metal bilayers

abstract: The thickness dependence of spin-orbit torque and magnetoresistance in ferromagnet/heavy-metal bilayers is studied using the first-principles non-equilibrium Green's function formalism combined with the Anderson disorder model. A systematic expansion in orthogonal vector spherical harmonics is used for the angular dependence of the torque. The damping-like torque in Co/Pt and Co/Au bilayers can be described as a sum of the spin-Hall contribution, which increases with thickness in agreement with the spin-diffusion model, and a comparable interfacial contribution. The magnetoconductance in the plane perpendicular to the current in Co/Pt bilayers is of the order of a conductance quantum per interfacial atom, exceeding the prediction of the spin-Hall model by more than an order of magnitude. This suggests that the "spin-Hall magnetoresistance," similarly to the damping-like torque, has a large interfacial contribution unrelated to the spin-Hall effect.


title: Rendering Point Clouds with Compute Shaders

abstract: We propose a compute shader based point cloud rasterizer with up to 10 times higher performance than classic point-based rendering with the GL_POINT primitive. In addition to that, our rasterizer offers 5 byte depth-buffer precision with uniform or customizable distribution, and we show that it is possible to implement a high-quality splatting method that blends together overlapping fragments while still maintaining higher frame-rates than the traditional approach.


title: Opposite skew left braces and applications

abstract: Given a skew left brace $\mathfrak{B}$, we introduce the notion of an "opposite" skew left brace $\mathfrak{B}'$, which is closely related to the concept of the opposite of a group, and provide several applications. Skew left braces are closely linked with both solutions to the Yang-Baxter Equation and Hopf-Galois structures on Galois field extensions. We show that the set-theoretic solution to the YBE given by $\mathfrak{B}'$ is the inverse to the solution given by $\mathfrak{B}$; this allows us to identify the group-like elements in the Hopf algebra providing the Hopf-Galois structure using only these solutions. We also show how left ideals of $\mathfrak{B}'$ correspond to the realizable intermediate fields of a certain Hopf-Galois extension of a Galois extension.


title: The Nine Axes of Merit for Technosignature Searches

abstract: The diverse methodologies and myriad orthogonal proposals for the best technosignatures to search for in SETI can make it difficult to develop an effective and balanced search strategy, especially from a funding perspective. Here I propose a framework to compare the relative advantages and disadvantages of various proposed technosignatures based on nine "axes of merit". This framework was first developed at the NASA Technosignatures Workshop in Houston in 2018 and published in that report. I give the definition and rationale behind the nine axes as well as the history of each axis in the SETI and technosignature literature. These axes are then applied to three example classes of technosignature searches as an illustration of their use. An open-source software tool is available to allow technosignature researchers to make their own version of the figure.


title: Bayesian Structure Learning in Graphical Models using Shrinkage priors

abstract: We consider the problem of learning the structure of a high dimensional precision matrix under sparsity assumptions. We propose to use a shrinkage prior, called the DL-graphical prior based on the Dirichlet-Laplace prior used for the Gaussian mean problem. A posterior sampling scheme based on Gibbs sampling is also provided along with theoretical guarantees of the method by obtaining the posterior convergence rate of the precision matrix.


title: Energy partition in two M-class circular-ribbon flares

abstract: In this paper, we investigate the energy partition of two homologous M1.1 circular-ribbon flares (CRFs) in active region (AR) 12434. They were observed by \textit{SDO}, \textit{GOES}, and \textit{RHESSI} on 2015 October 15 and 16, respectively. The peak thermal energy, nonthermal energy of flare-accelerated electrons, total radiative loss of hot plasma, and radiant energies in 1$-$8 {\AA} and 1$-$70 {\AA} of the flares are calculated. The two flares have similar energetics. The peak thermal energies are (1.94$\pm$0.13)$\times$10$^{30}$ erg. The nonthermal energies in flare-accelerated electrons are (3.9$\pm$0.7)$\times$10$^{30}$ erg. The radiative outputs of the flare loops in 1$-$70 {\AA}, which are $\sim$200 times greater than the outputs in 1$-$8 {\AA}, account for $\sim$62.5\% of the peak thermal energies. The radiative losses of SXR-emitting plasma are one order of magnitude lower than the peak thermal energies. Therefore, the total heating requirements of flare loops including radiative loss are (2.1$\pm$0.1)$\times$10$^{30}$ erg, which could sufficiently be supplied by nonthermal electrons.


title: Interpretable and Fine-Grained Visual Explanations for Convolutional Neural Networks

abstract: To verify and validate networks, it is essential to gain insight into their decisions, limitations as well as possible shortcomings of training data. In this work, we propose a post-hoc, optimization based visual explanation method, which highlights the evidence in the input image for a specific prediction. Our approach is based on a novel technique to defend against adversarial evidence (i.e. faulty evidence due to artefacts) by filtering gradients during optimization. The defense does not depend on human-tuned parameters. It enables explanations which are both fine-grained and preserve the characteristics of images, such as edges and colors. The explanations are interpretable, suited for visualizing detailed evidence and can be tested as they are valid model inputs. We qualitatively and quantitatively evaluate our approach on a multitude of models and datasets.


title: Towards 100,000-pixel microcalorimeter arrays using multi-absorber transition-edge sensors

abstract: We report on the development of multi-absorber transition edge sensors (TESs), referred to as hydras. A hydra consists of multiple x-ray absorbers each with a different thermal conductance to a TES. Position information is encoded in the pulse shape. With some trade-off in performance, hydras enable very large format arrays without the prohibitive increase in bias and read-out components associated with arrays of individual TESs. Hydras are under development for the next generation of space telescope such as Lynx. Lynx is a NASA concept under study that will combine a < 1 arcsecond angular resolution optic with 100,000-pixel microcalorimeter array with energy resolution of deltaE_FWHM ~ 3 eV in the soft x-ray energy range. We present first results from hydras with 25-pixels for Lynx. Designs with absorbers on a 25 micron and 50 micron pitch are studied. Arrays incorporate, for the first time, microstrip buried wiring layers of suitable pitch and density required to readout a full-scale Lynx array. The resolution from the coadded energy histogram including all 25-pixels was deltaE_FWHM = 1.66+/-0.02 eV and 3.34+/-0.06 eV at an energy of 1.5 keV for the 25 micron and 50 micron absorber designs respectively. Position discrimination is demonstrated from parameterization of the rise-time.


title: Construction of High Codimension Ancient Mean Curvature Flows

abstract: We construct a class of compact ancient solutions to the mean curvature flow in Euclidean space with high codimension. In particular, we construct higher codimensional ancient curve shortening flows. Moreover, we characterize the asymptotic behavior of these solutions. Add on remark: the construction in this paper has been discovered by Altschuler D.-Altschuler S.-Angenent-Wu in [AAAW13].


title: Cable-driven robotic interface for lower limb neuromechanics identification

abstract: This paper presents a versatile cable-driven robotic interface to investigate the single-joint joint neuromechanics of the hip, knee and ankle. This endpoint-based interface offers highly dynamic interaction and accurate position control, as is typically required for neuromechanics identification. It can be used with the subject upright, corresponding to natural posture during walking or standing, and does not impose kinematic constraints on a joint, in contrast to existing interfaces. Mechanical evaluations demonstrated that the interface yields a rigidity above 500N/m with low viscosity. Tests with a rigid dummy leg and linear springs show that it can identify the mechanical impedance of a limb accurately. A smooth perturbation is developed and tested with a human subject, which can be used to estimate the hip neuromechanics.


title: A 21 cm pilot survey for pulsars and transients using the Focal L-Band Array for the Green Bank Telescope

abstract: Phased Array Feed (PAF) receivers are at the forefront of modern day radio astronomy. PAFs are currently being developed for spectral line and radio continuum surveys and to search for pulsars and fast radio bursts. Here, we present results of the pilot survey for pulsars and fast radio bursts using the Focal plane L-band Array for the Green Bank Telescope (FLAG) receiver operating in the frequency range of 1.3--1.5 GHz. With a system temperature of $\sim$18 K, the receiver provided unprecedented sensitivity to the survey over an instantaneous field of view (FoV) of 0.1 deg$^{2}$. For the survey, we implemented both time and frequency domain search pipelines designed to find pulsars and fast radio bursts that were validated by test pulsar observations. Although no new sources were found, we were able to demonstrate the capability of this instrument from observations of known pulsars. We report an upper limit on the rate of fast radio bursts above a fluence of 0.36~Jy ms to be 1.3 $\times$ 10$^6$ events per day per sky. Using population simulations, we show that the FLAG will find a factor of 2--3 more pulsars in same survey duration compared to its single pixel counterpart at the Green Bank Telescope. We also demonstrate that the new phased array receiver, ALPACA for the Arecibo telescope, will be a superior survey instrument and will find pulsars at a higher rate than most contemporary receivers by a factor of 2--10.


title: Peering into the Anneal Process of a Quantum Annealer

abstract: Commercial adiabatic quantum annealers have the potential to solve important NP-hard optimization problems efficiently. The newest generation of those machines additionally allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the annealing. In this work we use the aforementioned feature of the D-Wave 2000Q to attempt to monitor how the anneal solution evolves during the anneal process. This process we call slicing: at each time slice during the anneal, we are able to obtain an approximate distribution of anneal solutions. We use our technique to obtain a variety of insights into the D-Wave 2000Q. For example, we observe when individual bits flip during the anneal process and when they stabilize, which allows us to determine the freeze-out point for each qubit individually. We highlight our results using both random QUBO (quadratic unconstrained binary optimization) instances and, for better visualization, instances which we specifically optimize (using our own genetic algorithm) to exhibit a pronounced evolution of its solution during the anneal.


title: Approximating the Convex Hull via Metric Space Magnitude

abstract: Magnitude of a finite metric space and the related notion of magnitude functions on metric spaces is an active area of research in algebraic topology. Magnitude originally arose in the context of biology, where it represents the number of effective species in an environment; when applied to a one-parameter family of metric spaces $tX$ with scale parameter $t$, the magnitude captures much of the underlying geometry of the space. Prior work has mostly focussed on properties of magnitude in a global sense; in this paper we restrict the sets to finite subsets of Euclidean space and investigate its individual components. We give an explicit formula for the corrected inclusion-exclusion principle, and define a quantity associated with each point, called the $\textit{moment}$ which gives an intrinsic ordering to the points. We exploit this in order to form an algorithm which approximates the convex hull.


title: New structurally unstable families of planar vector fields

abstract: We study global bifurcations in generic 3-parameter families of vector fields on $S^2$. In the recent article [arXiv:1506.06797], Ilyashenko, Kudryashov, and Schurov show that 3-parameter unfoldings of vector fields with the polycycle "tears of the heart" are structurally unstable. We consider 3-parameter unfoldings of vector fields with separatrix graphs "ears" and "glasses", and prove that these families are structurally unstable as well. We also study in more details the classical bifurcation of a saddle loop, and use it as a building block in our main example.


title: A distributed active subspace method for scalable surrogate modeling of function valued outputs

abstract: We present a distributed active subspace method for training surrogate models of complex physical processes with high-dimensional inputs and function valued outputs. Specifically, we represent the model output with a truncated Karhunen-Lo\`eve (KL) expansion, screen the structure of the input space with respect to each KL mode via the active subspace method, and finally form an overall surrogate model of the output by combining surrogates of individual output KL modes. To ensure scalable computation of the gradients of the output KL modes, needed in active subspace discovery, we rely on adjoint-based gradient computation. The proposed method combines benefits of active subspace methods for input dimension reduction and KL expansions used for spectral representation of the output field. We provide a mathematical framework for the proposed method and conduct an error analysis of the mixed KL active subspace approach. Specifically, we provide an error estimate that quantifies errors due to active subspace projection and truncated KL expansion of the output. We demonstrate the numerical performance of the surrogate modeling approach with an application example from biotransport.


title: Multi-Quark Colour Decompositions from Unitarity

abstract: Any loop QCD amplitude at full colour is constructed from kinematic and gauge-group building blocks. In a unitarity-based on-shell framework, both objects can be reconstructed from their respective counterparts in tree-level amplitudes. This procedure is at its most powerful when aligned with flexible colour decompositions of tree-level QCD amplitudes. In this note we derive such decompositions for amplitudes with an arbitrary number of quarks and gluons from the same principle that is used to bootstrap kinematics - unitarity factorisation. In the process we formulate new multi-quark bases and provide closed-form expressions for the new decompositions. We then elaborate upon their application in colour decompositions of loop multi-quark amplitudes.


title: Finsler metrics on surfaces admitting three projective vector fields

abstract: We show that in dimension 2 every Finsler metric with at least 3-dimensional Lie algebra of projective vector fields is locally projectively equivalent to a Randers metric. We give a short list of such Finsler metrics which is complete up to coordinate change and projective equivalence.


title: Physical geometry of the quasispherical Szekeres models

abstract: The quasispherical Szekeres metric is an exact solution to Einstein's equations describing an inhomogeneous and anisotropic cosmology. Though its governing equations are well-known, there are subtle, often-overlooked details in how the model's functions relate to its physical layout, including the shapes and relative positions of structures. We present an illustrated overview of the quasispherical Szekeres models and show exactly how the model functions relate to the physical shape and distribution of matter. In particular, we describe a shell rotation effect that has not previously been fully understood. We show how this effect relates to other known properties, and lay out some mathematical tools useful for constructing models and picturing them accurately.


title: Probing the Origin of Stellar Flares on M dwarfs Using TESS Data Sectors 1-3

abstract: Detailed studies of the Sun have shown that sunspots and solar flares are closely correlated. Photometric data from Kepler/K2 has allowed similar studies to be carried out on other stars. Here, we utilise TESS photometric 2-min cadence of 167 low mass stars from Sectors 1 - 3 to investigate the relationship between starspots and stellar flares. From our sample, 90 percent show clear rotational modulation likely due to the presence of a large, dominant starspot and we use this to determine a rotational period for each star. Additionally, each low mass star shows one or more flares in its lightcurve and using Gaia DR2 parallaxes and SkyMapper magnitudes we can estimate the energy of the flares in the TESS band-pass. Overall, we have 1834 flares from the 167 low mass stars with energies from $6.0\times 10^{29}$ - $2.4\times 10^{35}$~erg. We find none of the stars in our sample show any preference for rotational phase suggesting the lack of a correlation between the large, dominant star spot and flare number. We discuss this finding in greater detail and present further scenarios to account for the origin of flares on these low mass stars.


title: Search for dark photons in decays of Higgs bosons produced in association with Z bosons in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search is presented for a Higgs boson that is produced in association with a Z boson and that decays to an undetected particle together with an isolated photon. The search is performed by the CMS Collaboration at the Large Hadron Collider using a data set corresponding to an integrated luminosity of 137 fb$^{-1}$ recorded at a center-of-mass energy of 13 TeV. No significant excess of events above the expectation from the standard model background is found. The results are interpreted in the context of a theoretical model in which the undetected particle is a massless dark photon. An upper limit is set on the product of the cross section for associated Higgs and Z boson production and the branching fraction for such a Higgs boson decay, as a function of the Higgs boson mass. For a mass of 125 GeV, assuming the standard model production cross section, this corresponds to an observed (expected) upper limit on this branching fraction of 4.6 (3.6)% at 95% confidence level. These are the first limits on Higgs boson decays to final states that include an undetected massless dark photon.


title: Multiscatter capture of superheavy dark matter by Pop. III stars

abstract: If captured by the gravitational field of stars or other compact objects, dark matter can self-annihilate and produce a potentially detectable particle flux. In the case of superheavy dark matter ($ m_{X} \gtrsim 10^{8} GeV $), a large number of scattering events with nuclei inside stars are necessary to slow down the dark matter particles below the escape velocity of the stars, at which point the Dark Matter (DM) particle becomes trapped, or captured. Using the recently developed analytical formalism for multiscatter capture, combined with the latest results on the constraints of dark-matter-baryon scattering cross-section, we calculate upper bounds on the capture rates for superheavy dark matter particles by the first (Pop. III) stars. Assuming that a non-zero fraction of the products of captured superheavy dark matter (SHDM) annihilations can be trapped and thermalized inside the star we find that this additional heat source could influence the evolutionary phase of Pop. III stars. Moreover, requiring that Pop. III stars shine with sub-Eddington luminosity, we find upper bounds on the masses of the Pop. III stars. This implies a DM dependent cutoff on the initial mass function (IMF) of Pop. III stars, thus opening up the intriguing possibility of constraining DM properties using the IMF of extremely metal-poor stars.


title: Projectively equivalent Finsler metrics on surfaces of negative Euler characteristic

abstract: We proof that on a surface of negative Euler characteristic, two real-analytic Finsler metrics have the same unparametrized oriented geodesics, if and only if they differ by a scaling constant and addition of a closed 1-form.


title: Performance Comparison for Neuroscience Application Benchmarks

abstract: Researchers within the Human Brain Project and related projects have in the last couple of years expanded their needs for high-performance computing infrastructures. The needs arise from a diverse set of science challenges that range from large-scale simulations of brain models to processing of extreme-scale experimental data sets. The ICEI project, which is in the process of creating a distributed infrastructure optimised for brain research, started to build-up a set of benchmarks that reflect the diversity of applications in this field. In this paper we analyse the performance of some selected benchmarks on an IBM POWER8 and Intel Skylake based systems with and without GPUs.


title: Singlet Pathway to the Ground State of Ultracold Polar Molecules

abstract: Starting from weakly bound Feshbach molecules, we demonstrate a two-photon pathway to the dipolar ground state of bi-alkali molecules that involves only singlet-to-singlet optical transitions. This pathway eliminates the search for a suitable intermediate state with sufficient singlet-triplet mixing and the exploration of its hyperfine structure, as is typical for pathways starting from triplet dominated Feshbach molecules. By selecting a Feshbach state with a stretched singlet hyperfine component and controlling the polarization of the excitation laser, we assure coupling to only a single hyperfine component of the $\textrm{A}^{1}\Sigma^{+}$ excited potential, even if the hyperfine structure is not resolved. Similarly, we address a stretched hyperfine component of the $\textrm{X}^{1}\Sigma^{+}$ rovibrational ground state, and therefore an ideal three level system is established. We demonstrate this pathway with ${}^{6}\textrm{Li}{}^{40}\textrm{K}$ molecules. By exploring deeply bound states of the $\textrm{A}^{1}\Sigma^{+}$ potential, we are able to obtain large and balanced Rabi frequencies for both transitions. This method can be applied to other molecular species.


title: Unified Simulation and Test Platform for Control Systems of Unmanned Vehicles

abstract: Control systems on unmanned vehicles are safety-critical systems whose requirements on reliability and safety are ever-increasing. Currently, testing a complex autonomous control system is an expensive and time-consuming process, which requires massive repeated experimental testing during the whole development stage. This paper presents a unified simulation and test platform for vehicle autonomous control systems aiming to significantly improve the development speed and safety level of unmanned vehicles. First, a unified modular modeling framework compatible with different types of vehicles is proposed with methods to ensure modeling credibility. Then, the simulation software system is developed by the model-based design framework, whose modular programming methods and automatic code generation functions ensure the efficiency, credibility, and standardization of the system development process. Finally, an FPGA-based real-time hardware-in-the-loop simulation platform is proposed to ensure the comprehensiveness and credibility of the simulation and test results. In the end, the proposed platform is applied to a multicopter control system. By comparing with experimental results, the accuracy and credibility of the simulation testing results are verified by using the simulation credibility assessment method proposed in our previous work. To verify the practicability of the proposed platform, several successful applications are presented for the multicopter rapid prototyping, estimation algorithm verification, autonomous flight testing, and automatic safety testing with automatic fault injection and result evaluation of unmanned vehicles.


title: A non-nuclear $C^*$-algebra with the Weak Expectation Property and the Local Lifting Property

abstract: We construct the first example of a $C^*$-algebra $A$ with the properties in the title. This gives a new example of non-nuclear $A$ for which there is a unique $C^*$-norm on $A \otimes A^{op}$. This example is of particular interest in connection with the Connes-Kirchberg problem, which is equivalent to the question whether $C^*({\bb F}_2)$, which is known to have the LLP, also has the WEP. Our $C^*$-algebra $A$ has the same collection of finite dimensional operator subspaces as $C^*({\bb F}_2)$ or $C^*({\bb F}_\infty)$. In addition our example can be made to be quasidiagonal and of similarity degree (or length) 3. In the second part of the paper we reformulate our construction in the more general framework of a $C^*$-algebra that can be described as the \emph{limit both inductive and projective} for a sequence of $C^*$-algebras $(C_n)$ when each $C_n$ is a \emph{subquotient} of $C_{n+1}$. We use this to show that for certain local properties of injective (non-surjective) $*$-homomorphisms, there are $C^*$-algebras for which the identity map has the same properties as the $*$-homomorphisms.


title: Zero-Shot Deep Hashing and Neural Network Based Error Correction for Face Template Protection

abstract: In this paper, we present a novel architecture that integrates a deep hashing framework with a neural network decoder (NND) for application to face template protection. It improves upon existing face template protection techniques to provide better matching performance with one-shot and multi-shot enrollment. A key novelty of our proposed architecture is that the framework can also be used with zero-shot enrollment. This implies that our architecture does not need to be re-trained even if a new subject is to be enrolled into the system. The proposed architecture consists of two major components: a deep hashing (DH) component, which is used for robust mapping of face images to their corresponding intermediate binary codes, and a NND component, which corrects errors in the intermediate binary codes that are caused by differences in the enrollment and probe biometrics due to factors such as variation in pose, illumination, and other factors. The final binary code generated by the NND is then cryptographically hashed and stored as a secure face template in the database. The efficacy of our approach with zero-shot, one-shot, and multi-shot enrollments is shown for CMU-PIE, Extended Yale B, WVU multimodal and Multi-PIE face databases. With zero-shot enrollment, the system achieves approximately 85% genuine accept rates (GAR) at 0.01% false accept rate (FAR), and with one-shot and multi-shot enrollments, it achieves approximately 99.95% GAR at 0.01% FAR, while providing a high level of template security.


title: Highly accurate local basis sets for large-scale DFT calculations in CONQUEST

abstract: Given the widespread use of density functional theory (DFT), there is an increasing need for the ability to model large systems (beyond 1,000 atoms). We present a brief overview of the large-scale DFT code Conquest, which is capable of modelling such large systems, and discuss approaches to the generation of consistent, well-converged pseudo-atomic basis sets which will allow such large scale calculations. We present tests of these basis sets for a variety of materials, comparing to fully converged plane wave results using the same pseudopotentials and grids.


title: On double-membership graphs of models of Anti-Foundation

abstract: We answer some questions about graphs which are reducts of countable models of Anti-Foundation, obtained by considering the binary relation of double-membership $x\in y\in x$. We show that there are continuum-many such graphs, and study their connected components. We describe their complete theories and prove that each has continuum-many countable models, some of which are not reducts of models of Anti-Foundation.


title: A minimal core calculus for Solidity contracts

abstract: The Ethereum platform supports the decentralized execution of smart contracts, i.e. computer programs that transfer digital assets between users. The most common language used to develop these contracts is Solidity, a Javascript-like language which compiles into EVM bytecode, the language actually executed by Ethereum nodes. While much research has addressed the formalisation of the semantics of EVM bytecode, relatively little attention has been devoted to that of Solidity. In this paper we propose a minimal calculus for Solidity contracts, which extends an imperative core with a single primitive to transfer currency and invoke contract procedures. We build upon this formalisation to give semantics to the Ethereum blockchain. We show our calculus expressive enough to reason about some typical quirks of Solidity, like e.g. re-entrancy.


title: Maximum likelihood convolutional beamformer for simultaneous denoising and dereverberation

abstract: This article describes a probabilistic formulation of a Weighted Power minimization Distortionless response convolutional beamformer (WPD). The WPD unifies a weighted prediction error based dereverberation method (WPE) and a minimum power distortionless response beamformer (MPDR) into a single convolutional beamformer, and achieves simultaneous dereverberation and denoising in an optimal way. However, the optimization criterion is obtained simply by combining existing criteria without any clear theoretical justification. This article presents a generative model and a probabilistic formulation of a WPD, and derives an optimization algorithm based on a maximum likelihood estimation. We also describe a method for estimating the steering vector of the desired signal by utilizing WPE within the WPD framework to provide an effective and efficient beamformer for denoising and dereverberation.


title: I Bet You Are Wrong: Gambling Adversarial Networks for Structured Semantic Segmentation

abstract: Adversarial training has been recently employed for realizing structured semantic segmentation, in which the aim is to preserve higher-level scene structural consistencies in dense predictions. However, as we show, value-based discrimination between the predictions from the segmentation network and ground-truth annotations can hinder the training process from learning to improve structural qualities as well as disabling the network from properly expressing uncertainties. In this paper, we rethink adversarial training for semantic segmentation and propose to formulate the fake/real discrimination framework with a correct/incorrect training objective. More specifically, we replace the discriminator with a "gambler" network that learns to spot and distribute its budget in areas where the predictions are clearly wrong, while the segmenter network tries to leave no clear clues for the gambler where to bet. Empirical evaluation on two road-scene semantic segmentation tasks shows that not only does the proposed method re-enable expressing uncertainties, it also improves pixel-wise and structure-based metrics.


title: Saturated Free Algebras and Almost Indiscernible Theories

abstract: We extend the concept of "almost indiscernible theory" introduced by Pillay and Sklinos in [Bull. Symb. Log., 2015] (which was itself a modernization and expansion of Baldwin and Shelah [Algebra Universalis, 1983]), to uncountable languages and uncountable parameter sequences. Roughly speaking a theory $T$ is almost indiscernible if some saturated model is in the algebraic closure of an indiscernible set of sequences. We show that such a theory $T$ is nonmultidimensional, superstable, and stable in all cardinals $\ge |T |$ . We prove a structure theorem for sufficiently large $a$-models $M$: Theorem 2.10 which states that over a suitable base, $M$ is in the algebraic closure of an independent set of realizations of weight one types (in possibly infinitely many variables). We also explore further the saturated free algebras of Baldwin and Shelah in both the countable and uncountable context. We study in particular theories and varieties of $R$-modules, characterizing those rings $R$ for which the free $R$-module on $|R|^+$ generators is saturated (Theorem 3.15), and pointing out a counterexample to a conjecture from Pillay-Sklinos (Example 3.16).


title: Reference frames which separately store non-commuting conserved quantities

abstract: Even in the presence of conservation laws, one can perform arbitrary transformations on a system if given access to a suitable reference frame, since conserved quantities may be exchanged between the system and the frame. Here we explore whether these quantities can be separated into different parts of the reference frame, with each part acting as a `battery' for a distinct quantity. For systems composed of spin-$\frac12$ particles, we show that the components of angular momentum $S_x$, $S_y$ and $S_z$ (non-commuting conserved quantities) may be separated in this way, and also provide several extensions of this result. These results also play a key role in the quantum thermodynamics of non-commuting conserved quantities.


title: Relighting Humans: Occlusion-Aware Inverse Rendering for Full-Body Human Images

abstract: Relighting of human images has various applications in image synthesis. For relighting, we must infer albedo, shape, and illumination from a human portrait. Previous techniques rely on human faces for this inference, based on spherical harmonics (SH) lighting. However, because they often ignore light occlusion, inferred shapes are biased and relit images are unnaturally bright particularly at hollowed regions such as armpits, crotches, or garment wrinkles. This paper introduces the first attempt to infer light occlusion in the SH formulation directly. Based on supervised learning using convolutional neural networks (CNNs), we infer not only an albedo map, illumination but also a light transport map that encodes occlusion as nine SH coefficients per pixel. The main difficulty in this inference is the lack of training datasets compared to unlimited variations of human portraits. Surprisingly, geometric information including occlusion can be inferred plausibly even with a small dataset of synthesized human figures, by carefully preparing the dataset so that the CNNs can exploit the data coherency. Our method accomplishes more realistic relighting than the occlusion-ignored formulation.


title: Thermodynamics of spin-orbit coupled bosons in two dimensions from complex Langevin

abstract: We investigate the thermal properties of interacting spin-orbit coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state non-perturbatively in a range of spin-orbit couplings and chemical potentials. Our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudo-condensates. These have been observed to be destroyed by the spin-orbit interactions.


title: Compressive-strain induced enhancement of exchange interactions and short-range magnetic order in Sr$_2$IrO$_4$ investigated by Raman spectroscopy

abstract: We have carried out Raman spectroscopy experiments to investigate two-magnon excitations in epitaxial thin films of the quasi-two-dimensional antiferromagnetic Mott insulator Sr$_2$IrO$_4$ under in-plane misfit strain. With in-plane biaxial compression, the energy of the two-magnon peak increases, and the peak remains observable over a wider temperature range above the N\'eel temperature, indicating strain-induced enhancement of the superexchange interactions between $\it{J}_{eff}$ = 1/2 pseudospins. From density functional theory calculations, we have found an increase of the nearest-neighbor hopping parameter and exchange interaction with increasing biaxial compressive strain, in agreement with the experimental observations. Our experimental and theoretical results provide perspectives for systematic, theory-guided strain control of the primary exchange interactions in 5$\it{d}$ transition metal oxides.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.supr-con


title: MUSE library of stellar spectra

abstract: Empirical stellar spectral libraries have applications in both extragalactic and stellar studies, and they have an advantage over theoretical libraries because they naturally include all relevant chemical species and physical processes. During recent years we see a stream of new high quality sets of spectra, but increasing the spectral resolution and widening the wavelength coverage means resorting to multi-order echelle spectrographs. Assembling the spectra from many pieces results in lower fidelity of their shapes. We aim to offer the community a library of high signal-to-noise spectra with reliable continuum shapes. Furthermore, the using an integral field unit (IFU) alleviates the issue of slit losses. Our library was build with the MUSE (Multi-Unit Spectroscopic Explorer) IFU instrument. We obtained spectra over nearly the entire visual band (lambda~4800-9300 Ang). We assembled a library of 35 high-quality MUSE spectra for a subset of the stars from the X-shooter Spectral Library. We verified the continuum shape of these spectra with synthetic broad band colors derived from the spectra. We also report some spectral indices from the Lick system, derived from the new observations. We offer a high-fidelity set of stellar spectra that covers the Hertzsprung-Russell diagram. It can be useful for both extragalactic and stellar work and demonstrates that the IFUs are excellent tools for building reliable spectral libraries.


title: A Characterization of Mean Squared Error for Estimator with Bagging

abstract: Bagging can significantly improve the generalization performance of unstable machine learning algorithms such as trees or neural networks. Though bagging is now widely used in practice and many empirical studies have explored its behavior, we still know little about the theoretical properties of bagged predictions. In this paper, we theoretically investigate how the bagging method can reduce the Mean Squared Error (MSE) when applied on a statistical estimator. First, we prove that for any estimator, increasing the number of bagged estimators $N$ in the average can only reduce the MSE. This intuitive result, observed empirically and discussed in the literature, has not yet been rigorously proved. Second, we focus on the standard estimator of variance called unbiased sample variance and we develop an exact analytical expression of the MSE for this estimator with bagging. This allows us to rigorously discuss the number of iterations $N$ and the batch size $m$ of the bagging method. From this expression, we state that only if the kurtosis of the distribution is greater than $\frac{3}{2}$, the MSE of the variance estimator can be reduced with bagging. This result is important because it demonstrates that for distribution with low kurtosis, bagging can only deteriorate the performance of a statistical prediction. Finally, we propose a novel general-purpose algorithm to estimate with high precision the variance of a sample.


title: Properties and nature of Be stars 31. The binary nature, light variability, physical elements, and emission-line changes of HD~81357

abstract: Reliable determination of the basic physical properties of hot emission-line binaries with Roche-lobe filling secondaries is important for developing the theory of mass exchange in binaries. It is not easy, however, due to the presence of circumstellar matter. Here, we report the first detailed investigation of a new representative of this class of binaries, HD~81357, based on the analysis of spectra and photometry from several observatories. HD~81357 was found to be a double-lined spectroscopic binary and an ellipsoidal variable seen under an intermediate orbital inclination of $\sim(63\pm5)^\circ$, having an orbital period of 33\fd77445(41) and a~circular orbit. From an automated comparison of the observed and synthetic spectra, we estimate the component's effective temperatures to be 12930(540)~K and 4260(24)~K. The combined light-curve and orbital solutions, also constrained by a very accurate Gaia Data Release 2 parallax, give the following values of the basic physical properties: masses $3.36\pm0.15$ and $0.34\pm0.04$~\Mnom, radii $3.9\pm0.2$ and 13.97\pm0.05$~\Rnom, and a~mass ratio $10.0\pm0.5$. Evolutionary modelling of the system including the phase of mass transfer between the components indicated that HD~81357 is a~system observed in the final slow phase of the mass exchange after the mass-ratio reversal. Contrary to what has been seen for similar binaries like AU~Mon, no cyclic light variations were found on a~time scale an~order of magnitude longer than the orbital period. 243,1 15%


title: WFIRST and EUCLID: enabling the microlensing parallax measurement from space

abstract: The Wide Field Infrared Survey Telescope (WFIRST) is expected to detect hundreds of free-floating planets, but it will not be able to measure their masses. However, simultaneous microlensing observations by both Euclid and WFIRST spacecraft, separated by ~ 100, 000 km in orbits around the Sun-Earth L2 Lagrange point, will enable measurements of microlensing parallax for low-mass lenses such as free-floating planets. Using simple Fisher matrix estimates of the parallax measurement uncertainties, we show that high-cadence observations by Euclid could be used to measure ~ 1 free-floating planet microlens parallax per 6 days of simultaneous Euclid observations. Accounting for Euclid's pointing constraints, it could therefore potentially measure ~ 20 free-floating planet parallaxes with 120 days of observations split equally between Euclid's main mission and an extended mission, with a potential to increase this number if spacecraft pointing constraints can be relaxed after the end of the main mission. These Euclid observations would also provide additional mass measurements or cross-checks for larger numbers of WFIRST's bound planets, among other benefits to several science cases.


title: Faster Tensor Train Decomposition for Sparse Data

abstract: In recent years, the application of tensors has become more widespread in fields that involve data analytics and numerical computation. Due to the explosive growth of data, low-rank tensor decompositions have become a powerful tool to harness the notorious curse of dimensionality. The main forms of tensor decomposition include CP decomposition, Tucker decomposition, tensor train (TT) decomposition, etc. Each of the existing TT decomposition algorithms, including the TT-SVD and randomized TT-SVD, is successful in the field, but neither can both accurately and efficiently decompose large-scale sparse tensors. Based on previous research, this paper proposes a new quasi-best fast TT decomposition algorithm for large-scale sparse tensors with proven correctness and the upper bound of its complexity is derived. In numerical experiments, we verify that the proposed algorithm can decompose sparse tensors faster than the TT-SVD, and have more speed, precision and versatility than randomized TT-SVD, and it can be used to decomposes arbitrary high-dimensional tensor without losing efficiency when the number of non-zero elements is limited. The new algorithm implements a large-scale sparse matrix TT decomposition that was previously unachievable, enabling tensor decomposition based algorithms to be applied in larger-scale scenarios.


title: Integrable Geometric Flows for Curves in the Pseudoconformal 3-Sphere

abstract: We consider evolution equations for curves in the 3-dimensional sphere $S^3$ that are invariant under the group $SU(2,1)$ of pseudoconformal transformations, which preserves the standard contact structure on the sphere. In particular, we investigate how invariant evolutions of Legendrian and transverse curves induce well-known integrable systems and hierarchies at the level of their geometric invariants.


title: Advocacy Learning: Learning through Competition and Class-Conditional Representations

abstract: We introduce advocacy learning, a novel supervised training scheme for attention-based classification problems. Advocacy learning relies on a framework consisting of two connected networks: 1) $N$ Advocates (one for each class), each of which outputs an argument in the form of an attention map over the input, and 2) a Judge, which predicts the class label based on these arguments. Each Advocate produces a class-conditional representation with the goal of convincing the Judge that the input example belongs to their class, even when the input belongs to a different class. Applied to several different classification tasks, we show that advocacy learning can lead to small improvements in classification accuracy over an identical supervised baseline. Though a series of follow-up experiments, we analyze when and how such class-conditional representations improve discriminative performance. Though somewhat counter-intuitive, a framework in which subnetworks are trained to competitively provide evidence in support of their class shows promise, in many cases performing on par with standard learning approaches. This provides a foundation for further exploration into competition and class-conditional representations in supervised learning.


title: Envisioning Tool Support for Designing Privacy-Aware Internet of Thing Applications

abstract: The design and development process for Internet of Things (IoT) applications is more complicated than for desktop, mobile, or web applications. IoT applications require both software and hardware to work together across multiple different types of nodes (e.g., microcontrollers, system-on-chips, mobile phones, miniaturised single-board computers, and cloud platforms) with different capabilities under different conditions. IoT applications typically collect and analyse personal data that can be used to derive sensitive information about individuals. Without proper privacy protections in place, IoT applications could lead to serious privacy violations. Thus far, privacy concerns have not been explicitly considered in software engineering processes when designing and developing IoT applications, partly due to a lack of tools, technologies, and guidance. This paper presents a research vision that argues the importance of developing a privacy-aware IoT application design tool to address the challenges mentioned above. This tool should not only transform IoT application designs into privacy-aware application designs but also validate and verify them. First, we outline how this proposed tool should work in practice and its core functionalities. Then, we identify research challenges and potential directions towards developing the proposed tool. We anticipate that this proposed tool will save many engineering hours which engineers would otherwise need to spend on developing privacy expertise and applying it. We also highlight the usefulness of this tool towards privacy education and privacy compliance.


title: Towards closing the gap between the theory and practice of SVRG

abstract: Among the very first variance reduced stochastic methods for solving the empirical risk minimization problem was the SVRG method (Johnson & Zhang 2013). SVRG is an inner-outer loop based method, where in the outer loop a reference full gradient is evaluated, after which $m \in \mathbb{N}$ steps of an inner loop are executed where the reference gradient is used to build a variance reduced estimate of the current gradient. The simplicity of the SVRG method and its analysis have led to multiple extensions and variants for even non-convex optimization. We provide a more general analysis of SVRG than had been previously done by using arbitrary sampling, which allows us to analyse virtually all forms of mini-batching through a single theorem. Furthermore, our analysis is focused on more practical variants of SVRG including a new variant of the loopless SVRG (Hofman et al 2015, Kovalev et al 2019, Kulunchakov and Mairal 2019) and a variant of k-SVRG (Raj and Stich 2018) where $m=n$ and where $n$ is the number of data points. Since our setup and analysis reflect what is done in practice, we are able to set the parameters such as the mini-batch size and step size using our theory in such a way that produces a more efficient algorithm in practice, as we show in extensive numerical experiments.


title: Cascaded Revision Network for Novel Object Captioning

abstract: Image captioning, a challenging task where the machine automatically describes an image by sentences, has drawn significant attention in recent years. Despite the remarkable improvements of recent approaches, however, these methods are built upon a large set of training image-sentence pairs. The expensive labor efforts hence limit the captioning model to describe the wider world. In this paper, we present a novel network structure, Cascaded Revision Network, which aims at relieving the problem by equipping the model with out-of-domain knowledge. CRN first tries its best to describe an image using the existing vocabulary from in-domain knowledge. Due to the lack of out-of-domain knowledge, the caption may be inaccurate or include ambiguous words for the image with unknown (novel) objects. We propose to re-edit the primary captioning sentence by a series of cascaded operations. We introduce a perplexity predictor to find out which words are most likely to be inaccurate given the input image. Thereafter, we utilize external knowledge from a pre-trained object detection model and select more accurate words from detection results by the visual matching module. In the last step, we design a semantic matching module to ensure that the novel object is fit in the right position. By this novel cascaded captioning-revising mechanism, CRN can accurately describe images with unseen objects. We validate the proposed method with state-of-the-art performance on the held-out MSCOCO dataset as well as scale to ImageNet, demonstrating the effectiveness of this method.


title: Robust magnetotransport in disordered ferromagnetic kagome layers with quantum anomalous Hall effect

abstract: The magnetotransport properties of disordered ferromagnetic kagome layers are investigated numerically. We show that a large domain-wall magnetoresistance or negative magnetoresistance can be realized in kagome layered materials (e.g. Fe$_3$Sn$_2$, Co$_3$Sn$_2$S$_2$, and Mn$_3$Sn), which show the quantum anomalous Hall effect. The kagome layers show a strong magnetic anisotropy and a large magnetoresistance depending on their magnetic texture. These domain-wall magnetoresistances are expected to be robust against disorder and observed irrespective of the domain-wall thickness, in contrast to conventional domain-wall magnetoresistance in ferromagnetic metals.


title: Addressing Limited Weight Resolution in a Fully Optical Neuromorphic Reservoir Computing Readout

abstract: Using optical hardware for neuromorphic computing has become more and more popular recently due to its efficient high-speed data processing capabilities and low power consumption. However, there are still some remaining obstacles to realizing the vision of a completely optical neuromorphic computer. One of them is that, depending on the technology used, optical weighting elements may not share the same resolution as in the electrical domain. Moreover, noise and drift are important considerations as well. In this article, we investigate a new method for improving the performance of optical weighting, even in the presence of noise and in the case of very low resolution. Even with only 8 to 32 levels of resolution, the method can outperform the naive traditional low-resolution weighting by several orders of magnitude in terms of bit error rate and can deliver performance very close to full-resolution weighting elements, also in noisy environments.


title: Robust Learning with Jacobian Regularization

abstract: Design of reliable systems must guarantee stability against input perturbations. In machine learning, such guarantee entails preventing overfitting and ensuring robustness of models against corruption of input data. In order to maximize stability, we analyze and develop a computationally efficient implementation of Jacobian regularization that increases classification margins of neural networks. The stabilizing effect of the Jacobian regularizer leads to significant improvements in robustness, as measured against both random and adversarial input perturbations, without severely degrading generalization properties on clean data.


title: Cross Section for $\gamma n \to \pi^0 n$ measured at Mainz/A2

abstract: The $\gamma n \to \pi^0 n$ differential cross section evaluated for 27 energy bins span the photon-energy range 290-813 MeV (W = 1.195-1.553 GeV) and the pion c.m. polar production angles, ranging from 18 deg to 162 deg, making use of model-dependent nuclear corrections to extract pi0 production data on the neutron from measurements on the deuteron target. Additionally, the total photoabsorption cross section was measured. The tagged photon beam produced by the 883-MeV electron beam of the Mainz Microtron MAMI was used for the 0-meson production. Our accumulation of 3.6 x 10^6 $\gamma n \to \pi^0 n$ events allowed a detailed study of the reaction dynamics. Our data are in reasonable agreement with previous A2 measurements and extend them to lower energies. The data are compared to predictions of previous SAID, MAID, and BnGa partial-wave analyses and to the latest SAID fit MA19 that included our data. Selected photon decay amplitudes $N^* \to \gamma n$ at the resonance poles are determined for the first time.


title: The Relation between Composability and Splittability of Permutation Classes

abstract: A permutation class $C$ is said to be splittable if there exist two proper subclasses $A, B \subsetneq C$ such that any $\sigma \in C$ can be red-blue colored so that the red (respectively, blue) subsequence of $\sigma$ is order isomorphic to an element of $A$ (respectively, $B$). The class $C$ is said to be composable if there exists some number of proper subclasses $A_1, \dots, A_k \subsetneq C$ such that any $\sigma \in C$ can be written as $\alpha_1 \circ \dots \circ \alpha_k$ for some $\alpha_i \in A_i$. We answer a question of Karpilovskij by showing that there exists a composable permutation class that is not splittable. We also give a condition under which an infinite composable class must be splittable.


title: Correlations of multiplicative functions along deterministic and independent sequences

abstract: We study correlations of multiplicative functions taken along deterministic sequences and sequences that satisfy certain linear independence assumptions. The results obtained extend recent results of Tao and Ter\"av\"ainen and results of the author. Our approach is to use tools from ergodic theory in order to effectively exploit feedback from analytic number theory. The results on deterministic sequences crucially use structural properties of measure preserving systems associated with bounded multiplicative functions that were recently obtained by the author and Host. The results on independent sequences depend on multiple ergodic theorems obtained using the theory of characteristic factors and qualitative equidistribution results on nilmanifolds.


title: Force Balance Approach for Advanced Approximations in Density Functional Theories

abstract: We propose a systematic and constructive way to determine the exchange-correlation potentials of density-functional theories including vector potentials. The approach does not rely on energy or action functionals. Instead it is based on equations of motion of current quantities (force balance equations) and is feasible both in the ground-state and the time-dependent setting. This avoids, besides differentiability and causality issues, the optimized-effective-potential procedure of orbital-dependent functionals. We provide straightforward exchange-type approximations for different density functional theories that for a homogeneous system and no external vector potential reduce to the exchange-only local-density and Slater X$\alpha$ approximations.


title: Stochastic First-order Methods for Convex and Nonconvex Functional Constrained Optimization

abstract: Functional constrained optimization is becoming more and more important in machine learning and operations research. Such problems have potential applications in risk-averse machine learning, semisupervised learning, and robust optimization among others. In this paper, we first present a novel Constraint Extrapolation (ConEx) method for solving convex functional constrained problems, which utilizes linear approximations of the constraint functions to define the extrapolation (or acceleration) step. We show that this method is a unified algorithm that achieves the best-known rate of convergence for solving different functional constrained convex composite problems, including convex or strongly convex, and smooth or nonsmooth problems with a stochastic objective and/or stochastic constraints. Many of these rates of convergence were in fact obtained for the first time in the literature. In addition, ConEx is a single-loop algorithm that does not involve any penalty subproblems. Contrary to existing primal-dual methods, it does not require the projection of Lagrangian multipliers into a (possibly unknown) bounded set. Second, for nonconvex functional constrained problems, we introduce a new proximal point method that transforms the initial nonconvex problem into a sequence of convex problems by adding quadratic terms to both the objective and constraints. Under a certain MFCQ-type assumption, we establish the convergence and rate of convergence of this method to KKT points when the convex subproblems are solved exactly or inexactly. For large-scale and stochastic problems, we present a more practical proximal point method in which the approximate solutions of the subproblems are computed by the aforementioned ConEx method. To the best of our knowledge, most of these convergence and complexity results of the proximal point method for nonconvex problems also seem to be new in the literature.


title: Metric Learning With HORDE: High-Order Regularizer for Deep Embeddings

abstract: Learning an effective similarity measure between image representations is key to the success of recent advances in visual search tasks (e.g. verification or zero-shot learning). Although the metric learning part is well addressed, this metric is usually computed over the average of the extracted deep features. This representation is then trained to be discriminative. However, these deep features tend to be scattered across the feature space. Consequently, the representations are not robust to outliers, object occlusions, background variations, etc. In this paper, we tackle this scattering problem with a distribution-aware regularization named HORDE. This regularizer enforces visually-close images to have deep features with the same distribution which are well localized in the feature space. We provide a theoretical analysis supporting this regularization effect. We also show the effectiveness of our approach by obtaining state-of-the-art results on 4 well-known datasets (Cub-200-2011, Cars-196, Stanford Online Products and Inshop Clothes Retrieval).


title: Reconfigurable Liquid Crystal Reflectarray Metasurface for THz Communications

abstract: We present computational studies on a proposed 20 by 20 elements electronically reconfigurable liquid crystal (LC) based binary phase reflectarray metasurface, operational at 108 GHz. LC was modelled after Mer's GT3-23001, and full wave simulations have shown a phase difference of 177 degrees between ON and OFF states, while reflection amplitudes were both 0.88 for ON and OFF. We present preliminary full wave simulation results on the Genetic Algorithm (GA) optimised far-fields. We also present the basic design procedures and cross-platform implementations on optimisation routines involving Matlab, CST and VBA environments.


title: Black Hole and Cosmological Particle Production in Schwarzschild de Sitter

abstract: We compute the spectra and total fluxes of quantum mechanically produced particles crossing the black hole and cosmological horizons in Schwarzschild de Sitter (SdS). Particle states are defined with respect to well-behaved, Kruskal coordinates near the horizons, and as a consequence we find that these spectra are generally non-thermal. The non-thermal Bogoliubov coefficient for a vacuum fluctuation near the black hole horizon to produce a particle that crosses the cosmological horizon is shown to equal to the convolution of two thermal coefficients, one at the cosmological temperature and one at the black hole temperature, weighted by the transmission coefficient for wave propagation in static SdS coordinates. In this sense virtual thermal propagation underlies the production process. This representation leads to the useful result that the geometric optics approximation is reliable when used together with a low frequency cut-off determined by the transmission coefficient. The large black hole limit is a quasi-equilibrium situation as both temperatures approach the common value of zero, the particle spectra become equal, and both emissions are exponentially suppressed. Small black holes radiate as thermal bodies and absorb a tiny flux of cosmological particles. The behavior of the quantum fluctuations on the horizons is seen to be consistent with the Schottky anomaly behavior of classical gravitational fluctuations.


title: Learning Conditional Deformable Templates with Convolutional Networks

abstract: We develop a learning framework for building deformable templates, which play a fundamental role in many image analysis and computational anatomy tasks. Conventional methods for template creation and image alignment to the template have undergone decades of rich technical development. In these frameworks, templates are constructed using an iterative process of template estimation and alignment, which is often computationally very expensive. Due in part to this shortcoming, most methods compute a single template for the entire population of images, or a few templates for specific sub-groups of the data. In this work, we present a probabilistic model and efficient learning strategy that yields either universal or conditional templates, jointly with a neural network that provides efficient alignment of the images to these templates. We demonstrate the usefulness of this method on a variety of domains, with a special focus on neuroimaging. This is particularly useful for clinical applications where a pre-existing template does not exist, or creating a new one with traditional methods can be prohibitively expensive. Our code and atlases are available online as part of the VoxelMorph library at http://voxelmorph.csail.mit.edu.


title: Bayesian estimation of the functional spatial lag model

abstract: The spatial lag model (SLM) has been widely studied in the literature for spatialised data modeling in various disciplines such as geography, economics, demography, regional sciences, etc. This is an extension of the classical linear model that takes into account the proximity of spatial units in modeling. The extension of the SLM model in the functional framework (the FSLM model) as well as its estimation by the truncated maximum likelihood technique have been proposed by \cite{Ahmed}. In this paper, we propose a Bayesian estimation of the FSLM model. The Bayesian MCMC technique is used as estimation methods of the parameters of the model. A simulation study is conducted in order to compare the results of the Bayesian estimation method with the truncated maximum likelihood method. As an illustration, the proposed Bayesian method is used to establish a relationship between the unemployment rate and the curves of illiteracy rate observed in the 45 departments of Senegal.


title: Semigroup-theoretic approach to diffusion in thin layers separated by semi-permeable membranes

abstract: Using techniques of the theory of semigroups of linear operators we study the question of approximating solutions to equations governing diffusion in thin layers separated by a semi-permeable membrane. We show that as thickness of the layers converges to $0$, the solutions, which by nature are functions of $3$ variables, gradually lose dependence on the vertical variable and thus may be regarded as functions of $2$ variables. The limit equation describes diffusion on the lower and upper sides of a two-dimensional surface (the membrane) with jumps from one side to the other. The latter possibility is expressed as an additional term in the generator of the limit semigroup, and this term is build from permeability coefficients of the membrane featuring in the transmission conditions of the approximating equations (i.e., in the description of the domains of the generators of the approximating semigroups). We prove this convergence result in the spaces of square integrable and continuous functions, and study the way the choice of transmission conditions influences the limit.


title: Parallel Finger Search Structures

abstract: In this paper we present two versions of a parallel finger structure FS on p processors that supports searches, insertions and deletions, and has a finger at each end. This is to our knowledge the first implementation of a parallel search structure that is work-optimal with respect to the finger bound and yet has very good parallelism (within a factor of O( (log p)^2 ) of optimal). We utilize an extended implicit batching framework that transparently facilitates the use of FS by any parallel program P that is modelled by a dynamically generated DAG D where each node is either a unit-time instruction or a call to FS. The total work done by either version of FS is bounded by the finger bound F[L] (for some linearization L of D ), i.e. each operation on an item with finger distance r takes O( log r + 1 ) amortized work; it is cheaper for items closer to a finger. Running P using the simpler version takes O( ( T[1] + F[L] ) / p + T[inf] + d * ( (log p)^2 + log n ) ) time on a greedy scheduler, where T[1],T[inf] are the size and span of D respectively, and n is the maximum number of items in FS, and d is the maximum number of calls to FS along any path in D. Using the faster version, this is reduced to O( ( T[1] + F[L] ) / p + T[inf] + d * (log p)^2 + s[L] ) time, where s[L] is the weighted span of D where each call to FS is weighted by its cost according to F[L]. We also sketch how to extend FS to support a fixed number of movable fingers. The data structures in our paper fit into the dynamic multithreading paradigm, and their performance bounds are directly composable with other data structures given in the same paradigm. Also, the results can be translated to practical implementations using work-stealing schedulers.


title: How planets grow by pebble accretion II: Analytical calculations on the evolution of polluted envelopes

abstract: Proto-planets embedded in their natal disks acquire hot envelopes as they grow and accrete solids. This ensures that the material they accrete - pebbles, as well as (small) planetesimals - will vaporize to enrich their atmospheres. Enrichment modifies an envelope's structure and significantly alters its further evolution. Our aim is to describe the formation of planets with polluted envelopes from the moment that impactors begin to sublimate to beyond the disk's eventual dissipation. We constructed an analytical interior structure model, characterized by a hot and uniformly mixed high-Z vapor layer surrounding the core, located below the usual unpolluted radiative-convective regions. The evolution of planets with uniformly mixed polluted envelopes follows four potential phases. Initially, the central core grows directly through impacts and rainout until the envelope becomes hot enough to vaporize and absorb all incoming solids. We find that a planet reaches runaway accretion when the sum of its core and vapor mass exceeds a value that we refer to as the critical metal mass - a criterion that supersedes the traditional critical core mass. It scales positively with both the pollutant's evaporation temperature and with the planet's core mass. Hence, planets at shorter orbital separations require the accretion of more solids to reach runaway as they accrete less volatile materials. If the solids accretion rate dries up, we identify the decline of the mean molecular weight - dilution - as a mechanism to limit gas accretion during a polluted planet's embedded cooling phase. When the disk ultimately dissipates, the envelope's inner temperature declines and its vapor eventually rains out, augmenting the mass of the core. The energy release that accompanies this does not result in significant mass-loss, as it only occurs after the planet has substantially contracted.


title: Byzantine Approximate Agreement on Graphs

abstract: Consider a distributed system with $n$ processors out of which $f$ can be Byzantine faulty. In the approximate agreement task, each processor $i$ receives an input value $x_i$ and has to decide on an output value $y_i$ such that - the output values are in the convex hull of the non-faulty processors' input values, - the output values are within distance $d$ of each other. Classically, the values are assumed to be from an $m$-dimensional Euclidean space, where $m \ge 1$. In this work, we study the task in a discrete setting, where input values with some structure expressible as a graph. Namely, the input values are vertices of a finite graph $G$ and the goal is to output vertices that are within distance $d$ of each other in $G$, but still remain in the graph-induced convex hull of the input values. For $d=0$, the task reduces to consensus and cannot be solved with a deterministic algorithm in an asynchronous system even with a single crash fault. For any $d \ge 1$, we show that the task is solvable in asynchronous systems when $G$ is chordal and $n > (\omega+1)f$, where $\omega$ is the clique number of~$G$. In addition, we give the first Byzantine-tolerant algorithm for a variant of lattice agreement. For synchronous systems, we show tight resilience bounds for the exact variants of these and related tasks over a large class of combinatorial structures.


title: Green-Lazarsfeld Condition for Toric Edge Ideals of Bipartite Graphs

abstract: Previously, Ohsugi and Hibi gave a combinatorial description of bipartite graphs $G$ whose toric edge ideal $I_G$ is generated by quadrics, showing that every cycle of $G$ of length at least $6$ must have a chord. This corresponds to the Green-Lazarsfeld condition $\mathbf{N}_1$. In this paper, we investigate the higher syzygies of $I_G$ and give combinatorial descriptions of the Green-Lazarsfeld conditions $\mathbf{N}_p$ of toric edge ideals of bipartite graphs for all $p \ge 1$. In particular, we show that $I_G$ is linearly presented (i.e. satisfies condition $\mathbf{N}_2$) if and only if the bipartite complement of $G$ is a tree of diameter at most $3$. We also investigate the regularity of linearly presented toric edge ideals and give criteria for polyomino ideals to satisfy the Green-Lazarsfeld conditions.


title: Developing a Simple Model for Sand-Tool Interaction and Autonomously Shaping Sand

abstract: Autonomy for robots interacting with sand will enable a wide range of beneficial behaviors, from earth moving for construction and farming vehicles to navigating rough terrain for Mars rovers. The goal of this work is to shape sand into desired forms. Unlike other common autonomous tasks of achieving desired state of a robot, achieving a desired shape of a continuously deformable environment like sand is a much more challenging task. The state of robot can be described with a couple of states-x, y, z, roll, pitch, yaw-but the desired shape of sand can not be described with just a few values. Sand is an aggregation of billions of small particles. After simplifying the model of sand and tool interaction by looking only at the surface of the heightmap, we can formulate the problems into something that is still high dimensional (hundreds to thousands of state dimensions) but much more solvable. We show how this problem can be formulated into a graph search problem and solve it with the A-star algorithm and report preliminary results on using deep reinforcement learning methods like Deep Q-Network and Deep Deterministic Policy Gradient.


title: Measurement of the inclusive isolated-photon cross section in $pp$ collisions at $\sqrt{s}=13$ TeV using 36 fb$^{-1}$ of ATLAS data

abstract: The differential cross section for isolated-photon production in $pp$ collisions is measured at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC using an integrated luminosity of 36.1 fb$^{-1}$. The differential cross section is presented as a function of the photon transverse energy in different regions of photon pseudorapidity. The differential cross section as a function of the absolute value of the photon pseudorapidity is also presented in different regions of photon transverse energy. Next-to-leading-order QCD calculations from JETPHOX and SHERPA as well as next-to-next-to-leading-order QCD calculations from NNLOJET are compared with the measurement, using several parameterisations of the proton parton distribution functions. The predictions provide a good description of the data within the experimental and theoretical uncertainties.


title: Distributed Gradient Descent: Nonconvergence to Saddle Points and the Stable-Manifold Theorem

abstract: The paper studies a distributed gradient descent (DGD) process and considers the problem of showing that in nonconvex optimization problems, DGD typically converges to local minima rather than saddle points. The paper considers unconstrained minimization of a smooth objective function. In centralized settings, the problem of demonstrating nonconvergence to saddle points of gradient descent (and variants) is typically handled by way of the stable-manifold theorem from classical dynamical systems theory. However, the classical stable-manifold theorem is not applicable in distributed settings. The paper develops an appropriate stable-manifold theorem for DGD showing that convergence to saddle points may only occur from a low-dimensional stable manifold. Under appropriate assumptions (e.g., coercivity), this result implies that DGD typically converges to local minima and not to saddle points.


title: Deblending and Classifying Astronomical Sources with Mask R-CNN Deep Learning

abstract: We apply a new deep learning technique to detect, classify, and deblend sources in multi-band astronomical images. We train and evaluate the performance of an artificial neural network built on the Mask R-CNN image processing framework, a general code for efficient object detection, classification, and instance segmentation. After evaluating the performance of our network against simulated ground truth images for star and galaxy classes, we find a precision of 92% at 80% recall for stars and a precision of 98% at 80% recall for galaxies in a typical field with $\sim30$ galaxies/arcmin$^2$. We investigate the deblending capability of our code, and find that clean deblends are handled robustly during object masking, even for significantly blended sources. This technique, or extensions using similar network architectures, may be applied to current and future deep imaging surveys such as LSST and WFIRST. Our code, Astro R-CNN, is publicly available at https://github.com/burke86/astro_rcnn.


title: The largest angle bisection procedure

abstract: The {\it largest angle bisection} procedure is the operation which partitions a given triangle, $T$, into two smaller triangles by constructing the angle bisector of the largest angle of $T$. Applying the procedure to each of these two triangles produces a partition of $T$ into four smaller triangles. Continuing in this manner, after $n$ iterations, the initial triangle is divided into $2^n$ small triangles. We prove that as $n$ approaches infinity, the diameters of all these $2^n$ triangles tend to $0$, the smallest angle of all these triangles is bounded away from $0$, and that, with the exception of $T$ being an isosceles right triangle, the number of dissimilar triangles is unbounded.


title: A physics-informed reinforcement learning approach for the interfacial area transport in two-phase flow

abstract: The prediction of interfacial structure in two-phase flow systems is difficult and challenging. In this paper, a novel physics-informed reinforcement learning-aided framework (PIRLF) for the interfacial area transport is proposed. A Markov Decision Process that describes the bubble transport is established by assuming that the development of two-phase flow is a stochastic process with Markov property. The framework aims to capture the complexity of two-phase flow using the advantage of reinforcement learning (RL) in discovering complex patterns with the help of the physical model (Interfacial Area Transport Equation) as reference. The details of the framework design are described including the design of the environment and the algorithm used in solving the RL problem. The performance of the PIRLF is tested through experiments using the experimental database for vertical upward bubbly air-water flows. The result shows a good performance of PIRLF with rRMSE of 6.556%. The case studies on the PIRLF performance also show that the type of reward function that is related to the physical model can affect the framework performance. Based on the study, the optimal reward function is established. The approaches to extending the capability of PIRLF are discussed, which can be a reference for the further development of this methodology.

categories: physics.comp-ph cs.GT cs.LG physics.flu-dyn


title: Elliptical trajectories of a point on the elliptical 2-sphere

abstract: The focus of this work is to analyze the trajectories of a point on the ellipsoid $\mathbb{S}_{a_{1},a_{2},a_{3}}^{2}$ while it is under the influence of a Killing vector field $K$. For this purpose, we introduce the generalized Darboux frame and the variational vector fields of $\mathbb{S}_{a_{1},a_{2},a_{3}}^{2}$. Then, we determine the Killing equations in terms of the Darboux frame invariants along an ellipsoidal curve. The Killing equations make it possible for us to interpret the magnetic trajectory of a point on the ellipsoid $\mathbb{S}_{a_{1},a_{2},a_{3}}^{2}$. Then, we determine two special trajectories using the variational method. The first one is magnetic curves that are the trajectories produced by the Killing magnetic field $K$ are satisfied the following Lorentz force equation $F_{L} (t)=K\times _{E}t=\nabla _{T}t$, where $\times _{E}$ is elliptical cross product and $\nabla $ is the Levi-Civita connection of the ellipsoid $\mathbb{S}_{a_{1},a_{2},a_{3}}^{2}$. The second one is generalized magnetic helices that are trajectories described by the trajectory of a point on a great ellipse of the ellipsoid rolling without slipping on a fixed ellipse of the ellipsoid using the elliptical motion on the $\mathbb{S}_{a_{1},a_{2},a_{3}}^{2}$. Furthermore, we give various examples and visualized them with the program Mathematica.


title: On maximal multiplicities for Hamiltonians with separable variables

abstract: For $\mathbb N^*:=\mathbb N \setminus \{0\}$, we consider the collection $\mathfrak M(N)$ of all the $N$ rows, for which, for $n=1,\cdots,N$, the $n-th$ row consists of an increasing sequence $(a_j^n)_j$ of real numbers. For $\mathfrak A \in \mathfrak M(N)$, we define its spectrum $\sigma(\mathfrak A)$ by $\sigma(\mathfrak A)=\{\lambda\in \mathbb R \;|\; \lambda=\sum_{n=1}^Na_{j_n}^n\}\,,$ where $(j_1,j_2,\dots,j_N)\in (\mathbb N^*)^N$. This spectrum is discrete and consists of an infinite sequence that can be ordered as a strictly increasing sequence $\lambda_k(\mathfrak A)$. For $\lambda \in \sigma (\mathfrak A)$ we denote by $m(\lambda,\mathfrak A) $ the number of representations of such a $\lambda$, hence the multiplicity of $\lambda$.\\ In this paper we investigate for given $N\in \mathbb N^*$ and $k\in \mathbb N^*$ the highest possible multiplicity (denoted by $\mathfrak m_k(N)$) of $\lambda_k(\mathfrak A)$ for $\mathfrak A \in \mathfrak M(N)$. We give the exact result for $N=2$ and for $N=3$ prove a lower bound which appears, according to numerical experiments, as a "good" conjecture. For the general case, we give examples demonstrating that the problem is quite difficult. \\ This problem is equivalent to the analogue eigenvalue multiplicity questions for Schr\"odinger operators describing a system of N non-interacting one-dimensional particles.


title: Disorder in Twisted Bilayer Graphene

abstract: We develop a theory for a qualitatively new type of disorder in condensed matter systems arising from local twist-angle fluctuations in two strongly coupled van der Waals monolayers twisted with respect to each other to create a flat band moir\'e superlattice. The new paradigm of 'twist angle disorder' arises from the currently ongoing intense research activity in the physics of twisted bilayer graphene. In experimental samples of pristine twisted bilayer graphene, which are nominally free of impurities and defects, the main source of disorder is believed to arise from the unavoidable and uncontrollable non-uniformity of the twist angle across the sample. To address this new physics of twist-angle disorder, we develop a real-space, microscopic model of twisted bilayer graphene where the angle enters as a free parameter. In particular, we focus on the size of single-particle energy gaps separating the miniband from the rest of the spectrum, the Van Hove peaks, the renormalized Dirac cone velocity near charge neutrality, and the minibandwidth. We find that the energy gaps and minibandwidth are strongly affected by disorder while the renormalized velocity remains virtually unchanged. We discuss the implications of our results for the ongoing experiments on twisted bilayer graphene. Our theory is readily generalized to future studies of twist angle disorder effects on all electronic properties of moir\'e superlattices created by twisting two coupled van der Waals materials with respect to each other.


title: Quantum Overlapping Tomography

abstract: It is now experimentally possible to entangle thousands of qubits, and efficiently measure each qubit in parallel in a distinct basis. To fully characterize an unknown entangled state of $n$ qubits, one requires an exponential number of measurements in $n$, which is experimentally unfeasible even for modest system sizes. By leveraging (i) that single-qubit measurements can be made in parallel, and (ii) the theory of perfect hash families, we show that all $k$-qubit reduced density matrices of an $n$ qubit state can be determined with at most $e^{\mathcal{O}(k)} \log^2(n)$ rounds of parallel measurements. We provide concrete measurement protocols which realize this bound. As an example, we argue that with current experiments, the entanglement between every pair of qubits in a system of 1000 qubits could be measured and completely characterized in a few days. This corresponds to completely characterizing entanglement of nearly half a million pairs of qubits.


title: Quantum Penrose Inequality

abstract: The classical Penrose inequality specifies a lower bound on the total mass in terms of the area of certain trapped surfaces. This fails at the semiclassical level. We conjecture a Quantum Penrose Inequality: the mass at spatial infinity is lower-bounded by a function of the generalized entropy of the lightsheet of appropriate quantum trapped surfaces. This is the first relation between quantum information in quantum gravity, and the total energy.


title: Dark Matter, Dark Radiation and Gravitational Waves from Mirror Higgs Parity

abstract: An exact parity replicates the Standard Model giving a Mirror Standard Model, SM $\leftrightarrow$ SM$'$. This "Higgs Parity" and the mirror electroweak symmetry are spontaneously broken by the mirror Higgs, $\left\langle H'\right\rangle = v' \gg \left\langle H\right\rangle$, yielding the Standard Model Higgs as a Pseudo-Nambu-Goldstone Boson of an approximate $SU(4)$ symmetry, with a quartic coupling $\lambda_{SM}(v') \sim 10^{-3}$. Mirror electromagnetism is unbroken and dark matter is composed of $e'$ and $\bar{e}'$. Direct detection may be possible via the kinetic mixing portal, and in unified theories this rate is correlated with the proton decay rate. With a high reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-out followed by dilution from decays of mirror neutrinos, $\nu' \rightarrow \ell H$. Remarkably, this requires $v' \sim (10^8 - 10^{10})$ GeV, consistent with the Higgs mass, and a Standard Model neutrino mass of $(10^{-2} - 10^{-1})$ eV, consistent with observed neutrino masses. The mirror QCD sector exhibits a first order phase transition producing gravitational waves that may be detected by future observations. Mirror glueballs decay to mirror photons giving dark radiation with $\Delta N_{\rm eff} \sim 0.03 - 0.4$. With a low reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-in from the SM sector by either the Higgs or kinetic mixing portal.


title: Intrinsic superconducting instabilities of a solvable model for an incoherent metal

abstract: We construct a family of translationally invariant lattice models with a large number ($N$) of orbitals at every site coupled together via single electron tunneling. By tuning the relative strength of the electronic bandwidth and on-site interactions, that have a modified Sachdev-Ye-Kitaev (SYK) form, we demonstrate a number of unusual features at strong coupling and in the large$-N$ limit. We find examples of (i) an intrinsic non-BCS superconducting instability arising out of an incoherent non-Fermi liquid metal, and, (ii) an instability of an incipient heavy Fermi liquid metal to superconductivity with transition temperatures comparable to its renormalized bandwidth. At strong-coupling, these solvable models display pairing instabilities that are not driven by any special "nesting" properties associated with an underlying Fermi surface.


title: Fractionalized Degrees of Freedom at Infinite Coupling in large Nf QED in 2+1 dimensions

abstract: I consider quantum electrodynamics with many electrons in 2+1 space-time dimensions at finite temperature. The relevant dimensionless interaction parameter for this theory is the fine structure constant divided by the temperature. The theory is solvable at any value of the coupling, in particular for very weak (high temperature) and infinitely strong coupling (corresponding to the zero temperature limit). Concentrating on the photon, each of its physical degrees of freedom at infinite coupling only contributes half of the free-theory value to the entropy. These fractional degrees of freedom are reminiscent of what has been observed in other strongly coupled systems (such as N=4 SYM), and bear similarity to the fractional Quantum Hall effect, potentially suggesting connections between these phenomena. The results found for QED3 are fully consistent with the expectations from particle-vortex duality.


title: Reciprocal skin effect and its realization in a topolectrical circuit

abstract: A system is non-Hermitian when it exchanges energy with its environment and non-reciprocal when it behaves differently upon the interchange of input and response. Within the field of metamaterial research on synthetic topological matter, the skin effect describes the conspiracy of non-Hermiticity and non-reciprocity to yield extensive anomalous localization of all eigenmodes in a (quasi) one-dimensional geometry. Here, we introduce the reciprocal skin effect, which occurs in non-Hermitian but reciprocal systems in two or more dimensions: Eigenmodes with opposite longitudinal momentum exhibit opposite transverse anomalous localization. We experimentally demonstrate the reciprocal skin effect in a passive RLC circuit, suggesting convenient alternative implementations in optical, acoustic, mechanical, and related platforms. Skin mode localization brings forth potential applications in directional and polarization detectors for electromagnetic waves.


title: DoFun 3.0: Functional equations in Mathematica

abstract: We present version 3.0 of the Mathematica package DoFun for the derivation of functional equations. In this version, the derivation of equations for correlation functions of composite operators was added. In the update, the general workflow was slightly modified taking into account experience with the previous version. In addition, various tools were included to improve the usage experience and the code was partially restructured for easier maintenance.


title: Scale Invariant Entanglement Negativity at the Many-Body Localization Transition

abstract: The exact nature of the many-body localization transition remains an open question. An aspect which has been posited in various studies is the emergence of scale invariance around this point, however the direct observation of this phenomenon is still absent. Here we achieve this by studying the logarithmic negativity and mutual information between disjoint blocks of varying size across the many-body localization transition. The two length scales, block sizes and the distance between them, provide a clear quantitative probe of scale invariance across different length scales. We find that at the transition point, the logarithmic negativity obeys a scale invariant exponential decay with respect to the ratio of block separation to size, whereas the mutual information obeys a polynomial decay. The observed scale invariance of the quantum correlations in a microscopic model opens the direction to probe the fractal structure in critical eigenstates using tensor network techniques and provide constraints on the theory of the many-body localization transition.


title: Valence bond solid and possible deconfined quantum criticality in an extended kagome lattice Heisenberg antiferromagnet

abstract: We present numerical evidence for the emergence of an extended valence bond solid (VBS) phase at $T=0$ in the kagome $S=1/2$ Heisenberg antiferromagnet with ferromagnetic further-neighbor interactions. The VBS is located at the boundary between two magnetically ordered regions and extends close to the nearest-neighbor Heisenberg point. It exhibits a diamond-like singlet covering pattern with a $12$-site unit-cell. Our results suggest the possibility of a direct, possibly continuous, quantum phase transition from the neighboring $\mathbf{q}=0$ coplanar magnetically ordered phase into the VBS phase. Moreover, a second phase which breaks lattice symmetries, and is of likely spin-nematic type, is found close to the transition to the ferromagnetic phase. The results have been obtained using numerical Exact Diagonalization. We discuss implications of our results on the nature of nearest-neighbor Heisenberg antiferromagnet.