title: High-cooperativity coupling of a rare-earth spin ensemble to a superconducting resonator using yttrium orthosilicate as a substrate

abstract: Yttrium orthosilicate (Y$_2$SiO$_5$, or YSO) has proved to be a convenient host for rare-earth ions used in demonstrations of microwave quantum memories and optical memories with microwave interfaces, and shows promise for coherent microwave--optical conversion owing to its favourable optical and spin properties. The strong coupling required by such microwave applications could be achieved using superconducting resonators patterned directly on Y$_2$SiO$_5$, and hence we investigate here the use of Y$_2$SiO$_5$ as an alternative to sapphire or silicon substrates for superconducting hybrid device fabrication. A NbN resonator with frequency 6.008 GHz and low power quality factor $Q \approx 400000$ was fabricated on a Y$_2$SiO$_5$ substrate doped with isotopically enriched Nd$^{145}$. Measurements of dielectric loss yield a loss-tangent $\tan\delta = 4 \times 10^{-6}$, comparable to sapphire. Electron spin resonance (ESR) measurements performed using the resonator show the characteristic angular dependence expected from the anisotropic Nd$^{145}$ spin, and the coupling strength between resonator and electron spins is in the high cooperativity regime ($C = 30$). These results demonstrate Y$_2$SiO$_5$ as an excellent substrate for low-loss, high-Q microwave resonators, especially in applications for coupling to optically-accessible rare earth spins.


title: A quasi-robust discretization error estimate for discontinuous Galerkin Isogeometric Analysis

abstract: Isogeometric Analysis is a spline-based discretization method to partial differential equations which shows the approximation power of a high-order method. The number of degrees of freedom, however, is as small as the number of degrees of freedom of a low-order method. This does not come for free as the original formulation of Isogeometric Analysis requires a global geometry function. Since this is too restrictive for many kinds of applications, the domain is usually decomposed into patches, where each patch is parameterized with its own geometry function. In simpler cases, the patches can be combined in a conforming way. However, for non-matching discretizations or for varying coefficients, a non-conforming discretization is desired. An symmetric interior penalty discontinuous Galerkin (SIPG) method for Isogeometric Analysis has been previously introduced. In the present paper, we give error estimates that only depend poly-logarithmically on the spline degree. This opens the door towards the construction and the analysis of fast linear solvers, particularly multigrid solvers for non-conforming multipatch Isogeometric Analysis.


title: The Capacity Achieving Distribution for the Amplitude Constrained Additive Gaussian Channel: An Upper Bound on the Number of Mass Points

abstract: This paper studies an $n$-dimensional additive Gaussian noise channel with a peak-power-constrained input. It is well known that, in this case, when $n=1$ the capacity-achieving input distribution is discrete with finitely many mass points, and when $n>1$ the capacity-achieving input distribution is supported on finitely many concentric shells. However, due to the previous proof technique, neither the exact number of mass points/shells of the optimal input distribution nor a bound on it was available. This paper provides an alternative proof of the finiteness of the number mass points/shells of the capacity-achieving input distribution and produces the first firm bounds on the number of mass points and shells, paving an alternative way for approaching many such problems. Roughly, the paper consists of three parts. The first part considers the case of $n=1$. The first result, in this part, shows that the number of mass points in the capacity-achieving input distribution is within a factor of two from the downward shifted capacity-achieving output probability density function (pdf). The second result, by showing a bound on the number of zeros of the downward shifted capacity-achieving output pdf, provides a first firm upper on the number of mass points. Specifically, it is shown that the number of mass points is given by $O(\mathsf{A}^2)$ where $\mathsf{A}$ is the constraint on the input amplitude. The second part generalizes the results of the first part to the case of $n>1$. In particular, for every dimension $n>1$, it is shown that the number of shells is given by $O(\mathsf{A}^2)$ where $\mathsf{A}$ is the constraint on the input amplitude. Finally, the third part provides bounds on the number of points for the case of $n=1$ with an additional power constraint.


title: Event activity measurements and mid-rapidity correlations in 200 GeV p+Au collisions at STAR

abstract: These proceedings report preliminary measurements of correlations between mid-rapidity charged tracks and high-rapidity event activity (EA) at STAR for $\sqrt{s_\mathrm{NN}}= 200\, \mathrm{GeV}$ p+Au collisions taken in 2015. These correlations are intriguing because they inform the current debate over use of the Glauber model in `small' systems (here meaning p+A or d+A and denoted as `s+A') and have implications for calculating nuclear modification and quenching observables in these systems. The results support concerns about centrality binning in p+Au collisions, and as such motivate using ratios of semi-inclusive, as opposed to fully inclusive, jet spectra to look for jet enhancement or suppression.


title: Categories of Two-Colored Pair Partitions, Part II: Categories Indexed by Semigroups

abstract: Within the framework of unitary easy quantum groups, we study an analogue of Brauer's Schur-Weyl approach to the representation theory of the orthogonal group. We consider concrete combinatorial categories whose morphisms are formed by partitions of finite sets into disjoint subsets of cardinality two; the points of these sets are colored black or white. These categories correspond to "half-liberated easy" interpolations between the unitary group and Wang's quantum counterpart. We complete the classification of all such categories demonstrating that the subcategories of a certain natural halfway point are equivalent to additive subsemigroups of the natural numbers; the categories above this halfway point have been classified in a preceding article. We achieve this using combinatorial means exclusively. Our work reveals that the half-liberation procedure is quite different from what was previously known from the orthogonal case.


title: Nonlinear stability in nonlocal gravity

abstract: We address the stability issue of Ricci-flat and maximally symmetric spacetimes in nonlocal gravity to all perturbative orders in the gravitational perturbation. Assuming a potential at least cubic in curvature tensors but quadratic in the Ricci tensor, our proof consists on a mapping of the stability analysis in nonlocal gravity to the same problem in Einstein-Hilbert theory. One of the consequences is that only the graviton field can propagate and the theory is ghost-free at all perturbative orders. All the results known in Einstein gravity in vacuum with or without a cosmological constant can be exported to the case of nonlocal gravity: if a spacetime is stable at all perturbative orders in Einstein gravity, it is stable also in nonlocal gravity. Minkowski and de Sitter spacetimes are particular examples. We also study how the theory affects the propagation of gravitational waves in a cosmological background.


title: Thermal conductivity for III-V and II-VI semiconductor wurtzite and zinc-blende polytypes: the role of anharmonicity and phase space

abstract: We calculate the lattice thermal conductivity ($\kappa$) for cubic (zinc-blende) and hexagonal (wurtzite) phases for 8 semiconductors using $\textit{ab initio}$ calculations and solving the Phonon Boltzmann Transport Equation, explaining the different behavior of the ratio $\kappa_{\rm hex}/\kappa_{\rm cub}$ between the two phases. We show that this behavior depends on the relative importance of two antagonistic factors: anharmonicity, which we find to be always higher in the cubic phase; and the accessible phase space, which is higher for the less symmetric hexagonal phase. Based on that, we develop a method that predicts the most conducting phase---cubic or hexagonal---where other more heuristic approaches fail. We also present results for nanowires made of the same materials, showing the possibility to tune $\kappa_{\rm hex}/\kappa_{\rm cub}$ over a wide range by modifying their diameter, thus making them attractive materials for complex phononic and thermoelectric applications/systems.


title: Smoothing Spline Semiparametric Density Models

abstract: Density estimation plays a fundamental role in many areas of statistics and machine learning. Parametric, nonparametric and semiparametric density estimation methods have been proposed in the literature. Semiparametric density models are flexible in incorporating domain knowledge and uncertainty regarding the shape of the density function. Existing literature on semiparametric density models is scattered and lacks a systematic framework. In this paper, we consider a unified framework based on the reproducing kernel Hilbert space for modeling, estimation, computation and theory. We propose general semiparametric density models for both a single sample and multiple samples which include many existing semiparametric density models as special cases. We develop penalized likelihood based estimation methods and computational methods under different situations. We establish joint consistency and derive convergence rates of the proposed estimators for both the finite dimensional Euclidean parameters and an infinite-dimensional functional parameter. We validate our estimation methods empirically through simulations and an application.


title: Scheduling in distributed systems: A cloud computing perspective

abstract: Scheduling is essentially a decision-making process that enables resource sharing among a number of activities by determining their execution order on the set of available resources. The emergence of distributed systems brought new challenges on scheduling in computer systems, including clusters, grids, and more recently clouds. On the other hand, the plethora of research makes it hard for both newcomers researchers to understand the relationship among different scheduling problems and strategies proposed in the literature, which hampers the identification of new and relevant research avenues. In this paper we introduce a classification of the scheduling problem in distributed systems by presenting a taxonomy that incorporates recent developments, especially those in cloud computing. We review the scheduling literature to corroborate the taxonomy and analyze the interest in different branches of the proposed taxonomy. Finally, we identify relevant future directions in scheduling for distributed systems.


title: Integrating Blocking and Non-Blocking MPI Primitives with Task-Based Programming Models

abstract: In this paper we present the Task-Aware MPI library (TAMPI) that integrates both blocking and non-blocking MPI primitives with task-based programming models. The TAMPI library leverages two new runtime APIs to improve both programmability and performance of hybrid applications. The first API allows to pause and resume the execution of a task depending on external events. This API is used to improve the interoperability between blocking MPI communication primitives and tasks. When an MPI operation executed inside a task blocks, the task running is paused so that the runtime system can schedule a new task on the core that became idle. Once the blocked MPI operation is completed, the paused task is put again on the runtime system's ready queue, so eventually it will be scheduled again and its execution will be resumed. The second API defers the release of dependencies associated with a task completion until some external events are fulfilled. This API is composed only of two functions, one to bind external events to a running task and another function to notify about the completion of external events previously bound. TAMPI leverages this API to bind non-blocking MPI operations with tasks, deferring the release of their task dependencies until both task execution and all its bound MPI operations are completed. Our experiments reveal that the enhanced features of TAMPI not only simplify the development of hybrid MPI+OpenMP applications that use blocking or non-blocking MPI primitives but they also naturally overlap computation and communication phases, which improves application performance and scalability by removing artificial dependencies across communication tasks.


