title: Forecasting Future Murders of Mr. Boddy by Numerical Weather Prediction

abstract: Despite a previous description of his state as a stable fixed point, just past midnight this morning Mr. Boddy was murdered again. In fact, over 70 years Mr. Boddy has been reported murdered $10^6$ times, while there exist no documented attempts at intervention. Using variational data assimilation, we train a model of Mr. Boddy's dynamics on the time series of observed murders, to forecast future murders. The parameters to be estimated include instrument, location, and murderer. We find that a successful estimation requires three additional elements. First, to minimize the effects of selection bias, generous ranges are placed on parameter searches, permitting values such as the Cliff, the Poisoned Apple, and the Wife. Second, motive, which was not considered relevant to previous murders, is added as a parameter. Third, Mr. Boddy's little-known asthmatic condition is considered as an alternative cause of death. Following this morning's event, the next local murder is forecast for 17:19:03 EDT this afternoon, with a standard deviation of seven hours, at The Kitchen at 4330 Katonah Avenue, Bronx, NY, 10470, with either the Lead Pipe or the Lead Bust of Washington Irving. The motive is: Case of Mistaken Identity, and there was no convergence upon a murderer. Testing of the procedure's predictive power will involve catching the D train to 205th Street and a few transfers over to Katonah Avenue, and sitting around waiting with our eyes peeled. We discuss the problem of identifying a global solution - that is, the best reason for murder on a landscape riddled with pretty-decent reasons. We also discuss the procedure's assumption of Gaussian-distributed errors, which will under-predict rare events. This under-representation of highly improbable events may be offset by the fact that the training data, after all, consists of multiple murders of a single person.


title: Stable, Concurrent Controller Composition for Multi-Objective Robotic Tasks

abstract: Robotic systems often need to consider multiple tasks concurrently. This challenge calls for controller synthesis algorithms that fulfill multiple control specifications while maintaining the stability of the overall system. In this paper, we decompose multi-objective tasks into subtasks, where individual subtask controllers are designed independently and then combined to generate the overall control policy. In particular, we adopt Riemannian Motion Policies (RMPs), a recently proposed controller structure in robotics, and, RMPflow, its associated computational framework for combining RMP controllers. We re-establish and extend the stability results of RMPflow through a rigorous Control Lyapunov Function (CLF) treatment. We then show that RMPflow can stably combine individually designed subtask controllers that satisfy certain CLF constraints. This new insight leads to an efficient CLF-based computational framework to generate stable controllers that consider all the subtasks simultaneously. Compared with the original usage of RMPflow, our framework provides users the flexibility to incorporate design heuristics through nominal controllers for the subtasks. We validate the proposed computational framework through numerical simulation and robotic implementation.


title: The ALMA-PILS survey: Gas dynamics in IRAS 16293$-$2422 and the connection between its two protostars

abstract: [Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We selected molecular gas line transitions of CO, H2CO, HCN, CS, SiO, and CCH from the ALMA-PILS spectral imaging survey (329-363 GHz) and used them as tracers of kinematics, density, and temperature in the IRAS16293-2422 system. The angular resolution of the PILS data set allows us to study these quantities at a resolution of 0.5 arcsec (60 au [..]). Line-of-sight velocity maps of both optically thick and optically thin molecular lines reveal: (i) new manifestations of previously known outflows emanating from protostar A; (ii) a kinematically quiescent bridge of dust and gas spanning between the two protostars, with an inferred density between 4 10^4 and 3 10^7 cm^-3; and (iii) a separate, straight filament seemingly connected to protostar B seen only in CCH, with a flat kinematic signature. Signs of various outflows, all emanating from source A, are evidence of high-density and warmer gas; none of them coincide spatially and kinematically with the bridge. We hypothesize that the bridge arc is a remnant of filamentary substructure in the protostellar envelope material from which protostellar sources A and B have formed. One particular morphological structure appears to be due to outflowing gas impacting the quiescent bridge material. The continuing lack of clear outflow signatures unambiguously associated to protostar B and the vertically extended shape derived for its disk-like structure lead us to conclude that source B may be in an earlier evolutionary stage than source A.


title: The Power Light Cone of the Discrete Bak-Sneppen, Contact and other local processes

abstract: We consider a class of random processes on graphs that include the discrete Bak-Sneppen (DBS) process and the several versions of the contact process (CP), with a focus on the former. These processes are parametrized by a probability $0\leq p \leq 1$ that controls a local update rule. Numerical simulations reveal a phase transition when $p$ goes from 0 to 1. Analytically little is known about the phase transition threshold, even for one-dimensional chains. In this article we consider a power-series approach based on representing certain quantities, such as the survival probability or the expected number of steps per site to reach the steady state, as a power-series in $p$. We prove that the coefficients of those power series stabilize as the length $n$ of the chain grows. This is a phenomenon that has been used in the physics community but was not yet proven. We show that for local events $A,B$ of which the support is a distance $d$ apart we have $\mathrm{cor}(A,B) = \mathcal{O}(p^d)$. The stabilization allows for the (exact) computation of coefficients for arbitrary large systems which can then be analyzed using the wide range of existing methods of power series analysis.


title: How the power spectrum of dust continuum images may hide the presence of a characteristic filament width

abstract: Herschel observations of interstellar clouds support a paradigm for star formation in which filaments play a central role. One of the foundations of this paradigm is the finding, based on studies of the radial column density profiles, that nearby filaments share a common inner width of ~0.1 pc. The existence of a characteristic filament width has been questioned, however, as it seems inconsistent with the scale-free nature of the power spectrum of cloud images. In an effort to clarify the origin of this apparent discrepancy, we examined the power spectra of the Herschel 250 micron images of the Polaris, Aquila, and Taurus clouds and performed a number of simple numerical experiments by injecting synthetic filaments in both the Herschel images and synthetic background images. We constructed several populations of synthetic filaments of 0.1 pc width with realistic area filling factors ($A_{\rm fil}$) and distributions of column density contrasts ($\delta_c$). After adding synthetic filaments to the Herschel images, we re-computed the image power spectra and compared the results with the original, essentially scale-free power spectra. We used the $\chi^2$ variance of the residuals between the best power-law fit and the output power spectrum as a diagnostic of the presence of a significant departure from a scale-free power spectrum. We found that the $\chi^2$ variance depends primarily on the combined parameter $\delta_c^2\, A_{\rm fil}$. Analysis of the real Herschel data shows that the observed filamentary structure contributes only 1/5 of the power in the power spectra at angular frequencies where an effect of the characteristic filament width is expected. We conclude that the scale-free power spectra of Herschel images remain consistent with the existence of a characteristic filament width ~0.1 pc and do not invalidate the conclusions drawn from studies of the filament profiles.


title: Methods for suspensions of passive and active filaments

abstract: Flexible filaments and fibres are essential components of important complex fluids that appear in many biological and industrial settings. Direct simulations of these systems that capture the motion and deformation of many immersed filaments in suspension remain a formidable computational challenge due to the complex, coupled fluid--structure interactions of all filaments, the numerical stiffness associated with filament bending, and the various constraints that must be maintained as the filaments deform. In this paper, we address these challenges by describing filament kinematics using quaternions to resolve both bending and twisting, applying implicit time-integration to alleviate numerical stiffness, and using quasi-Newton methods to obtain solutions to the resulting system of nonlinear equations. In particular, we employ geometric time integration to ensure that the quaternions remain unit as the filaments move. We also show that our framework can be used with a variety of models and methods, including matrix-free fast methods, that resolve low Reynolds number hydrodynamic interactions. We provide a series of tests and example simulations to demonstrate the performance and possible applications of our method. Finally, we provide a link to a MATLAB/Octave implementation of our framework that can be used to learn more about our approach and as a tool for filament simulation.


title: Knife-edge based measurement of the 4D transverse phase space of electron beams with picometer-scale emittance

abstract: Precise manipulation of high brightness electron beams requires detailed knowledge of the particle phase space shape and evolution. As ultrafast electron pulses become brighter, new operational regimes become accessible with emittance values in the picometer range, with enormous impact on potential scientific applications. Here we present a new characterization method for such beams and demonstrate experimentally its ability to reconstruct the 4D transverse beam matrix of strongly correlated electron beams with sub-nanometer emittance and sub-micrometer spot size, produced with the HiRES beamline at LBNL. Our work extends the reach of ultrafast electron accelerator diagnostics into picometer-range emittance values, opening the way to complex nanometer-scale electron beam manipulation techniques.


title: Information content of queries in training Parameterized Quantum Circuits

abstract: Parameterized quantum circuits (PQC, aka, variational quantum circuits) are among the proposals for a computational advantage over classical computation of near-term (not fault tolerant) digital quantum computers. PQCs have to be "trained" -- i.e., the expectation value function has to be maximized over the space of parameters. This paper deals with the number of samples (or "runs" of the quantum computer) which are required to train the PQC, and approaches it from an information theoretic viewpoint. The main take-away is a disparity in the large amount of information contained in a single exact evaluation of the expectation value, vs the exponentially small amount contained in the random sample obtained from a single run of the quantum circuit.


title: Topology of the Stokes phenomenon

abstract: Several intrinsic topological ways to encode connections on vector bundles on smooth complex algebraic curves will be described. In particular the notion of {\em Stokes decompositions} will be formalised, as a convenient intermediate category between the Stokes filtrations and the Stokes local systems/wild monodromy representations. The main result establishes a new simple characterisation of the Stokes decompositions.


title: Estimating Total Open Heliospheric Magnetic Flux

abstract: Over the solar-activity cycle, there are extended periods where significant discrepancies occur between the spacecraft-observed total (unsigned) open magnetic flux and that determined from coronal models. In this article, the total open heliospheric magnetic flux is computed using two different methods and then compared with results obtained from in-situ interplanetary magnetic-field observations. The first method uses two different types of photospheric magnetic-field maps as input to the Wang Sheeley Arge (WSA) model: i) traditional Carrington or diachronic maps, and ii) Air Force Data Assimilative Photospheric Flux Transport model synchronic maps. The second method uses observationally derived helium and extreme-ultraviolet coronal-hole maps overlaid on the same magnetic-field maps in order to compute total open magnetic flux. The diachronic and synchronic maps are both constructed using magnetograms from the same source, namely the National Solar Observatory Kitt Peak Vacuum Telescope and Vector Spectromagnetograph. The results of this work show that the total open flux obtained from observationally derived coronal holes agrees remarkably well with that derived from WSA, especially near solar minimum. This suggests that, on average, coronal models capture well the observed large-scale coronal-hole structure over most of the solar cycle. Both methods show considerable deviations from total open flux deduced from spacecraft data, especially near solar maximum, pointing to something other than poorly determined coronal-hole area specification as the source of these discrepancies.


