title: Challenges in Collaborative HRI for Remote Robot Teams

abstract: Collaboration between human supervisors and remote teams of robots is highly challenging, particularly in high-stakes, distant, hazardous locations, such as off-shore energy platforms. In order for these teams of robots to truly be beneficial, they need to be trusted to operate autonomously, performing tasks such as inspection and emergency response, thus reducing the number of personnel placed in harm's way. As remote robots are generally trusted less than robots in close-proximity, we present a solution to instil trust in the operator through a `mediator robot' that can exhibit social skills, alongside sophisticated visualisation techniques. In this position paper, we present general challenges and then take a closer look at one challenge in particular, discussing an initial study, which investigates the relationship between the level of control the supervisor hands over to the mediator robot and how this affects their trust. We show that the supervisor is more likely to have higher trust overall if their initial experience involves handing over control of the emergency situation to the robotic assistant. We discuss this result, here, as well as other challenges and interaction techniques for human-robot collaboration.


title: The Game of Triangles

abstract: We present STR (Star-Triangle Relations), a Mathematica package designed to solve Feynman integrals by means of the method of uniqueness in any Euclidean spacetime dimension. We provide a set of tools to draw Feynman diagrams and interact with them only by the use of the mouse. Throughout the use of a graphic interface, the package should be easily accessible to users with little or no previous experience on diagrams computation.


title: Symmetry-protected metallic and topological phases in penta-materials

abstract: We analyze the symmetry and topological features of a family of materials closely related to penta-graphene, derived from it by adsorption or substitution of different atoms. Our description is based on a novel approach, called topological quantum chemistry, that allows to characterize the topology of the electronic bands, based on the mapping between real and reciprocal space. In particular, by adsorption of alkaline (Li or Na) atoms we obtain a nodal line metal at room temperature, with a continuum of Dirac points around the perimeter of the Brillouin zone. This behavior is also observed in some substitutional derivatives of penta-graphene, such as penta-PC$_2$. Breaking of time-reversal symmetry can be achieved by the use of magnetic atoms; we study penta-MnC$_2$, which also presents spin-orbit coupling and reveals a topological insulator phase. We find that for this family of materials, symmetry is the source of protection for metallic and nontrivial topological phases that can be associated to the presence of fractional band filling, spin-orbit coupling and time-reversal symmetry breaking.


title: Merging versus Ensembling in Multi-Study Prediction: Theoretical Insight from Random Effects

abstract: A critical decision point when training predictors using multiple studies is whether these studies should be combined or treated separately. We compare two multi-study learning approaches in the presence of potential heterogeneity in predictor-outcome relationships across datasets. We consider 1) merging all of the datasets and training a single learner, and 2) multi-study ensembling, which involves training a separate learner on each dataset and combining the predictions resulting from each learner. In a linear regression setting, we show analytically and confirm via simulation that merging yields lower prediction error than ensembling when the predictor-outcome relationships are relatively homogeneous across studies. However, as cross-study heterogeneity increases, there exists a transition point beyond which ensembling outperforms merging. We provide analytic expressions for the transition point in various scenarios, study asymptotic properties, and illustrate how transition point theory can be used for deciding when studies should be combined with an application from metabolomics.


title: Superconvergent HDG methods for Maxwell's equations via the $M$-decomposition

abstract: The concept of the $M$-decomposition was introduced by Cockburn et al.\ in Math. Comp.\ vol.\ 86 (2017), pp.\ 1609-1641 {to provide criteria to guarantee optimal convergence rates for the Hybridizable Discontinuous Galerkin (HDG) method for coercive elliptic problems}. In that paper they systematically constructed superconvergent hybridizable discontinuous Galerkin (HDG) methods to approximate the solutions of elliptic PDEs on unstructured meshes. In this paper, we use the $M$-decomposition to construct HDG methods for the Maxwell's equations on unstructured meshes in two dimension. In particular, we show the any choice of spaces having an $M$-decomposition, together with sufficiently rich auxiliary spaces, has an optimal error estimate and superconvergence even though the problem is not in general coercive. Unlike the elliptic case, we obtain a superconvergent rate for the curl of the solution, not the solution, and this is confirmed by our numerical experiments.


title: Ultra-high Purcell factor using long-range modes in asymmetric plasmonic waveguides

abstract: For integrated optical devices, realistic utilization of the superior wave-matter interaction offered by plasmonics is typically impeded by optical losses, which increase rapidly with mode volume reduction. Although coupled-mode plasmonic structures has demonstrated effective alleviation of the loss-confinement trade-off, stringent symmetry requirements must be enforced for such reduction to prevail. In this work, we report an asymmetric plasmonic waveguide that is not only capable of guiding subwavelength optical mode with long-range propagation, but is also unrestricted by structural, material, or modal symmetry. In these composite hybrid plasmonic waveguides (CHPWs), the versatility afforded by coupling dissimilar plasmonic structures, within the same waveguide, allow better fabrication tolerance and provide more degrees of design freedom to simultaneously optimize various device attributes. The CHPWs used to demonstrated the concept in this work, exhibit propagation loss and mode area of only 0.03 dB/{\mu}m and 0.002 {\mu}m^2 respectively, corresponding to the smallest combination amongst experimentally demonstrated long-range plasmonic structures to-date. Moreover, CHPW micro-rings were realized with record in/out coupling excitation efficiency (71%), extinction ratio (29 dB), as well as Purcell factor (1.5\times10^4).


title: Representation Learning on Visual-Symbolic Graphs for Video Understanding

abstract: Events in natural videos typically arise from spatio-temporal interactions between actors and objects and involve multiple co-occurring activities and object classes. To capture this rich visual and semantic context, we propose using two graphs: (1) an attributed spatio-temporal visual graph whose nodes correspond to actors and objects and whose edges encode different types of interactions, and (2) a symbolic graph that models semantic relationships. We further propose a graph neural network for refining the representations of actors, objects and their interactions on the resulting hybrid graph. Our model goes beyond current approaches that assume nodes and edges are of the same type, operate on graphs with fixed edge weights and do not use a symbolic graph. In particular, our framework: a) has specialized attention-based message functions for different node and edge types; b) uses visual edge features; c) integrates visual evidence with label relationships; and d) performs global reasoning in the semantic space. Experiments on challenging video understanding tasks, such as temporal action localization on the Charades dataset, show that the proposed method leads to state-of-the-art performance.


title: Extreme Covering Systems of the Integers

abstract: It is proved that if the least modulus of a distinct covering system is 4, its largest modulus is at least 60; also if the least modulus is 3, the LCM of the moduli is at least 120; finally, if the least modulus is 4, the LCM of the moduli is at least 360. The constants 60, 120, and 360 can not be replaced by larger constants.


title: POPQORN: Quantifying Robustness of Recurrent Neural Networks

abstract: The vulnerability to adversarial attacks has been a critical issue for deep neural networks. Addressing this issue requires a reliable way to evaluate the robustness of a network. Recently, several methods have been developed to compute $\textit{robustness quantification}$ for neural networks, namely, certified lower bounds of the minimum adversarial perturbation. Such methods, however, were devised for feed-forward networks, e.g. multi-layer perceptron or convolutional networks. It remains an open problem to quantify robustness for recurrent networks, especially LSTM and GRU. For such networks, there exist additional challenges in computing the robustness quantification, such as handling the inputs at multiple steps and the interaction between gates and states. In this work, we propose $\textit{POPQORN}$ ($\textbf{P}$ropagated-$\textbf{o}$ut$\textbf{p}$ut $\textbf{Q}$uantified R$\textbf{o}$bustness for $\textbf{RN}$Ns), a general algorithm to quantify robustness of RNNs, including vanilla RNNs, LSTMs, and GRUs. We demonstrate its effectiveness on different network architectures and show that the robustness quantification on individual steps can lead to new insights.


title: Avalanches, thresholds, and diffusion in meso-scale amorphous plasticity

abstract: We present results on a meso-scale model for amorphous matter in athermal, quasi-static (a-AQS), steady state shear flow. In particular, we perform a careful analysis of the scaling with the lateral system size, $L$, of: i) statistics of individual relaxation events in terms of stress relaxation, $S$, and individual event mean-squared displacement, $M$, and the subsequent load increments, $\Delta \gamma$, required to initiate the next event; ii) static properties of the system encoded by $x=\sigma_y-\sigma$, the distance of local stress values from threshold; and iii) long-time correlations and the emergence of diffusive behavior. For the event statistics, we find that the distribution of $S$ is similar to, but distinct from, the distribution of $M$. We find a strong correlation between $S$ and $M$ for any particular event, with $S\sim M^{q}$ with $q\approx 0.65$. $q$ completely determines the scaling exponents for $P(M)$ given those for $P(S)$. For the distribution of local thresholds, we find $P(x)$ is analytic at $x=0$, and has a value $\left. P(x)\right|_{x=0}=p_0$ which scales with lateral system length as $p_0\sim L^{-0.6}$. Extreme value statistics arguments lead to a scaling relation between the exponents governing $P(x)$ and those governing $P(S)$. Finally, we study the long-time correlations via single-particle tracer statistics. The value of the diffusion coefficient is completely determined by $\langle \Delta \gamma \rangle$ and the scaling properties of $P(M)$ (in particular from $\langle M \rangle$) rather than directly from $P(S)$ as one might have naively guessed. Our results: i) further define the a-AQS universality class, ii) clarify the relation between avalanches of stress relaxation and diffusive behavior, iii) clarify the relation between local threshold distributions and event statistics.


title: Online Distributed Estimation of Principal Eigenspaces

abstract: Principal components analysis (PCA) is a widely used dimension reduction technique with an extensive range of applications. In this paper, an online distributed algorithm is proposed for recovering the principal eigenspaces. We further establish its rate of convergence and show how it relates to the number of nodes employed in the distributed computation, the effective rank of the data matrix under consideration, and the gap in the spectrum of the underlying population covariance matrix. The proposed algorithm is illustrated on low-rank approximation and $\boldsymbol{k}$-means clustering tasks. The numerical results show a substantial computational speed-up vis-a-vis standard distributed PCA algorithms, without compromising learning accuracy.


