title: Unsupervised Probabilistic Deformation Modeling for Robust Diffeomorphic Registration

abstract: We propose a deformable registration algorithm based on unsupervised learning of a low-dimensional probabilistic parameterization of deformations. We model registration in a probabilistic and generative fashion, by applying a conditional variational autoencoder (CVAE) network. This model enables to also generate normal or pathological deformations of any new image based on the probabilistic latent space. Most recent learning-based registration algorithms use supervised labels or deformation models, that miss important properties such as diffeomorphism and sufficiently regular deformation fields. In this work, we constrain transformations to be diffeomorphic by using a differentiable exponentiation layer with a symmetric loss function. We evaluated our method on 330 cardiac MR sequences and demonstrate robust intra-subject registration results comparable to two state-of-the-art methods but with more regular deformation fields compared to a recent learning-based algorithm. Our method reached a mean DICE score of 78.3% and a mean Hausdorff distance of 7.9mm. In two preliminary experiments, we illustrate the model's abilities to transport pathological deformations to healthy subjects and to cluster five diseases in the unsupervised deformation encoding space with a classification performance of 70%.


title: Efficient Interpolation for the Theory of Arrays

abstract: Existing techniques for Craig interpolation for the quantifier-free fragment of the theory of arrays are inefficient for computing sequence and tree interpolants: the solver needs to run for every partitioning $(A, B)$ of the interpolation problem to avoid creating $AB$-mixed terms. We present a new approach using Proof Tree Preserving Interpolation and an array solver based on Weak Equivalence on Arrays. We give an interpolation algorithm for the lemmas produced by the array solver. The computed interpolants have worst-case exponential size for extensionality lemmas and worst-case quadratic size otherwise. We show that these bounds are strict in the sense that there are lemmas with no smaller interpolants. We implemented the algorithm and show that the produced interpolants are useful to prove memory safety for C programs.


title: Optimizing thermoelectric efficiency of superlattice nanowires at room temperature

abstract: It is known that the figure of merit ($ZT$) of thin nanowires can be significantly enhanced at room temperature due to the reduction of phonon thermal conductance arising from the increase of boundary scattering of phonons. It is expected that the phonon thermal conductance of nanowires filled with quantum dots (QDs) will be further reduced. Here we consider a superlattice nanowire (SLNW) modeled by a linear chain of strongly coupled QDs connected to electrodes. We study the dependence of $ZT$ on the QD energy level ($E_0$) (relative to the Fermi level $E_F$ in the electrodes), inter-dot coupling strength ($t_c$), tunneling rate ($\Gamma$), and temperature $T$ in order to optimize the design. It is found that at room temperature the maximum power factor occurs when $(E_0-E_F)/k_BT \approx 2.4$ and $\Gamma=t_c$, a result almost independent of the number of QDs in SLNW as long as $t_c/k_BT <0.5$. By using reasonable physical parameters we show that thin SLNW with cross-sectional width near $3~nm$ has a potential to achieve $ZT\ge3$.


title: Existence of weak solutions to first-order stationary mean-field games with Dirichlet conditions

abstract: In this paper, we study first-order stationary monotone mean-field games (MFGs) with Dirichlet boundary conditions. While for Hamilton--Jacobi equations Dirichlet conditions may not be satisfied, here, we establish the existence of solutions of MFGs that satisfy those conditions. To construct these solutions, we introduce a monotone regularized problem. Applying Schaefer's fixed-point theorem and using the monotonicity of the MFG, we verify that there exists a unique weak solution to the regularized problem. Finally, we take the limit of the solutions of the regularized problem and using Minty's method, we show the existence of weak solutions to the original MFG.


title: Some implications between Grothendieck's anabelian conjectures

abstract: Grothendieck gave two forms of his "main conjecture of anabelian geometry", i.e. the section conjecture and the hom conjecture. He stated that these two forms are equivalent and that if they hold for hyperbolic curves then they hold for elementary anabelian varieties too. We state a stronger form of Grothendieck's conjecture (equivalent in the case of curves) and prove that Grothendieck's statements hold for our form of the conjecture. We work with DM stacks, rather than schemes. If $X$ is a DM stack over a field $k\subseteq\mathbb{C}$ finitely generated over $\mathbb{Q}$, we prove that whether $X$ satisfies the conjecture or not depends only on $X_{\mathbb{C}}$. We prove that the section conjecture for hyperbolic orbicurves stated by Borne and Emsalem follows from the conjecture for hyperbolic curves.


title: Recognizing Birds from Sound - The 2018 BirdCLEF Baseline System

abstract: Reliable identification of bird species in recorded audio files would be a transformative tool for researchers, conservation biologists, and birders. In recent years, artificial neural networks have greatly improved the detection quality of machine learning systems for bird species recognition. We present a baseline system using convolutional neural networks. We publish our code base as reference for participants in the 2018 LifeCLEF bird identification task and discuss our experiments and potential improvements.


title: Vehicle Communication Strategies for Simulated Highway Driving

abstract: Interest in emergent communication has recently surged in Machine Learning. The focus of this interest has largely been either on investigating the properties of the learned protocol or on utilizing emergent communication to better solve problems that already have a viable solution. Here, we consider self-driving cars coordinating with each other and focus on how communication influences the agents' collective behavior. Our main result is that communication helps (most) with adverse conditions.


title: Data-Driven Analysis of Pareto Set Topology

abstract: When and why can evolutionary multi-objective optimization (EMO) algorithms cover the entire Pareto set? That is a major concern for EMO researchers and practitioners. A recent theoretical study revealed that (roughly speaking) if the Pareto set forms a topological simplex (a curved line, a curved triangle, a curved tetrahedron, etc.), then decomposition-based EMO algorithms can cover the entire Pareto set. Usually, we cannot know the true Pareto set and have to estimate its topology by using the population of EMO algorithms during or after the runtime. This paper presents a data-driven approach to analyze the topology of the Pareto set. We give a theory of how to recognize the topology of the Pareto set from data and implement an algorithm to judge whether the true Pareto set may form a topological simplex or not. Numerical experiments show that the proposed method correctly recognizes the topology of high-dimensional Pareto sets within reasonable population size.


title: Fifth forces, Higgs portals and broken scale invariance

abstract: We study the relationship between the strength of fifth forces and the origin of scale breaking in the Standard Model (SM) of particle physics. We start with a light scalar field that is conformally coupled to a toy SM matter sector through a Weyl rescaling of the metric. After appropriately normalizing the fields, the conformally coupled scalar only interacts directly with the would-be Higgs field through kinetic-mixing and Higgs-portal terms. Thus, for the first time, we describe the equivalence of conformally coupled scalar-tensor modifications of gravity and Higgs-portal theories, and we find that the usual tree-level fifth forces only emerge if there is mass mixing between the conformally coupled scalar and the Higgs field. The strength of the fifth force, mediated by the light scalar, then depends on whether the mass of the Higgs arises from an explicit symmetry-breaking term or a spontaneous mechanism of scale breaking. Solar System tests of gravity and the non-observation of fifth forces therefore have the potential to provide information about the structure of the Higgs sector and the origin of its symmetry breaking, setting an upper bound on the magnitude of any explicit scale-breaking terms. These results demonstrate the phenomenological importance (both for cosmology and high-energy physics) of considering how scalar-tensor modifications of gravity are embedded within extensions of the SM.


title: Beam Selection for MmWave Massive MIMO Systems Under Hybrid Transceiver Architecture

abstract: Hybrid analog/digital precoding is a promising technology to cut down hardware cost in massive multi-input multi-output (MIMO) systems at millimeter wave frequencies. In this letter, we propose a joint beam selection scheme for analog domain precoding under discrete lens array. By considering channel correlation among users, the proposed scheme is able to avoid inter-user interference and maximize system sum rate. Moreover, a bio-inspired ant colony optimization-based algorithm is proposed to obtain a near-optimal solution with dramatically reduced computational complexity. Finally, simulations show the advantages of the proposed scheme in improving system sum rate.


title: Beyond-proximity-force-approximation Casimir force between two spheres at finite temperature II: plasma versus Drude modeling, grounded versus isolated spheres

abstract: A recent experiment [J. L. Garrett, D. A. T. Somers, and J. N. Munday, Phys. Rev. Lett {\bf 120}, 040401 (2018)] measured for the first time the gradient of the Casimir force between two gold spheres in vacuum at room temperature, and placed a bound on the magnitude of the deviation of the measured force from the proximity force approximation (PFA). The present work extends a previous theoretical analysis of this experiment [G. Bimonte, Phys. Rev. D {\bf 97}, 085011 (2018)], by analyzing in detail how the magnitude of the deviation from PFA is affected by the inclusion or neglect of ohmic dissipation at zero frequency, a much debated issue in the current Casimir literature, which goes by the name of the Drude vs plasma controversy. We analyze as well the effect of connecting the plates to charge reservoirs, which is the standard configuration used in Casimir experiments. We describe a simple and effective decimation procedure, allowing for a faster computation of the Casimir force for large aspect ratios of the system.


title: Industrial Symbiotic Networks as Coordinated Games

abstract: We present an approach for implementing a specific form of collaborative industrial practices-called Industrial Symbiotic Networks (ISNs)-as MC-Net cooperative games and address the so called ISN implementation problem. This is, the characteristics of ISNs may lead to inapplicability of fair and stable benefit allocation methods even if the collaboration is a collectively desired one. Inspired by realistic ISN scenarios and the literature on normative multi-agent systems, we consider regulations and normative socioeconomic policies as two elements that in combination with ISN games resolve the situation and result in the concept of coordinated ISNs.


title: Joint MMSE Transceiver Design for Downlink Heterogeneous Network

