title: Time-integrated Neutrino Source Searches with 10 years of IceCube Data

abstract: This paper presents the results from point-like neutrino source searches using ten years of IceCube data collected between Apr.~6, 2008 and Jul.~10, 2018. We evaluate the significance of an astrophysical signal from a point-like source looking for an excess of clustered neutrino events with energies typically above $\sim1\,$TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the Northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of $2.9\,\sigma$ after accounting for statistical trials. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the Northern catalog are inconsistent with background at 3.3$\,\sigma$ significance. These results, all based on searches for a cumulative neutrino signal integrated over the ten years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.


title: Privacy-preserving Federated Bayesian Learning of a Generative Model for Imbalanced Classification of Clinical Data

abstract: In clinical research, the lack of events of interest often necessitates imbalanced learning. One approach to resolve this obstacle is data integration or sharing, but due to privacy concerns neither is practical. Therefore, there is an increasing demand for a platform on which an analysis can be performed in a federated environment while maintaining privacy. However, it is quite challenging to develop a federated learning algorithm that can address both privacy-preserving and class imbalanced issues. In this study, we introduce a federated generative model learning platform for generating samples in a data-distributed environment while preserving privacy. We specifically propose approximate Bayesian computation-based Gaussian Mixture Model called 'Federated ABC-GMM', which can oversample data in a minor class by estimating the posterior distribution of model parameters across institutions in a privacy-preserving manner. PhysioNet2012, a dataset for prediction of mortality of patients in an Intensive Care Unit (ICU), was used to verify the performance of the proposed method. Experimental results show that our method boosts classification performance in terms of F1 score up to nearly an ideal situation. It is believed that the proposed method can be a novel alternative to solving class imbalance problems.


title: The approximate variation to pointwise selection principles

abstract: Let $T\subset\mathbb{R}$, $M$ be a metric space with metric $d$, and $M^T$ be the set of all functions mapping $T$ into $M$. Given $f\in M^T$, we study the properties of the approximate variation $\{V_\varepsilon(f)\}_{\varepsilon>0}$, where $V_\varepsilon(f)$ is the greatest lower bound of Jordan variations $V(g)$ of functions $g\in M^T$ such that $d(f(t),g(t))\le\varepsilon$ for all $t\in T$. The notion of $\varepsilon$-variation $V_\varepsilon(f)$ was introduced by Fra\'nkov\'a [Math. Bohem. 116 (1991), 20-59] for intervals $T=[a,b]$ in $\mathbb{R}$ and $M=\mathbb{R}^N$ and extended to the general case by Chistyakov and Chistyakova [Studia Math. 238 (2017), 37-57]. We prove directly the following basic pointwise selection principle: If a sequence of functions $\{f_j\}_{j=1}^\infty$ from $M^T$ is such that the closure in $M$ of the set $\{f_j(t):j\in\mathbb{N}\}$ is compact for all $t\in T$ and $\limsup_{j\to\infty}V_\varepsilon(f_j)$ is finite for all $\varepsilon>0$, then it contains a subsequence, which converges pointwise on $T$ to a bounded regulated function $f\in M^T$. We establish several variants of this result for sequences of regulated and nonregulated functions, for functions with values in reflexive separable Banach spaces, for the almost everywhere convergence and weak pointwise convergence of extracted subsequences, and comment on the necessity of assumptions in the selection principles. The sharpness of all assertions is illustrated by examples.


title: Weakly stationary stochastic processes valued in a separable Hilbert space: Gramian-Cram\'er representations and applications

abstract: The spectral theory for weakly stationary processes valued in a separable Hilbert space has known renewed interest in the past decade. Here we follow earlier approaches which fully exploit the normal Hilbert module property of the time domain. The key point is to build the Gramian-Cram\'er representation as an isomorphic mapping from the modular spectral domain to the modular time domain. We also discuss the general Bochner theorem and provide useful results on the composition and inversion of lag-invariant linear filters. Finally, we derive the Cram\'er-Karhunen-Lo\`eve decomposition and harmonic functional principal component analysis, which are established without relying on additional assumptions.


title: Invariant Gibbs measures and global strong solutions for nonlinear Schr\"odinger equations in dimension two

abstract: We consider the defocusing nonlinear Schr\"odinger equation on $\mathbb{T}^2$ with Wick ordered power nonlinearity, and prove almost sure global well-posedness with respect to the associated Gibbs measure. The heart of the matter is the uniqueness of the solution as limit of solutions to canonically truncated systems. The invariance of the Gibbs measure under the global dynamics follows as a consequence. The proof relies on the novel idea of random averaging operators.


title: Design and operation of the air-cooled beam dump for the extraction line of CERN's Proton Synchrotron Booster (PSB)

abstract: A new beam dump has been designed, built, installed and operated to withstand the future proton beam extracted from the Proton Synchrotron Booster (PSB) in the framework of the LHC Injector Upgrade (LIU) Project at CERN, consisting of up to 1E14 protons per pulse at 2 GeV, foreseen after the machine upgrades planned for CERN's Long Shutdown 2 (2019-2020). In order to be able to efficiently dissipate the heat deposited by the primary beam, the new dump was designed as a cylindrical block assembly, made out of a copper alloy and cooled by forced airflow. In order to determine the energy density distribution deposited by the beam in the dump, Monte Carlo simulations were performed using the FLUKA code, and thermo-mechanical analyses were carried out by importing the energy density into ANSYS. In addition, Computational Fluid Dynamics (CFD) simulations of the airflow were performed in order to accurately estimate the heat transfer convection coefficient on the surface of the dump. This paper describes the design process, highlights the constraints and challenges of integrating a new dump for increased beam power into the existing facility and provides data on the operation of the dump.


title: DLB: Deep Learning Based Load Balancing

abstract: In this paper, we introduce DLB, a Deep Learning based load Balancing mechanism, to effectively address the data skew problem. The key idea of DLB is to replace hash functions in the load balancing mechanisms with deep learning models, which are trained to be able to map different distributions of workloads and data to the servers in a uniform manner. We implemented DLB and deployed it on a practical Cloud environment using CloudSim. Experimental results using both synthetic and real-world data sets show that compared with traditional hash function-based load balancing methods, DLB is able to achieve more balanced mappings, especially when the workload is highly skewed.


title: Logical Performance of 9 Qubit Compass Codes in Ion Traps with Crosstalk Errors

abstract: We simulate four quantum error correcting codes under error models inspired by realistic noise sources in near-term ion trap quantum computers: $T_2$ dephasing, gate overrotation, and crosstalk. We use this data to find preferred codes for given error parameters along with logical error biases and a pseudothreshold which compares the physical and logical gate failure rates for a CNOT gate. Using these results we conclude that Bacon-Shor-13 is the most promising near term candidate as long as the impact of crosstalk can be mitigated through other means.


title: On the supersymmetric solutions of the Heterotic Superstring effective action

abstract: We consider the effective action of the Heterotic Superstring to first order in alpha prime and derive the necessary and sufficient conditions that a field configuration has to satisfy in order to admit at least one Killing spinor using the spinor bilinear method and making minimal coordinate and frame choices. As a previous step in this derivation, we compute the complete spinor bilinear algebra using the Fierz identities, obtaining as a by-product the algebra satisfied by the Spin(7) structure contained in the bilinears in an arbitrary basis. We find the relations existing between the left-hand-sides of the bosonic equations of motion evaluated on supersymmetric field configurations using the Killing Spinor Identities instead of the (far more complicated) integrability conditions of the Killing Spinor Equations as it is common in the literature. We show how to include the Kalb-Ramond's Bianchi identity in the Killing Spinor Identities.


title: A dynamical approach towards Collatz conjecture

abstract: The present work focuses on the study of the renowned Collatz conjecture, also known as the $3x +1$ problem. The distinguished analysis approach lies on the dynamics of an iterative map in binary form. A new estimation of the enlargement of iterated numbers is given. Within the associated iterative map, characteristic periods for periodic orbits are identified.


title: A Benchmark Set of Highly-efficient CUDA and OpenCL Kernels and its Dynamic Autotuning with Kernel Tuning Toolkit

abstract: Autotuning of performance-relevant source-code parameters allows to automatically tune applications without hard coding optimizations and thus helps with keeping the performance portable. In this paper, we introduce a benchmark set of ten autotunable kernels for important computational problems implemented in OpenCL or CUDA. Using our Kernel Tuning Toolkit, we show that with autotuning most of the kernels reach near-peak performance on various GPUs and outperform baseline implementations on CPUs and Xeon Phis. Our evaluation also demonstrates that autotuning is key to performance portability. In addition to offline tuning, we also introduce dynamic autotuning of code optimization parameters during application runtime. With dynamic tuning, the Kernel Tuning Toolkit enables applications to re-tune performance-critical kernels at runtime whenever needed, for example, when input data changes. Although it is generally believed that autotuning spaces tend to be too large to be searched during application runtime, we show that it is not necessarily the case when tuning spaces are designed rationally. Many of our kernels reach near peak-performance with moderately sized tuning spaces that can be searched at runtime with acceptable overhead. Finally we demonstrate, how dynamic performance tuning can be integrated into a real-world application from cryo-electron microscopy domain.


title: Sirius Digital Low-Level RF Status