title: Theoretical and experimental electron diffraction intensity maps for single crystal silicon from an ultrafast source

abstract: Electron diffraction through a thin patterned silicon membrane can be used to create complex spatial modulations in electron distributions by varying the intensity of different reflections using parameters such as crystallographic orientation and wafer thickness, then selecting specific spots in the diffraction plane using apertures. The patterned electron beams can be used to control phase and amplitude of subsequent x-ray emission, enabling novel coherent x-ray methods. The electrons themselves can also be used for femtosecond time resolved diffraction and microscopy. As a first step toward patterned beams, we demonstrate experimentally and through simulation the ability to accurately predict and control diffraction spot intensities. We simulate MeV transmission electron diffraction patterns using the multislice method for various crystallographic orientations of a single crystal Si(100) membrane near beam normal. The resulting intensity maps of the Bragg reflections are compared to experimental results obtained at the Accelerator Structure Test Area Ultrafast Electron Diffraction (ASTA UED) facility at SLAC. Furthermore, the fraction of inelastic and elastic scattering of the initial charge is estimated along with the absorption of the membrane to determine the contrast that would be seen in a patterned version of the Si(100) membrane.


title: Beyond-Luttinger-Liquid thermodynamics of a one-dimensional contact repulsive Bose gas

abstract: We present a thorough study of the thermodynamics of a one-dimensional repulsive Bose gas, focusing in particular on corrections beyond the Luttinger-liquid description. We compute the chemical potential, the pressure and the contact, as a function of temperature and gas parameter with exact thermal Bethe-Ansatz. In addition, we provide interpretations of the main features in the analytically tractable regimes, based on a variety of approaches (Bogoliubov, hard-core, Sommerfeld and virial). The beyond Luttinger-liquid thermodynamic effects are found to be non-monotonic as a function of gas parameter. Such behavior is explained in terms of non-linear dispersion and ``negative excluded volume'' effects, for weak and strong repulsion respectively, responsible for the opposite sign corrections in the thermal next-to-leading term of the thermodynamic quantities at low temperatures. Our predictions can be applied to other systems including super Tonks-Girardeau gases, dipolar and Rydberg atoms, helium, quantum liquid droplets in bosonic mixtures and impurities in a quantum bath.


title: An embedded boundary approach for efficient simulations of viscoplastic fluids in three dimensions

abstract: We present a methodology for simulating three-dimensional flow of incompressible viscoplastic fluids modelled by generalised Newtonian rheological equations. It is implemented in a highly efficient framework for massively parallelisable computations on block-structured grids. In this context, geometric features are handled by the embedded boundary approach, which requires specialised treatment only in cells intersecting or adjacent to the boundary. This constitutes the first published implementation of an embedded boundary algorithm for simulating flow of viscoplastic fluids on structured grids. The underlying algorithm employs a two-stage Runge-Kutta method for temporal discretisation, in which viscous terms are treated semi-implicitly and projection methods are utilised to enforce the incompressibility constraint. We augment the embedded boundary algorithm to deal with the variable apparent viscosity of the fluids. Since the viscosity depends strongly on the strain rate tensor, special care has been taken to approximate the components of the velocity gradients robustly near boundary cells, both for viscous wall fluxes in cut cells and for updates of apparent viscosity in cells adjacent to them. After performing convergence analysis and validating the code against standard test cases, we present the first ever fully three-dimensional simulations of creeping flow of Bingham plastics around translating objects. Our results shed new light on the flow fields around these objects.


title: Efficient color code decoders in $d\geq 2$ dimensions from toric code decoders

abstract: We introduce an efficient decoder of the color code in $d\geq 2$ dimensions, the Restriction Decoder, which uses any $d$-dimensional toric code decoder combined with a local lifting procedure to find a recovery operation. We prove that the Restriction Decoder successfully corrects errors in the color code if and only if the corresponding toric code decoding succeeds. We also numerically estimate the Restriction Decoder threshold for the color code in two and three dimensions against the bit-flip and phase-flip noise with perfect syndrome extraction. We report that the 2D color code threshold $p_{\textrm{2D}} \approx 10.2\%$ on the square-octagon lattice is on a par with the toric code threshold on the square lattice.


title: MiSC: Mixed Strategies Crowdsourcing

abstract: Popular crowdsourcing techniques mostly focus on evaluating workers' labeling quality before adjusting their weights during label aggregation. Recently, another cohort of models regard crowdsourced annotations as incomplete tensors and recover unfilled labels by tensor completion. However, mixed strategies of the two methodologies have never been comprehensively investigated, leaving them as rather independent approaches. In this work, we propose $\textit{MiSC}$ ($\textbf{Mi}$xed $\textbf{S}$trategies $\textbf{C}$rowdsourcing), a versatile framework integrating arbitrary conventional crowdsourcing and tensor completion techniques. In particular, we propose a novel iterative Tucker label aggregation algorithm that outperforms state-of-the-art methods in extensive experiments.


title: Interaction between Northern Coal Sack in the Cyg OB 7 cloud complex and the multiple super nova remnants including HB 21

abstract: We report possible interaction between multiple super nova remnants (SNRs) and Northern Coal Sack (NCS) which is a massive clump (~1000 Mo) in the Cyg OB 7 cloud complex and is forming a massive Class 0 object. We performed molecular observations of the 12CO(J=1-0), 13CO(J=1-0), and C18O(J=1-0) emission lines using the 45m telescope at the Nobeyama Radio Observatory, and we found that there are mainly four velocity components at Vlsr=-20, -6, -4, and 10 km/s. The -6 and -4 km/s components correspond to the systemic velocities of NCS and the Cygnus OB 7 complex, respectively, and the other velocity components originate from distinct smaller clouds. Interestingly, there are apparent correlations and anti-correlations among the spatial distributions of the four components, suggesting that they are physically interacting with one another. On a larger scale, we find that a group of small clouds belonging to the -20 and 10 km/s components are located along two different arcs around some SNRs including HB 21 which has been suggested to be interacting with the Cyg OB 7 cloud complex, and we also find that NCS is located right at the interface of the arcs. The small clouds are likely to be the gas swept up by the stellar wind of the massive stars which created the SNRs. We suggest that the small clouds alined along the two arcs recently encountered NCS and the massive star formation in NCS was triggered by the strong interaction with the small clouds.


title: Maximum Likelihood Estimation of Toric Fano Varieties

abstract: We study the maximum likelihood estimation problem for several classes of toric Fano models. We start by exploring the maximum likelihood degree for all $2$-dimensional Gorenstein toric Fano varieties. We show that the ML degree is equal to the degree of the surface in every case except for the quintic del Pezzo surface with two ordinary double points and provide explicit expressions that allow one to compute the maximum likelihood estimate in closed form whenever the ML degree is less than 5. We then explore the reasons for the ML degree drop using $A$-discriminants and intersection theory. Finally, we show that toric Fano varieties associated to 3-valent phylogenetic trees have ML degree one and provide a formula for the maximum likelihood estimate. We prove it as a corollary to a more general result about the multiplicativity of ML degrees of codimension zero toric fiber products, and it also follows from a connection to a recent result about staged trees.


title: Blockchain Mining Games with Pay Forward

abstract: We study the strategic implications that arise from adding one extra option to the miners participating in the bitcoin protocol. We propose that when adding a block, miners also have the ability to pay forward an amount to be collected by the first miner who successfully extends their branch, giving them the power to influence the incentives for mining. We formulate a stochastic game for the study of such incentives and show that with this added option, smaller miners can guarantee that the best response of even substantially more powerful miners is to follow the expected behavior intended by the protocol designer.


title: Compensation of polarization dependent loss using noiseless amplification and attenuation

abstract: Polarization dependent loss (PDL) is a serious problem that hinders the transfer of polarization qubits through quantum networks. Recently it has been shown that the detrimental effects of PDL on qubit fidelity can be compensated for with the introduction of an additional passive PDL element that rebalances the polarization modes of the transmitted qubit. This procedure works extremely well when the output of the system is postselected on photon detection. However, in cases where the qubit might be needed for further analysis this procedure introduces unwanted vacuum terms into the state. Here we present procedures for the compensation of the effects of PDL using noiseless amplification and attenuation. Each of these techniques introduces a heralding signal into the correction procedure that significantly reduces the vacuum terms in the final state. When detector inefficiency and dark counts are included in the analysis noiseless amplification remains superior, in terms of the fidelity of the final state, to both noiseless attenuation and passive PDL compensation for detector efficiencies greater than 40%.


title: A Compact Millimeter-Wavelength Fourier-Transform Spectrometer

abstract: We have constructed a Fourier-transform spectrometer (FTS) operating between 50 and 330 GHz with minimum volume (355 x260 x64 mm) and weight (13 lbs) while maximizing optical throughput (100 $\mathrm{mm}^2$ sr) and optimizing the spectral resolution (4 GHz). This FTS is designed as a polarizing Martin-Puplett interferometer with unobstructed input and output in which both input polarizations undergo interference. The instrument construction is simple with mirrors milled on the box walls and one motorized stage as the single moving element. We characterize the performance of the FTS, compare the measurements to an optical simulation, and discuss features that relate to details of the FTS design. The simulation is also used to determine the tolerance of optical alignments for the required specifications. We detail the FTS mechanical design and provide the control software as well as the analysis code online.


title: Sparsity-Promoting Optimal Control of Cyber-Physical Systems over Shared Communication Networks

abstract: Recent years have seen several new directions in the design of sparse control of cyber-physical systems (CPSs) driven by the objective of reducing communication cost. One common assumption made in these designs is that the communication happens over a dedicated network. For many practical applications, however, communication must occur over shared networks, leading to two critical design challenges, namely - time-delays in the feedback and fair sharing of bandwidth among users. In this paper, we present a set of sparse H2 control designs under these two design constraints. An important aspect of our design is that the delay itself can be a function of sparsity, which leads to an interesting pattern of trade-offs in the H2 performance. We present three distinct algorithms. The first algorithm preconditions the assignable bandwidth to the network and produces an initial guess for a stabilizing controller. This is followed by our second algorithm, which sparsifies this controller while simultaneously adapting the feedback delay and optimizing the H2 performance using alternating directions method of multipliers (ADMM). The third algorithm extends this approach to a multiple user scenario where optimal number of communication links, whose total sum is fixed, is distributed fairly among users by minimizing the variance of their H2 performances. The problem is cast as a difference-of-convex (DC) program with mixed-integer linear program (MILP) constraints. We provide theorems to prove convergence of these algorithms, followed by validation through numerical simulations.