abstract: In this paper, we propose a minimum mean square error (MMSE) based transceiver design scheme for a downlink multiple-input multiple-output (MIMO) two-tier heterogeneous network with general linear equality per-cell power constraints. Three practical channel models with both perfect and imperfect channel state information are used in simulations. In each channel model, we consider two system configurations, two data transmission schemes and two cellular cooperation scenarios. Our study shows that the proposed MMSE scheme is more flexible than interference alignment (IA) based scheme. For the cases where the IA-type scheme is applicable, the proposed scheme generally outperforms IA-type scheme in terms of average sum rate and bit error rate (BER), but is computationally more complex than the IA-type scheme.


title: Do neutrons publish? A neutron publication survey 2005-2015

abstract: Publication in scientific journals is the main product of scientific research. The amount of papers published, their placement in high impact journals, and their citations are used as a measure of the productivity of individual scientists, institutes or fields of science. To give a profound basis on the publication record and the quality of the publication efforts in neutron scattering, a survey has been done following the approach to use bibliographic databases. Questions to be addressed by this survey are: Is the productivity of research with neutrons changing over the years? Which is the geographic distribution in this field of research? Which ones are leading facilities? Is the quality of publications changing? The main results found are presented.


title: Constraints on axion-like particle properties with very high energy gamma-ray observations of Galactic sources

abstract: Axion-like particles (ALPs) can oscillate to photons and vice versa in electromagnetic fields. The photon-ALP oscillation provides an attractive solution to the apparent transparency of the Universe to TeV photons. The allowed parameter regions for the ALP mass $m_{\rm a}\leq 10^{-7}$ eV have been tightly constrained by the Fermi-LAT and H.E.S.S observations of some extragalactic sources. In this work we show for the first time that the H.E.S.S observations of some TeV sources in the Galactic plane exclude the highest ALP mass region (i.e., $m_{\rm a}\sim {\rm a ~few\times 10^{-7}}$ eV) that accounts for the TeV transparency of the Universe. The upcoming CTA observations of the Galactic TeV sources are shown to be able to improve the constraints significantly.


title: Motion Fused Frames: Data Level Fusion Strategy for Hand Gesture Recognition

abstract: Acquiring spatio-temporal states of an action is the most crucial step for action classification. In this paper, we propose a data level fusion strategy, Motion Fused Frames (MFFs), designed to fuse motion information into static images as better representatives of spatio-temporal states of an action. MFFs can be used as input to any deep learning architecture with very little modification on the network. We evaluate MFFs on hand gesture recognition tasks using three video datasets - Jester, ChaLearn LAP IsoGD and NVIDIA Dynamic Hand Gesture Datasets - which require capturing long-term temporal relations of hand movements. Our approach obtains very competitive performance on Jester and ChaLearn benchmarks with the classification accuracies of 96.28% and 57.4%, respectively, while achieving state-of-the-art performance with 84.7% accuracy on NVIDIA benchmark.


title: Future Decline of African Dust: Insights from the Recent Past and Paleo-records

abstract: African dust exhibits strong variability on a range of time scales. Here we show that the interhemispheric contrast in Atlantic SST (ICAS) drives African dust variability on interannual, multidecadal, and millennial timescales, and a strong anthropogenic decline of African dust in the future can be expected due to the projected increase of the ICAS. During the recent past, the ICAS is found to significantly correlate with various dust observations and proxies that extend as far back as 1851. Physically, positive ICAS anomalies induce large-scale circulation changes that push the Intertropical Convergence Zone (ITCZ) northward and decrease surface wind speed over African dust source regions, which reduces dust emission and transport to the tropical Atlantic. The ICAS-ITCZ-dust relationship also finds robust support from paleo-climate observations that span the last 17,000 years. The ICAS drive of African dust variability is consistent with documented relationships between dust activity and Sahel precipitation6, the North Atlantic Oscillation7, and time series of a surface wind speed pattern over Northern Africa8, and offers a unified framework to understand them. The ICAS-dust connection implies that human activities that change ICAS through emitting greenhouse gases and pollutions have affected and will continue to affect African dust. Climate models project that anthropogenic increase of the ICAS can push the ICAS value to surpass the highest level attained during the Holocene by the end of this century and decrease African dust activity by as much as 60% off its current level, which has broad consequences for aspects of the climate in the North Atlantic region and beyond.


title: Slow relaxation dynamics in binary glasses during stress-controlled, tension-compression cyclic loading

abstract: The effect of cyclic loading on relaxation dynamics and mechanical properties of metallic glasses is studied using molecular dynamics simulations. We consider the Kob-Andersen three-dimensional binary mixture rapidly cooled across the glass transition and subjected to thousands of tension-compression cycles in the elastic range. It was found that during cyclic loading at constant pressure, the system is relocated to progressively lower levels of the potential energy, thus promoting greater densification and higher strength. Furthermore, with increasing stress amplitude, the average glass density increases and the minimum of the potential energy becomes deeper, while the elastic modulus is reduced. The typical size of clusters of atoms with large nonaffine displacements becomes smaller over consecutive cycles, which correlates with the gradual decrease in the potential energy. These results are important for thermomechanical processing of metallic glasses with improved mechanical properties.


title: Reliable IoT Storage: Minimizing Bandwidth Use in Storage Without Newcomer Nodes

abstract: This letter characterizes the optimal policies for bandwidth use and storage for the problem of distributed storage in Internet of Things (IoT) scenarios, where lost nodes cannot be replaced by new nodes as is typically assumed in Data Center and Cloud scenarios. We develop an information flow model that captures the overall process of data transmission between IoT devices, from the initial preparation stage (generating redundancy from the original data) to the different repair stages with fewer and fewer devices. Our numerical results show that in a system with 10 nodes, the proposed optimal scheme can save as much as 10.3% of bandwidth use, and as much as 44% storage use with respect to the closest suboptimal approach.


title: Phenomenological Three-Orbital Spin-Fermion Model for Cuprates

abstract: A spin-fermion model that captures the charge-transfer properties of Cu-based high critical temperature superconductors is introduced and studied via Monte Carlo simulations. The strong Coulomb repulsion among $d$-electrons in the Cu orbitals is phenomenologically replaced by an exchange coupling between the spins of the itinerant electrons and localized spins at the Cu sites, formally similar to double-exchange models for manganites. This interaction induces a charge-transfer insulator gap in the undoped case (five electrons per unit cell). Adding a small antiferromagnetic Heisenberg coupling between localized spins reinforces the global tendency towards antiferromagnetic order. To perform numerical calculations the localized spins are considered classical, as in previous related efforts. In this first study, undoped and doped $8\times 8$ clusters are analyzed in a wide range of temperatures. The numerical results reproduce experimental features in the one-particle spectral function and the density-of-states such as $(i)$ the formation of a Zhang-Rice-like band with a dispersion of order $\sim 0.5$ eV and with rotational symmetry about wavevector $(\pi/2,\pi/2)$ at the top of the band, and $(ii)$ the opening of a pseudogap at the chemical potential upon doping. We also observed incipient tendencies towards spin incommensurability. This simple model offers a formalism intermediate between standard mean-field approximations, that fail at finite temperatures in regimes with short-range order, and sophisticated many-body techniques such as Quantum Monte Carlo, that suffer sign problems.


title: Multiple factor analysis of distributional data

abstract: In the framework of Symbolic Data Analysis (SDA), distribution-variables are a particular case of multi-valued variables: each unit is represented by a set of distributions (e.g. histograms, density functions or quantile functions), one for each variable. Factor analysis (FA) methods are primary exploratory tools for dimension reduction and visualization. In the present work, we use Multiple Factor Analysis (MFA) approach for the analysis of data described by distributional variables. Each distributional variable induces a set new numeric variable related to the quantiles of each distribution. We call these new variables as \textit{quantile variables} and the set of quantile variables related to a distributional one is a block in the MFA approach. Thus, MFA is performed on juxtaposed tables of quantile variables. \\ We show that the criterion decomposed in the analysis is an approximation of the variability based on a suitable metrics between distributions: the squared $L_2$ Wasserstein distance. \\ Applications on simulated and real distributional data corroborate the method. The interpretation of the results on the factorial planes is performed by new interpretative tools that are related to the several characteristics of the distributions (location, scale and shape).


title: Lipschitz Continuity in Model-based Reinforcement Learning

abstract: We examine the impact of learning Lipschitz continuous models in the context of model-based reinforcement learning. We provide a novel bound on multi-step prediction error of Lipschitz models where we quantify the error using the Wasserstein metric. We go on to prove an error bound for the value-function estimate arising from Lipschitz models and show that the estimated value function is itself Lipschitz. We conclude with empirical results that show the benefits of controlling the Lipschitz constant of neural-network models.


title: Juggling with light

abstract: We discovered that when a pair of small particles is optically levitated, the particles execute a dance whose motion resembles the orbits of balls being juggled. This motion lies in a plane perpendicular to the polarization of the incident light. We ascribe the dance to a mechanism by which the dominant force on each particle cyclically alternates between radiation pressure and gravity as each particle takes turns eclipsing the other. We explain the plane of motion by considering the anisotropic scattering of polarized light at a curved interface.

categories: physics.optics cond-mat.soft physics.comp-ph physics.flu-dyn


title: Accelerating Hybrid Monte Carlo simulations of the Hubbard model on the hexagonal lattice

abstract: We present different methods to increase the performance of Hybrid Monte Carlo simulations of the Hubbard model in two-dimensions. Our simulations concentrate on a hexagonal lattice, though can be easily generalized to other lattices. It is found that best results can be achieved using a flexible GMRES solver for matrix inversions and the second order Omelyan integrator with Hasenbusch acceleration on different time scales for molecular dynamics. We demonstrate how an arbitrary number of Hasenbusch mass terms can be included into this geometry and find that the optimal speed depends weakly on the choice of the number of Hasenbusch masses and their values. As such, the tuning of these masses is amenable to automization and we present an algorithm for this tuning that is based on the knowledge of the dependence of solver time and forces on the Hasenbusch masses. We benchmark our algorithms to systems where direct numerical diagonalization is feasible and find excellent agreement. We also simulate systems with hexagonal lattice dimensions up to $102\times 102$ and $N_t=64$. We find that the Hasenbusch algorithm leads to a speed up of more than an order of magnitude.


