title: On the Local Lipschitz Stability of Bayesian Inverse Problems

abstract: In this note we consider the stability of posterior measures occuring in Bayesian inference w.r.t. perturbations of the prior measure and the log-likelihood function. This extends the well-posedness analysis of Bayesian inverse problems. In particular, we prove a general local Lipschitz continuous dependence of the posterior on the prior and the log-likelihood w.r.t. various common distances of probability measures. These include the total variation, Hellinger, and Wasserstein distance and the Kullback-Leibler divergence. We only assume the boundedness of the likelihoods and measure their perturbations in an $L^p$-norm w.r.t. the prior. The obtained stability yields under mild assumptions the well-posedness of Bayesian inverse problems, in particular, a well-posedness w.r.t. the Wasserstein distance. Moreover, our results indicate an increasing sensitivity of Bayesian inference as the posterior becomes more concentrated, e.g., due to more or more accurate data. This confirms and extends previous observations made in the sensitivity analysis of Bayesian inference.


title: Torsion points and isogenies on CM elliptic curves

abstract: Let $\mathcal{O}$ be an order in the imaginary quadratic field $K$. For positive integers $M \mid N$, we determine the least degree of an $\mathcal{O}$-CM point on the modular curve $X(M,N)_{/K(\zeta_M)}$ and also on the modular curve $X(M,N)_{/\mathbb{Q}(\zeta_M)}$: that is, we treat both the case in which the complex multiplication is rationally defined and the case in which we do not assume that the complex multiplication is rationally defined. To prove these results we establish several new theorems on rational cyclic isogenies of CM elliptic curves. In particular, we extend a result of Kwon that determines the set of positive integers $N$ for which there is an $\mathcal{O}$-CM elliptic curve $E$ admitting a cyclic, $\mathbb{Q}(j(E))$-rational $N$-isogeny.


title: Hierarchical Soft Actor-Critic: Adversarial Exploration via Mutual Information Optimization

abstract: We describe a novel extension of soft actor-critics for hierarchical Deep Q-Networks (HDQN) architectures using mutual information metric. The proposed extension provides a suitable framework for encouraging explorations in such hierarchical networks. A natural utilization of this framework is an adversarial setting, where meta-controller and controller play minimax over the mutual information objective but cooperate on maximizing expected rewards.


title: L\'evy flights and hydrodynamic superdiffusion on the Dirac cone of Graphene

abstract: We show that hydrodynamic collision processes of graphene at the neutrality point can be described in terms of a Fokker-Planck equation with fractional derivative, corresponding to a L\'evy flight in momentum space. Thus, electron-electron collisions give rise to frequent small-angle scattering processes that are interrupted by rare large-angle events. The latter give rise to superdiffusive dynamics of collective excitations. We argue that such superdiffusive dynamics is of more general importance to the out-of-equilibrium dynamics of quantum-critical systems.


title: Slicing the IO execution with ReLayTracer

abstract: Analyzing IO performance anomalies is a crucial task in various computing environments, ranging from large-scale cloud applications to desktop applications. However, the IO stack of modern operating systems is complicated, making it hard to understand the performance anomalies with existing tools. Kernel IO executions are frequently interrupted by internal kernel activities, requiring a sophisticated IO profile tool to deal with the noises. Furthermore, complicated interactions of concurrent IO requests cause different sources of tail latencies in kernel IO stack. As a consequence, developers want to know fine-grained latency profile across IO layers, which may differ in each IO requests. To meet the requirements, this paper suggests ReLayTracer, a per-request, per-layer IO profiler. ReLayTracer enables a detailed analysis to identify root causes of IO performance anomalies by providing per-layer latency distributions of each IO request, hardware performance behavior, and time spent by kernel activities such as an interrupt.


title: Replacing the do-calculus with Bayes rule

abstract: The concept of causality has a controversial history. The question of whether it is possible to represent and address causal problems with probability theory, or if fundamentally new mathematics such as the do calculus is required has been hotly debated, e.g. Pearl (2001) states "the building blocks of our scientific and everyday knowledge are elementary facts such as "mud does not cause rain" and "symptoms do not cause disease" and those facts, strangely enough, cannot be expressed in the vocabulary of probability calculus". This has lead to a dichotomy between advocates of causal graphical modeling and the do calculus, and researchers applying Bayesian methods. In this paper we demonstrate that, while it is critical to explicitly model our assumptions on the impact of intervening in a system, provided we do so, estimating causal effects can be done entirely within the standard Bayesian paradigm. The invariance assumptions underlying causal graphical models can be encoded in ordinary Probabilistic graphical models, allowing causal estimation with Bayesian statistics, equivalent to the do calculus. Elucidating the connections between these approaches is a key step toward enabling the insights provided by each to be combined to solve real problems.


title: A Global Solution Method for Decentralized Multi-Area SCUC and Savings Allocation Based on MILP Value Functions

abstract: To address the issue that Lagrangian dual function based algorithms cannot guarantee convergence and global optimality for decentralized multi-area security constrained unit commitment (M-SCUC) problems, a novel decomposition and coordination method using MILP (mixed integer linear programming) value functions is proposed in this paper. Each regional system operator sets the tie-line power injections as variational parameters in its regional SCUC model, and utilizes a finite algorithm to generate a MILP value function, which returns the optimal generation cost for any given interchange scheduling. With the value functions available from all system operators, theoretically, a coordinator is able to derive a globally optimal interchange scheduling. Since power exchanges may alter the financial position of each area considerably from what it would have been via scheduling independently, we then propose a fair savings allocation method using the values functions derived above and the Shapley value in cooperative game theory. Numerical experiments on a two-area 12-bus system and a three-area 457-bus system are carried out. The validness of the value functions based method is verified for the decentralized M-SCUC problems. The outcome of savings allocation is compared with that of the locational marginal cost based method.


title: Danger of using fully homomorphic encryption: A look at Microsoft SEAL

abstract: Fully homomorphic encryption is a promising crypto primitive to encrypt your data while allowing others to compute on the encrypted data. But there are many well-known problems with fully homomorphic encryption such as CCA security and circuit privacy problem. Despite these problems, there are still many companies are currently using or preparing to use fully homomorphic encryption to build data security applications. It seems that the full homomorphic encryption is very close to practicality and these problems can be easily mitigated in implementation. Although the those problems are well known in theory, there is no public discussion of their actual impact on real application. Our research shows that there are many security pitfalls in fully homomorphic encryption from the perspective of practical application. The security problems of a fully homomorphic encryption in a real application is more severe than imagined. In this paper, we will take Microsoft SEAL as an examples to introduce the security pitfalls of fully homomorphic encryption from the perspective of implementation and practical application


title: Weak Geodesics for the deformed Hermitian-Yang-Mills equation

abstract: We study weak geodesics in the space of potentials for the deformed Hermitian-Yang-Mills equation. The geodesic equation can be formulated as a degenerate elliptic equation, allowing us to employ nonlinear Dirichlet duality theory, as developed by Harvey-Lawson. By exploiting the convexity of the level sets of the Lagrangian angle operator in the highest branch, we are able to construct continuous solutions of the associated Dirichlet problem.


title: A family of non-weight modules over the super-Virasoro algebras

abstract: In this paper, we construct a family of non-weight modules over the super-Virasoro algebras. Those modules when regarded as modules of the Ramond algebra and further restricted as modules over the Cartan subalgebra $\mathfrak{h}$ are free of rank $1$, while when regarded as modules of the Neveu-Schwarz algebra and further restricted as modules over the Cartan subalgebra $\mathfrak{H}$ are free of rank $2$. We obtain a sufficient and necessary condition for such modules to be simple. Moreover, we determine the isomorphism classes of these modules. Finally, we show that these modules constitute a complete classification of free $U(\mathfrak{h})$-modules of rank $1$ over the super-Virasoro algebra of Ramond type, and also constitute a complete classification of free $U(\mathfrak{H})$-modules of rank $2$ over the super-Virasoro algebra of Neveu-Schwarz type.


title: Explicit and compact representations for the Green's function and the Solution of Linear Difference Equations with variable coefficients