title: Long-term microstructural evolution of tungsten under heat and neutron loads

abstract: In nuclear fusion reactors, tungsten will be exposed to high neutron loads at high temper atures (>900 C). The evolution and degradation of the mechanical properties under these conditions is uncertain and therefore constitutes a major risk. Here, the microstructural evolution of tungsten under combined heat and neutron loads is studied, using a multi-scale approach incorporating clusters dynamics and a mean-field recrystallization model. The mean-field recrystallization model contains both nucleation in the bulk and at the grain boundaries. The cluster dynamics model includes the incorporation of loops in the dynamics of the dislocation network as a mechanism. The effects of bulk nucleation on the microstructural evolution are explored. The simulations predict a cyclically occuring neutron-induced recrystallization at all studied temperatures. Furthermore, the evolution of the irradiation hardening during neutron-induced recrystallization is assessed from the simulated microstructures.


title: Fractals and the monadic second order theory of one successor

abstract: We show that if $X$ is virtually any classical fractal subset of $\mathbb{R}^n$, then $(\mathbb{R},<,+,X)$ interprets the monadic second-order theory of $(\mathbb{N},+1)$. This result is sharp in the sense that the standard model of the monadic second-order theory of $(\mathbb{N},+1)$ is known to interpret $(\mathbb{R},<,+,X)$ for various classical fractals $X$ including the middle-thirds Cantor set and the Sierpinski carpet. Let $X \subseteq \mathbb{R}^n$ be closed and nonempty. We show that if the $C^k$-smooth points of $X$ are not dense in $X$ for some $k \geq 1$, then $(\mathbb{R},<,+,X)$ interprets the monadic second-order theory of $(\mathbb{N},+1)$. The same conclusion holds if the packing dimension of $X$ is strictly greater than the topological dimension of $X$ and $X$ has no affine points.


title: Towards an Algebraic Network Information Theory: Simultaneous Joint Typicality Decoding

abstract: Consider a receiver in a multi-user network that wishes to decode several messages. Simultaneous joint typicality decoding is one of the most powerful techniques for determining the fundamental limits at which reliable decoding is possible. This technique has historically been used in conjunction with random i.i.d. codebooks to establish achievable rate regions for networks. Recently, it has been shown that, in certain scenarios, nested linear codebooks in conjunction with "single-user" or sequential decoding can yield better achievable rates. For instance, the compute-forward problem examines the scenario of recovering $L \le K$ linear combinations of transmitted codewords over a $K$-user multiple-access channel (MAC), and it is well established that linear codebooks can yield higher rates. Here, we develop bounds for simultaneous joint typicality decoding used in conjunction with nested linear codebooks, and apply them to obtain a larger achievable region for compute-forward over a $K$-user discrete memoryless MAC. The key technical challenge is that competing codeword tuples that are linearly dependent on the true codeword tuple introduce statistical dependencies, which requires careful partitioning of the associated error events.


title: Penalized estimation of flexible hidden Markov models for time series of counts

abstract: Hidden Markov models are versatile tools for modeling sequential observations, where it is assumed that a hidden state process selects which of finitely many distributions generates any given observation. Specifically for time series of counts, the Poisson family often provides a natural choice for the state-dependent distributions, though more flexible distributions such as the negative binomial or distributions with a bounded range can also be used. However, in practice, choosing an adequate class of (parametric) distributions is often anything but straightforward, and an inadequate choice can have severe negative consequences on the model's predictive performance, on state classification, and generally on inference related to the system considered. To address this issue, we propose an effectively nonparametric approach to fitting hidden Markov models to time series of counts, where the state-dependent distributions are estimated in a completely data-driven way without the need to select a distributional family. To avoid overfitting, we add a roughness penalty based on higher-order differences between adjacent count probabilities to the likelihood, which is demonstrated to produce smooth probability mass functions of the state-dependent distributions. The feasibility of the suggested approach is assessed in a simulation experiment, and illustrated in two real-data applications, where we model the distribution of i) major earthquake counts and ii) acceleration counts of an oceanic whitetip shark (Carcharhinus longimanus) over time.


title: Complementary Quality Measures for Noisy Braiding Operations

abstract: Topological quantum computing with non-abelian anyons in a network of one-dimensional chains relies on braiding operations. In real devices, a noisy environment may compromise these braiding operations. In order to assess the failure acquired during braiding with noisy parameters, I define three quality measures. To keep the results as general as possible, I study these quality measures in a model with minimal assumptions that still allows for different kinds of noise.


title: Multi-Parameter Regression Survival Modelling: An Alternative to Proportional Hazards

abstract: It is standard practice for covariates to enter a parametric model through a single distributional parameter of interest, for example, the scale parameter in many standard survival models. Indeed, the well-known proportional hazards model is of this kind. In this paper we discuss a more general approach whereby covariates enter the model through more than one distributional parameter simultaneously (e.g., scale and shape parameters). We refer to this practice as "multi-parameter regression" (MPR) modelling and explore its use in a survival analysis context. We find that multi-parameter regression leads to more flexible models which can offer greater insight into the underlying data generating process. To illustrate the concept, we consider the two-parameter Weibull model which leads to time-dependent hazard ratios, thus relaxing the typical proportional hazards assumption and motivating a new test of proportionality. A novel variable selection strategy is introduced for such multi-parameter regression models. It accounts for the correlation arising between the estimated regression coefficients in two or more linear predictors -- a feature which has not been considered by other authors in similar settings. The methods discussed have been implemented in the mpr package in R.


title: Region Proposal by Guided Anchoring

abstract: Region anchors are the cornerstone of modern object detection techniques. State-of-the-art detectors mostly rely on a dense anchoring scheme, where anchors are sampled uniformly over the spatial domain with a predefined set of scales and aspect ratios. In this paper, we revisit this foundational stage. Our study shows that it can be done much more effectively and efficiently. Specifically, we present an alternative scheme, named Guided Anchoring, which leverages semantic features to guide the anchoring. The proposed method jointly predicts the locations where the center of objects of interest are likely to exist as well as the scales and aspect ratios at different locations. On top of predicted anchor shapes, we mitigate the feature inconsistency with a feature adaption module. We also study the use of high-quality proposals to improve detection performance. The anchoring scheme can be seamlessly integrated into proposal methods and detectors. With Guided Anchoring, we achieve 9.1% higher recall on MS COCO with 90% fewer anchors than the RPN baseline. We also adopt Guided Anchoring in Fast R-CNN, Faster R-CNN and RetinaNet, respectively improving the detection mAP by 2.2%, 2.7% and 1.2%. Code will be available at https://github.com/open-mmlab/mmdetection.


title: FireLedger: A High Throughput Blockchain Consensus Protocol

abstract: Blockchains are distributed secure ledgers to which transactions are issued continuously and each block of transactions is tightly coupled to its predecessors. Permissioned blockchains place special emphasis on transactions throughput. In this paper we present FireLedger, which leverages the iterative nature of blockchains in order to improve their throughput in optimistic execution scenarios. FireLedger trades latency for throughput in the sense that in FireLedger the last f + 1 blocks of each node's blockchain are considered tentative, i.e., they may be rescinded in case one of the last f + 1 blocks proposers was Byzantine. Yet, when optimistic assumptions are met, a new block is decided in each communication step, which consists of a proposer that sends only its proposal and all other participants are sending a single bit each. Our performance study demonstrates that in a single Amazon data-center, FireLedger running on 10 mid-range Amazon nodes obtains a throughput of up to 160K transactions per second for (typical Bitcoin size) 512 bytes transactions. In a 10 nodes Amazon geo-distributed setting with 512 bytes transactions, FireLedger obtains a throughput of 30K tps. Moreover, on higher end Amazon machines, FireLedger obtains $20%-600%$ better throughput than state of the art protocols like HotStuff and BFT-SMaRt, depending on the exact configuration.


title: Nematic twist-bend phase of a bent liquid crystal dimer: field-induced deformations of the helical structure, macroscopic polarization and fast switching speeds

abstract: The twist-bend nematic (Ntb) phase is a recent addition to the nematic (N) phases of liquid crystals (LCs). A net polar order in the Ntb phase under an external electric field was predicted in several recent theoretical studies but yet to be experimentally realized. We investigated the polar nature, dielectric properties, electro-optical switching and optical transmission properties of a bent LC dimer CB7CB. The LC showed a relatively high-temperature nematic (N) phase and a lower-temperature nematic (NX) phase (also called Ntb in literature). A threshold-dependent polarization current response with large polarization values was obtained in the entire mesophase range. The associated switching times were found in sub-millisecond region. This ferroelectric-like polarization resulted from collective reorientation of polar cybotactic clusters. In NX phase, electric field-induced deformation of twisted helical structures also contributed to net polarization. Dielectric measurements confirmed the presence of cybotactic clusters via relaxation processes with large activation energies. Deformation of the NX helical structure under external electric field was corroborated by polarized optical microscopy and optical transmission studies. The field-induced deformations, net polar order and fast switching will contribute towards greater understanding of the NX (or Ntb) phase dynamics. It may also find applications in next-generation electro-optic devices.


title: Multispectral and Hyperspectral Image Fusion by MS/HS Fusion Net