title: Equilibrium magnetization of a quasispherical cluster of single-domain particles

abstract: Equilibrium magnetization curve of a rigid finite-size spherical cluster of single-domain particles is investigated both numerically and analytically. The spatial distribution of particles within the cluster is random. Dipole-dipole interactions between particles are taken into account. The particles are monodisperse. It is shown, using the stochastic Landau-Lifshitz-Gilbert equation that the magnetization of such clusters is generally lower than predicted by the classical Langevin model. In a broad range of dipolar coupling parameters and particle volume fractions, the cluster magnetization in the weak field limit can be successfully described by the modified mean-field theory, which was originally proposed for the description of concentrated ferrofluids. In moderate and strong fields, the theory overestimates the cluster magnetization. However, predictions of the theory can be improved by adjusting the corresponding mean-field parameter. If magnetic anisotropy of particles is additionally taken into account and if the distribution of the particles' easy axes is random and uniform, then the cluster equilibrium response is even weaker. The decrease of the magnetization with increasing anisotropy constant is more pronounced at large applied fields. The phenomenological generalization of the modified mean-field theory, that correctly describes this effect for small coupling parameters, is proposed.


title: Spectral estimates for Dirichlet Laplacian on tubes with exploding twisting velocity

abstract: We study the spectrum of the Dirichlet Laplacian on an unbounded twisted tube with twisting velocity exploding to infinity. If the tube cross section does not intersect the axis of rotation, then its spectrum is purely discrete under some additional conditions on the twisting velocity (D.Krejcirik, 2015). In the current work we prove a Berezin type upper bound for the eigenvalue moments.


title: Breaking Chiral Symmetry with Microfluidics

abstract: A robust route for the biased production of single-handed chiral structures has been found in generating non-spherical, multi-component double emulsions using microfluidics. The specific type of handedness is determined by the final packing geometry of four different inner drops inside an ultra-thin sheath of oil. Before three-dimensional chiral structures are formed, the quasi-one-dimensional chain re-arranges in two dimensions into either checkerboard or stripe patterns. We derive an analytical model predicting which pattern is more likely and assembles in the least amount of time. Moreover, our dimensionless model accurately predicts our experimental results and is based on local bending dynamics, rather than global surface energy minimization. This better reflects the underlying self-assembly process which will not, in general, reach a global energy minimum. In summary, using glass microfluidic techniques for channeling aqueous fluids through narrow orifices of multi-bore injection capillaries while encapsulating these fluids as drops inside an ultra-thin sheath of oil is sufficient to produce single-handed chiral structures.


title: Probing Ferromagnetic Order in Few-Fermion Correlated Spin-Flip Dynamics

abstract: We unravel the dynamical stability of a fully polarized one-dimensional ultracold few-fermion spin-1/2 gas subjected to inhomogeneous driving of the itinerant spins. Despite the unstable character of the total spin-polarization the existence of an interaction regime is demonstrated where the spin-correlations lead to almost maximally aligned spins throughout the dynamics. The resulting ferromagnetic order emerges from the build up of superpositions of states of maximal total spin. They comprise a decaying spin-polarization and a dynamical evolution towards an almost completely unpolarized NOON-like state. Via single-shot simulations we demonstrate that our theoretical predictions can be detected in state-of-the-art ultracold experiments.


title: Silence measurements and measures for ET: characterisation of long term seismic noise in the M\'atra Mountains

abstract: The analysis of long term seismological data collected underground in the M\'atra Mountains, Hungary, using the facilities of the M\'atra Gravitational and Geophysical Laboratory (MGGL) is reported. The laboratory is situated inside the Gy\"ongy\"osoroszi mine, Hungary, 88m below the surface. This study focuses on the requirements of the Einstein Telescope (ET), one of the planned third generation gravitational wave observatories, which is designed for underground operation. After a short introduction of the geophysical environment the evaluation of the collected long term data follows including the comparison of a two-week measurement campaign deeper in the mine. Based on our analysis and considering the specialities of long term data collection, refinements of the performance and evaluation criteria are suggested as well as performance estimation of a possible M\'atra site.


title: Anonymous Single-Sign-On for n designated services with traceability

abstract: Anonymous Single-Sign-On authentication schemes have been proposed to allow users to access a service protected by a verifier without revealing their identity which has become more important due to the introduction of strong privacy regulations. In this paper we describe a new approach whereby anonymous authentication to different verifiers is achieved via authorisation tags and pseudonyms. The particular innovation of our scheme is authentication can only occur between a user and its designated verifier for a service, and the verification cannot be performed by any other verifier. The benefit of this authentication approach is that it prevents information leakage of a user's service access information, even if the verifiers for these services collude which each other. Our scheme also supports a trusted third party who is authorised to de-anonymise the user and reveal her whole services access information if required. Furthermore, our scheme is lightweight because it does not rely on attribute or policy-based signature schemes to enable access to multiple services. The scheme's security model is given together with a security proof, an implementation and a performance evaluation.


title: Wiggling skyrmion propagation under parametric pumping

abstract: We address the problem of how magnetic skyrmions can propagate along a guided direction by parametric pumping. As evidenced by our micromagnetic simulations, skyrmions can hardly be driven by either a static electric field or a static magnetic field alone. Although the magnetic anisotropy can be modified by an electric field, parametric pumping with an oscillating electric field can only excite the breathing modes. On the other hand, a static magnetic field can break rotational symmetry through the Zeeman interaction, but it cannot serve as an energy source for propelling a skyrmion. Here we found that the combination of a perpendicularly oscillating electric field and an in-plane static magnetic field can drive a skyrmion undergoing a wiggling motion along a well defined trajectory. The most efficient driving occurs when the frequency of the oscillating field is close to that of the breathing motion. The physics is revealed in a generalized Thiele equation where a net spin current excited by the parametric pumping can drive the skyrmion propagation through angular momentum transfer. Compared with other alternative proposals, our results open new possibilities for manipulating skyrmions in both metals and insulators with low-power consumption. The oscillating skyrmion motion can also be a microwave generator for future spintronic applications such as an nano-tool on a diamond Nitrogen-Vacancy center.


title: The Hardness of Conditional Independence Testing and the Generalised Covariance Measure

abstract: It is a common saying that testing for conditional independence, i.e., testing whether whether two random vectors $X$ and $Y$ are independent, given $Z$, is a hard statistical problem if $Z$ is a continuous random variable (or vector). In this paper, we prove that conditional independence is indeed a particularly difficult hypothesis to test for. Valid statistical tests are required to have a size that is smaller than a predefined significance level, and different tests usually have power against a different class of alternatives. We prove that a valid test for conditional independence does not have power against any alternative. Given the non-existence of a uniformly valid conditional independence test, we argue that tests must be designed so their suitability for a particular problem may be judged easily. To address this need, we propose in the case where $X$ and $Y$ are univariate to nonlinearly regress $X$ on $Z$, and $Y$ on $Z$ and then compute a test statistic based on the sample covariance between the residuals, which we call the generalised covariance measure (GCM). We prove that validity of this form of test relies almost entirely on the weak requirement that the regression procedures are able to estimate the conditional means $X$ given $Z$, and $Y$ given $Z$, at a slow rate. We extend the methodology to handle settings where $X$ and $Y$ may be multivariate or even high-dimensional. While our general procedure can be tailored to the setting at hand by combining it with any regression technique, we develop the theoretical guarantees for kernel ridge regression. A simulation study shows that the test based on GCM is competitive with state of the art conditional independence tests. Code is available as the R package GeneralisedCovarianceMeasure on CRAN.


title: Self-enforcing Game Theory-based Resource Allocation for LoRaWAN Assisted Public Safety Communications

abstract: Public safety networks avail to disseminate information during emergency situations through its dedicated servers. Public safety networks accommodate public safety communication (PSC) applications to track the location of its utilizers and enable to sustain transmissions even in the crucial scenarios. Despite that, if the traditional setups responsible for PSCs are unavailable, it becomes prodigiously arduous to handle any of the safety applications, which may cause havoc in the society. Dependence on a secondary network may assist to solve such an issue. But, the secondary networks should be facilely deployable and must not cause exorbitant overheads in terms of cost and operation. For this, LoRaWAN can be considered as an ideal solution as it provides low power and long-range communication. However, an excessive utilization of the secondary network may result in high depletion of its own resources and can lead to a complete shutdown of services, which is a quandary at hand. As a solution, this paper proposes a novel network model via a combination of LoRaWAN and traditional public safety networks, and uses a self-enforcing agreement based game theory for allocating resources efficiently amongst the available servers. The proposed approach adopts memory and energy constraints as agreements, which are satisfied through Nash equilibrium. The numerical results show that the proposed approach is capable of efficiently allocating the resources with sufficiently high gains for resource conservation, network sustainability, resource restorations and probability to continue at the present conditions even in the complete absence of traditional Access Points (APs) compared with a baseline scenario with no failure of nodes.


title: Statistical features of rapidly rotating decaying turbulence: enstrophy and energy spectra, and coherent structures