abstract: Leibniz' combinatorial formula for determinants is modified to establish a condensed and easily handled compact representation for Hessenbergians, referred to here as Leibnizian representation. Alongside, the elements of a fundamental solution set associated with linear difference equations with variable coefficients of order $p$ are explicitly represented by $p$ banded Hessenbergian solutions, built up solely of the variable coefficients. This yields banded Hessenbergian representations for the elements both of the product of companion matrices and of the determinant ratio formula of the one-sided Green's function (Green's function for short). Combining the above results, the elements of the foregoing notions are endowed with compact representations formulated here by Leibnizian and nested sum representations. We show that the elements of the fundamental solution set can be expressed in terms of the first banded Hessenbergian fundamental solution, called principal determinant function. We also show that the Green's function coincides with the principal determinant function, when both functions are restricted to a fairly large domain. These results yield, an explicit and compact representation of the Green's function restriction along with an explicit and compact solution representation of the previously stated type of difference equations in terms of the variable coefficients, the initial conditions and the forcing term. The equivalence of the Green's function solution representation and the well known single determinant solution representation is derived from first principles. Algorithms and automated software are employed to illustrate the main results of this paper.


title: Propagation of the gravitational waves in a cosmological background

abstract: We investigate the propagation of the gravitational waves in a cosmological background. Based on the framework of spatially covariant gravity, we derive the general quadratic action for the gravitational waves. The spatial derivatives of the extrinsic curvature and the parity-violating terms are systematically introduced. Special attention is paid to the propagation speed of the gravitational waves. We find that it is possible to make the two polarization modes propagate in the same speed, which may differ from that of the light, in the presence of parity-violating terms in the action. In particular, we identify a large class of spatially covariant gravity theories with parity violation, in which both the polarization modes propagate in the speed of light. Our results imply that there are more possibilities in the framework of spatially covariant gravity in light of the propagation speed of the gravitational waves.


title: Avoiding Reasoning Shortcuts: Adversarial Evaluation, Training, and Model Development for Multi-Hop QA

abstract: Multi-hop question answering requires a model to connect multiple pieces of evidence scattered in a long context to answer the question. In this paper, we show that in the multi-hop HotpotQA (Yang et al., 2018) dataset, the examples often contain reasoning shortcuts through which models can directly locate the answer by word-matching the question with a sentence in the context. We demonstrate this issue by constructing adversarial documents that create contradicting answers to the shortcut but do not affect the validity of the original answer. The performance of strong baseline models drops significantly on our adversarial evaluation, indicating that they are indeed exploiting the shortcuts rather than performing multi-hop reasoning. After adversarial training, the baseline's performance improves but is still limited on the adversarial evaluation. Hence, we use a control unit that dynamically attends to the question at different reasoning hops to guide the model's multi-hop reasoning. We show that this 2-hop model trained on the regular data is more robust to the adversaries than the baseline model. After adversarial training, this 2-hop model not only achieves improvements over its counterpart trained on regular data, but also outperforms the adversarially-trained 1-hop baseline. We hope that these insights and initial improvements will motivate the development of new models that combine explicit compositional reasoning with adversarial training.


title: Active Generative Adversarial Network for Image Classification

abstract: Sufficient supervised information is crucial for any machine learning models to boost performance. However, labeling data is expensive and sometimes difficult to obtain. Active learning is an approach to acquire annotations for data from a human oracle by selecting informative samples with a high probability to enhance performance. In recent emerging studies, a generative adversarial network (GAN) has been integrated with active learning to generate good candidates to be presented to the oracle. In this paper, we propose a novel model that is able to obtain labels for data in a cheaper manner without the need to query an oracle. In the model, a novel reward for each sample is devised to measure the degree of uncertainty, which is obtained from a classifier trained with existing labeled data. This reward is used to guide a conditional GAN to generate informative samples with a higher probability for a certain label. With extensive evaluations, we have confirmed the effectiveness of the model, showing that the generated samples are capable of improving the classification performance in popular image classification tasks.


title: Pre-Calabi-Yau algebras and double Poisson brackets

abstract: We give an explicit formula showing how the double Poisson algebra introduced in \cite{VdB} appears as a particular part of a pre-Calabi-Yau structure, i.e. cyclically invariant, with respect to the natural inner form, solution of the Maurer-Cartan equation on $A\oplus A^*$. Specific part of this solution is described, which is in one-to-one correspondence with the double Poisson algebra structures. The result holds for any associative algebra $A$ and emphasizes the special role of the fourth component of a pre-Calabi-Yau structure in this respect. As a consequence we have that appropriate pre-Calabi-Yau structures induce a Poisson brackets on representation spaces $({\rm Rep}_n A)^{Gl_n}$ for any associative algebra $A$.


title: Non-Gaussian Formation of Primordial Black Holes: Effects on the Threshold

abstract: Primordial black holes could have been formed in the early universe from sufficiently large cosmological perturbations re-entering the horizon when the Universe is still radiation dominated. These originate from the spectrum of curvature perturbations generated during inflation at small-scales. Because of the non-linear relation between the curvature perturbation $\zeta$ and the overdensity $\delta\rho$, the formation of the primordial black holes is affected by intrinsic non-Gaussianity even though the curvature perturbation is Gaussian. We investigate the impact of this non-Gaussianity on the critical threshold $\delta_c$ which measures the excess of mass of the perturbation, finding a relative change with respect to the value obtained using a linear relation between $\zeta$ and $\delta\rho$, of a few percent. This shows that the value of the critical threshold is rather robust against non-linearities. The same holds also when cosmologically interesting values of local primordial non-Gaussianity are added to the curvature perturbation.


title: A Bayesian Solution to the M-Bias Problem

abstract: It is common practice in using regression type models for inferring causal effects, that inferring the correct causal relationship requires extra covariates are included or ``adjusted for''. Without performing this adjustment erroneous causal effects can be inferred. Given this phenomenon it is common practice to include as many covariates as possible, however such advice comes unstuck in the presence of M-bias. M-Bias is a problem in causal inference where the correct estimation of treatment effects requires that certain variables are not adjusted for i.e. are simply neglected from inclusion in the model. This issue caused a storm of controversy in 2009 when Rubin, Pearl and others disagreed about if it could be problematic to include additional variables in models when inferring causal effects. This paper makes two contributions to this issue. Firstly we provide a Bayesian solution to the M-Bias problem. The solution replicates Pearl's solution, but consistent with Rubin's advice we condition on all variables. Secondly the fact that we are able to offer a solution to this problem in Bayesian terms shows that it is indeed possible to represent causal relationships within the Bayesian paradigm, albeit in an extended space. We make several remarks on the similarities and differences between causal graphical models which implement the do-calculus and probabilistic graphical models which enable Bayesian statistics. We hope this work will stimulate more research on unifying Pearl's causal calculus using causal graphical models with traditional Bayesian statistics and probabilistic graphical models.


title: Rising from Ashes or Dying Flash? Mega Outburst of Small Comet 289P/Blanpain in 2013

abstract: Jupiter-family comet 289P/Blanpain was first discovered in 1819 and was then lost for $\sim200$ years, only to be rediscovered in 2003 as a small, weakly active comet. The comet is associated with the Phoenicids, an otherwise minor meteor shower that produced significant outbursts in 1956 and 2014. The shower points to the existence of significant mass-loss events of P/Blanpain in recent history. P/Blanpain was recovered during an apparent large outburst in 2013 July at an appreciable heliocentric distance of 3.9 au, with brightness increase of 9 mag, making it one of the largest comet outbursts ever observed. Here we present an analysis of archival data taken by several telescopes. We find that the 2013 outburst has produced $\sim10^{8}$ kg of dust, which accounts for a modest fraction ($\sim1\%$) of the mass of P/Blanpain's nucleus as measured in 2004. Based on analysis of long-term lightcurve and modeling of coma morphology, we conclude that the 2013 outburst was most likely driven by the crystallization of amorphous water ice triggered by a spin-up disruption of the nucleus. Dust dynamical model shows that a small fraction of the dust ejecta will reach the Earth in 2036 and 2041, but are only expected to produce minor enhancements to the Phoenicid meteor shower. The 2013 outburst of P/Blanpain, though remarkable for a comet at small sizes, does not necessary imply a catastrophic disruption of the nucleus. The upcoming close encounter of P/Blanpain in 2020 January will provide an opportunity to examine the current state of the comet.


title: Machine-Assisted Map Editing

abstract: Mapping road networks today is labor-intensive. As a result, road maps have poor coverage outside urban centers in many countries. Systems to automatically infer road network graphs from aerial imagery and GPS trajectories have been proposed to improve coverage of road maps. However, because of high error rates, these systems have not been adopted by mapping communities. We propose machine-assisted map editing, where automatic map inference is integrated into existing, human-centric map editing workflows. To realize this, we build Machine-Assisted iD (MAiD), where we extend the web-based OpenStreetMap editor, iD, with machine-assistance functionality. We complement MAiD with a novel approach for inferring road topology from aerial imagery that combines the speed of prior segmentation approaches with the accuracy of prior iterative graph construction methods. We design MAiD to tackle the addition of major, arterial roads in regions where existing maps have poor coverage, and the incremental improvement of coverage in regions where major roads are already mapped. We conduct two user studies and find that, when participants are given a fixed time to map roads, they are able to add as much as 3.5x more roads with MAiD.


title: Computation of intrinsic spin Hall conductivities from first principles using maximally-localized Wannier functions

abstract: We present a method to compute the intrinsic spin Hall conductivity from first principles using an interpolation scheme based on maximally-localized Wannier functions. After obtaining the relevant matrix elements among the ab initio Bloch states calculated on a coarse k-point mesh, we Fourier transform them to find the corresponding matrix elements between Wannier states. We then perform an inverse Fourier transform to interpolate the velocity and spin-current matrix elements onto a dense k-point mesh, and use them to evaluate the spin Hall conductivity as a Brillouin-zone integral. This strategy has a much lower computational cost than a direct ab initio calculation, without sacrificing the accuracy. We demonstrate that the spin Hall conductivities of platinum and doped gallium arsenide, computed with our interpolation scheme as a function of the Fermi energy, are in good agreement with those obtained in previous first-principles studies. We also discuss certain approximations that can be made, in the spirit of the tight-binding method, to simplify the calculation of the velocity and spin-current matrix elements in the Wannier representation.


title: Evidence of orbit-selective electronic kagome lattice with planar flat-band in correlated paramagnetic YCr6Ge6

abstract: Electronic properties of kagome lattice have drawn great attention recently. In associate with flat-band induced by destructive interference and Dirac cone-type dispersion, abundant exotic phenomena have been theoretically discussed. The material realization of electronic kagome lattice is a crucial step towards comprehending kagome physics and achieving novel quantum phases. Here, combining angle-resolved photoemission spectroscopy, transport measurements and first-principle calculations, we expose a planar flat-band in paramagnetic YCr6Ge6 as a typical signature of electronic kagome lattice. We unearth that the planar flat-band arises from the d_(z^2 ) electrons with intra-kagome-plane hopping forbidden by destructive interference. On the other hand, the destructive interference and flatness of the d_(x^2-y^2 ) and d_xy bands are decomposed possibly due to additional in-plane hopping terms, but the Dirac cone-type dispersion is reserved near chemical potential. We explicitly unveil that orbital character plays an essential role to realize electronic kagome lattice in bulk materials with transition metal kagome layers. Paramagnetic YCr6Ge6 provides an opportunity to comprehend intrinsic properties of electronic kagome lattice as well as its interplays with spin orbit coupling and electronic correlation of Cr-3d electrons, and be free from complications induced by strong local moment of ions in kagome planes.


title: Impact of HTTP Cookie Violations in Web Archives

abstract: Certain HTTP Cookies on certain sites can be a source of content bias in archival crawls. Accommodating Cookies at crawl time, but not utilizing them at replay time may cause cookie violations, resulting in defaced composite mementos that never existed on the live web. To address these issues, we propose that crawlers store Cookies with short expiration time and archival replay systems account for values in the Vary header along with URIs.


title: Interfacial Dzyaloshinskii-Moriya interaction and chiral magnetic textures in a ferrimagnetic insulator

abstract: The interfacial Dzyaloshinskii-Moriya interaction (DMI) in multilayers of heavy metal and ferromagnetic metals enables the stabilization of novel chiral spin structures such as skyrmions. Magnetic insulators, on the other hand can exhibit enhanced dynamics and properties such as lower magnetic damping and therefore it is of interest to combine the properties enabled by interfacial DMI with insulating systems. Here, we demonstrate the presence of interfacial DMI in heterostructures that include insulating magnetic layers. We use a bilayer of perpendicularly magnetized insulating thulium iron garnet (TmIG) and the heavy metal platinum, and find a surprisingly strong interfacial DMI that, combined with spin-orbit torque results, in efficient switching. The interfacial origin is confirmed through thickness dependence measurements of the DMI, revealing the characteristic 1/thickness dependence with one order of magnitude longer decay length compared to metallic layers. We combine chiral spin structures and spin-orbit torques for efficient switching and identify skyrmions that allow us to establish the GGG/TmIG interface as the origin of the DMI.


title: Cosmology Inference from Biased Tracers using the EFT-based Likelihood

abstract: The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-order $n$-point functions of biased tracers, as well as their covariance, and comparing with measurements. Recently, we presented an EFT-based expression for the conditional probability of the density field of a biased tracer given the matter density field, which in principle combines information from arbitrarily high order $n$-point functions. Here, we report results based on this likelihood by applying it to halo catalogs in real space, specifically on the inference of the power spectrum normalization $\sigma_8$. We include bias terms up to second order as well as the leading higher-derivative term. For a cutoff value of $\Lambda = 0.1 h\,{\rm Mpc}^{-1}$, we recover the ground-truth value of $\sigma_8$ to within 95% CL for different halo samples and redshifts. We discuss possible sources for the remaining systematic bias in $\sigma_8$ as well as future developments.


title: Is vortex stretching the main cause of the turbulent energy cascade?

abstract: In three dimensional turbulence there is on average a cascade of kinetic energy from the largest to the smallest scales of the flow. While the dominant idea is that the cascade occurs through the physical process of vortex stretching, evidence for this is debated. In the framework of the Karman-Howarth equation for the two point turbulent kinetic energy, we derive a new result for the average flux of kinetic energy between two points in the flow that reveals the role of vortex stretching. However, the result shows that vortex stretching is in fact not the main contributor to the average energy cascade; the main contributor is the self-amplification of the strain-rate field. We emphasize the need to correctly distinguish and not conflate the roles of vortex stretching and strain-self amplification in order to correctly understand the physics of the cascade, and also resolve a paradox regarding the differing role of vortex stretching on the mechanisms of the energy cascade and energy dissipation rate. Direct numerical simulations are used to confirm the results, as well as provide further results and insights on vortex stretching and strain-self amplification at different scales in the flow. Interestingly, the results imply that while vortex stretching plays a sub-leading role in the average cascade, it may play a leading order role during large fluctuations of the energy cascade about its average behavior.


title: Modeling Music Modality with a Key-Class Invariant Pitch Chroma CNN

abstract: This paper presents a convolutional neural network (CNN) that uses input from a polyphonic pitch estimation system to predict perceived minor/major modality in music audio. The pitch activation input is structured to allow the first CNN layer to compute two pitch chromas focused on different octaves. The following layers perform harmony analysis across chroma and time scales. Through max pooling across pitch, the CNN becomes invariant with regards to the key class (i.e., key disregarding mode) of the music. A multilayer perceptron combines the modality activation output with spectral features for the final prediction. The study uses a dataset of 203 excerpts rated by around 20 listeners each, a small challenging data size requiring a carefully designed parameter sharing. With an R2 of about 0.71, the system clearly outperforms previous systems as well as individual human listeners. A final ablation study highlights the importance of using pitch activations processed across longer time scales, and using pooling to facilitate invariance with regards to the key class.


title: Interpolating between promotion and the long cycle

abstract: We give a new proof of the cyclic sieving phenomena for promotion on rectangular standard tableaux. This uses an action of the cactus groups in the seminormal bases of the irreducible representations of the Hecke algebras.


title: Cusp transitivity in hyperbolic 3-manifolds

abstract: Let $M$ be a cusped finite-volume hyperbolic three-manifold with isometry group $G$. Then $G$ induces a $k$-transitive action by permutation on the cusps of $M$ for some integer $k\ge 0$. Generically $G$ is trivial and $k=0$, but $k>0$ does occur in special cases. We show examples with $k=1,2,4$. An interesting question concerns the possible number of cusps for a fixed $k$. Our main result provides an answer for $k=2$ by constructing a family of manifolds having no upper bound on the number of cusps.


title: CheckNet: Secure Inference on Untrusted Devices

abstract: We introduce CheckNet, a method for secure inference with deep neural networks on untrusted devices. CheckNet is like a checksum for neural network inference: it verifies the integrity of the inference computation performed by untrusted devices to 1) ensure the inference has actually been performed, and 2) ensure the inference has not been manipulated by an attacker. CheckNet is completely transparent to the third party running the computation, applicable to all types of neural networks, does not require specialized hardware, adds little overhead, and has negligible impact on model performance. CheckNet can be configured to provide different levels of security depending on application needs and compute/communication budgets. We present both empirical and theoretical validation of CheckNet on multiple popular deep neural network models, showing excellent attack detection (0.88-0.99 AUC) and attack success bounds.