abstract: Hyperspectral imaging can help better understand the characteristics of different materials, compared with traditional image systems. However, only high-resolution multispectral (HrMS) and low-resolution hyperspectral (LrHS) images can generally be captured at video rate in practice. In this paper, we propose a model-based deep learning approach for merging an HrMS and LrHS images to generate a high-resolution hyperspectral (HrHS) image. In specific, we construct a novel MS/HS fusion model which takes the observation models of low-resolution images and the low-rankness knowledge along the spectral mode of HrHS image into consideration. Then we design an iterative algorithm to solve the model by exploiting the proximal gradient method. And then, by unfolding the designed algorithm, we construct a deep network, called MS/HS Fusion Net, with learning the proximal operators and model parameters by convolutional neural networks. Experimental results on simulated and real data substantiate the superiority of our method both visually and quantitatively as compared with state-of-the-art methods along this line of research.


title: An essential singularity of the cotangent of the Coulomb-nuclear phase shift, and a finite limit of the nuclear part of the effective-range function derived at zero energy

abstract: The Coulomb-nuclear phase shift $\delta^{(cs)}_l$, $\cot\delta^{(cs)}_l$ and a finite limit of the nuclear part $\Delta_l(k)$ of the effective-range function (ERF) are derived for an arbitrary orbital momentum $l$ when energy $E\rightarrow0$. It is proved that $\cot\delta^{(cs)}_l$ has an essential singularity at zero energy, but $\Delta_l(k)$ does not. The explicit finite limit of $\Delta_l(0)$ is found. The property of $\Delta_l(k)$ as a meromorphic function makes possible the analytical continuation of a re-normalized scattering amplitude from the physical energy region to a bound state pole. Then the asymptotic normalization coefficients (ANC) can be deduced from experimental phase-shift data and applied to radiative capture processes which are important in nuclear astrophysics for new elements creation. Our results are in agreement with the results published for $S$ wave scattering.


title: Bayesian calibration and sensitivity analysis for a karst aquifer model using active subspaces

abstract: In this article, we perform a parameter study for a recently developed karst hydrological model. The study consists of a high-dimensional Bayesian inverse problem and a global sensitivity analysis. For the first time in karst hydrology, we use the active subspace method to find directions in the parameter space that dominate the Bayesian update from the prior to the posterior distribution in order to effectively reduce the dimension of the problem and for computational efficiency. Additionally, the calculated active subspace can be exploited to construct sensitivity metrics on each of the individual parameters and be used to construct a natural model surrogate. The model consists of 21 parameters to reproduce the hydrological behavior of spring discharge in a karst aquifer located in the Kerschbaum spring recharge area at Waidhofen a.d. Ybbs in Austria. The experimental spatial and time series data for the inference process were collected by the water works in Waidhofen. We show that this case study has implicit low-dimensionality, and we run an adjusted Markov chain Monte Carlo algorithm in a low-dimensional subspace to construct samples of the posterior distribution. The results are visualized and verified by plots of the posterior's push-forward distribution displaying the uncertainty in predicting discharge values due to the experimental noise in the data. Finally, a discussion provides hydrological interpretation of these results for the Kerschbaum area.


title: Near-UV and optical spectroscopy of comets using the ISIS spectrograph on the WHT

abstract: We present an analysis of long-slit cometary spectroscopy using the dual-arm ISIS spectrograph on the 4.2 m WHT. Eleven comets were observed over two nights in 2016 March and we detected the OH (0-0) emission band at 3085 {\AA} in the spectra of five of these comets. Emission bands of the species NH, CN, C$_{3}$, C$_{2}$, NH$_{2}$ and [OI] were also detected. We used Haser modelling to determine molecular production rates and abundance ratios for the observed species. We found that our average abundances relative to OH and CN were generally consistent with those measured in previous studies.


title: Protograph-Based LDPC Code Design for Probabilistic Shaping with On-Off Keying

abstract: This work investigates protograph-based LDPC codes for the AWGN channel with OOK modulation. A non-uniform distribution of the OOK modulation symbols is considered to improve the power efficiency especially for low SNRs. To this end, a specific transmitter architecture based on time sharing is proposed that allows probabilistic shaping of (some) OOK modulation symbols. Tailored protograph-based LDPC code designs outperform standard schemes with uniform signaling and off-the-shelf codes by 1.1 dB for a transmission rate of 0.25 bits/channel use.


title: Spin-wave Confinement and Coupling in Organic-Based Magnetic Nanostructures

abstract: Vanadium tetracyanoethylene (V[TCNE]$_\text{x}$) is an organic-based ferrimagnet that exhibits robust magnetic ordering (T$_\text{C}$ of over 600 K), high quality-factor (high-Q) microwave resonance (Q up to 3,500), and compatibility with a wide variety of substrates and encapsulation technologies. Here, we substantially expand the potential scope and impact of this emerging material by demonstrating the ability to produce engineered nanostructures with tailored magnetic anisotropy that serve as a platform for the exploration of cavity magnonics, revealing strongly coupled quantum confined standing wave modes that can be tuned into and out of resonance with an applied magnetic field. Specifically, time-domain micromagnetic simulations of these nanostructures faithfully reproduce the experimentally measured spectra, including the quasi-uniform mode and higher-order spin-wave (magnon) modes. Finally, when the two dominant magnon modes present in the spectra are brought into resonance by varying the orientation of the in-plane magnetic field, we observe anti-crossing behavior indicating strong coherent coupling between these two magnon modes at room temperature. These results position V[TCNE]$_\text{x}$ as a leading candidate for the development of coherent magnonics, with potential applications ranging from microwave electronics to quantum information.


title: Monitoring induced distributed double-couple sources using Marchenko-based virtual receivers

abstract: We aim to monitor and characterize signals in the subsurface by combining these passive signals with recorded reflection data at the surface of the Earth. To achieve this, we propose a method to create virtual receivers from reflection data using the Marchenko method. By applying homogeneous Green's function retrieval, these virtual receivers are then used to monitor the responses from subsurface sources. We consider monopole point sources with a symmetric source signal, where the full wavefield without artefacts in the subsurface can be obtained. Responses from more complex source mechanisms, such as double-couple sources, can also be used and provide results with comparable quality as the monopole responses. If the source signal is not symmetric in time, our technique that is based on homogeneous Green's function retrieval provides an incomplete signal, with additional artefacts. The duration of these artefacts is limited and they are only present when the source of the signal is located above the virtual receiver. For sources along a fault rupture, this limitation is also present and more severe due to the source activating over a longer period of time. Part of the correct signal is still retrieved, as well as the source location of the signal. These artefacts do not occur in another method which creates virtual sources as well as receivers from reflection data at the surface. This second method can be used to forecast responses to possible future induced seismicity sources (monopoles, double-couple sources and fault ruptures). This method is applied to field data, where similar results to synthetic data are achieved, which shows the potential for the application on real data signals.


title: On the Origin of Otho-Gardenhose Heliospheric Flux

abstract: Parker spiral theory predicts the heliospheric magnetic field (HMF) will have components of opposite polarity radially toward the Sun and tangentially antiparallel to the solar rotation (i.e. in GSE coordinates, with Bx/By<0). This theory explains the average orientation of the HMF very well but does not predict orthogardenhose (OGH) flux with BX/BY>0 that is often observed. We here study the occurrence and structure of OGH flux, as seen in near-Earth space (r=1AU) and by the Helios spacecraft at 0.29<r<1AU, in order to evaluate the contributions to OGH flux generation of the various mechanisms and factors that are not accounted for in Parker spiral theory. We study the loss of OGH flux with increasing averaging timescale between 16 seconds and 100 hours and so determine its spectrum of scale sizes. OGH flux at sunspot minimum is shown to be more common than at sunspot maximum and caused by smaller-scale structure in the HMF, with a mode temporal scale at a fixed point of 10hrs compared to 40hrs. OGH generated by rotating the HMF through the radial direction is also shown to differ in its spectrum of scale sizes from that generated by rotating the HMF through the tangential direction: the former does not contribute to the excess open heliospheric flux at a given r but the latter does. We show that roughly half of the HMF deflection from the ideal Parker spiral giving OGH at Earth occurs at r below 0.3AU. We highlight some questions which can be addressed by the Parker Solar Probe mission which will study the HMF down to r=0.046AU.


title: Gender Differences in Injury Severity Risk of Single-Vehicle Crashes in Virginia: A Nested Logit Analysis of Heterogeneity

abstract: This paper explores gender differences in injury severity risk using a comprehensive crash dataset including the driver, vehicle, environment, and roadway characteristics. For the purpose of this study, only single vehicle crashes that occurred in the Commonwealth of Virginia, collected in a five year period from 2011 to 2015, were used. This study contributes to the literature of crash analysis by incorporating an extensive dataset with normalized attributes. The dataset used for model estimation integrated from two different data sources. These data sources include Virginia Traffic Records Electronic Data System (TREDS) to incorporate crash information and Virginia Base Mapping Program to incorporate roadway characteristics. A two level nested logit model is developed for each gender, in order to relax the Independence of Irrelevant Alternatives (IIA) assumption. Several crash determinants including but not limited to, drivers characteristics, drivers behavior, vehicle condition, weather conditions, roadway geometry and surface conditions were found to have an impact on both male and female drivers. This study is an assessment of influence of different factors on crash severity segmented by gender. A few interesting results include; (i) female fatality risk exceeds male fatality risk when driving under the influence of alcohol, (ii) speeding contributes to lower fatality risk among female drivers compared to their male counterparts, (iii) poor vehicle conditions have no impact on injury risk to female drivers while having an increased injury risk among male drivers, and (iv) while driving in a work zone area increases the risk of property damage crashes, the impact is higher for male drivers compared to female drivers.


title: Tautological relations for stable maps to a target variety

abstract: We define tautological relations for the moduli space of stable maps to a target variety. Using the double ramification cycle formula for target varieties of Janda-Pandharipande-Pixton-Zvonkine, we construct nontrivial tautological relations parallel to Pixton's double ramification cycle relations for the moduli of curves. Examples and applications are discussed.


title: Monomial multiplicities in explicit form