abstract: In this paper we investigate the properties of rapidly rotating decaying turbulence using numerical simulations and phenomenological modelling. We find that as the turbulent flow evolves in time, the Rossby number decreases to $\sim 10^{-3}$, and the flow becomes quasi-two-dimensional with strong coherent columnar structures arising due to the inverse cascade of energy. We establish that a major fraction of energy is confined in Fourier modes $(\pm1,0,0)$ and $(0,\pm1,0)$ that correspond to the largest columnar structure in the flow. For wavenumbers ($k$) greater than the enstrophy dissipation wavenumber ($k_d$), our phenomenological arguments and numerical study show that the enstrophy flux and spectrum of a horizontal cross-section perpendicular to the axis of rotation are given by $\epsilon_\omega\exp(-C(k/k_d)^2)$ and $C\epsilon_\omega^{2/3}k^{-1}\exp(-C(k/k_d)^2)$ respectively; for this 2D flow, $\epsilon_\omega$ is the enstrophy dissipation rate, and $C$ is a constant. Using these results, we propose a new form for the energy spectrum of rapidly rotating decaying turbulence: $E(k)=C\epsilon_\omega^{2/3}k^{-3}\exp(-C(k/k_d)^2)$. This model of the energy spectrum is based on wavenumber-dependent enstrophy flux, and it deviates significantly from power law energy spectrum reported earlier.


title: Functional bosonization of a Dirac field in 2+1 dimensions, in the presence of a boundary

abstract: We apply the functional bosonization procedure to a massive Dirac field defined on a 2+1 dimensional spacetime which has a non-trivial boundary. We find the form of the bosonized current both for the bulk and boundary modes, showing that the gauge field in the bosonized theory contains a perfect-conductor boundary condition on the worldsheet spanned by the boundary. We find the bononized action for the corresponding boundary modes.


title: The Effect of Mechanical Cold Work on the Magnetic Flux Expulsion of Niobium

abstract: Expulsion of ambient flux has been shown to be crucial to obtain high quality factors in bulk niobium SRF cavities. However, there remain many questions as to what properties of the niobium material determine its flux expulsion behavior. In this paper, we present first results from a new study of two cavities that were specially fabricated to study flux expulsion. Both cavities were made from large grain ingot niobium slices, one of which had its slices rolled prior to fabrication, and none these slices were annealed prior to measurement. Expulsion measurements indicate that a dense network of grain boundaries is not necessary for a cavity to have near-complete flux trapping behavior up to large thermal gradients. The results also contribute to a body of evidence that cold work is a strong determinant of flux expulsion behavior in SRF-grade niobium.


title: A stochastic model for the evolution of species with random fitness

abstract: We generalize the evolution model introduced by Guiol, Machado and Schinazi (2010). In our model at odd times a random number X of species is created. Each species is endowed with a random fitness with arbitrary distribution on $[0, 1]$. At even times a random number Y of species is removed, killing the species with lower fitness. We show that there is a critical fitness $f_c$ below which the number of species hits zero i.o. and above of which this number goes to infinity. We prove uniform convergence for the distribution of surviving species and describe the phenomena which could not be observed in previous works with uniformly distributed fitness.


title: NAIS-Net: Stable Deep Networks from Non-Autonomous Differential Equations

abstract: This paper introduces Non-Autonomous Input-Output Stable Network(NAIS-Net), a very deep architecture where each stacked processing block is derived from a time-invariant non-autonomous dynamical system. Non-autonomy is implemented by skip connections from the block input to each of the unrolled processing stages and allows stability to be enforced so that blocks can be unrolled adaptively to a pattern-dependent processing depth. NAIS-Net induces non-trivial, Lipschitz input-output maps, even for an infinite unroll length. We prove that the network is globally asymptotically stable so that for every initial condition there is exactly one input-dependent equilibrium assuming $tanh$ units, and incrementally stable for ReL units. An efficient implementation that enforces the stability under derived conditions for both fully-connected and convolutional layers is also presented. Experimental results show how NAIS-Net exhibits stability in practice, yielding a significant reduction in generalization gap compared to ResNets.


title: Investigating Negative Interactions in Multiplex Networks: A Mutual Information Approach

abstract: Many interesting real-world systems are represented as complex networks with multiple types of interactions and complicated dependency structures between layers. These interactions can be encoded as having a valence with positive links marking interactions such as trust and friendship and negative links denoting distrust or hostility. Extracting information from these negative interactions is challenging since standard topological metrics are often poor predictors of negative link formation, particularly across network layers. In this paper, we introduce a method based on mutual information which enables us to predict both negative and positive relationships. Our experiments show that SMLP (Signed Multiplex Link Prediction) can leverage negative relationship layers in multiplex networks to improve link prediction performance.


title: Nucleon Mass Sum Rule Violation in QCD and Confinement

abstract: In this paper we show that for a nucleon at rest if the quarks and/or antiquarks are in motion inside the nucleon producing chromo-magnetic field then the mass sum rule in QCD is violated when the confinement potential energy at large distance $r$ rises linearly with $r$ (or faster). Hence we find that the mass of the nucleon at rest is not equal to the mass-energy of all the quarks plus antiquarks plus gluons inside the nucleon if there exists chromo-magnetic field inside the nucleon and the confinement potential energy at large distance $r$ rises linearly with $r$ (or faster).


title: Learning Disentangled Representations of Texts with Application to Biomedical Abstracts

abstract: We propose a method for learning disentangled representations of texts that code for distinct and complementary aspects, with the aim of affording efficient model transfer and interpretability. To induce disentangled embeddings, we propose an adversarial objective based on the (dis)similarity between triplets of documents with respect to specific aspects. Our motivating application is embedding biomedical abstracts describing clinical trials in a manner that disentangles the populations, interventions, and outcomes in a given trial. We show that our method learns representations that encode these clinically salient aspects, and that these can be effectively used to perform aspect-specific retrieval. We demonstrate that the approach generalizes beyond our motivating application in experiments on two multi-aspect review corpora.


title: A refined convergence analysis of pDCA$_e$ with applications to simultaneous sparse recovery and outlier detection

abstract: We consider the problem of minimizing a difference-of-convex (DC) function, which can be written as the sum of a smooth convex function with Lipschitz gradient, a proper closed convex function and a continuous possibly nonsmooth concave function. We refine the convergence analysis in [38] for the proximal DC algorithm with extrapolation (pDCA$_e$) and show that the whole sequence generated by the algorithm is convergent when the objective is level-bounded, {\em without} imposing differentiability assumptions in the concave part. Our analysis is based on a new potential function and we assume such a function is a Kurdyka-{\L}ojasiewicz (KL) function. We also establish a relationship between our KL assumption and the one used in [38]. Finally, we demonstrate how the pDCA$_e$ can be applied to a class of simultaneous sparse recovery and outlier detection problems arising from robust compressed sensing in signal processing and least trimmed squares regression in statistics. Specifically, we show that the objectives of these problems can be written as level-bounded DC functions whose concave parts are {\em typically nonsmooth}. Moreover, for a large class of loss functions and regularizers, the KL exponent of the corresponding potential function are shown to be 1/2, which implies that the pDCA$_e$ is locally linearly convergent when applied to these problems. Our numerical experiments show that the pDCA$_e$ usually outperforms the proximal DC algorithm with nonmonotone linesearch [24, Appendix A] in both CPU time and solution quality for this particular application.


title: An Action for the Infrared Regime of Gauge Theories and the Problem of Color Transformations

abstract: It has been known for a while that there is spontaneous breaking of Lorentz symmetry in the nonzero charged sectors of quantum electrodynamics due to the infrared problem of soft photons. More recently, it has also been suggested that similar results hold for color transformations in a nonabelian gauge theory. Here we show that an action where a diffeomorphism has been carried out for the part describing hard gauge particles and matter fields can be used to analyze these issues. In addition to rederiving old results in this formalism, we also show that color transformations cannot be unitarily implemented on perturbative gluon states if gluon fields of arbitrarily low energy are allowed. Implications for confinement and mass gap are briefly commented upon.


title: On the mechanical properties of novamene: a fully atomistic molecular dynamics and DFT investigation

abstract: We have investigated through fully atomistic reactive molecular dynamics and DFT simulations, the mechanical properties and fracture dynamics of novamene, a new 3D carbon allotrope structure recently proposed. Our results showed that novamene is an anisotropic structure with relation to tensile deformation. Although novamente shares some mechanical features with other carbon allotropes, it also exhibits distinct ones, such as, extensive structural reconstructions (self-healing effect). Novamene presents ultimate strength (~ 100 GPa) values lower than other carbon allotropes, but it has the highest ultimate strain along the z-direction (~ 22.5%). Although the Young's modulus (~ 600 GPa) and ultimate strength values are smaller than for other carbon allotropes, they still outperform other materials, such as for example silicon, steel or titanium alloys. With relation to the fracture dynamics, novamene is again anisotropic with the fracture/crack propagation originating from deformed heptagons and pentagons for x and y directions and broken sp3 bonds connecting structural planes. Another interesting feature is the formation of multiple and long carbon linear chains in the final fracture stages.


title: The Protons in Primary Cosmic Rays in the Energy Range 10^15-10^17 eV According to Data from the PAMIR Experiment

abstract: Adjusted data on the fraction of protons in the mass composition of primary cosmic rays (PCRs) in the energy range of 10^15-10^17 eV are presented. Adjustments are made according to detailed calculations of the response of the X-ray emulsion chamber in the PAMIR experiment. It is demonstrated that the fraction of protons in a PCR is 16-18% for E0 = 10^15-10^16 eV and does not change within the error for E0 = 10^16-10^17 eV.

categories: astro-ph.HE hep-ex physics.ins-det physics.space-ph


title: The GALAH Survey: Velocity fluctuations in the Milky Way using red clump giants

abstract: If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.


title: Quantum phase transitions in highly crystalline two-dimensional superconductors