title: Spin-valley density wave in moir\'e materials

abstract: We introduce and study a minimum two-orbital Hubbard model on a triangular lattice, which captures the key features of both the trilayer ABC-stacked graphene-boron nitride heterostructure and twisted transition metal dichalcogenides in a broad parameter range. Our model comprises first- and second-nearest neighbor hoppings with valley-contrasting flux that accounts for trigonal warping in the band structure. For the strong-coupling regime with one electron per site, we derive a spin-orbital exchange Hamiltonian and find the semiclassical ground state to be a spin-valley density wave. We show that a relatively small second-neighbor exchange interaction is sufficient to stabilize the ordered state against quantum fluctuations. Effects of spin- and valley Zeeman fields as well as thermal fluctuations are also examined.


title: Superconductivity in the $\eta$-carbide-type oxides Zr4Rh2Ox

abstract: We report on the synthesis and the superconductivity of Zr$_4$Rh$_2$O$_{x}$ ($x$ = 0.4, 0.5, 0.6, 0.7, 1.0). These compounds crystallize in the $\eta$-carbide structure, which is a filled version of the complex intermetallic Ti$_2$Ni structure. We find that in the system Zr$_4$Rh$_2$O$_{x}$, already a small amount ($x$ $\geq$ 0.4) of oxygen addition stabilizes the $\eta$-carbide structure over the more common intermetallic CuAl$_2$ structure-type, in which Zr$_2$Rh crystallizes. We show that Zr$_4$Rh$_2$O$_{0.7}$ and Zr$_4$Rh$_2$O are bulk superconductors with critical temperatures of $T_c \approx$ 2.8 K and 4.7 K in the resistivity, respectively. Our analysis of the superconducting properties reveal both compounds to be strongly type-II superconductors with critical fields up to $\mu_0 H_{c1}$(0) $\approx$ 8.8 mT and $\mu_0 H_{c2}$(0) $\approx$ 6.08 T. Our results support that the $\eta$-carbides are a versatile family of compounds for the investigation of the interplay of interstitial doping on physical properties, especially for superconductivity.


title: Relating the curvature of De Sitter Universe to Open Quantum Lamb Shift Spectroscopy

abstract: In this paper, we explore the connection between the curvature of the background De Sitter space-time with the spectroscopic study of entanglement of two atoms. Our set up is in the context of an Open Quantum System (OQS), where the two atoms, each having two energy levels and represented by Pauli spin tensor operators projected along any arbitrary direction. The system mimics the role of a pair of freely falling Unruh De-Witt detectors, which are allowed to non-adiabatically interact with a conformally coupled massless probe scalar field which has the role of background thermal bath. The effective dynamics of this combined system takes into account of the non-adiabatic interaction, which is commonly known as the Resonant Casimir Polder Interaction (RCPI) with the thermal bath. Our analysis reveals that the RCPI of two stable entangled atoms in the quantum vacuum states in OQS depends on the de Sitter space-time curvature relevant to the temperature of the thermal bath felt by the static observer. We also find that, in OQS, RCPI produces a new significant contribution appearing in the effective Hamiltonian of the total system and thermal bath under consideration. We find that the Lamb Shift is characterized by a decreasing inverse square power-law behavior, $L^{-2}$, when inter atomic Euclidean distance, $L$, is much larger than a characteristic length scale, $k$, which is the inverse surface gravity of the background De Sitter space. If the background space-time would have been Minkowskian this shift decreases as, $L^{-1}$, and is independent of temperature. Thus, we establish a connection between the curvature of the De Sitter space-time with the Lamb Shift spectroscopy.


title: Eigenvector of a matrix in $SO_3(\mathbb{R})$

abstract: Let $A=[a_{ij}]\in O_3(\mathbb{R})$. We give several different proofs of the fact that the vector $$ V:=\left[\begin{array}{ccc} \displaystyle \frac{1}{a_{23}+a_{32}} & \displaystyle \frac{1}{a_{13}+a_{31}} & \displaystyle \frac{1}{a_{12}+a_{21}} \end{array}\right]^T, $$ if it exists, is an eigenvector of $A$ corresponding to the eigenvalue $1$.


title: Quantum dissipation in a scalar field theory with gapped momentum states

abstract: Understanding quantum dissipation is important from both theoretical perspective and applications. Here, we show how to describe dissipation in a scalar field theory. We treat dissipation non-perturbatively, represent it by a bilinear term in the Lagrangian and quantize the theory. We find that dissipation promotes a gap in momentum space and reduces the particle energy. As a result, particle mass becomes dressed by dissipation due to self-interaction. The underlying mechanism is similar to that governing the propagation of transverse collective modes in liquids. We discuss the interplay between the dissipative and mass terms, the associated different regimes of field dynamics and the emergence of ultraviolet and infrared cutoffs due to dissipation.


title: First fully kinetic three-dimensional simulation of the AWAKE baseline scenario

abstract: The "Advanced Proton Driven Plasma Wakefield Acceleration Experiment" (AWAKE) aims to accelerate leptons via proton-beam-driven wakefield acceleration. It comprises extensive numerical studies as well as experiments at the CERN laboratory. The baseline scenario incorporates a plasma volume of approximately $62\,\mathrm{cm}^3$. The plasma wavelength is about $1.25\,\mathrm{mm}$ and needs to be adequately resolved, using a minimum of $130$ points per plasma wavelength, in order to accurately reproduce the physics. The baseline scenario incorporates the proton beam micro-bunching, the concurrent non-linear wakefield growth as well as the off-axis electron beam injection, trapping and acceleration. We present results for the first three-dimensional simulation of this baseline scenario with a full model, using a sufficient resolution. The simulation consumed about $22\,\mathrm{Mch}$ of computer resources and scaled up to 32768 cores, thanks to a multitude of adaptions, improvements and optimization of the simulation code PSC. Through this large-scale simulation effort we were able to verify the results of reduced-model simulations as well as identify important novel effects during the electron injection process.


title: Neural Network Astronomy as a New Tool for Observing Bright and Compact Objects

abstract: We propose a new method for solving an important problem of astronomy that arises in observations with ultrahigh-angular-resolution interferometers. This method is based on the application of the theory of artificial neural networks. We propose and compute a multiparameter model for a celestial object like Sgr A*. For this model we have numerically constructed a number of probable images for neural network training. After neural network training on these images, the quality of its operation has been tested on another series of images from the same model. We have proven that a neural network can recognize and classify celestial objects (also obtained from interferometers) virtually no worse than can be done by a human.


title: Functorial Question Answering

abstract: Distributional compositional (DisCo) models are functors that compute the meaning of a sentence from the meaning of its words. We show that DisCo models in the category of sets and relations correspond precisely to relational databases. As a consequence, we get complexity-theoretic reductions from semantics and entailment of a fragment of natural language to evaluation and containment of conjunctive queries, respectively. Finally, we define question answering as an NP-complete problem.


title: Nematic topological semimetal and insulator in magic angle bilayer graphene at charge neutrality

abstract: We report on a fully self-consistent Hartree-Fock calculation of interaction effects on the Moir\'e flat bands of twisted bilayer graphene, assuming that valley U(1) symmetry is respected. We use realistic band structures and interactions and focus on the charge neutrality point, where experiments have variously reported either insulating or semimetallic behavior. Restricting the search to orders for which the valley U(1) symmetry remains unbroken, we find three types of self-consistent solutions with competitive ground state energy (i) insulators that break $C_2 {\mathcal T}$ symmetry, including valley Chern insulators (ii) spin or valley polarized insulators and (iii) rotation $C_3$ symmetry breaking semimetals whose gaplessness is protected by the topology of the Moir\'e flat bands. We find that the relative stability of these states can be tuned by weak strains that break $C_3$ rotation. The nematic semimetal and also, somewhat unexpectedly, the $C_2 {\mathcal T}$ breaking insulators, are stabilized by weak strain. These ground states may be related to the semi-metallic and insulating behaviors seen at charge neutrality, and the sample variability of their observation. We also compare with the results of STM measurements near charge neutrality.


title: Cosmic Inference: Constraining Parameters With Observations and Highly Limited Number of Simulations

abstract: Cosmological probes pose an inverse problem where the measurement result is obtained through observations, and the objective is to infer values of model parameters which characterize the underlying physical system -- our Universe. Modern cosmological probes increasingly rely on measurements of the small-scale structure, and the only way to accurately model physical behavior on those scales, roughly 65 Mpc/h or smaller, is via expensive numerical simulations. In this paper, we provide a detailed description of a novel statistical framework for obtaining accurate parameter constraints by combining observations with a very limited number of cosmological simulations. The proposed framework utilizes multi-output Gaussian process emulators that are adaptively constructed using Bayesian optimization methods. We compare several approaches for constructing multi-output emulators that enable us to take possible inter-output correlations into account while maintaining the efficiency needed for inference. Using Lyman alpha forest flux power spectrum, we demonstrate that our adaptive approach requires considerably fewer --- by a factor of a few in Lyman alpha P(k) case considered here --- simulations compared to the emulation based on Latin hypercube sampling, and that the method is more robust in reconstructing parameters and their Bayesian credible intervals.


title: Subleading Power Factorization with Radiative Functions