title: Formation of Incommensurate Charge Density Waves in Cuprates

abstract: Although charge density waves (CDWs) are omnipresent in cuprate high-temperature superconductors, they occur at significantly different wavevectors, confounding efforts to understand their formation mechanism. Here, we use resonant inelastic x-ray scattering to investigate the doping- and temperature-dependent CDW evolution in La2-xBaxCuO4 (x=0.115-0.155). We discovered that the CDW develops in two stages with decreasing temperature. A precursor CDW with quasi-commensurate wavevector emerges first at high-temperature. This doping-independent precursor CDW correlation originates from the CDW phase mode coupled with a phonon and "seeds" the low-temperature CDW with strongly doping dependent wavevector. Our observation reveals the precursor CDW and its phase mode as the building blocks of the highly intertwined electronic ground state in the cuprates.


title: A high-order compact finite difference scheme and precise integration method based on modified Hopf-Cole transformation for numerical simulation of n-dimensional Burgers' system

abstract: This paper modifies a n-dimensional Hopf-Cole transformation to the n-dimensional Burgers' system. We obtain the n-dimensional heat conduction equation through the modification of the Hopf-Cole transformation. Then the fourth-order precise integration method (PIM) in combination with a spatially global sixth-order compact finite difference (CFD) scheme is presented to solve the equation with high accuracy. Moreover, coupling with the Strang splitting method, the scheme is extended to multi-dimensional (two, three-dimensional) Burgers' system. Numerical results show that the proposed method appreciably improves the computational accuracy compared with the existing numerical method.Moreover, the two-dimensional and three-dimensional examples demonstrate excellent adaptability, and the numerical simulation results also have very high accuracy in medium Reynolds numbers.


title: Evaluation of Mater Bi and Polylactic Acid as materials for biodegradable innovative mini-radiosondes to track small scale fluctuations within clouds

abstract: Turbulence plays an important part in determining the chemical and physical processes, on both the micro- and macro-scales, whereby clouds are formed and behave. However, exactly how these are linked together and how turbulence impacts each of these processes is not yet fully understood. This is partly due to a lack of in-situ small scale fluctuation measurements due to a limitation in the available technology. It is in this context that the radiosondes, for which the material characterisation is presented in this paper, are being developed to generate a Lagrangian set of data which can be used to improve the ever-expanding knowledge of atmospheric processes and, in particular, the understanding of the interaction between turbulence and micro-physical phenomenologies inside clouds (www.complete-h2020network.eu). Specifically, the materials developed for the balloons are discussed in further detail within this paper. Mater Bi and polylactic acid are the two common biodegradable thermoplastics that were used initially to make the balloons. To tailor their properties, the balloons were then coated with carnauba wax blended with either pine resin or SiO_2 nanoparticles. The properties such as hydrophobicity, toughness, elasticity and helium gas permeability are investigated and improved in order to keep the density of the radiosondes as constant as possible for a couple of hours. This will allow them to float inside and outside clouds on an isopycnic surface, to measure various properties such as velocity, temperature, pressure and humidity by means of solid state sensors and to transmit them to receivers on Earth. Tests have been made under a rigorous metrological approach comparing the 6 new materials with two reference balloon materials, latex and mylar. It was found that Mater Bi with the two carnaubua wax coatings is the most suited..


title: Geometric and conventional contribution to superfluid weight in twisted bilayer graphene