abstract: Superconductor-insulator transition is one of the remarkable phenomena driven by quantum fluctuation in two-dimensional (2D) systems. Such a quantum phase transition (QPT) was investigated predominantly on highly disordered thin films with amorphous or granular structures using scaling law with constant exponents. Here, we provide a totally different view of QPT in highly crystalline 2D superconductors. According to the magneto-transport measurements in 2D superconducting ZrNCl and MoS2, we found that the quantum metallic state commonly observed at low magnetic fields is converted via the quantum Griffiths state to the weakly localized metal at high magnetic fields. The scaling behavior, characterized by the diverging dynamical critical exponent (Griffiths singularity), indicates that the quantum fluctuation manifests itself as superconducting puddles, in marked contrast with the thermal fluctuation. We suggest that an evolution from the quantum metallic to the quantum Griffiths state is generic nature in highly crystalline 2D superconductors with weak pinning potentials.


title: Loschmidt echo driven by hyperfine and electric-quadrupole interactions in nanoscale nuclear spin baths

abstract: The nuclear spin bath (NSB) dynamics and its quantum control are of importance for the storage and processing of quantum information within a semiconductor environment. In the presence of a carrier spin, primarily it is the hyperfine interaction that rules the high frequency NSB characteristics. Here, we first study the overall coherence decay and rephasings in a hyperfine-driven NSB through the temporal and spectral behaviors of the so-called Loschmidt echo (LE). Its dependence on the NSB size, initial polarization, and coupling inhomogeneity are separately investigated, which leads to a simple phenomenological expression that can accommodate all of these attributes. Unlike the prevailing emphasis on spin 1/2, the NSBs with larger spin quantum numbers are equally considered. For this case, additionally the effect of nuclear electric quadrupole interaction is taken into account where its biaxiality term is influential on the decoherence. The insights gained from model systems are then put to use for two generic realistic semiconductor systems, namely, a donor center and a quantum dot that represent small and large nanoscale NSB examples, respectively. The spectrum of LE for large quantum dots can reach the 100~MHz range, whereas, for donor centers, it reduces to a few MHz, making them readily amenable for dynamical decoupling techniques. The effect of quadrupole interaction on LE is seen to be negligible for large quantum dots, while it becomes significant for donor centers, most notably in the form of depolarizing a polarized NSB.


title: Randomized opinion dynamics over networks: influence estimation from partial observations

abstract: In this paper, we propose a technique for the estimation of the influence matrix in a sparse social network, in which $n$ individual communicate in a gossip way. At each step, a random subset of the social actors is active and interacts with randomly chosen neighbors. The opinions evolve according to a Friedkin and Johnsen mechanism, in which the individuals updates their belief to a convex combination of their current belief, the belief of the agents they interact with, and their initial belief, or prejudice. Leveraging recent results of estimation of vector autoregressive processes, we reconstruct the social network topology and the strength of the interconnections starting from \textit{partial observations} of the interactions, thus removing one of the main drawbacks of finite horizon techniques. The effectiveness of the proposed method is shown on randomly generation networks.


title: A Decomposition-based Approach towards the Control of Boolean Networks (Technical Report)

abstract: We study the problem of computing a minimal subset of nodes of a given asynchronous Boolean network that need to be controlled to drive its dynamics from an initial steady state (or attractor) to a target steady state. Due to the phenomenon of state-space explosion, a simple global approach that performs computations on the entire network, may not scale well for large networks. We believe that efficient algorithms for such networks must exploit the structure of the networks together with their dynamics. Taking such an approach, we derive a decomposition-based solution to the minimal control problem which can be significantly faster than the existing approaches on large networks. We apply our solution to both real-life biological networks and randomly generated networks, demonstrating promising results.


title: Eclipsing damped Ly$\alpha$ systems in the Sloan Digital Sky Survey Data Release 12

abstract: We present the results of our automatic search for proximate damped Ly$\alpha$ absorption (PDLA) systems in the quasar spectra from the Sloan Digital Sky Survey Data Release 12. We constrain our search to those PDLAs lying within 1500 km s$^{-1}$ from the quasar to make sure that the broad DLA absorption trough masks most of the strong Ly$\alpha$ emission from the broad line region (BLR) of the quasar. When the Ly$\alpha$ emission from the BLR is blocked by these so-called eclipsing DLAs, narrow Ly$\alpha$ emission from the host galaxy could be revealed as a narrow emission line (NEL) in the DLA trough. We define a statistical sample of 399 eclipsing DLAs with log$N$(HI)$\ge$21.10. We divide our statistical sample into three subsamples based on the strength of the NEL detected in the DLA trough. By studying the stacked spectra of these subsamples, we found that absorption from high ionization species are stronger in DLAs with stronger NEL in their absorption core. Moreover, absorption from the excited states of species like SiII are also stronger in DLAs with stronger NEL. We also found no correlation between the luminosity of the Ly$\alpha$ NEL and the quasar luminosity. These observations are consistent with a scenario in which the DLAs with stronger NEL are denser and physically closer to the quasar. We propose that these eclipsing DLAs could be the product of the interaction between infalling and outflowing gas. High resolution spectroscopic observation would be needed to shed some light on the nature of these eclipsing DLAs.


title: Phase Transition of the 2-Choices Dynamics on Core-Periphery Networks

abstract: Consider the following process on a network: Each agent initially holds either opinion blue or red; then, in each round, each agent looks at two random neighbors and, if the two have the same opinion, the agent adopts it. This process is known as the 2-Choices dynamics and is arguably the most basic non-trivial opinion dynamics modeling voting behavior on social networks. Despite its apparent simplicity, 2-Choices has been analytically characterized only on restricted network classes---under assumptions on the initial configuration that establish it as a fast majority consensus protocol. In this work, we aim at contributing to the understanding of the 2-Choices dynamics by considering its behavior on a class of networks with core-periphery structure, a well-known topological assumption in social networks. In a nutshell, assume that a densely-connected subset of agents, the core, holds a different opinion from the rest of the network, the periphery. Then, depending on the strength of the cut between the core and the periphery, a phase-transition phenomenon occurs: Either the core's opinion rapidly spreads among the rest of the network, or a metastability phase takes place, in which both opinions coexist in the network for superpolynomial time. The interest of our result is twofold. On the one hand, by looking at the 2-Choices dynamics as a simplistic model of competition among opinions in social networks, our theorem sheds light on the influence of the core on the rest of the network, as a function of the core's connectivity toward the latter. On the other hand, we provide one of the first analytical results which shows a heterogeneous behavior of a simple dynamics as a function of structural parameters of the network. Finally, we validate our theoretical predictions with extensive experiments on real networks.


title: Morse-Smale systems without heteroclinic submanifolds on codimension one separatrices

abstract: We study a topological structure of a closed $n$-manifold $M^n$ ($n\geq 3$) which admits a Morse-Smale diffeomorphism such that codimension one separatrices of saddles periodic points have no heteroclinic intersections different from heteroclinic points. Also we consider gradient like flow on $M^n$ such that codimension one separatices of saddle singularities have no intersection at all. We show that $M^n$ is either an $n$-sphere $S^n$, or the connected sum of a finite number of copies of $S^{n-1}\otimes S^1$ and a finite number of special manifolds $N^n_i$ admitting polar Morse-Smale systems. Moreover, if some $N^n_i$ contains a single saddle, then $N^n_i$ is projective-like (in particular, $n\in\{4,8,16\}$, and $N^n_i$ is a simply-connected and orientable manifold). Given input dynamical data, one constructs a supporting manifold $M^n$. We give a formula relating the number of sinks, sources and saddle periodic points to the connected sum for $M^n$. As a consequence, we obtain conditions for the existence of heteroclinic intersections for Morse-Smale diffeomorphisms and a periodic trajectory for Morse-Smale flows.


title: Examples of genuine QM abelian surfaces which are modular

abstract: Let $K$ be an imaginary quadratic field. Modular forms for GL(2) over $K$ are known as Bianchi modular forms. Standard modularity conjectures assert that every weight 2 rational Bianchi newform has either an associated elliptic curve over $K$ or an associated abelian surface with quaternionic multiplication over $K$. We give explicit evidence in the way of examples to support this conjecture in the latter case. Furthermore, the quaternionic surfaces given correspond to genuine Bianchi newforms, which answers a question posed by J. Cremona as to whether this phenomenon can happen.


title: Resonance energy transfer between two atoms in a conducting cylindrical waveguide

abstract: We consider the energy transfer process between two identical atoms placed inside a perfectly conducting cylindrical waveguide. We first introduce a general analytical expression of the energy transfer amplitude in terms of the electromagnetic Green's tensor; we then evaluate it in the case of a cylindrical waveguide made of a perfect conductor, for which analytical forms of the Green's tensor exist. We numerically analyse the energy transfer amplitude when the radius of the waveguide is such that the transition frequency of both atoms is below the lower cutoff frequency of the waveguide, so that the resonant photon exchange is strongly suppressed. We consider both cases of atomic dipoles parallel and orthogonal to the axis of the guide. In both cases, we find that the energy transfer is modified by the presence of the waveguide. In the near zone, that is when the atomic separation is smaller than the atomic transition wavelength, the change, with respect to the free-space case, is small for axial dipoles, while it is larger for radial dipoles; it grows when the intermediate region between near and far zone is approached. In the far zone, we find that the energy transfer amplitude is strongly suppressed by the waveguide, becoming virtually zero. A physical interpretation of these results is discussed. Finally, we discuss the resonance interaction energy and force between two identical correlated atoms in the waveguide, one excited and the other in the ground state, prepared in their symmetric or antisymmetric superposition.


title: A $\mathrm{G}_2$-period of a Fourier coefficient of an Eisenstein series on $\mathrm{E}_6$

abstract: We calculate a $\mathrm{G}_2$-period of a Fourier coefficient of a cuspidal Eisenstein series on the split simply-connected group $\mathrm{E}_6$, and relate this period to the Ginzburg-Rallis period of cusp forms on $\mathrm{GL}_6$. This gives us a relation between the Ginzburg-Rallis period and the central value of the exterior cube L-function of $\mathrm{GL}_6$