abstract: In this article we give explicit descriptions of the multiplicities of some classes of monomial ideals. For instance, we give a formula for the multiplicities of all codimension 1 monomial ideals, and another formula for the multiplicities of almost complete intersections. We also introduce a new class of ideals that extends the family of monomial complete intersections, and give a formula for their multiplicity, as well as a visual interpretation of this invariant.


title: Spinning Toroidal Brane Cosmology; A Classical and Quantum Survey

abstract: We construct a cosmological model based on a free particle model which is constrained on an embedded toroidal brane, with a general rotation around a specific axis in the bulk space. Some related issues such as the rotation axis of the brane, the presence of gravitomagnetic background and its relation to the general angular velocity of the brane, and its quantum mechanics and related issues such as minimal length and minimal momentum of the quantum model in the $\mathbb{T}^3$ brane are studied. Also, some cosmological features such as the constraint which is imposed upon the toroidal universe by the isotropy and homogeneity conditions, the corresponding Hubble law, and accelerating expansion for the spinning toroidal model without considering a cosmological constant are also studied.


title: Chiral phase structure and sixteen meson states in SU(3) Polyakov linear-sigma model at finite temperature and chemical potential in strong magnetic field

abstract: In characterizing the chiral phase-structure of pseudoscalars ($J^{pc}=0^{-+}$), scalars ($J^{pc}=0^{++}$), vectors ($J^{pc}=1^{--}$) and axial-vectors ($J^{pc}=1^{++}$) meson states and their dependence on temperature, chemical potential, and magnetic fields, we utilize SU($3$) Polyakov linear-sigma model (PLSM) in mean-field approximation. We first determine the chiral (non)strange quark condensates, $\sigma_l$ and $\sigma_s$ and the corresponding deconfinement order parameters, $\phi$ and $\phi^*$, respectively, in thermal and dense (finite chemical potential) medium and finite magnetic field. The temperature and the chemical potential characteristics of nonet meson states normalized to the lowest {\it bosonic} Matsubara frequencies are analyzed. We noticed that all normalized meson masses become temperature independent at different {\it critical} temperatures. We observe that the chiral and deconfinement phase transitions are shifted to lower {\it quasicritical} temperatures with increasing chemical potential and magnetic field. Thus, we conclude that the magnetic field seems to have almost the same effect as that of the chemical potential, especially on accelerating the phase transition, i.e. inverse magnetic catalysis. We also find that increasing chemical potential enhances the mass degeneracy of the various meson masses, while increasing the magnetic field seems to reduce the critical chemical potential, at which the chiral phase transition takes place. Our mass spectrum calculations agree well with the recent PDG compilations and PNJL, lattice QCD calculations, and QMD/UrQMD simulations.


title: New wavelength calibration of the HARPS spectrograph

abstract: (abridged) Even if the HARPS spectrograph has been operational for more than 15 years and it provides among the most precise Doppler measurements, improvements are still possible. One known problem, for instance, is the non-fully regular block-stitching of the CCDs, which introduces, in some cases, one-year period parasitic signals in the measured radial velocity. The aim is to improve the wavelength calibration of HARPS to push further its planet-detection capabilities. The properties of the CCD stitching-induced pixel-size anomalies are determined with LED flat-field frames, and then a physical, gap-corrected map of the CCDs is used for the fitting model of the spectral orders. We also use a new thorium line list, based on much higher-accuracy measurements than the one used up to now. We derive new wavelength solutions for the 15 years of HARPS data, both before and after the 2015 fibre upgrade. We demonstrate that we correct the gap anomalies by computing the wavelength solutions of laser frequency comb exposures, both with and without taking the gap correction into account. By comparing the rms of the most stable stars of the HARPS sample, we show that we globally decrease the radial velocity dispersion of the data, especially for the data acquired after the change of fibres. Finally, the comparative analysis of several individual systems shows that we manage to attenuate the periodogram power at one year in most cases. The analysis of the RVs derived from individual stellar lines also shows that we correct the stitching-induced RV variation. This improved calibration of the HARPS spectrograph allows to go deeper in the search for low-amplitude radial-velocity signals. It will be further improved by combining the thorium calibration spectra with laser frequency comb and Fabry-Perot calibration spectra, and not only for HARPS but notably also for HARPS-N and ESPRESSO.


title: Generalization Studies of Neural Network Models for Cardiac Disease Detection Using Limited Channel ECG

abstract: Acceleration of machine learning research in healthcare is challenged by lack of large annotated and balanced datasets. Furthermore, dealing with measurement inaccuracies and exploiting unsupervised data are considered to be central to improving existing solutions. In particular, a primary objective in predictive modeling is to generalize well to both unseen variations within the observed classes, and unseen classes. In this work, we consider such a challenging problem in machine learning driven diagnosis: detecting a gamut of cardiovascular conditions (e.g. infarction, dysrhythmia etc.) from limited channel ECG measurements. Though deep neural networks have achieved unprecedented success in predictive modeling, they rely solely on discriminative models that can generalize poorly to unseen classes. We argue that unsupervised learning can be utilized to construct effective latent spaces that facilitate better generalization. This work extensively compares the generalization of our proposed approach against a state-of-the-art deep learning solution. Our results show significant improvements in F1-scores.


title: Prestress tuning of negative refraction and wave channeling from flexural sources

abstract: The quest for wave channeling and manipulation has driven a strong research effort on topological and architected materials, capable of propagating localized electromagnetical or mechanical signals. With reference to an elastic structural grid, which elements can sustain both axial and flexural deformations, it is shown that material interfaces can be created with structural properties tuned by prestress states to achieve total reflection, negative refraction, and strongly localized signal channeling. The achievement of a flat lens and topologically localized modes is demonstrated and tunability of the system allows these properties to hold for a broad range of wavelengths. An ingredient to obtain these effects is the use, suggested here and never attempted before, of concentrated pulsating moments. The important aspect of the proposed method is that states of prestress can be easily removed or changed to tune with continuity the propagational characteristics of the medium, so that a new use of vibration channeling and manipulation is envisaged for elastic materials.


title: Robotics Evolution: from Remote Brain to Cloud

abstract: Robotic systems have been evolving since decades and touching almost all aspects of life, either for leisure or critical applications. Most of traditional robotic systems operate in well-defined environments utilizing pre-configured on-board processing units. However, modern and foreseen robotic applications ask for complex processing requirements that exceed the limits of on-board computing power. Cloud computing and the related technologies have high potential to overcome on-board hardware restrictions and can improve the performance efficiency. This research highlights the advancements in robotic systems with focus on cloud robotics as an emerging trend. There exists an extensive amount of effort to leverage the potentials of robotic systems and to handle arising shortcomings. Moreover, there are promising insights for future breed of intelligent, flexible, and autonomous robotic systems in the Internet of Things era.


title: Automatic detection of passable roads after floods in remote sensed and social media data

abstract: This paper addresses the problem of floods classification and floods aftermath detection utilizing both social media and satellite imagery. Automatic detection of disasters such as floods is still a very challenging task. The focus lies on identifying passable routes or roads during floods. Two novel solutions are presented, which were developed for two corresponding tasks at the MediaEval 2018 benchmarking challenge. The tasks are (i) identification of images providing evidence for road passability and (ii) differentiation and detection of passable and non-passable roads in images from two complementary sources of information. For the first challenge, we mainly rely on object and scene-level features extracted through multiple deep models pre-trained on the ImageNet and Places datasets. The object and scene-level features are then combined using early, late and double fusion techniques. To identify whether or not it is possible for a vehicle to pass a road in satellite images, we rely on Convolutional Neural Networks and a transfer learning-based classification approach. The evaluation of the proposed methods are carried out on the large-scale datasets provided for the benchmark competition. The results demonstrate significant improvement in the performance over the recent state-of-art approaches.


title: An Analysis of the Accuracy of the P300 BCI

abstract: The P300 Brain-Computer Interface (BCI) is a well-established communication channel for severely disabled people. The P300 event-related potential is mostly characterized by its amplitude or its area, which correlate with the spelling accuracy of the P300 speller. Here, we introduce a novel approach for estimating the efficiency of this BCI by considering the P300 signal-to-noise ratio (SNR), a parameter that estimates the spatial and temporal noise levels and has a significantly stronger correlation with spelling accuracy. Furthermore, we suggest a Gaussian noise model, which utilizes the P300 event-related potential SNR to predict spelling accuracy under various conditions for LDA-based classification. We demonstrate the utility of this analysis using real data and discuss its potential applications, such as speeding up the process of electrode selection.


title: On Emergence and Complexity of Ergodic Decompositions

abstract: A concept of emergence was recently introduced in the paper [Berger] in order to quantify the richness of possible statistical behaviors of orbits of a given dynamical system. In this paper, we develop this concept and provide several new definitions, results, and examples. We introduce the notion of topological emergence of a dynamical system, which essentially evaluates how big the set of all its ergodic probability measures is. On the other hand, the metric emergence of a particular reference measure (usually Lebesgue) quantifies how non-ergodic this measure is. We prove fundamental properties of these two emergences, relating them with classical concepts such as Kolmogorov's $\epsilon$-entropy of metric spaces and quantization of measures. We also relate the two types of emergences by means of a variational principle. Furthermore, we provide several examples of dynamics with high emergence. First, we show that the topological emergence of some standard classes of hyperbolic dynamical systems is essentially the maximal one allowed by the ambient. Secondly, we construct examples of smooth area-preserving diffeomorphisms that are extremely non-ergodic in the sense that the metric emergence of the Lebesgue measure is essentially maximal. These examples confirm that super-polynomial emergence indeed exists, as conjectured in the paper [Berger]. Finally, we prove that such examples are locally generic among smooth diffeomorphisms.


title: Joint Power Splitting and Relay Selection in Energy-Harvesting Communications for IoT Networks