abstract: The study of amplitudes and cross sections in the soft and collinear limits allows for an understanding of their all orders behavior, and the identification of universal structures. At leading power soft emissions are eikonal, and described by Wilson lines. Beyond leading power the eikonal approximation breaks down, soft fermions must be added, and soft radiation resolves the nature of the energetic partons from which they were emitted. For both subleading power soft gluon and quark emissions, we use the soft collinear effective theory (SCET) to derive an all orders gauge invariant bare factorization, at both amplitude and cross section level. This yields universal multilocal matrix elements, which we refer to as radiative functions. These appear from subleading power Lagrangians inserted along the lightcone which dress the leading power Wilson lines. The use of SCET enables us to determine the complete set of radiative functions that appear to $\mathcal{O}(\lambda^2)$ in the power expansion, to all orders in $\alpha_s$. For the particular case of event shape observables in $e^+e^-\to$ dijets we derive how the radiative functions contribute to the factorized cross section to $\mathcal{O}(\lambda^2)$.


title: Probing and dressing magnetic impurities in a superconductor

abstract: We propose a method to probe and control the interactions between an ensemble of magnetic impurities in a superconductor via microwave radiation. Our method relies upon the presence of sub-gap Yu-Shiba-Rusinov (YSR) states associated with the impurities. Depending on the sign of the detuning, radiation generates either a ferro- or antiferromagnetic contribution to the exchange interaction. This contribution can bias the statistics of the random exchange constants stemming from the RKKY interaction. Moreover, by measuring the microwave response at the YSR resonance, one gains information about the magnetic order of the impurities. To this end, we estimate the absorption coefficient as well as the achievable strength of the microwave-induced YSR-interactions using off-resonant radiation. The ability to utilize microwave fields to both probe and control impurity spins in a superconducting host may open new paths to studying metallic spin glasses.


title: Holographic integral geometry with time dependence

abstract: We write down Crofton formulas--expressions that compute lengths of spacelike curves in asymptotically AdS$_3$ geometries as integrals over kinematic space--which apply when the curve and/or the background spacetime is time-dependent. Relative to their static predecessor, the time-dependent Crofton formulas display several new features, whose origin is the local null rotation symmetry of the bulk geometry. In pure AdS$_3$ where null rotations are global symmetries, the Crofton formulas simplify and become integrals over the null planes, which intersect the bulk curve.


title: On the Spectrum of Electrons Accelerated in Supernova Remnants

abstract: Using a semi-analytic model of non-linear diffusive shock acceleration, we model the spectrum of cosmic ray (CR) electrons accelerated by supernova remnants (SNRs). Because electrons experience synchrotron losses in the amplified magnetic fields characteristic of SNRs, they exhibit substantially steeper spectra than protons. In particular, we find that the difference between the electron and proton spectral index (power law slope) ranges from 0.1 to 0.4. Our findings must be reckoned with theories of Galactic CR transport, which often assume that electrons and protons are injected with the same slope, and may especially have implications for the observed "positron excess."


title: Power-law Temperature Dependence of the Penetration Depth in a Topological Superconductor due to Surface States

abstract: We study the temperature dependence of the magnetic penetration depth in a 3D topological superconductor (TSC), incorporating the paramagnetic current due to the surface states. A TSC is predicted to host a gapless 2D surface Majorana fluid. In addition to the bulk-dominated London response, we identify a $T^3$ power-law-in-temperature contribution from the surface, valid in the low-temperature limit. Our system is fully gapped in the bulk, and should be compared to bulk nodal superconductivity, which also exhibits power-law behavior. Power-law temperature dependence of the penetration depth can be one indicator of topological superconductivity.


title: Influence of ground motion duration on the structural response at multiple seismic intensity levels

abstract: This paper aims to investigate the effects of motion duration on the structural seismic demands, seeking potential correlations between motion durations and structural responses at several seismic intensity levels. Three seismic intensity levels with 100years, 475years, and 2475years earthquake return periods (RPs) are first considered for correlation computations. Spectrally matched ground motions are employed to isolate the contribution of duration from the effects of ground motion amplitudes and response spectral shape. Four single degree of freedom systems derived from four real reinforced concrete structures are studied, where both degrading and non-degrading equivalent SDOF systems are included for structural modeling. Results show a low positive correlation between motion duration and structural displacement demand, but this correlation increases with an increase in earthquake RP. It is also investigated whether or not this insignificant positive correlation has an impact on the incremental dynamic analysis curves. The spectrally matched ground motions are divided into two distinct groups in this case: short and long duration ground motions. The comparison of incremental dynamic analysis of these two groups at the collapse limit reveals that long-duration ground motions can cause up to a 20 percent decrease in the spectral acceleration demand of considered structural systems.


title: Spectral weight in Chern-Simons theory with symmetry breaking

abstract: We calculate the low-energy spectral weight of a holographic superfluid coupled to a Chern-Simons term in IR radial scaling geometries characterized by a parameter $\eta$. This work was motivated by previous results where an unexpected low-energy spectral weight and a region of instability were seen, both at finite momentum, for the holographic superfluid. We characterize the effect of varying the Chern-Simons coupling $\alpha$ and condensate charge parameter $\zeta$ on these regions supporting low-energy spectral weight or a finite momentum instability. We show that $\eta$, $\alpha$ and $\zeta$ each plays a unique role in shaping these regions. We find a surface $\alpha_{\text{crit}}(\eta, \zeta)$ above which the theory is unstable. In the longitudinal channel we extend our analysis to general dimension $d$. We briefly analyze the Einstein-Maxwell-dilaton theory and find that Fermi shells exist for $d>4$.


title: Event shape and Multiplicity dependence of Freeze-out Scenario and System Thermodynamics in Proton+Proton Collisions at $\sqrt{s}$ = 13 TeV using PYTHIA8

abstract: Recent observations of QGP-like conditions in high-multiplicity pp collisions from ALICE experiment at the LHC warrants an introspection whether to use pp collisions as a baseline measurement to characterize heavy-ion collisions for possible formation of a Quark-Gluon Plasma. A double differential study of the particle spectra and thermodynamics of the produced system as a function of charged-particle multiplicity and transverse spherocity in pp collisions would shed light into the underlying event dynamics. Transverse spherocity, one of the event shape observables, allows to separate the events in terms of jetty and isotropic events. We analyse the identified particle transverse momentum ($p_{\rm T}$) spectra as a function of charged-particle multiplicity and transverse spherocity using Tsallis non-extensive statistics and Boltzmann-Gibbs Blastwave (BGBW) model in pp collisions at $\sqrt{s}$ = 13 TeV using PYTHIA8 event generator. The extracted parameters such as temperature ($T$), radial flow ($\beta$) and non-extensive parameter ($q$) are shown as a function of charged-particle multiplicity for different spherocity classes. We observe that the isotropic events approach to thermal equilibrium while the jetty ones remain far from equilibrium. We argue that, while studying the QGP-like conditions in small systems, one should separate the isotropic events from the spherocity-integrated events, as the production dynamics are different.


title: Dynamic Vision Sensor integration on FPGA-based CNN accelerators for high-speed visual classification

abstract: Deep-learning is a cutting edge theory that is being applied to many fields. For vision applications the Convolutional Neural Networks (CNN) are demanding significant accuracy for classification tasks. Numerous hardware accelerators have populated during the last years to improve CPU or GPU based solutions. This technology is commonly prototyped and tested over FPGAs before being considered for ASIC fabrication for mass production. The use of commercial typical cameras (30fps) limits the capabilities of these systems for high speed applications. The use of dynamic vision sensors (DVS) that emulate the behavior of a biological retina is taking an incremental importance to improve this applications due to its nature, where the information is represented by a continuous stream of spikes and the frames to be processed by the CNN are constructed collecting a fixed number of these spikes (called events). The faster an object is, the more events are produced by DVS, so the higher is the equivalent frame rate. Therefore, these DVS utilization allows to compute a frame at the maximum speed a CNN accelerator can offer. In this paper we present a VHDL/HLS description of a pipelined design for FPGA able to collect events from an Address-Event-Representation (AER) DVS retina to obtain a normalized histogram to be used by a particular CNN accelerator, called NullHop. VHDL is used to describe the circuit, and HLS for computation blocks, which are used to perform the normalization of a frame needed for the CNN. Results outperform previous implementations of frames collection and normalization using ARM processors running at 800MHz on a Zynq7100 in both latency and power consumption. A measured 67% speedup factor is presented for a Roshambo CNN real-time experiment running at 160fps peak rate.


title: Engineering First-Order Quantum Phase Transitions for Weak Signal Detection

abstract: The quantum critical detector (QCD), recently introduced for weak-signal amplification [Opt. Express 27, 10482 (2019)], functions by exploiting high sensitivity near the phase transition point of first-order quantum phase transitions. We contrast the behavior of the first-order as well as the second-order quantum phase transitions (QPTs) in the detector. We find that the giant sensitivity to a weak input signal, which can be utilized for quantum amplification, only exists in first-order QPTs. We define two new magnetic order parameters to quantitatively characterize the first-order QPT of the interacting spins in the detector. We also introduce the Husimi $Q$-functions as a powerful tool to show the fundamental change in the ground-state wave function of the detector during the QPTs and especially, the intrinsic dynamical change within the detector during a quantum critical amplification. We explicitly show the high figures of merit of the QCD via the quantum gain and signal-to-quantum noise ratio. Specifically, we predict the existence of a universal first-order QPT in the interacting spin system resulting from two competing ferromagnetic orders. Our results motivate new designs of weak signal detectors by engineering first-order QPTs, which are of fundamental significance in the search for new particles, quantum metrology, and information science.


title: Boltzmann stochastic thermodynamics

abstract: The Boltzmann kinetic equation is obtained from an integro-differential master equation that describes a stochastic dynamics in phase space of an isolated thermodynamic system. The stochastic evolution yields a generation of entropy, leading to an increase of Gibbs entropy, in contrast to a Hamiltonian dynamics, described by the Liouville equation, for which the entropy is constant in time. By considering transition rates corresponding to collisions of two particles, the Boltzmann equation is attained. When the angle of the scattering produced by collisions is small, the master equation is shown to be reduced to a differential equation of the Fokker-Planck type. When the dynamics is of the Hamiltonian type, the master equation reduces to the Liouville equation. The present approach is understood as a stochastic interpretation of the reasonings employed by Maxwell and Boltzmann in the kinetic theory of gases regarding the microscopic time evolution.


title: Supernova Photometric Classification Pipelines Trained on Spectroscopically Classified Supernovae from the Pan-STARRS1 Medium-Deep Survey