abstract: By tuning the angle between graphene layers to specific "magic angles" the lowest energy bands of twisted bilayer graphene (TBLG) can be made flat. The flat nature of the bands favors the formation of collective ground states and, in particular, TBLG has been shown to support superconductivity. When the energy bands participating in the superconductivity are well-isolated, the superfluid weight scales inversely with the effective mass of such bands. For flat-band systems one would therefore conclude that even if superconducting pairing is present most of the signatures of the superconducting state should be absent. This conclusion is at odds with the experimental observations for TBLG. We calculate the superfluid weight for TBLG taking into account both the conventional contribution and the contribution arising from the quantum geometry of the bands. We find that both contributions are larger than one would expect treating the bands as well-isolated, that at the magic angle the geometric contribution is larger than the conventional one, and that for small deviations away from the magic angle the conventional contribution is larger than the geometric one. Our results show that, despite the flatness of the bands the superfluid weight in TBLG is finite and consistent with experimental observations. We also show how the superfluid weight can be tuned by varying the chemical potential and the twist angle opening the possibility to tune the nature of the superconducting transition between the standard BCS transition and the Berezinskii-Kosterlitz-Thouless transition.


title: Adversarial attacks on Copyright Detection Systems

abstract: It is well-known that many machine learning models are susceptible to adversarial attacks, in which an attacker evades a classifier by making small perturbations to inputs. This paper discusses how industrial copyright detection tools, which serve a central role on the web, are susceptible to adversarial attacks. We discuss a range of copyright detection systems, and why they are particularly vulnerable to attacks. These vulnerabilities are especially apparent for neural network based systems. As a proof of concept, we describe a well-known music identification method, and implement this system in the form of a neural net. We then attack this system using simple gradient methods. Adversarial music created this way successfully fools industrial systems, including the AudioTag copyright detector and YouTube's Content ID system. Our goal is to raise awareness of the threats posed by adversarial examples in this space, and to highlight the importance of hardening copyright detection systems to attacks.


title: Machine Learning Software Engineering in Practice: An Industrial Case Study

abstract: SAP is the market leader in enterprise software offering an end-to-end suite of applications and services to enable their customers worldwide to operate their business. Especially, retail customers of SAP deal with millions of sales transactions for their day-to-day business. Transactions are created during retail sales at the point of sale (POS) terminals and then sent to some central servers for validations and other business operations. A considerable proportion of the retail transactions may have inconsistencies due to many technical and human errors. SAP provides an automated process for error detection but still requires a manual process by dedicated employees using workbench software for correction. However, manual corrections of these errors are time-consuming, labor-intensive, and may lead to further errors due to incorrect modifications. This is not only a performance overhead on the customers' business workflow but it also incurs high operational costs. Thus, automated detection and correction of transaction errors are very important regarding their potential business values and the improvement in the business workflow. In this paper, we present an industrial case study where we apply machine learning (ML) to automatically detect transaction errors and propose corrections. We identify and discuss the challenges that we faced during this collaborative research and development project, from three distinct perspectives: Software Engineering, Machine Learning, and industry-academia collaboration. We report on our experience and insights from the project with guidelines for the identified challenges. We believe that our findings and recommendations can help researchers and practitioners embarking into similar endeavors.


title: MMDetection: Open MMLab Detection Toolbox and Benchmark

abstract: We present MMDetection, an object detection toolbox that contains a rich set of object detection and instance segmentation methods as well as related components and modules. The toolbox started from a codebase of MMDet team who won the detection track of COCO Challenge 2018. It gradually evolves into a unified platform that covers many popular detection methods and contemporary modules. It not only includes training and inference codes, but also provides weights for more than 200 network models. We believe this toolbox is by far the most complete detection toolbox. In this paper, we introduce the various features of this toolbox. In addition, we also conduct a benchmarking study on different methods, components, and their hyper-parameters. We wish that the toolbox and benchmark could serve the growing research community by providing a flexible toolkit to reimplement existing methods and develop their own new detectors. Code and models are available at https://github.com/open-mmlab/mmdetection. The project is under active development and we will keep this document updated.


title: Probing the origin of the off-pulse emission from the pulsars B0525+21 and B2045-16

abstract: Pulsars typically exhibit radio emission in the form of narrow pulses originated from confined regions of their magnetospheres. A potential presence of magnetospherically originated emission outside this region, the so-called off-pulse emission, would challenge the existing theories. Detection of significant off-pulse emission has been reported so far from only two pulsars, B0525+21 and B2045-16, at 325 and 610 MHz. However, the nature of this newly uncovered off-pulse emission remains unclear. To probe its origin we conducted very high resolution radio observations of B0525+21 and B2045-16 with the European VLBI Network (EVN) at 1.39 GHz. Whereas the pulsed emission is detected at a level consistent with previous observations, we report absence of any off-pulse emission above $42$ and $96\ \mathrm{\mu Jy\ beam^{-1}}$ (three times the rms noise levels) for B0525+21 and B2045-16, respectively. Our stringent upper limits imply the off-pulse emission to be less than $0.4$ and $0.3\%$ of the period-averaged pulsed flux density, i.e., much fainter than the previously suggested values of $1$-$10\%$. Since the EVN data are most sensitive to extremely compact angular scales, our results suggest a non-magnetospheric origin for the previously reported off-pulse emission. Presence of extended emission that is resolved out on these milliarcsecond scales still remains plausible. In this case, we constrain the emission to arise from structures with sizes of $\sim (0.61$-$19) \times 10^3\ \mathrm{au}$ for B0525+21 and $\sim (0.48$-$8.3) \times 10^3\ \mathrm{au}$ for B2045-16. These constraints might indicate that the two pulsars are accompanied by compact bow-shock pulsar wind nebulae. Future observations probing intermediate angular scales ($\sim 0.1$-$5\ \mathrm{arcsec}$) will help in clarifying the actual origin of the off-pulse emission.


title: The suspended pinch point and SU($2$)$\times$U($1$) gauge theories

abstract: We show that the suspended pinch point can be seen as an elliptically fibered variety with singular fibers of type I$_2$ over codimension-one points of the base and a torsionless Mordell--Weil group of rank one. In the F-theory algorithm, this corresponds to a Lie group $SU(2)\times U(1)$. We also identify the matter content as given by the direct sum of the adjoint representation (with zero U($1$)-charge) and the fundamental representation with U($1$)-charge $\pm 1$. We then study the geometry of an SU($2$)$\times$U($1$)-model given by a compact elliptically fibered variety with the singularities of a suspended pinch point. We describe in detail the crepant resolutions and the network of flops of this geometry. We compute topological invariants including the Euler characteristic and Hodge numbers. We also study the weak coupling limit of this geometry and show that it corresponds to an orientifold theory with an Sp($1$)-stack transverse to the orientifold and two brane-image-branes wrapping the orientifold.


title: A Note on Sequences of Lattices

abstract: We investigate the relation between the convergence of a sequence of lattices and the set-theoretic convergence of their corresponding Voronoi cells sequence. We prove that if a sequence of full rank lattices converges to a full rank lattice, then the closures of the limit infimum and limit supremum of the Voronoi cells converges to the corresponding Voronoi cell. It remains an open question if the converse is also true.


title: vGraph: A Generative Model for Joint Community Detection and Node Representation Learning

abstract: This paper focuses on two fundamental tasks of graph analysis: community detection and node representation learning, which capture the global and local structures of graphs, respectively. In the current literature, these two tasks are usually independently studied while they are actually highly correlated. We propose a probabilistic generative model called vGraph to learn community membership and node representation collaboratively. Specifically, we assume that each node can be represented as a mixture of communities, and each community is defined as a multinomial distribution over nodes. Both the mixing coefficients and the community distribution are parameterized by the low-dimensional representations of the nodes and communities. We designed an effective variational inference algorithm which regularizes the community membership of neighboring nodes to be similar in the latent space. Experimental results on multiple real-world graphs show that vGraph is very effective in both community detection and node representation learning, outperforming many competitive baselines in both tasks. We show that the framework of vGraph is quite flexible and can be easily extended to detect hierarchical communities.


title: Enforcing temporal consistency in Deep Learning segmentation of brain MR images

abstract: Longitudinal analysis has great potential to reveal developmental trajectories and monitor disease progression in medical imaging. This process relies on consistent and robust joint 4D segmentation. Traditional techniques are dependent on the similarity of images over time and the use of subject-specific priors to reduce random variation and improve the robustness and sensitivity of the overall longitudinal analysis. This is however slow and computationally intensive as subject-specific templates need to be rebuilt every time. The focus of this work to accelerate this analysis with the use of deep learning. The proposed approach is based on deep CNNs and incorporates semantic segmentation and provides a longitudinal relationship for the same subject. The proposed approach is based on deep CNNs and incorporates semantic segmentation and provides a longitudinal relationship for the same subject. The state of art using 3D patches as inputs to modified Unet provides results around ${0.91 \pm 0.5}$ Dice and using multi-view atlas in CNNs provide around the same results. In this work, different models are explored, each offers better accuracy and fast results while increasing the segmentation quality. These methods are evaluated on 135 scans from the EADC-ADNI Harmonized Hippocampus Protocol. Proposed CNN based segmentation approaches demonstrate how 2D segmentation using prior slices can provide similar results to 3D segmentation while maintaining good continuity in the 3D dimension and improved speed. Just using 2D modified sagittal slices provide us a better Dice and longitudinal analysis for a given subject. For the ADNI dataset, using the simple UNet CNN technique gives us ${0.84 \pm 0.5}$ and while using modified CNN techniques on the same input yields ${0.89 \pm 0.5}$. Rate of atrophy and RMS error are calculated for several test cases using various methods and analyzed.


title: Steering Eco-Evolutionary Games Dynamics with Manifold Control