title: SOAP: A generalised application of the Viterbi algorithm to searches for continuous gravitational-wave signals

abstract: All-sky and wide parameter space searches for continuous gravitational waves are generally template-matching schemes which test a bank of signal waveforms against data from a gravitational wave detector. Such searches can offer optimal sensitivity for a given computing cost and signal model, but are highly-tuned to specific signal types and are computationally expensive, even for semi-coherent searches. We have developed a search method based on the well-known Viterbi algorithm which is model-agnostic and has a computational cost several orders of magnitude lower than template methods, with a modest reduction in sensitivity. In particular, this method can search for signals which have an unknown frequency evolution. We test the algorithm on three simulated and real data sets: gapless Gaussian noise, Gaussian noise with gaps and real data from the final run of initial LIGO (S6). We show that at 95% efficiency, with a 1% false alarm rate, the algorithm has a depth sensitivity of $\sim 33$, $10$ and $13$ ,Hz$^{-1/2}$ with corresponding SNRs of $\sim 60$, $72$ and $74$ in these datasets. we discuss the use of this algorithm for detecting a wide range of quasi-monochromatic gravitational wave signals and instrumental lines.


title: Encoding an oscillator into many oscillators

abstract: An outstanding challenge for quantum information processing using bosonic systems is Gaussian errors such as excitation loss and added thermal noise errors. Thus, bosonic quantum error correction (QEC) is essential. Most bosonic QEC schemes encode a finite-dimensional logical qubit or qudit into noisy bosonic oscillator modes. In this case, however, the infinite-dimensional bosonic nature of the physical system is lost at the error-corrected logical level. On the other hand, there are several proposals for encoding an oscillator mode into many noisy oscillator modes. However, these oscillator-into-oscillators encoding schemes are in the class of Gaussian quantum error correction. Therefore, these codes cannot correct practically relevant Gaussian errors due to the established no-go theorems which state that Gaussian errors cannot be corrected by using only Gaussian resources. Here, we circumvent these no-go results and show that it is possible to correct Gaussian errors by using Gottesman-Kitaev-Preskill (GKP) states as non-Gaussian resources. In particular, we propose a non-Gaussian oscillator-into-oscillators code, the two-mode GKP-repetition code, and demonstrate that it can correct additive Gaussian noise errors. In addition, we generalize the two-mode GKP-repetition code to an even broader class of non-Gaussian oscillator codes, namely, GKP-stabilizer codes. Specifically, we show that there exists a highly hardware-efficient GKP-stabilizer code, the GKP-two-mode-squeezing code, that can quadratically suppress additive Gaussian noise errors in both the position and momentum quadratures up to a small logarithmic correction. Moreover, for any GKP-stabilizer code, we show that logical Gaussian operations can be readily implemented by using only physical Gaussian operations. We also show that our oscillator codes can correct practically relevant excitation loss and thermal noise errors.


title: Activity Classification Using Smartphone Gyroscope and Accelerometer Data

abstract: Activities, such as walking and sitting, are commonly used in biomedical settings either as an outcome or covariate of interest. Researchers have traditionally relied on surveys to quantify activity levels of subjects in both research and clinical settings, but surveys are not objective in nature and have many known limitations, such as recall bias. Smartphones provide an opportunity for unobtrusive objective measurement of various activities in naturalistic settings, but their data tends to be noisy and needs to be analyzed with care. We explored the potential of smartphone accelerometer and gyroscope data to distinguish between five different types of activity: walking, sitting, standing, ascending stairs, and descending stairs. We conducted a study in which four participants followed a study protocol and performed a sequence of various activities with one phone in their front pocket and another phone in their back pocket. The subjects were filmed throughout, and the obtained footage was annotated to establish ground truth activity. We applied the so-called movelet method to classify their activity. Our results demonstrate the promise of smartphones for activity detection in naturalistic settings, but they also highlight common challenges in this field of research.


title: Some Experimental Results of Relieving Discomfort in Virtual Reality by Disturbing Feedback Loop in Human Brain

abstract: Recently, great progress has been made in virtual reality(VR) research and application. However, virtual reality faces a big problem since its appearance, i.e. discomfort (nausea, stomach awareness, etc). Discomfort can be relieved by increasing hardware (sensor, cpu and display) speed. But this will increase cost. This paper gives another low cost solution. The phenomenon of cybersickness is explained with the control theory: discomfort arises if feedback scene differs from expectation, so it can be relieved by disturbing feedback loop in human brain. A hardware platform is build to test this explanation. The VR display on a Samsung S6 is blurred while head movement is detected. The effect is evaluated by comparing responses to the Simulated Sickness Questionnaire (SSQ) between a control and experimental condition. Experimental results show that the new method can ease discomfort remarkably with little extra cost. As a result, VR may be used more widely in teaching (like foreign language, medicine). It's also reasonable to expect likewise merits in other VR applications.


title: Optimal power and efficiency of single quantum dot heat engines: theory and experiment

abstract: Quantum dots (QDs) can serve as near perfect energy filters and are therefore of significant interest for the study of thermoelectric energy conversion close to thermodynamic efficiency limits. Indeed, recent experiments in [Nat. Nano. 13, 920 (2018)] realized a QD heat engine with performance near these limits and in excellent agreement with theoretical predictions. However, these experiments also highlighted a need for more theory to help guide and understand the practical optimization of QD heat engines, in particular regarding the role of tunnel couplings on the performance at maximum power and efficiency for QDs that couple seemingly weakly to electronic reservoirs. Furthermore, these experiments also highlighted the critical role of the external load when optimizing the performance of a QD heat engine in practice. To provide further insight into the operation of these engines we use the Anderson impurity model together with a Master equation approach to perform power and efficiency calculations up to co-tunneling order. This is combined with additional thermoelectric experiments on a QD embedded in a nanowire where the power is measured using two methods. We use the measurements to present an experimental procedure for efficiently finding the external load $R_P$ which should be connected to the engine to optimize power output. Our theoretical estimates of $R_P$ show a good agreement with the experimental results, and we show that second order tunneling processes and non-linear effects have little impact close to maximum power, allowing us to derive a simple analytic expression for $R_P$. In contrast, we find that the electron contribution to the thermoelectric efficiency is significantly reduced by second order tunneling processes, even for rather weak tunnel couplings.


title: A diameter gap for quotients of the unit sphere

abstract: We prove that for any isometric action of a group on a unit sphere of dimension larger than one, the quotient space has diameter zero or larger than a universal dimension-independent positive constant.


title: Alternating Weak Automata from Universal Trees