abstract: Sirius is a Synchrotron Light Source Facility based on a 4th generation low emittance storage ring. The facility is presently in the final assembly phase in Campinas, Brazil, and comprises a 3 GeV electron storage ring, a full energy Booster and a 150 MeV Linac. The Booster RF system is based on a single 5-cell RF cavity driven by a 45 kW amplifier at 500 MHz and is designed to operate at 2 Hz repetition rate. The storage ring will initially use one normal conducting 7-cell RF cavity and, in the final configuration, will comprise two superconducting cavities, each one driven by a 240 kW RF amplifier. A digital LLRF system based on ALBA DLLRF has been designed and commissioned to control both normal conducting and future superconducting RF configurations. The first DLLRF system is already in operation in the injector Booster. This paper presents an overview of the performance of the control loops in high power operation and the results of the Booster DLLRF commissioning. It also shows the integration of DLLRF with the RF system as well as other utilities of the system such as automatic start-up, conditioning, fast interlocks and post-mortem analysis.


title: Strong-field-driven dynamics and high-harmonic generation in interacting 1D systems

abstract: We explore the roles of electronic band structure and Coulomb interactions in solid-state HHG by studying the optical response of linear atomic chains and carbon nanotubes to intense ultrashort pulses. Specifically, we simulate electron dynamics by solving the single-particle density matrix equation of motion in the presence of intense ultrafast optical fields, incorporating tight-binding electronic states and a self-consistent electron-electron interaction. While linear atomic chains constitute an idealized system, our realistic 1D model readily provides insight on the temporal evolution of electronic states in reciprocal space, both in the absence or presence of electron interactions, which we demonstrate to play an important role in the HHG yield. This model further predicts that doped semiconductors generate high harmonics more efficiently than their metallic and undoped counterparts. To complement this idealized system we also show results for HHG in more realistic quasi-1D structures such as carbon nanotubes, the behavior of which is found to be in good qualitative agreement with the atomic chains. Our findings apply directly to extreme nonlinear optical phenomena in atoms on surfaces, carbon-based structures, linear arrays of dopant atoms in semiconductors, and linear molecules, such as polycyclic aromatic hydrocarbon chains, and can be straightforwardly extended to optimize existing platforms for HHG or identify new solid-state alternatives in the context of nonlinear plasmonics.


title: Classification of spherical objects based on the form function of acoustic echoes

abstract: One way to recognise an object is to study how the echo has been shaped during the interaction with the target. Wideband sonar allows the study of the energy distribution for a large range of frequencies. The frequency distribution contains information about an object, including its inner structure. This information is a key for automatic recognition. The scattering by a target can be quantitatively described by its Form Function. The Form Function can be calculated based on the data of the initial pulse, reflected pulse and parameters of a medium where the pulse is propagating. In this work spherical objects are classified based on their filler material - water or air. We limit the study to spherical 2 layered targets immersed in water. The Form Function is used as a descriptor and fed into a Neural Network classifier, Multilayer Perceptron (MLP). The performance of the classifier is compared with Support Vector Machine (SVM) and the Form Function descriptor is examined in contrast to the Time and Frequency Representation of the echo.


title: End-to-End Speech Recognition: A review for the French Language

abstract: Recently, end-to-end ASR based either on sequence-to-sequence networks or on the CTC objective function gained a lot of interest from the community, achieving competitive results over traditional systems using robust but complex pipelines. One of the main features of end-to-end systems, in addition to the ability to free themselves from extra linguistic resources such as dictionaries or language models, is the capacity to model acoustic units such as characters, subwords or directly words; opening up the capacity to directly translate speech with different representations or levels of knowledge depending on the target language. In this paper we propose a review of the existing end-to-end ASR approaches for the French language. We compare results to conventional state-of-the-art ASR systems and discuss which units are more suited to model the French language.


title: Proof complexity of systems of (non-deterministic) decision trees and branching programs

abstract: This paper studies propositional proof systems in which lines are sequents of decision trees or branching programs - deterministic and nondeterministic. The systems LDT and LNDT are propositional proof systems in which lines represent deterministic or non-deterministic decision trees. Branching programs are modeled as decision dags. Adding extension to LDT and LNDT gives systems eLDT and eLNDT in which lines represent deterministic and non-deterministic branching programs, respectively. Deterministic and non-deterministic branching programs correspond to log-space (L) and nondeterministic log-space (NL). Thus the systems eLDT and eLNDT are propositional proof systems that reason with (nonuniform) L and NL properties. The main results of the paper are simulation and non-simulation results for tree-like and dag-like proofs in the systems LDT, LNDT, eLDT, and eLNDT. These systems are also compared with Frege systems, constantdepth Frege systems and extended Frege systems


title: Generalized Punctual Hilbert Schemes and $\mathfrak{g}$-complex structures

abstract: We define and analyze various generalizations of the punctual Hilbert scheme of the plane, associated to complex or real Lie algebras. Out of these, we construct new geometric structures on surfaces whose moduli spaces share multiple properties with Hitchin components, and which are conjecturally homeomorphic to them. For simple complex Lie algebras, this generalizes the higher complex structure. For real Lie algebras, this should give an alternative description of the Hitchin-Kostant-Rallis section.


title: The REDTOP experiment

abstract: The $\eta$ and $\eta$' mesons are almost unique in the particle universe since they are Goldstone and the dynamics of their decays are strongly constrained. Because the $\eta$ has no charge, decays that violate conservation laws can occur without interfering with a corresponding current. The integrated \eta-meson samples collected in earlier experiments amount to $\sim10^{9}$ events, dominated by the WASA-at-Cosy experiment. A new experiment, REDTOP (Rare Eta Decays with a TPC for Optical Photons), is being proposed, either at the proton Delivery Ring of Fermilab or at the Proton Synchrotron of CERN, with the intent of collecting in excess of $10^{12} \eta/yr$ ($10^{10}$ $\eta'/year$) for studying rare decays. Such statistics are sufficient for investigating several symmetry violations, and for searching particles beyond the Standard Model. Recent studies have indicated that REDTOP has very good sensitivity to processes that couple the Standard Model to new physics through three of the so-called portals: the vector, the scalar and the axion portal. An upgraded version of the detector will be proposed, at a later stage, for a run at the PIP-II accelerator complex, where the $\eta$ mesons will be produced in tagged mode. The physics program, the accelerator systems and the detector for REDTOP are discussed.


title: Adaptive Partitioning for Template Functions on Persistence Diagrams

abstract: As the field of Topological Data Analysis continues to show success in theory and in applications, there has been increasing interest in using tools from this field with methods for machine learning. Using persistent homology, specifically persistence diagrams, as inputs to machine learning techniques requires some mathematical creativity. The space of persistence diagrams does not have the desirable properties for machine learning, thus methods such as kernel methods and vectorization methods have been developed. One such featurization of persistence diagrams by Perea, Munch and Khasawneh uses continuous, compactly supported functions, referred to as "template functions," which results in a stable vector representation of the persistence diagram. In this paper, we provide a method of adaptively partitioning persistence diagrams to improve these featurizations based on localized information in the diagrams. Additionally, we provide a framework to adaptively select parameters required for the template functions in order to best utilize the partitioning method. We present results for application to example data sets comparing classification results between template function featurizations with and without partitioning, in addition to other methods from the literature.


title: Wiles defect for Hecke algebras that are not complete intersections

abstract: In his work on modularity theorems, Wiles proved a numerical criterion for a map of rings $R\to T$ to be an isomorphism of complete intersections. He used this to show that certain deformation rings and Hecke algebras associated to a mod $p$ Galois representation at non-minimal level were isomorphic and complete intersections, provided the same was true at minimal level. In this paper we study Hecke algebras acting on cohomology of Shimura curves arising from maximal orders in indefinite quaternion algebras over the rationals localized at a semistable irreducible mod $p$ Galois representation $\overline \rho$. If $\overline \rho$ is scalar at some primes dividing the discriminant of the quaternion algebra, then Hecke algebra is still isomorphic to the deformation ring, but is not a complete intersection, or even Gorenstein, so the Wiles numerical criterion cannot apply. We consider a weight 2 newform $f$ which contributes to the cohomology of the Shimura curve and gives rise to an augmentation $\lambda_f$ of the Hecke algebra. We quantify the failure of the Wiles numerical criterion at $\lambda_f$ by computing the associated {\it Wiles defect} purely in terms of the local behavior at primes dividing the discriminant of the global Galois representation $\rho_f$ which $f$ gives rise to by the Eichler--Shimura construction. One of the main tools used in the proof is Taylor--Wiles--Kisin patching.


title: Band edge localization beyond regular Floquet eigenvalues

abstract: We prove that localization near band edges of multi-dimensional ergodic random Schr\"odinger operators with periodic background potential in $L^2(\mathbb{R}^d)$ is universal. By this we mean that localization in its strongest dynamical form holds without extra assumptions on the random variables and independently of regularity or degeneracy of the Floquet eigenvalues of the background operator. The main novelty is an initial scale estimate the proof of which avoids Floquet theory altogether and uses instead an interplay between quantitative unique continuation and large deviation estimates. Furthermore, our reasoning is sufficiently flexible to prove this initial scale estimate in a non-ergodic setting, which promises to be an ingredient for understanding band edge localization also in these situations.


title: The need for adequate sampling in a well-functioning market surveillance system

abstract: Adequate sampling is essential for the well-functioning of a market surveillance system. As small as possible statistically significant sample size is the main factor that determines the costs of market surveillance actions. This paper studies various possibilities for calculation of the size of the sample with an emphasis on the method based on the binomial distribution. Examples, comparisons, and conclusions are provided.


title: Analysis of Nakamoto Consensus, Revisited

abstract: In the Bitcoin white paper, Nakamoto proposed a very simple Byzantine fault tolerant consensus algorithm that is also known as Nakamoto consensus. Despite its simplicity, some existing analysis of Nakamoto consensus appears to be long and involved. In this technical report, we aim to make such analysis simple and transparent so that we can teach senior undergraduate students and graduate students in our institutions. This report is largely based on a 3-hour tutorial given by one of the authors in June 2019.


title: Extreme eigenvalue statistics of $m$-dependent heavy-tailed matrices