abstract: In this paper, we consider an energy-harvesting (EH) relay system consisting of a source, a destination, and multiple EH decode-and-forward (DF) relays that can harvest the energy from their received radio signals. A power-splitting ratio is employed at an EH DF relay to characterize a trade-off between an energy used for decoding the source signal received at the relay and the remaining energy harvested for retransmitting the decode outcome. We propose an optimal power splitting and relay selection (OPS-RS) framework and also consider the conventional equal power splitting and relay selection (EPS-RS) for comparison purposes. We derive closed-form expressions of outage probability for both the OPS-RS and EPS-RS schemes and characterize their diversity gains through an asymptotic outage analysis in the high signal-to-noise ratio region. We further examine an extension of our OPS-RS framework to an energy-harvesting battery (EHB) enabled cooperative relay scenario, where the EH relays are equipped with batteries used to store their harvested energies. We propose an amplify-and-forward (AF) based EHB-OPS-RS and a DF based EHB-OPS-RS schemes for AF and DF relay networks, respectively. Numerical results show that the proposed OPS-RS always outperforms the EPS-RS in terms of outage probability. Moreover, the outage probabilities of AF and DF based EHB-OPS-RS schemes are much smaller than that of OPS-RS and EPS-RS methods, demonstrating the benefit of exploiting the batteries in EH relays. Additionally, the DF based EHB-OPS-RS substantially outperforms the AF based EHB-OPS-RS and such an outage advantage becomes more significant, as the number of EH relays increases.


title: Catalog of Energy Patterns for Mobile Applications

abstract: Software engineers make use of design patterns for reasons that range from performance to code comprehensibility. Several design patterns capturing the body of knowledge of best practices have been proposed in the past, namely creational, structural and behavioral patterns. However, with the advent of mobile devices, it becomes a necessity a catalog of design patterns for energy efficiency. In this work, we inspect commits, issues and pull requests of 1027 Android and 756 iOS apps to identify common practices when improving energy efficiency. This analysis yielded a catalog, available online, with 22 design patterns related to improving the energy efficiency of mobile apps. We argue that this catalog might be of relevance to other domains such as Cyber-Physical Systems and Internet of Things. As a side contribution, an analysis of the differences between Android and iOS devices shows that the Android community is more energy-aware.


title: Stability and Controllability results for a Timoshenko system

abstract: In this paper, we study the indirect boundary stability and exact controllability of a one-dimensional Timoshenko system. In the first part of the paper, we consider the Timoshenko system with only one boundary fractional damping. We first show that the system is strongly stable but not uniformly stable. Hence, we look for a polynomial decay rate for smooth initial data. Using frequency domain arguments combined with the multiplier method, we prove that the energy decay rate depends on coefficients appearing in the PDE and on the order of the fractional damping. Moreover, under the equal speed propagation condition, we obtain the optimal polynomial energy decay rate. In the second part of this paper, we study the indirect boundary exact controllability of the Timoshenko system with mixed Dirichlet-Neumann boundary conditions and boundary control. Using non-harmonic analysis, we first establish a weak observability inequality, which depends on the ratio of the waves propagation speeds. Next, using the HUM method, we prove that the system is exactly controllable in appropriate spaces and that the control time can be small.


title: Risk of Cascading Blackouts Given Correlated Component Outages

abstract: Cascading blackouts typically occur when nearly simultaneous outages occur in k out of N components in a power system, triggering subsequent failures that propagate through the network and cause significant load shedding. While large cascades are rare, their impact can be catastrophic, so quantifying their risk is important for grid planning and operation. A common assumption in previous approaches to quantifying such risk is that the $k$ initiating component outages are statistically independent events. However, when triggered by a common exogenous cause, initiating outages may actually be correlated. Here, copula analysis is used to quantify the impact of correlation of initiating outages on the risk of cascading failure. The method is demonstrated on two test cases; a 2383-bus model of the Polish grid under varying load conditions and a synthetic 10,000-bus model based on the geography of the Western US. The large size of the Western US test case required development of new approaches for bounding an estimate of the total number of N-3 blackout-causing contingencies. The results suggest that both risk of cascading failure, and the relative contribution of higher order contingencies, increase as a function of spatial correlation in component failures.


title: On the Finiteness Problem for classes of modular lattices

abstract: The Finiteness Problem is shown to be unsolvable for any sufficiently large class of modular lattices.


title: Almost Interior Points in Ordered Banach Spaces and the Long--Term Behaviour of Strongly Positive Operator Semigroups

abstract: The first part of this article is a brief survey of the properties of so-called almost interior points in ordered Banach spaces. Those vectors can be seen as a generalization of ``functions which are strictly positive almost everywhere'' on $L^p$-spaces and of ``quasi-interior points'' in Banach lattices. In the second part we study the long--term behaviour of strongly positive operator semigroups on ordered Banach spaces; these are semigroups which, in a sense, map every non-zero positive vector to an almost interior point. Using the Jacobs--de Leeuw--Glicksberg decomposition together with the theory presented in the first part of the paper we deduce sufficiency criteria for such semigroups to converge (strongly or in operator norm) as time tends to infinity. This generalises known results for semigroups on Banach lattices as well as on normally ordered Banach spaces with unit.


title: Sclera Force Control in Robot-assisted Eye Surgery: Adaptive Force Control vs. Auditory Feedback

abstract: Surgeon hand tremor limits human capability during microsurgical procedures such as those that treat the eye. In contrast, elimination of hand tremor through the introduction of microsurgical robots diminishes the surgeon's tactile perception of useful and familiar tool-to-sclera forces. While the large mass and inertia of eye surgical robot prevents surgeon microtremor, loss of perception of small scleral forces may put the sclera at risk of injury. In this paper, we have applied and compared two different methods to assure the safety of sclera tissue during robot-assisted eye surgery. In the active control method, an adaptive force control strategy is implemented on the Steady-Hand Eye Robot in order to control the magnitude of scleral forces when they exceed safe boundaries. This autonomous force compensation is then compared to a passive force control method in which the surgeon performs manual adjustments in response to the provided audio feedback proportional to the magnitude of sclera force. A pilot study with three users indicate that the active control method is potentially more efficient.


title: Lower bounds for rainbow Tur\'{a}n numbers of paths and other trees

abstract: For a fixed graph $F$, we would like to determine the maximum number of edges in a properly edge-colored graph on $n$ vertices which does not contain a rainbow copy of $F$, that is, a copy of $F$ all of whose edges receive a different color. This maximum, denoted by $ex^*(n, F)$, is the rainbow Tur\'{a}n number of $F$. We show that $ex^*(n,P_k)\geq \frac{k}{2}n + O(1)$ where $P_k$ is a path on $k\geq 3$ edges, generalizing a result by Maamoun and Meyniel and by Johnston, Palmer and Sarkar. We show similar bounds for brooms on $2^s-1$ edges and diameter $\leq 10$ and a few other caterpillars of small diameter.


title: A Universal Density Matrix Functional from Molecular Orbital-Based Machine Learning: Transferability across Organic Molecules

abstract: We address the degree to which machine learning can be used to accurately and transferably predict post-Hartree-Fock correlation energies. Refined strategies for feature design and selection are presented, and the molecular-orbital-based machine learning (MOB-ML) method is applied to several test systems. Strikingly, for the MP2, CCSD, and CCSD(T) levels of theory, it is shown that the thermally accessible (350 K) potential energy surface for a single water molecule can be described to within 1 millihartree using a model that is trained from only a single reference calculation at a randomized geometry. To explore the breadth of chemical diversity that can be described, MOB-ML is also applied to a new dataset of thermalized (350 K) geometries of 7211 organic models with up to seven heavy atoms. In comparison with the previously reported $\Delta$-ML method, MOB-ML is shown to reach chemical accuracy with three-fold fewer training geometries. Finally, a transferability test in which models trained for seven-heavy-atom systems are used to predict energies for thirteen-heavy-atom systems reveals that MOB-ML reaches chemical accuracy with 36-fold fewer training calculations than $\Delta$-ML (140 versus 5000 training calculations).


title: Low-Latency Gravitational Wave Alerts for Multi-Messenger Astronomy During the Second Advanced LIGO and Virgo Observing Run

abstract: Advanced LIGO's second observing run (O2), conducted from November 30, 2016 to August 25, 2017, combined with Advanced Virgo's first observations in August 2017 witnessed the birth of gravitational-wave multi-messenger astronomy. The first ever gravitational-wave detection from the coalescence of two neutron stars, GW170817, and its gamma-ray counterpart, GRB 170817A, led to an electromagnetic follow-up of the event at an unprecedented scale. Several teams from across the world searched for EM/neutrino counterparts to GW170817, paving the way for the discovery of optical, X-ray, and radio counterparts. In this article, we describe the online identification of gravitational-wave transients and the distribution of gravitational-wave alerts by the LIGO and Virgo collaborations during O2. We also describe the gravitational-wave observables which were sent in the alerts to enable searches for their counterparts. Finally, we give an overview of the online candidate alerts shared with observing partners during O2. Alerts were issued for 14 candidates, six of which have been confirmed as gravitational-wave events associated with the merger of black holes or neutron stars. Eight of the 14 alerts were issued less than an hour after data acquisition.


title: Surrogate-assisted reliability-based design optimization: a survey and a new general framework