abstract: An improved translation from alternating parity automata on infinite words to alternating weak automata is given. The blow-up of the number of states is related to the size of the smallest universal ordered trees and hence it is quasi-polynomial, and only polynomial if the asymptotic number of priorities is logarithmic in the number of states. This is an exponential improvement on the translation of Kupferman and Vardi (2001) and a quasi-polynomial improvement on the translation of Boker and Lehtinen (2018). Any slightly better such translation would (if---like all presently known such translations---it is efficiently constructive) lead to algorithms for solving parity games that are asymptotically faster in the worst case than the current state of the art (Calude, Jain, Khoussainov, Li, and Stephan, 2017; Jurdzi\'nski and Lazi\'c, 2017; and Fearnley, Jain, Schewe, Stephan, and Wojtczak, 2017), and hence it would yield a significant breakthrough.


title: The Complexity Geometry of a Single Qubit

abstract: The computational complexity of a quantum state quantifies how hard it is to make. `Complexity geometry', first proposed by Nielsen, is an approach to defining computational complexity using the tools of differential geometry. Here we demonstrate many of the attractive features of complexity geometry using the example of a single qubit, which turns out to be rich enough to be illustrative but simple enough to be illuminating.


title: Asymptotic behaviour of the one-dimensional "rock-paper-scissors" cyclic cellular automaton

abstract: The one-dimensional three-state cyclic cellular automaton is a simple spatial model with three states in a cyclic "rock-paper-scissors" prey-predator relationship. Starting from a random configuration, similar states gather in increasingly large clusters; asymptotically, any finite region is filled with a uniform state that is, after some time, driven out by its predator, each state taking its turn in dominating the region (heteroclinic cycles). We consider the situation where each site in the initial configuration is chosen independently at random with a different probability for each state. We prove that the asymptotic probability that a state dominates a finite region corresponds to the initial probability of its prey. The proof methods are based on discrete probability tools, mainly particle systems and random walks.


title: Update on decaying and annihilating heavy dark matter with the 6-year IceCube HESE data

abstract: In view of the IceCube's 6-year high-energy starting events (HESE) sample, we revisit the possibility that the updated data may be better explained by a combination of neutrino fluxes from dark matter decay and an isotropic astrophysical power-law than purely by the latter. We find that the combined two-component flux qualitatively improves the fit to the observed data over a purely astrophysical one, and discuss how these updated fits compare against a similar analysis done with the 4-year HESE data. We also update fits involving dark matter decay via multiple channels, without any contribution from the astrophysical flux. We find that a DM-only explanation is not excluded by neutrino data alone. Finally, we also consider the possibility of a signal from dark matter annihilations and perform analogous analyses to the case of decays, commenting on its implications.


title: Numerical study of the Lorentzian Engle-Pereira-Rovelli-Livine spin foam amplitude

abstract: The Lorentzian EPRL spin foam amplitude for loop quantum gravity is a multi-dimensional non-compact integral of highly oscillating functions. Using a method based on the decomposition of Clebsch-Gordan coefficients for the unitary infinite-dimensional representations of SL(2,C) in terms of those of SU(2), we are able to provide for the first time numerical evaluations of the vertex amplitude. The values obtained support the asymptotic formula obtained by Barrett and collaborators with a saddle point approximation, showing, in particular, a power-law decay and oscillations related to the Regge action. The comparison offers a test of the efficiency of the method. Truncating the decomposition to the first few terms provides a qualitative matching of the power-law decay and oscillations. For vector and Euclidean Regge boundary data, a qualitative matching is obtained with just the first term, which corresponds to the simplified EPRL model. We comment on future developments for the numerics and extension to higher vertices. We complete our work with some analytic results: We provide an algorithm and explicit configurations for the different geometries that can arise as boundary data, and explain the geometric consequences of the decomposition used.


title: Nematic and Antiferromagnetic Quantum Criticality in a Multi-Orbital Hubbard Model for Iron Pnictides

abstract: The extent to which quantum criticality drives the physics of iron pnictides is a central question in the field. Earlier theoretical considerations were based on an effective field theory, and the proposed realization in P-doped iron arsenides has received extensive experimental evidence. To connect the quantum critical behavior with the underlying electronic physics, it is important to analyze it within microscopic models. Here, we do so for a multi-orbital model containing both Hubbard and Hund's interactions, by a variational Monte Carlo method based on Jastrow-Slater wave functions that allow for a non-perturbative treatment of electron correlations. We find strong evidence for the existence of a unique quantum critical point, where both nematic and $(\pi,0)$ antiferromagnetic orders develop together, in the bad-metal regime of the phase diagram. Implications of our results for the iron-based superconductivity are discussed.


title: Integrating Semantic Knowledge to Tackle Zero-shot Text Classification

abstract: Insufficient or even unavailable training data of emerging classes is a big challenge of many classification tasks, including text classification. Recognising text documents of classes that have never been seen in the learning stage, so-called zero-shot text classification, is therefore difficult and only limited previous works tackled this problem. In this paper, we propose a two-phase framework together with data augmentation and feature augmentation to solve this problem. Four kinds of semantic knowledge (word embeddings, class descriptions, class hierarchy, and a general knowledge graph) are incorporated into the proposed framework to deal with instances of unseen classes effectively. Experimental results show that each and the combination of the two phases achieve the best overall accuracy compared with baselines and recent approaches in classifying real-world texts under the zero-shot scenario.


title: $W$ radiative decay to heavy-light mesons in HQET factorization through ${\cal O}(\alpha_s)$

abstract: Analogous to NRQCD factorization for heavy quarkonium exclusive production, in this work we propose to employ the heavy-quark-effective-theory (HQET) factorization, which has been predominantly applied to account for exclusive $B$ decays, to study the exclusive production of the heavy-flavored mesons. We take $W\to B(D_s)+\gamma$ as a prototype process. The validity of the HQET factorization rests upon the presumed scale hierarchy: $m_W\sim m_b\gg \Lambda_{\rm QCD}$. Through an explicit analysis at next-to-leading order in $\alpha_s$ yet at leading order in $1/m_b$, we verify that the decay form factors can indeed be expressed as the convolution between perturbatively calculable hard-scattering kernel and the $B$ meson light-cone distribution amplitude (LCDA) defined in HQET. It is observed that the factorization scale dependence becomes reduced after incorporating the NLO perturbative correction. An interesting future investigation is to identify and resum large collinear logarithms of $m_W/m_b$ that arise ubiquitously in the fixed-order expressions of the hard-scattering kernel in HQET factorization.


title: Magnetorotational Mechanism of the Explosion of Core-Collapse Supernovae

abstract: The idea of the magnetorotational explosion mechanism is that the energy of rotation of the neutron star formed in the course of a collapse is transformed into the energy of an expanding shock wave by means of a magnetic field. In the two-dimensional case, the time of this transformation depends weakly on the initial strength of the poloidal magnetic field because of the development of a magnetorotational instability. Differential rotation leads to the twisting and growth of the toroidal magnetic-field component, which becomes much stronger than the poloidal component. In the case where the initial configuration of the magnetic field is close to a dipole configuration, the ejection of matter has a jet character, whereas, in the case of a quadrupole configuration, there arises an equatorial ejection. In either case, the energy release is sufficient for explaining the observed average energy of supernova explosion. Neutrinos are emitted as the collapse and the formation of a rapidly rotating neutron star proceeds. In addition, neutrino radiation arises in the process of magnetorotational explosion owing to additional rotational-energy losses. In order to explain an interval of 4.5 hours between the two observed neutrino signals from SN 1987A, it is necessary to assume a weakening of the magnetorotional instability and a small initial magnetic field (10(9)-10(10)G) in the newly formed rotating neutron star. The existence of a black hole in the SN 1987A remnant could explain the absence of any visible pointlike source at the center of the explosion.


title: A resource theory of Maxwell's demons

abstract: Motivated by recent progress on the motive power of information in quantum thermodynamics, we put forth an operational resource theory of Maxwell's demons. We show that the resourceful ({\em daemonic}) states can be partitioned into at most nine irreducible subsets. The sets can be classified by a rank akin to the Schmidt rank for entanglement theory. Moreover, we show that there exists a natural monotone, called the wickedness, which quantifies the multilevel resource content of the states. The present resource theory is shown to share deep connections with the resource theory of thermodynamics. In particular, the nine irreducible sets are found to be characterized by well defined temperatures which, however, are not monotonic in the wickedness. This result, as we demonstrate, is found to have dramatic consequences for Landauer's erasure principle. Our analysis therefore settles a longstanding debate surrounding the identity of Maxwell's demons and the operational significance of other related fundamental thermodynamic entities.


title: Quantum differential ghost microscopy

abstract: Quantum correlations become formidable tools for beating classical capacities of measurement. Preserving these advantages in practical systems, where experimental imperfections are unavoidable, is a challenge of the utmost importance. Here we propose and realize a quantum ghost imaging protocol able to compensate for the detrimental effect of detection noise and losses. This represents an important improvement as quantum correlations allow low brightness imaging, desirable for reducing the absorption dose. In particular, we develop a comprehensive model starting from a ghost imaging scheme elaborated for bright thermal light, known as differential ghost imaging and particularly suitable in the relevant case of faint or sparse objects. We perform the experiment using SPDC light in microscopic configuration. The image is reconstructed exploiting non-classical intensity correlation rather than photon pairs detection coincidences. On one side we validate the theoretical model and on the other we show the applicability of this technique by reconstructing a biological object with 5 micrometers resolution.


title: Differential Observation Techniques for the SZE-21cm and radio sources

abstract: The SZE-21cm has been proposed as an alternative probe for the Dark Ages (DA) and the Epoch of Reionization (EoR). The effect is produced when photons of the 21cm background are inverse Compton (up-)scattered by electrons residing in hot plasma of cosmic structures and can be studied through differential observations of the large-scale structure, towards and away a region of the radio sky containing the cosmic structure of interest. This work makes use of the 21cmFAST code to simulate low frequency radio observations of a galaxy cluster from which we then extract the SZE-21cm signal. We further explore uses of the differential technique to radio observations of active radio galaxies using radio observations from the MWA GLEAM survey. The differential techniques this work studies enable us to extract the SZE-21cm signal from simulated galaxy cluster observations, which can then be used to obtain the global features of the 21cm signal during the DA and EoR. The technique shows further benefits in source extraction and characterisation, and feature enhancement in radio observations particularly for low-surface brightness and extended radio sources.


title: Critical excitation-inhibition balance in dense neural networks

abstract: The "edge of chaos" phase transition in artificial neural networks is of renewed interest in light of recent evidence for criticality in brain dynamics. Statistical mechanics traditionally studied this transition with connectivity $k$ as the control parameter and an exactly balanced excitation-inhibition ratio. While critical connectivity has been found to be low in these model systems, typically around $k=2$, which is unrealistic for natural systems, a recent study utilizing the excitation-inhibition ratio as the control parameter found a new, nearly degree independent, critical point when connectivity is large. However, the new phase transition is accompanied by an unnaturally high level of activity in the network. Here we study random neural networks with the additional properties of (i) a high clustering coefficient and (ii) neurons that are solely either excitatory or inhibitory, a prominent property of natural neurons. As a result, we observe an additional critical point for networks with large connectivity, regardless of degree distribution, which exhibits low activity levels that compare well with neuronal brain networks.


title: The representation of asteroid shapes: a test for the inversion of Gaia photometry

abstract: It is common practice nowadays to derive spins and 3D shapes of asteroids from the inversion of photometric light curves. However, this method requires a good number of photometric points and dedicated observing sessions. On the other hand, the photometric observations carried out by the Gaia mission will be sparse and their number relatively small. For this reason, a multi-parametric shape described by a large number of elementary facets cannot probably be derived from Gaia data alone. Therefore, the Data Processing and Analysis Consortium (DPAC), implemented a simpler solution as an unattended data reduction pipeline which relies on three axial ellipsoids for the shape representation. However, overall accuracy of such triaxial shape solutions has to be assessed. How adequate is an ellipsoidal approximation to represent the overall properties of an irregular body? Which error is made on the volume in comparison a more accurate model with irregular shape? To answer these questions, we have implemented numerical procedures for comparing ellipsoids to more complex and irregular shapes, and we performed a full simulation of the photometric signal from these objects, using both shape representations. Implementing the same software algorithm that will be used for the analysis of Gaia asteroid photometry, rotation period, spin axis orientation and ellipsoidal shape were derived from simulated observations of selected Main Belt asteroids assuming a geometric scattering model (work is in progress for more complex scattering models). We found that the ellipsoids by photometry inversion are closely similar to the best-fitting ellipsoids of the simulated complex shapes and that the error on the asteroid volume (relative to a complex shape) is generally low, usually around $10\%$.


title: Birkhoff polytopes of different type and the orthant-lattice property

abstract: The Birkhoff polytope, defined to be the convex hull of $n\times n$ permutation matrices, is a well studied polytope in the context of the Ehrhart theory. This polytope is known to have many desirable properties, such as the Gorenstein property and existence of regular, unimodular triangulations. In this paper, we study analogues of the Birkhoff polytope for finite irreducible Coxeter groups of other types. We focus on a type-$B$ Birkhoff polytope $\mathcal{BB}(n)$ arising from signed permutation matrices and prove that it and its dual polytope are reflexive, and hence Gorenstein, and also possess regular, unimodular triangulations. Noting that our triangulation proofs do not rely on the combinatorial structure of $\mathcal{BB}(n)$, we define the notion of an orthant-lattice property polytope and use this to prove more general results for the existence of regular, unimodular triangulations and unimodular covers for a significant family of reflexive polytopes. We conclude by remarking on some connections to Gale-duality, Birkhoff polytopes of other types, and possible applications of orthant-lattice property.


title: Monopole production via photon fusion at the LHC

abstract: The existence of magnetic monopoles, also predicted in some GUT theories, would symmetrise Maxwell equations and explain the charge quantisation. Searches for them are being performed in cosmic telescopes as well as in collider experiments, such as MoEDAL and ATLAS. In this report, we focus on the, least explored, photon-fusion mechanism, yet Drell-Yan results are discussed, too. Cross sections for monopoles of spin 0, 1/2 and 1 for an effective monopole-velocity-dependent magnetic charge are presented. For spin-1/2 and spin-1 monopoles, a magnetic-moment term is included, which is treated as a new phenomenological parameter and, together with the velocity-dependent coupling, allows for a perturbative treatment of the cross-section calculation. We present an appropriate implementation of photon-fusion and Drell-Yan processes into MadGraph UFO models, aimed to serve as a useful tool in monopole searches at LHC, especially for photon fusion, given that it has not been considered by experimental collaborations recently. Moreover, the experimental implications of such perturbatively reliable monopole searches are discussed.


title: Recovery of zeroth order coefficients in non-linear wave equations

abstract: This paper is concerned with the resolution of an inverse problem related to the recovery of a scalar (potential) function $V$ from the source to solution map, of the semi-linear equation $(\Box_{g}+V)u+u^3=0$ on a globally hyperbolic Lorentzian manifold $(M,g)$. We first study the simpler model problem where the geometry is the Minkowski space and prove the uniqueness of $V$ through the use of geometric optics and a three-fold wave interaction arising from the cubic non-linearity. Subsequently, the result is generalized to globally hyperbolic Lorentzian manifolds by using Gaussian beams.


title: Temporal correlation between the optical and {\gamma}-ray flux variations in the blazar 3C 454.3

abstract: Blazars show optical and $\gamma$-ray flux variations that are generally correlated, although there are exceptions. Here we present anomalous behaviour seen in the blazar 3C 454.3 based on an analysis of quasi-simultaneous data at optical, UV, X-ray and $\gamma$-ray energies, spanning about 9 years from August 2008 to February 2017. We have identified four time intervals (epochs), A, B, D and E, when the source showed large-amplitude optical flares. In epochs A and B the optical and $\gamma$-ray flares are correlated, while in D and E corresponding flares in $\gamma$-rays are weak or absent. In epoch B the degree of optical polarization strongly correlates with changes in optical flux during a short-duration optical flare superimposed on one of long duration. In epoch E the optical flux and degree of polarization are anti-correlated during both the rising and declining phases of the optical flare. We carried out broad-band spectral energy distribution (SED) modeling of the source for the flaring epochs A,B, D and E, and a quiescent epoch, C. Our SED modeling indicates that optical flares with absent or weak corresponding $\gamma$-ray flares in epochs D and E could arise from changes in a combination of parameters, such as the bulk Lorentz factor, magnetic field and electron energy density, or be due to changes in the location of the $\gamma$-ray emitting regions.


title: Cosmic-ray propagation around the Sun: investigating the influence of the solar magnetic field on the cosmic-ray Sun shadow

abstract: The cosmic-ray Sun shadow, which is caused by high-energy charged cosmic rays being blocked and deflected by the Sun and its magnetic field, has been observed by various experiments, such as Argo-YBJ, HAWC, Tibet, and IceCube. Most notably, the shadow's size and depth was recently shown to correlate with the 11-year solar cycle. The interpretation of such measurements, which help to bridge the gap between solar physics and high-energy particle astrophysics, requires a solid theoretical understanding of cosmic-ray propagation in the coronal magnetic field. It is the aim of this paper to establish theoretical predictions for the cosmic-ray Sun shadow in order to identify observables that can be used to study this link in more detail. To determine the cosmic-ray Sun shadow, we numerically compute trajectories of charged cosmic rays in the energy range of 5-316 TeV for five different mass numbers. We present and analyze the resulting shadow images for protons and iron, as well as for typically measured cosmic-ray compositions. We confirm the observationally established correlation between the magnitude of the shadowing effect and both the mean sunspot number and the polarity of the magnetic field during the solar cycle. We also show that during low solar activity, the Sun's shadow behaves similarly to that of a dipole, for which we find a non-monotonous dependence on energy. In particular, the shadow can become significantly more pronounced than the geometrical disk expected for a totally unmagnetized Sun. For times of high solar activity, we instead predict the shadow to depend monotonously on energy, and to be generally weaker than the geometrical shadow for all tested energies. These effects should become visible in energy-resolved measurements of the Sun shadow, and may in the future become an independent measure for the level of disorder in the solar magnetic field.


title: Optical Supersymmetry in the Time Domain

abstract: Originally emerged within the context of string and quantum field theory, and later fruitfully extrapolated to photonics, the algebraic transformations of quantum-mechanical supersymmetry were conceived in the space realm. Here, we introduce a paradigm shift, demonstrating that Maxwell's equations also possess an underlying supersymmetry in the time domain. As a result, we obtain a simple analytic relation between the scattering coefficients of a large variety of time-varying optical systems and uncover a wide new class of reflectionless, three dimensional, all-dielectric, isotropic, omnidirectional, polarization-independent, non-complex media. Temporal supersymmetry is also shown to arise in dispersive media supporting temporal bound states, which allows engineering their momentum spectra and dispersive properties. These unprecedented features define a promising design platform for free-space and integrated photonics, enabling the creation of a number of novel reconfigurable reflectionless devices, such as frequency-selective, polarization-independent and omnidirectional invisible materials, compact frequency-independent phase shifters, broadband isolators, and versatile pulse-shape transformers.


title: A new metric for statistical properties of long time behaviors

abstract: Let $(X,T)$ be a topological dynamical system with metric $d$. We define a new function $\overline{F}(x,y)=\limsup\limits_{n \to +\infty} \inf\limits_{\sigma \in S_n} \frac 1n \sum\limits_{k=1}^n d(T^k x,T^{\sigma(k)} y)$ by using permutation group $S_n$. It's shown $F(x,y)=\lim\limits_{n \to +\infty} \inf\limits_{\sigma \in S_n} \frac 1n \sum\limits_{k=1}^n d(T^k x,T^{\sigma(k)} y)$ exists when $x,y \in X$ are generic points. Applying this function, we prove $(X,T)$ is uniquely ergodic if and only if $\overline{F}(x,y)=0$ for any $x,y \in X$. The characterizations of ergodic measures and physical measures by $\overline{F}(x,y)$ are given. We introduce the notion of weak mean equicontinuity and prove that $(X,T)$ is weak mean equicontinuous if and only if the time averages $f^{*}(x)=\lim\limits_{n \to +\infty}\frac 1n \sum\limits_{k=1}^n f(T^k x)$ exist and are continuous for all $f \in C(X)$.


title: Connected max cut is polynomial for graphs without $K_5\backslash e$ as a minor

abstract: Given a graph $G=(V, E)$, a connected cut $\delta (U)$ is the set of edges of E linking all vertices of U to all vertices of $V\backslash U$ such that the induced subgraphs $G[U]$ and $G[V\backslash U]$ are connected. Given a positive weight function $w$ defined on $E$, the connected maximum cut problem (CMAX CUT) is to find a connected cut $\Omega$ such that $w(\Omega)$ is maximum among all connected cuts. CMAX CUT is NP-hard even for planar graphs. In this paper, we prove that CMAX CUT is polynomial for graphs without $K_5\backslash e$ as a minor. We deduce a quadratic time algorithm for the minimum cut problem in the same class of graphs without computing the maximum flow.


title: The "Sphered Cube": A New Method for the Solution of Partial Differential Equations in Cubical Geometry

abstract: A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere. By choosing to embed the cube inside of the sphere, one obtains a discretization that is free of any sharp edges or corners. Taking full advantage of the simple geometry of the sphere, spectral bases based on spin-weighted spherical harmonics and Jacobi polynomials, which properly capture the regularity of scalar, vector and tensor components in spherical coordinates, can be applied to obtain moderately efficient and accurate numerical solutions of partial differential equations in the cube. This technique demonstrates the advantages of these bases over other methods for solving PDEs in spherical coordinates. We present results for a test case of incompressible hydrodynamics in cubical geometry: Rayleigh-B\'enard convection with fully Dirichlet boundary conditions. Analysis of the simulations provides what is, to our knowledge, the first result on the scaling of the heat flux with the thermal forcing for this type of convection in a cube in a sphere.

categories: physics.comp-ph astro-ph.IM math.NA physics.flu-dyn


title: The Marshland Conjecture

abstract: We posit the existence of the Marshland within string theory. This region is the boundary between the landscape of consistent low-energy limits of quantum gravity, and the swampland of theories that cannot be embedded within string theory because they violate certain trendy and obviously uncontroversial conjectures. The Marshland is probably fractal, and we show some pretty pictures of fractals that will be useful in talks. We further show that the Marshland contains theories with a large number of light axions, allowing us to cite lots of our own papers. We show that the Marshland makes up most of the volume of the landscape, and admits a novel, weakly broken $\mathbb{Z}_2$ Marshymmetry that we find strong evidence for by considering a carefully crafted example.


title: Optimising maintenance: What are the expectations for Cyber Physical Systems

abstract: The need for maintenance is based on the wear of components of machinery. If this need can be defined reliably beforehand so that no unpredicted failures take place then the maintenance actions can be carried out economically with mini-mum disturbances to production. There are two basic challenges in solving the above. First understanding the development of wear and failures, and second managing the measurement and diagnosis of such parameters that can reveal the development of wear. In principle the development of wear and failures can be predicted through monitoring time, load or wear as such. Moni-toring time is not very efficient, as there are only limited numbers of components that suffer from aging which as such is the result of chemical wear i.e. changes in the material. In most cases the loading of components influences their wear. In principle the loading can be stable or varying in nature. Of these two cases the varying load case is much more challenging than the stable one. The monitoring of wear can be done either directly e.g. optical methods or indirectly e.g. vibration. Monitoring actual wear is naturally the most reliable approach, but it often means that additional investments are needed. The paper discusses how the monitoring of wear and need for maintenance can be done based on the use of Cyber Physical Systems.


title: Complex manifolds with negative curvature operator

abstract: We prove that compact complex manifolds with admitting metrics with negative Chern curvature operator either admit a $dd^c$-exact positive (1,1) current, or are K\"ahler with ample canonical bundle. In the case of complex surfaces we obtain a complete classification. The proofs rely on a global existence and convergence result for the pluriclosed flow.


title: On the Bohr inequality with a fixed zero coefficient

abstract: In this paper, we introduce the study of the Bohr phenomenon for a quasi-subordination family of functions, and establish the classical Bohr's inequality for the class of quasisubordinate functions. As a consequence, we improve and obtain the exact version of the classical Bohr's inequality for bounded analytic functions and also for $K$-quasiconformal harmonic mappings by replacing the constant term by the absolute value of the analytic part of the given function. We also obtain the Bohr radius for the subordination family of odd analytic functions.


title: Derived categories of (one-sided) exact categories and their localizations

abstract: We consider the quotient of an exact or one-sided exact category $\mathcal{E}$ by a so-called percolating subcategory $\mathcal{A}$. For exact categories, such a quotient is constructed in two steps. Firstly, one localizes $\mathcal{E}$ at a suitable class $S_\mathcal{A} \subseteq \operatorname{Mor}(\mathcal{E})$ of morphisms. The localization $\mathcal{E}[S_\mathcal{A}^{-1}]$ need not be an exact category, but will be a one-sided exact category. Secondly, one constructs the exact hull $\mathcal{E}{/\mkern-6mu/} \mathcal{A}$ of $\mathcal{E}[S_\mathcal{A}^{-1}]$ and shows that this satisfies the 2-universal property of a quotient amongst exact categories. In this paper, we show that this quotient $\mathcal{E} \to \mathcal{E} {/\mkern-6mu/} \mathcal{A}$ induces a Verdier localization $\mathbf{D}^b(\mathcal{E}) \to \mathbf{D}^b(\mathcal{E} {/\mkern-6mu/} \mathcal{A})$ of bounded derived categories. Specifically, (i) we study the derived category of a one-sided exact category, (ii) we show that the localization $\mathcal{E} \to \mathcal{E}[S_\mathcal{A}^{-1}]$ induces a Verdier quotient $\mathbf{D}^b(\mathcal{E}) \to \mathbf{D}^b(\mathcal{E}[S^{-1}_\mathcal{A}])$, and (iii) we show that the natural embedding of a one-sided exact category $\mathcal{F}$ into its exact hull $\overline{\mathcal{F}}$ lifts to a derived equivalence $\mathbf{D}^b(\mathcal{F}) \to \mathbf{D}^b(\overline{\mathcal{F}})$. We furthermore show that the Verdier localization is compatible with several enhancements of the bounded derived category, so that the above Verdier localization can be used in the study of localizing invariants, such as non-connective $K$-theory.


title: Overcoming Catastrophic Forgetting with Unlabeled Data in the Wild

abstract: Lifelong learning with deep neural networks is well-known to suffer from catastrophic forgetting: the performance on previous tasks drastically degrades when learning a new task. To alleviate this effect, we propose to leverage a large stream of unlabeled data easily obtainable in the wild. In particular, we design a novel class-incremental learning scheme with (a) a new distillation loss, termed global distillation, (b) a learning strategy to avoid overfitting to the most recent task, and (c) a confidence-based sampling method to effectively leverage unlabeled external data. Our experimental results on various datasets, including CIFAR and ImageNet, demonstrate the superiority of the proposed methods over prior methods, particularly when a stream of unlabeled data is accessible: our method shows up to 15.8% higher accuracy and 46.5% less forgetting compared to the state-of-the-art method. The code is available at https://github.com/kibok90/iccv2019-inc.


title: The effects of dust evolution on disks in the mid-IR

abstract: In this paper, we couple together the dust evolution code two-pop-py with the thermochemical disk modelling code ProDiMo. We create a series of thermochemical disk models that simulate the evolution of dust over time from 0.018 Myr to 10 Myr, including the radial drift, growth, and settling of dust grains. We examine the effects of this dust evolution on the mid-infrared gas emission, focussing on the mid-infrared spectral lines of C2H2, CO2, HCN, NH3, OH, and H2O that are readily observable with Spitzer and the upcoming E-ELT and JWST. The addition of dust evolution acts to increase line fluxes by reducing the population of small dust grains. We find that the spectral lines of all species except C2H2 respond strongly to dust evolution, with line fluxes increasing by more than an order of magnitude across the model series as the density of small dust grains decreases over time. The C2H2 line fluxes are extremely low due to a lack of abundance in the infrared line-emitting regions, despite C2H2 being commonly detected with Spitzer, suggesting that warm chemistry in the inner disk may need further investigation. Finally, we find that the CO2 flux densities increase more rapidly than the other species as the dust disk evolves. This suggests that the flux ratios of CO2 to other species may be lower in disks with less-evolved dust populations.


title: Yet Another Accelerated SGD: ResNet-50 Training on ImageNet in 74.7 seconds

abstract: There has been a strong demand for algorithms that can execute machine learning as faster as possible and the speed of deep learning has accelerated by 30 times only in the past two years. Distributed deep learning using the large mini-batch is a key technology to address the demand and is a great challenge as it is difficult to achieve high scalability on large clusters without compromising accuracy. In this paper, we introduce optimization methods which we applied to this challenge. We achieved the training time of 74.7 seconds using 2,048 GPUs on ABCI cluster applying these methods. The training throughput is over 1.73 million images/sec and the top-1 validation accuracy is 75.08%.


title: Adiabatic population transfer of dressed spin states with quantum optimal control

abstract: We report theoretical studies of adiabatic population transfer using dressed spin states. Quantum optimal control using the algorithm of Chopped Random Basis (CRAB) has been implemented in a negatively charged diamond nitrogen vacancy center that is coupled to a strong and resonant microwave field. We show that the dressed spin states are highly effective in suppressing effects of spin dephasing on adiabatic population transfer. The numerical simulation also demonstrates that CRAB-based quantum optimal control can enable an efficient and robust adiabatic population transfer.


title: Novel effective ergodicity breaking phase transition in a driven-dissipative system

abstract: We show that the Olami-Feder-Christensen model exhibits an effective ergodicity breaking transition as the noise is varied. Above the critical noise, the average stress on each site converges to the global average. Below the critical noise, the stress on individual sites becomes trapped in different limit cycles. We use ideas from the study of dynamical systems and compute recurrence plots and the recurrence rate. We identify the order parameter as the recurrence rate averaged over all sites and find numerical evidence that the transition can be characterized by exponents that are consistent with hyperscaling.


title: Fooling the Parallel Or Tester with Probability $8/27$

abstract: It is well-known that the higher-order language PCF is not fully abstract: there is a program - the so-called parallel or tester, meant to test whether its input behaves as a parallel or - which never terminates on any input, operationally, but is denotationally non-trivial. We explore a probabilistic variant of PCF, and ask whether the parallel or tester exhibits a similar behavior there. The answer is no: operationally, one can feed the parallel or tester an input that will fool it into thinking it is a parallel or. We show that the largest probability of success of such would-be parallel ors is exactly $8/27$. The bound is reached by a very simple probabilistic program. The difficult part is to show that that bound cannot be exceeded.


title: Structural modification of thin Bi(1 1 1) films by passivation and native oxide model

abstract: The structure of thin terminated Bi(1 1 1) films of approximately 1 nm thickness is investigated from first principles. Our density functional theory calculations show that covalent bonds to the surface can change the orientation of the films completely. For thicker films, the effect is limited to the surface only. Based on these observations, we further present a simple model structure for the native oxide and chemically similar oxides, which form a protective capping layer, leaving the orientation of the films unchanged. The advantages of this energetically favorable layered termination are discussed in the context of the films' technological exploitation in nanoelectronic devices.


title: Engineering problems in machine learning systems

abstract: Fatal accidents are a major issue hindering the wide acceptance of safety-critical systems that employ machine learning and deep learning models, such as automated driving vehicles. In order to use machine learning in a safety-critical system, it is necessary to demonstrate the safety and security of the system through engineering processes. However, thus far, no such widely accepted engineering concepts or frameworks have been established for these systems. The key to using a machine learning model in a deductively engineered system is decomposing the data-driven training of machine learning models into requirement, design, and verification, particularly for machine learning models used in safety-critical systems. Simultaneously, open problems and relevant technical fields are not organized in a manner that enables researchers to select a theme and work on it. In this study, we identify, classify, and explore the open problems in engineering (safety-critical) machine learning systems --- that is, in terms of requirement, design, and verification of machine learning models and systems --- as well as discuss related works and research directions, using automated driving vehicles as an example. Our results show that machine learning models are characterized by a lack of requirements specification, lack of design specification, lack of interpretability, and lack of robustness. We also perform a gap analysis on a conventional system quality standard SQuARE with the characteristics of machine learning models to study quality models for machine learning systems. We find that a lack of requirements specification and lack of robustness have the greatest impact on conventional quality models.


title: Longitudinal spin polarization in a thermal model

abstract: We use a thermal model with single freeze-out to determine longitudinal polarization of $\Lambda$ hyperons emitted from a hot and rotating hadronic medium. We consider the top RHIC energies and use the model parameters determined in the previous analyses of particle spectra and elliptic flow. Using a direct connection between the spin polarization tensor and thermal vorticity, we reproduce earlier results which indicate a quadrupole structure of the longitudinal component of the polarization three-vector with an opposite sign compared to that found in the experiment. We further use only the spatial components of the thermal vorticity in the laboratory system to define polarization and show that this leads to the correct sign and magnitude of the quadrupole structure. This procedure resembles a non-relativistic connection between the polarization three-vector and vorticity employed in other works. In general, our results bring further evidence that the spin polarization dynamics in heavy-ion collisions may be not directly related to the thermal vorticity. The additional material explains the construction of the hydrodynamicaly consistent gradients of fluid velocity and temperature in thermal models with the help of the perfect-fluid equations of motion.


title: Firsthand Opiates Abuse on Social Media: Monitoring Geospatial Patterns of Interest Through a Digital Cohort

abstract: In the last decade drug overdose deaths reached staggering proportions in the US. Besides the raw yearly deaths count that is worrisome per se, an alarming picture comes from the steep acceleration of such rate that increased by 21% from 2015 to 2016. While traditional public health surveillance suffers from its own biases and limitations, digital epidemiology offers a new lens to extract signals from Web and Social Media that might be complementary to official statistics. In this paper we present a computational approach to identify a digital cohort that might provide an updated and complementary view on the opioid crisis. We introduce an information retrieval algorithm suitable to identify relevant subspaces of discussion on social media, for mining data from users showing explicit interest in discussions about opioid consumption in Reddit. Moreover, despite the pseudonymous nature of the user base, almost 1.5 million users were geolocated at the US state level, resembling the census population distribution with a good agreement. A measure of prevalence of interest in opiate consumption has been estimated at the state level, producing a novel indicator with information that is not entirely encoded in the standard surveillance. Finally, we further provide a domain specific vocabulary containing informal lexicon and street nomenclature extracted by user-generated content that can be used by researchers and practitioners to implement novel digital public health surveillance methodologies for supporting policy makers in fighting the opioid epidemic.


title: Tweeting MPs: Digital Engagement between Citizens and Members of Parliament in the UK

abstract: Disengagement and disenchantment with the Parliamentary process is an important concern in today's Western democracies. Members of Parliament (MPs) in the UK are therefore seeking new ways to engage with citizens, including being on digital platforms such as Twitter. In recent years, nearly all (579 out of 650) MPs have created Twitter accounts, and have amassed huge followings comparable to a sizable fraction of the country's population. This paper seeks to shed light on this phenomenon by examining the volume and nature of the interaction between MPs and citizens. We find that although there is an information overload on MPs, attention on individual MPs is focused during small time windows when something topical may be happening relating to them. MPs manage their interaction strategically, replying selectively to UK-based citizens and thereby serving in their role as elected representatives, and using retweets to spread their party's message. Most promisingly, we find that Twitter opens up new avenues with substantial volumes of cross-party interaction, between MPs of one party and citizens who support (follow) MPs of other parties.


title: Multiplicity dependence for the production of strange hadrons and charged particles in proton-proton collisions

abstract: In this contribution, the production rates and the transverse momentum distributions of strange hadrons are reported as a function of charged particle multiplicity. In this analysis, the data collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector at the LHC are used. It is found that the production rate of $K_{S}^{0}$, $\Lambda$, $\Xi^{\pm}$, and $\Omega$ increases with multiplicity faster than that for charged particles. The higher the strangeness content of the hadron, the more pronounced is the increase. Moreover, the energy and multiplicity dependence of charged particle production in pp collisions are presented and the results are compared to predictions from Monte Carlo (MC) event generators. It turns out that the average multiplicity density increases steeply with center-of-mass energy for high multiplicity classes.


title: The Dynamical equations for $\mathfrak{gl}(n\vert m)$

abstract: In this note we propose a compatible set of equations which commutes with the Knizhnik-Zamolodchikov equations based on the $\mathfrak{{gl}}(n\vert m)$ symmetry algebra and establish the Matsuo-Cherednik correspondence in this context.


title: News vs Information

abstract: We consider the relative entropy between the vacuum state and a coherent state in linearized quantum gravity around a stationary black hole spacetime. Combining recent results by Casini et al. and Longo with the Raychaudhuri equation, the following result is obtained: Let $\frak A$ be the algebra of observables assoiciated with a region that is the causal future of some compact set in the interior of the spacetime. Let $S$ be the relative entropy with respect to this algebra, $A$ the area of the horizon cross section defined by the region, computed to second order in the gravitational perturbation. If the region is time-translated by the Killing parameter $t$, then $\frac{d}{d t}(S+A/4)= 2\pi F$, with $F$ the flux of the gravitational/matter radiation (integrated squared news tensor) emitted towards the future of the region.


title: Quadrotor Manipulation System: Development of a Robust Contact Force Estimation and Impedance Control Scheme Based on DOb and FTRLS

abstract: The research on aerial manipulation systems has been increased rapidly in recent years. These systems are very attractive for a wide range of applications due to their unique features. However, dynamics, control and manipulation tasks of such systems are quite challenging because they are naturally unstable, have very fast dynamics, have strong nonlinearities, are very susceptible to parameters variations due to carrying a payload besides the external disturbances, and have complex inverse kinematics. In addition, the manipulation tasks require estimating (applying) a certain force of (at) the end-effector as well as the accurate positioning of it. Thus, in this article, a robust force estimation and impedance control scheme is proposed to address these issues. The robustness is achieved based on the Disturbance Observer (DOb) technique. Then, a tracking and performance low computational linear controller is used. For teleoperation purpose, the contact force needs to be identified. However, the current developed techniques for force estimation have limitations because they are based on ignoring some dynamics and/or requiring of an indicator of the environment contact. Unlike these techniques, we propose a technique based on linearization capabilities of DOb and a Fast Tracking Recursive Least Squares (FTRLS) algorithm. The complex inverse kinematics problem of such a system is solved by a Jacobin based algorithm. The stability analysis of the proposed scheme is presented. The algorithm is tested to achieve tracking of task space reference trajectories besides the impedance control. The efficiency of the proposed technique is enlightened via numerical simulation.


title: Reconstructing Rational Functions with $\texttt{FireFly}$

abstract: We present the open-source $\texttt{C++}$ library $\texttt{FireFly}$ for the reconstruction of multivariate rational functions over finite fields. We discuss the involved algorithms and their implementation. As an application, we use $\texttt{FireFly}$ in the context of integration-by-parts reductions and compare runtime and memory consumption to a fully algebraic approach with the program $\texttt{Kira}$.


title: Deconfined quantum critical point in one dimension

abstract: We perform a numerical study of a spin-1/2 model with $\mathbb{Z}_2 \times \mathbb{Z}_2$ symmetry in one dimension which demonstrates an interesting similarity to the physics of two-dimensional deconfined quantum critical points (DQCP). Specifically, we investigate the quantum phase transition between Ising ferromagnetic and valence bond solid (VBS) symmetry-breaking phases. Working directly in the thermodynamic limit using uniform matrix product states, we find evidence for a direct continuous phase transition that lies outside of the Landau-Ginzburg-Wilson paradigm. In our model, the continuous transition is found everywhere on the phase boundary. We find that the magnetic and VBS correlations show very close power law exponents, which is expected from the self-duality of the parton description of this DQCP. Critical exponents vary continuously along the phase boundary in a manner consistent with the predictions of the field theory for this transition. We also find a regime where the phase boundary splits, as suggested by the theory, introducing an intermediate phase of coexisting ferromagnetic and VBS order parameters. Interestingly, we discover a transition involving this coexistence phase which is similar to the DQCP, being also disallowed by Landau-Ginzburg-Wilson symmetry-breaking theory.


title: Measuring Star-Formation Histories, Distances, and Metallicities with Pixel Color-Magnitude Diagrams I: Model Definition and Mock Tests

abstract: We present a comprehensive study of the applications of the pixel color-magnitude diagram (pCMD) technique for measuring star formation histories (SFHs) and other stellar population parameters of galaxies, and demonstrate that the technique can also constrain distances. SFHs have previously been measured through either the modeling of resolved-star CMDs or of integrated-light SEDs, yet neither approach can easily be applied to galaxies in the "semi-resolved regime". The pCMD technique has previously been shown to have the potential to measure stellar populations and star formation histories in semi-resolved galaxies. Here we present Pixel Color-Magnitude Diagrams with Python (pcmdpy), a GPU-accelerated package that makes significant computational improvements to the original code and including more realistic physical models. These advances include the simultaneous fitting of distance, modeling a Gaussian metallicity-distribution function, and an observationally-motivated dust model. GPU-acceleration allows these more realistic models to be fit roughly 7x faster than the simpler models in the original code. We present results from a suite of mock tests, showing that with proper model assumptions, the code can simultaneously recover SFH, [Fe/H], distance, and dust extinction. Our results suggest the code, applied to observations with HST-like resolution, should constrain these properties with high precision within 10 Mpc and can be applied to systems out to as far as 100 Mpc. pCMDs open a new window to studying the stellar populations of many galaxies that cannot be readily studied through other means.


title: A massive white-dwarf merger product before final collapse

abstract: Gravitational wave emission can lead to the coalescence of close pairs of compact objects orbiting each other. For the case of neutron stars such mergers may yield masses above the Tolman-Oppenheimer-Volkoff limit, leading to the formation of black holes. For the case of white dwarfs the merger product may exceed the Chandrasekhar limit, leading either to a thermonuclear explosion as Type Ia supernova, or to a collapse forming a neutron star. If a Type Ia supernova explosion is avoided, the merger of two massive white dwarfs is expected to form a hydrogen- and helium-free circumstellar nebula with a hot and luminous, rapidly rotating and highly magnetized central star for several 10,000 yr before its final collapse. Here we report the discovery of a hot star with an emission line dominated spectrum in the centre of a circular mid-infrared nebula. Both the star and the nebula appear to be free of hydrogen and helium. Our tailored stellar atmosphere and wind models indicate a stellar surface temperature of about 200,000 K, and a record outflow velocity of 16,000 km/s. This extreme velocity, together with the derived mass outflow rate, imply rapid stellar rotation and a strong magnetic field aiding the wind acceleration. The Gaia distance of the star leads to a luminosity of 10^{4.5} Lsun, which matches models of the post-merger evolution of super-Chandrasekhar mass white dwarfs. The high stellar temperature and the nebular size argue for a short remaining lifetime of the star, which will produce a bright optical and high-energy transient upon collapse. Our observations indicate that super-Chandrasekhar mass white dwarf mergers can indeed avoid a thermonuclear explosion as Type Ia supernova, and provide empirical evidence for magnetic field generation in stellar mergers.


title: Magnus Hall Effect

abstract: A new type of linear response Hall effect is predicted in time-reversal-invariant systems with built-in electric field at zero magnetic field. The Hall response results from a quantum Magnus effect where a self-rotating Bloch electron wavepacket moving under electric field develops an anomalous velocity in the transverse direction. We show that in the ballistic limit the Magnus Hall conductance measures the distribution of Berry curvature on the Fermi surface.


title: RAPID: Early Classification of Explosive Transients using Deep Learning

abstract: We present RAPID (Real-time Automated Photometric IDentification), a novel time-series classification tool capable of automatically identifying transients from within a day of the initial alert, to the full lifetime of a light curve. Using a deep recurrent neural network with Gated Recurrent Units (GRUs), we present the first method specifically designed to provide early classifications of astronomical time-series data, typing 12 different transient classes. Our classifier can process light curves with any phase coverage, and it does not rely on deriving computationally expensive features from the data, making RAPID well-suited for processing the millions of alerts that ongoing and upcoming wide-field surveys such as the Zwicky Transient Facility (ZTF), and the Large Synoptic Survey Telescope (LSST) will produce. The classification accuracy improves over the lifetime of the transient as more photometric data becomes available, and across the 12 transient classes, we obtain an average area under the receiver operating characteristic curve of 0.95 and 0.98 at early and late epochs, respectively. We demonstrate RAPID's ability to effectively provide early classifications of observed transients from the ZTF data stream. We have made RAPID available as an open-source software package (https://astrorapid.readthedocs.io) for machine learning-based alert-brokers to use for the autonomous and quick classification of several thousand light curves within a few seconds.


title: Entanglement entropy of the long-range Dyson hierarchical model

abstract: We study the ground state entanglement entropy of the quantum Dyson hierarchical spin chain in which the interaction decays algebraically with the distance as $r^{-1-\sigma}$. We exploit the real-space renormalisation group solution which gives the ground-state wave function in the form of a tree tensor network and provides a manageable recursive expression for the reduced density matrix of the renormalised ground state. Surprisingly, we find that at criticality the entanglement entropy obeys an area law, as opposite to the logarithmic scaling of short-range critical systems and of other non-hierarchical long-range models. We provide also some analytical results in the limit of large and small $\sigma$ that are tested against the numerical solution of the recursive equations.


title: Photon pair condensation by engineered dissipation

abstract: Dissipation can usually induce detrimental decoherence in a quantum system. However, engineered dissipation can be used to prepare and stabilize coherent quantum many-body states. Here, we show that by engineering dissipators containing photon pair operators, one can stabilize an exotic dark state, which is a condensate of photon pairs with a phase-nematic order. In this system, the usual superfluid order parameter, i.e. single-photon correlation, is absent, while the photon pair correlation exhibits long-range order. Although the dark state is not unique due to multiple parity sectors, we devise an additional type of dissipators to stabilize the dark state in a particular parity sector via a diffusive annihilation process which obeys Glauber dynamics in an Ising model. Furthermore, we propose an implementation of these photon-pair dissipators in circuit-QED architecture.


title: Bootstrapping MN and Tetragonal CFTs in Three Dimensions

abstract: Conformal field theories (CFTs) with MN and tetragonal global symmetry in $d=2+1$ dimensions are relevant for structural, antiferromagnetic and helimagnetic phase transitions in a wide class of materials. The study of these theories with the nonperturbative numerical conformal bootstrap is initiated in this work. Bounds for operator dimensions are obtained and they are found to possess sharp kinks in the MN case, suggesting the existence of full-fledged CFTs. In the tetragonal case, no new kinks are found. Estimates for critical exponents are provided for a few cases describing phase transitions in various materials. In two particular MN cases, corresponding to theories with global symmetry groups $O(2)^2\rtimes S_2$ and $O(2)^3\rtimes S_3$, a second kink is found. In the $O(2)^2\rtimes S_2$ case it is argued to be saturated by a CFT that belongs to a new universality class relevant for the structural phase transition of NbO$_2$ and paramagnetic-helimagnetic transitions of the rare-earth metals Ho and Dy. In the $O(2)^3\rtimes S_3$ case it is suggested that the CFT that saturates the second kink belongs to a new universality class relevant for the paramagnetic-antiferromagnetic phase transition of the rare-earth metal Nd.


title: Nuclear Parton Distributions from Lepton-Nucleus Scattering and the Impact of an Electron-Ion Collider

abstract: We present a first determination of the nuclear parton distribution functions (nPDF) based on the NNPDF methodology: nNNPDF1.0. This analysis is based on neutral-current deep-inelastic structure function data and is performed up to NNLO in QCD calculations with heavy quark mass effects. For the first time in the NNPDF fits, the $\chi^2$ minimization is achieved using stochastic gradient descent with reverse-mode automatic differentiation (backpropagation). We validate the robustness of the fitting methodology through closure tests, assess the perturbative stability of the resulting nPDFs, and compare them with other recent analyses. The nNNPDF1.0 distributions satisfy the boundary condition whereby the NNPDF3.1 proton PDF central values and uncertainties are reproduced at $A=1$, which introduces important constraints particularly for low-$A$ nuclei. We also investigate the information that would be provided by an Electron-Ion Collider (EIC), finding that EIC measurements would significantly constrain the nPDFs down to $x\simeq 5\times 10^{-4}$. Our results represent the first-ever nPDF determination obtained using a Monte Carlo methodology consistent with that of state-of-the-art proton PDF fits, and provide the foundation for a subsequent global nPDF analyses including also proton-nucleus data.


title: Making Trotters Sprint: A Variational Imaginary Time Ansatz for Quantum Many-body Systems

abstract: We introduce a variational wavefunction for many-body ground states that involves imaginary time evolution with two different Hamiltonians in an alternating fashion with variable time intervals. We successfully apply the ansatz on the one- and two-dimensional transverse-field Ising model and systematically study its scaling for the one-dimensional model at criticality. We find the total imaginary time required scales logarithmically with system size, in contrast to the linear scaling in conventional Quantum Monte Carlo. We suggest this is due to unique dynamics permitted by alternating imaginary time evolution, including the exponential growth of bipartite entanglement. For generic models, the superior scaling of our ansatz potentially mitigates the negative sign problem at the expense of having to optimize variational parameters.


title: Searches for other vacua I: bubbles in our universe

abstract: We discuss models in which vacua other than our own can be directly observed in the present universe. Models with density-dependent vacuum structure can give rise to `non-lethal'-vacua: vacua with lower energy-density than our vacuum, but only in regions with finite Standard Model densities. These models provide an explicit example of a bubble which is confined to a finite region of space and produces potentially detectable signatures, unlike standard Coleman tunneling events where bubbles expand at the speed of light and are never directly observable. We study the expansion and contraction of a confined bubble created after a core-collapse supernova, focusing on energy deposition that may be observable in the vicinity of a supernova remnant due to the formation and evolution of a confined bubble.


title: Emergent Symmetry and Conserved Current at a One Dimensional Incarnation of Deconfined Quantum Critical Point

abstract: The deconfined quantum critical point (DQCP) was originally proposed as a continuous transition between two spontaneous symmetry breaking phases in 2D spin-1/2 systems. While great efforts have been spent on the DQCP for 2D systems, both theoretically and numerically, ambiguities among the nature of the transition are still not completely clarified. Here we shift the focus to a recently proposed 1D incarnation of DQCP in a spin-1/2 chain. By solving it with the variational matrix product state in the thermodynamic limit, a continuous transition between a valence-bond solid phase and a ferromagnetic phase is discovered. The scaling dimensions of various operators are calculated and compared with those from field theoretical description. At the critical point, two emergent O(2) symmetries are revealed, and the associated conserved current operators with exact integer scaling dimensions are determined with scrutiny. Our findings provide the low-dimensional analog of DQCP where unbiased numerical results are in perfect agreement with the controlled field theoretical predictions and have extended the realm of the unconventional phase transition as well as its identification with the advanced numerical methodology.


title: Charting the coming synergy between lattice QCD and high-energy phenomenology

abstract: Building upon the PDFSense framework developed in Ref. [1], we perform a comprehensive analysis of the sensitivity of present and future high-energy data to a number of quantities commonly evaluated in lattice gauge theory, with a particular focus on the integrated Mellin moments of nucleon parton distribution functions (PDFs), such as $\langle x \rangle_{u^+ - d^+}$ and $\langle x \rangle_{g}$, as well as $x$-dependent quark quasi-distributions -- in particular, that of the isovector combination. Our results demonstrate the potential for lattice calculations and phenomenological quark distributions informed by high-energy experimental data to cooperatively improve the picture of the nucleon's collinear structure. This will increasingly be the case as computational resources for lattice calculations further expand, and QCD global analyses continue to grow in sophistication. Our sensitivity analysis suggests that a future lepton-hadron collider would be especially instrumental in providing phenomenological constraints to lattice observables.


title: The evolution of the quenching of star formation in cluster galaxies since z $\sim$ 1

abstract: We study the star-forming (SF) population of galaxies within a sample of 209 IR-selected galaxy clusters at 0.3$\,\leq\,z\,\leq\,$1.1 in the ELAIS-N1 and XMM-LSS fields, exploiting the first HSC-SSP data release. The large area and depth of these data allows us to analyze the dependence of the SF fraction, $f_{SF}$, on stellar mass and environment separately. Using $R/R_{200}$ to trace environment, we observe a decrease in $f_{SF}$ from the field towards the cluster core, which strongly depends on stellar mass and redshift. The data show an accelerated growth of the quiescent population within the cluster environment: the $f_{SF}$ vs. stellar mass relation of the cluster core ($R/R_{200}\,\leq\,$0.4) is always below that of the field (4$\,\leq\,R/R_{200}\,<\,$6). Finally, we find that environmental and mass quenching efficiencies depend on galaxy stellar mass and distance to the center of the cluster, demonstrating that the two effects are not separable in the cluster environment. We suggest that the increase of the mass quenching efficiency in the cluster core may emerge from an initial population of galaxies formed ``in situ.'' The dependence of the environmental quenching efficiency on stellar mass favors models in which galaxies exhaust their reservoir of gas through star formation and outflows, after new gas supply is truncated when galaxies enter the cluster.


title: BAO Modulation as a Probe of Compensated Isocurvature Perturbations

abstract: Compensated isocurvature perturbations (CIPs) are opposite spatial fluctuations in the baryon and dark matter density. They can be generated for example in the curvaton model in the early Universe but are difficult to observe because their gravitational imprint nearly cancels. We therefore propose a new measurement method by searching for a spatial modulation of the baryon acoustic oscillation (BAO) scale that CIPs induce. We find that for a Euclid-like survey the sensitivity is marginally better than the WMAP cosmic microwave background (CMB) constraint, which exploits the CIP-induced modulation of the CMB sound horizon. For a cosmic-variance limited BAO survey using emission-line galaxies up to $z\sim7$ the sensitivity is between stage 3 and stage 4 CMB experiments. These results include using CIP-galaxy cross-correlations, which improves the sensitivity by a factor of $\sim2-3$ for correlated CIPs. The method could be further improved with an optimal estimator, similarly to the CMB, and could provide a useful cross-check of other CIP probes. Finally, if CIPs exist, they can bias cosmological measurements made assuming no CIPs. In particular, they can act as a super-sample fluctuation of the baryon density and bias measurements of the BAO scale. For modern BAO surveys, the largest 2$\sigma$ CIP fluctuation allowed by Planck's 95\% bound could bias BAO measurements of $H(z)$ by 2.2\%, partially reducing the tension with the local $H_0$ measurements from 3.1$\sigma$ to 2.3$\sigma$.


title: Exploring A Photospheric Radius Correction to Model Secondary Eclipse Spectra for Transiting Exoplanets

abstract: We highlight a physical effect that is often not considered that impacts the calculation of model spectra of planets at secondary eclipse, affecting both emission and reflection spectra. The radius of the emitting surface of the planet is not merely one value measured from a transit light curve, but is itself a function of wavelength, yet it is not directly measurable. At high precision, a similar effect is well-known in transit "transmission spectroscopy" but this related effect also impacts emission and reflection. As is well-appreciated, the photospheric radius can vary across $\sim$4-8 atmospheric scale heights, depending on atmospheric opacity and spectral resolution. This effect leads to a decreased weighting in model calculations at wavelengths where atmospheric opacity is low, and one sees more deeply into the atmosphere, to a smaller radius. The overall effect serves to mute emission spectra features for atmospheres with no thermal inversion but to enhance features for atmospheres with a thermal inversion. While this effect can be ignored for current \emph{Hubble} observations, it can lead to wavelength-dependent 10-20\% changes in planet-to-star flux ratios in the infrared at $R\sim~200-1000$ (readily achievable for JWST) for low-gravity hot Jupiters, although values of 5\% are more typical for the population. The effect is mostly controlled by the ratio of the atmospheric scale height to the planet radius, and can be important at any planetary temperature. Of known planets, the effect is largest for the cool "super-puffs" at very low surface gravity, where it can alter calculated flux ratios by over 100\%. We discuss complexities of including this photospheric radius effect in 1D and 3D atmosphere models.


title: A Golden Probe of Nonlinear Higgs Dynamics

abstract: The most salient generic feature of a composite Higgs boson resides in the nonlinearity of its dynamics, which arises from degenerate vacua associated with the pseudo-Nambu-Goldstone (PNGB) nature of the Higgs boson. It has been shown that the nonlinear Higgs dynamics is universal in the IR and controlled only by a single parameter $f$, the decay constant of the PNGB Higgs. In this work we perform a fit, for the first time, to Wilson coefficients of ${\cal O}(p^4)$ operators in the nonlinear Lagrangian using the golden H $\to$ 4L decay channel. By utilizing both the "rate" information in the signal strength and the "shape" information in the fully differential spectra, we provide limits on the Goldstone decay constant $f$, as well as ${\cal O}(p^4)$ Wilson coefficients, using Run 2 data at the LHC. In rate measurements alone, the golden channel prefers a negative $\xi =v^2/f^2$ corresponding to a non-compact coset structure. Including the shape information, we identify regions of parameter space where current LHC constraint on $f$ is still weak, allowing for $\xi \lesssim 0.5$ or $\xi \gtrsim -0.5$. We also comment on future sensitivity at the high-luminosity upgrade of the LHC which could allow for simultaneous fits to multiple Wilson coefficients.


title: The Mass of Stirring Bodies in the AU Mic Debris Disk Inferred from Resolved Vertical Structure

abstract: The vertical distribution of dust in debris disks is sensitive to the number and size of large planetesimals dynamically stirring the disk, and is therefore well-suited for constraining the prevalence of otherwise unobservable Uranus and Neptune analogs. Information regarding stirring bodies has previously been inferred from infrared and optical observations of debris disk vertical structure, but theoretical works predict that the small particles traced by short-wavelength observations will be `puffed up' by radiation pressure, yielding only upper limits. The large grains that dominate the disk emission at millimeter wavelengths are much less sensitive to the effects of stellar radiation or stellar winds, and therefore trace the underlying mass distribution more directly. Here we present ALMA 1.3 mm dust continuum observations of the debris disk around the nearby M star AU Mic. The 3 au spatial resolution of the observations, combined with the favorable edge-on geometry of the system, allows us to measure the vertical thickness of the disk. We report a scale height-to-radius aspect ratio of $h = 0.031_{-0.004}^{+0.005}$ between radii of $\sim 23$ au and $\sim 41$ au. Comparing this aspect ratio to a theoretical model of size-dependent velocity distributions in the collisional cascade, we find that the perturbing bodies embedded in the local disk must be larger than about 400 km, and the largest perturbing body must be smaller than roughly 1.8 M$_\odot$. These measurements rule out the presence of a gas giant or Neptune analog near the $\sim 40$ au outer edge of the debris ring, but are suggestive of large planetesimals or an Earth-sized planet stirring the dust distribution.


title: Tensor network investigation of the double layer Kagome compound Ca$_{10}$Cr$_7$O$_{28}$

abstract: Quantum spin liquids are exotic quantum phases of matter that do not order even at zero temperature. While there are several toy models and simple Hamiltonians that could host a quantum spin liquid as their ground state, it is very rare to find actual, realistic materials that exhibits their properties. At the same time, the classical simulation of such instances of strongly correlated systems is intricate and reliable methods are scarce. In this work, we investigate the quantum magnet Ca$_{10}$Cr$_7$O$_{28}$ that has recently been discovered to exhibit properties of a quantum spin liquid in inelastic neutron scattering experiments. This compound has a distorted bilayer Kagome lattice crystal structure consisting of Cr$^{5+}$ ions with spin-$1/2$ moments. Coincidentally, the lattice structure renders a tensor network algorithm in 2D applicable that can be seen as a new variant of a projected entangled simplex state algorithm in the thermodynamic limit. In this first numerical investigation of this material that takes into account genuine quantum correlations, good agreement with the experimental findings is found. We argue that this is one of the very first studies of physical materials in the laboratory with tensor network methods, contributing to uplifting tensor networks from conceptual tools to methods to describe real two-dimensional quantum materials.


title: Parametric identification of the dynamics of inter-sectoral balance: modelling and forecasting

abstract: This work is devoted to modelling and identification of the dynamics of the inter-sectoral balance of a macroeconomic system. An approach to the problem of specification and identification of a weakly formalized dynamical system is developed. A matching procedure for parameters of a linear stationary Cauchy problem with a decomposition of its upshot trend and a periodic component, is proposed. Moreover, an approach for detection of significant harmonic waves, which are inherent to real macroeconomic dynamical systems, is developed.


title: Thermodynamics and Phase Transitions of NUTty Dyons

abstract: Picking up the threads on a recent proposal [arXiv:1903.08668], we show how to formulate consistent thermodynamics of Lorentzian Taub-NUT spacetimes in the presence of electric and magnetic charges. Namely, with an entropy identified with a quarter of the horizon area and no Misner time periodicity condition imposed, we show that a new pair of conjugate quantities can be introduced so that the NUT parameter can be independently varied and the corresponding first law and Smarr relation can be consistently formulated. Moreover, we show that (contrary to the statements in the literature) the electric and magnetic parameters need not be proportional to one another and a full cohomogeneity first law including variations of both charges can be written down, provided one charge is considered on the horizon and the other at infinity. The corresponding phase transitions are also briefly discussed.