title: Ferromagnetic phases due to competing short- and long-range interactions in the spin-S Blume-Capel model

abstract: We broaden the study of the statistical physics of the spin-$S$ Blume-Capel model with ferromagnetic mean-field interactions $J$ in competition with short-range antiferromagnetic interactions $K$ in a linear chain in the thermodynamic limit. This work describes the critical behavior of the model when the $S$ takes a half integer and an integer value. In both cases the phase diagrams exhibit new ferromagnetic phases (for certain values of $K$) enclosed by branches emerging from the first-order frontiers of the pure ferromagnetic model. For finite temperatures the complex topologies were obtained by numerical minimization of the free energy.


title: Identification of meteorite source regions in the Solar System

abstract: Over the past decade there has been a large increase in the number of automated camera networks that monitor the sky for fireballs. One of the goals of these networks is to provide the necessary information for linking meteorites to their pre-impact, heliocentric orbits and ultimately to their source regions in the solar system. We re-compute heliocentric orbits for the 25 meteorite falls published to date from original data sources. Using these orbits, we constrain their most likely escape routes from the main asteroid belt and the cometary region by utilizing a state-of-the-art orbit model of the near-Earth-object population, which includes a size-dependence in delivery efficiency. While we find that our general results for escape routes are comparable to previous work, the role of trajectory measurement uncertainty in escape-route identification is explored for the first time. Moreover, our improved size-dependent delivery model substantially changes likely escape routes for several meteorite falls, most notably Tagish Lake which seems unlikely to have originated in the outer main belt as previously suggested. We find that reducing the uncertainty of fireball velocity measurements below $\sim0.1$ km/s does not lead to reduced uncertainties in the identification of their escape routes from the asteroid belt and, further, their ultimate source regions. This analysis suggests that camera networks should be optimized for the largest possible number of meteorite recoveries with measured speed precisions of order 0.1 km/s.


title: Uniqueness of two-convex closed ancient solutions to the mean curvature flow

abstract: In this paper we consider closed non-collapsed ancient solutions to the mean curvature flow ($n \ge 2$) which are uniformly two-convex. We prove that any two such ancient solutions are the same up to translations and scaling. In particular, they must coincide up to translations and scaling with the rotationally symmetric closed ancient non-collapsed solution constructed by Brian White in (2000), and by Robert Haslhofer and Or Hershkovits in (2016).


title: Stationarily ordered types and the number of countable models

abstract: We introduce notions of stationarily ordered types and theories; the latter generalizes weak o-minimality and the first is a relaxed version of weak o-minimality localized at the locus of a single type. We show that forking, as a binary relation on elements realizing stationarily ordered types, is an equivalence relation and that each stationarily ordered type in a model determines some order-type as an invariant of the model. We study weak and forking non-orthogonality of stationarily ordered types, show that they are equivalence relations and prove that invariants of non-orthogonal types are closely related. The developed techniques are applied to prove that in the case of a binary, stationarily ordered theory with fewer than $2^{\aleph_0}$ countable models, the isomorphism type of a countable model is determined by a certain sequence of invariants of the model. In particular, we confirm Vaught's conjecture for binary, stationarily ordered theories.


title: A third strike against perfect phylogeny

abstract: Perfect phylogenies are fundamental in the study of evolutionary trees because they capture the situation when each evolutionary trait emerges only once in history; if such events are believed to be rare, then by Occam's Razor such parsimonious trees are preferable as a hypothesis of evolution. A classical result states that 2-state characters permit a perfect phylogeny precisely if each subset of 2 characters permits one. More recently, it was shown that for 3-state characters the same property holds but for size-3 subsets. A long-standing open problem asked whether such a constant exists for each number of states. More precisely, it has been conjectured that for any fixed integer $r$, there exists a constant $f(r)$ such that a set of $r$-state characters $C$ has a perfect phylogeny if and only if every subset of at most $f(r)$ characters has a perfect phylogeny. In this paper, we show that this conjecture is false. In particular, we show that for any constant $t$, there exists a set $C$ of $8$-state characters such that $C$ has no perfect phylogeny, but there exists a perfect phylogeny for every subset of $t$ characters. This negative result complements the two negative results ("strikes") of Bodlaender et al. We reflect on the consequences of this third strike, pointing out that while it does close off some routes for efficient algorithm development, many others remain open.


title: Evaluation of IEEE 802.11ad for mmWave V2V Communications

abstract: Autonomous vehicles can construct a more accurate perception of their surrounding environment by exchanging rich sensor data with nearby vehicles. Such exchange can require larger bandwidths than currently provided by ITS-G5/DSRC and Cellular V2X. Millimeter wave (mmWave) communications can provide higher bandwidth and could complement current V2X standards. Recent studies have started investigating the potential of IEEE 802.11ad to support high bandwidth vehicular communications. This paper introduces the first performance evaluation of the IEEE 802.11ad MAC (Medium Access Control) and beamforming mechanism for mmWave V2V communications. The study highlights existing opportunities and shortcomings that should guide the development of mmWave communications for V2V communications.


title: Limited Evaluation Cooperative Co-evolutionary Differential Evolution for Large-scale Neuroevolution

abstract: Many real-world control and classification tasks involve a large number of features. When artificial neural networks (ANNs) are used for modeling these tasks, the network architectures tend to be large. Neuroevolution is an effective approach for optimizing ANNs; however, there are two bottlenecks that make their application challenging in case of high-dimensional networks using direct encoding. First, classic evolutionary algorithms tend not to scale well for searching large parameter spaces; second, the network evaluation over a large number of training instances is in general time-consuming. In this work, we propose an approach called the Limited Evaluation Cooperative Co-evolutionary Differential Evolution algorithm (LECCDE) to optimize high-dimensional ANNs. The proposed method aims to optimize the pre-synaptic weights of each post-synaptic neuron in different subpopulations using a Cooperative Co-evolutionary Differential Evolution algorithm, and employs a limited evaluation scheme where fitness evaluation is performed on a relatively small number of training instances based on fitness inheritance. We test LECCDE on three datasets with various sizes, and our results show that cooperative co-evolution significantly improves the test error comparing to standard Differential Evolution, while the limited evaluation scheme facilitates a significant reduction in computing time.


title: Correlation between tides and seismicity in Northwestern South America: the case of Colombia

abstract: We present the first systematic exploration of earth tides-seismicity correlation in northwestern South America, with a special emphasis in Colombia. For this purpose, we use a dataset of ~167,000 earthquakes, gathered by the Colombian Seismological Network between 1993 and 2017. Most of the events are intermediate-depth earthquakes from the Bucaramanga seismic nest and the Cauca seismic cluster. For this purpose, we implemented a novel approach for the calculation of tidal phases that considers the relative positions of the Earth-Moon-Sun system at the time of the events. After applying the standard Schuster test to the whole dataset and to several earthquake samples (classified by time, location, magnitude and depth), we found strong correlation anomalies with the diurnal and monthly components of the tide (global log(p) values around -7.0 for the diurnal constituent and -12.1 for the monthly constituent), especially for the intermediate depth events. These anomalies suggest that around 16% of the deep earthquakes in Colombia may be triggered by tides, especially when the monthly phase is between 350$^\circ$-10$^\circ$. We attribute our positive results, which favor the tidal-triggering hypothesis, in contrast to previous negative ones to: 1) the size of our dataset, and 2) the method we used to calculate tidal phases. Anyone willing to reproduce our results or to apply our methodology to custom datasets can use the public information system tQuakes that we developed for this work.


title: Multiprecision Arithmetic for Cryptology in C++ - Compile-Time Computations and Beating the Performance of Hand-Optimized Assembly at Run-Time

abstract: We describe a new C++ library for multiprecision arithmetic for numbers in the order of 100--500 bits, i.e., representable with just a few limbs. The library is written in "optimizing-compiler-friendly" C++, with an emphasis on the use of fixed-size arrays and particular function-argument-passing styles (including the avoidance of naked pointers) to allow the limbs to be allocated on the stack or even in registers. Depending on the particular functionality, we get close to, or significantly beat the performance of existing libraries for multiprecision arithmetic that employ hand-optimized assembly code. Most functions in the library are constant-time, which is a necessity for secure implementations of cryptographic protocols. Beyond the favorable runtime performance, our library is, to the best of the author's knowledge, the first library that offers big-integer computations during compile-time. For example, when implementing finite-field arithmetic with a fixed modulus, this feature enables the automatic precomputation (at compile time) of the special modulus-dependent constants required for Barrett and Montgomery reduction. Another application is to parse (at compile-time) a base-10-encoded big-integer literal.


title: Multi-view Hybrid Embedding: A Divide-and-Conquer Approach

abstract: We present a novel cross-view classification algorithm where the gallery and probe data come from different views. A popular approach to tackle this problem is the multi-view subspace learning (MvSL) that aims to learn a latent subspace shared by multi-view data. Despite promising results obtained on some applications, the performance of existing methods deteriorates dramatically when the multi-view data is sampled from nonlinear manifolds or suffers from heavy outliers. To circumvent this drawback, motivated by the Divide-and-Conquer strategy, we propose Multi-view Hybrid Embedding (MvHE), a unique method of dividing the problem of cross-view classification into three subproblems and building one model for each subproblem. Specifically, the first model is designed to remove view discrepancy, whereas the second and third models attempt to discover the intrinsic nonlinear structure and to increase discriminability in intra-view and inter-view samples respectively. The kernel extension is conducted to further boost the representation power of MvHE. Extensive experiments are conducted on four benchmark datasets. Our methods demonstrate overwhelming advantages against the state-of-the-art MvSL based cross-view classification approaches in terms of classification accuracy and robustness.


