title: On commensurability of some right-angled Artin groups II: RAAGs defined by paths

abstract: In this paper we continue the study of right-angled Artin groups up to commensurability initiated in [CKZ]. We show that RAAGs defined by different paths of length greater than 3 are not commensurable. We also characterise which RAAGs defined by paths are commensurable to RAAGs defined by trees of diameter 4. More precisely, we show that a RAAG defined by a path of length $n>4$ is commensurable to a RAAG defined by a tree of diameter 4 if and only if $n$ is 2 modulo 4. These results follow from the connection that we establish between the classification of RAAGs up to commensurability and linear integer-programming.


title: Trade-off between information and disturbance in qubit thermometry

abstract: We address the trade-off between information and disturbance in qubit thermometry from the perspective of quantum estimation theory. Given a quantum measurement, we quantify information via the Fisher information of the measurement and disturbance via four different figures of merit, which capture different aspects (statistical, thermodynamical, geometrical) of the trade-off. For each disturbance measure, the efficient measurements, i.e. the measurements that introduce a disturbance not greater than any other measurement extracting the same amount of information, are determined explicitly. The family of efficient measurements varies with the choice of the disturbance measure. On the other hand, commutativity between the elements of the probability operator-valued measure (POVM) and the equilibrium state of the thermometer is a necessary condition for efficiency with respect to any figure of disturbance.


title: Series solution of Laplace problems

abstract: At the ANZIAM conference in Hobart in February, 2018, there were several talks on the solution of Laplace problems in multiply connected domains by means of conformal mapping. It appears to be not widely known that such problems can also be solved by the elementary method of series expansions with coefficients determined by least-squares fitting on the boundary. (These are not convergent series; the coefficients depend on the degree of the approximation.) Here we give a tutorial introduction to this method, which converges at an exponential rate if the boundary data are sufficiently well-behaved. The mathematical foundations go back to Runge in 1885 and Walsh in 1929. One of our examples involves an approximate Cantor set with up to 2048 components.


title: Hashing with Mutual Information

abstract: Binary vector embeddings enable fast nearest neighbor retrieval in large databases of high-dimensional objects, and play an important role in many practical applications, such as image and video retrieval. We study the problem of learning binary vector embeddings under a supervised setting, also known as hashing. We propose a novel supervised hashing method based on optimizing an information-theoretic quantity: mutual information. We show that optimizing mutual information can reduce ambiguity in the induced neighborhood structure in the learned Hamming space, which is essential in obtaining high retrieval performance. To this end, we optimize mutual information in deep neural networks with minibatch stochastic gradient descent, with a formulation that maximally and efficiently utilizes available supervision. Experiments on four image retrieval benchmarks, including ImageNet, confirm the effectiveness of our method in learning high-quality binary embeddings for nearest neighbor retrieval.


title: Entangled Pauli Principles: the DNA of Quantum Hall Fluids

abstract: A formalism is developed for the rigorous study of solvable fractional quantum Hall parent Hamiltonians with Landau level mixing. The idea of organization through "generalized Pauli principles" is expanded to allow for root level entanglement, giving rise to "entangled Pauli principles". Through the latter, aspects of the effective field theory description become ingrained in exact microscopic solutions for a great wealth of phases for which no similar single Landau level description is known. We discuss in detail braiding statistic, edge theory, and rigorous zero mode counting for the Jain-221 state as derived from a microscopic Hamiltonian. The relevant root-level entanglement is found to feature an AKLT-type MPS structure associated with an emergent SU(2) symmetry.


title: Sparse power-law network model for reliable statistical predictions based on sampled data

abstract: A projective network model is a model that enables predictions to be made based on a subsample of the network data, with the predictions remaining unchanged if a larger sample is taken into consideration. An exchangeable model is a model that does not depend on the order in which nodes are sampled. Despite a large variety of non-equilibrium (growing) and equilibrium (static) sparse complex network models that are widely used in network science, how to reconcile sparseness (constant average degree) with the desired statistical properties of projectivity and exchangeability is currently an outstanding scientific problem. Here we propose a network process with hidden variables which is projective and can generate sparse power-law networks. Despite the model not being exchangeable, it can be closely related to exchangeable uncorrelated networks as indicated by its information theory characterization and its network entropy. The use of the proposed network process as a null model is here tested on real data, indicating that the model offers a promising avenue for statistical network modelling.


title: Collective effects in Casimir-Polder forces

abstract: We study cooperative phenomena in the fluctuation-induced forces between a surface and a system of neutral two-level quantum emitters prepared in a coherent collective state, showing that the total Casimir-Polder force on the emitters can be modified via their mutual correlations. Particularly, we find that a collection of emitters prepared in a super- or subradiant state experiences an enhanced or suppressed collective vacuum-induced force, respectively. The collective nature of dispersion forces can be understood as resulting from the interference between the different processes contributing to the surface-modified resonant dipole-dipole interaction. Such cooperative fluctuation forces depend singularly on the surface response at the resonance frequency of the emitters, thus being easily maneuverable. Our results demonstrate the potential of collective phenomena as a new tool to selectively tailor vacuum forces.


title: Control of Battery Storage Systems for the Simultaneous Provision of Multiple Services

abstract: In this paper, we propose a control framework for a battery energy storage system to provide simultaneously multiple services to the electrical grid. The objective is to maximise the battery exploitation from these services in the presence of uncertainty (load, stochastic distributed generation, grid frequency). The framework is structured in two phases. In a period-ahead phase, we solve an optimization problem that allocates the battery power and energy budgets to the different services. In the subsequent real-time phase the control set-points for the deployment of such services are calculated separately and superimposed. The control framework is first formulated in a general way and then casted in the problem of providing dispatchability of a medium voltage feeder in conjunction to primary frequency control. The performance of the proposed framework are validated by simulations and real-scale experi- ments, performed with a grid-connected 560 kWh/720 kVA Li-ion battery energy storage system.


title: A lower bound for the number of elastic collisions

abstract: We prove by example that the number of elastic collisions of $n$ balls of equal mass and equal size in $d$-dimensional space can be greater than $n^3/27$ for $n\geq 3$ and $d\geq 2$. The previously known lower bound was of order $n^2$.


title: Estimation of Poisson arrival processes under linear models

abstract: In this paper we consider the problem of estimating the parameters of a Poisson arrival process where the rate function is assumed to lie in the span of a known basis. Our goal is to estimate the basis expansions coefficients given a realization of this process. We establish novel guarantees concerning the accuracy achieved by the maximum likelihood estimate. Our initial result is near-optimal, with the exception of an undesirable dependence on the dynamic range of the rate function. We then show how to remove this dependence through a process of "noise regularization", which results in an improved bound. We conjecture that a similar guarantee should be possible when using a more direct (deterministic) regularization scheme. We conclude with a discussion of practical applications and an empirical examination of the proposed regularization schemes.


title: On non-positive curvature properties of the Hilbert metric

abstract: In this paper, we consider different types of non-positive curvature properties of the Hilbert metric of a convex domain in R^n. First, we survey the relationships among the concepts and prove that in the case of Hilbert metric some of them are equivalent. Furthermore, we show some condition which implies the rigidity feature: if the Hilbert metric is Berwald, i.e., its Finslerian Chern connection reduces to a linear one, then the domain is an ellipsoid and the metric is Riemannian.


title: Excitons in few-layer hexagonal boron nitride: Davydov splitting and surface localization

abstract: Hexagonal boron nitride (hBN) has been attracting great attention because of its strong excitonic effects. Taking into account few-layer systems, we investigate theoretically the effects of the number of layers on quasiparticle energies, absorption spectra, and excitonic states, placing particular focus on the Davydov splitting of the lowest bound excitons. We describe how the inter-layer interaction as well as the variation in electronic screening as a function of layer number $N$ affects the electronic and optical properties. Using both \textit{ab initio} simulations and a tight-binding model for an effective Hamiltonian describing the excitons, we characterize in detail the symmetry of the excitonic wavefunctions and the selection rules for their coupling to incoming light. We show that for $N > 2$, one can distinguish between surface excitons that are mostly localized on the outer layers and inner excitons, leading to an asymmetry in the energy separation between split excitonic states. In particular, the bound surface excitons lie lower in energy than their inner counterparts. Additionally, this enables us to show how the layer thickness affects the shape of the absorption spectrum.


title: Power Control and Channel Allocation for D2D Underlaid Cellular Networks

abstract: Device-to-Device (D2D) communications underlaying cellular networks is a viable network technology that can potentially increase spectral utilization and improve power efficiency for proximitybased wireless applications and services. However, a major challenge in such deployment scenarios is the interference caused by D2D links when sharing the same resources with cellular users. In this work, we propose a channel allocation (CA) scheme together with a set of three power control (PC) schemes to mitigate interference in a D2D underlaid cellular system modeled as a random network using the mathematical tool of stochastic geometry. The novel aspect of the proposed CA scheme is that it enables D2D links to share resources with multiple cellular users as opposed to one as previously considered in the literature. Moreover, the accompanying distributed PC schemes further manage interference during link establishment and maintenance. The first two PC schemes compensate for large-scale path-loss effects and maximize the D2D sum rate by employing distance-dependent pathloss parameters of the D2D link and the base station, including an error estimation margin. The third scheme is an adaptive PC scheme based on a variable target signal-to-interference-plus-noise ratio, which limits the interference caused by D2D users and provides sufficient coverage probability for cellular users. Closed-form expressions for the coverage probability of cellular links, D2D links, and sum rate of D2D links are derived in terms of the allocated power, density of D2D links, and path-loss exponent. The impact of these key system parameters on network performance is analyzed and compared with previous work. Simulation results demonstrate an enhancement in cellular and D2D coverage probabilities, and an increase in spectral and power efficiency.


title: Topological Edge Modes by Smart Patterning