abstract: Reliability-based design optimization (RBDO) is an active field of research with an ever increasing number of contributions. Numerous methods have been proposed for the solution of RBDO, a complex problem that combines optimization and reliability analysis. Classical approaches are based on approximation methods and have been classified in review papers. In this paper, we first review classical approaches based on approximation methods such as FORM, and also more recent methods that rely upon surrogate modelling and Monte Carlo simulation. We then propose a general framework for the solution of RBDO problems that includes three independent blocks, namely adaptive surrogate modelling, reliability analysis and optimization. These blocks are non-intrusive with respect to each other and can be plugged independently in the framework. After a discussion on numerical considerations that require attention for the framework to yield robust solutions to various types of problems, the latter is applied to three examples (using two analytical functions and a finite element model). Kriging and support vector machines together with their own active learning schemes are considered in the surrogate model block. In terms of reliability analysis, the proposed framework is illustrated using both crude Monte Carlo and subset simulation. Finally, the covariance-matrix adaptation - evolution scheme (CMA-ES), a global search algorithm, or sequential quadratic programming (SQP), a local gradient-based method, are used in the optimization block. The comparison of the results to benchmark studies show the effectiveness and efficiency of the proposed framework.


title: Dark Photon Dark Matter from a Network of Cosmic Strings

abstract: We study the production of ultra-light dark photons from a network of near-global, Abelian-Higgs cosmic strings. We find that dark photons produced in this way are nonrelativistic today and can make up all of the dark matter for dark photon masses as small as $m_A \sim 10^{-22} \, \mathrm{eV}$.


title: Mechanization of Separation in Generic Extensions

abstract: We mechanize, in the proof assistant Isabelle, a proof of the axiom-scheme of Separation in generic extensions of models of set theory by using the fundamental theorems of forcing. We also formalize the satisfaction of the axioms of Extensionality, Foundation, Union, and Powerset. The axiom of Infinity is likewise treated, under additional assumptions on the ground model. In order to achieve these goals, we extended Paulson's library on constructibility with renaming of variables for internalized formulas, improved results on definitions by recursion on well-founded relations, and sharpened hypotheses in his development of relativization and absoluteness.


title: Uniformly perfect finitely generated simple left orderable groups

abstract: We show that the finitely generated simple left orderable groups $G_{\rho}$ constructed by the first two authors in arXiv:1807.06478 are uniformly perfect - each element in the group can be expressed as a product of three commutators of elements in the group. This implies that the group does not admit any homogeneous quasimorphism. Moreover, any nontrivial action of the group on the circle, which lifts to an action on the real line, admits a fixed point. Most strikingly, it follows that the groups are examples of left orderable monsters, which means that any faithful action on the real line without a global fixed point is globally contracting. This answers Question 4 from the 2018 ICM proceedings article of the third author. (This question has also been answered simultaneously and independently, using completely different methods, by Matte Bon and Triestino in arXiv:1811.12256.) To prove our results, we provide a certain characterisation of elements of the group $G_{\rho}$ which is a useful new tool in the study of these examples.


title: Automated Synthesis of Safe Digital Controllers for Sampled-Data Stochastic Nonlinear Systems

abstract: We present a new method for the automated synthesis of digital controllers with formal safety guarantees for systems with nonlinear dynamics, noisy output measurements, and stochastic disturbances. Our method derives digital controllers such that the corresponding closed-loop system, modeled as a sampled-data stochastic control system, satisfies a safety specification with probability above a given threshold. The proposed synthesis method alternates between two steps: generation of a candidate controller pc, and verification of the candidate. pc is found by maximizing a Monte Carlo estimate of the safety probability, and by using a non-validated ODE solver for simulating the system. Such a candidate is therefore sub-optimal but can be generated very rapidly. To rule out unstable candidate controllers, we prove and utilize Lyapunov's indirect method for instability of sampled-data nonlinear systems. In the subsequent verification step, we use a validated solver based on SMT (Satisfiability Modulo Theories) to compute a numerically and statistically valid confidence interval for the safety probability of pc. If the probability so obtained is not above the threshold, we expand the search space for candidates by increasing the controller degree. We evaluate our technique on three case studies: an artificial pancreas model, a powertrain control model, and a quadruple-tank process.


title: Compactification of the space of Hamiltonian stationary Lagrangian submanifolds with bounded total extrinsic curvature and volume

abstract: For a sequence of immersed connected closed Hamiltonian stationary Lagrangian submaniolds in $\mathbb{C}^{n}$ with uniform bounds on their volumes and the total extrinsic curvatures, we prove that a subsequence converges either to a point or to a Hamiltonian stationary Lagrangian $n$-varifold locally uniformly in $C^{k}$ for any nonnegative integer $k$ away from a finite set of points, and the limit is Hamiltonian stationary in ${\mathbb{C}}^{n}$. We also obtain a theorem on extending Hamiltonian stationary Lagrangian submanifolds $L$ across a compact set $N$ of Hausdorff codimension at least 2 that is locally noncollapsing in volumes matching its Hausdorff dimension, provided the mean curvature of $L$ is in $L^{n}$ and a condition on local volume of $L$ near $N$ is satisfied.


title: Accelerated Flow for Probability Distributions

abstract: This paper presents a methodology and numerical algorithms for constructing accelerated gradient flows on the space of probability distributions. In particular, we extend the recent variational formulation of accelerated gradient methods in (wibisono, et. al. 2016) from vector valued variables to probability distributions. The variational problem is modeled as a mean-field optimal control problem. The maximum principle of optimal control theory is used to derive Hamilton's equations for the optimal gradient flow. The Hamilton's equation are shown to achieve the accelerated form of density transport from any initial probability distribution to a target probability distribution. A quantitative estimate on the asymptotic convergence rate is provided based on a Lyapunov function construction, when the objective functional is displacement convex. Two numerical approximations are presented to implement the Hamilton's equations as a system of $N$ interacting particles. The continuous limit of the Nesterov's algorithm is shown to be a special case with $N=1$. The algorithm is illustrated with numerical examples.


title: Thermodynamics of majority-logic decoding in information erasure

abstract: We investigate the performance of majority-logic decoding in both reversible and finite-time information erasure processes performed on macroscopic bits that contain $N$ microscopic binary units. While we show that for reversible erasure protocols single-unit transformations are more efficient than majority-logic decoding, the latter is found to offer several benefits for finite-time erasure processes: Both the minimal erasure duration for a given erasure and the minimal erasure error for a given erasure duration are reduced, if compared to a single unit. Remarkably, the majority-logic decoding is also more efficient in both the small erasure error and fast erasure region. These benefits are also preserved under the optimal erasure protocol that minimizes the dissipated heat. Our work therefore shows that majority-logic decoding can lift the precision-speed-efficiency trade-off in information erasure processes.


title: Skeletonisation Algorithms with Theoretical Guarantees for Unorganised Point Clouds with High Levels of Noise

abstract: Data Science aims to extract meaningful knowledge from unorganised data. Real datasets usually come in the form of a cloud of points with only pairwise distances. Numerous applications require to visualise an overall shape of a noisy cloud of points sampled from a non-linear object that is more complicated than a union of disjoint clusters. The skeletonisation problem in its hardest form is to find a 1-dimensional skeleton that correctly represents a shape of the cloud. This paper compares several algorithms that solve the above skeletonisation problem for any point cloud and guarantee a successful reconstruction. For example, given a highly noisy point sample of an unknown underlying graph, a reconstructed skeleton should be geometrically close and homotopy equivalent to (has the same number of independent cycles as) the underlying graph. One of these algorithm produces a Homologically Persistent Skeleton (HoPeS) for any cloud without extra parameters. This universal skeleton contains sub-graphs that provably represent the 1-dimensional shape of the cloud at any scale. Other subgraphs of HoPeS reconstruct an unknown graph from its noisy point sample with a correct homotopy type and within a small offset of the sample. The extensive experiments on synthetic and real data reveal for the first time the maximum level of noise that allows successful graph reconstructions.


title: Half-filled Stripes in a Hole-Doped Three-Orbital Spin-Fermion Model for Cuprates

abstract: Using Monte Carlo techniques, we study a three-orbital CuO$_2$ spin-fermion model for copper-based high critical temperature superconductors that captures the charge-transfer properties of these compounds. Our studies reveal the presence of spin order in the parent compound and, more importantly, stripe spin and charge order under hole doping. Due to the $p$-$d$ orbital hybridizati on, the added holes are approximately equally distributed among the two $p$ orbitals of the oxygen atoms and the $d$ orbital of the copper atoms in the unit cell. In rectangular clusters of dimension $16\times 4$ %that break the $\pi/2$ lattice rotational symmetry, {\it half-filled} stripes are observed upon hole doping, namely when $N_h=2n$ holes are introduced in the system then $n$ stripes of length 4 are formed along the short direction. The original antiferromagnetic order observed in the parent compound develops a $\pi-$shift across each stripe and the magnetic structure factor has a peak at wavevector ${\bf k}=(\pi-\delta,\pi)$ with $\delta=2\pi N_h/N=\pi N_h/2L$, where $L=16$. The electronic charge is also modulated and the charge structure factor is maximized at ${\bf k}=(2\delta,0)$. As electrons are removed from the system, intracell orbital nematicity with $\langle n_{p_x}\rangle-\langle n_{p_y}\rangle\ne 0$ dev elops in the oxygen sector, as well as intercell magnetic nematicity with $\langle S^z_{{\bf i},d}(S^z_{{\bf i}+{\bf x},d}-S^z_{{\bf i}+{\ bf y},d})\rangle\ne 0$ in the spin copper sector, in the standard notation. This occurs not only in rectangular but also in square 8$\times$8 lattices. Overall, our results suggest that the essence of the stripe spin and charge distribution experimentally observed in hole-doped cuprates are captured by unbiased Monte Carlo studies of a simple hole-doped charge-transfer insulator CuO$_2$ spin-fermion model.


title: Standard sirens with a running Planck mass