title: Shortest Dubins Path to a Circle

abstract: The Dubins path problem had enormous applications in path planning for autonomous vehicles. In this paper, we consider a generalization of the Dubins path planning problem, which is to find a shortest Dubins path that starts from a given initial position and heading, and ends on a given target circle with the heading in the tangential direction. This problem has direct applications in Dubins neighborhood traveling salesman problem, obstacle avoidance Dubins path planning problem etc. We characterize the length of the four CSC paths as a function of angular position on the target circle, and derive the conditions which to find the shortest Dubins path to the target circle.


title: Parsimonious Volterra System Identification

abstract: In this short paper, we aim at developing algorithms for sparse Volterra system identification when the system to be identified has infinite impulse response. Assuming that the impulse response is represented as a sum of exponentials and given input-output data, the problem of interest is to find the "simplest" nonlinear Volterra model which is compatible with the a priori information and the collected data. By simplest, we mean the model whose impulse response has the least number of exponentials. The algorithms provided are able to handle both fragmented data and measurement noise. Academic examples at the end of paper show the efficacy of proposed approach.


title: A sequential sampling strategy for extreme event statistics in nonlinear dynamical systems

abstract: We develop a method for the evaluation of extreme event statistics associated with nonlinear dynamical systems, using a small number of samples. From an initial dataset of design points, we formulate a sequential strategy that provides the 'next-best' data point (set of parameters) that when evaluated results in improved estimates of the probability density function (pdf) for a scalar quantity of interest. The approach utilizes Gaussian process regression to perform Bayesian inference on the parameter-to-observation map describing the quantity of interest. We then approximate the desired pdf along with uncertainty bounds utilizing the posterior distribution of the inferred map. The 'next-best' design point is sequentially determined through an optimization procedure that selects the point in parameter space that maximally reduces uncertainty between the estimated bounds of the pdf prediction. Since the optimization process utilizes only information from the inferred map it has minimal computational cost. Moreover, the special form of the metric emphasizes the tails of the pdf. The method is practical for systems where the dimensionality of the parameter space is of moderate size, i.e. order O(10). We apply the method to estimate the extreme event statistics for a very high-dimensional system with millions of degrees of freedom: an offshore platform subjected to three-dimensional irregular waves. It is demonstrated that the developed approach can accurately determine the extreme event statistics using limited number of samples.


title: Unparticle Decay of Neutrinos and its Possible Signatures at a ${\rm Km}^2$ Detector for (3+1) Flavour Framework

abstract: We consider a scenario where ultra high energy neutrinos undergo unparticle decay during its passage from its cosmological source to Earth. The idea of unparticle had been first proposed by Georgi by considering the possible existence of an unknown scale invariant sector at high energies and the unparticles in this sector manifest itself below a dimensional transmutation scale $\Lambda_{\cal U}$. We then explore the possible signature of such decaying neutrinos to unparticles at a square kilometer detector such as IceCube.


title: Dynamical gluon mass at non-zero temperature in instanton vacuum model

abstract: In the framework of the instanton liquid model (ILM), we consider thermal modifications of the gluon properties in different scenarios of temperature $T$ dependence of the average instanton size $\bar{\rho}(T)$ and the instanton density $n(T)$ known from the literature. Due to interactions with instantons, the gluons acquire the dynamical temperature dependent "electric" gluon mass $M_{el}(q,T).$ We found that at small momenta and zero temperature $M_{el}(0,0)\approx362\,{\rm MeV}$ at the phenomenological values of $\bar{\rho}(0)=1/3\,{\rm fm}$ and $n(0)=1\,{\rm fm}^{-4}$, however the $T$-dependence of the mass is very sensitive to the temperature dependence of the instanton vacuum parameters $\bar{\rho}(T),\,n(T)$: it is very mild in case of the lattice-motivated dependence and decreases steeply in the whole range with theoretical parametrization. We see that in region $0<T<T_{c}$ ILM is able to reproduce lattice results for the dynamical "electric" gluon mass.


title: The boundary algebra of a GL$_m$-dimer

abstract: We consider GL$_m$-dimers of triangulations of regular convex $n$-gons, which give rise to a dimer model with boundary $Q$ and a dimer algebra $\Lambda_Q$. Let $e_b$ be the sum of the idempotents of all the boundary vertices, and $\mathcal{B}_Q:= e_b \Lambda_Q e_b$ the associated boundary algebra. In this article we show that given two different triangulations $T_1$ and $T_2$ of the $n$-gon, the boundary algebras are isomorphic, i.e. $e_b \Lambda_{Q_{T_1}} e_b \cong e_b \Lambda_{Q_{T_2}} e_b$.


title: Odd-frequency Pairing in Conventional Josephson Junctions

abstract: Using a simple theoretical model, we demonstrate the emergence of odd-frequency pair amplitudes in conventional Josephson junctions both in the absence of a voltage (DC effect) and in the presence of a finite voltage (AC effect). In both cases, we find that odd-frequency interlead pairing emerges whenever a Josephson current is expected to flow. Additionally, we show that the interlead spin-susceptibility is directly influenced by the presence of the odd-frequency pair amplitudes. Specifically, we find that the spin-susceptibility is suppressed when the odd-frequency component is the largest. By establishing a novel link between the physics of Josephson junctions and odd-frequency pairing, this work demonstrates the importance of odd-frequency pairing for understanding conventional superconducting systems.


title: Meteor showers from active asteroids and dormant comets in near-Earth space: a review

abstract: Small bodies in the solar system are conventionally classified into asteroids and comets. However, it is recently found that a small number of objects can exhibit properties of both asteroids and comets. Some are more consistent with asteroids despite episodic ejections and are labeled as "active asteroids", while some might be aging comets with depleting volatiles. Ejecta produced by active asteroids and/or dormant comets are potentially detectable as meteor showers at the Earth if they are in Earth-crossing orbits, allowing us to retrieve information about the historic activities of these objects. Meteor showers from small bodies with low and/or intermittent activities are usually weak, making shower confirmation and parent association challenging. We show that statistical tests are useful for identifying likely parent-shower pairs. Comprehensive analyses of physical and dynamical properties of meteor showers can lead to deepen understanding on the history of their parents. Meteor outbursts can trace to recent episodic ejections from the parents, and "orphan" showers may point to historic disintegration events. The flourish of NEO and meteor surveys during the past decade has produced a number of high-confidence parent-shower associations, most have not been studied in detail. More work is needed to understand the formation and evolution of these parent-shower pairs.


title: A fourth-order maximum principle preserving operator splitting scheme for three-dimensional fractional Allen-Cahn equations

abstract: In this paper, by using Strang's second-order splitting method, the numerical procedure for the three-dimensional (3D) space fractional Allen-Cahn equation can be divided into three steps. The first and third steps involve an ordinary differential equation, which can be solved analytically. The intermediate step involves a 3D linear fractional diffusion equation, which is solved by the Crank-Nicolson alternating directional implicit (ADI) method. The ADI technique can convert the multidimensional problem into a series of one-dimensional problems, which greatly reduces the computational cost. A fourth-order difference scheme is adopted for discretization of the space fractional derivatives. Finally, Richardson extrapolation is exploited to increase the temporal accuracy. The proposed method is shown to be unconditionally stable by Fourier analysis. Another contribution of this paper is to show that the numerical solutions satisfy the discrete maximum principle under reasonable time step constraint. For fabricated smooth solutions, numerical results show that the proposed method is unconditionally stable and fourth-order accurate in both time and space variables. In addition, the discrete maximum principle is also numerically verified.


title: Semantic Text Analysis for Detection of Compromised Accounts on Social Networks

abstract: Compromised accounts on social networks are regular user accounts that have been taken over by an entity with malicious intent. Since the adversary exploits the already established trust of a compromised account, it is crucial to detect these accounts to limit the damage they can cause. We propose a novel general framework for semantic analysis of text messages coming out from an account to detect compromised accounts. Our framework is built on the observation that normal users will use language that is measurably different from the language that an adversary would use when the account is compromised. We propose to use the difference of language models of users and adversaries to define novel interpretable semantic features for measuring semantic incoherence in a message stream. We study the effectiveness of the proposed semantic features using a Twitter data set. Evaluation results show that the proposed framework is effective for discovering compromised accounts on social networks and a KL-divergence-based language model feature works best.


title: A family of random sup-measures with long-range dependence

abstract: A family of self-similar and translation-invariant random sup-measures with long-range dependence are investigated. They are shown to arise as the limit of the empirical random sup-measure of a stationary heavy-tailed process, inspired by an infinite urn scheme, where same values are repeated at several random locations. The random sup-measure reflects the long-range dependence nature of the original process, and in particular characterizes how locations of extremes appear as long-range clusters represented by random closed sets. A limit theorem for the corresponding point-process convergence is established.


title: Conformal spacelike-timelike correspondence in QCD

abstract: This paper is a study of a spacelike-timelike conformal correspondence in QCD. When the times at vertices are fixed in an A_+ = 0 gauge calculation the distribution of gluons in a highly virtual decay have an exact correspondence with the gluons in the lightcone wavefunction of a high energy dipole with the identification of angles in the timelike case and transverse coordinates in the lightcone wavefunction. Divergences show up when the time integrals are done. A procedure for dropping these divergences, analogous to the Gell-Mann Low procedure in QED, allows one to define a conformal QCD, at least through NLO. Possible uses of such a conformal QCD are discussed.


title: Random Tilings with the GPU

abstract: We present GPU accelerated implementations of Markov chain algorithms to sample random tilings, dimers, and the six-vertex model.


title: The Algorithm of Accumulated Mutual Influence of The Vertices in Semantic Networks