abstract: The research in topological materials and meta-materials reached maturity and is now gradually entering the phase of practical applications and devices. However, scaling down the experimental demonstrations definitely presents a challenge. In this work, we study coupled identical resonators whose collective dynamics is fully determined by the pattern in which the resonators are arranged. We call a pattern topological if boundary resonant modes fully fill all existing spectral gaps whenever the pattern is halved. This is a characteristic of the pattern and is entirely independent of the structure of the resonators and the details of the couplings. Existence of such patterns is proven using $K$-theory and exemplified using a novel experimental platform based on magnetically coupled spinners. Topological meta-materials built on these principles can be easily engineered at any scale, providing a practical platform for applications and devices.

categories: physics.app-ph cond-mat.mtrl-sci cond-mat.other math-ph math.MP


title: Hybrid Model For Word Prediction Using Naive Bayes and Latent Information

abstract: Historically, the Natural Language Processing area has been given too much attention by many researchers. One of the main motivation beyond this interest is related to the word prediction problem, which states that given a set words in a sentence, one can recommend the next word. In literature, this problem is solved by methods based on syntactic or semantic analysis. Solely, each of these analysis cannot achieve practical results for end-user applications. For instance, the Latent Semantic Analysis can handle semantic features of text, but cannot suggest words considering syntactical rules. On the other hand, there are models that treat both methods together and achieve state-of-the-art results, e.g. Deep Learning. These models can demand high computational effort, which can make the model infeasible for certain types of applications. With the advance of the technology and mathematical models, it is possible to develop faster systems with more accuracy. This work proposes a hybrid word suggestion model, based on Naive Bayes and Latent Semantic Analysis, considering neighbouring words around unfilled gaps. Results show that this model could achieve 44.2% of accuracy in the MSR Sentence Completion Challenge.


title: Niching an Archive-based Gaussian Estimation of Distribution Algorithm via Adaptive Clustering

abstract: As a model-based evolutionary algorithm, estimation of distribution algorithm (EDA) possesses unique characteristics and has been widely applied to global optimization. However, traditional Gaussian EDA (GEDA) may suffer from premature convergence and has a high risk of falling into local optimum when dealing with multimodal problem. In this paper, we first attempts to improve the performance of GEDA by utilizing historical solutions and develops a novel archive-based EDA variant. The use of historical solutions not only enhances the search efficiency of EDA to a large extent, but also significantly reduces the population size so that a faster convergence could be achieved. Then, the archive-based EDA is further integrated with a novel adaptive clustering strategy for solving multimodal optimization problems. Taking the advantage of the clustering strategy in locating different promising areas and the powerful exploitation ability of the archive-based EDA, the resultant algorithm is endowed with strong capability in finding multiple optima. To verify the efficiency of the proposed algorithm, we tested it on a set of well-known niching benchmark problems and compared it with several state-of-the-art niching algorithms. The experimental results indicate that the proposed algorithm is competitive.


title: Rashba Sandwiches with Topological Superconducting Phases

abstract: We introduce a versatile heterostructure harboring various topological superconducting phases characterized by the presence of helical, chiral, or unidirectional edge states. Changing parameters, such as an effective Zeeman field or chemical potential, one can tune between these three topological phases in the same setup. Our model relies only on conventional non-topological ingredients. The bilayer setup consists of an s-wave superconductor sandwiched between two two-dimensional electron gas layers with strong Rashba spin-orbit interaction. The interplay between two different pairing mechanisms, proximity induced direct and crossed Andreev superconducting pairings, gives rise to multiple topological phases. In particular, helical edge states occur if crossed Andreev superconducting pairing is dominant. In addition, an in-plane Zeeman field leads to a 2D gapless topological phase with unidirectional edge states, which were previously predicted to exist only in non-centrosymmetric superconductors. If the Zeeman field is tilted out of the plane, the system is in a topological phase hosting chiral edge states.


title: Cantor spectrum of graphene in magnetic fields

abstract: We consider a quantum graph as a model of graphene in magnetic fields and give a complete analysis of the spectrum, for all constant fluxes. In particular, we show that if the reduced magnetic flux $\Phi/2\pi$ through a honeycomb is irrational, the continuous spectrum is an unbounded Cantor set of Lebesgue measure zero.


title: Secure and Privacy-Aware Data Dissemination for Cloud-Based Applications

abstract: In this paper we propose a data dissemination platform that supports data security and different privacy levels even when the platform and the data are hosted by untrusted infrastructures. The proposed system aims at enabling an application ecosystem that uses off-the-shelf trusted platforms (in this case, Intel SGX), so that users may allow or disallow third parties to access the live data stream with a specific sensitivity-level. Moreover, this approach does not require users to manage the encryption keys directly. Our experiments show that such an approach is indeed practical for medium scale systems, where participants disseminate small volumes of data at a time, such as in smart grids and IoT environments.


title: Controlling the Deformation of Metamaterials: Corner Modes via Topology

abstract: Topological metamaterials have invaded the mechanical world, demonstrating acoustic cloaking and waveguiding at finite frequencies and variable, tunable elastic response at zero frequency. Zero frequency topological states have previously relied on the Maxwell condition, namely that the system has equal numbers of degrees of freedom and constraints. Here, we show that otherwise rigid periodic mechanical structures are described by a map with a nontrivial topological degree (a generalization of the winding number introduced by Kane and Lubensky) that creates, directs and protects modes on their boundaries. We introduce a model system consisting of rigid quadrilaterals connected via free hinges at their corners in a checkerboard pattern. This bulk structure generates a topological linear deformation mode exponentially localized in one corner, as investigated numerically and via experimental prototype. Unlike the Maxwell lattices, these structures select a single desired mode, which controls variable stiffness and mechanical amplification that can be incorporated into devices at any scale.


title: Dissipative hydrodynamics with higher-form symmetry

abstract: A theory of parity-invariant dissipative fluids with $q$-form symmetry is formulated to first order in a derivative expansion. The fluid is anisotropic with symmetry $\text{SO}(D-1-q)\times\text{SO}(q)$ and carries dissolved $q$-dimensional charged objects that couple to a $(q+1)$-form background gauge field. The case $q=1$ for which the fluid carries string charge is related to magnetohydrodynamics in $D=4$ spacetime dimensions. We identify $q$+7 parity-even independent transport coefficients at first order in derivatives for $q>1$. In particular, compared to the $q=1$ case under the assumption of parity and charge conjugation invariance, fluids with $q>1$ are characterised by $q$ extra transport coefficients with the physical interpretation of shear viscosity in the $\text{SO}(q)$ sector and current resistivities. We discuss certain issues related to the existence of a hydrostatic sector for fluids with higher-form symmetry for any $q\ge1$. We extend these results in order to include an interface separating different fluid phases and study the dispersion relation of capillary waves finding clear signatures of anisotropy. The formalism developed here can be easily adapted to study hydrodynamics with multiple higher-form symmetries.


title: Label Sanitization against Label Flipping Poisoning Attacks

abstract: Many machine learning systems rely on data collected in the wild from untrusted sources, exposing the learning algorithms to data poisoning. Attackers can inject malicious data in the training dataset to subvert the learning process, compromising the performance of the algorithm producing errors in a targeted or an indiscriminate way. Label flipping attacks are a special case of data poisoning, where the attacker can control the labels assigned to a fraction of the training points. Even if the capabilities of the attacker are constrained, these attacks have been shown to be effective to significantly degrade the performance of the system. In this paper we propose an efficient algorithm to perform optimal label flipping poisoning attacks and a mechanism to detect and relabel suspicious data points, mitigating the effect of such poisoning attacks.


title: Supernova 1987A Constraints on Sub-GeV Dark Sectors, Millicharged Particles, the QCD Axion, and an Axion-like Particle

abstract: We consider the constraints from Supernova 1987A on particles with small couplings to the Standard Model. We discuss a model with a fermion coupled to a dark photon, with various mass relations in the dark sector; millicharged particles; dark-sector fermions with inelastic transitions; the hadronic QCD axion; and an axion-like particle that couples to Standard Model fermions with couplings proportional to their mass. In the fermion cases, we develop a new diagnostic for assessing when such a particle is trapped at large mixing angles. Our bounds for a fermion coupled to a dark photon constrain small couplings and masses <200 MeV, and do not decouple for low fermion masses. They exclude parameter space that is otherwise unconstrained by existing accelerator-based and direct-detection searches. In addition, our bounds are complementary to proposed laboratory searches for sub-GeV dark matter, and do not constrain several "thermal" benchmark-model targets. For a millicharged particle, we exclude charges between 10^(-9) to a few times 10^(-6) in units of the electron charge; this excludes parameter space to higher millicharges and masses than previous bounds. For the QCD axion and an axion-like particle, we apply several updated nuclear physics calculations and include the energy dependence of the optical depth to accurately account for energy loss at large couplings. We rule out a hadronic axion of mass between 0.1 and a few hundred eV, or equivalently bound the PQ scale between a few times 10^4 and 10^8 GeV, closing the hadronic axion window. For an axion-like particle, our bounds disfavor decay constants between a few times 10^5 GeV up to a few times 10^8 GeV. In all cases, our bounds differ from previous work by more than an order of magnitude across the entire parameter space. We also provide estimated systematic errors due to the uncertainties of the progenitor.

categories: hep-ph astro-ph.CO astro-ph.GA astro-ph.HE hep-ex


title: A complex network description of thermal quantum states in the Ising spin chain