abstract: Feedback loops between population dynamics of individuals and their ecological environment are ubiquitously found in nature, and have shown profound effects on the resulting eco-evolutionary dynamics. Incorporating linear environmental feedback law into replicator dynamics of two-player games, recent theoretical studies shed light on understanding the oscillating dynamics of social dilemma. However, detailed effects of more general nonlinear feedback loops in multi-player games, which is more common especially in microbial systems, remain unclear. Here, we focus on ecological public goods games with environmental feedbacks driven by nonlinear selection gradient. Unlike previous models, multiple segments of stable and unstable equilibrium manifolds can emerge from the population dynamical systems. We find that a larger relative asymmetrical feedback speed for group interactions centered on cooperators not only accelerates the convergence of stable manifolds, but also increases the attraction basin of these stable manifolds. Furthermore, our work offers an innovative manifold control approach: by designing appropriate switching control laws, we are able to steer the eco-evolutionary dynamics to any desired population states. Our mathematical framework is an important generalization and complement to coevolutionary game dynamics, and also fills the theoretical gap in guiding the widespread problem of population state control in microbial experiments.


title: MQTTg: An Android Implementation

abstract: The Internet of Things (IoT) age is upon us. As we look to build larger networks with more devices connected to the Internet, the need for lightweight protocols that minimize the use of both energy and computation gain popularity. One such protocol is Message Queue Telemetry Transport (MQTT). Since its introduction in 1999, it has slowly increased in use cases and gained a huge spike in popularity since it was used in the popular messaging application Facebook Messenger. In our previous works, we focused on adding geolocation to MQTT, to help modernize the protocol into the IoT age. In this paper, we build off our previous work on MQTTg and build an IoT Android Application that can pull geolocation information from the Operating System. We then use the geolocation data to create geofences to help further tailor the use cases of MQTTg.


title: Exotic Spin-Dependent Interaction Searches at Indiana University

abstract: The axion is a hypothesized particle appearing in various theories beyond the Standard Model. It is a light spin-0 boson initially postulated to solve the strong CP problem and is also a strong candidate for dark matter. If the axion or an axion-like particle exists, it would mediate a P-odd and T-odd spin-dependent interaction. We describe two experiments under development at Indiana University-Bloomington to search for such an interaction.


title: When Risks and Uncertainties Collide: Mathematical Finance for Arbitrage Markets in a Quantum Mechanical View

abstract: Geometric arbitrage theory reformulates a generic asset model possibly allowing for arbitrage by packaging all asset and their forward dynamics into a stochastic principal fibre bundle, with a connection whose parallel transport encodes discounting and portfolio rebalancing, and whose curvature measures, in this geometric language, the instantaneous arbitrage capability generated by the market itself. The asset and market portfolio dynamics have a quantum mechanical description, which is constructed by quantizing the deterministic version of the stochastic Lagrangian system describing a market allowing for arbitrage. Results, obtained by solving the Schroedinger equation, coincide with those obtained by solving the stochastic Euler Lagrange equations derived by a variational principle and providing therefore consistency.


title: High Speed and High Dynamic Range Video with an Event Camera

abstract: Event cameras are novel sensors that report brightness changes in the form of a stream of asynchronous "events" instead of intensity frames. They offer significant advantages with respect to conventional cameras: high temporal resolution, high dynamic range, and no motion blur. While the stream of events encodes in principle the complete visual signal, the reconstruction of an intensity image from a stream of events is an ill-posed problem in practice. Existing reconstruction approaches are based on hand-crafted priors and strong assumptions about the imaging process as well as the statistics of natural images. In this work we propose to learn to reconstruct intensity images from event streams directly from data instead of relying on any hand-crafted priors. We propose a novel recurrent network to reconstruct videos from a stream of events, and train it on a large amount of simulated event data. During training we propose to use a perceptual loss to encourage reconstructions to follow natural image statistics. We further extend our approach to synthesize color images from color event streams. Our network surpasses state-of-the-art reconstruction methods by a large margin in terms of image quality (> 20%), while comfortably running in real-time. We show that the network is able to synthesize high framerate videos (> 5,000 frames per second) of high-speed phenomena (e.g. a bullet hitting an object) and is able to provide high dynamic range reconstructions in challenging lighting conditions. We also demonstrate the effectiveness of our reconstructions as an intermediate representation for event data. We show that off-the-shelf computer vision algorithms can be applied to our reconstructions for tasks such as object classification and visual-inertial odometry and that this strategy consistently outperforms algorithms that were specifically designed for event data.


title: From harmonic mappings to Ricci flow due to the Bochner technique

abstract: The present paper is devoted to the study a global aspect of the geometry of harmonic mappings and, in particular, infinitesimal harmonic transformations, and represents the application of our results to the theory of Ricci solutions and the Ricci flow. These results will be obtained using the methods of Geometric analysis and, in particular, due to theorems of Yau, Li and Schoen on the connections between the geometry of a complete smooth manifold and the global behavior of its subharmonic functions.


title: Microscope objective for imaging atomic strontium with 0.63 micrometer resolution

abstract: Imaging and manipulating individual atoms with submicrometer separation can be instrumental for quantum simulation of condensed matter Hamiltonians and quantum computation with neutral atoms. Quantum gas microscope experiments in most cases rely on quite costly solutions. Here we present an open-source design of a microscope objective for atomic strontium consisting solely of off-the-shelf lenses that is diffraction-limited for 461${\,}$nm light. A prototype built with a simple stacking design is measured to have a resolution of 0.63(4)${\,\mu}$m, which is in agreement with the predicted value. This performance, together with the near diffraction-limited performance for 532${\,}$nm light makes this design useful for both quantum gas microscopes and optical tweezer experiments with strontium. Our microscope can easily be adapted to experiments with other atomic species such as erbium, ytterbium, and dysprosium, as well as Rydberg experiments with rubidium.


title: Deep-Fill: Deep Learning Based Sinogram Domain Gap Filling in Positron Emission Tomography

abstract: One of the major challenges in design and developing of PET, scanners are the presence of inactive areas between the detector blocks which degrade the image spatial resolution and leads to streaking artifacts especially when we employ analytical image reconstruction. The aim of this study is to assess the feasibility of generating the gap-free PET image using the deep convolutional encoder-decoder in sinogram space. The gap-corrupted sinograms of simulated HRRT scanner, sinograms without gaps as ideal/ground truth and predicted sinograms owing to the implemented our deep-fill method were quantitatively compared. In total, 1293 phantom images divided into three main sets of training 1000, 150 Validation, and 143 test set. The 1000 training images were augmented using affine transformations with various sub-transforms including rotation (rotate randomly), translation to 12000 Image with 6 frequencies of 2 main methods. The Deep-Fill architecture consists of an encoder and a decoder part, and it is composed of convolution operation, max pooling, ReLU activation, concatenation, and up convolution layers. The gap image is going through the network along with all possible paths then the gap-free image was generated by decoder part of the network. The quality of the generated images was quantitatively assessed by different quality metrics in both sinogram space and reconstruction images. We demonstrated that deep learning based approaches applied to inter-detector gap filling can recover the missing data in sinogram with high quantitative accuracy and have the potential to significantly improve the reconstructed image quality and prevent degradation of PET image quantification.


title: The limiting distribution of the hook length of a randomly chosen cell in a random Young diagram

abstract: Let $p(n)$ be the number of all integer partitions of the positive integer $n$ and let $\lambda$ be a partition, selected uniformly at random from among all such $p(n)$ partitions. It is known that each partition $\lambda$ has a unique graphical representation, composed by $n$ non-overlapping cells in the plane called Young diagram. As a second step of our sampling experiment, we select a cell $c$ uniformly at random from among all $n$ cells of the Young diagram of the partition $\lambda$. For large $n$, we study the asymptotic behavior of the hook length $Z_n=Z_n(\lambda,c)$ of the cell $c$ of a random partituion $\lambda$. This two-step sampling procedure suggests a product probability measure, which assigns the probability $1/np(n)$ to each pair $(\lambda,c)$. With respect to this probability measure, we show that the random variable $\pi Z_n/\sqrt{6n}$ converges weakly, as $n\to\infty$, to a random variable whose probability density function equals $6y/\pi^2 (e^y-1)$ if $0<y<\infty$, and zero elsewhere.


title: Helicity Maximization Below the Diffraction Limit

abstract: Optimally-chiral electromagnetic fields with maximized helicity density, recently introduced in [1], enable chirality characterization of optically small nanoparticles. Here, we demonstrate a technique to obtain optimally-chiral nearfields that leads to the maximization of helicity density, under the constraint of constant energy density, beyond the diffraction limit. We show how optimally-chiral illumination induces balanced electric and magnetic dipole moments in an achiral dielectric nanoantenna which leads to generating optimally-chiral scattered and total nearfield. In particular, we explore helicity and energy densities in nearfield of a spherical dielectric nanoantenna illuminated by an optimally-chiral combination of azimuthally and radially polarized beams that generates parallel induced electric and magnetic dipole moments that in turn also generate optimally-chiral scattered field with the same handedness of the incident field. The application of helicity maximization to nearfields results in helicity enhancement at nanoscale which is of great advantage in the detection of nanoscale chiral samples, microscopy, and optical manipulation of chiral nanoparticles.


title: The brightest point in accretion disk and black hole spin: implication to the image of black hole M87*

abstract: We propose the simple new method for extracting the value of the black hole spin from the direct high-resolution image of black hole by using a thin accretion disk model. In this model the observed dark region on the first image of the supermassive black hole in the galaxy M87, obtained by the Event Horizon Telescope, is a silhouette of the black hole event horizon. The outline of this silhouette is the equator of the event horizon sphere. The dark silhouette of the black hole event horizon is placed within the expected position of the black hole shadow, which is not revealed on the first image. We calculated numerically the relation between the observed position of the black hole silhouette and the brightest point in the thin accretion disk, depending on the black hole spin. From this relation we derive the spin of the supermassive black hole M87*, $a=0.75\pm0.15$.


title: Equivariant neural networks and equivarification

abstract: We provide a process to modify a neural network to an equivariant one, which we call equivarification. As an illustration, we build an equivariant neural network for image classification by equivarifying a convolutional neural network.