abstract: We analyze the largest eigenvalue statistics of m-dependent heavy-tailed Wigner matrices as well as the associated sample covariance matrices having entry-wise regularly varying tail distributions with parameter $0<\alpha<4$. Our analysis extends results in the previous literature for the corresponding random matrices with independent entries above the diagonal, by allowing for m-dependence between the entries of a given matrix. We prove that the limiting point process of extreme eigenvalues is a Poisson cluster process.


title: Learning the piece-wise constant graph structure of a varying Ising model

abstract: This work focuses on the estimation of multiple change-points in a time-varying Ising model that evolves piece-wise constantly. The aim is to identify both the moments at which significant changes occur in the Ising model, as well as the underlying graph structures. For this purpose, we propose to estimate the neighborhood of each node by maximizing a penalized version of its conditional log-likelihood. The objective of the penalization is twofold: it imposes sparsity in the learned graphs and, thanks to a fused-type penalty, it also enforces them to evolve piece-wise constantly. Using few assumptions, we provide two change-points consistency theorems. Those are the first in the context of unknown number of change-points detection in time-varying Ising model. Finally, experimental results on several synthetic datasets and a real-world dataset demonstrate the performance of our method.


title: Rotating a supersolid dipolar gas

abstract: Distintictive features of supersolids show up in their rotational properties. We calculate the moment of inertia of a harmonically trapped dipolar Bose-Einstein condensed gas as a function of the tunable scattering length parameter, providing the transition from the (fully) superfluid to the supersolid phase and eventually to an incoherent crystal of self-bound droplets. The transition from the superfluid to the supersolid phase is characterized by a jump in the moment on inertia, revealing its first order nature. In the case of elongated trapping in the plane of rotation we show that the the moment of inertia determines the value of the frequency of the scissors mode, which is significantly affected by the reduction of superfluidity in the supersolid phase. The case of isotropic trapping is instead well suited to study the formation of quantized vortices, which are shown to be characterized, in the supersolid phase, by a sizeable deformed core, caused by the presence of the sorrounding density peaks.


title: Monodromy relations from twisted homology

abstract: We reformulate the monodromy relations of open-string scattering amplitudes as boundary terms of twisted homologies on the configuration spaces of Riemann surfaces of arbitrary genus. This allows us to write explicit linear relations involving loop integrands of open-string theories for any number of external particles and, for the first time, to arbitrary genus. In the non-planar sector, these relations contain seemingly unphysical contributions, which we argue clarify mismatches in previous literature. The text is mostly self-contained and presents a concise introduction to twisted homologies. As a result of this powerful formulation, we can propose estimates on the number of independent loop integrands based on Euler characteristics of the relevant configuration spaces, leading to a higher-genus generalization of the famous $(n-3)!$ result at genus zero.


title: Single and Cross-Dimensional Feature Detection and Description: An Evaluation

abstract: Three-dimensional local feature detection and description techniques are widely used for object registration and recognition applications. Although several evaluations of 3D local feature detection and description methods have already been published, these are constrained in a single dimensional scheme, i.e. either 3D or 2D methods that are applied onto multiple projections of the 3D data. However, cross-dimensional (mixed 2D and 3D) feature detection and description has yet to be investigated. Here, we evaluated the performance of both single and cross-dimensional feature detection and description methods on several 3D datasets and demonstrated the superiority of cross-dimensional over single-dimensional schemes.


title: Lattice study of electromagnetic conductivity of quark-gluon plasma in external magnetic field

abstract: We study the electromagnetic (e.m.) conductivity of QGP in a magnetic background by lattice simulations with $N_f = 2+1$ dynamical rooted staggered fermions at the physical point. We study the correlation functions of the e.m.~currents at $T=200,\,250$\,MeV and use the Tikhonov approach to extract the conductivity. This is found to rise with the magnetic field in the direction parallel to it and to decrease in the transverse direction, giving evidence for both the Chiral Magnetic Effect and the magnetoresistance phenomenon in QGP. We also estimate the chiral charge relaxation time in QGP.


title: Cluster Editing parameterized above modification-disjoint $P_3$-packings

abstract: Given a graph $G=(V,E)$ and an integer $k$, the Cluster Editing problem asks whether we can transform $G$ into a union of vertex-disjoint cliques by at most $k$ modifications (edge deletions or insertions). In this paper, we study the following variant of Cluster Editing. We are given a graph $G=(V,E)$, a packing $\cal H$ of modification-disjoint induced $P_3$s (no pair of $P_3$s in $\cal H$ share an edge or non-edge) and an integer $\ell$. The task is to decide whether $G$ can be transformed into a union of vertex-disjoint cliques by at most $\ell+|\cal H|$ modifications (edge deletions or insertions). We show that this problem is NP-hard even when $\ell=0$ (in which case the problem asks to turn $G$ into a disjoint union of cliques by performing exactly one edge deletion or insertion per element of $\cal H$) and when each vertex is in at most 23 $P_3$s of the packing. This answers negatively a question of van Bevern, Froese, and Komusiewicz (CSR 2016, ToCS 2018), repeated by C. Komusiewicz at Shonan meeting no. 144 in March 2019. We then initiate the study to find the largest integer $c$ such that the problem remains tractable when restricting to packings such that each vertex is in at most $c$ packed $P_3$s. Here packed $P_3$s are those belonging to the packing $\cal H$. Van Bevern et al. showed that the case $c = 1$ is fixed-parameter tractable with respect to $\ell$ and we show that the case $c = 2$ is solvable in $|V|^{2\ell + O(1)}$ time.


title: Pumping lemmas for classes of languages generated by folding systems

abstract: Geometric folding processes are ubiquitous in natural systems ranging from protein biochemistry to patterns of insect wings and leaves. In a previous study, a folding operation between strings of formal languages was introduced as a model of such processes. The operation was then used to define a folding system (F-system) as a construct consisting of a core language, containing the strings to be folded, and a folding procedure language, which defines how the folding is done. This paper reviews main definitions associated with F-systems and next it determines necessary conditions for a language to belong to classes generated by such systems. The conditions are stated in the form of pumping lemmas and four classes are considered, in which the core and folding procedure languages are both regular, one of them is regular and the other context-free, or both are context-free. Full demonstrations of the lemmas are provided, and the analysis is illustrated with examples.


title: Texture Bias Of CNNs Limits Few-Shot Classification Performance

abstract: Accurate image classification given small amounts of labelled data (few-shot classification) remains an open problem in computer vision. In this work we examine how the known texture bias of Convolutional Neural Networks (CNNs) affects few-shot classification performance. Although texture bias can help in standard image classification, in this work we show it significantly harms few-shot classification performance. After correcting this bias we demonstrate state-of-the-art performance on the competitive miniImageNet task using a method far simpler than the current best performing few-shot learning approaches.


title: Optimization Hierarchy for Fair Statistical Decision Problems

abstract: Data-driven decision-making has drawn scrutiny from policy makers due to fears of potential discrimination, and a growing literature has begun to develop fair statistical techniques. However, these techniques are often specialized to one model context and based on ad-hoc arguments, which makes it difficult to perform theoretical analysis. This paper develops an optimization hierarchy, which is a sequence of optimization problems with an increasing number of constraints, for fair statistical decision problems. Because our hierarchy is based on the framework of statistical decision problems, this means it provides a systematic approach for developing and studying fair versions of hypothesis testing, decision-making, estimation, regression, and classification. We use the insight that qualitative definitions of fairness are equivalent to statistical independence between the output of a statistical technique and a random variable that measures attributes for which fairness is desired. We use this insight to construct an optimization hierarchy that lends itself to numerical computation, and we use tools from variational analysis and random set theory to prove that higher levels of this hierarchy lead to consistency in the sense that it asymptotically imposes this independence as a constraint in corresponding statistical decision problems. We demonstrate numerical effectiveness of our hierarchy using several data sets, and we use our hierarchy to fairly perform automated dosing of morphine.


title: Fast Local Planning and Mapping in Unknown Off-Road Terrain

abstract: In this paper, we present a fast, on-line mapping and planning solution for operation in unknown, off-road, environments. We combine obstacle detection along with a terrain gradient map to make simple and adaptable cost map. This map can be created and updated at 10 Hz. An A* planner finds optimal paths over the map. Finally, we take multiple samples over the control input space and do a kinematic forward simulation to generated feasible trajectories. Then the most optimal trajectory, as determined by the cost map and proximity to A* path, is chosen and sent to the controller. Our method allows real time operation at rates of 30 Hz. We demonstrate the efficiency of our method in various off-road terrain at high speed.


title: A super-Earth and a mini-Neptune around Kepler-59

abstract: We characterize the radii and masses of the star and planets in the Kepler-59 system, as well as their orbital parameters. The star parameters are determined through a standard spectroscopic analysis, resulting in a mass of $1.359\pm 0.155\,M_\odot$ and a radius of $1.367\pm 0.078\,R_\odot$. The planetary radii obtained are $1.5\pm 0.1\,R_\oplus$ for the inner and $2.2\pm 0.1\,R_\oplus$ for the outer planet. The orbital parameters and the planetary masses are determined by the inversion of Transit Timing Variations (TTV) signals. For this, we consider two different data sets, one provided by Holczer et al. 2016, with TTVs only for the planet Kepler-59c, and the other provided by Rowe et al. 2015, with TTVs signals for both planets. The inversion method is carried out by applying an algorithm of Bayesian inference (MultiNest) combined with an efficient N-body integrator (Swift). For each of the data sets, two possible solutions are found, both having the same probability according to their corresponding Bayesian evidences. All four solutions appear to be indistinguishable within their 2-$\sigma$ uncertainties. Nevertheless, statistical analyses show that the solutions from Rowe et al. 2015 data better characterize the data. The first and second solutions identify masses of $5_{-2}^{+4}~M_{\mathrm{\oplus}}$ and $4.6_{-2.0}^{+3.6}~M_{\mathrm{\oplus}}$, and $3.0^{+0.8}_{-0.8}~M_{\mathrm{\oplus}}$ and $2.6^{+1.9}_{-0.8}~M_{\mathrm{\oplus}}$ for the inner and outer planet, respectively. This points to a system with an inner super-Earth and an outer mini-Neptune. Dynamical studies show the planets have almost co-planar orbits with small eccentricities ($e<0.1$), close but not into the 3:2 mean motion resonance. Stability analysis indicates that this configuration is stable over million years of evolution.