abstract: We use network analysis to describe and characterize an archetypal quantum system - an Ising spin chain in a transverse magnetic field. We analyze weighted networks for this quantum system, with link weights given by various measures of spin-spin correlations such as the von Neumann and Renyi mutual information, concurrence, and negativity. We analytically calculate the spin-spin correlations in the system at an arbitrary temperature by mapping the Ising spin chain to fermions, as well as numerically calculate the correlations in the ground state using matrix product state methods, and then analyze the resulting networks using a variety of network measures. We demonstrate that the network measures show some traits of complex networks already in this spin chain, arguably the simplest quantum many-body system. The network measures give insight into the phase diagram not easily captured by more typical quantities, such as the order parameter or correlation length. For example, the network structure varies with transverse field and temperature, and the structure in the quantum critical fan is different from the ordered and disordered phases.


title: Universal Behavior of the Thermoelectric Figure of Merit, zT, vs. Quality Factor

abstract: To increase the performance of thermoelectric materials, the electronic parameters in the figure of merit must be improved. In this paper, we use full, numerical band structures and solve the Boltzmann equation in the relaxation time approximation using energy-dependent scattering times informed by first principles simulations. By varying the strength of the electron-phonon coupling or the lattice thermal conductivity, we compute the thermoelectric figure of merit, zT, vs. a generalized thermoelectric quality factor. More than a dozen different complex electronic structures are examined. Surprisingly, we find that at a given quality factor, none provides a better figure of merit than that of a material with a simple, parabolic band and acoustic deformation potential scattering. A qualitative argument for this unexpected finding is presented. This apparent universal behavior suggests that even for complex electronic band structures, the thermoelectric figure of merit depends solely on the ratio of electrical to thermal conductivity; the Seebeck coefficient and Lorenz number need not be considered. This observation should simplify the search for promising new materials, but if exceptions to this behavior can be identified, new paths for increasing thermoelectric material performance will open up.


title: Towards 2+1 photon tomography: Energy-based selection of two 511 keV photons and a prompt photon with the J-PET scanner

abstract: The possibility to separate signals caused by 511 keV photons created in annihilation of electron-positron pairs and the so-called prompt photons from nuclei de- excitation is investigated. It could potentially be used to improve the quality of reconstructed images in the J-PET scanner in 2+1 photon tomography. Firstly, a research is conducted for several radioisotopes that decay via \b{eta}+ decay followed by de-excitation of an excited nucleus. Efficiency, purity and false positive rate are calculated for each isotope as a function of energy deposited threshold, with a hypothesis that signals caused by 511 keV photons deposit smaller values of energy than 1 z 13the selected threshold, while prompt photons deposit larger energy than the threshold. Analysis of the results accompanied with physical properties of radioisotopes suggests using 44 Sc, which is the most promising candidate for medical applications. With the use of GATE and J-POS simulation software, in-phantom scattering was introduced and the best energy deposited threshold value was estimated to be approximately 375 keV. It corresponds to almost 100% efficiency for 511 keV signals, 75% purity for 511 keV photons, and approximately 70% efficiency and purity for prompt photons.


title: A computational perspective of the role of Thalamus in cognition