abstract: We consider the effect of a time-varying Planck mass on the propagation of gravitational waves (GWs). A running Planck mass arises naturally in several modified gravity theories, and here we focus on those that carry an additional dark energy field responsible for the late-time accelerated expansion of the universe, yet--like general relativity (GR)--propagate only two GW polarizations, both traveling at the speed of light. Because a time-varying Planck mass affects the amplitude of the GWs and therefore the inferred distance to the source, standard siren measurements of $H_0$ are degenerate with the parameter $c_M$ characterizing the time-varying Planck mass, where $c_M=0$ corresponds to GR with a constant Planck mass. The effect of non-zero $c_M$ will have a noticeable impact on GWs emitted by binary neutron stars (BNSs) at the sensitivities and distances observable by ground-based GW detectors such as advanced LIGO and A+, implying that standard siren measurements can provide joint constraints on $H_0$ and $c_M$. Assuming a $\Lambda$CDM evolution of the universe and taking Planck's measurement of $H_0$ as a prior, we find that GW170817 constrains $c_M = -9^{+21}_{-28}$ ($68.3\%$ credibility). We also discuss forecasts, finding that if we assume $H_0$ is known independently, then 100 BNS events detected by advanced LIGO can constrain $c_M$ to within $\pm0.9$. This is comparable to the current best constraints from cosmology. Similarly, for 100 LIGO A+ BNS detections, it is possible to constrain $c_M$ to $\pm0.5$. When analyzing joint $H_0$ and $c_M$ constraints we find that $\sim 400$ LIGO A+ events are needed to constrain $H_0$ to $1\%$ accuracy. Finally, we discuss the possibility of a nonzero value of $c_M$ biasing standard siren $H_0$ measurements from 100 LIGO A+ detections, and find that $c_M=+1.35$ could bias $H_0$ by 3-4$\sigma$ too low if we incorrectly assume $c_M=0$.


title: Quantifying magic for multi-qubit operations

abstract: The development of a framework for quantifying "non-stabiliserness" of quantum operations is motivated by the magic state model of fault-tolerant quantum computation, and by the need to estimate classical simulation cost for noisy intermediate-scale quantum (NISQ) devices. The robustness of magic was recently proposed as a well-behaved magic monotone for multi-qubit states and quantifies the simulation overhead of circuits composed of Clifford+T gates, or circuits using other gates from the Clifford hierarchy. Here we present a general theory of the "non-stabiliserness" of quantum operations rather than states, which are useful for classical simulation of more general circuits. We introduce two magic monotones, called channel robustness and magic capacity, which are well-defined for general n-qubit channels and treat all stabiliser-preserving CPTP maps as free operations. We present two complementary Monte Carlo-type classical simulation algorithms with sample complexity given by these quantities and provide examples of channels where the complexity of our algorithms is exponentially better than previous known simulators. We present additional techniques that ease the difficulty of calculating our monotones for special classes of channels.


title: Comparing the geometry of the basins of attraction, the speed and the efficiency of several numerical methods

abstract: We use simple equations in order to compare the basins of attraction on the complex plane, corresponding to a large collection of numerical methods, of several order. Two cases are considered, regarding the total number of the roots, which act as numerical attractors. For both cases we use the iterative schemes for performing a thorough and systematic classification of the nodes on the complex plane. The distributions of the required iterations as well as the probability and their correlations with the corresponding basins of convergence are also discussed. Our numerical calculations suggest that most of the iterative schemes provide relatively similar convergence structures on the complex plane. In addition, several aspects of the numerical methods are compared in an attempt to obtain general conclusions regarding their speed and efficiency. Moreover, we try to determine how the complexity of the each case influences the main characteristics of the numerical methods.


title: Quasi-derivations of Lie-Yamaguti algebras

abstract: The concept of derivation for Lie-Yamaguti algebras is generalized in this paper. A quasi-derivation of an LY-algebra is embedded as derivation in a larger LY-algebra. The relationship between quasi-derivations and robustness of Lie-Yamaguti algebras has been studied.


title: Laboratory photo-chemistry of covalently bonded fluorene clusters: observation of an interesting PAH bowl-forming mechanism

abstract: The fullerene C$_{60}$, one of the largest molecules identified in the interstellar medium (ISM), has been proposed to form top-down through the photo-chemical processing of large (more than 60 C-atoms) polycyclic aromatic hydrocarbon (PAH) molecules. In this article, we focus on the opposite process, investigating the possibility that fullerenes form from small PAHs, in which bowl-forming plays a central role. We combine laboratory experiments and quantum chemical calculations to study the formation of larger PAHs from charged fluorene clusters. The experiments show that with visible laser irradiation, the fluorene dimer cation - [C$_{13}$H$_{9}$$-$C$_{13}$H$_{9}$]$^+$ - and the fluorene trimer cation - [C$_{13}$H$_{9}$$-$C$_{13}$H$_{8}$$-$C$_{13}$H$_{9}$]$^+$ - undergo photo-dehydrogenation and photo-isomerization resulting in bowl structured aromatic cluster-ions, C$_{26}$H$_{12}$$^+$ and C$_{39}$H$_{20}$$^+$, respectively. To study the details of this chemical process, we employ quantum chemistry that allows us to determine the structures of the newly formed cluster-ions, to calculate the hydrogen loss dissociation energies, and to derive the underlying reaction pathways. These results demonstrate that smaller PAH clusters (with less than 60 C-atoms) can convert to larger bowled geometries that might act as building blocks for fullerenes, as the bowl-forming mechanism greatly facilitates the conversion from dehydrogenated PAHs to cages. Moreover, the bowl-forming induces a permanent dipole moment that - in principle - allows to search for such species using radio astronomy.


title: High Throughput Computation of Reference Ranges of Biventricular Cardiac Function on the UK Biobank Population Cohort

abstract: The exploitation of large-scale population data has the potential to improve healthcare by discovering and understanding patterns and trends within this data. To enable high throughput analysis of cardiac imaging data automatically, a pipeline should comprise quality monitoring of the input images, segmentation of the cardiac structures, assessment of the segmentation quality, and parsing of cardiac functional indexes. We present a fully automatic, high throughput image parsing workflow for the analysis of cardiac MR images, and test its performance on the UK Biobank (UKB) cardiac dataset. The proposed pipeline is capable of performing end-to-end image processing including: data organisation, image quality assessment, shape model initialisation, segmentation, segmentation quality assessment, and functional parameter computation; all without any user interaction. To the best of our knowledge,this is the first paper tackling the fully automatic 3D analysis of the UKB population study, providing reference ranges for all key cardiovascular functional indexes, from both left and right ventricles of the heart. We tested our workflow on a reference cohort of 800 healthy subjects for which manual delineations, and reference functional indexes exist. Our results show statistically significant agreement between the manually obtained reference indexes, and those automatically computed using our framework.


title: A New Tensioning Method using Deep Reinforcement Learning for Surgical Pattern Cutting

abstract: Surgeons normally need surgical scissors and tissue grippers to cut through a deformable surgical tissue. The cutting accuracy depends on the skills to manipulate these two tools. Such skills are part of basic surgical skills training as in the Fundamentals of Laparoscopic Surgery. The gripper is used to pinch a point on the surgical sheet and pull the tissue to a certain direction to maintain the tension while the scissors cut through a trajectory. As the surgical materials are deformable, it requires a comprehensive tensioning policy to yield appropriate tensioning direction at each step of the cutting process. Automating a tensioning policy for a given cutting trajectory will support not only the human surgeons but also the surgical robots to improve the cutting accuracy and reliability. This paper presents a multiple pinch point approach to modelling an autonomous tensioning planner based on a deep reinforcement learning algorithm. Experiments on a simulator show that the proposed method is superior to existing methods in terms of both performance and robustness.


title: Modified Jaccard Index Analysis and Adaptive Feature Selection for Location Fingerprinting with Limited Computational Complexity

abstract: We propose an approach for fingerprinting-based positioning which reduces the data requirements and computational complexity of the online positioning stage. It is based on a segmentation of the entire region of interest into subregions, identification of candidate subregions during the online-stage, and position estimation using a preselected subset of relevant features. The subregion selection uses a modified Jaccard index which quantifies the similarity between the features observed by the user and those available within the reference fingerprint map. The adaptive feature selection is achieved using an adaptive forward-backward greedy search which determines a subset of features for each subregion, relevant with respect to a given fingerprinting-based positioning method. In an empirical study using signals of opportunity for fingerprinting the proposed subregion and feature selection reduce the processing time during the online-stage by a factor of about 10 while the positioning accuracy does not deteriorate significantly. In fact, in one of the two study cases the 90th percentile of the circular error increased by 7.5% while in the other study case we even found a reduction of the corresponding circular error by 30%.


title: BrailleBand: Blind Support Haptic Wearable Band for Communication using Braille Language

abstract: Visually impaired people are neglected from many modern communication and interaction procedures. Assistive technologies such as text-to-speech and braille displays are the most commonly used means of connecting such visually impaired people with mobile phones and other smart devices. Both these solutions face usability issues, thus this study focused on developing a user friendly wearable solution called the "BrailleBand" with haptic technology while preserving affordability. The "BrailleBand" enables passive reading using the Braille language. Connectivity between the BrailleBand and the smart device (phone) is established using Bluetooth protocol. It consists of six nodes in three bands worn on the arm to map the braille alphabet, which are actuated to give the sense of touch corresponding to the characters. Three mobile applications were developed for training the visually impaired and to integrate existing smart mobile applications such as navigation and short message service (SMS) with the device BrailleBand. The adaptability, usability and efficiency of reading was tested on a sample of blind users which reflected progressive results. Even though, the reading accuracy depends on the time duration between the characters (character gap) an average Character Transfer Rate of 0.4375 characters per second can be achieved with a character gap of 1000 ms.


title: Exploratory analysis of high-resolution power interruption data reveals spatial and temporal heterogeneity in electric grid reliability