title: Magnetic field dependence of the thermopower of Kondo-correlated quantum dots

abstract: Recent experiments have measured the signatures of the Kondo effect in the zero-field thermopower of strongly correlated quantum dots [Svilans {\em et al.,} Phys. Rev. Lett. {\bf 121}, 206801 (2018); Dutta {\em et al.,} Nano Lett. {\bf 19}, 506 (2019)]. They confirm the predicted Kondo-induced sign change in the thermopower, upon increasing the temperature through a gate-voltage dependent value $T_{1}\gtrsim T_{\rm K}$, where $T_{\rm K}$ is the Kondo temperature. Here, we use the numerical renormalization group (NRG) method to investigate the effect of a finite magnetic field $B$ on the thermopower of such quantum dots. We show that, for fields $B$ exceeding a gate-voltage dependent value $B_{0}$, an additional sign change takes place in the Kondo regime at a temperature $T_{0}(B\geq B_{0})>0$ with $T_0<T_1$. The field $B_{0}$ is comparable to, but larger than, the field $B_{c}$ at which the zero-temperature spectral function splits in a magnetic field. The validity of the NRG results for $B_{0}$ are checked by comparison with asymptotically exact higher-order Fermi-liquid calculations [Oguri {\em et al.,} Phys. Rev. B {\bf 97}, 035435 (2018)]. Our calculations clarify the field-dependent signatures of the Kondo effect in the thermopower of Kondo-correlated quantum dots and explain the recently measured trends in the $B$-field dependence of the thermoelectric response of such systems [Svilans {\em et al.,} Phys. Rev. Lett. {\bf 121}, 206801 (2018)].


title: The full $\mathit{Kepler}$ phase curve of the eclipsing hot white dwarf binary system KOI-964

abstract: We analyze the full $\mathit{Kepler}$ phase curve of KOI-964, a binary system consisting of a hot white dwarf on an eclipsing orbit around an A-type host star. Using all 18 quarters of long-cadence photometry, we carry out a joint light curve fit and obtain improved phase curve amplitudes, occultation depths, orbital parameters, and transit ephemeris over the previous results of Carter et al. 2011. A periodogram of the residuals from the phase curve fit reveals an additional stellar variability signal from the host star with a characteristic period of $0.620276\pm0.000011$ days and a full amplitude of $24\pm2$ ppm. We also present new Keck/HIRES radial velocity observations which we use to measure the orbit and obtain a mass ratio of $q=0.106\pm0.012$. Combining this measurement with the results of a stellar isochrone analysis, we find that the masses of the host star and white dwarf companion are $2.23\pm0.12\,M_{\odot}$ and $0.236^{+0.028}_{-0.027}\,M_{\odot}$, respectively. The effective temperatures of the two components are $9940^{+260}_{-230}$ K and $15080\pm400$ K, respectively, and we determine the age of the system to be $0.21^{+0.11}_{-0.08}$ Gyr. We use the measured system properties to compute predicted phase curve amplitudes and find that while the measured Doppler boosting and mutual illumination components agree well with theory, the ellipsoidal distortion amplitude is significantly underestimated. We detail possible explanations for this discrepancy, including interactions between the dynamical tide of the host star and the tidal bulge and possible non-synchronous rotation of the host star.


title: MARTHE: Scheduling the Learning Rate Via Online Hypergradients

abstract: We study the problem of fitting task-specific learning rate schedules from the perspective of hyperparameter optimization, aiming at good generalization. We describe the structure of the gradient of a validation error w.r.t. the learning rate schedule -- the hypergradient. Based on this, we introduce MARTHE, a novel online algorithm guided by cheap approximations of the hypergradient that uses past information from the optimization trajectory to simulate future behaviour. It interpolates between two recent techniques, RTHO (Franceschi et al., 2017) and HD (Baydin et al. 2018), and is able to produce learning rate schedules that are more stable leading to models that generalize better.


title: A flexible integer linear programming formulation for scheduling clinician on-call service in hospitals

abstract: Scheduling of personnel in a hospital environment is vital to improving the service provided to patients and balancing the workload assigned to clinicians. Many approaches have been tried and successfully applied to generate efficient schedules in such settings. However, due to the computational complexity of the scheduling problem in general, most approaches resort to heuristics to find a non-optimal solution in a reasonable amount of time. We designed an integer linear programming formulation to find an optimal schedule in a clinical division of a hospital. Our formulation mitigates issues related to computational complexity by minimizing the set of constraints, yet retains sufficient flexibility so that it can be adapted to a variety of clinical divisions. We then conducted a case study for our approach using data from the Infectious Diseases division at St. Michael's Hospital in Toronto, Canada. We analyzed and compared the results of our approach to manually-created schedules at the hospital, and found improved adherence to departmental constraints and clinician preferences. We used simulated data to examine the sensitivity of the runtime of our linear program for various parameters and observed reassuring results, signifying the practicality and generalizability of our approach in different real-world scenarios.


title: Masked Gradient-Based Causal Structure Learning

abstract: This paper studies the problem of learning causal structures from observational data. We reformulate the Structural Equation Model (SEM) with additive noises in a form parameterized by binary graph adjacency matrix and show that, if the original SEM is identifiable, then the binary adjacency matrix can be identified up to super-graphs of the true causal graph under mild conditions. We then utilize the reformulated SEM to develop a causal structure learning method that can be efficiently trained using gradient-based optimization, by leveraging a smooth characterization on acyclicity and the Gumbel-Softmax approach to approximate the binary adjacency matrix. It is found that the obtained entries are typically near zero or one and can be easily thresholded to identify the edges. We conduct experiments on synthetic and real datasets to validate the effectiveness of the proposed method, and show that it readily includes different smooth model functions and achieves a much improved performance on most datasets considered.


title: Atomistic theory of electronic and optical properties of InAsP/InP nanowire quantum dots

abstract: We present here an atomistic theory of the electronic and optical properties of hexagonal InAsP quantum dots in InP nanowires in the wurtzite phase. These self-assembled quantum dots are unique in that their heights, shapes, and diameters are well known. Using a combined valence-force-field, tight-binding, and configuration-interaction approach we perform atomistic calculations of single-particle states and excitonic, biexcitonic and trion complexes as well as emission spectra as a function of the quantum dot height, diameter and As versus P concentration. The atomistic tight-binding parameters for InAs and InP in the wurtzite crystal phase were obtained by ab initio methods corrected by empirical band gaps. The low energy electron and hole states form electronic shells similar to parabolic or cylindrical quantum confinement, only weakly affected by hexagonal symmetry and As fluctuations. The relative alignment of the emission lines from excitons, trions and biexcitons agrees with that for InAs/InP dots in the zincblende phase in that biexcitons and positive trions are only weakly bound. The random distribution of As atoms leads to dot-to-dot fluctuations of a few meV for the single-particle states and the spectral lines. Due to the high symmetry of hexagonal InAsP nanowire quantum dots the exciton fine structure splitting is found to be small, of the order a few $\mu$eV with significant random fluctuations in accordance with experiments.


title: Prescribed-time Control for Perturbed Euler-Lagrange Systems with Obstacle Avoidance

abstract: This paper introduces a class of time-varying controllers for Euler-Lagrange systems such that the convergence occurs at an arbitrary finite time, independently of initial conditions, and free of chattering. The proposed controller is based on a mapping technique and is designed in two steps: First, a conventional (obstacle avoidance) asymptotically stable controller is specified for the nominal system; then, by a simple substitution, a prescribed-time (obstacle avoidance) controller is achievable for the perturbed system. It is proved that the proposed scheme is uniformly prescribed-time stable for unperturbed systems and prescribed-time attractive for perturbed systems as it rejects matched disturbances with unknown upper bounds without disturbance observation. As an example, a two-link robot manipulator is considered for numerical simulations.


title: Vapor drainage in the protolunar disk as the cause for the depletion in volatile elements of the Moon

abstract: Lunar rocks are severely depleted in moderately volatile elements such as Rb, K, and Zn relative to Earth. Identifying the cause of this depletion is important for understanding how the Earth-Moon system evolved in the aftermath of the Moon-forming giant impact. We measured the Rb isotopic compositions of lunar and terrestrial rocks to understand why moderately volatile elements are depleted in the Moon. Combining our new measurements with previous data reveals that the Moon has an 87Rb/85Rb ratio higher than Earth by +0.16 permil. This isotopic composition is consistent with evaporation of Rb into a vapor medium that was ~99 percent saturated. Evaporation under this saturation can also explain the previously documented isotopic fractionations of K, Ga, Cu and Zn of lunar rocks relative to Earth. We show that a possible setting for achieving the same saturation upon evaporation of elements with such diverse volatilities is through viscous drainage of a partially vaporized protolunar disk onto Earth. In the framework of an alpha-disk model, the alpha-viscosity needed to explain the ~99% saturation calculated here is 0.001 to 0.01, which is consistent with a vapor disk where viscosity is controlled by magnetorotational instability.