abstract: Thalamus has traditionally been considered as only a relay source of cortical inputs, with hierarchically organized cortical circuits serially transforming thalamic signals to cognitively-relevant representations. Given the absence of local excitatory connections within the thalamus, the notion of thalamic `relay' seemed like a reasonable description over the last several decades. Recent advances in experimental approaches and theory provide a broader perspective on the role of the thalamus in cognitively-relevant cortical computations, and suggest that only a subset of thalamic circuit motifs fit the relay description. Here, we discuss this perspective and highlight the potential role for the thalamus -- and specifically mediodorsal (MD) nucleus -- in dynamic selection of cortical representations through a combination of intrinsic thalamic computations and output signals that change cortical network functional parameters. We suggest that through the contextual modulation of cortical computation, thalamus and cortex jointly optimize the information/cost tradeoff in an emergent fashion. We emphasize that coordinated experimental and theoretical efforts will provide a path to understanding the role of the thalamus in cognition, along with an understanding to augment cognitive capacity in health and disease.


title: Vu Ngoc's Conjecture on focus-focus singular fibers with multiple pinched points

abstract: We classify, up to fiberwise symplectomorphisms, a saturated neighborhood of a singular fiber of an integrable system (which is proper onto its image and has connected fibers) containing $k > 1$ focus-focus critical points. Our result shows that there is a one-to-one correspondence between such neighborhoods and $k$ formal power series, up to a $(\mathbb{Z}_2 \times D_k)$-action, where $D_k$ is the $k$-th dihedral group. The $k$ formal power series determine the dynamical behavior of the Hamiltonian vector fields associated to the components of the momentum map on the symplectic manifold $(M,\omega)$ near the singular fiber containing the $k$ focus-focus critical points. This proves a conjecture of San Vu Ngoc from 2002.


title: Spinon Fermi surface in a cluster Mott insulator model on a triangular lattice and possible application to 1T-TaS$_2$

abstract: 1T-TaS$_2$ is a cluster Mott insulator on the triangular lattice with 13 Ta atoms forming a star of David cluster as the unit cell. We derive a two dimensional XXZ spin-1/2 model with four-spin ring exchange term to describe the effective low energy physics of a monolayer 1T-TaS$_2$, where the effective spin-1/2 degrees of freedom arises from the Kramers degenerate spin-orbital states on each star of David. A large scale density matrix renormalization group simulation is further performed on this effective model and we find a gapless spin liquid phase with spinon Fermi surface at moderate to large strength region of four-spin ring exchange term. All peaks in the static spin structure factor are found to be located on the "$2k_F$" surface of half-filled spinon on the triangular lattice. Experiments to detect the spinon Fermi surface phase in 1T-TaS$_2$ are discussed.


title: Cognitive Radars: A Reality?

abstract: This paper describes some key ideas and applications of cognitive radars, highlighting the limits and the path forward. Cognitive radars are systems based on the perception-action cycle of cognition that sense the environment, learn from it relevant information about the target and the background, then adapt the radar sensor to optimally satisfy the needs of their mission according to a desired goal. The concept of cognitive radar was introduced originally for active radar only. In this paper we describe how this paradigm can be applied also to passive radar. In particular, we describe (i) cognitive active radars that work in a spectrally dense environment and change the transmitted waveform on-the-fly to avoid interference with the primary user of the channel, such as broadcast or communication systems, (ii) cognitive active radars that adjust transmit waveform parameters to achieve a specified level of target tracking performance, and (iii) cognitive passive radars, that contrary to the active radars cannot directly change the transmitted waveforms, but can instead select the best source of opportunity to improve detection and tracking performance.


title: The ASAS-SN Catalog of Variable Stars I: The Serendipitous Survey

abstract: The All-Sky Automated Survey for Supernovae (ASAS-SN) is the first optical survey to routinely monitor the whole sky with a cadence of $\sim2-3$ days down to V$\lesssim17$ mag. ASAS-SN has monitored the whole sky since 2014, collecting $\sim100-500$ epochs of observations per field. The V-band light curves for candidate variables identified during the search for supernovae are classified using a random forest classifier and visually verified. We present a catalog of 66,179 bright, new variable stars discovered during our search for supernovae, including 27,479 periodic variables and 38,700 irregular variables. V-band light curves for the ASAS-SN variables are available through the ASAS-SN variable stars database (https://asas-sn.osu.edu/variables). The database will begin to include the light curves of known variable stars in the near future along with the results for a systematic, all-sky variability survey.


title: Majorana Superconducting Qubit

abstract: We propose a platform for universal quantum computation that uses conventional $s$-wave superconducting leads to address a topological qubit stored in spatially separated Majorana bound states in a multi-terminal topological superconductor island. Both the manipulation and read-out of this "Majorana superconducting qubit" are realized by tunnel couplings between Majorana bound states and the superconducting leads. The ability of turning on and off tunnel couplings on-demand by local gates enables individual qubit addressability while avoiding cross-talk errors. By combining the scalability of superconducting qubit and the robustness of topological qubits, the Majorana superconducting qubit may provide a promising and realistic route towards quantum computation.


title: Modeling microlensing events with MulensModel

abstract: We introduce MulensModel, a software package for gravitational microlensing modeling. The package provides a framework for calculating microlensing model magnification curves and goodness-of-fit statistics for microlensing events with single and binary lenses as well as a variety of higher-order effects: extended sources with limb-darkening, annual microlensing parallax, satellite microlensing parallax, and binary lens orbital motion. The software could also be used for analysis of the planned microlensing survey by the NASA flag-ship WFIRST satellite. MulensModel is available at https://github.com/rpoleski/MulensModel/.


title: VISION - Vienna Survey in Orion II. Infrared extinction in Orion A

abstract: We have investigated the shape of the extinction curve in the infrared up to ~25 {\mu}m for the Orion A star-forming complex. The basis of this work is near-infrared data acquired with VISTA, in combination with Pan-STARRS and mid-infrared Spitzer photometry. We obtain colour excess ratios for eight passbands by fitting a series of colour-colour diagrams. The fits are performed using Markov chain Monte Carlo methods, together with a linear model under a Bayesian formalism. The resulting colour excess ratios are directly interpreted as a measure of the extinction law. We show that the Orion A molecular cloud is characterized by flat mid-infrared extinction, similar to many other recently studied sightlines. Moreover, we find statistically significant evidence that the extinction law from ~1 {\mu}m to at least ~6 {\mu}m varies across the cloud. In particular, we find a gradient along galactic longitude, where regions near the ONC show a different extinction law compared to L1641 and L1647, the low-mass star-forming sites in the cloud complex. These variations are of the order of only 3% and are most likely caused by the influence of the massive stars on their surrounding medium. While the observed general trends in our measurements are in agreement with model predictions, both well-established and new dust grain models are not able to fully reproduce our infrared extinction curve. We also present a new extinction map featuring a resolution of 1 arcmin and revisit the correlation between extinction and dust optical depth. This analysis shows that cloud substructure, which is not sampled by background sources, affects the conversion factor between these two measures. In conclusion, we argue that specific characteristics of the infrared extinction law are still not well understood, but Orion A can serve as an unbiased template for future studies.


title: Production of Cool Gas in Thermally-Driven Outflows

abstract: Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the CGOLS (Cholla Galactic OutfLow Simulations) suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ($v \sim 1000$ km/s), cool ($10^4$ K) gas observed in absorption line studies of outflows around some star-forming galaxies. By varying the mass-loading and geometry of the simulated winds, we identify a region of parameter space that leads to cool gas in outflows. In particular, when using an analytically-motivated central feedback model, we find that cooling flows can be produced with reasonable mass-loading rates ($\dot{M}_{wind} / \dot{M}_{SFR} \sim 0.5$), provided the star formation rate surface density is high. When a more realistic clustered feedback model is applied, destruction of high density clouds near the disk and interactions between different outflow regions indicate that lower mass-loading rates of the hot gas within the feedback region may still produce multiphase outflows. These results suggest an origin for fast-moving cool gas in outflows that does not rely on directly accelerating cool gas from the interstellar medium. These cooling flows may additionally provide an explanation for the multiphase gas ubiquitously observed in the halos of star-forming galaxies at low redshift.


title: Identification of the Central Compact Object in the young supernova remnant 1E0102.2-7219

abstract: Oxygen-rich young supernova remnants are valuable objects for probing the outcome of nucleosynthetic processes in massive stars, as well as the physics of supernova explosions. Observed within a few thousand years after the supernova explosion, these systems contain fast-moving oxygen-rich and hydrogen-poor filaments visible at optical wavelengths: fragments of the progenitor's interior expelled at a few 1000 km/s during the supernova explosion. Here we report the first identification of the compact object in 1E0102.2-7219 in reprocessed Chandra X-ray Observatory data, enabled via the discovery of a ring-shaped structure visible primarily in optical recombination lines of Ne I and O I. The optical ring, discovered in integral field spectroscopy observations from the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope, has a radius of $(2.10\pm0.35)$ arcsec = $(0.63\pm0.11)$ pc, and is expanding at a velocity of $90.5_{-30}^{+40}$ km/s. It surrounds an X-ray point source with an intrinsic X-ray luminosity $L_{i}$ (1.2--2.0 keV)= $(1.4\pm0.2)\times10^{33}$ erg/s. The energy distribution of the source indicates that this object is an isolated neutron star: a Central Compact Object akin to those present in the Cas A and Puppis A supernova remnants, and the first of its kind to be identified outside of our Galaxy.


title: Sowing black hole seeds: Direct collapse black hole formation with realistic Lyman-Werner radiation in cosmological simulations

abstract: We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation.


title: Introducing CGOLS: The Cholla Galactic OutfLow Simulation Suite

abstract: We present the Cholla Galactic OutfLow Simulations (CGOLS) suite, a set of extremely high resolution global simulations of isolated disk galaxies designed to clarify the nature of multiphase structure in galactic winds. Using the GPU-based code Cholla, we achieve unprecedented resolution in these simulations, modeling galaxies over a 20 kpc region at a constant resolution of 5 pc. The simulations include a feedback model designed to test the effects of different mass- and energy-loading factors on galactic outflows over kiloparsec scales. In addition to describing the simulation methodology in detail, we also present the results from an adiabatic simulation that tests the frequently-adopted analytic galactic wind model of Chevalier & Clegg. Our results indicate that the Chevalier & Clegg model is a good fit to nuclear starburst winds in the non-radiative region of parameter space. Lastly, we investigate the role of resolution and convergence in large-scale simulations of multiphase galactic winds. While our largest-scale simulations do show convergence of observable features like soft X-ray emission, our tests demonstrate that simulations of this kind with resolutions greater than 10 pc are not yet converged, confirming the need for extreme resolution in order to study the structure of winds and their effects on the circumgalactic medium.


title: Multi-messenger constraints to the local emission of cosmic-ray electrons

abstract: The data on the inclusive flux of cosmic positrons and electrons ($e^++e^{-}$) have been recently collected from GeV to tens of TeV energies by several experiments with unprecedented precision. In addition, the Fermi-LAT Collaboration has provided a new energy spectrum for the upper bounds on the $e^++e^{-}$ dipole anisotropy. This observable can bring information on the emission from local Galactic sources, notably measured with high precision at radio frequencies. We develop a framework in which $e^-$ and $e^+$ measured at Earth from GeV up to tens of TeV energies have a composite origin. A dedicated analysis is deserved to Vela YZ and Cygnus Loop Supernova Remnants (SNRs), for which we consider two different models for the injection of $e^-$. We investigate the consistency of these models using the three physical observables: the radio flux from Vela YZ and Cygnus Loop at all the available frequencies, the $e^++e^-$ flux from five experiments from the GeV to tens of TeV energy, the $e^++e^-$ dipole anisotropy upper limits from 50 GeV to about 1 TeV. We find that the radio flux for these nearby SNRs strongly constraints the properties of the injection electron spectrum, partially compatible with the looser constraints derived from the $e^+ + e^-$ flux data. We also perform a multi-wavelength multi-messenger analysis by fitting simultaneously the radio flux on Vela YZ and Cygnus Loop and the $e^+ + e^-$ flux, and checking the outputs against the $e^+ + e^-$ dipole anisotropy data. Remarkably, we find a model which is compatible with all the $e^++e^-$ flux data, the radio data for Vela YZ and Cygnus Loop, and with the anisotropy upper bounds. We show the severe constraints imposed by the most recent data on the $e^+ + e^-$ dipole anisotropy.


title: Quantized Rabi Oscillations and Circular Dichroism in Quantum Hall Systems

abstract: The dissipative response of a quantum system upon a time-dependent drive can be exploited as a probe of its geometric and topological properties. In this work, we explore the implications of such phenomena in the context of two-dimensional gases subjected to a uniform magnetic field. It is shown that a filled Landau level exhibits a quantized circular dichroism, which can be traced back to its underlying non-trivial topology. Based on selection rules, we find that this quantized circular dichroism can be suitably described in terms of Rabi oscillations, whose frequencies satisfy simple quantization laws. Moreover, we discuss how these quantized dissipative responses can be probed locally, both in the bulk and at the boundaries of the quantum Hall system. This work suggests alternative forms of topological probes in quantum systems based on circular dichroism.


title: Response and phase transition of a Kitaev spin liquid in a local magnetic field

abstract: We study the response of the Kitaev spin liquid (KSL) to a local magnetic field perpendicular to the Kitaev honeycomb lattice. The local magnetic field induces a dynamical excitation of a flux pair in the spin liquid and the system can be described by a generally particle-hole asymmetric interacting resonant level model. The dynamical excitation of the flux pair closes the flux gap in the spectrum of the spin correlation function locally for the gapless KSL even from the perturbative response to a weak magnetic field. Beyond the perturbative regime, the p-h asymmetry competes with the magnetic field and results in a rich phase diagram. Moreover, the magnetic field breaks the gauge equivalence of the ferromagnetic and anti-ferromagnetic Kitaev couplings of the ground state and leads to very different behaviors for the two cases. The anti-ferromagnetic case experiences a first order phase transition to the polarized state during magnetization whereas the ferromagnetic case does not. This study can be generalized to the Kitaev model in a uniform magnetic field and may help understand issues in recent experiments on KSL candidates.


title: Signatures of Plexcitonic States in Molecular Electroluminescence

abstract: We develop a quantum master equation (QME) approach to investigate the electroluminesence (EL) of molecules confined between metallic electrodes and coupled to quantum plasmonic modes. Within our general state-based framework, we describe electronic tunneling, vibrational damping, environmental dephasing, and the quantum coherent dynamics of coupled quantum electromagnetic field modes. As an example, we calculate the STM-induced spontaneous emission of a tetraphenylporphyrin (TPP) molecule coupled to a nanocavity plasmon. In the weak molecular exciton-plasmon coupling regime we find excellent agreement with experiments, including above-threshold hot luminescence, an effect not described by previous semiclassical calculations. In the strong coupling regime, we analyze the spectral features indicative of the formation of plexcitonic states.


title: An Overview of Robust Subspace Recovery

abstract: This paper will serve as an introduction to the body of work on robust subspace recovery. Robust subspace recovery involves finding an underlying low-dimensional subspace in a dataset that is possibly corrupted with outliers. While this problem is easy to state, it has been difficult to develop optimal algorithms due to its underlying nonconvexity. This work emphasizes advantages and disadvantages of proposed approaches and unsolved problems in the area.


title: Measurement of a Superconducting Qubit with a Microwave Photon Counter

abstract: Fast, high-fidelity measurement is a key ingredient for quantum error correction. Conventional approaches to the measurement of superconducting qubits, involving linear amplification of a microwave probe tone followed by heterodyne detection at room temperature, do not scale well to large system sizes. Here we introduce an alternative approach to measurement based on a microwave photon counter. We demonstrate raw single-shot measurement fidelity of 92%. Moreover, we exploit the intrinsic damping of the counter to extract the energy released by the measurement process, allowing repeated high-fidelity quantum non-demolition measurements. Crucially, our scheme provides access to the classical outcome of projective quantum measurement at the millikelvin stage. In a future system, counter-based measurement could form the basis for a scalable quantum-to-classical interface.


title: The Dirac equation as a quantum walk over the honeycomb and triangular lattices

abstract: A discrete-time Quantum Walk (QW) is essentially an operator driving the evolution of a single particle on the lattice, through local unitaries. Some QWs admit a continuum limit, leading to well-known physics partial differential equations, such as the Dirac equation. We show that these simulation results need not rely on the grid: the Dirac equation in $(2+1)$--dimensions can also be simulated, through local unitaries, on the honeycomb or the triangular lattice. The former is of interest in the study of graphene-like materials. The latter, we argue, opens the door for a generalization of the Dirac equation to arbitrary discrete surfaces.

categories: quant-ph cond-mat.mes-hall cs.DC hep-lat


title: Model-Free Control for Distributed Stream Data Processing using Deep Reinforcement Learning

abstract: In this paper, we focus on general-purpose Distributed Stream Data Processing Systems (DSDPSs), which deal with processing of unbounded streams of continuous data at scale distributedly in real or near-real time. A fundamental problem in a DSDPS is the scheduling problem with the objective of minimizing average end-to-end tuple processing time. A widely-used solution is to distribute workload evenly over machines in the cluster in a round-robin manner, which is obviously not efficient due to lack of consideration for communication delay. Model-based approaches do not work well either due to the high complexity of the system environment. We aim to develop a novel model-free approach that can learn to well control a DSDPS from its experience rather than accurate and mathematically solvable system models, just as a human learns a skill (such as cooking, driving, swimming, etc). Specifically, we, for the first time, propose to leverage emerging Deep Reinforcement Learning (DRL) for enabling model-free control in DSDPSs; and present design, implementation and evaluation of a novel and highly effective DRL-based control framework, which minimizes average end-to-end tuple processing time by jointly learning the system environment via collecting very limited runtime statistics data and making decisions under the guidance of powerful Deep Neural Networks. To validate and evaluate the proposed framework, we implemented it based on a widely-used DSDPS, Apache Storm, and tested it with three representative applications. Extensive experimental results show 1) Compared to Storm's default scheduler and the state-of-the-art model-based method, the proposed framework reduces average tuple processing by 33.5% and 14.0% respectively on average. 2) The proposed framework can quickly reach a good scheduling solution during online learning, which justifies its practicability for online control in DSDPSs.


title: A limit theorem for a class of stationary increments L\'{e}vy moving average process with multiple singularities

abstract: In this paper we present some new limit theorems for power variations of stationary increment L\'{e}vy driven moving average processes. Recently, such asymptotic results have been investigated in [Ann. Probab. 45(6B) (2017), 4477--4528, Festschrift for Bernt {\O}ksendal, Stochastics 81(1) (2017), 360--383] under the assumption that the kernel function potentially exhibits a singular behaviour at $0$. The aim of this work is to demonstrate how some of the results change when the kernel function has multiple singularity points. Our paper is also related to the article [Stoch. Process. Appl. 125(2) (2014), 653--677] that studied the same mathematical question for the class of Brownian semi-stationary models.


title: The EBEX Balloon Borne Experiment - Detectors and Readout

abstract: EBEX was a long-duration balloon-borne experiment to measure the polarization of the cosmic microwave background. The experiment had three frequency bands centered at 150, 250, and 410 GHz and was the first to use a kilo-pixel array of transition edge sensor (TES) bolometers aboard a balloon platform; shortly after reaching float we operated 504, 342, and 109 TESs at each of the bands, respectively. We describe the design and characterization of the array and the readout system. We give the distributions of measured thermal conductances, normal resistances, and transition temperatures. With the exception of the thermal conductance at 150 GHz. We measured median low-loop-gain time constants $\tau_{0}=$ 88, 46, and 57 ms and compare them to predictions. Two measurements of bolometer absorption efficiency show high ($\sim$0.9) efficiency at 150 GHz and medium ($\sim$0.35, and $\sim$0.25) at the two higher bands, respectively. We measure a median total optical load of 3.6, 5.3 and 5.0 pW absorbed at the three bands, respectively. EBEX pioneered the use of the digital version of the frequency domain multiplexing (FDM) system which multiplexed the bias and readout of 16 bolometers onto two wires. We present accounting of the measured noise equivalent power. The median per-detector noise equivalent temperatures referred to a black body with a temperature of 2.725 K are 400, 920, and 14500 $\mu$K$\sqrt{s}$ for the three bands, respectively. We compare these values to our pre-flight predictions and to a previous balloon payload, discuss the sources of excess noise, and the path for a future payload to make full use of the balloon environment.


title: Semidiscrete approximation to Benjamin-type equations

abstract: In this paper a semidiscrete Fourier pseudospectral method for approximating Benjamin-type equations is introduced and analyzed. A study of convergence is presented.


title: Bohmian mechanics and Fisher information for $q$-deformed Schr\"odinger equation

abstract: We discuss the Bohmian mechanics by means of the deformed Schr\"odinger equation for position dependent mass, in the context of a $q$-algebra inspired by nonextensive statistics. A deduction of the Bohmian quantum formalism is performed by means of a deformed Fisher information functional, from which a deformed Cram\'er-Rao bound is derived. Lagrangian and Hamiltonian formulations, inherited by the $q$-algebra, are also developed. Then, we illustrate the results with a particle confined in an infinite square potential well. The preservation of the deformed Cram\'er-Rao bound for the stationary states shows the role played by the $q$-algebraic structure.


title: Noise-induced backscattering in a quantum-spin-Hall edge

abstract: Time-reversal symmetry suppresses electron backscattering in a quantum-spin-Hall edge, yielding quantized conductance at zero temperature. Understanding the dominant corrections in finite-temperature experiments remains an unsettled issue. We study a novel mechanism for conductance suppression: backscattering caused by incoherent electromagnetic noise. Specifically, we show that an electric potential fluctuating randomly in time can backscatter electrons inelastically without constraints faced by electron-electron interactions. We quantify noise-induced corrections to the dc conductance in various regimes and propose an experiment to test this scenario.


title: Programming Quantum Computers Using Design Automation

abstract: Recent developments in quantum hardware indicate that systems featuring more than 50 physical qubits are within reach. At this scale, classical simulation will no longer be feasible and there is a possibility that such quantum devices may outperform even classical supercomputers at certain tasks. With the rapid growth of qubit numbers and coherence times comes the increasingly difficult challenge of quantum program compilation. This entails the translation of a high-level description of a quantum algorithm to hardware-specific low-level operations which can be carried out by the quantum device. Some parts of the calculation may still be performed manually due to the lack of efficient methods. This, in turn, may lead to a design gap, which will prevent the programming of a quantum computer. In this paper, we discuss the challenges in fully-automatic quantum compilation. We motivate directions for future research to tackle these challenges. Yet, with the algorithms and approaches that exist today, we demonstrate how to automatically perform the quantum programming flow from algorithm to a physical quantum computer for a simple algorithmic benchmark, namely the hidden shift problem. We present and use two tool flows which invoke RevKit. One which is based on ProjectQ and which targets the IBM Quantum Experience or a local simulator, and one which is based on Microsoft's quantum programming language Q$\#$.


title: Geodesic Orbit Riemannian Structures on $\mathbf{R}^n$

abstract: A geodesic orbit manifold is a complete Riemannian manifold all of whose geodesics are orbits of one-parameter groups of isometries. We give both a geometric and an algebraic characterization of geodesic orbit manifolds that are diffeomorphic to $\mathbf{R}^n$. Along the way, we establish various structural properties of more general geodesic orbit manifolds.


title: PRESISTANT: Learning based assistant for data pre-processing

abstract: Data pre-processing is one of the most time consuming and relevant steps in a data analysis process (e.g., classification task). A given data pre-processing operator (e.g., transformation) can have positive, negative or zero impact on the final result of the analysis. Expert users have the required knowledge to find the right pre-processing operators. However, when it comes to non-experts, they are overwhelmed by the amount of pre-processing operators and it is challenging for them to find operators that would positively impact their analysis (e.g., increase the predictive accuracy of a classifier). Existing solutions either assume that users have expert knowledge, or they recommend pre-processing operators that are only "syntactically" applicable to a dataset, without taking into account their impact on the final analysis. In this work, we aim at providing assistance to non-expert users by recommending data pre-processing operators that are ranked according to their impact on the final analysis. We developed a tool PRESISTANT, that uses Random Forests to learn the impact of pre-processing operators on the performance (e.g., predictive accuracy) of 5 different classification algorithms, such as J48, Naive Bayes, PART, Logistic Regression, and Nearest Neighbor. Extensive evaluations on the recommendations provided by our tool, show that PRESISTANT can effectively help non-experts in order to achieve improved results in their analytical tasks.


title: Derivations and differential operators on rings and fields

abstract: Let $R$ be an integral domain of characteristic zero. We prove that a function $D\colon R\to R$ is a derivation of order $n$ if and only if $D$ belongs to the closure of the set of differential operators of degree $n$ in the product topology of $R^R$, where the image space is endowed with the discrete topology. In other words, $f$ is a derivation of order $n$ if and only if, for every finite set $F\subset R$, there is a differential operator $D$ of degree $n$ such that $f=D$ on $F$. We also prove that if $d_1, \dots, d_n$ are nonzero derivations on $R$, then $d_1 \circ \ldots \circ d_n$ is a derivation of exact order $n$.


title: An optical and terahertz instrumentation system at the FAST linac at Fermilab

abstract: FAST is a facility at Fermilab that consists of a photoinjector, two superconducting capture cavities, one superconducting ILC-style cryomodule, and a small ring for studying non-linear, integrable beam optics called IOTA. This paper discusses the layout for the optical transport system that provides optical radiation to an externally located streak camera for bunch length measurements, and THz radiation to a Martin-Puplett interferometer, also for bunch length measurements. It accepts radiation from two synchrotron radiation ports in a chicane bunch compressor and a diffraction/transition radiation screen downstream of the compressor. It also has the potential to access signal from a transition radiation screen or YAG screen after the spectrometer magnet for measurements of energy-time correlations. Initial results from both the streak camera and Martin-Puplett will be presented.


title: Red Misfits in the Sloan Digital Sky Survey: Properties of Star-forming Red Galaxies

abstract: We study Red Misfits, a population of red, star-forming galaxies in the local Universe. We classify galaxies based on inclination-corrected optical colours and specific star formation rates derived from the Sloan Digital Sky Survey Data Release 7. Although the majority of blue galaxies are star-forming and most red galaxies exhibit little to no ongoing star formation, a small but significant population of galaxies ($\sim$11 per cent at all stellar masses) are classified as red in colour yet actively star-forming. We explore a number of properties of these galaxies and demonstrate that Red Misfits are not simply dusty or highly-inclined blue cloud galaxies or quiescent red galaxies with poorly-constrained star formation. The proportion of Red Misfits is nearly independent of environment and this population exhibits both intermediate morphologies and an enhanced likelihood of hosting an AGN. We conclude that Red Misfits are a transition population, gradually quenching on their way to the red sequence and this quenching is dominated by internal processes rather than environmentally-driven processes. We discuss the connection between Red Misfits and other transition galaxy populations, namely S0's, red spirals and green valley galaxies.


title: Investigating transverse Hall conductances using two-terminal setups

abstract: In this paper we present a method to numerically study transverse Hall conductances using a two-terminal setup. Using nonlinear transport concepts we find that the Hall voltage dependence on the model parameters can be investigated from the difference between the injectivities from each terminal. The method is suitable to work with non-equilibrium Green's functions as well as for scattering matrix approaches. We illustrate the proposed idea by studying the quantum spin Hall effect in graphene with disordered spin-orbit scattering centers induced by adatoms. We use two distinct models, a finite difference implementation of the Dirac Hamiltonian and a tight-binding Hamiltonian combined with the scattering matrix approach and the non-equilibrium Green's functions approach, respectively.


title: Tuning thermal transport in graphene via combinations of molecular antiresonances

abstract: We propose a method to engineer the phonon thermal transport properties of low dimensional systems. The method relies on introducing a predetermined combination of molecular adsorbates, which give rise to antiresonances at frequencies specific to the molecular species. Despite their dissimilar transmission spectra, thermal resistances due to individual molecules remain almost the same for all species. On the other hand, thermal resistance due to combinations of different species are not additive and show large differences depending on the species. Using a toy model, the physics underlying the violation of resistance summation rule is investigated. It is demonstrated that equivalent resistance of two scatterers having the same resistances can be close to the sum of the constituents or $\sim$70\% of it depending on the relative positions of the antiresonances. The relative positions of the antiresonances determine the net change in transmission, therefore the equivalent resistance. Since the entire spectrum is involved in phonon spectrum changes in different parts of the spectrum become important. Performing extensive first-principles based computations, we show that these distinctive attributes of phonon transport can be useful to tailor the thermal transport through low dimensional materials, especially for thermoelectric and thermal management applications.


title: On Mittag-Leffler d-orthogonal polynomials

abstract: This paper presents a first result of a long term research project dealing with the construction of d-orthogonal polynomials with Hahn's property. We shall show that the latter class could be characterized by expanding a polynomial as a finite sum of first derivatives of the elements of the sequence and we shall explain how this characterization could be used to construct Hahn-classical d-orthogonal polynomials as well. In this paper we look for solutions of linear combinations of the first derivatives of two consecutive elements of the sequence by considering the derivative operator and Delta (discrete) operator. The resulting polynomials constitute a particular class of Laguerre d-orthogonal polynomials and a generalization of Mittag-Leffler polynomials, respectively.


title: Unlimited parity alternating partitions

abstract: We introduce a new type of partitions that consists of partitions whose different parts alternate in parity (e.g., $3+2+2+1+1$). Various properties of this partition function are studied. In particular, we obtain its asymptotic behavior by employing Ingham's Tauberian theorem.


title: Drift parameter estimation for nonlinear stochastic differential equations driven by fractional Brownian motion

abstract: We derive the strong consistency of the least squares estimator for the drift coefficient of a fractional stochastic differential system. The drift coeffcient is one-sided dissipative Lipschitz and the driving noise is additive and fractional with Hurst parameter $H \in (\frac{1}{4}, 1)$. We assume that continuous observation is possible. The main tools are ergodic theorem and Malliavin calculus. As a by-product, we derive a maximum inequality for Skorohod integrals, which plays an important role to obtain the strong consistency of the least squares estimator.


title: Commissioning and first results of the Electron Beam Profiler in the Main Injector at Fermilab

abstract: The planned neutrino program at Fermilab requires large proton beam intensities in excess of 2 MW. Measuring the transverse profiles of these high intensity beams is challenging and often depends on non-invasive techniques. One such technique involves measuring the deflection of a probe beam of electrons with a trajectory perpendicular to the proton beam. A device such as this is already in use at the Spallation Neutron Source at ORNL and a similar device has been installed in the Main Injector at Fermilab. Commissioning of the device is in progress with the goal of having it operational by the end of the year. The status of the commissioning and initial results will be presented.


title: Control Contraction Metrics on Finsler Manifolds

abstract: Control Contraction Metrics (CCMs) provide a nonlinear controller design involving an offline search for a Riemannian metric and an online search for a shortest path between the current and desired trajectories. In this paper, we generalize CCMs to Finsler geometry, allowing the use of non-Riemannian metrics. We provide open loop and sampled data controllers. The sampled data control construction presented here does not require real time computation of globally shortest paths, simplifying computation.


title: Machine Learning Entanglement Freedom Or: How I Learned to Stop Worrying and Love Linear Regression

abstract: Quantum many-body systems realise many different phases of matter characterised by their exotic emergent phenomena. While some simple versions of these properties can occur in systems of free fermions, their occurrence generally implies that the physics is dictated by an interacting Hamiltonian. The interaction distance has been successfully used to quantify the effect of interactions in a variety of states of matter via the entanglement spectrum [Nat. Commun. 8, 14926 (2017), arXiv:1705.09983]. The computation of the interaction distance reduces to a global optimisation problem whose goal is to search for the free-fermion entanglement spectrum closest to the given entanglement spectrum. In this work, we employ techniques from machine learning in order to perform this same task. In a supervised learning setting, we use labelled data obtained by computing the interaction distance and predict its value via linear regression. Moving to a semi-supervised setting, we train an auto-encoder to estimate an alternative measure to the interaction distance, and we show that it behaves in a similar manner.


title: The temperature dependence of the shell corrections

abstract: We have examined the dependence of the shell correction to the nuclear liquid drop energy at finite excitations on the excitation energy (temperature). For this we have calculated the shell correction to the energy and free energy in very broad region of nuclei and deformations starting directly from their formal definitions. We have found out that the dependence of the shell corrections on the excitation energy differ substantially from the widely used approximation $\delta E(E^*)=\delta E(0)\exp(-E^*/E_d)$ both at small and large excitations. In particular, below the critical temperature at which the pairing effects vanish, the shell correction to the free energy is rather insensitive to the excitation energy. We suggest a more accurate approximation for the temperature dependence of the shell correction to the energy and free energy that is expressed in terms of the shell correction to the energy of independent particles and the shell correction to the pairing energy at T=0 and few fitted constants.


title: Nuclear quantum effects enter the mainstream

abstract: Over the past decades, atomistic simulations of chemical, biological and materials systems have become increasingly precise and predictive thanks to the development of accurate and efficient techniques that describe the quantum mechanical behavior of electrons. However, the overwhelming majority of such simulations still assume that the nuclei behave as classical particles. While historically this approximation could sometimes be justified due to complexity and computational overhead, the lack of nuclear quantum effects has become one of the biggest sources of error when systems containing light atoms are treated using current state-of-the-art descriptions of chemical interactions. Over the past decade, this realization has spurred a series of methodological advances that have led to dramatic reductions in the cost of including these important physical effects in the structure and dynamics of chemical systems. Here we show how these developments are now allowing nuclear quantum effects to become a mainstream feature of molecular simulations. These advances have led to new insights into chemical processes in the condensed phase and open the door to many exciting future opportunities.


title: Primordial gravitational waves, precisely: The role of thermodynamics in the Standard Model

abstract: In this paper, we revisit the estimation of the spectrum of primordial gravitational waves originated from inflation, particularly focusing on the effect of thermodynamics in the Standard Model of particle physics. By collecting recent results of perturbative and non-perturbative analysis of thermodynamic quantities in the Standard Model, we obtain the effective degrees of freedom including the corrections due to non-trivial interaction properties of particles in the Standard Model for a wide temperature interval. The impact of such corrections on the spectrum of primordial gravitational waves as well as the damping effect due to free-streaming particles is investigated by numerically solving the evolution equation of tensor perturbations in the expanding universe. It is shown that the reevaluation of the effects of free-streaming photons and neutrinos gives rise to some additional damping features overlooked in previous studies. We also observe that the continuous nature of the QCD crossover results in a smooth spectrum for modes that reenter the horizon at around the epoch of the QCD phase transition. Furthermore, we explicitly show that the values of the effective degrees of freedom remain smaller than the commonly used value 106.75 even at temperature much higher than the critical temperature of the electroweak crossover, and that the amplitude of primordial gravitational waves at a frequency range relevant to direct detection experiments becomes $\mathcal{O}(1)\,\%$ larger than previous estimates that do not include such corrections. This effect can be relevant to future high-sensitivity gravitational wave experiments such as ultimate DECIGO. Our results on the temperature evolution of the effective degrees of freedom are made available as tabulated data and fitting functions, which can also be used in the analysis of other cosmological relics.


title: Stepanov-like Weighted Pseudo-Almost Automorphic Solutions on Time Scales for a Novel High-order BAM Neural Network with Mixed Time-varying Delays in the Leakage Terms

abstract: We first propose the concept of Stepanov-like weighted pseudo-almost automorphic on time-space scales and we apply this type of oscillation to high-order BAM neural networks with mixed delays. Then, we study the existence and exponential stability of Sp-weighted pseudo almost automorphic on time-space scales solutions for the suggested system. Some criteria assuring the convergence are proposed. Our method is mainly based on the Banachs fixed point theorem, the theory of calculus on time scales and the Lyapunov-Krasovskii functional method. Moreover, a numerical example is given to show the effectiveness of the main results.


title: Potentials for $\mathcal{A}$-quasiconvexity

abstract: We show that each constant rank operator $\mathcal{A}$ admits an exact potential $\mathbb{B}$ in frequency space. We use this fact to show that the notion of $\mathcal{A}$-quasiconvexity can be tested against compactly supported fields. We also show that $\mathcal{A}$-free Young measures are generated by sequences $\mathbb{B}u_j$, modulo shifts by the barycentre.


title: Spin-wave induced phonon resonance in multiferroic BiFeO$_3$

abstract: We report the direct observation of a "resonance" mode in the lowest-energy optic phonon very near the zone center around (111) in the multiferroic BiFeO$_3$ using neutron scattering methods. The phonon scattering intensity is enhanced when antiferromagnetic (AFM) order sets in at T$_N = 640$~K, and it increases on cooling. This "resonance" is confined to a very narrow region in energy-momentum space where no spin-wave excitation intensity is expected, and it can be modified by an external magnetic field. Our results suggest the existence of a novel coupling between the lattice and spin fluctuations in this multiferroic system in which the spin-wave excitations are mapped onto the lattice vibrations via the Dzyaloshinskii-Moriya (DM) interaction.


title: On the Optimal Shape of Tree Roots and Branches

abstract: This paper introduces two classes of variational problems, determining optimal shapes for tree roots and branches. Given a measure $\mu$, describing the distribution of leaves, we introduce a sunlight functional $\S(\mu)$ computing the total amount of light captured by the leaves. On the other hand, given a measure $\mu$ describing the distribution of root hair cells, we consider a harvest functional $\H(\mu)$ computing the total amount of water and nutrients gathered by the roots. In both cases, we seek to maximize these functionals subject to a ramified transportation cost, for transporting nutrients from the roots to the trunk and from the trunk to the leaves. The main results establish various properties of these functionals, and the existence of optimal distributions. In particular, we prove the upper semicontinuity of $\S$ and $\H$, together with a priori estimates on the support of optimal distributions.


title: Building a Telescope to Look Into High-Dimensional Image Spaces

abstract: An image pattern can be represented by a probability distribution whose density is concentrated on different low-dimensional subspaces in the high-dimensional image space. Such probability densities have an astronomical number of local modes corresponding to typical pattern appearances. Related groups of modes can join to form macroscopic image basins that represent pattern concepts. Recent works use neural networks that capture high-order image statistics to learn Gibbs models capable of synthesizing realistic images of many patterns. However, characterizing a learned probability density to uncover the Hopfield memories of the model, encoded by the structure of the local modes, remains an open challenge. In this work, we present novel computational experiments that map and visualize the local mode structure of Gibbs densities. Efficient mapping requires identifying the global basins without enumerating the countless modes. Inspired by Grenander's jump-diffusion method, we propose a new MCMC tool called Attraction-Diffusion (AD) that can capture the macroscopic structure of highly non-convex densities by measuring metastability of local modes. AD involves altering the target density with a magnetization potential penalizing distance from a known mode and running an MCMC sample of the altered density to measure the stability of the initial chain state. Using a low-dimensional generator network to facilitate exploration, we map image spaces with up to 12,288 dimensions (64 $\times$ 64 pixels in RGB). Our work shows: (1) AD can efficiently map highly non-convex probability densities, (2) metastable regions of pattern probability densities contain coherent groups of images, and (3) the perceptibility of differences between training images influences the metastability of image basins.


title: Multi-Agent Imitation Learning for Driving Simulation

abstract: Simulation is an appealing option for validating the safety of autonomous vehicles. Generative Adversarial Imitation Learning (GAIL) has recently been shown to learn representative human driver models. These human driver models were learned through training in single-agent environments, but they have difficulty in generalizing to multi-agent driving scenarios. We argue these difficulties arise because observations at training and test time are sampled from different distributions. This difference makes such models unsuitable for the simulation of driving scenes, where multiple agents must interact realistically over long time horizons. We extend GAIL to address these shortcomings through a parameter-sharing approach grounded in curriculum learning. Compared with single-agent GAIL policies, policies generated by our PS-GAIL method prove superior at interacting stably in a multi-agent setting and capturing the emergent behavior of human drivers.


title: Quantitatively Evaluating GANs With Divergences Proposed for Training

abstract: Generative adversarial networks (GANs) have been extremely effective in approximating complex distributions of high-dimensional, input data samples, and substantial progress has been made in understanding and improving GAN performance in terms of both theory and application. However, we currently lack quantitative methods for model assessment. Because of this, while many GAN variants are being proposed, we have relatively little understanding of their relative abilities. In this paper, we evaluate the performance of various types of GANs using divergence and distance functions typically used only for training. We observe consistency across the various proposed metrics and, interestingly, the test-time metrics do not favour networks that use the same training-time criterion. We also compare the proposed metrics to human perceptual scores.


title: A Class of Multiparameter Oscillatory Singular Integral Operators: Endpoint Hardy Space Bounds

abstract: We establish endpoint bounds on a Hardy space $H^1$ for a natural class of multiparameter singular integral operators which do not decay away from the support of rectangular atoms. Hence the usual argument via a Journ\'e-type covering lemma to deduce bounds on product $H^1$ is not valid. We consider the class of multiparameter oscillatory singular integral operators given by convolution with the classical multiple Hilbert transform kernel modulated by a general polynomial oscillation. Various characterisations are known which give $L^2$ (or more generally $L^p, 1<p<\infty$) bounds. Here we initiate an investigation of endpoint bounds on the rectangular Hardy space $H^1$ in two dimensions; we give a characterisation when bounds hold which are uniform over a given subspace of polynomials and somewhat surprisingly, we discover that the Hardy space and $L^p$ theories for these operators are very different.


title: Unsupervised Odometry and Depth Learning for Endoscopic Capsule Robots

abstract: In the last decade, many medical companies and research groups have tried to convert passive capsule endoscopes as an emerging and minimally invasive diagnostic technology into actively steerable endoscopic capsule robots which will provide more intuitive disease detection, targeted drug delivery and biopsy-like operations in the gastrointestinal(GI) tract. In this study, we introduce a fully unsupervised, real-time odometry and depth learner for monocular endoscopic capsule robots. We establish the supervision by warping view sequences and assigning the re-projection minimization to the loss function, which we adopt in multi-view pose estimation and single-view depth estimation network. Detailed quantitative and qualitative analyses of the proposed framework performed on non-rigidly deformable ex-vivo porcine stomach datasets proves the effectiveness of the method in terms of motion estimation and depth recovery.


title: Magnetic-Visual Sensor Fusion-based Dense 3D Reconstruction and Localization for Endoscopic Capsule Robots

abstract: Reliable and real-time 3D reconstruction and localization functionality is a crucial prerequisite for the navigation of actively controlled capsule endoscopic robots as an emerging, minimally invasive diagnostic and therapeutic technology for use in the gastrointestinal (GI) tract. In this study, we propose a fully dense, non-rigidly deformable, strictly real-time, intraoperative map fusion approach for actively controlled endoscopic capsule robot applications which combines magnetic and vision-based localization, with non-rigid deformations based frame-to-model map fusion. The performance of the proposed method is demonstrated using four different ex-vivo porcine stomach models. Across different trajectories of varying speed and complexity, and four different endoscopic cameras, the root mean square surface reconstruction errors 1.58 to 2.17 cm.


title: Classical Transitions

abstract: We introduce the calculus of Classical Transitions (CT), which extends the research line on the relationship between linear logic and processes to labelled transitions. The key twist from previous work is registering parallelism in typing judgements, by generalising linear logic judgements from one sequents to many (hypersequents). This allows us to bridge the gap between the structures of operators used as proof terms in previous work and those of the standard {\pi}-calculus (in particular parallel operator and restriction). The proof theory of CT allows for new proof transformations, which we show correspond to a labelled transition system (LTS) for processes. We prove that CT enjoys subject reduction and progress.


title: SN 2017ein and the Possible First Identification of a Type Ic Supernova Progenitor

abstract: We have identified a progenitor candidate in archival Hubble Space Telescope (HST) images for the Type Ic SN 2017ein in NGC 3938, pinpointing the candidate's location via HST Target-of-Opportunity imaging of the SN itself. This would be the first identification of a stellar-like object as a progenitor candidate for any Type Ic supernova to date. We also present observations of SN 2017ein during the first ~49 days since explosion. We find that SN 2017ein most resembles the well-studied Type Ic SN 2007gr. We infer that SN 2017ein experienced a total visual extinction of A_V~1.0--1.9 mag, predominantly because of dust within the host galaxy. Although the distance is not well known, if this object is the progenitor, it was likely of high initial mass, ~47--48 M_sun if a single star, or ~60--80 M_sun if in a binary system. However, we also find that the progenitor candidate could be a very blue and young compact cluster, further implying a very massive (>65 M_sun) progenitor. Furthermore, the actual progenitor might not be associated with the candidate at all and could be far less massive. From the immediate stellar environment, we find possible evidence for three different populations; if the SN progenitor was a member of the youngest population, this would be consistent with an initial mass of ~57 M_sun. After it has faded, the SN should be reobserved at high spatial resolution and sensitivity, to determine whether the candidate is indeed the progenitor.


title: Bounding $\chi$ by a fraction of $\Delta$ for graphs without large cliques

abstract: The greedy coloring algorithm shows that a graph of maximum degree at most $\Delta$ has chromatic number at most $\Delta + 1$, and this is tight for cliques. Much attention has been devoted to improving this "greedy bound" for graphs without large cliques. Brooks famously proved that this bound can be improved by one if $\Delta \geq 3$ and the graph contains no clique of size $\Delta + 1$. Reed's Conjecture states that the "greedy bound" can be improved by $k$ if the graph contains no clique of size $\Delta + 1 - 2k$. Johansson proved that the "greedy bound" can be improved by a factor of $\Omega(\ln(\Delta)^{-1})$ or $\Omega\left(\frac{\ln(\ln(\Delta))}{\ln(\Delta)}\right)$ for graphs with no triangles or no cliques of any fixed size, respectively. Notably missing is a $\textit{linear}$ improvement on the "greedy bound" for graphs without large cliques. In this paper, we prove that for sufficiently large $\Delta$, if $G$ is a graph with maximum degree at most $\Delta$ and no clique of size $\omega$, then $$\chi(G) \leq 72\Delta\sqrt{\frac{\ln(\omega)}{\ln(\Delta)}}.$$ This implies that for sufficiently large $\Delta$, if $\omega^{(72c)^2} \leq \Delta$ then $\chi(G) \leq \Delta/c$. This bound actually holds for the list-chromatic and even the correspondence-chromatic number (also known as the DP-chromatic number). In fact, we prove what we call a "local version" of it, a result implying the existence of a coloring when the number of available colors for each vertex depends on local parameters, like the degree and the clique number of its neighborhood. Our result simultaneously implies the linear improvement over the "greedy bound" and the two aforementioned results of Johansson.


title: The Binary Dwarf Carbon Star SDSS J125017.90+252427.6

abstract: Although dwarf carbon (dC) stars are thought universally to be binaries to explain the presence of $C_2$ in their spectra while still near main sequence luminosity, direct observational evidence for binarity is remarkably scarce. Here we report the detection of a 2.92 d periodicity in both photometry and radial velocity of SDSS J125017.90+252427.6, an $r=16.4$ dC star. This is the first photometric binary dC, and only the second dC spectroscopic binary. The relative phase of the photometric period to the spectroscopic observations suggests that the photometric variations are a reflection effect due to heating from an unseen companion. The observed radial velocity amplitude of the dC component ($K = 98.8\pm10.7$ km/s) is consistent with a white dwarf companion, presumably the evolved star that earlier donated the carbon to the dC, although substantial orbital evolution must have occurred. Large synoptic photometric surveys such as the Palomar Transient Factory, used for this work, may prove useful for identifying binaries among the shorter period dC stars.


title: The Geometry of Most Probable Trajectories in Noise-Driven Dynamical Systems

abstract: This paper presents a heuristic derivation of a geometric minimum action method that can be used to determine most-probable transition paths in noise-driven dynamical systems. Particular attention is focused on systems that violate detailed balance, and the role of the stochastic vorticity tensor is emphasized. The general method is explored through a detailed study of a two-dimensional quadratic shear flow which exhibits bifurcating most-probable transition pathways.


title: Long lifetime of thermally-excited magnons in bulk yttrium iron garnet

abstract: Spin currents are generated within the bulk of magnetic materials due to heat flow, an effect called intrinsic spin-Seebeck. This bulk bosonic spin current consists of a diffusing thermal magnon cloud, parametrized by the magnon chemical potential ($\mu_{m}$), with a diffusion length of several microns in yttrium iron garnet (YIG). Transient opto-thermal measurements of the spin-Seebeck effect (SSE) as a function of temperature reveal the time evolution of $\mu_{m}$ due to intrinsic SSE in YIG. The interface SSE develops at times < 2 ns while the intrinsic SSE signal continues to evolve at times > 500 $\mu$s, dominating the temperature dependence of SSE in bulk YIG. Time-dependent SSE data are fit to a multi-temperature model of coupled spin/heat transport using finite element method (FEM), where the magnon spin lifetime ($\tau$) and magnon-phonon thermalization time ($\tau_{mp}$) are used as fit parameters. From 300 K to 4 K, $\tau_{mp}$ varies from 1 to 10 ns, whereas $\tau$ varies from 2 to 60 $\mu$s with the spin lifetime peaking at 90 K. At low temperature, a reduction in $\tau$ is observed consistent with impurity relaxation reported in ferromagnetic resonance measurements. These results demonstrate that the thermal magnon cloud in YIG contains extremely low frequency magnons (~10 GHz) providing spectral insight to the microscopic scattering processes involved in magnon spin/heat diffusion.


title: On k-11-representable graphs

abstract: Distinct letters $x$ and $y$ alternate in a word $w$ if after deleting in $w$ all letters but the copies of $x$ and $y$ we either obtain a word of the form $xyxy\cdots$ (of even or odd length) or a word of the form $yxyx\cdots$ (of even or odd length). A simple graph $G=(V,E)$ is word-representable if there exists a word $w$ over the alphabet $V$ such that letters $x$ and $y$ alternate in $w$ if and only if $xy$ is an edge in $E$. Thus, edges of $G$ are defined by avoiding the consecutive pattern 11 in a word representing $G$, that is, by avoiding $xx$ and $yy$. In 2017, Jeff Remmel has introduced the notion of a $k$-$11$-representable graph for a non-negative integer $k$, which generalizes the notion of a word-representable graph. Under this representation, edges of $G$ are defined by containing at most $k$ occurrences of the consecutive pattern $11$ in a word representing $G$. Thus, word-representable graphs are precisely $0$-$11$-representable graphs. Our key result in this paper is showing that any graph is $2$-$11$-representable by a concatenation of permutations, which is rather surprising taking into account that concatenation of permutations has limited power in the case of $0$-$11$-representation. Also, we show that the class of word-representable graphs, studied intensively in the literature, is contained strictly in the class of $1$-$11$-representable graphs. Another result that we prove is the fact that the class of interval graphs is precisely the class of $1$-$11$-representable graphs that can be represented by uniform words containing two copies of each letter. This result can be compared with the known fact that the class of circle graphs is precisely the class of $0$-$11$-representable graphs that can be represented by uniform words containing two copies of each letter.


title: Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination

abstract: Despite the extreme diversity of T cell repertoires, many identical T-cell receptor (TCR) sequences are found in a large number of individual mice and humans. These widely-shared sequences, often referred to as `public', have been suggested to be over-represented due to their potential immune functionality or their ease of generation by V(D)J recombination. Here we show that even for large cohorts the observed degree of sharing of TCR sequences between individuals is well predicted by a model accounting for by the known quantitative statistical biases in the generation process, together with a simple model of thymic selection. Whether a sequence is shared by many individuals is predicted to depend on the number of queried individuals and the sampling depth, as well as on the sequence itself, in agreement with the data. We introduce the degree of publicness conditional on the queried cohort size and the size of the sampled repertoires. Based on these observations we propose a public/private sequence classifier, `PUBLIC' (Public Universal Binary Likelihood Inference Classifier), based on the generation probability, which performs very well even for small cohort sizes.