abstract: Modern grid monitoring equipment enables utilities to collect detailed records of power interruptions. These data are aggregated to compute publicly reported metrics describing high-level characteristics of grid performance. The current work explores the depth of insights that can be gained from public data, and the implications of losing visibility into heterogeneity in grid performance through aggregation. We present an exploratory analysis examining three years of high-resolution power interruption data collected by archiving information posted in real-time on the public-facing website of a utility in the Western United States. We report on the size, frequency and duration of individual power interruptions, and on spatio-temporal variability in aggregate reliability metrics. Our results show that metrics of grid performance can vary spatially and temporally by orders of magnitude, revealing heterogeneity that is not evidenced in publicly reported metrics. We show that limited access to granular information presents a substantive barrier to conducting detailed policy analysis, and discuss how more widespread data access could help to answer questions that remain unanswered in the literature to date. Given open questions about whether grid performance is adequate to support societal needs, we recommend establishing pathways to make high-resolution power interruption data available to support policy research.


title: Joint Two-Dimensional Resummation in $q_T$ and $0$-Jettiness at NNLL

abstract: We consider Drell-Yan production $pp \to Z/\gamma^* \to \ell^+\ell^-$ with the simultaneous measurement of the $Z$-boson transverse momentum $q_T$ and $0$-jettiness $\mathcal{T}_0$. Since both observables resolve the initial-state QCD radiation, the double-differential cross section in $q_T$ and $\mathcal{T}_0$ contains Sudakov double logarithms of both $q_T/Q$ and $\mathcal{T}_0/Q$, where $Q \sim m_Z$ is the dilepton invariant mass. We simultaneously resum the logarithms in $q_T$ and $\mathcal{T}_0$ to next-to-next-to-leading logarithmic order (NNLL) matched to next-to-leading fixed order (NLO). Our results provide the first genuinely two-dimensional analytic Sudakov resummation for initial-state radiation. Integrating the resummed double-differential spectrum with an appropriate scale choice over either $\mathcal{T}_0$ or $q_T$ recovers the corresponding single-differential resummation for the remaining variable. We discuss in detail the required effective field theory setups and their combination using two-dimensional resummation profile scales. We also introduce a new method to perform the $q_T$ resummation where the underlying resummation is carried out in impact-parameter space, but is consistently turned off depending on the momentum-space target value for $q_T$. Our methods apply at any order and for any color-singlet production process, such that our results can be systematically extended when the relevant perturbative ingredients become available.


title: Type Ibn Supernovae May not all Come from Massive Stars

abstract: Because core-collapse supernovae are the explosions of massive stars, which have relatively short lifetimes, they occur almost exclusively in galaxies with active star formation. On the other hand, the Type Ibn supernova PS1-12sk exploded in an environment much more typical of thermonuclear (Type Ia) supernovae: on the outskirts of the brightest elliptical galaxy in a galaxy cluster. The lack of any obvious star formation at that location presented a challenge to models of Type Ibn supernovae as the explosions of very massive Wolf-Rayet stars. Here we present a supplementary search for star formation at the site of PS1-12sk, now that the supernova has faded, via deep ultraviolet imaging of the host cluster with the Hubble Space Telescope. We do not detect any ultraviolet emission within 1 kpc of the supernova location, which allows us deepen the limit on star formation rate by an order of magnitude compared to the original study on this event. In light of this new limit, we discuss whether the progenitors of Type Ibn supernovae can be massive stars and what reasonable alternatives have been proposed.


title: Phase Transitions and Gravitational Wave Tests of Pseudo-Goldstone Dark Matter in the Softly Broken U(1) Scalar Singlet Model

abstract: We study phase transitions in a softly broken $U(1)$ complex singlet scalar model in which the dark matter is the pseudo-scalar part of a singlet whose direct detection coupling to matter is strongly suppressed. Our aim is to find ways to test this model with the stochastic gravitational wave background from the scalar phase transition. We find that the phase transition which induces vacuum expectation values for both the Higgs boson and the singlet - necessary to provide a realistic dark matter candidate - is always of the second order. If the stochastic gravitational wave background characteristic to a first order phase transition will be discovered by interferometers, the soft breaking of $U(1)$ cannot be the explanation to the suppressed dark matter-baryon coupling, providing a conclusive negative test for this class of singlet models.


title: MUSE sneaks a peek at extreme ram-pressure stripping events - IV. Hydrodynamic and gravitational interactions in the Blue Infalling Group

abstract: We report new wide-field (~ 4 x 4 square arcmin) MUSE observations of the Blue Infalling Group (BIG), a compact group of galaxies located at a projected distance of ~ 150 kpc from the X-Ray centre of the A1367 cluster at z = 0.021. Our MUSE observations map in detail the extended ionized gas, primarily traced by H$\alpha$ emission, in between the members of the group. The gas morphology and its kinematics appear consistent with a tidal origin due to galaxy encounters, as also supported by the disturbed kinematics visible in one of the group members and the presence of tidal dwarf systems. A diffuse tail extending in the direction opposite to the cluster centre is also detected, hinting at a global ram-pressure stripping of the intra-group material as BIG falls inside A1367. Based on the analysis of spatially-resolved emission line maps, we identify multiple ionization mechanisms for the diffuse gas filaments, including in situ photoionization from embedded HII regions and shocks. Combining spatially resolved kinematics and line ratios, we rule out the association of the most massive galaxy, CGCG097-120, with the group as this system appears to be decoupled from the intragroup medium and subject to strong ram pressure as it falls into A1367. Through our new analysis, we conclude that BIG is shaped by pre-processing produced by gravitational interactions in the local group environment combined with ram pressure stripping by the global cluster halo.


title: Collisional unfolding of quantum Darwinism

abstract: We examine the emergence of objectivity via quantum Darwinism through the use of a collision model, i.e. where the dynamics is modeled through sequences of unitary interactions between the system and the individual constituents of the environment, termed "ancillas". By exploiting versatility of this framework, we show that one can transition from a "Darwinistic" to an "encoding" environment by simply tuning their interaction. Furthermore we establish that in order for a setting to exhibit quantum Darwinism we require a mutual decoherence to occur between the system and environmental ancillas, thus showing that system decoherence alone is not sufficient. Finally, we demonstrate that the observation of quantum Darwinism is sensitive to a non-uniform system-environment interaction.


title: Disruption of satellite galaxies in simulated groups and clusters: the roles of accretion time, baryons, and pre-processing

abstract: We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes.


title: The unWISE Catalog: Two Billion Infrared Sources from Five Years of WISE Imaging

abstract: We present the unWISE Catalog, containing the positions and fluxes of roughly two billion objects observed by the Wide-field Infrared Survey Explorer (WISE) over the full sky. The unWISE Catalog has two advantages over the existing WISE catalog (AllWISE): first, it is based on significantly deeper imaging, and second, it features improved modeling of crowded regions. The deeper imaging used in the unWISE Catalog comes from the coaddition of all publicly available 3$-$5 micron WISE imaging, including that from the ongoing NEOWISE-Reactivation mission, thereby increasing the total exposure time by a factor of 5 relative to AllWISE. At these depths, even at high Galactic latitudes many sources are blended with their neighbors; accordingly, the unWISE analysis simultaneously fits thousands of sources to obtain accurate photometry. Our new catalog detects sources at 5-sigma roughly 0.7 magnitudes fainter than the AllWISE catalog and more accurately models millions of faint sources in the Galactic plane, enabling a wealth of Galactic and extragalactic science. In particular, relative to AllWISE, unWISE doubles the number of galaxies detected between redshifts 0 and 1 and triples the number between redshifts 1 and 2, cataloging more than half a billion galaxies over the whole sky.


title: A new view on Auger data and cosmogenic neutrinos in light of different nuclear disintegration and air-shower models

abstract: We study the implications of Ultra-High Energy Cosmic Ray (UHECR) data from the Pierre Auger Observatory for potential accelerator candidates and cosmogenic neutrino fluxes for different combinations of nuclear disintegration and air-shower models. We exploit the most recent spectral and mass composition data (2017) with a new, computationally very efficient simulation code PriNCe. We extend the systematic framework originally developed by the Pierre Auger Collaboration with the cosmological source evolution as an additional free parameter. In this framework, an ensemble of generalized UHECR accelerators is characterized by a universal spectral index (equal for all injection species), a maximal rigidity, and the normalizations for five nuclear element groups. We find that the 2017 data favor a small but constrained contribution of heavy elements (iron) at the source. We demonstrate that the results moderately depend on the nuclear disintegration (PSB, Peanut, or Talys) model, and more strongly on the air-shower (EPOS-LHC, Sibyll-2.3, or QGSjet-II-04) model. Variations of these models result in different source evolutions and spectral indices, limiting the interpretation in terms of a particular class of cosmic accelerators. Better constrained parameters include the maximal rigidity and the mass composition at the source. Hence, the cosmogenic neutrino flux can be robustly predicted, since it originates from interactions with the cosmic infrared background and peaks at $10^8 \, \mathrm{GeV}$. Depending on the source evolution at high redshifts the flux is likely out of reach of future neutrino observatories in most cases, and a minimal cosmogenic neutrino flux cannot be claimed from data without assuming a cosmological distribution of the sources.


title: The Effect of General Relativistic Precession on Tidal Disruption Events from Eccentric Nuclear Disks

abstract: An eccentric nuclear disk consists of stars moving on apsidally-aligned orbits around a central black hole. The secular gravitational torques that dynamically stabilize these disks can also produce tidal disruption events (TDEs) at very high rates in Newtonian gravity. General relativity, however, is known to quench secular torques via rapid apsidal precession. Here we show that for a disk to black hole mass ratio $\gtrsim 10^{-3}$, the system is in the full loss cone regime. The magnitude of the torque per orbital period acting on a stellar orbit means that general relativistic precession does not have a major effect on the dynamics. Thus we find that TDE rates from eccentric nuclear disks are not affected by general relativistic precession. Furthermore, we show that orbital elements between successive TDEs from eccentric nuclear disks are correlated, potentially resulting in unique observational signatures.