title: Healthy... Distress... Default

abstract: We discuss a simple, exactly solvable model of stochastic stock dynamics that incorporates regime switching between healthy and distressed regimes. Using this model, which is analytically tractable, we discuss a way of extracting expected returns for stocks from realized CDS spreads, essentially, the CDS market sentiment about future stock returns. This alpha/signal could be useful in a cross-sectional (statistical arbitrage) context for equities trading.


title: Extracting Dimensional Parameters of Gratings Produced with Self-Aligned Multiple Patterning Using GISAXS

abstract: Background: To ensure consistent and high-quality semiconductor production at future logic nodes, additional metrology tools are needed. For this purpose, grazing-incidence small-angle X-ray scattering (GISAXS) is being considered because measurements are fast with a proven capability to reconstruct average grating line profiles with high accuracy. Aim: GISAXS measurements of grating line shapes should be extended to samples with pitches smaller than 50 nm and their defects. The method's performance should be evaluated. Approach: A series of gratings with 32 nm pitch and deliberately introduced pitchwalk is measured using GISAXS. The grating line profiles with associated uncertainties are reconstructed using a Maxwell solver and Markov-Chain Monte Carlo (MCMC) sampling combined with a simulation library approach. Results: The line shape and the pitchwalk are generally in agreement with previously published transmission small-angle X-ray scattering (SAXS) results; however the line height and line width show deviations of (1.0 +/- 0.2) nm and (2.0 +/- 0.7) nm, respectively. The complex data evaluation leads to relatively high pitchwalk uncertainties between 0.5 nm and 2 nm. Conclusions: GISAXS shows great potential as a metrology tool for small-pitch line gratings with complex line profiles. Faster simulation methods would enable more accurate results.


title: A blinding solution for inference from astronomical data

abstract: This paper presents a joint blinding and deblinding strategy for inference of physical laws from astronomical data. The strategy allows for up to three blinding stages, where the data may be blinded, the computations of theoretical physics may be blinded, and --assuming Gaussianly distributed data-- the covariance matrix may be blinded. We found covariance blinding to be particularly effective, as it enables the blinder to determine close to exactly where the blinded posterior will peak. Accordingly, we present an algorithm which induces posterior shifts in predetermined directions by hiding untraceable biases in a covariance matrix. The associated deblinding takes the form of a numerically lightweight post-processing step, where the blinded posterior is multiplied with deblinding weights. We illustrate the blinding strategy for cosmic shear from KiDS-450, and show that even though there is no direct evidence of the KiDS-450 covariance matrix being biased, the famous cosmic shear tension with Planck could easily be induced by a mischaracterization of correlations between $\xi_-$ at the highest redshift and all lower redshifts. The blinding algorithm illustrates the increasing importance of accurate uncertainty assessment in astronomical inferences, as otherwise involuntary blinding through biases occurs.

categories: astro-ph.CO astro-ph.GA gr-qc physics.data-an


title: Many Faces of Feature Importance: Comparing Built-in and Post-hoc Feature Importance in Text Classification

abstract: Feature importance is commonly used to explain machine predictions. While feature importance can be derived from a machine learning model with a variety of methods, the consistency of feature importance via different methods remains understudied. In this work, we systematically compare feature importance from built-in mechanisms in a model such as attention values and post-hoc methods that approximate model behavior such as LIME. Using text classification as a testbed, we find that 1) no matter which method we use, important features from traditional models such as SVM and XGBoost are more similar with each other, than with deep learning models; 2) post-hoc methods tend to generate more similar important features for two models than built-in methods. We further demonstrate how such similarity varies across instances. Notably, important features do not always resemble each other better when two models agree on the predicted label than when they disagree.


title: Isogeometric analysis with piece-wise constant test functions

abstract: We focus on the finite element method computations with higher-order C1 continuity basis functions that preserve the partition of unity. We show that the rows of the system of linear equations can be combined, and the test functions can be sum up to 1 using the partition of unity property at the quadrature points. Thus, the test functions in higher continuity IGA can be set to piece-wise constants. This formulation is equivalent to testing with piece-wise constant basis functions, with supports span over some parts of the domain. The resulting method is a Petrov-Galerkin formulation with piece-wise constant test functions. This observation has the following consequences. The numerical integration cost can be reduced because we do not need to evaluate the test functions since they are equal to 1. This observation is valid for any basis functions preserving the partition of unity property. It is independent of the problem dimension and geometry of the computational domain. It also can be used in time-dependent problems, e.g., in the explicit dynamics computations, where we can reduce the cost of generation of the right-hand side. This summation of test functions can be performed for an arbitrary linear differential operator resulting from the Galerkin method applied to a PDE where we discretize with C1 continuity basis functions. The resulting method is equivalent to a linear combination of the collocations at points and with weights resulting from applied quadrature over the spans defined by supports of the piece-wise constant test functions.


title: LanCe: A Comprehensive and Lightweight CNN Defense Methodology against Physical Adversarial Attacks on Embedded Multimedia Applications

abstract: Recently, adversarial attacks can be applied to the physical world, causing practical issues to various Convolutional Neural Networks (CNNs) powered applications. Most existing physical adversarial attack defense works only focus on eliminating explicit perturbation patterns from inputs, ignoring interpretation to CNN's intrinsic vulnerability. Therefore, they lack the expected versatility to different attacks and thereby depend on considerable data processing costs. In this paper, we propose LanCe -- a comprehensive and lightweight CNN defense methodology against different physical adversarial attacks. By interpreting CNN's vulnerability, we find that non-semantic adversarial perturbations can activate CNN with significantly abnormal activations and even overwhelm other semantic input patterns' activations. We improve the CNN recognition process by adding a self-verification stage to detect the potential adversarial input with only one CNN inference cost. Based on the detection result, we further propose a data recovery methodology to defend the physical adversarial attacks. We apply such defense methodology into both image and audio CNN recognition scenarios and analyze the computational complexity for each scenario, respectively. Experiments show that our methodology can achieve an average 91% successful rate for attack detection and 89% accuracy recovery. Moreover, it is at most 3x faster compared with the state-of-the-art defense methods, making it feasible to resource-constrained embedded systems, such as mobile devices.


title: Normal Estimation for 3D Point Clouds via Local Plane Constraint and Multi-scale Selection

abstract: In this paper, we propose a normal estimation method for unstructured 3D point clouds. In this method, a feature constraint mechanism called Local Plane Features Constraint (LPFC) is used and then a multi-scale selection strategy is introduced. The LPEC can be used in a single-scale point network architecture for a more stable normal estimation of the unstructured 3D point clouds. In particular, it can partly overcome the influence of noise on a large sampling scale compared to the other methods which only use regression loss for normal estimation. For more details, a subnetwork is built after point-wise features extracted layers of the network and it gives more constraints to each point of the local patch via a binary classifier in the end. Then we use multi-task optimization to train the normal estimation and local plane classification tasks simultaneously.Also, to integrate the advantages of multi-scale results, a scale selection strategy is adopted, which is a data-driven approach for selecting the optimal scale around each point and encourages subnetwork specialization. Specifically, we employed a subnetwork called Scale Estimation Network to extract scale weight information from multi-scale features. More analysis is given about the relations between noise levels, local boundary, and scales in the experiment. These relationships can be a better guide to choosing particular scales for a particular model. Besides, the experimental result shows that our network can distinguish the points on the fitting plane accurately and this can be used to guide the normal estimation and our multi-scale method can improve the results well. Compared to some state-of-the-art surface normal estimators, our method is robust to noise and can achieve competitive results.


title: Quasinormal modes of five-dimensional black holes in non-commutative geometry

abstract: We compute the spectrum of quasinormal frequencies of five-dimensional black holes obtained in noncommutative geometry. In particular, we study scalar perturbations of a massive scalar field adopting the 6th order WKB approximation. We investigate in detail the impact of the mass of the scalar field, the angular degree and the overtone number on the spectrum. All modes are found to be stable.


title: On semigroup orbits of polynomials and multiplicative orders

abstract: We show, under some natural restrictions, that some semigroup orbits of polynomials cannot contain too many elements of small multiplicative order modulo a large prime $p$, extending previous work of Shparlinski (2017).


title: Semi-supervised Learning using Adversarial Training with Good and Bad Samples

abstract: In this work, we investigate semi-supervised learning (SSL) for image classification using adversarial training. Previous results have illustrated that generative adversarial networks (GANs) can be used for multiple purposes. Triple-GAN, which aims to jointly optimize model components by incorporating three players, generates suitable image-label pairs to compensate for the lack of labeled data in SSL with improved benchmark performance. Conversely, Bad (or complementary) GAN, optimizes generation to produce complementary data-label pairs and force a classifier's decision boundary to lie between data manifolds. Although it generally outperforms Triple-GAN, Bad GAN is highly sensitive to the amount of labeled data used for training. Unifying these two approaches, we present unified-GAN (UGAN), a novel framework that enables a classifier to simultaneously learn from both good and bad samples through adversarial training. We perform extensive experiments on various datasets and demonstrate that UGAN: 1) achieves state-of-the-art performance among other deep generative models, and 2) is robust to variations in the amount of labeled data used for training.


title: Joint Beamforming for Intelligent Reflecting Surface-Assisted Millimeter Wave Communications