abstract: Photometric classification of supernovae (SNe) is imperative as recent and upcoming optical time-domain surveys, such as the Large Synoptic Survey Telescope (LSST), overwhelm the available resources for spectrosopic follow-up. Here we develop a range of light curve classification pipelines, trained on 518 spectroscopically-classified SNe from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS): 357 Type Ia, 93 Type II, 25 Type IIn, 21 Type Ibc, and 17 Type I SLSNe. We present a new parametric analytical model that can accommodate a broad range of SN light curve morphologies, including those with a plateau, and fit this model to data in four PS1 filters (griz). We test a number of feature extraction methods, data augmentation strategies, and machine learning algorithms to predict the class of each SN. Our best pipelines result in 90% average accuracy, 70% average purity, and 80% average completeness for all SN classes, with the highest success rates for Type Ia SNe and SLSNe and the lowest for Type Ibc SNe. Despite the greater complexity of our classification scheme, the purity of our Type Ia SN classification, 95%, is on par with methods developed specifically for Type Ia versus non-Type Ia binary classification. As the first of its kind, this study serves as a guide to developing and training classification algorithms for a wide range of SN types with a purely empirical training set, particularly one that is similar in its characteristics to the expected LSST main survey strategy. Future work will implement this classification pipeline on ~3000 PS1/MDS light curves that lack spectroscopic classification.


title: Strain-Based Room-Temperature Non-Volatile MoTe$_2$ Ferroelectric Phase Change Transistor

abstract: The primary mechanism of operation of almost all transistors today relies on electric-field effect in a semiconducting channel to tune its conductivity from the conducting 'on'-state to a non-conducting 'off'-state. As transistors continue to scale down to increase computational performance, physical limitations from nanoscale field-effect operation begin to cause undesirable current leakage that is detrimental to the continued advancement of computing. Using a fundamentally different mechanism of operation, we show that through nanoscale strain engineering with thin films and ferroelectrics (FEs) the transition metal dichalcogenide (TMDC) MoTe$_2$ can be reversibly switched with electric-field induced strain between the 1T'-MoTe$_2$ (semimetallic) phase to a semiconducting MoTe$_2$ phase in a field effect transistor geometry. This alternative mechanism for transistor switching sidesteps all the static and dynamic power consumption problems in conventional field-effect transistors (FETs). Using strain, we achieve large non-volatile changes in channel conductivity (G$_{on}$/G$_{off}$~10$^7$ vs. G$_{on}$/G$_{off}$~0.04 in the control device) at room temperature. Ferroelectric devices offer the potential to reach sub-ns nonvolatile strain switching at the attojoule/bit level, having immediate applications in ultra-fast low-power non-volatile logic and memory while also transforming the landscape of computational architectures since conventional power, speed, and volatility considerations for microelectronics may no longer exist.


title: Galaxy Zoo: Probabilistic Morphology through Bayesian CNNs and Active Learning

abstract: We use Bayesian convolutional neural networks and a novel generative model of Galaxy Zoo volunteer responses to infer posteriors for the visual morphology of galaxies. Bayesian CNN can learn from galaxy images with uncertain labels and then, for previously unlabelled galaxies, predict the probability of each possible label. Our posteriors are well-calibrated (e.g. for predicting bars, we achieve coverage errors of 11.8% within a vote fraction deviation of 0.2) and hence are reliable for practical use. Further, using our posteriors, we apply the active learning strategy BALD to request volunteer responses for the subset of galaxies which, if labelled, would be most informative for training our network. We show that training our Bayesian CNNs using active learning requires up to 35-60% fewer labelled galaxies, depending on the morphological feature being classified. By combining human and machine intelligence, Galaxy Zoo will be able to classify surveys of any conceivable scale on a timescale of weeks, providing massive and detailed morphology catalogues to support research into galaxy evolution.


title: Abelian dualities of $\mathcal{N}=(0,4)$ boundary conditions

abstract: We propose dual pairs of $\mathcal{N}=(0,4)$ half-BPS boundary conditions for 3d $\mathcal{N}=4$ Abelian gauge theories related to mirror symmetry and S-duality by showing the matching of boundary 't Hooft anomalies and supersymmetric indices. We find simple $\mathcal{N}=(0,4)$ mirror symmetry between 2d $\mathcal{N}=(0,4)$ Abelian gauge theories and free Fermi multiplets that generalizes $\mathcal{N}=(0,2)$ Abelian duality. We also propose a prescription for computing half-index of enriched Neumann boundary condition including 2d boundary bosonic matters by gauging the 2d boundary flavor symmetry of Dirichlet boundary condition. By coupling $\mathcal{N}=(0,4)$ half-BPS boundary configurations of 3d $\mathcal{N}=4$ gauge theories to quarter-BPS corner configurations of 4d $\mathcal{N}=4$ Super Yang-Mills theories, we further obtain a new type of 4d-3d-2d duality that may involve 3d non-Abelian gauge symmetry.


title: Error analysis for a fractional-derivative parabolic problem on quasi-graded meshes using barrier functions

abstract: An initial-boundary value problem with a Caputo time derivative of fractional order $\alpha\in(0,1)$ is considered, solutions of which typically exhibit a singular behaviour at an initial time. For this problem, we give a simple and general numerical-stability analysis using barrier functions, which yields sharp pointwise-in-time error bounds on quasi-graded temporal meshes with arbitrary degree of grading. L1-type and Alikhanov-type discretization in time are considered. In particular, those results imply that milder (compared to the optimal) grading yields optimal convergence rates in positive time. Semi-discretizations in time and full discretizations are addressed. The theoretical findings are illustrated by numerical experiments.


title: Multilinear Time Invariant System Theory

abstract: In biological and engineering systems, structure, function and dynamics are highly coupled. Such interactions can be naturally and compactly captured via tensor based state space dynamic representations. However, such representations are not amenable to the standard system and controls framework which requires the state to be in the form of a vector. In order to address this limitation, recently a new class of multiway dynamical systems has been introduced in which the states, inputs and outputs are tensors. We propose a new form of multilinear time invariant (MLTI) systems based on the Einstein product and even-order paired tensors. We extend classical linear time invariant (LTI) system notions including stability, reachability and observability for the new MLTI system representation by leveraging recent advances in tensor algebra.


title: An exact algorithm for biobjective integer programming problems

abstract: We propose an exact algorithm for solving biobjective integer programming problems, which arise in various applications of operations research. The algorithm is based on solving Pascoletti-Serafini scalarizations to search specified regions (boxes) in the objective space and returns the set of nondominated points. We implement the algorithm with different strategies, where the choices of the scalarization model parameters and splitting rule differ. We then derive bounds on the number of scalarization models solved; and demonstrate the performances of the variants through computational experiments both as exact algorithms and as solution approaches under time restriction. The experiments demonstrate that different strategies have advantages in different aspects: while some are quicker in finding the whole set of nondominated solutions, others return good-quality solutions in terms of representativeness when run under time restriction. We also compare the proposed approach with existing algorithms. The results of our experiments show the satisfactory behaviour of our algorithm, especially when run under time limit, as it achieves better coverage of the whole frontier with a smaller number of solutions compared to the existing algorithms.


title: T-structures and twisted complexes on derived injectives

abstract: In the paper "Deformation theory of abelian categories", the last two authors proved that an abelian category with enough injectives can be reconstructed as the category of finitely presented modules over the category of its injective objects. We show a generalization of this to pretriangulated dg-categories with a left bounded non-degenerate t-structure with enough derived injectives, the latter being derived enhancements of the injective objects in the heart of the t-structure. Such dg-categories (with an additional hypothesis of closure under suitable products) can be completely described in terms of left bounded twisted complexes of their derived injectives.


title: Modified transverse Ising model for the dielectric properties of SrTiO$_3$ films and interfaces

abstract: The transverse Ising model (TIM), with pseudospins representing the lattice polarization, is often used as a simple description of ferroelectric materials. However, we demonstrate that the TIM, as it is usually formulated, provides an incorrect description of SrTiO$_{3}$ films and interfaces because of its inadequate treatment of spatial inhomogeneity. We correct this deficiency by adding a pseudospin anisotropy to the model. We demonstrate the physical need for this term by comparison of the TIM to a typical Landau-Ginzburg-Devonshire model. We then demonstrate the physical consequences of the modification for two model systems: a ferroelectric thin film, and a metallic LaAlO$_{3}$/SrTiO$_{3}$ interface. We show that, in both cases, the modified TIM has a substantially different polarization profile than the conventional TIM. In particular, at low temperatures the formation of quantized states at LaAlO$_{3}$/SrTiO$_{3}$ interfaces only occurs in the modified TIM.


title: Rational Groups and a Characterization of a Class of Permutation Groups

abstract: We prove that a finite group is rational if and only if it has a set of permutation characters which separate conjugacy classes. It follows from this that a finite group is rational if and only if it has a representation as a permutation group in which any two elements fixing the same number of letters are conjugate.


title: Evaluation of 4D Light Field Compression Methods

abstract: Light field data records the amount of light at multiple points in space, captured e.g. by an array of cameras or by a light-field camera that uses microlenses. Since the storage and transmission requirements for such data are tremendous, compression techniques for light fields are gaining momentum in recent years. Although plenty of efficient compression formats do exist for still and moving images, only a little research on the impact of these methods on light field imagery is performed. In this paper, we evaluate the impact of state-of-the-art image and video compression methods on quality of images rendered from light field data. The methods include recent video compression standards, especially AV1 and XVC finalised in 2018. To fully exploit the potential of common image compression methods on four-dimensional light field imagery, we have extended these methods into three and four dimensions. In this paper, we show that the four-dimensional light field data can be compressed much more than independent still images while maintaining the same visual quality of a perceived picture. We gradually compare the compression performance of all image and video compression methods, and eventually answer the question, "What is the best compression method for light field data?".


title: Critical Network Cascades with Re-excitable nodes: Why tree-like approximations usually work, when they breakdown, and how to correct them