title: Canonical sphere bundles of the Grassmann manifold

abstract: For a given Hilbert space $\mathcal H$, consider the space of self-adjoint projections $\mathcal P(\mathcal H)$. In this paper we study the differentiable structure of a canonical sphere bundle over $\mathcal P(\mathcal H)$ given by $$ \mathcal R=\{\, (P,f)\in \mathcal P(\mathcal H)\times \mathcal H \, : \, Pf=f , \, \|f\|=1\, \}. $$ We establish the smooth action on $\mathcal R$ of the group of unitary operators of $\mathcal H$, therefore $\mathcal R$ is an homogeneous space. Then we study the metric structure of $\mathcal R$ by endowing it first with the uniform quotient metric, which is a Finsler metric, and we establish minimality results for the geodesics. These are given by certain one-parameter groups of unitary operators, pushed into $\mathcal R$ by the natural action of the unitary group. Then we study the restricted bundle $\mathcal R_2^+$ given by considering only the projections in the restricted Grassmannian, locally modelled by Hilbert-Schmidt operators. Therefore we endow $\mathcal R_2^+$ with a natural Riemannian metric that can be obtained by declaring that the action of the group is a Riemannian submersion. We study the Levi-Civita connection of this metric and establish a Hopf-Rinow theorem for $\mathcal R_2^+$, again obtaining a characterization of the geodesics as the image of certain one-parameter groups with special speeds.