abstract: Millimeter wave (mmWave) communications are evolving as a promising technology to meet the ever increasing data rate requirements. However, high directivity and severe path loss make it vulnerable to blockages, which could be frequent in indoor or urban environments. To address this issue, intelligent reflecting surfaces (IRSs) are introduced to provide additional adjustable reflected paths. Most prior works assume that elements of IRSs have an infinite phase resolution, which is difficult to be realized in practical systems. In this paper, IRSs with low-resolution phase shifters are considered. We aim to maximize the receive signal power at the user by jointly optimizing discrete phase shifts of IRSs and the transmit beamforming vector at the base station for mmWave downlink systems. An analytical near-optimal solution is developed by exploiting some important characteristics of mmWave channels. Our theoretical analysis reveals that low-resolution phase shifters can still achieve a receive signal power that increases quadratically with the number of reflecting elements. Simulation results are provided to corroborate our analysis and show the effectiveness of the proposed solution.


title: Circuit QED: single-step realization of a multiqubit controlled phase gate with one microwave photonic qubit simultaneously controlling $n-1$ microwave photonic qubits

abstract: We present a novel method to realize a multi-target-qubit controlled phase gate with one microwave photonic qubit simultaneously controlling $n-1$ target microwave photonic qubits. This gate is implemented with $n$ microwave cavities coupled to a superconducting flux qutrit. Each cavity hosts a microwave photonic qubit, whose two logic states are represented by the vacuum state and the single photon state of a single cavity mode, respectively. During the gate operation, the qutrit remains in the ground state and thus decoherence from the qutrit is greatly suppressed. This proposal requires only a single-step operation and thus the gate implementation is quite simple. The gate operation time is independent of the number of the qubits. In addition, this proposal does not need applying classical pulse or any measurement. Numerical simulations demonstrate that high-fidelity realization of a controlled phase gate with one microwave photonic qubit simultaneously controlling two target microwave photonic qubits is feasible with current circuit QED technology. The proposal is quite general and can be applied to implement the proposed gate in a wide range of physical systems, such as multiple microwave or optical cavities coupled to a natural or artificial $\Lambda$-type three-level atom.


title: Minimal automaton for multiplying and translating the Thue-Morse set

abstract: The Thue-Morse set $\mathcal{T}$ is the set of those non-negative integers whose binary expansions have an even number of $1$. The name of this set comes from the fact that its characteristic sequence is given by the famous Thue-Morse word ${\tt abbabaabbaababba\cdots}$, which is the fixed point starting with ${\tt a}$ of the word morphism ${\tt a\mapsto ab,b\mapsto ba}$. The numbers in $\mathcal{T}$ are commonly called the {\em evil numbers}. We obtain an exact formula for the state complexity of the set $m\mathcal{T}+r$ (i.e.\ the number of states of its minimal automaton) with respect to any base $b$ which is a power of $2$. Our proof is constructive and we are able to explicitly provide the minimal automaton of the language of all $2^p$-expansions of the set of integers $m\mathcal{T}+r$ for any positive integers $p$ and $m$ and any remainder $r\in\{0,\ldots,m-1\}$. The proposed method is general for any $b$-recognizable set of integers. As an application, we obtain a decision procedure running in quadratic time for the problem of deciding whether a given $2^p$-recognizable set is equal to a set of the form $m\mathcal{T}+r$.


title: Secular effects of Ultralight Dark Matter on Binary Pulsars

abstract: Dark matter (DM) can consist of very light bosons behaving as a classical scalar field that experiences coherent oscillations. The presence of this DM field would perturb the dynamics of celestial bodies, either because the (oscillating) DM stress tensor modifies the gravitational potentials of the galaxy, or if DM is directly coupled to the constituents of the body. We study secular variations of the orbital parameters of binary systems induced by such perturbations. Two classes of effects are identified. Effects of the first class appear if the frequency of DM oscillations is in resonance with the orbital motion; these exist for general DM couplings including the case of purely gravitational interaction. Effects of the second class arise if DM is coupled quadratically to the masses of the binary system members and do not require any resonant condition. The exquisite precision of binary pulsar timing can be used to constrain these effects. Current observations are not sensitive to oscillations in the galactic gravitational field, though a discovery of pulsars in regions of high DM density may improve the situation. For DM with direct coupling to ordinary matter, the current timing data are already competitive with other existing constraints in the range of DM masses $\sim 10^{-22}-10^{-18}\,{\rm eV}$. Future observations are expected to increase the sensitivity and probe new regions of parameters.


title: Escape of a forced-damped particle from weakly nonlinear truncated potential well

abstract: Escape from a potential well is an extreme example of transient behavior. We consider the escape of the harmonically forced particle under viscous damping from the benchmark truncated weakly nonlinear potential well. Main attention is paid to most interesting case of primary 1:1 resonance. The treatment is based on multiple-scales analysis and exploration of the slow-flow dynamics. Contrary to Hamiltonian case described in earlier works, in the case with damping the slow-flow equations are not integrable. However, if the damping is small enough, it is possible to analyze the perturbed slow-flow equations. The effect of the damping on the escape threshold is evaluated in the explicit analytic form. Somewhat unexpectedly, the escape mechanisms in terms of the slow flow are substantially different for the linear and weakly nonlinear cases.


title: Automatic sequences are also non-uniformly morphic

abstract: It is well-known that there exist infinite sequences that are the fixed point of non-uniform morphisms, but not $k$-automatic for any $k$. In this note we show that every $k$-automatic sequence is the image of a fixed point of a {\it non-uniform\/} morphism.


title: CDAG: A Serialized blockDAG for Permissioned Blockchain

abstract: Blockchain is maintained as a global log between a network of nodes and uses cryptographic distributed protocols to synchronize the updates. As adopted by Bitcoin and Ethereum these update operations to the ledger are serialized, and executed in batches. To safeguard the system against the generation of conflicting sets of updates and maintain the consistency of the ledger, the frequency of the updates is controlled, which severely affects the performance of the system. This paper presents Converging Directed Acyclic Graph (CDAG), as a substitute for the chain and DAG structures used in other blockchain protocols. CDAG allows multiple parallel updates to the ledger and converges them at the next step providing finality to the blocks. It partitions the updates into non-intersecting buckets of transactions to prevent the generation of conflicting blocks and divide the time into slots to provide enough time for them to propagate in the network. Multiple simultaneous updates improve the throughput of CDAG, and the converging step helps to finalize them faster, even in the presence of conflicts. Moreover, CDAG provides a total order among the blocks of the ledger to support smart contracts, unlike some of the other blockDAG protocols. We evaluate the performance of CDAG on Google Cloud Platform using Google Kubernetes Engine, simulating a real-time network. Experimental results show that CDAG achieves a throughput of more than 2000 transactions per second and confirms them well in under 2 minutes. Also, the protocol scales well in comparison to other permissioned protocols, and the capacity of the network only limits the performance.


title: An introduction to multiple orthogonal polynomials and Hermite-Pad\'e approximation

abstract: We present a brief introduction to the theory of multiple orthogonal polynomials on the basis of known results for an important class of measures known as Nikishin systems. For type I and type II multiple orthogonal polynomials with respect to such systems of measures, we describe some of their most relevant properties regarding location and distribution of zeros as well as their weak and ratio asymptotic behavior.


title: Towards Learning Cross-Modal Perception-Trace Models

abstract: Representation learning is a key element of state-of-the-art deep learning approaches. It enables to transform raw data into structured vector space embeddings. Such embeddings are able to capture the distributional semantics of their context, e.g. by word windows on natural language sentences, graph walks on knowledge graphs or convolutions on images. So far, this context is manually defined, resulting in heuristics which are solely optimized for computational performance on certain tasks like link-prediction. However, such heuristic models of context are fundamentally different to how humans capture information. For instance, when reading a multi-modal webpage (i) humans do not perceive all parts of a document equally: Some words and parts of images are skipped, others are revisited several times which makes the perception trace highly non-sequential; (ii) humans construct meaning from a document's content by shifting their attention between text and image, among other things, guided by layout and design elements. In this paper we empirically investigate the difference between human perception and context heuristics of basic embedding models. We conduct eye tracking experiments to capture the underlying characteristics of human perception of media documents containing a mixture of text and images. Based on that, we devise a prototypical computational perception-trace model, called CMPM. We evaluate empirically how CMPM can improve a basic skip-gram embedding approach. Our results suggest, that even with a basic human-inspired computational perception model, there is a huge potential for improving embeddings since such a model does inherently capture multiple modalities, as well as layout and design elements.


title: Theoretical study of deeply virtual Compton scattering off $^4$He

abstract: An interesting breakthrough in understanding the elusive inner content of nuclear systems in terms of partonic degrees of freedom is represented by deeply virtual Compton scattering processes. In such a way, tomographic view of nuclei and bound nucleons in coordinate space could be achieved for the first time. Moreover, nowadays experimental results for such a process considering $^4$He targets recently released at Jefferson Lab are available. In this talk, the recent results of our rigorous Impulse Approximation for DVCS off $^4$He, in terms of state-of-the-art models of the nuclear spectral function and of the parton structure of the bound proton, able to explain present data, has been shown.


title: Quantum Matrix Algebras of BMW type: Structure of the Characteristic Subalgebra

abstract: A notion of quantum matrix (QM-) algebra generalizes and unifies two famous families of algebras from the theory of quantum groups: the RTT-algebras and the reflection equation (RE-) algebras. These algebras being generated by the components of a `quantum' matrix $M$ possess certain properties which resemble structure theorems of the ordinary matrix theory. It turns out that such structure results are naturally derived in a more general framework of the QM-algebras. In this work we consider a family of Birman-Murakami-Wenzl (BMW) type QM-algebras. These algebras are defined with the use of R-matrix representations of the BMW algebras. Particular series of such algebras include orthogonal and symplectic types RTT- and RE- algebras, as well as their super-partners. For a family of BMW type QM-algebras, we investigate the structure of their `characteristic subalgebras' --- the subalgebras where the coefficients of characteristic polynomials take values. We define three sets of generating elements of the characteristic subalgebra and derive recursive Newton and Wronski relations between them. We also define an associative $\star$-product for the matrix $M$ of generators of the QM-algebra which is a proper generalization of the classical matrix multiplication. We determine the set of all matrix `descendants' of the quantum matrix $M$, and prove the $\star$-commutativity of this set in the BMW type.