abstract: Network science is a rapidly expanding field, with a large and growing body of work on network-based dynamical processes. Most theoretical results in this area rely on the so-called \emph{locally tree-like approximation}. This is, however, usually an `uncontrolled' approximation, in the sense that the magnitudes of the error are typically unknown, although numerical results show that this error is often surprisingly small. In this paper, we place this approximation on more rigorous footing by calculating the magnitude of deviations away from tree-based theories in the context of discrete-time critical network cascades with re-excitable nodes. We discuss the conditions under which tree-like approximations give good results for calculating network criticality, and also explain the reasons for deviation from this approximation, in terms of the density of certain kinds of network motifs. Using this understanding, we derive results for network criticality that apply to general networks that explicitly do not satisfy the locally tree-like approximation. In particular, we focus on the bi-parallel motif, the smallest motif relevant to the failure of a tree-based theory in this context, and we derive the corrections due to such motifs on the conditions for criticality. We verify our claims on computer-generated networks, and we confirm that our theory accurately predicts the observed deviations from criticality. Using our theory, we explain why numerical simulations often show that deviations from a tree-based theory are surprisingly small. More specifically, we show that these deviations are negligible for networks whose average degree is even modestly large compared to one, justifying why tree-based theories appear to work well for most real-world networks.


title: A Coordinated Search Strategy for Multiple Solitary Robots: An Extension

abstract: The problem of coordination without a priori information about the environment is important in robotics. Applications vary from formation control to search and rescue. This paper considers the problem of search by a group of solitary robots: self-interested robots without a priori knowledge about each other, and with restricted communication capacity. When the capacity of robots to communicate is limited, they may obliviously search in overlapping regions (i.e. be subject to interference). Interference hinders robot progress, and strategies have been proposed in the literature to mitigate interference [1], [2]. Interaction of solitary robots has attracted much interest in robotics, but the problem of mitigating interference when time for search is limited remains an important area of research. We propose a coordination strategy based on the method of cellular decomposition [3] where we employ the concept of soft obstacles: a robot considers cells assigned to other robots as obstacles. The performance of the proposed strategy is demonstrated by means of simulation experiments. Simulations indicate the utility of the strategy in situations where a known upper bound on the search time precludes search of the entire environment.


title: Alpha MAML: Adaptive Model-Agnostic Meta-Learning

abstract: Model-agnostic meta-learning (MAML) is a meta-learning technique to train a model on a multitude of learning tasks in a way that primes the model for few-shot learning of new tasks. The MAML algorithm performs well on few-shot learning problems in classification, regression, and fine-tuning of policy gradients in reinforcement learning, but comes with the need for costly hyperparameter tuning for training stability. We address this shortcoming by introducing an extension to MAML, called Alpha MAML, to incorporate an online hyperparameter adaptation scheme that eliminates the need to tune meta-learning and learning rates. Our results with the Omniglot database demonstrate a substantial reduction in the need to tune MAML training hyperparameters and improvement to training stability with less sensitivity to hyperparameter choice.


title: A Dynamical Systems Perspective on Nesterov Acceleration

abstract: We present a dynamical system framework for understanding Nesterov's accelerated gradient method. In contrast to earlier work, our derivation does not rely on a vanishing step size argument. We show that Nesterov acceleration arises from discretizing an ordinary differential equation with a semi-implicit Euler integration scheme. We analyze both the underlying differential equation as well as the discretization to obtain insights into the phenomenon of acceleration. The analysis suggests that a curvature-dependent damping term lies at the heart of the phenomenon. We further establish connections between the discretized and the continuous-time dynamics.


title: BEST potential in testing the eV-scale sterile neutrino explanation of reactor antineutrino anomalies

abstract: Baksan Experiment on Sterile Neutrino (BEST) is presently at the stage of production of the artificial neutrino source Cr-51, the gallium exposure will start in July and proceed for three months. While aiming specifically at investigating the Gallium neutrino anomaly (SAGE and GALLEX experiments), BEST can do more and it is tempting to estimate its ability in testing sterile neutrino explanation of antineutrino (reactor) anomalies. We observe a moderate sensitivity to the region in model parameter space (sterile neutrino mass and mixing with active electron neutrino) outlined by the old reactor antineutrino anomaly and the best fit of DANSS experiment, while the Neutrino-4 favorite region falls right in the BEST ballpark. In particular, by analyzing SAGE+GALLEX and Neutrino-4 $\chi^2$ distributions we find that Neutrino-4 results are fully consistent with the Gallium anomaly, the significance of the combined anomaly almost reaches 4$\sigma$ level. If the BEST confirms the Neutrino-4 results, the joint analysis will indicate more than 5$\sigma$ evidence for the sterile neutrino of eV-scale mass.


title: A Fast and Scalable Implementation Method for Competing Risks Data with the R Package fastcmprsk

abstract: Advancements in medical informatics tools and high-throughput biological experimentation make large-scale biomedical data routinely accessible to researchers. Competing risks data are typical in biomedical studies where individuals are at risk to more than one cause (type of event) which can preclude the others from happening. The Fine-Gray model is a popular and well-appreciated model for competing risks data and is currently implemented in a number of statistical software packages. However, current implementations are not computationally scalable for large-scale competing risks data. We have developed an R package, fastcmprsk, that uses a novel forward-backward scan algorithm to significantly reduce the computational complexity for parameter estimation by exploiting the structure of the subject-specific risk sets. Numerical studies compare the speed and scalability of our implementation to current methods for unpenalized and penalized Fine-Gray regression and show impressive gains in computational efficiency.


title: Randomization of Approximate Bilinear Computation for Matrix Multiplication

abstract: We present a method for randomizing formulas for bilinear computation of matrix products. We consider the implications of such randomization when there are two sources of error: One due to the formula itself only being approximately correct, and one due to using floating point arithmetic. Our theoretical results and numerical experiments indicate that our method can improve performance when each of these error sources are present individually, as well as when they are present at the same time.


title: Accelerating lattice quantum Monte Carlo simulation using artificial neural networks: an application to the Holstein model

abstract: Monte Carlo (MC) simulations are essential computational approaches with widespread use throughout all areas of science. We present a method for accelerating lattice MC simulations using fully connected and convolutional artificial neural networks that are trained to perform local and global moves in configuration space, respectively. Both networks take local spacetime MC configurations as input features and can, therefore, be trained using samples generated by conventional MC runs on smaller lattices before being utilized for simulations on larger systems. This new approach is benchmarked for the case of determinant quantum Monte Carlo (DQMC) studies of the two-dimensional Holstein model. We find that both artificial neural networks are capable of learning an unspecified effective model that accurately reproduces the MC configuration weights of the original Hamiltonian and achieve an order of magnitude speedup over the conventional DQMC algorithm. Our approach is broadly applicable to many classical and quantum lattice MC algorithms.


title: Peak sets and boundary interpolation sets for the unit disc: a survey

abstract: This paper surveys results concerning peak sets and boundary interpolation sets for the unit disc. It includes hitherto unpublished results proved by David C. Ullrich on peak sets for the Zygmund class.


title: Transport-Based Neural Style Transfer for Smoke Simulations

abstract: Artistically controlling fluids has always been a challenging task. Optimization techniques rely on approximating simulation states towards target velocity or density field configurations, which are often handcrafted by artists to indirectly control smoke dynamics. Patch synthesis techniques transfer image textures or simulation features to a target flow field. However, these are either limited to adding structural patterns or augmenting coarse flows with turbulent structures, and hence cannot capture the full spectrum of different styles and semantically complex structures. In this paper, we propose the first Transport-based Neural Style Transfer (TNST) algorithm for volumetric smoke data. Our method is able to transfer features from natural images to smoke simulations, enabling general content-aware manipulations ranging from simple patterns to intricate motifs. The proposed algorithm is physically inspired, since it computes the density transport from a source input smoke to a desired target configuration. Our transport-based approach allows direct control over the divergence of the stylization velocity field by optimizing incompressible and irrotational potentials that transport smoke towards stylization. Temporal consistency is ensured by transporting and aligning subsequent stylized velocities, and 3D reconstructions are computed by seamlessly merging stylizations from different camera viewpoints.


title: AutoDispNet: Improving Disparity Estimation With AutoML

abstract: Much research work in computer vision is being spent on optimizing existing network architectures to obtain a few more percentage points on benchmarks. Recent AutoML approaches promise to relieve us from this effort. However, they are mainly designed for comparatively small-scale classification tasks. In this work, we show how to use and extend existing AutoML techniques to efficiently optimize large-scale U-Net-like encoder-decoder architectures. In particular, we leverage gradient-based neural architecture search and Bayesian optimization for hyperparameter search. The resulting optimization does not require a large-scale compute cluster. We show results on disparity estimation that clearly outperform the manually optimized baseline and reach state-of-the-art performance.


title: Percival: Making In-Browser Perceptual Ad Blocking Practical With Deep Learning

abstract: In this paper we present Percival, a browser-embedded, lightweight, deep learning-powered ad blocker. Percival embeds itself within the browser's image rendering pipeline, which makes it possible to intercept every image obtained during page execution and to perform blocking based on applying machine learning for image classification to flag potential ads. Our implementation inside both Chromium and Brave browsers shows only a minor rendering performance overhead of 4.55%, demonstrating the feasibility of deploying traditionally heavy models (i.e. deep neural networks) inside the critical path of the rendering engine of a browser. We show that our image-based ad blocker can replicate EasyList rules with an accuracy of 96.76%. To show the versatility of the Percival's approach we present case studies that demonstrate that Percival 1) does surprisingly well on ads in languages other than English; 2) Percival also performs well on blocking first-party Facebook ads, which have presented issues for other ad blockers. Percival proves that image-based perceptual ad blocking is an attractive complement to today's dominant approach of block lists


title: Search for anomalous electroweak production of vector boson pairs in association with two jets in proton-proton collisions at 13 TeV

abstract: A search for anomalous electroweak production of WW, WZ, and ZZ boson pairs in association with two jets in proton-proton collisions at $\sqrt{s} =$ 13 TeV at the LHC is reported. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ collected with the CMS detector. Events are selected by requiring two jets with large rapidity separation and invariant mass, one or two leptons (electrons or muons), and a W or Z boson decaying hadronically. No excess of events with respect to the standard model background predictions is observed and constraints on the structure of quartic vector boson interactions in the framework of dimension-8 effective field theory operators are reported. Stringent limits on parameters of the effective field theory operators are obtained. The observed 95% confidence level limits for the S0, M0, and T0 operators are $-$2.7 $<$ f$_{\mathrm{S0}}/ \Lambda^{4}$ $<$ 2.7, $-$1.0 $<$ f$_{\mathrm{M0}}/ \Lambda^{4}$ $<$ 1.0, and $-$0.17 $<$ f$_{\mathrm{T0}}/ \Lambda^{4}$ $<$ 0.16, in units of TeV$^{-4}$. Constraints are also reported on the product of the cross section and branching fraction for vector boson fusion production of charged Higgs bosons as a function of mass from 600 to 2000 GeV. The results are interpreted in the context of the Georgi-Machacek model.


title: Limitations and Biases in Facial Landmark Detection -- An Empirical Study on Older Adults with Dementia

abstract: Accurate facial expression analysis is an essential step in various clinical applications that involve physical and mental health assessments of older adults (e.g. diagnosis of pain or depression). Although remarkable progress has been achieved toward developing robust facial landmark detection methods, state-of-the-art methods still face many challenges when encountering uncontrolled environments, different ranges of facial expressions, and different demographics of the population. A recent study has revealed that the health status of individuals can also affect the performance of facial landmark detection methods on front views of faces. In this work, we investigate this matter in a much greater context using seven facial landmark detection methods. We perform our evaluation not only on frontal faces but also on profile faces and in various regions of the face. Our results shed light on limitations of the existing methods and challenges of applying these methods in clinical settings by indicating: 1) a significant difference between the performance of state-of-the-art when tested on the profile or frontal faces of individuals with vs. without dementia; 2) insights on the existing bias for all regions of the face; and 3) the presence of this bias despite re-training/fine-tuning with various configurations of six datasets.


title: REPLAB: A Reproducible Low-Cost Arm Benchmark Platform for Robotic Learning

abstract: Standardized evaluation measures have aided in the progress of machine learning approaches in disciplines such as computer vision and machine translation. In this paper, we make the case that robotic learning would also benefit from benchmarking, and present the "REPLAB" platform for benchmarking vision-based manipulation tasks. REPLAB is a reproducible and self-contained hardware stack (robot arm, camera, and workspace) that costs about 2000 USD, occupies a cuboid of size 70x40x60 cm, and permits full assembly within a few hours. Through this low-cost, compact design, REPLAB aims to drive wide participation by lowering the barrier to entry into robotics and to enable easy scaling to many robots. We envision REPLAB as a framework for reproducible research across manipulation tasks, and as a step in this direction, we define a template for a grasping benchmark consisting of a task definition, evaluation protocol, performance measures, and a dataset of 92k grasp attempts. We implement, evaluate, and analyze several previously proposed grasping approaches to establish baselines for this benchmark. Finally, we also implement and evaluate a deep reinforcement learning approach for 3D reaching tasks on our REPLAB platform. Project page with assembly instructions, code, and videos: https://goo.gl/5F9dP4.


title: Shortest Path Algorithms between Theory and Practice

abstract: Utilizing graph algorithms is a common activity in computer science. Algorithms that perform computations on large graphs are not always efficient. This work investigates the Single-Source Shortest Path (SSSP) problem, which is considered to be one of the most important and most studied graph problems. This thesis contains a review of the SSSP problem in both theory and practice. In addition, it discusses a new single-source shortest-path algorithm that achieves the same $O(n \cdot m)$ time bound as the traditional Bellman-Ford-Moore algorithm but outperforms it and other state-of-the-art algorithms in practice. The work is comprised of three parts. The first discusses some basic shortest-path and negative-cycle-detection algorithms in literature from the theoretical and practical point of view. The second contains a discussion of a new algorithm for the single-source shortest-path problem that outperforms most state-of-the-art algorithms for several well-known families of graphs. The main idea behind the proposed algorithm is to select the fewest most-effective vertices to scan. We also propose a discussion of correctness, termination, and the proof of the worst-case time bound of the proposed algorithm. This section also suggests two different implementations for the proposed algorithm, the first runs faster while the second performs a fewer number of operations. Finally, an extensive computational study of the different shortest paths algorithms is conducted. The results are proposed using a new evaluation metric for shortest-path algorithms. A discussion of outcomes, strengths, and weaknesses of the various shortest path algorithms are also included in this work.


title: Universal thermodynamic bounds on nonequilibrium response with biochemical applications

abstract: Diverse physical systems are characterized by their response to small perturbations. Near thermodynamic equilibrium, the fluctuation-dissipation theorem provides a powerful theoretical and experimental tool to determine the nature of response by observing spontaneous equilibrium fluctuations. In this spirit, we derive here a collection of equalities and inequalities valid arbitrarily far from equilibrium that constrain the response of nonequilibrium steady states in terms of the strength of nonequilibrium driving. Our work opens new avenues for characterizing nonequilibrium response. As illustrations, we show how our results rationalize the energetic requirements of two common biochemical motifs.


title: Transport and Interface: an Uncertainty Principle for the Wasserstein distance

abstract: Let $f: [0,1]^d \rightarrow \mathbb{R}$ be a continuous function with zero mean and interpret $f_{+} = \max(f, 0)$ and $f_{-} = -\min(f, 0)$ as the densities of two measures. We prove that if the cost of transport from $f_{+}$ to $f_{-}$ is small (in terms of the Wasserstein distance $W^1$), then the nodal set $\left\{x \in (0,1)^d: f(x) = 0 \right\}$ has to be large (`if it is always easy to buy milk, there must be many supermarkets'). More precisely, we show that $$ W_1(f_+, f_-) \cdot \mathcal{H}^{d-1}\left\{x \in (0,1)^d: f(x) = 0 \right\} \gtrsim_{d} \left( \frac{\|f\|_{L^1}}{\|f\|_{L^{\infty}}} \right)^{4 - \frac1d} \|f\|_{L^1} \, .$$ We apply this ``uncertainty principle" to the metric Sturm-Liouville theory in higher dimensions to show that a linear combination of eigenfunctions of an elliptic operator cannot have an arbitrarily small zero set.