title: Riccati parametric deformations of the Cornu spiral

abstract: A parametric deformation of the Cornu spiral is introduced. The parameter is an integration constant which appears in the general solution of the Riccati equation related to the Fresnel integrals. Argand plots of the deformed spirals are presented and a supersymmetric (Darboux) structure of the deformation is revealed through the factorization approach.


title: Perpetual Orbit Coupled Simulated Annealing for Continuous Optimization

abstract: This paper proposes a novel technique called Perpetual Orbit (PO) to control the generation temperature of the Coupled Simulated Annealing (CSA) algorithm. The new algorithm, PO-CSA, takes advantage of the PO technique combined to the CSA's automatic quasi-optimal scheduling of the acceptance temperature to make optimization of continuous variables more robust w.r.t. the initialization parameters. The PO-CSA was tested with 14 optimization problems, compared to five reference algorithms, including the original CSA with two different initialization procedures. While the PO-CSA performed better than its competitors for the majority of the cost functions, its control of the generation temperature also proved to be more effective than the original CSA with exhaustively tuned initialization parameters. The new parameter-free method also proved more effective in the use of an extra budget of cost function evaluations when compared to the exhaustively tuned case.


title: Fragmentation of charm to charmonium in $e^+e^-$ and $pp$ collisions

abstract: We perform numerical comparison of the fragmentation mechanism of charmonium production ($g\,g\to c\,\bar{c}$ followed by $c\to\psi\,c$) with the full leading order calculation ($g\,g\to\psi\,c\,\bar{c}$ at ${\O} (\alpha_s^4)$). We conclude that the non-fragmentation contributions remain important up to $J/\psi$ transverse momenta about as large as 50 GeV, thus making questionable the applicability of the fragmentation approximation at smaller transverse momenta.