title: Early-type galaxy density profiles from IllustrisTNG: II. Evolutionary trend of the total density profile

abstract: We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $\gamma^{\prime}$ of the total density profile for their main progenitors. We find that their $\gamma^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $\gamma^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $\gamma^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $\gamma^{\prime}$ and maintains the overall mass distribution close to isothermal.


title: Breaking a dark degeneracy: The gamma-ray signature of early matter domination

abstract: The Universe's early thermal history is poorly constrained, and it is possible that it underwent a period of early matter domination driven by a heavy particle or an oscillating scalar field that decayed into radiation before the onset of Big Bang nucleosynthesis. The entropy sourced by this particle's decay reduces the cross section required for thermal-relic dark matter to achieve the observed abundance. This degeneracy between dark matter properties and the thermal history vastly widens the field of viable dark matter candidates, undermining efforts to constrain dark matter's identity. Fortunately, an early matter-dominated era also amplifies density fluctuations at small scales and leads to early microhalo formation, boosting the dark matter annihilation rate and bringing smaller cross sections into the view of existing indirect-detection probes. Employing several recently developed models of microhalo formation and evolution, we develop a procedure to derive indirect-detection constraints on dark matter annihilation in cosmologies with early matter domination. This procedure properly accounts for the unique morphology of microhalo-dominated signals. While constraints depend on dark matter's free-streaming scale, the microhalos make it possible to obtain upper bounds as small as $\langle\sigma v\rangle \lesssim 10^{-32}$ cm$^3$s$^{-1}$ using Fermi-LAT observations of the isotropic gamma-ray background and the Draco dwarf galaxy.


title: Induced metal-free star formation around a massive black hole seed

abstract: The direct formation of a massive black hole is a potential seeding mechanism of the earliest observed supermassive black holes. We investigate how the existence of a massive black hole seed impacts the ionization and thermal state of its pre-galactic host halo and subsequent star formation. We show that its X-ray radiation ionizes and heats the medium, enhancing $\rm{H}_2$ formation in shielded regions, within the nuclear region in the span of a million years. The enhanced molecular cooling triggers the formation of a $\sim 10^4~{\rm M}_\odot$ metal-free stellar cluster at a star formation efficiency of $\sim 0.1\%$ in a single event. Star formation occurs near the edges of the H II region that is partially ionized by X-rays, thus the initial size depends on the black hole properties and surrounding environment. The simulated metal-free galaxy has an initial half-light radius of $\sim 10$ pc but expands to $\sim 50$ pc after 10 million years because of the outward velocities of their birth clouds. Supernova feedback then quenches any further star formation for tens of millions of years, allowing the massive black hole to dominate the spectrum once the massive metal-free stars die.


title: Supersymmetric Yang-Mills, Spherical Branes, and Precision Holography

abstract: Using supersymmetric localization we compute the free energy and BPS Wilson loop vacuum expectation values for planar maximally supersymmetric Yang-Mills theory on $S^d$ in the strong coupling limit for $2\leq d<6$. The same calculation can also be performed in supergravity using the recently found spherical brane solutions. We find excellent agreement between the two sets of results. This constitutes a non-trivial precision test of holography in a non-conformal setting. The free energy of maximal SYM on $S^6$ diverges in the strong coupling limit which might signify the onset of little string theory. We show how this divergence can be regularized both in QFT and in supergravity. We also consider $d=7$ with a small negative 't Hooft coupling and show that the free energy and Wilson loop vacuum expectation value agree with the results from supergravity after addressing some subtleties.


title: Electron-phonon and electron-electron interaction effects in twisted bilayer graphene

abstract: By comparing with recently available experimental data from several groups, we critically discuss the manifestation of continuum many body interaction effects in twisted bilayer graphene (tBLG) with small twist angles and low carrier densities, which arise naturally within the Dirac cone approximation for the non-interacting band structure. We provide two specific examples of such continuum many body theories: one involving electron-phonon interaction and one involving electron-electron interaction. In both cases, the experimental findings are only partially quantitatively consistent with rather clear-cut leading-order theoretical predictions based on well-established continuum many body theories. We provide a critical discussion, based mainly on the currently available tBLG experimental data, on possible future directions for understanding many body renormalization involving electron-phonon and electron-electron interactions in the system. One definitive conclusion based on the comparison between theory and experiment is that the leading order 1-loop perturbative renormalization group theory completely fails to account for the electron-electron interaction effects in the strong-coupling limit of flatband moir\'e tBLG system near the magic twist angle even at low doping where the Dirac cone approximation should apply. By contrast, approximate nonperturbative theoretical results based on Borel-Pad\'e resummation or $1/N$ expansion seems to work well compared with experiments, indicating rather small interaction corrections to Fermi velocity or carrier effective mass. For electron-phonon interactions, however, the leading-order continuum theory works well except when van Hove singularities in the density of states come into play.


title: Optimal in situ electromechanical sensing of molecular species

abstract: We investigate protocols for optimal molecular detection with electromechanical nanoscale sensors in ambient conditions. Our models are representative of suspended graphene nanoribbons, which due to their piezoelectric and electronic properties, provide responsive and versatile sensors. In particular, we analytically account for the corrections in the electronic transmission function and signal-to-noise ratio originating in environmental perturbations, such as thermal fluctuations and solvation effects. We also investigate the role of the sampling time in the current statistics. As a result, we formulate a protocol for optimal sensing based on the modulation of the Fermi level at fixed bias, and provide approximate forms for the current, linear susceptibility, and current fluctuations. We show how the algebraic tails in the thermally broadened transmission function affect the behavior of the signal-to-noise ratio and optimal sensing. These results provide further insights into the operation of graphene deflectometers and other techniques for electromechanical sensing.


title: CP violation and circular polarisation in neutrino radiative decay

abstract: The radiative decay of neutral fermions has been studied for decades but $CP$ violation induced within such a paradigm has evaded attention. $CP$ violation in these processes can produce an asymmetry between circularly polarised directions of the radiated photons and produces an important source of net circular polarisation in particle and astroparticle physics observables. The results presented in this work outlines the general connection between $CP$ violation and circular polarisation for both Dirac and Majorana fermions and can be used for any class of models that produce such radiative decays. The total $CP$ violation is calculated based on a widely studied Yukawa interaction considered in both active and sterile neutrino radiative decay scenarios as well as searches for dark matter via direct detection and collider signatures. Finally, the phenomenological implications of the formalism on keV sterile neutrino decay, leptogenesis-induced right-handed neutrino radiative decay and IceCube-driven heavy dark matter decay are discussed.


title: Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS). I. A Successful Machine Learning Selection of Metal-Poor Galaxies and the Discovery of a Galaxy with M*<10^6 M_sun and 0.016 Z_sun

abstract: We have initiated a new survey for local extremely metal-poor galaxies (EMPGs) with Subaru/Hyper Suprime-Cam (HSC) large-area (~500 deg^2) optical images reaching a 5 sigma limit of ~26 magnitude, about 100 times deeper than the Sloan Digital Sky Survey (SDSS). To select Z/Z_sun<0.1 EMPGs from ~40 million sources detected in the Subaru images, we first develop a machine-learning (ML) classifier based on a deep neural network algorithm with a training data set consisting of optical photometry of galaxy, star, and QSO models. We test our ML classifier with SDSS objects having spectroscopic metallicity measurements, and confirm that our ML classifier accomplishes 86%-completeness and 46%-purity EMPG classifications with photometric data. Applying our ML classifier to the photometric data of the Subaru sources as well as faint SDSS objects with no spectroscopic data, we obtain 27 and 86 EMPG candidates from the Subaru and SDSS photometric data, respectively. We conduct optical follow-up spectroscopy for 10 out of our EMPG candidates with Magellan/LDSS-3+MagE, Keck/DEIMOS, and Subaru/FOCAS, and find that the 10 EMPG candidates are star-forming galaxies at z=0.007-0.03 with large H_beta equivalent widths of 104-265 A, stellar masses of log(M*/M_sun)=5.0-7.1, and high specific star-formation rates of ~300 Gyr^{-1}, which are similar to those of early galaxies at z>6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun<0.1, one of which is HSC J1631+4426, the most metal-poor galaxy with Z/Z_sun=0.016 reported ever.


title: The impact of planet wakes on the location and shape of the water iceline in a protoplanetary disk

abstract: Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.


title: Repurposing lattice QCD results for composite phenomenology

abstract: A number of proposed extensions of the Standard Model include new strongly interacting dynamics, in the form of SU(N) gauge fields coupled to various numbers of fermions. Often, these extensions allow N = 3 as a plausible choice, or even require N = 3, such as in twin Higgs models, where the new dynamics is a "copy" of QCD. However, the fermion masses in such a sector are typically different from (often heavier than) the ones of real-world QCD, relative to the confinement scale. Many of the strong interaction masses and matrix elements for SU(3) at heavy fermion masses have already been computed on the lattice, typically as a byproduct of the approach to the physical point of real QCD. We provide a summary of these relevant results for the phenomenological community.


title: MiSTree: a Python package for constructing and analysing Minimum Spanning Trees