title: Graph-based Semi-Supervised & Active Learning for Edge Flows

abstract: We present a graph-based semi-supervised learning (SSL) method for learning edge flows defined on a graph. Specifically, given flow measurements on a subset of edges, we want to predict the flows on the remaining edges. To this end, we develop a computational framework that imposes certain constraints on the overall flows, such as (approximate) flow conservation. These constraints render our approach different from classical graph-based SSL for vertex labels, which posits that tightly connected nodes share similar labels and leverages the graph structure accordingly to extrapolate from a few vertex labels to the unlabeled vertices. We derive bounds for our method's reconstruction error and demonstrate its strong performance on synthetic and real-world flow networks from transportation, physical infrastructure, and the Web. Furthermore, we provide two active learning algorithms for selecting informative edges on which to measure flow, which has applications for optimal sensor deployment. The first strategy selects edges to minimize the reconstruction error bound and works well on flows that are approximately divergence-free. The second approach clusters the graph and selects bottleneck edges that cross cluster-boundaries, which works well on flows with global trends.


title: The generalized Hadamard product of polynomials and its stability

abstract: For two polynomials of degrees $n$ and $m$ ($n\geq m$) $$ f\left( s\right) =a_{0}+a_{1}s+\ldots+a_{n-1}s^{n-1}+a_{n}s^{n}$$ and $$g\left( s\right) =b_{0}+b_{1}s+\ldots+b_{m-1}s^{m-1}+b_{m}s^{m}$$ we define a set of polynomials $f\bullet g =\left\{ F_{0},\ldots,F_{n-m}\right\} $, where \[ F_{j}\left( s\right) =a_{j}b_{0}+a_{j+1}b_{1}s+\ldots+a_{j+m}b_{m}s^{m}, \] for $j=0,\ldots,n-m$, and call it \textit{a generalized Hadamard product of $f$ and $g$}. We give sufficient conditions for the Hurwitz stability of $f\bullet g$. The obtained results show that the famous Garloff--Wagner theorem on the Hurwitz stability of the Hadamard product of polynomials is a special case of a more general fact. We also show that for every polynomial with positive coefficients (even not necessarily stable) one can find a polynomial such that their generalized Hadamard product is stable. Some connections with polynomials admitting the Hadamard factorization are also given. Numerical examples complete and illustrate the considerations.


title: Search for a light charged Higgs boson decaying to a W boson and a CP-odd Higgs boson in final states with e$\mu\mu$ or $\mu\mu\mu$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search for a light charged Higgs boson (H$^+$) decaying to a W boson and a CP-odd Higgs boson (A) in final states with e$\mu\mu$ or $\mu\mu\mu$ is performed using data from pp collisions at $\sqrt{s}=$ 13 TeV, recorded by the CMS detector at the LHC and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. In this search, it is assumed that the H$^+$ boson is produced in decays of top quarks, and the A boson decays to two oppositely charged muons. The presence of signals for H$^+$ boson masses between 100 and 160 GeV and A boson masses between 15 and 75 GeV is investigated. No evidence for the production of the H$^+$ boson is found. Upper limits at 95% confidence level are obtained on the combined branching fraction for the decay chain t $\to$ bH$^+$ $\to$ bW$^+$A $\to$ bW$^+\mu^+\mu^-$, of 1.9 $\times$ 10$^{-6}$ to 8.6 $\times$ 10$^{-6}$, depending on the masses of the H$^+$ and A bosons. These are the first limits for these decay modes of the H$^+$ and A bosons.


title: Braiding properties of worldline configurations in hardcore lattice bosons

abstract: In this manuscript, we study braiding properties of worldline configurations for a variety of ground-states of hardcore Bose-Hubbard models in two dimensions. Configurations are collections of particle paths and result from the path-integral formulation of statistical mechanics. For hard-core bosons, configurations can be seen as geometric braids and therefore can be assigned a certain topological structure, i.e. a way to classify braiding events among worldlines. By means of Monte Carlo calculations, we study superfluid phase and a variety of insulating phases and observe that ground-states of different quantum phases correspond to different probability distributions of braiding properties.