title: Semigroup properties of solutions of SDEs driven by L{\'e}vy processes with independent coordinates

abstract: We study the stochastic differential equation $dX_t = A(X_{t-}) \, dZ_t$, $ X_0 = x$, where $Z_t = (Z_t^{(1)},\ldots,Z_t^{(d)})^T$ and $Z_t^{(1)}, \ldots, Z_t^{(d)}$ are independent one-dimensional L{\'e}vy processes with characteristic exponents $\psi_1, \ldots, \psi_d$. We assume that each $\psi_i$ satisfies a weak lower scaling condition WLSC($\alpha,0,\underline{C}$), a weak upper scaling condition WUSC($\beta,1,\overline{C}$) (where $0< \alpha \le \beta < 2$) and some additional regularity properties. We consider two mutually exclusive assumptions: either (i) all $\psi_1, \ldots, \psi_d$ are the same and $\alpha, \beta$ are arbitrary, or (ii) not all $\psi_1, \ldots, \psi_d$ are the same and $\alpha > (2/3)\beta$. We also assume that the determinant of $A(x) = (a_{ij}(x))$ is bounded away from zero, and $a_{ij}(x)$ are bounded and Lipschitz continuous. In both cases (i) and (ii) we prove that for any fixed $\gamma \in (0,\alpha) \cap (0,1]$ the semigroup $P_t$ of the process $X$ satisfies $|P_t f(x) - P_t f(y)| \le c t^{-\gamma/\alpha} |x - y|^{\gamma} ||f||_\infty$ for arbitrary bounded Borel function $f$. We also show the existence of a transition density of the process $X$.


title: Sound Reception System by an Acoustic Luneburg Lens

abstract: Applying the directionality of a Luneburg lens, a sound reception system for ships is designed. A two-dimensional acoustic Luneburg lens which is made of 253 acrylic pipes of different radii which control the refractive index is used. It has a large disk shape; a diameter of 180 cm and a height of 20 cm. The 8 microphones were installed around the edges of the lens. The sound reception ability of the lens was tested on the sea using the horn of a university training ship as an emitter and the lens loaded on a small boat as a receiver. Both are moving at 20 ~ 40 km/h at a distance of between 300 ~ 700 m. The lens could focus to recognize the direction of the input sound in the frequency range between 300 ~ 900 Hz.


title: A Secure Contained Testbed for Analyzing IoT Botnets

abstract: Many security issues have come to the fore with the increasingly widespread adoption of Internet-of-Things (IoT) devices. The Mirai attack on Dyn DNS service, in which vulnerable IoT devices such as IP cameras, DVRs and routers were infected and used to propagate large-scale DDoS attacks, is one of the more prominent recent examples. IoT botnets, consisting of hundreds-of-thousands of bots, are currently present ``in-the-wild'' at least and are only expected to grow in the future, with the potential to cause significant network downtimes and financial losses to network companies. We propose, therefore, to build testbeds for evaluating IoT botnets and design suitable mitigation techniques against them. A DETERlab-based IoT botnet testbed is presented in this work. The testbed is built in a secure contained environment and includes ancillary services such as DHCP, DNS as well as botnet infrastructure including CnC and scanListen/loading servers. Developing an IoT botnet testbed presented us with some unique challenges which are different from those encountered in non-IoT botnet testbeds and we highlight them in this paper. Further, we point out the important features of our testbed and illustrate some of its capabilities through experimental results.


title: PolSAR Image Classification based on Polarimetric Scattering Coding and Sparse Support Matrix Machine

abstract: POLSAR image has an advantage over optical image because it can be acquired independently of cloud cover and solar illumination. PolSAR image classification is a hot and valuable topic for the interpretation of POLSAR image. In this paper, a novel POLSAR image classification method is proposed based on polarimetric scattering coding and sparse support matrix machine. First, we transform the original POLSAR data to get a real value matrix by the polarimetric scattering coding, which is called polarimetric scattering matrix and is a sparse matrix. Second, the sparse support matrix machine is used to classify the sparse polarimetric scattering matrix and get the classification map. The combination of these two steps takes full account of the characteristics of POLSAR. The experimental results show that the proposed method can get better results and is an effective classification method.


title: (f)RFCDE: Random Forests for Conditional Density Estimation and Functional Data

abstract: Random forests is a common non-parametric regression technique which performs well for mixed-type unordered data and irrelevant features, while being robust to monotonic variable transformations. Standard random forests, however, do not efficiently handle functional data and runs into a curse-of dimensionality when presented with high-resolution curves and surfaces. Furthermore, in settings with heteroskedasticity or multimodality, a regression point estimate with standard errors do not fully capture the uncertainty in our predictions. A more informative quantity is the conditional density p(y | x) which describes the full extent of the uncertainty in the response y given covariates x. In this paper we show how random forests can be efficiently leveraged for conditional density estimation, functional covariates, and multiple responses without increasing computational complexity. We provide open-source software for all procedures with R and Python versions that call a common C++ library.


title: Throttling numbers for adversaries on connected graphs

abstract: In this paper, we answer two open problems from [Breen et al., Throttling for the game of Cops and Robbers on graphs, Discrete Math., 341 (2018) 2418-2430]. The throttling number $th_c(G)$ of a graph $G$ is the minimum possible value of $k + capt_k(G)$ over all positive integers $k$, where $capt_k(G)$ is the number of rounds needed for $k$ cops to capture the robber on $G$. One of the problems from [Breen et al., 2018] was to determine whether there exists a family of trees $T$ of order $n$ for which $th_c(T)$ is asymptotically equal to $2 \sqrt{n}$. We show that such a family cannot exist by improving the upper bound on $\displaystyle \max_{T} th_c(T)$ for all trees $T$ of order $n$ from $2 \sqrt{n}$ to $\frac{\sqrt{14}}{2} \sqrt{n} + O(1)$. We prove this bound by deriving a more general throttling bound for connected graphs that applies to multiple graph adversaries, including the robber and the gambler. This also improves the best known upper bounds on $th_c(G)$ for chordal graphs and unicyclic graphs $G$, as well as throttling numbers for positive semidefinite (PSD) zero forcing on trees. In addition to the results about cop versus robber, we use our general throttling bound to improve previous upper bounds on throttling numbers for the cop versus gambler game on connected graphs. Another open problem from [Breen et al., 2018] was to obtain a bound on $th_c(G)$ for cactus graphs $G$. We prove an $O(\sqrt{n})$ bound for all cactus graphs $G$ of order $n$. Furthermore, we exhibit a family of trees $T$ of order $n$ that have $th_c(T) > 1.4502 \sqrt{n}$ for all $n$ sufficiently large, improving on the previous lower bound of $\lceil \sqrt{2n}-\frac{1}{2} \rceil + 1$ on $\displaystyle \max_{T} th_c(T)$ for trees $T$ of order $n$.


title: Leavitt path algebras over a poset of fields

abstract: Let $E$ be a finite directed graph, and let $I$ be the poset obtained as the antisymmetrization of its set of vertices with respect to a pre-order $\le$ that satisfies $v\le w$ whenever there exists a directed path from $w$ to $v$. Assuming that $I$ is a tree, we define a poset of fields over $I$ as a family $\mathbf K = \{ K_i :i\in I \}$ of fields $K_i$ such that $K_i\subseteq K_j$ if $j\le i$. We define the concepts of a Leavitt path algebra $L_{\mathbf K} (E)$ and a regular algebra $Q_{\mathbf K}(E)$ over the poset of fields $\mathbf K$, and we show that $Q_{\mathbf K}(E)$ is a hereditary von Neumann regular ring, and that its monoid $\mathcal V (Q_{\mathbf K}(E))$ of isomorphism classes of finitely generated projective modules is canonically isomorphic to the graph monoid $M(E)$ of $E$.


title: Search for triple and quadruple beta decay of $^{150}$Nd

abstract: Triple beta decay of $^{150}$Nd to the ground and excited states of $^{150}$Eu and quadruple beta decay of $^{150}$Nd to the excited states of $^{150}$Gd have been studied using a 400 cm$^3$ low-background HPGe detector and an external source consisting of 3046 g of natural Nd$_2$O$_3$ powder. A half-life limit for the quadruple beta decay to the 0$_1^+$ state of $^{150}$Gd was found to be $T_{1/2} (0\nu+4\nu) > 8.7\times10^{20}$ yr (90\% C.L.). For other (0$\nu$+4$\nu$) transitions to the 2$_1^+$, 3$_1^-$, 4$^+_1$, 2$_2^+$ and 2$^+_3$ excited states limits for the first time have been obtained at the level of $(6.1-9.5)\times10^{20}$ yr (90\% C.L.). We report here also the results for the first search for triple beta decay to the ground and excited final states of $^{150}$Eu. We find no evidence of this decay and set lower limits on the half-life in the range $T_{1/2} (0\nu+3\nu) > (0.04-4.8)\times10^{20}$ yr (90\% C.L.).


title: Learning Execution through Neural Code Fusion

abstract: As the performance of computer systems stagnates due to the end of Moore's Law, there is a need for new models that can understand and optimize the execution of general purpose code. While there is a growing body of work on using Graph Neural Networks (GNNs) to learn representations of source code, these representations do not understand how code dynamically executes. In this work, we propose a new approach to use GNNs to learn fused representations of general source code and its execution. Our approach defines a multi-task GNN over low-level representations of source code and program state (i.e., assembly code and dynamic memory states), converting complex source code constructs and complex data structures into a simpler, more uniform format. We show that this leads to improved performance over similar methods that do not use execution and it opens the door to applying GNN models to new tasks that would not be feasible from static code alone. As an illustration of this, we apply the new model to challenging dynamic tasks (branch prediction and prefetching) from the SPEC CPU benchmark suite, outperforming the state-of-the-art by 26% and 45% respectively. Moreover, we use the learned fused graph embeddings to demonstrate transfer learning with high performance on an indirectly related task (algorithm classification).


title: Are there any challenges in the charmonia production and polarization at the LHC?

abstract: We analyze the inclusive prompt production of $\psi^\prime$, $\chi_c$, $J/\psi$ and $\eta_c$ mesons at the LHC using the $k_T$-factorization approach. Our consideration is based on the off-shell production amplitudes for hard partonic subprocesses, nonrelativistic QCD formalism for the formation of bound states and transverse momentum dependent (or unintegrated) gluon densities in a proton derived from Catani-Ciafaloni-Fiorani-Marchesini equation. The nonperturbative color octet transitions are treated in terms of the multipole radiation theory. We extract the corresponding long-distance matrix elements from a combined fit to transverse momentum distributions measured at various LHC experiments. We make predictions for the polarization parameters $\lambda_\theta$, $\lambda_\phi$, $\lambda_{\theta \phi}$ and the frame-independent parameter $\lambda^*$ and compare them to the available $\psi^\prime$ and $J/\psi$ data. Finally, we present a universal set of parameters that provides a reasonable simultaneous description for the whole body of the LHC data (on the $p_T$ distributions, relative production rates and polarization observables) for the whole charmonium family.


title: Eye Gaze Metrics and Analysis of AOI for Indexing Working Memory towards Predicting ADHD

abstract: ADHD is being recognized as a diagnosis which persists into adulthood impacting economic, occupational, and educational outcomes. There is an increased need to accurately diagnose and recommend interventions for this population. One consideration is the development and implementation of reliable and valid outcome measures which reflect core diagnostic criteria. For example, adults with ADHD have reduced working memory capacity when compared to their peers (Michalek et al., 2014). A reduction in working memory capacity indicates attentional control deficits which align with many symptoms outlined on behavioral checklists used to diagnose ADHD. Using computational methods, such as eye tracking technology, to generate a relationship between ADHD and measures of working memory capacity would be useful to advancing our understanding and treatment of the diagnosis in adults. This chapter will outline a feasibility study in which eye tracking was used to measure eye gaze metrics during a working memory capacity task for adults with and without ADHD and machine learning algorithms were applied to generate a feature set unique to the ADHD diagnosis. The chapter will summarize the purpose, methods, results, and impact of this study.


title: Spectropolarimetry of stripped envelope core collapse supernovae and their progenitors

abstract: Over the past 30 years, spectropolarimetry has proven to be a great tool to probe the 3D geometry of core collapse supernovae (CCSNe). The number of high-quality multi-epoch spectropolarimetric data sets for CCSNe remains quite low, however, due to the challenging nature of such observations. In this Thesis we present and analyse the data of two Type IIb SNe (SN 2008aq - two epochs, SN 2011hs - seven epochs) and one Type Ic-bl (SN 2014ad - seven epochs). The latter is the most complete spectropolarimetric data set for a Type Ic-bl to date. We also provide a preliminary re-analysis of the spectropolarimetric data of SN 1993J. The new datasets demonstrate that the spectropolarimetry of stripped envelope supernovae is complex and varied, and although early studies of type IIb SNe showed similar behaviour to type IIP, it is now clear that this early trend is broken by a growing sample of SNe. Additionally, we identified potential signatures of a jet in SN 2014ad. Furthermore, a toy model to simulate spectropolarimetry from simple geometries (similar to Reilly et al. 2016) was created. A systematic approach to parameter space exploration was developped; this showed that unfortunately these models suffer extreme degeneracies and may not lead to accurate interpretations. Finally, spectropolarimetry can also be used to probe asymmetries caused by fast rotation in the potential progenitors of stripped envelope CCSNe. Although numerous Wolf-Rayet stars of type WN and WC have been observed with spectropolarimetry, no single WO stars have been studied to date. In our data of two Galactic WO stars, we found no line effect and could not sufficiently constrain the rotational velocity to exclude a collapsar scenario. Therefore, the absence of a line effect may not necessarily equate to a low rotational velocity.


title: A Dynamic Game Approach to Strategic Design of Secure and Resilient Infrastructure Network

abstract: Infrastructure networks are vulnerable to both cyber and physical attacks. Building a secure and resilient networked system is essential for providing reliable and dependable services. To this end, we establish a two-player three-stage game framework to capture the dynamics in the infrastructure protection and recovery phases. Specifically, the goal of the infrastructure network designer is to keep the network connected before and after the attack, while the adversary aims to disconnect the network by compromising a set of links. With costs for creating and removing links, the two players aim to maximize their utilities while minimizing the costs. In this paper, we use the concept of subgame perfect equilibrium (SPE) to characterize the optimal strategies of the network defender and attacker. We derive the SPE explicitly in terms of system parameters. We further investigate the resilience planning of the defender and the strategic timing of attack of the adversary. Finally, we use case studies of UAV-enabled communication networks for disaster recovery to corroborate the obtained analytical results.


title: Efficient computation of the cumulative distribution function of a linear mixture of independent random variables

abstract: For a variant of the algorithm in [Pit19] (arXiv:1903.10816) to compute the approximate density or distribution function of a linear mixture of independent random variables known by a finite sample, it is presented a proof of the functional correctness, i.e. the convergence of the computed distribution function towards the true distribution function (given the observations) as the algorithm resolution is increased to infinity. The algorithm (like its predecessor version) bears elements which are closely related to early known methods for numerical inversion of the characteristic function of a probability distribution, however here efficiently computes the complete distribution function. Possible applications are in computing the distribution of the bootstrap estimate in any linear bootstrap method (e.g. in the block bootstrap for the mean as parameter of interest, or residual bootstrap in linear regression with fixed design), or in elementary analysis-of-variance hypothesis testing.


title: Jet fragmentation functions for $Z$-tagged jets

abstract: Recently the LHCb collaboration has measured both longitudinal and transverse momentum distribution of hadrons produced inside $Z$-tagged jets in proton-proton collisions at the Large Hadron Collider. These distributions are commonly referred to as jet fragmentation functions and are characterized by the longitudinal momentum fraction $z_h$ of the jet carried by the hadron and the transverse momentum $j_\perp$ with respect to the jet direction. We derive a QCD formalism within Soft-Collinear Effective Theory to describe these distributions and find that the $z_h$-dependence provides information on standard collinear fragmentation functions, while $j_\perp$-dependence probes transverse momentum dependent (TMD) fragmentation functions. We perform theoretical calculations and compare our results with the LHCb data. We find good agreement for the intermediate $z_h$ region. For $j_\perp$-dependence, we suggest binning in both $z_h$ and $j_\perp$, which would lead to a more direct probing of TMD fragmentation functions.


title: Parafermions in Hierarchical Fractional Quantum Hall States

abstract: Motivated by recent theoretical progress demonstrating the existence of non-Abelian parafermion zero modes in domain walls on interfaces between two dimensional Abelian topological phases of matter, we investigate the properties of gapped interfaces of hierarchical fractional quantum Hall states, in the lowest Landau level, characterized by the Hall conductance $\sigma_{xy}(m,p) = \frac{p}{2mp+1} \frac{e^{2}}{h}$, for integer numbers $(m,p)$ with $m,p \geq 1$. The case $m=1$ corresponds to the experimentally well established sequence of fractional quantum Hall states with $\sigma_{xy} = \frac{1}{3}\frac{e^{2}}{h}, \frac{2}{5}\frac{e^{2}}{h}, \frac{3}{7}\frac{e^{2}}{h}, ...\, $, which has been observed in many two dimensional electron gases. Exploring the mechanism by which the $(m,p+1)$ hierarchical state is generated from the condensation of quasiparticles of the ``parent" state $(m,p)$, we uncover a remarkably rich sequence of parafermions in hierarchical interfaces whose quantum dimension $d_{m,p}$ depends both upon the total quantum dimension $\mathcal{D}_{m,p} = \sqrt{2mp+1}$ of the bulk Abelian phase, as well as on the parity of the ``hierarchy level" $p$, which we associate with the $\mathbb{Z}_2$ stability of Majorana zero modes in one dimensional topological superconductors. We show that these parafermions reside on domain walls separating segments of the interface where the low energy modes are gapped by two distinct mechanisms: (1) a charge neutral backscattering process or (2) an interaction that breaks $U(1)$ charge conservation symmetry and stabilizes a condensate whose charge depends on $p$. Remarkably, this charge condensate corresponds to a clustering of quasiparticles of fractional charge $\frac{p}{2mp+1}\,e$, allowing us to draw a correspondence between these fractionalized condensates and Read-Rezayi non-Abelian fractional quantum Hall cluster states.


title: The Effect of Distant Encounters on Black Hole Binaries in Globular Clusters: Systematic Increase of In-cluster Mergers in the LISA band

abstract: We present the first systematic study on how distant weak interactions impact the dynamical evolution of merging binary black holes (BBHs) in dense stellar clusters. Recent studies indicate that dense clusters are likely to significantly contribute to the rate of merging BBHs observable through gravitational waves (GWs), and that many of these mergers will appear with notable eccentricities measurable in the LISA and LIGO sensitivity bands. This is highly interesting, as eccentricity can be used to distinguish between different astrophysical merger channels. However, all of these recent studies are based on various Monte Carlo (MC) techniques that only include strong interactions for the dynamical evolution of BBHs, whereas any binary generally undergoes orders-of-magnitude more weak interactions than strong. It is well known that weak interactions primarily lead to a change in the binary's eccentricity, which for BBHs implies that weak interactions can change their GW inspiral time and thereby their merger probability. With this motivation, we perform MC simulations of BBHs evolving in dense clusters under the influence of both weak and strong interactions. We find that including weak interactions leads to a notable increase in the number of BBHs that merge inside their cluster, which correspondingly leads to a higher number of eccentric LISA sources. These preliminary results illustrate the importance of including weak interactions for accurately modeling how BBHs merge in clusters, and how to link their emitted GW signals to their astrophysical environment.


title: Extended Thermodynamics and Complexity in Gravitational Chern-Simons Theory

abstract: We study several aspects of the extended thermodynamics of BTZ black holes with thermodynamic mass $M=\alpha m + \gamma \frac{j}{\ell}$ and angular momentum $J = \alpha j + \gamma \ell m$, for general values of the parameters $(\alpha, \gamma)$ ranging from regular ($\alpha=1, \gamma=0$) to exotic ($\alpha=0, \gamma=1$). We show that there exist two distinct behaviours for the black holes, one when $\alpha > \gamma$ ("mostly regular"), and the other when $\gamma < \alpha$ ("mostly exotic"). We find that the Smarr formula holds for all $(\alpha, \gamma)$. We derive the corresponding thermodynamic volumes, which we find to be positive provided $\alpha$ and $\gamma$ satisfy a certain constraint. The dependence of pressure on volume is unremarkable and strictly decreasing when $\alpha > \gamma$, but a maximum volume emerges for large $J\gg T$ when $\gamma > \alpha$; consequently an exotic black hole of a given horizon circumference and temperature can exist in two distinct anti de Sitter backgrounds. We compute the reverse isoperimetric ratio, and study the Gibbs free energy and criticality conditions for each. Finally we investigate the complexity growth of these objects and find that they are all proportional to the complexity of the BTZ black hole. Somewhat surprisingly, purely exotic BTZ black holes have vanishing complexity growth.


title: Interacting Valley Chern Insulator and its Topological Imprint on Moir\'{e} Superconductors

abstract: One salient feature of systems with Moir\'{e} superlattice is that, Chern number of "minibands" originating from each valley of the original graphene Brillouin zone becomes a well-defined quantized number because the miniband from each valley can be isolated from the rest of the spectrum due to the Moir\'{e} potential. Then a Moir\'{e} system with a well-defined valley Chern number can become a nonchiral topological insulator with $U(1) \times Z_3 $ symmetry and a $\mathbb{Z}$ classification at the free fermion level. Here we demonstrate that the strongly interacting nature of the Moir\'{e} system reduces the classification of the valley Chern insulator from $\mathbb{Z}$ to $\mathbb{Z}_3$, and it is topologically equivalent to a bosonic symmetry protected topological state made of local boson operators. We also demonstrate that, even if the system becomes a superconductor when doped away from the valley Chern insulator, the valley Chern insulator still leaves a topological imprint as the localized Majorana fermion zero mode in certain geometric configuration.


title: A new energy upper bound for AdS black holes inspired by free field theory

abstract: We consider the toroidally compactified planar AdS-Schwarzschild solution to 4-dimensional gravity with negative cosmological constant. This has a flat torus conformal boundary metric. We show that if the spatial part of the boundary metric is deformed, keeping it static and the temperature and area fixed, then assuming a static bulk solution exists, its energy is less than that of the AdS-Schwarzschild solution. The proof is non-perturbative in the metric deformation. While we expect the same holds for the free energy for black hole solutions we are so far are not able to prove it. In the context of AdS-CFT this implies a 3-dimensional holographic CFT on a flat spatial torus whose bulk dual is AdS-Schwarzschild has a greater energy than if the spatial geometry is deformed in any way that preserves temperature and area. This work was inspired by previous results in free field theory, where scalars and fermions in 3-dimensions have been shown to energetically disfavour flat space.


title: Orbital Stability and Precession Effects in the Kepler-89 System

abstract: Among the numerous discoveries resulting from the {\it Kepler} mission are a plethora of compact planetary systems that provide deep insights into planet formation theories. The architecture of such compact systems also produces unique opportunities to study orbital dynamics in compact environments and the subsequent evolution of orbital parameters. One of the compact {\it Kepler} systems is Kepler-89; a system for which the radial velocity follow-up observations place strong upper limits on the masses of the planets and their Keplerian orbital elements. The potential for non-circular orbits in this system make it a compelling system to study dynamical constraints on the measured orbital parameters. We present a dynamical analysis of the system that demonstrates the stability of the circular model and shows the eccentric model of the system is not stable. The analysis indicates that planets c and d, although close to the 2:1 secular resonance, do not permanently occupy the 2:1 resonance configuration. We explore regions of orbital parameter space to identify the upper bounds of orbital eccentricity for the planets. We further show how the dynamics in the compact system leads to significant periastron precession of the innermost planets. Finally, we quantify the effect of the periastron precession on the transit times of the planets compared with the cyclic variations expected from transit timing variations.


title: Finding Your Voice: The Linguistic Development of Mental Health Counselors

abstract: Mental health counseling is an enterprise with profound societal importance where conversations play a primary role. In order to acquire the conversational skills needed to face a challenging range of situations, mental health counselors must rely on training and on continued experience with actual clients. However, in the absence of large scale longitudinal studies, the nature and significance of this developmental process remain unclear. For example, prior literature suggests that experience might not translate into consequential changes in counselor behavior. This has led some to even argue that counseling is a profession without expertise. In this work, we develop a computational framework to quantify the extent to which individuals change their linguistic behavior with experience and to study the nature of this evolution. We use our framework to conduct a large longitudinal study of mental health counseling conversations, tracking over 3,400 counselors across their tenure. We reveal that overall, counselors do indeed change their conversational behavior to become more diverse across interactions, developing an individual voice that distinguishes them from other counselors. Furthermore, a finer-grained investigation shows that the rate and nature of this diversification vary across functionally different conversational components.


title: AdS$_4$/CFT$_3$ from Weak to Strong String Coupling

abstract: We consider the four-point function of operators in the stress tensor multiplet of the $U(N)_k\times U(N)_{-k}$ ABJM theory, in the limit where $N$ is taken to infinity while $N/k^{5}$ is held fixed. In this limit, ABJM theory is holographically dual to type IIA string theory on $AdS_4\times \mathbb{CP}^3$ at finite string coupling $g_s \sim (N/k^5)^{1/4}$. While at leading order in $1/N$, the stress tensor multiplet four-point function can be computed from type IIA supergravity, in this work we focus on the first subleading correction, which comes from tree level Witten diagrams with an $R^4$ interaction vertex. Using superconformal Ward identities, bulk locality, and the mass deformed sphere free energy previously computed to all orders in $1/N$ from supersymmetric localization, we determine this $R^4$ correction as a function of $N/k^5$. Taking its flat space limit, we recover the known $R^4$ contribution to the type IIA S-matrix and reproduce the fact that it only receives perturbative contributions in $g_s$ from genus zero and genus one string worldsheets. This is the first check of AdS/CFT at finite $g_s$ for local operators. Our result for the four-point correlator interpolates between the large $N$, large 't Hooft coupling limit and the large $N$ finite $k$ limit. From the bulk perspective, this is an interpolation between type IIA string theory on $AdS_4\times \mathbb{CP}^3$ at small string coupling and M-theory on $AdS_4\times S^7/\mathbb{Z}_k$.


title: The CARMENES search for exoplanets around M dwarfs. Two temperate Earth-mass planet candidates around Teegarden's Star

abstract: Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.


title: Ly$\alpha$ Halos Around $z\sim6$ Quasars

abstract: We present deep MUSE observations of five quasars within the first Gyr of the Universe ($z\gtrsim6$), four of which display extended Ly$\alpha$ halos. After PSF-subtraction, we reveal halos surrounding two quasars for the first time, as well as confirming the presence of two more halos for which tentative detections exist in long-slit spectroscopic observations and narrow-band imaging. The four Ly$\alpha$ halos presented here are diverse in morphology and size, they each display spatial asymmetry, and none are centred on the position of the quasar. Spectra of the diffuse halos demonstrate that none are dramatically offset in velocity from the systemic redshift of the quasars ($\Delta$ v $< 200$ kms$^{-1}$), however each halo shows a broad Ly$\alpha$ line, with a velocity width of order $\sim1000$ kms$^{-1}$. Total Ly$\alpha$ luminosities range between $\sim$ $2 \times 10^{43}$ erg s$^{-1}$ and $\sim$ $2 \times 10^{44}$ erg s$^{-1}$, reaching maximum radial extents of $13 - 30$ pkpc from the quasar positions. We find larger sizes and higher Ly$\alpha$ luminosities than previous literature results at this redshift, but find no correlation between the quasar properties and the Ly$\alpha$ halo, suggesting that the detected emission is most closely related to the physical properties of the circum-galactic medium


title: Renormalons in quantum mechanics

abstract: We present a nonrelativistic one-particle quantum mechanics whose perturbative S-matrix exhibits a renormalon divergence that we explicitely compute. The potential of our model is the sum of the 2d Dirac $\delta$-potential -- known to require renormalization -- and a 1d Dirac $\delta$-potential tilted at an angle. We argue that renormalons are not specific to this example and exist for a much wider class of potentials. The ambiguity in the Borel summation of the perturbative series due to the renormalon pole is resolved by the physical condition of causality through careful consideration of the $i\epsilon$ prescription. The suitably summed perturbative result coincides with the exact answer obtained through the operator formalism for scattering.


title: Notes on 8 Majorana Fermions

abstract: Eight Majorana fermions in $d=1+1$ dimensions enjoy a triality that permutes the representation of the $SO(8)$ global symmetry in which the fermions transform. This triality plays an important role in the quantization of the superstring, and in the analysis of interacting topological insulators and the associated phenomenon of symmetric mass generation. The purpose of these notes is to provide an introduction to the triality and its applications, with careful attention paid to various ${\bf Z}_2$ global and gauge symmetries and their coupling to background spin structures.


title: A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks

abstract: Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modeling, and future directions are discussed.