abstract: In this article the algorithm for calculating a mutual influence of the vertices in cognitive maps is introduced. It has been shown, that in the proposed algorithm, there is no problem in comparing with a widely used method - the impulse method, as the proposed algorithm always gives a result regardless of whether impulse process, which corresponds to the weighted directed graph, is a stable or not. Also the result of calculation according to the proposed algorithm does not depend on the initial impulse, and vice versa the initial values of the weights of the vertices influence on the result of calculation. Unlike the impulse method, the proposed algorithm for calculating a mutual influence of the vertices does not violate the scale invariance after increasing the elements of the adjacent matrix, which corresponds to the cognitive map, in the same value. Also in this article the advantages of the proposed algorithm on numerous examples of analysis of cognitive maps are presented.


title: Stability of calcium and magnesium carbonates at lower mantle thermodynamic conditions

abstract: We present a theoretical investigation, based on ab initio calculations and the quasi-harmonic approximation, on the stability properties of magnesium (MgCO$_3$) and calcium (CaCO$_3$) carbonates at high temperatures and pressures. The results indicate that those carbonates should be stable in Earth's lower mantle, instead of dissociating into other minerals, in chemical environments with excesses of SiO$_2$, MgO, or MgSiO$_3$. Therefore, considering the lower mantle chemical composition, filled with major minerals such as MgSiO$_3$ and MgO, calcium and magnesium carbonates are the primary candidates as carbon hosts in that region. For the thermodynamic conditions of the mantle, the results also indicate that carbon should be primarily hosted on MgCO$_3$. Finally, the results indicate that carbon, in the form of free CO$_2$, is unlikely in the lower mantle.


title: Helping or Hurting? Predicting Changes in Users' Risk of Self-Harm Through Online Community Interactions

abstract: In recent years, online communities have formed around suicide and self-harm prevention. While these communities offer support in moment of crisis, they can also normalize harmful behavior, discourage professional treatment, and instigate suicidal ideation. In this work, we focus on how interaction with others in such a community affects the mental state of users who are seeking support. We first build a dataset of conversation threads between users in a distressed state and community members offering support. We then show how to construct a classifier to predict whether distressed users are helped or harmed by the interactions in the thread, and we achieve a macro-F1 score of up to 0.69.


title: Performance Analysis of Massive MIMO Multi-Way Relay Networks with Low-Resolution ADCs

abstract: High power consumption and hardware cost are two barriers for practical massive multiple-input multiple-output (mMIMO) systems. A promising solution is to employ low-resolution analog-to-digital converters (ADCs). In this paper, we consider a general mMIMO multi-way relaying system with a multi-level mixed-ADC architecture, in which each antenna is connected to an ADC pair of an arbitrary resolution. By leveraging on Bussgang's decomposition theorem and Lloyd-Max algorithm for quantization, tight closed-form approximations are derived for the average achievable rates of zero-forcing (ZF) relaying considering both perfect and imperfect channel state information (CSI). To conquer the challenges caused by multi-way relaying, the complicated ZF beam-forming matrix, and the general mixed-ADC structure, we develop a novel method for the achievable rate analysis using the singular-value decomposition (SVD) for Gaussian matrices, distributions of the singular values of Gaussian matrices, and properties of Haar matrices. The results explicitly show the achievable rate behavior in terms of the user and relay transmit powers and the numbers of relay antennas and users. Most importantly, it quantifies the performance degradation caused by low-resolution ADCs and channel estimation error. We demonstrate that the average achievable rate has an almost linear relation with the square of the average of quantization coefficients pertaining to the ADC resolution profile.


title: Scalar Dark Matter Clumps with Angular Momentum

abstract: The behavior of light scalar dark matter has been a subject of much interest recently as it can lead to interesting small scale behavior. In particular, this can lead to the formation of gravitationally bound clumps for these light scalars, including axions. In Ref. [1] we analyzed the behavior of these clumps, assuming spherical symmetry, allowing for both attractive and repulsive self-interactions. There is a maximum allowed mass for the clumps in the case of attractive interactions, and a minimum radius for the clumps in the case of repulsive interactions, which is saturated at large mass. Here we extend this work to include non-spherically symmetric clumps. Since the system tries to re-organize into a BEC ground state, we consider configurations with a conserved non-zero angular momentum, and construct minimum energy configurations at fixed particle number and fixed angular momentum. We find generalizations of the previous spherically symmetric results. In particular, while there is still a maximum mass for the attractive case, its value increases with angular momentum. Also, the minimum radius in the repulsive case is raised to higher radii. We clarify how a recent claim in the literature of an upper bound on angular momentum is due to inaccurate numerics. In a forthcoming paper we shall investigate the possibility of resonance of axion clumps into both visible and hidden sector photons, and analyze how the altered mass and radius from non-zero angular momentum affects the resonance.

categories: hep-ph astro-ph.CO astro-ph.GA gr-qc hep-th


title: Relic Density of Asymmetric Dark Matter in Modified Cosmological Scenarios

abstract: We discuss the relic abundance of asymmetric Dark Matter particles in modified cosmological scenarios where the Hubble rate is changed with respect to the standard cosmological scenario. The modified Hubble rate leaves its imprint on the relic abundance of asymmetric Dark Matter particles if the asymmetric Dark Matter particles freeze--out in this era. For generality we parameterize the modification of the Hubble rate and then calculate the relic abundance of asymmetric Dark Matter particles and anti--particles. We found the abundances for the Dark Matter particles and anti--particles are enhanced in the modified cosmological models. The indirect detection signal is possible for the asymmetric Dark Matter particles due to the increased annihilation rate in the modified cosmological models. Applying Planck data, we find the constraints on the parameters of the modified cosmological models.


title: A Survey on Graph Drawing Beyond Planarity

abstract: Graph Drawing Beyond Planarity is a rapidly growing research area that classifies and studies geometric representations of non-planar graphs in terms of forbidden crossing configurations. Aim of this survey is to describe the main research directions in this area, the most prominent known results, and some of the most challenging open problems.


title: Nonlinear system modeling based on constrained Volterra series estimates

abstract: A simple nonlinear system modeling algorithm designed to work with limited \emph{a priori }knowledge and short data records, is examined. It creates an empirical Volterra series-based model of a system using an $l_{q}$-constrained least squares algorithm with $q\geq 1$. If the system $m\left( \cdot \right) $ is a continuous and bounded map with a finite memory no longer than some known $\tau$, then (for a $D$ parameter model and for a number of measurements $N$) the difference between the resulting model of the system and the best possible theoretical one is guaranteed to be of order $\sqrt{N^{-1}\ln D}$, even for $D\geq N$. The performance of models obtained for $q=1,1.5$ and $2$ is tested on the Wiener-Hammerstein benchmark system. The results suggest that the models obtained for $q>1$ are better suited to characterize the nature of the system, while the sparse solutions obtained for $q=1$ yield smaller error values in terms of input-output behavior.


title: Storing single photons emitted by a quantum memory on a highly excited Rydberg state

abstract: Strong interaction between two single photons is a long standing and important goal in quantum photonics. This would enable a new regime of nonlinear optics and unlock several applications in quantum information science, including photonic quantum gates and deterministic Bell-state measurements. In the context of quantum networks, it would be important to achieve interactions between single photons from independent photon pairs storable in quantum memories. So far, most experiments showing nonlinearities at the single-photon level have used weak classical input light. Here, we demonstrate the storage and retrieval of a paired single photon emitted by an ensemble quantum memory in a strongly nonlinear medium based on highly excited Rydberg atoms. We show that nonclassical correlations between the two photons persist after retrieval from the Rydberg ensemble. Our result is an important step towards deterministic photon-photon interactions, and may enable deterministic Bell-state measurements with multimode quantum memories.


title: Magnetic and pulsational variability of Przybylski's star (HD101065)

abstract: Since its discovery more than half a century ago Przybylski's star (HD101065) continues to excite the astronomical community by the unusual nature of its spectrum, exhibiting exotic element abundances. This star was also the first magnetic chemically peculiar A-type star for which the presence of rapid oscillations was established. Our analysis of newly acquired and historic longitudinal magnetic field measurements indicates that Przybylski's star is also unusual with respect to its extremely slow rotation. Adopting a dipolar structure for the magnetic field and using a sine wave fit to all reported longitudinal magnetic field values over the last 43yr, we find a probable rotation period P_rot ~ 188yr, which however has to be considered tentative as it does not represent a unique solution and has to be verified by future observations. Additionally, based on our own spectropolarimetric material obtained with HARPSpol, we discuss the impact of the anomalous structure of its atmosphere, in particular of the non-uniform horizontal and vertical distributions of chemical elements on the magnetic field measurements and the pulsational variability. Anomalies related to the vertical abundance stratification of Pr and Nd are for the first time used to establish the presence of a radial magnetic field gradient.


title: MontePython 3: boosted MCMC sampler and other features

abstract: MontePython is a parameter inference package for cosmology. We present the latest development of the code over the past couple of years. We explain, in particular, two new ingredients both contributing to improve the performance of Metropolis-Hastings sampling: an adaptation algorithm for the jumping factor, and a calculation of the inverse Fisher matrix, which can be used as a proposal density. We present several examples to show that these features speed up convergence and can save many hundreds of CPU-hours in the case of difficult runs, with a poor prior knowledge of the covariance matrix. We also summarise all the functionalities of MontePython in the current release, including new likelihoods and plotting options.


title: Nonparametric Bayesian label prediction on a large graph using truncated Laplacian regularization

abstract: This article describes an implementation of a nonparametric Bayesian approach to solving binary classification problems on graphs. We consider a hierarchical Bayesian approach with a prior that is constructed by truncating a series expansion of the soft label function using the graph Laplacian eigenfunctions as basis functions. We compare our truncated prior to the untruncated Laplacian based prior in simulated and real data examples to illustrate the improved scalability in terms of size of the underlying graph.