title: Speeding up the Karatsuba algorithm

abstract: This paper describes an $\sim {\cal O}(n)$ pre-compute technique to speed up the Karatsuba algorithm for multiplying two numbers.


title: Optimizing Interim Analysis Timing for Bayesian Adaptive Commensurate Designs

abstract: In developing products for rare diseases, statistical challenges arise due to the limited number of patients available for participation in drug trials and other clinical research. Bayesian adaptive clinical trial designs offer the possibility of increased statistical efficiency, reduced development cost and ethical hazard prevention via their incorporation of evidence from external sources (historical data, expert opinions, and real-world evidence), and flexibility in the specification of interim looks. In this paper, we propose a novel Bayesian adaptive commensurate design that borrows adaptively from historical information and also uses a particular payoff function to optimize the timing of the study's interim analysis. The trial payoff is a function of how many samples can be saved via early stopping and the probability of making correct early decisions for either futility or efficacy. We calibrate our Bayesian algorithm to have acceptable long-run frequentist properties (Type I error and power) via simulation at the design stage. We illustrate our approach using a pediatric trial design setting testing the effect of a new drug for a rare genetic disease. The optimIA R package available at https://github.com/wxwx1993/Bayesian_IA_Timing provides an easy-to-use implementation of our approach.


title: Overfitting in Synthesis: Theory and Practice (Extended Version)

abstract: In syntax-guided synthesis (SyGuS), a synthesizer's goal is to automatically generate a program belonging to a grammar of possible implementations that meets a logical specification. We investigate a common limitation across state-of-the-art SyGuS tools that perform counterexample-guided inductive synthesis (CEGIS). We empirically observe that as the expressiveness of the provided grammar increases, the performance of these tools degrades significantly. We claim that this degradation is not only due to a larger search space, but also due to overfitting. We formally define this phenomenon and prove no-free-lunch theorems for SyGuS, which reveal a fundamental tradeoff between synthesizer performance and grammar expressiveness. A standard approach to mitigate overfitting in machine learning is to run multiple learners with varying expressiveness in parallel. We demonstrate that this insight can immediately benefit existing SyGuS tools. We also propose a novel single-threaded technique called hybrid enumeration that interleaves different grammars and outperforms the winner of the 2018 SyGuS competition (Inv track), solving more problems and achieving a $5\times$ mean speedup.


title: Neural Metric Learning for Fast End-to-End Relation Extraction

abstract: Relation extraction (RE) is an indispensable information extraction task in several disciplines. RE models typically assume that named entity recognition (NER) is already performed in a previous step by another independent model. Several recent efforts, under the theme of end-to-end RE, seek to exploit inter-task correlations by modeling both NER and RE tasks jointly. Earlier work in this area commonly reduces the task to a table-filling problem wherein an additional expensive decoding step involving beam search is applied to obtain globally consistent cell labels. In efforts that do not employ table-filling, global optimization in the form of CRFs with Viterbi decoding for the NER component is still necessary for competitive performance. We introduce a novel neural architecture utilizing the table structure, based on repeated applications of 2D convolutions for pooling local dependency and metric-based features, that improves on the state-of-the-art without the need for global optimization. We validate our model on the ADE and CoNLL04 datasets for end-to-end RE and demonstrate $\approx 1\%$ gain (in F-score) over prior best results with training and testing times that are seven to ten times faster --- the latter highly advantageous for time-sensitive end user applications.


title: Bidirectional Information Flow and the Roles of Privacy Masks in Cloud-Based Control

abstract: We consider a cloud-based control architecture for a linear plant with Gaussian process noise, where the state of the plant contains a client's sensitive information. We assume that the cloud tries to estimate the state while executing a designated control algorithm. The mutual information between the client's actual state and the cloud's estimate is adopted as a measure of privacy loss. We discuss the necessity of uplink and downlink privacy masks. After observing that privacy is not necessarily a monotone function of the noise levels of privacy masks, we discuss the joint design procedure for uplink and downlink privacy masks. Finally, the trade-off between privacy and control performance is explored.


title: Twisted complexes and simplicial homotopies

abstract: In this paper we consider the dg-category of twisted complexes over simplicial ringed spaces. It is clear that a simplicial map $f: (\mathcal{U},\mathcal{R})\to (\mathcal{V}, \mathcal{S})$ between simplicial ringed spaces induces a dg-functor $f^*: \text{Tw}(\mathcal{V}, \mathcal{S})\to \text{Tw}(\mathcal{U}, \mathcal{R})$ where $\text{Tw}(\mathcal{U}, \mathcal{R})$ denotes the dg-category of twisted complexes on $(\mathcal{U},\mathcal{R})$. In this paper we prove that for simplicial homotopic maps $f$ and $g$, there exists an $A_{\infty}$-natural transformation $\Phi: f^*\Rightarrow g^*$ between induced dg-functors. Moreover the $0$th component of $\Phi$ is an objectwise weak equivalence. If we restrict ourselves to the full dg-subcategory of twisted perfect complexes, then we prove that $\Phi$ admits an $A_{\infty}$-quasi-inverse when $(\mathcal{U},\mathcal{R})$ satisfies some additional conditions.


title: Polynomial Time Algorithms for Estimating Spectra of Adiabatic Hamiltonians

abstract: Much research regarding quantum adiabatic optimization has focused on stoquastic Hamiltonians with Hamming symmetric potentials, such as the well studied "spike" example. Due to the large amount of symmetry in these potentials such problems are readily open to analysis both analytically and computationally. However, more realistic potentials do not have such a high degree of symmetry and may have many local minima. Here we present a somewhat more realistic class of problems consisting of many individually Hamming symmetric potential wells. For two or three such wells we demonstrate that such a problem can be solved exactly in time polynomial in the number of qubits and wells. For greater than three wells, we present a tight binding approach with which to efficiently analyze the performance of such Hamiltonians in an adiabatic computation. We provide several basic examples designed to highlight the usefulness of this toy model and to give insight into using the tight binding approach to examining it, including: (1) adiabatic unstructured search with a transverse field driver and a prior guess to the marked item and (2) a scheme for adiabatically simulating the ground states of small collections of strongly interacting spins, with an explicit demonstration for an Ising model Hamiltonian.


title: Fluctuation diagnostics of the finite temperature quasi-antiferromagnetic regime of the 2D Hubbard model

abstract: We study the finite temperature Fermi-liquid to non-Fermi-liquid crossover in the 2D Hubbard model for a range of dopings using the self-consistent ladder dual fermion method. We consider relatively high temperatures where we identify a suppression of the density of states near the Fermi level caused by a quasi-antiferromagnetic behaviour that is itself characterized by a long, but finite, correlation length scale. We perform fluctuation diagnostics to decompose the single-particle self energy into scattering $q$-vector and bosonic frequency contributions. Within this framework we find that the key contributions to the single-particle self energy that give non-Fermi-liquid character, even at weak coupling, are caused by relatively sharp $q=(\pi,\pi)$ spin fluctuations, while the decomposition in the bosonic frequency channel shows a complicated dependence on the relative strengths of zero, positive and negative frequency contributions. Finally, variation in density suggests that the tendency towards non-Fermi-liquid behavior is not substantially different for electron or hole doped systems.


title: The formal affine Demazure algebra and real finite reflection groups

abstract: In this paper, we generalize the formal affine Demazure algebra of Hoffnung-Malag\'on-L\'opez-Savage-Zainoulline to all real finite reflection groups. We begin by generalizing the formal group ring of Calm\`es-Petrov-Zainoulline to all real finite reflection groups. We then define and study the formal Demazure operators that act on the formal group ring. Using these results and constructions, we define and study the formal affine Demazure algebra for all real finite reflection groups. Finally, we compute several structure coefficients that appear in a braid relation among the formal Demazure elements, and we conclude this paper by computing all structure coefficients for the reflection groups $I_2(5)$, $I_2(7)$, $H_3$, and $H_4$.


title: A Multi-Task Learning Framework for Extracting Drugs and Their Interactions from Drug Labels

abstract: Preventable adverse drug reactions as a result of medical errors present a growing concern in modern medicine. As drug-drug interactions (DDIs) may cause adverse reactions, being able to extracting DDIs from drug labels into machine-readable form is an important effort in effectively deploying drug safety information. The DDI track of TAC 2018 introduces two large hand-annotated test sets for the task of extracting DDIs from structured product labels with linkage to standard terminologies. Herein, we describe our approach to tackling tasks one and two of the DDI track, which corresponds to named entity recognition (NER) and sentence-level relation extraction respectively. Namely, our approach resembles a multi-task learning framework designed to jointly model various sub-tasks including NER and interaction type and outcome prediction. On NER, our system ranked second (among eight teams) at 33.00% and 38.25% F1 on Test Sets 1 and 2 respectively. On relation extraction, our system ranked second (among four teams) at 21.59% and 23.55% on Test Sets 1 and 2 respectively.


title: Optimizing Sequential Medical Treatments with Auto-Encoding Heuristic Search in POMDPs

abstract: Health-related data is noisy and stochastic in implying the true physiological states of patients, limiting information contained in single-moment observations for sequential clinical decision making. We model patient-clinician interactions as partially observable Markov decision processes (POMDPs) and optimize sequential treatment based on belief states inferred from history sequence. To facilitate inference, we build a variational generative model and boost state representation with a recurrent neural network (RNN), incorporating an auxiliary loss from sequence auto-encoding. Meanwhile, we optimize a continuous policy of drug levels with an actor-critic method where policy gradients are obtained from a stablized off-policy estimate of advantage function, with the value of belief state backed up by parallel best-first suffix trees. We exploit our methodology in optimizing dosages of vasopressor and intravenous fluid for sepsis patients using a retrospective intensive care dataset and evaluate the learned policy with off-policy policy evaluation (OPPE). The results demonstrate that modelling as POMDPs yields better performance than MDPs, and that incorporating heuristic search improves sample efficiency.