abstract: The minimum spanning tree (MST), a graph constructed from a distribution of points, draws lines between pairs of points so that all points are linked in a single skeletal structure that contains no loops and has minimal total edge length. The MST has been used in a broad range of scientific fields such as particle physics (to distinguish classes of events in collider collisions), in astronomy (to detect mass segregation in star clusters) and cosmology (to search for filaments in the cosmic web). Its success in these fields has been driven by its sensitivity to the spatial distribution of points and the patterns within. MiSTree, a public Python package, allows a user to construct the MST in a variety of coordinates systems, including Celestial coordinates used in astronomy. The package enables the MST to be constructed quickly by initially using a k-nearest neighbour graph (kNN, rather than a matrix of pairwise distances) which is then fed to Kruskal's algorithm to construct the MST. MiSTree enables a user to measure the statistics of the MST and provides classes for binning the MST statistics (into histograms) and plotting the distributions. Applying the MST will enable the inclusion of high-order statistics information from the cosmic web which can provide additional information to improve cosmological parameter constraints. This information has not been fully exploited due to the computational cost of calculating N-point statistics. MiSTree was designed to be used in cosmology but could be used in any field which requires extracting non-Gaussian information from point distributions. The source code for MiSTree is available on GitHub at https://github.com/knaidoo29/mistree


title: Analytic Functional Bootstrap for CFTs in $d>1$

abstract: We introduce analytic functionals which act on the crossing equation for CFTs in arbitrary spacetime dimension. The functionals fully probe the constraints of crossing symmetry on the first sheet, and are in particular sensitive to the OPE, (double) lightcone and Regge limits. Compatibility with the crossing equation imposes constraints on the functional kernels which we study in detail. We then introduce two simple classes of functionals. The first class has a simple action on generalized free fields and their deformations and can be used to bootstrap AdS contact interactions in general dimension. The second class is obtained by tensoring holomorphic and antiholomorphic copies of $d=1$ functionals which have been considered recently. They are dual to simple solutions to crossing in $d=2$ which include the energy correlator of the Ising model. We show how these functionals lead to optimal bounds on the OPE density of $d=2$ CFTs and argue that they provide an equivalent rewriting of the $d=2$ crossing equation which is better suited for numeric computations than current approaches.


title: Role of magnetically dominated disc-outflow symbiosis on bright hard-state black hole sources: ultra-luminous X-ray sources to quasars

abstract: We present optically thin solutions for magnetized, advective disc-outflow symbiosis around black holes (BHs). The main objective is to explain the bright, hard state observations of accreting systems with stellar-mass to supermassive BHs. We include the effects of magnetic fields and radiation counterpart in entropy gradient based on the first law of thermodynamics to represent energy advection. The cooling process includes bremsstrahlung, synchrotron radiation, and inverse Comptonization process. One of our main ventures is to explain some long-standing issues of ultra-luminous X-ray sources (ULXs). The existing physical scenarios to explain their unusual high luminosity are either the existence of the missing class of intermediate-mass BH (IMBH), or super-Eddington accretion around a stellar-mass BH. However, most ULXs with steep power-law spectrum can be well explained through super-Eddington accretion, while the existence of IMBH is indeed disputed extensively. Nevertheless, the interpretation of ULXs with a hard power-law dominated state remains mysterious. Here we show that our magnetically dominated disc-outflow symbiosis around rapidly spinning stellar-mass BHs can achieve such large luminosity even for sub-Eddington accretion rate. The magnetic field at the outer zone of the advective flow is more than the corresponding Eddington limit. Such field becomes dynamically dominant near the BH through continuous accretion process due to flux freezing, but maintaining its Eddington limit. This unique field configurations enhance the synchrotron, and synchrotron self-Comptonization process to achieve very large luminosity. Through the same mechanism, our solutions for supermassive BHs can explain the unusual large luminosity of ultra-luminous quasars.


title: Radiation-Hydrodynamical Models of X-ray Photoevaporation in Carbon Depleted Circumstellar Discs

abstract: The so-called transition discs provide an important tool to probe various mechanisms that might influence the evolution of protoplanetary discs and therefore the formation of planetary systems. One of these mechanisms is photoevaporation due to energetic radiation from the central star, which can in principal explain the occurrence of discs with inner cavities like transition discs. Current models, however, fail to reproduce a subset of the observed transition discs, namely objects with large measured cavities and vigorous accretion. For these objects the presence of (multiple) giant planets is often invoked to explain the observations. In our work we explore the possibility of X-ray photoevaporation operating in discs with different gas-phase depletion of carbon and show that the influence of photoevaporation can be extended in such low-metallicity discs. As carbon is one of the main contributors to the X-ray opacity, its depletion leads to larger penetration depths of X-rays in the disc and results in higher gas temperatures and stronger photoevaporative winds. We present radiation-hydrodynamical models of discs irradiated by internal X-ray+EUV radiation assuming Carbon gas-phase depletions by factors of 3,10 and 100 and derive realistic mass-loss rates and profiles. Our analysis yields robust temperature prescriptions as well as photoevaporative mass-loss rates and profiles which may be able to explain a larger fraction of the observed diversity of transition discs.


title: Modelling the emergence of phenotypic heterogeneity in vascularised tumours

abstract: We present a mathematical study of the emergence of phenotypic heterogeneity in vascularised tumours. Our study is based on formal asymptotic analysis and numerical simulations of a system of non-local parabolic equations that describes the phenotypic evolution of tumour cells and their nonlinear dynamic interactions with the oxygen, which is released from the intratumoural vascular network. Numerical simulations are carried out both in the case of arbitrary distributions of intratumour blood vessels and in the case where the intratumoural vascular network is reconstructed from clinical images obtained using dynamic optical coherence tomography. The results obtained support a more in-depth theoretical understanding of the eco-evolutionary process which underpins the emergence of phenotypic heterogeneity in vascularised tumours. In particular, our results offer a theoretical basis for empirical evidence indicating that the phenotypic properties of cancer cells in vascularised tumours vary with the distance from the blood vessels, and establish a relation between the degree of tumour tissue vascularisation and the level of intratumour phenotypic heterogeneity.


title: User Manual CAI version-1.0: An Open-Source Toolbox for Computer-Aided Investigation on the Fundamental Limits of Information Systems

abstract: We provide an open source toolbox on https://github.com/ct2641/CAI/releases/tag/1.0 to conduct computer-aided investigation on the fundamental limits of information systems. The toolbox relies on either Gurobi or Cplex as the linear program solving engine. The program can read a problem description file, and then fulfill the following tasks: 1) compute a bound for a given linear combination of information measures and provide the value of information measures at the optimal solution; 2) efficiently compute a polytope tradeoff outer bound between two information quantities; 3) produce a proof (as a weighted sum of known information inequalities; and 4) provide the range for information quantities between which the optimal value does not change (sensitivity analysis). This technical report provides an overview of this toolbox, a detailed description of the syntax of the problem description file, and a few example use cases.


title: Radio VLBA polarization and multi-band monitoring of the high-redshift quasar S5 0836+710 during a high activity period

abstract: We report on results of a multi-band monitoring campaign from radio to gamma rays of the high-redshift flat spectrum radio quasar S5 0836+710 during a high activity period detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Two major flares were detected, in 2015 August and November. In both episodes, the apparent isotropic gamma-ray luminosity exceeds 10^50 erg/s, with a doubling time scale of about 3 hours. The high gamma-ray activity may be related to a superluminal knot that emerged from the core in 2015 April at the peak of the radio activity and is moving downstream along the jet. The low variability observed in X-rays may indicate that X-ray emission is produced by the low-energy tail of the same electron population that produces the gamma-ray emission. The analysis of full-polarization pc-scale radio observations suggests the presence of a limb-brightened polarization structure at about 1 mas from the core in which a rotation measure gradient with a sign change is observed transverse to the jet direction. These characteristics are consistent with a scenario in which Faraday rotation is produced by a sheath of thermal electrons with a toroidal magnetic field surrounding the emitting jet.


title: Renormalization and Matching for the Collins-Soper Kernel from Lattice QCD

abstract: The Collins-Soper kernel, which governs the energy evolution of transverse-momentum dependent parton distribution functions (TMDPDFs), is required to accurately predict Drell-Yan like processes at small transverse momentum, and is a key ingredient for extracting TMDPDFs from experiment. Earlier we proposed a method to calculate this kernel from ratios of the so-called quasi-TMDPDFs determined with lattice QCD, which are defined as hadronic matrix elements of staple-shaped Euclidean Wilson line operators. Here we provide the one-loop renormalization of these operators in a regularization-independent momentum subtraction (RI$^\prime$/MOM) scheme, as well as the conversion factor from the RI$^\prime$/MOM-renormalized quasi-TMDPDF to the $\overline{\rm MS}$ scheme. We also propose a procedure for calculating the Collins-Soper kernel directly from position space correlators, which simplifies the lattice determination.


title: Magnetic order and single-ion anisotropy in Tb$_3$Ga$_5$O$_{12}$

abstract: Terbium gallium garnet (TGG), Tb$_3$Ga$_5$O$_{12}$, is well known for its applications in laser optics, but also exhibits complex low-temperature magnetism that is not yet fully understood. Its low-temperature magnetic order is determined by means of time-of-flight neutron powder diffraction. It is found to be a multiaxial antiferromagnet with magnetic Tb$^{3+}$ ions forming six sublattices of magnetic moments aligned parallel and anti-parallel to the $\langle100\rangle$ crystallographic directions of the cubic unit cell. The structure displays strong easy-axis anisotropy with respect to a two-fold axis of symmetry in the local orthorhombic environment of the Tb$^{3+}$ sites. The crystal-field splitting within the single-ion ground-state manifold is investigated by inelastic neutron scattering on powder samples. A strong temperature dependence of the quasidoublet ground-state is observed and revised parameters of the crystal-field Hamiltonian are given. The results of bulk magnetic susceptibility and magnetisation measurements are in good agreement with values based on the crystal-field model down to 20~K, where the onset of magnetic correlations is observed.