title: Multi-GPU Distributed Parallel Bayesian Differential Topic Modelling

abstract: There is an explosion of data, documents, and other content, and people require tools to analyze and interpret these, tools to turn the content into information and knowledge. Topic modeling have been developed to solve these problems. Topic models such as LDA [Blei et. al. 2003] allow salient patterns in data to be extracted automatically. When analyzing texts, these patterns are called topics. Among numerous extensions of LDA, few of them can reliably analyze multiple groups of documents and extract topic similarities. Recently, the introduction of differential topic modeling (SPDP) [Chen et. al. 2012] performs uniformly better than many topic models in a discriminative setting. There is also a need to improve the sampling speed for topic models. While some effort has been made for distributed algorithms, there is no work currently done using graphical processing units (GPU). Note the GPU framework has already become the most cost-efficient platform for many problems. In this thesis, I propose and implement a scalable multi-GPU distributed parallel framework which approximates SPDP. Through experiments, I have shown my algorithms have a gain in speed of about 50 times while being almost as accurate, with only one single cheap laptop GPU. Furthermore, I have shown the speed improvement is sublinearly scalable when multiple GPUs are used, while fairly maintaining the accuracy. Therefore on a medium-sized GPU cluster, the speed improvement could potentially reach a factor of a thousand. Note SPDP is just a representative of other extensions of LDA. Although my algorithm is implemented to work with SPDP, it is designed to be a general enough to work with other topic models. The speed-up on smaller collections (i.e., 1000s of documents), means that these more complex LDA extensions could now be done in real-time, thus opening up a new way of using these LDA models in industry.


title: $H\to Z\gamma$ in the gauge-Higgs unification

abstract: The decay rate of the Higgs decay $H \to Z \gamma$ is evaluated at the one-loop level in the $SO(5)\times U(1)$ gauge-Higgs unification. Although an infinite number of loops with Kaluza-Klein states contribute to the decay amplitude, there appears the cancellation among the loops, and the decay rate is found to be finite and non-zero. It is found that the decay rate is well approximated by the decay rate in the standard model multiplied by $\cos^2\theta_H$, where $\theta_H$ is the Aharonov-Bohm phase induced by the vacuum expectation value of an extra-dimensional component of the gauge field.


title: Nanowire-based thermoelectric ratchet in the hopping regime

abstract: We study a thermoelectric ratchet consisting of an array of disordered nanowires arranged in parallel on top of an insulating substrate, and contacted asymmetrically to two electrodes. Transport is investigated in the Mott hopping regime, when localized electrons can propagate through the nanowires via thermally assisted hops. When the electronic temperature in the nanowires is different from the phononic one in the substrate, we show that a finite electrical current is generated even in the absence of driving forces between the electrodes. We discuss the device performance both as an energy harvester, when an excess heat from the substrate is converted into useful power, and as a refrigerator, when an external power is supplied to cool down the substrate.


title: Exact observability and controllability for linear neutral type systems

abstract: The problem of exact observability is analyzed for a wide class of neutral type systems by an infinite dimensional approach. The duality with the exact controllabil-ity problem is the main tool. It is based on an explicit expression of a neutral type system which corresponding to the abstract adjoint system. A nontrivial relation is obtained between the initial neutral system and the system obtained via the adjoint abstract state operator. The characterization of the duality between controllability and observability is deduced, and then observability conditions are obtained.


title: Observability analysis and state estimation of lithium-ion batteries in the presence of sensor biases

abstract: This paper investigates the observability of one of the most commonly used equivalent circuit models (ECMs) for lithium-ion batteries and presents a method to estimate the state of charge (SOC) in the presence of sensor biases, highlighting the importance of observability analysis for choosing appropriate state estimation algorithms. Using a differential geometric approach, necessary and sufficient conditions for the nonlinear ECM to be observable are derived and are shown to be different from the conditions for the observability of the linearised model. It is then demonstrated that biases in the measurements, due to sensor ageing or calibration errors, can be estimated by applying a nonlinear Kalman filter to an augmented model where the biases are incorporated into the state vector. Experiments are carried out on a lithium-ion pouch cell and three types of nonlinear filters, the first-order extended Kalman filter (EKF), the second-order EKF and the unscented Kalman filter (UKF) are applied using experimental data. The different performances of the filters are explained from the point of view of observability.


title: Numerical study of the spin-3/2 Ashkin-Teller model

abstract: The study of the Ashkin-Teller model (ATM) of spin-3/2 on a hypercubic lattice is undertaken via Monte Carlo simulation. The phase diagrams are displayed and discussed in the physical parameter space. Rich physical properties are recovered, namely the second order transition and multicritical points. The phase diagrams have been obtained by varying the strength describing the four spin interaction and the single ion potential. This model shows a new high temperature partially ordered phase, called $\langle S\rangle$ and a new Baxter 3/2 ground state which do not exist either in the spin-1/2 ATM or in the spin-1 ATM.


title: On numerical Landau damping for splitting methods applied to the Vlasov-HMF model

abstract: We consider time discretizations of the Vlasov-HMF (Hamiltonian Mean-Field) equation based on splitting methods between the linear and non-linear parts. We consider solutions starting in a small Sobolev neighborhood of a spatially homogeneous state satisfying a linearized stability criterion (Penrose criterion). We prove that the numerical solutions exhibit a scattering behavior to a modified state, which implies a nonlinear Landau damping effect with polynomial rate of damping. Moreover, we prove that the modified state is close to the continuous one and provide error estimates with respect to the time stepsize.


title: Sharp thresholds between finite spread and uniform convergence for a reaction-diffusion equation with oscillating initial data

abstract: We investigate the large-time dynamics of solutions of multi-dimensional reaction-diffusion equations with ignition type nonlinearities. We consider solutions which are in some sense locally persistent at large time and initial data which asymptotically oscillate around the ignition threshold. We show that, as time goes to infinity, any solution either converges uniformly in space to a constant state, or spreads with a finite speed uniformly in all directions. Furthermore, the transition between these two behaviors is sharp with respect to the period vector of the asymptotic profile of the initial data. We also show the convergence to planar fronts when the initial data are asymptotically periodic in one direction.


title: Maxwell's equal area law for Lovelock Thermodynamics

abstract: We present the construction of Maxwell's equal area law for the Guass-Bonnet AdS black holes in $d=5,6$ and third order Lovelock AdS black holes in $d=7,8$. The equal area law can be used to find the number and location of the points of intersection in the plots of Gibbs free energy, so that we can get the thermodynamically preferred solution which corresponds to the first order phase transition. We obtain the radius of the small and large black holes in the phase transition which share the same Gibbs free energy. The case with two critical points is explored in much more details. The latent heat is also studied.


title: Transport through a quantum dot coupled to two Majorana bound states

abstract: We investigate electron transport inside a ring system composed of a quantum dot (QD) coupled to two Majorana bound states confined at the ends of a one-dimensional topological superconductor nanowire. By tuning the magnetic flux threading through the ring, the model system we consider can be switched into states with or without zero-energy modes when the nanowire is in its topological phase. We find that the Fano profile in the conductance spectrum due to the interference between bound and continuum states exhibits markedly different features for these two different situations, which consequently can be used to detect the Majorana zero-energy mode. Most interestingly, as a periodic function of magnetic flux, the conductance shows $2\pi$ periodicity when the two Majorana bound states are nonoverlapping (as in an infinitely long nanowire) but displays $4\pi$ periodicity when the overlapping becomes nonzero (as in a finite length nanowire). We map the model system into a QD--Kitaev ring in the Majorana fermion representation and affirm these different characteristics by checking the energy spectrum.


title: Non-uniqueness results for the anisotropic Calderon problem with data measured on disjoint sets

abstract: In this paper, we give some simple counterexamples to uniqueness for the Calderon problem on Riemannian manifolds with boundary when the Dirichlet and Neumann data are measured on disjoint sets of the boundary. We provide counterexamples in the case of two and three dimensional Riemannian manifolds with boundary having the topology of circular cylinders in dimension two and toric cylinders in dimension three. The construction could be easily extended to higher dimensional Riemannian manifolds.


title: Wave equations on the linear mass Vaidya metric

abstract: We discuss the near singularity region of the linear mass Vaidya metric. In particular we investi- gate the structure in the numerical solutions for the scattering of scalar and electromagnetic metric perturbations from the singularity. In addition to directly integrating the full wave-equation, we use the symmetry of the metric to reduce the problem to that of an ODE. We observe that, around the total evaporation point, quasi-normal like oscillations appear, indicating that this may be an interesting model for the description of the end-point of black hole evaporation.


title: High-order control for symplectic maps

abstract: We revisit the problem of introducing an a priori control for devices that can be modeled via a symplectic map in a neighborhood of an elliptic equilibrium. Using a technique based on Lie transform methods we produce a normal form algorithm that avoids the usual step of interpolating the map with a flow. The formal algorithm is completed with quantitative estimates that bring into evidence the asymptotic character of the normal form transformation. Then we perform an heuristic analysis of the dynamical behavior of the map using the invariant function for the normalized map. Finally, we discuss how control terms of different orders may be introduced so as to increase the size of the stable domain of the map. The numerical examples are worked out on a two dimensional map of H\'enon type.


title: Effects of correlated hopping on electronic ferroelectricity in the extended Falicov-Kimball model in two dimensions

abstract: We use the Hartree-Fock (HF) approximation with the charge-density-wave (CDW) instability to study the effect of correlated hopping on the stability of electronic ferroelectricity in the extended Falicov-Kimball model (FKM) in two dimensions. It is shown that the effects of correlated hopping term are very strong, especially for negative values of the $f$-level energy $E_f$, where they lead to: (i) stabilization of the ferroelectric ground state with a spontaneous hybridization $P_{df}=\langle d^+f \rangle$ for positive values of correlated hopping parameter $t'_d$, (ii) stabilization of the antiferroelectric phase for $t'_d < t'_{d,crit.} < 0$ and (iii) suppression of the ferroelectric ground state for $t'_{d,crit.} \leqslant t'_d < 0$, The effects of correlated hopping on valence transitions are also discussed.


title: Impact of supernova and cosmic-ray driving on the surface brightness of the galactic halo in soft X-rays

abstract: The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.


title: Pseudoscalar decays into lepton pairs from rational approximants

abstract: The pseudoscalar decays into lepton pairs $P\rightarrow\overline{\ell}\ell$ are analyzed with the machinery of Canterbury approximants, an extension of Pad\'e approximants to bivariate functions. This framework provides an ideal model-independent approach to implement all our knowledge of the pseudoscalar transition form factors driving these decays, can be used for data analysis, and allows to include experimental data and theoretical constraints in an easy way, and determine a systematic error. We find that previous theoretical estimates for these branching ratios have underestimated their theoretical uncertainties. From our updated results, the existing experimental discrepancies for $\pi^0\rightarrow e^+e^-$ and $\eta\rightarrow \mu^+\mu^-$ channels cannot be explained unless the doubly-virtual transition form factors behavior -not yet measured- are out of theoretical expectations, which is an interesting result both for anomalous magnetic moments of leptons, and for physics beyond the standard model.


title: Towards a spin-ensemble quantum memory for superconducting qubits

abstract: This article reviews efforts to build a new type of quantum device, which combines an ensemble of electronic spins with long coherence times, and a small-scale superconducting quantum processor. The goal is to store over long times arbitrary qubit states in orthogonal collective modes of the spin-ensemble, and to retrieve them on-demand. We first present the protocol devised for such a multi-mode quantum memory. We then describe a series of experimental results using NV center spins in diamond, which demonstrate its main building blocks: the transfer of arbitrary quantum states from a qubit into the spin ensemble, and the multi-mode retrieval of classical microwave pulses down to the single-photon level with a Hahn-echo like sequence. A reset of the spin memory is implemented in-between two successive sequences using optical repumping of the spins.


title: The harmonic transvector algebra in two vector variables

abstract: The decomposition of polynomials of one vector variable into irreducible modules for the orthogonal group is a crucial result in harmonic analysis which makes use of the Howe duality theorem and leads to the study of spherical harmonics. The aim of the present paper is to describe a decomposition of polynomials in two vector variables and to obtain projection operators on each of the irreducible components. To do so, a particular transvector algebra will be used as a new dual partner for the orthogonal group leading to a generalisation of the classical Howe duality. The results are subsequently used to obtain explicit projection operators and formulas for integration of polynomials over the associated Stiefel manifold.


title: Generalized conditional gradient: analysis of convergence and applications

abstract: The objectives of this technical report is to provide additional results on the generalized conditional gradient methods introduced by Bredies et al. [BLM05]. Indeed , when the objective function is smooth, we provide a novel certificate of optimality and we show that the algorithm has a linear convergence rate. Applications of this algorithm are also discussed.


title: Feasibility of Diagrammatic Monte-Carlo based on weak-coupling expansion in asymptotically free theories: case study of $O(N)$ sigma-model in the large-$N$ limit

abstract: We discuss the feasibility of applying Diagrammatic Monte-Carlo algorithms to the weak-coupling expansions of asymptotically free quantum field theories, taking the large-$N$ limit of the $O(N)$ sigma-model as the simplest example where exact results are available. We use stereographic mapping from the sphere to the real plane to set up the perturbation theory, which results in a small bare mass term proportional to the coupling $\lambda$. Counting the powers of coupling associated with higher-order interaction vertices, we arrive at the double-series representation for the dynamically generated mass gap in powers of both $\lambda$ and $\log(\lambda)$, which converges quite quickly to the exact non-perturbative answer. We also demonstrate that it is feasible to obtain the coefficients of these double series by a Monte-Carlo sampling in the space of Feynman diagrams. In particular, the sign problem of such sampling becomes milder at small $\lambda$, that is, close to the continuum limit.


title: BRST Analysis of the Supersymmetric Higher Spin Field Models

abstract: We develop the BRST approach for all massless integer and half-integer higher spins in 4D Minkowski space, using the two component spinor nota- tion and develop the Lagrangian formulation for supersymmetric higher spin models. It is shown that the problem of second class constraints disappears and the BRST procedure becomes much more simple than in tensorial nota- tion. Furthermore, we demonstrate that the BRST procedure automatically provides extra auxiliary components that belong in the set of supersymmetry auxiliary components. Finally, we demonstrate how supersymmetry transfo- rmations are realized in such an approach. As a result, we conclude that the BRST approach to higher spin supersymmetric theories allows to derive both the Lagrangian and the supersymmetry transformations. Although most part of the work is devoted to massless component supersymmetric models, we also discuss generalization for massive component supersymmetric models and for superfield models.


title: States of $\rho B^* \bar{B}^*$ with $J=3$ within the Fixed Center Approximation to Faddeev equations

abstract: In this work we study the $\rho B^*\bar{B}^*$ three-body system solving the Faddeev equations in the fixed center approximation. We assume the $B^*\bar{B}^*$ system forming a cluster, and in terms of the two-body $\rho B^*$ unitarized scattering amplitudes in the local Hidden Gauge approach we find a new $I(J^{PC})=1(3^{--})$ state. The mass of the new state corresponds to a two particle invariant mass of the $\rho B^\ast$ system close to the resonant energy of the $B^\ast_2(5747)$, indicating that the role of this $J=2$ resonance is important in the dynamical generation of the new state.


title: Monolayer charge-neutral graphene on platinum with extremely weak electron-phonon coupling

abstract: Epitaxial growth of graphene on transition metal substrates is an important route for obtaining large scale graphene. However, the interaction between graphene and the substrate often leads to multiple orientations, distorted graphene band structure, large doping and strong electron-phonon coupling. Here we report the growth of monolayer graphene with high crystalline quality on Pt(111) substrate by using a very low concentration of an internal carbon source with high annealing temperature. The controlled growth leads to electronically decoupled graphene: it is nearly charge neutral and has extremely weak electron-phonon coupling (coupling strength $\lambda$ $\approx$ 0.056) as revealed by angle-resolved photoemission spectroscopic measurements. The thermodynamics and kinetics of the carbon diffusion process is investigated by DFT calculation. Such graphene with negligible graphene-substrate interaction provides an important platform for fundamental research as well as device applications when combined with a nondestructive sample transfer technique.


title: Radio Resource Allocation in LTE-Advanced Cellular Networks with M2M Communications

abstract: Machine-to-machine (M2M) communications are expected to provide ubiquitous connectivity between machines without the need of human intervention. To support such a large number of autonomous devices, the M2M system architecture needs to be extremely power and spectrally efficient. This article thus briefly reviews the features of M2M services in the third generation (3G) long-term evolution and its advancement (LTE-Advanced) networks. Architectural enhancements are then presented for supporting M2M services in LTE-Advanced cellular networks. To increase spectral efficiency, the same spectrum is expected to be utilized for human-to-human (H2H) communications as well as M2M communications. We therefore present various radio resource allocation schemes and quantify their utility in LTE-Advanced cellular networks. System-level simulation results are provided to validate the performance effectiveness of M2M communications in LTE-Advanced cellular networks.


title: Generalized Alexander Polynomial Invariants

abstract: We propose an algorithm which allows to derive the generalized Alexander polynomial invariants of knots and links with the help of the q,p-numbers, appearing in bosonic two-parameter quantum algebra. These polynomials turn into HOMFLY ones by applying special parametrization. The Jones polynomials can be also obtained by using this algorithm.


title: Random graphs and Lindstrom quantifiers for natural graph properties

abstract: We study zero-one laws for random graphs. We focus on the following question that was asked by many: Given a graph property P, is there a language of graphs able to express P while obeying the zero-one law? Our results show that on the one hand there is a (regular) language able to express connectivity and k-colorability for any constant k and still obey the zero-one law. On the other hand we show that in any (semiregular) language strong enough to express Hamiltonicity one can interpret arithmetic and thus the zero-one law fails miserably. This answers a question of Blass and Harary.


title: JWST observations of stellar occultations by solar system bodies and rings

abstract: In this paper we investigate the opportunities provided by the James Webb Space Telescope (JWST) for significant scientific advances in the study of solar system bodies and rings using stellar occultations. The strengths and weaknesses of the stellar occultation technique are evaluated in light of JWST's unique capabilities. We identify several possible JWST occultation events by minor bodies and rings, and evaluate their potential scientific value. These predictions depend critically on accurate a priori knowledge of the orbit of JWST near the Sun-Earth Lagrange-point 2 (L2). We also explore the possibility of serendipitous stellar occultations by very small minor bodies as a by-product of other JWST observing programs. Finally, to optimize the potential scientific return of stellar occultation observations, we identify several characteristics of JWST's orbit and instrumentation that should be taken into account during JWST's development.


title: Chemical Tracers of Pre-Brown Dwarf Cores Formed Through Turbulent Fragmentation

abstract: A gas-grain time dependent chemical code, UCL\_CHEM, has been used to investigate the possibility of using chemical tracers to differentiate between the possible formation mechanisms of brown dwarfs. In this work, we model the formation of a pre-brown dwarf core through turbulent fragmentation by following the depth-dependent chemistry in a molecular cloud through the step change in density associated with an isothermal shock and the subsequent freefall collapse once a bound core is produced. Trends in the fractional abundance of molecules commonly observed in star forming cores are then explored to find a diagnostic for identifying brown dwarf mass cores formed through turbulence. We find that the cores produced by our models would be bright in CO and NH$_3$ but not in HCO$^+$. This differentiates them from models using purely freefall collapse as such models produce cores that would have detectable transitions from all three molecules.


title: A strong Dixmier-Moeglin equivalence for quantum Schubert cells

abstract: Dixmier and Moeglin gave an algebraic condition and a topological condition for recognising the primitive ideals among the prime ideals of the universal enveloping algebra of a finite-dimensional complex Lie algebra; they showed that the primitive, rational, and locally closed ideals coincide. In modern terminology, they showed that the universal enveloping algebra of a finite-dimensional complex Lie algebra satisfies the Dixmier-Moeglin equivalence. We define quantities which measure how "close" an arbitrary prime ideal of a noetherian algebra is to being primitive, rational, and locally closed; if every prime ideal is equally "close" to each of these three properties, then we say that the algebra satisfies the strong Dixmier-Moeglin equivalence. Using the example of the universal enveloping algebra of sl_2(C), we show that the strong Dixmier-Moeglin equivalence is stronger than the Dixmier-Moeglin equivalence. For a simple complex Lie algebra g, a non root of unity q\neq 0 in an infinite field K, and an element w of the Weyl group of g, De Concini, Kac, and Procesi have constructed a subalgebra U_q[w] of the quantised enveloping K-algebra U_q(g). These quantum Schubert cells U_q[w] are known to satisfy the Dixmier-Moeglin equivalence and we show that they in fact satisfy the strong Dixmier-Moeglin equivalence. Along the way, we show that commutative affine domains, uniparameter quantum tori, and uniparameter quantum affine spaces satisfy the strong Dixmier-Moeglin equivalence.


title: Securing the MIMO Wiretap Channel with Polar Codes and Encryption

abstract: Polar codes have been proven to be capacity achieving for any binary-input discrete memoryless channel, while at the same time they can reassure secure and reliable transmission over the single-input single-output wireless channel. However, the use of polar codes to secure multiple-antenna transmission and reception has not yet been reported in the open literature. In this paper, we assume a multiple-input multiple-output wiretap channel, where the legitimate receiver and the eavesdropper are equipped with the same number of antennas. We introduce a protocol that exploits the properties of both physical and media access control layer security by employing polar coding and encryption techniques in a hybrid manner in order to guarantee secure transmission. A novel security technique is also proposed, where a cryptographic key is generated based on the information transmitted and renewed every transmission block without the need for a separate key exchange method. Finally, to illustrate the effectiveness of the proposed protocol, we prove the weak and strong security conditions, and we provide a practical method to achieve computational security for the cases where these conditions cannot be established.


title: Soft-Defined Heterogeneous Vehicular Network: Architecture and Challenges

abstract: Heterogeneous Vehicular NETworks (HetVNETs) can meet various quality-of-service (QoS) requirements for intelligent transport system (ITS) services by integrating different access networks coherently. However, the current network architecture for HetVNET cannot efficiently deal with the increasing demands of rapidly changing network landscape. Thanks to the centralization and flexibility of the cloud radio access network (Cloud-RAN), soft-defined networking (SDN) can conveniently be applied to support the dynamic nature of future HetVNET functions and various applications while reducing the operating costs. In this paper, we first propose the multi-layer Cloud RAN architecture for implementing the new network, where the multi-domain resources can be exploited as needed for vehicle users. Then, the high-level design of soft-defined HetVNET is presented in detail. Finally, we briefly discuss key challenges and solutions for this new network, corroborating its feasibility in the emerging fifth-generation (5G) era.


title: A twisted bicanonical system with base points

abstract: By a theorem of Reider, a twisted bicanonical system, that means a linear system of divisors numerically equivalent to a bicanonical divisor, on a minimal surface of general type, is base point free if $K^2_S \geq 5$. Twisted bicanonical systems with base points are known in literature only for $K^2=1,2$. We prove in this paper that all surfaces in a family of surfaces with $K^2=3$ constructed in a previous paper with G. Bini and J. Neves have a twisted bicanonical system (different from the bicanonical system) with two base points. We show that the map induced by the above twisted bicanonical system is birational, and describe in detail the closure of its image and its singular locus. Inspired by this description, we reduced the problem of constructing a minimal surface of general type with $K^2=3$ whose bicanonical system has base points, under some reasonable assumptions, to the problem of constructing a curve in $\mathbb{P}^3$ with certain properties.


title: An extension of the Ehrenfeucht-Fraisse game for first order logics augmented with Lindstrom quantifiers

abstract: We propose an extension of the Ehrenfeucht-Fraisse game able to deal with logics augmented with Lindstrom quantifiers. We describe three different games with varying balance between simplicity and ease of use.


title: Collective Prediction of Individual Mobility Traces with Exponential Weights

abstract: We present and test a sequential learning algorithm for the short-term prediction of human mobility. This novel approach pairs the Exponential Weights forecaster with a very large ensemble of experts. The experts are individual sequence prediction algorithms constructed from the mobility traces of 10 million roaming mobile phone users in a European country. Average prediction accuracy is significantly higher than that of individual sequence prediction algorithms, namely constant order Markov models derived from the user's own data, that have been shown to achieve high accuracy in previous studies of human mobility prediction. The algorithm uses only time stamped location data, and accuracy depends on the completeness of the expert ensemble, which should contain redundant records of typical mobility patterns. The proposed algorithm is applicable to the prediction of any sufficiently large dataset of sequences.


title: Characters of real special 2-groups

abstract: It is well-known that special 2-groups can be described in terms of quadratic maps over fields of characteristic 2. In this article we develop methods to compute conjugacy classes, complex representations and characters of a real special $2$-group using quadratic maps alone.


title: The eccentric short-period orbit of the supergiant fast X-ray transient HD 74194 (=LM Vel)

abstract: Aims. We present the first orbital solution for the O-type supergiant star HD 74194, which is the optical counterpart of the supergiant fast X-ray transient IGR J08408-4503. Methods. We measured the radial velocities in the optical spectrum of HD 74194, and we determined the orbital solution for the first time. We also analysed the complex H{\alpha} profile. Results. HD 74194 is a binary system composed of an O-type supergiant and a compact object in a short-period ($P=9.5436\pm0.0002$ d) and high-eccentricity ($e=0.63\pm0.03$) orbit. The equivalent width of the H{\alpha} line is not modulated entirely with the orbital period, but seems to vary in a superorbital period ($P=285\pm10$ d) nearly 30 times longer than the orbital one.


title: Execution of Compound Multi-Kernel OpenCL Computations in Multi-CPU/Multi-GPU Environments

abstract: Current computational systems are heterogeneous by nature, featuring a combination of CPUs and GPUs. As the latter are becoming an established platform for high-performance computing, the focus is shifting towards the seamless programming of these hybrid systems as a whole. The distinct nature of the architectural and execution models in place raises several challenges, as the best hardware configuration is behaviour and workload dependent. In this paper, we address the execution of compound, multi-kernel, OpenCL computations in multi-CPU/multi-GPU environments. We address how these computations may be efficiently scheduled onto the target hardware, and how the system may adapt itself to changes in the workload to process and to fluctuations in the CPU's load. An experimental evaluation attests the performance gains obtained by the conjoined use of the CPU and GPU devices, when compared to GPU-only executions, and also by the use of data-locality optimizations in CPU environments.


title: On location of zeros of the first derivative

abstract: Let p(z) be a complex polynomial of degree n. Let C be a circle containing its n-1 zeros, having its center in the centroid of these zeros. We prove that C must contain at least int((n-1):2) zeros of its derivative.


title: Approximating Strategic Abilities under Imperfect Information: a Naive Approach

abstract: Alternating-time temporal logic (ATL) allows to specify requirements on abilities that different agents should (or should not) possess in a multi-agent system. However, model checking ATL specifications in realistic systems is computationally hard. In particular, if the agents have imperfect information about the global state of the system, the complexity ranges from Delta2P to undecidable, depending on the syntactic and semantic details. The problem is also hard in practice, as evidenced by several recent attempts to tackle it. On the other hand, model checking of alternating epistemic mu-calculus can have a distinctly lower computational complexity. In this work, we look at the idea of approximating the former problem by the verification of its "naive" translations to the latter. In other words, we look at what happens when one uses the (incorrect) fixpoint algorithm to verify formulae of ATL with imperfect information.


title: Recoller pour s\'eparer

abstract: We introduce the notion of a separator for a morphism of schemes f:T\to S; in particular, it is universal among morphisms from T to separated S-schemes. A separator is a local isomorphism; this property conveys the intuition of gluing some affine covering more, in order to make the scheme separated. When f is quasi-separated, its separator exists if and only if the schematic closure of the diagonal projects on both factors by flat morphisms of finite type. In particular, f admits a separator if T is Noetherian Dedekind and S=Spec(Z), or if f is \'etale of finite presentation and S is normal. Any normal scheme of finite type over a Noetherian ring admits an open subset containing all the points of codimension 1, which has a separator. A contrario, we give several examples of morphisms f that do not admit a separator. As an application, we attach to every smooth scheme T over a normal base S a morphism to a separated \'etale S-scheme of finite presentation, which is universal (a kind of separated alternative for "scheme of connected components of the fibres"). This simultaneously generalizes the classical case where the base is a field, and the case of a smooth and proper morphism (Stein factorisation).


title: Anomalous longitudinal relaxation of nuclear spins in CaF$_2$

abstract: We consider the effect of non-secular resonances for interacting nuclear spins in solids which were predicted theoretically to exist in the presence of strong static and strong radio-frequency magnetic fields. These resonances imply corrections to the standard secular approximation for the nuclear spin-spin interaction in solids, which, in turn, should lead to an anomalous longitudinal relaxation in nuclear magnetic resonance experiments. In this article, we investigate the feasibility of the experimental observation of this anomalous longitudinal relaxation in calcium fluoride (CaF$_2$) and conclude that such an observation is realistic.


title: Work fluctuations in a nonlinear micromechanical oscillator driven far from thermal equilibrium

abstract: We explore fluctuation relations in a periodically driven micromechanical torsional oscillator. In the linear regime where the modulation is weak, we verify that the ratio of the work variance to the mean work is constant, consistent with conventional fluctuation theorems. We then increase the amplitude of the periodic drive so that the response becomes nonlinear and two nonequilibrium oscillation states coexist. Due to interstate transitions, the work variance exhibits a peak at the driving frequency at which the occupation of the two states is equal. Moreover, the work fluctuations depend exponentially on the inverse noise intensity. Our data are consistent with recent theories on systems driven into bistability that predict generic behaviors different from conventional fluctuation theorems.


title: Airglow and magnetic field disturbances over Brazilian region during Chile tsunami (2015)

abstract: In this work, we present first report on disturbances over Brazilian atmosphere on 16--17 September, 2015 following the Chile tsunami occurrence. Using all-sky imager and magnetometer located at 2330 km away from the epicenter, the presence of disturbances is noted 1--3 hours after the tsunami beginning time and during time which seismic tremor was also felt in the region. We argue that their presence towards continent at 2000-3000 km away from the epicenter offers another example of similar atmospheric response as those observed during Tohoku-Oki tsunami, 2011. This similarity and their appearance during seismic tremor over the region classify them to be of tsunamigenic and/or seismogenic nature.


title: Exact two-dimensionalization of rapidly rotating large-Reynolds-number flows

abstract: We consider the flow of a Newtonian fluid in a three-dimensional domain, rotating about a vertical axis and driven by a vertically invariant horizontal body-force. This system admits vertically invariant solutions that satisfy the 2D Navier-Stokes equation. At high Reynolds number and without global rotation, such solutions are usually unstable to three-dimensional perturbations. By contrast, for strong enough global rotation, we prove rigorously that the 2D (and possibly turbulent) solutions are stable to vertically dependent perturbations: the flow becomes 2D in the long-time limit. These results shed some light on several fundamental questions of rotating turbulence: for arbitrary Reynolds number and small enough Rossby number, the system is attracted towards purely 2D flow solutions, which display no energy dissipation anomaly and no cyclone-anticyclone asymmetry. Finally, these results challenge the applicability of wave turbulence theory to describe stationary rotating turbulence in bounded domains.


title: Shaft Inflation and the Planck satellite observations

abstract: A new family of inflation models is introduced and studied. The models are characterised by a scalar potential which, far from the origin, approaches an inflationary plateau in a power-law manner, while near the origin becomes monomial, as in chaotic inflation. The models are obtained in the context of global supersymmetry starting with a superpotential, which interpolates from a generalised monomial to an O'Raifearteagh form for small to large values of the inflaton field respectively. It is demonstrated that the observables obtained, such as the scalar spectral index, its running and the tensor to scalar ratio, are in excellent agreement with the latest observations, without any fine-tuning. Moreover, by widening mildly the shaft in field space, it is shown that sizable tensors can be generated, which may well be observable in the near future.


title: A comprehensive study of noble gases and nitrogen in Hypatia, a diamond-rich pebble from SW Egypt

abstract: This is a follow-up study of a work by Kramers et al. (2013) on an unusual diamond-rich rock found in the SW side of the Libyan Desert Glass strewn field. This pebble, called Hypatia, is composed of almost pure carbon. Transmission Electron Microscopy and X-ray diffraction results reveal that Hypatia is made of defect-rich diamond containing lonsdaleite and deformation bands. These characteristics are compatible with an impact origin on Earth and/or in space. We analyzed concentrations and isotopic compositions of all five noble gases and nitrogen in several mg sized Hypatia samples. These data confirm that Hypatia is extra-terrestrial. The sample is rich in trapped noble gases with an isotopic composition close to the meteoritic Q component. 40Ar/36Ar ratios in individual steps are as low as 0.4. Concentrations of cosmic-ray produced 21Ne correspond to a nominal cosmic-ray exposure age of ca. 0.1 Myr if produced in a typical m-sized meteoroid. Such an atypically low nominal exposure age suggests high shielding in a considerably larger body. In addition to the Xe-Q composition, an excess of radiogenic 129Xe (from the decay of extinct 129I) is observed (129Xe/132Xe = 1.18 +/- 0.03). Two N components are present, an isotopically heavy component ({\delta}15N = +20 permil) released at low temp. and a major light component ({\delta}15N = -110 permil) at higher temp. This disequilibrium in N suggests that the diamonds in Hypatia were formed in space. Our data are broadly consistent with concentrations and isotopic compositions of noble gases in at least three different types of carbon-rich meteoritic materials. However, Hypatia does not seem to be related to any of these materials, but may have sampled a similar cosmochemical reservoir. Our study does not confirm the presence of exotic noble gases that led Kramers et al. to propose that Hypatia is a remnant of a comet that impacted the Earth.


title: Efficient Unsupervised Temporal Segmentation of Motion Data

abstract: We introduce a method for automated temporal segmentation of human motion data into distinct actions and compositing motion primitives based on self-similar structures in the motion sequence. We use neighbourhood graphs for the partitioning and the similarity information in the graph is further exploited to cluster the motion primitives into larger entities of semantic significance. The method requires no assumptions about the motion sequences at hand and no user interaction is required for the segmentation or clustering. In addition, we introduce a feature bundling preprocessing technique to make the segmentation more robust to noise, as well as a notion of motion symmetry for more refined primitive detection. We test our method on several sensor modalities, including markered and markerless motion capture as well as on electromyograph and accelerometer recordings. The results highlight our system's capabilities for both segmentation and for analysis of the finer structures of motion data, all in a completely unsupervised manner.


title: SN 2014J at M82: I. A middle-class type Ia supernova by all spectroscopic metrics

abstract: We present the intensive spectroscopic follow up of the type Ia supernova (SN Ia) 2014J in the starburst galaxy M82. Twenty-seven optical spectra have been acquired from January 22nd to September 1st 2014 with the Isaac Newton (INT) and William Herschel (WHT) Telescopes. After correcting the observations for the recession velocity of M82 and for Milky Way and host galaxy extinction, we measured expansion velocities from spectral line blueshifts and pseudo-equivalent width of the strongest features in the spectra, which gives an idea on how elements are distributed within the ejecta. We position SN 2014J in the Benetti (2005), Branch et al. (2006) and Wang et al. (2009) diagrams. These diagrams are based on properties of the Si II features and provide dynamical and chemical information about the SN ejecta. The nearby SN 2011fe, which showed little evidence for reddening in its host galaxy, is shown as a reference for comparisons. SN 2014J is a border-line object between the Core-normal (CN) and Broad-line (BL) groups, which corresponds to an intermediate position between Low Velocity Gradient (LVG) and High Velocity Gradient (HVG) objects. SN 2014J follows the R(Si II)-\Delta m15 correlation, which confirms its classification as a relatively normal SN Ia. Our description of the SN Ia in terms of the evolution of the pseudo-equivalent width of various ions as well as the position in the various diagrams put this specific SN Ia into the overall sample of SN Ia.


title: Spacings - An Example for Universality in Random Matrix Theory

abstract: Universality of local eigenvalue statistics is one of the most striking phenomena of Random Matrix Theory, that also accounts for a lot of the attention that the field has attracted over the past 15 years. In this paper we focus on the empirical spacing distribution and its Kolmogorov distance from the universal limit. We describe new results, some analytical, some numerical, that are contained in [27]. A large part of the paper is devoted to explain basic definitions and facts of Random Matrix Theory, culminating in a sketch of the proof of a weak version of convergence for the empirical spacing distribution $\sigma_N$.


title: Rectangular Young tableaux and the Jacobi ensemble

abstract: It has been shown by Pittel and Romik that the random surface associated with a large rectangular Young tableau converges to a deterministic limit. We study the fluctuations from this limit along the edges of the rectangle. We show that in the corner, these fluctuations are gaussian wheras, away from the corner and when the rectangle is a square, the fluctuations are given by the Tracy-Widom distribution. Our method is based on a connection with the Jacobi ensemble.


title: Mean-field predictions for a dipolar Bose-Einstein condensate with $^{164}$Dy

abstract: We study dipolar Bose-Einstein condensates for a realistic set of parameters close to actual experimental setups with dysprosium. Our analysis is based on the extended Gross-Pitaevskii equation, which we solve numerically exact on a three-dimensional grid. We present stability and phase diagrams, and study the expansion dynamics of dipolar condensates. Our calculations show signatures of the dipole-dipole interaction in terms of structured states and a deviation of the well-known inversion of the aspect ratio of the cloud during a time of flight.


title: Depth-Resolved Composition and Electronic Structure of Buried Layers and Interfaces in a LaNiO$_3$/SrTiO$_3$ Superlattice from Soft- and Hard- X-ray Standing-Wave Angle-Resolved Photoemission

abstract: LaNiO$_3$ (LNO) is an intriguing member of the rare-earth nickelates in exhibiting a metal-insulator transition for a critical film thickness of about 4 unit cells [Son et al., Appl. Phys. Lett. 96, 062114 (2010)]; however, such thin films also show a transition to a metallic state in superlattices with SrTiO$_3$ (STO) [Son et al., Appl. Phys. Lett. 97, 202109 (2010)]. In order to better understand this transition, we have studied a strained LNO/STO superlattice with 10 repeats of [4 unit-cell LNO/3 unit-cell STO] grown on an (LaAlO$_3$)$_{0.3}$(Sr$_2$AlTaO$_6$)$_{0.7}$ substrate using soft x-ray standing-wave-excited angle-resolved photoemission (SWARPES), together with soft- and hard- x-ray photoemission measurements of core levels and densities-of-states valence spectra. The experimental results are compared with state-of-the-art density functional theory (DFT) calculations of band structures and densities of states. Using core-level rocking curves and x-ray optical modeling to assess the position of the standing wave, SWARPES measurements are carried out for various incidence angles and used to determine interface-specific changes in momentum-resolved electronic structure. We further show that the momentum-resolved behavior of the Ni 3d eg and t2g states near the Fermi level, as well as those at the bottom of the valence bands, is very similar to recently published SWARPES results for a related La$_{0.7}$Sr$_{0.3}$MnO$_3$/SrTiO$_3$ superlattice that was studied using the same technique (Gray et al., Europhysics Letters 104, 17004 (2013)), which further validates this experimental approach and our conclusions. Our conclusions are also supported in several ways by comparison to DFT calculations for the parent materials and the superlattice, including layer-resolved density-of-states results.


title: Shielding $^2\Sigma$ ultracold dipolar molecular collisions with electric fields

abstract: The prospects for shielding ultracold, paramagnetic, dipolar molecules from inelastic and chemical collisions are investigated. Molecules placed in their first rotationally excited states are found to exhibit effective long-range repulsion for applied electric fields above a certain critical value, as previously shown for non-paramagnetic molecules. This repulsion can safely allow the molecules to scatter while reducing the risk of inelastic or chemically reactive collisions. Several molecular species of $^2\Sigma$ molecules of experimental interest -- RbSr, SrF, BaF, and YO -- are considered, and all are shown to exhibit orders of magnitude suppression in quenching rates in a sufficiently strong laboratory electric field. It is further shown that, for these molecules described by Hund's coupling case b, electronic and nuclear spins play the role of spectator with respect to the shielding.


title: Uniform asymptotic behaviour of Jacobi-$\sn$ near a singular point. The Lost formula from handbooks for elliptic functions

abstract: In this work we construct uniform asymptotic expansion of $\sn(t|m)$ - Jacobi when $m\to1-0$. The constructed expansion is valid over more than a half of period. The turning point is included into the interval of validity for the approximation. In addition we obtain the asymptotic formula for elliptic integral of the first kind and discuss the differences with the same formula from a handbook.


title: Fuel-Optimal Centralized Coordination of Truck Platooning Based on Shortest Paths

abstract: Platooning is a way to significantly reduce fuel consumption of trucks. Vehicles that drive at close inter-vehicle distance assisted by automatic controllers experience substantially lower air-drag. In this paper, we deal with the problem of coordinating the formation and the breakup of platoons in a fuel-optimal way. We formulate an optimization problem which accounts for routing, speed-dependent fuel consumption, and platooning decisions. An algorithm to obtain an approximate solution to the problem is presented. It first determines the shortest path for each truck. Then, possible platoon configurations are identified. For a certain platoon configuration the optimal speed profile is the solution to a convex program. The algorithm is illustrated by a realistic example.


title: Internal structure of Pluto and Charon with an iron core

abstract: Pluto has been observed by the New Horizons space probe to have some relatively fresh ice on the old ices covering most of the surface. Pluto was thought to consist of only a rocky core below the ice. Here I show that Pluto can have an iron core, as can also its companion Charon, which has recently been modelled to have one. The presence of an iron core means the giant impact origin calculations should be redone to include iron and thus higher temperatures. An iron core leads to the possibility of a different geology. An originally molten core becomes solid later, with contraction and a release of latent heat. The space vacated allows the upper rock layers to flow downwards at some locations at the surface of the core, and some of the ice above the rock to descend, filling the spaces left by the rock motion downwards. These phenomena can lead to the forces recently deforming the icy surface of Pluto, and in a lesser way, of Charon.


title: Efficient computation of capillary-gravity generalized solitary waves

abstract: This paper is devoted to the computation of capillary-gravity solitary waves of the irrotational incompressible Euler equations with free surface. The numerical study is a continuation of a previous work in several points: an alternative formulation of the Babenko-type equation for the wave profiles, a detailed description of both the numerical resolution and the analysis of the internal flow structure under a solitary wave. The numerical code used in this study is provided in open source for those interested readers.

categories: physics.flu-dyn math.NA nlin.PS physics.class-ph physics.comp-ph


title: Hecke algebras for $\mathrm{GL}_n$ over local fields

abstract: We study the local Hecke algebra $\mathcal{H}_{G}(K)$ for $G = \mathrm{GL}_n$ and $K$ a non-archimedean local field of characteristic zero. We show that for $G = \mathrm{GL}_2$ and any two such fields $K$ and $L$, there is a Morita equivalence $\mathcal{H}_{G}(K) \sim_M \mathcal{H}_{G}(L)$, by using the Bernstein decomposition of the Hecke algebra and determining the intertwining algebras that yield the Bernstein blocks up to Morita equivalence. By contrast, we prove that for $G = \mathrm{GL}_n$, there is an algebra isomorphism $\mathcal{H}_{G}(K) \cong \mathcal{H}_{G}(L)$ which is an isometry for the induced $L^1$-norm if and only if there is a field isomorphism $K \cong L$.


title: Reliable and Efficient Autonomous Driving: the Need for Heterogeneous Vehicular Networks

abstract: Autonomous driving technology has been regarded as a promising solution to reduce road accidents and traffic congestion, as well as to optimize the usage of fuel and lane. Reliable and high efficient Vehicle-to-Vehicle (V2V) and Vehicle-to-Infrastructure (V2I) communications are essential to let commercial autonomous driving vehicles be on the road before 2020. The current paper firstly presents the concept of Heterogeneous Vehicular NETworks (HetVNETs) for autonomous driving, in which an improved protocol stack is proposed to satisfy the communication requirements of not only safety but also non-safety services. We then consider and study in detail several typical scenarios for autonomous driving. In order to tackle the potential challenges raised by the autonomous driving vehicles in HetVNETs, new techniques from transmission to networking are proposed as potential solutions.


title: A Bramble like Witness for Large Branch-Width

abstract: Using the new framework of graph separations developed recently by Diestel and Oum, we show that there is a natural witness to large branch-width which can be defined in a similar way to a bramble. We also apply our methods to a more recent object of study in the area, a common generalization of tangles and blocks, called profiles.


title: Photon noise from chaotic and coherent millimeter-wave sources measured with horn-coupled, aluminum lumped-element kinetic inductance detectors

abstract: We report photon-noise limited performance of horn-coupled, aluminum lumped-element kinetic inductance detectors at millimeter wavelengths. The detectors are illuminated by a millimeter-wave source that uses an active multiplier chain to produce radiation between 140 and 160 GHz. We feed the multiplier with either amplified broadband noise or a continuous-wave tone from a microwave signal generator. We demonstrate that the detector response over a 40 dB range of source power is well-described by a simple model that considers the number of quasiparticles. The detector noise-equivalent power (NEP) is dominated by photon noise when the absorbed power is greater than approximately 1 pW, which corresponds to $\mathrm{NEP} \approx 2 \times 10^{-17} \, \mathrm{W} \, \mathrm{Hz}^{-1/2}$, referenced to absorbed power. At higher source power levels we observe the relationships between noise and power expected from the photon statistics of the source signal: $\mathrm{NEP} \propto P$ for broadband (chaotic) illumination and $\mathrm{NEP} \propto P^{1/2}$ for continuous-wave (coherent) illumination.


title: Unveiling the optical properties of a metamaterial synthesized by electron-beam-induced deposition

abstract: The direct writing using a focused electron beam allows for fabricating truly three-dimensional structures of sub-wavelength dimensions in the visible spectral regime. The resulting sophisticated geometries are perfectly suited for studying light-matter interaction at the nanoscale. Their overall optical response will strongly depend not only on geometry but also on the optical properties of the deposited material. In case of the typically used metal-organic precursors, the deposits show a substructure of metallic nanocrystals embedded in a carbonaceous matrix. Since gold-containing precursor media are especially interesting for optical applications, we experimentally determine the effective permittivity of such an effective material. Our experiment is based on spectroscopic measurements of planar deposits. The retrieved permittivity shows a systematic dependence on the gold particle density and cannot be sufficiently described using the common Maxwell-Garnett approach for effective medium.


title: Hopf and Homoclinic Loop Bifurcations on a DC-DC Boost Converter under a SMC Strategy

abstract: In this paper, a dc-dc boost converter with sliding mode control and washout filter is analysed. This device is modelled as a three-dimensional Filippov system, characterized by the existence of sliding movement and restricted to the switching manifold. The operating point of the boost converter is a pseudo-equilibrium, and it, undergoes a subcritical Hopf bifurcation. Such a bifurcation occurs in the sliding vector field and creates, in this field, an unstable limit cycle. The limit cycle is confined to the switching manifold and disappears when it touches the visible-invisible two-fold point, resulting in a homoclinic loop which itself closes in this two-fold point.


title: Light curves of the latest FUor: Indication of a close binary

abstract: We monitored the recent FUor 2MASS J06593158-0405277 (V960 Mon) since November 2009 at various observatories and multiple wavelengths. After the outburst by nearly 2.9 mag in $r$ around September 2014 the brightness gently fades until April 2015 by nearly 1 mag in $U$ and 0.5 mag in $z$. Thereafter the brightness at $\lambda>5000 \AA$ was constant until June 2015 while the shortest wavelengths ($U, B$) indicate a new rise, similar to that seen for the FUor V2493 Cyg (HBC722). Our near-infrared (NIR) monitoring between December 2014 and April 2015 shows a smaller outburst amplitude ($\sim$2 mag) and a smaller (0.2 $-$ 0.3 mag) post-outburst brightness decline. Optical and NIR color-magnitude diagrams indicate that the brightness decline is caused by growing extinction. The post-outburst light curves are modulated by an oscillating color-neutral pattern with a period of about 17 days and an amplitude declining from $\sim$0.08 mag in October 2014 to $\sim$0.04 mag in May 2015. The properties of the oscillating pattern lead us to suggest the presence of a close binary with eccentric orbit.


title: Maximal $L^2$ regularity for Ornstein-Uhlenbeck equation in convex sets of Banach spaces

abstract: We study the elliptic equation $\lambda u-L^{\Omega}u=f$ in an open convex subset $\Omega$ of an infinite dimensional separable Banach space $X$ endowed with a centered non-degenerate Gaussian measure $\gamma$, where $L^\Omega$ is the Ornstein-Uhlenbeck operator. We prove that for $\lambda>0$ and $f\in L^2(\Omega,\gamma)$ the weak solution $u$ belongs to the Sobolev space $W^{2,2}(\Omega,\gamma)$. Moreover we prove that $u$ satisfies the Neumann boundary condition in the sense of traces at the boundary of $\Omega$. This is done by finite dimensional approximation.


title: Quantum capacitance in monolayers of silicene and related buckled materials

abstract: Silicene and related buckled materials are distinct from both the conventional two dimensional electron gas and the famous graphene due to strong spin orbit coupling and the buckled structure. These materials have potential to overcome limitations encountered for graphene, in particular the zero band gap and weak spin orbit coupling. We present a theoretical realization of quantum capacitance which has advantages over the scattering problems of traditional transport measurements. We derive and discuss quantum capacitance as a function of the Fermi energy and temperature taking into account electron-hole puddles through a Gaussian broadening distribution. Our predicted results are very exciting and pave the way for future spintronic and valleytronic devices.


title: Quasi-hermitian lattices with imaginary zero-range interactions

abstract: We study a general class of PT-symmetric tridiagonal Hamiltonians with purely imaginary interaction terms in the quasi-hermitian representation of quantum mechanics. Our general Hamiltonian encompasses many previously studied lattice models as special cases. We provide numerical results regarding domains of observability and exceptional points, and discuss the possibility of explicit construction of general metric operators (which in turn determine all the physical Hilbert spaces). The condition of computational simplicity for the metrics motivates the introduction of certain one-parametric special cases, which consequently admit closed-form extrapolation patterns of the low-dimensional results.


title: Anomalous transport and thermoelectric performances of CuAgSe compounds

abstract: The copper silver selenide has two phases: the low-temperature semimetal phase ({\alpha}-CuAgSe) and high-temperature phonon-glass superionic phase (\b{eta}-CuAgSe). In this work, the electric transport and thermoelectric properties of the two phases are investigated. It is revealed that the \b{eta}-CuAgSe is a p-type semiconductor and exhibits low thermal conductivity while the {\alpha}-CuAgSe shows metallic conduction with dominant n-type carriers and low electrical resistivity. The thermoelectric figure of merit zT of the polycrystalline \b{eta}-CuAgSe at 623 K is ~0.95, suggesting that superionic CuAgSe can be a promising thermoelectric candidate in the intermediate temperature range.


title: The role of membrane curvature for the wrapping of nanoparticles

abstract: Cellular internalization of nanoparticles requires the full wrapping of the nanoparticles by the cell membrane. This wrapping process can occur spontaneously if the adhesive interactions between the nanoparticles and the membranes are sufficiently strong to compensate for the cost of membrane bending. In this article, we show that the membrane curvature prior to wrapping plays a key role for the wrapping process, besides the size and shape of the nanoparticles that have been investigated in recent years. For membrane segments that initially bulge away from nanoparticles by having a mean curvature of the same sign as the mean curvature of the particle surface, we find strongly stable partially wrapped states that can prevent full wrapping. For membrane segments that initially bulge towards the nanoparticles, in contrast, partially wrapped states can constitute a significant energetic barrier for the wrapping process.


title: Determinantal Sampling Designs

abstract: In this article, recent results about point processes are used in sampling theory. Precisely, we define and study a new class of sampling designs: determinantal sampling designs. The law of such designs is known, and there exists a simple selection algorithm. We compute exactly the variance of linear estimators constructed upon these designs by using the first and second order inclusion probabilities. Moreover, we obtain asymptotic and finite sample theorems. We construct explicitly fixed size determinantal sampling designs with given first order inclusion probabilities. We also address the search of optimal determinantal sampling designs.


title: RBF DQ Method for Solving Nonlinear Differential Equations of Lane Emden type

abstract: The Lane-Emden type equations are employed in the modelling of several phenomena in the areas of mathematical physics and astrophysics . In this paper a new numerical method is applied to investigate some well-known classes of Lane-Emden type equations which are nonlinear ordinary differential equations on the semi-infinite domain. We will apply a mesh-less method based on radial basis function differential quadrature method. In RBFs-DQ the derivative value of function with respect to a point is directly approximated by a linear combination of all functional values in the global domain . The main aim of this method is the determination of weight coefficients. Here we concentrate on Gaussian(GS) as a radial function for approximating the solution of the mentioned equation. The comparison of the results with the other numerical methods shows the efficiency and accuracy of this method.


title: Herschel HIFI Observations of the Sgr A +50 km/s Cloud. Deep Searches for O2 in Emission and Foreground Absorption

abstract: To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).


title: Machine learning for metagenomics: methods and tools

abstract: Owing to the complexity and variability of metagenomic studies, modern machine learning approaches have seen increased usage to answer a variety of question encompassing the full range of metagenomic NGS data analysis. We review here the contribution of machine learning techniques for the field of metagenomics, by presenting known successful approaches in a unified framework. This review focuses on five important metagenomic problems: OTU-clustering, binning, taxonomic profling and assignment, comparative metagenomics and gene prediction. For each of these problems, we identify the most prominent methods, summarize the machine learning approaches used and put them into perspective of similar methods. We conclude our review looking further ahead at the challenge posed by the analysis of interactions within microbial communities and different environments, in a field one could call "integrative metagenomics".


title: A surface of maximal canonical degree

abstract: It has been conjectured that the optimal canonical degree of a minimal surface of general type is 36, from a work in the 70's of Beauville who proved that 36 was an upper bound. The highest canonical degree known for the problem was 16 by Persson in 1978. The purpose of this article is to confirm the conjecture by providing an explicit surface.


title: Commitment and Oblivious Transfer in the Bounded Storage Model with Errors

abstract: The bounded storage model restricts the memory of an adversary in a cryptographic protocol, rather than restricting its computational power, making information theoretically secure protocols feasible. We present the first protocols for commitment and oblivious transfer in the bounded storage model with errors, i.e., the model where the public random sources available to the two parties are not exactly the same, but instead are only required to have a small Hamming distance between themselves. Commitment and oblivious transfer protocols were known previously only for the error-free variant of the bounded storage model, which is harder to realize.


title: Sigma-convergence of semilinear stochastic wave equations

abstract: We address the homogenization of a semilinear hyperbolic stochastic partial differential equation with highly oscillating coefficients, in the context of ergodic algebras with mean value. To achieve our goal, we use a suitable variant of the sigma-convergence concept that takes into account both the random and deterministic behaviours of the phenomenon modelled by the underlying problem. We also provide an appropriate scheme for the approximation of the effective coefficients. To illustrate our approach, we work out some concrete problems such as the periodic homogenization problem, the almost periodic and the asymptotically almost periodic ones, and many more besides.


title: Ratchet transport of a two-dimensional electron gas at cyclotron resonance

abstract: The driving of charge carriers confined in a quantum well lacking the center of space inversion by an alternating electric field leads to the formation of a direct electric current. We develop a microscopic theory of such a quantum ratchet effect for quantum wells subjected to a static magnetic field. We show that the ratchet current emerges for a linearly polarized alternating electric field as well as a rotating electric field and drastically increases at the cyclotron resonance conditions. For the magnetic field tilted with respect to the quantum well normal, the ratchet current contains an additional resonance at the first subharmonic of the cyclotron resonance.


title: $a_0-f_0$ mixing in the Khuri-Treiman equations for $\eta\to 3\pi$

abstract: A reliable determination of the isospin breaking double quark mass ratio from precise experimental data on $\eta\to 3\pi$ decays should be based on the chiral expansion of the amplitude supplemented with a Khuri-Treiman type dispersive treatment of the final-state interactions. We discuss an extension of this formalism which allows to estimate the effects of the $a_0(980)$ and $f_0(980)$ resonances and their mixing on the $\eta\to 3\pi$ amplitudes. Matrix generalisations of the equations describing elastic $\pi\pi$ rescattering with $I=0,\,2$ are introduced which accomodate both $\pi\pi/K\bar{K}$ and $\eta\pi/K\bar{K}$ coupled-channel rescattering. Isospin violation induced by the physical $K^+-K^0$ mass difference and by direct $u-d$ mass difference effects are both accounted for in the dispersive integrals. Numerical solutions are constructed which illustrate how the large resonance effects at 1 GeV propagate down to low energies. They remain small in the physical region of the decay, due to the matching constraints with the NLO chiral amplitude, but they are not negligible and go in the sense of further improving the agreement with experiment for the Dalitz plot parameters.


title: Generating domain walls between topologically ordered phases using quantum double models

abstract: Transitions between different topologically ordered phases have been studied by artificially creating boundaries between these gapped phases and thus studying their effects relating to condensation and tunneling of particles from one phase to the other. In this work we introduce exactly solvable models which are similar to the quantum double models (QDM) of Kitaev where such domain walls are dynamically generated making them a part of the spectrum. These systems have a local symmetry and may or may not have a global symmetry leading to the possibility of a ground state degeneracy arising from both, global symmetry breaking and a topological degeneracy. The domain wall states now separate the different topologically ordered phases belonging to the different sectors controlled by the global symmetry, when present. They have interesting properties including fusion with the deconfined anyons of the topologically ordered phase, to either create new domain walls, or their annihilation. They can also act like a scatterer permuting anyons of the topologically ordered phases on either side of the domain wall. Thus these domain wall states can be thought of as a synthetic scatterer which can be created in any part of the lattice. We show these effects for the simplest case of the $\mathbb{Z}_2$ toric code phase decorated with a global $\mathbb{Z}_2$ symmetry and then make remarks about the case when we have the QDM based on two arbitrary groups $G_1$ and $G_2$ in which case we may no longer have a global symmetry.


title: Quantum Cosmological Perturbations of Multiple Fluids

abstract: The formalism to treat quantization and evolution of cosmological perturbations of multiple fluids is described. We first construct the Lagrangian for both the gravitational and matter parts, providing the necessary relevant variables and momenta leading to the quadratic Hamiltonian describing linear perturbations. The final Hamiltonian is obtained without assuming any equations of motions for the background variables. This general formalism is applied to the special case of two fluids, having in mind the usual radiation and matter mix which made most of our current Universe history. Quantization is achieved using an adiabatic expansion of the basis functions. This allows for an unambiguous definition of a vacuum state up to the given adiabatic order. Using this basis, we show that particle creation is well defined for a suitable choice of vacuum and canonical variables, so that the time evolution of the corresponding quantum fields is unitary. This provides constraints for setting initial conditions for an arbitrary number of fluids and background time evolution. We also show that the common choice of variables for quantization can lead to an ill-defined vacuum definition. Our formalism is not restricted to the case where the coupling between fields is small, but is only required to vary adiabatically with respect to the ultraviolet modes, thus paving the way to consistent descriptions of general models not restricted to single-field (or fluid).


title: Using persistent homology to reveal hidden information in neural data

abstract: We propose a method, based on persistent homology, to uncover topological properties of a priori unknown covariates of neuron activity. Our input data consist of spike train measurements of a set of neurons of interest, a candidate list of the known stimuli that govern neuron activity, and the corresponding state of the animal throughout the experiment performed. Using a generalized linear model for neuron activity and simple assumptions on the effects of the external stimuli, we infer away any contribution to the observed spike trains by the candidate stimuli. Persistent homology then reveals useful information about any further, unknown, covariates.


title: Hitting probabilities of random covering sets in tori and metric spaces

abstract: We provide sharp lower and upper bounds for the Hausdorff dimension of the intersection of a typical random covering set with a fixed analytic set both in Ahlfors regular metric spaces and in the $d$-dimensional torus. In metric spaces, we consider covering sets generated by balls and, in the torus, we deal with general analytic generating sets.


title: Kinematic signatures of AGN feedback in moderately powerful radio galaxies at z~2 observed with SINFONI

abstract: Most successful galaxy formation scenarios now postulate that the intense star formation in massive, high-redshift galaxies during their major growth period was truncated when powerful AGNs launched galaxy-wide outflows of gas that removed large parts of the interstellar medium. The most powerful radio galaxies at z~2 show clear signatures of such winds, but are too rare to be good representatives of a generic phase in the evolution of all massive galaxies at high redshift. Here we present SINFONI imaging spectroscopy of 12 radio galaxies at z~2 that are intermediate between the most powerful radio and vigorous starburst galaxies in radio power, and common enough to represent a generic phase in the early evolution of massive galaxies. The kinematic properties are diverse, with regular velocity gradients with amplitudes of Delta v=200-400 km s^-1 as in rotating disks as well as irregular kinematics with multiple velocity jumps of a few 100 km s^-1. Line widths are generally high, typically around FWHM=800 km s^-1, consistent with wind velocities in hydrodynamic models. A broad H-alpha line in one target implies a black hole mass of a few 10^9 M$_sun. The ratio of line widths, sigma, to bulk velocity, v, is so large that even the gas in galaxies with regular velocity fields is unlikely to be gravitationally bound. It is unclear, however, whether the large line widths are due to turbulence or unresolved, local outflows as are sometimes observed at low redshifts. Comparison of the kinetic energy with the energy supply from the AGN through jet and radiation pressure suggests that the radio source still plays a dominant role for feedback, consistent with low-redshift radio-loud quasars.


title: Molecular Anisotropic Magnetoresistance

abstract: Using density functional theory calculations, we demonstrate that the effect of anisotropic magnetoresistance (AMR) can be enhanced by orders of magnitude with respect to conventional bulk ferromagnets in junctions containing molecules sandwiched between ferromagnetic leads. We study ballistic transport in metal-benzene complexes contacted by $3d$ transition-metal wires. We show that the gigantic AMR can arise from spin-orbit coupling effects in the leads, drastically enhanced by orbital-symmetry filtering properties of the molecules. We further discuss how this molecular anisotropic magnetoresistance (MAMR) can be tuned by proper choice of materials and their electronic properties.


title: Asymmetric acoustic propagation of wave packets via the self-demodulation effect

abstract: This article presents the experimental characterization of nonreciprocal elastic wave transmission in a single-mode elastic waveguide. This asymmetric system is obtained by coupling a selection layer with a conversion layer: the selection component is provided by a phononic crystal, while the conversion is achieved by a nonlinear self-demodulation effect in a 3D unconsolidated granular medium. A quantitative experimental study of this acoustic rectifier indicates a high rectifying ratio, up to $10^6$, with wide band (10 kHz) and an audible effect. Moreover, this system allows for wave-packet rectification and extends the future applications of asymmetric systems.


title: Centralized Power Control in Cognitive Radio Networks Using Modulation and Coding Classification Feedback

abstract: In this paper, a centralized Power Control (PC) scheme and an interference channel learning method are jointly tackled to allow a Cognitive Radio Network (CRN) access to the frequency band of a Primary User (PU) operating based on an Adaptive Coding and Modulation (ACM) protocol. The learning process enabler is a cooperative Modulation and Coding Classification (MCC) technique which estimates the Modulation and Coding scheme (MCS) of the PU. Due to the lack of cooperation between the PU and the CRN, the CRN exploits this multilevel MCC sensing feedback as implicit channel state information (CSI) of the PU link in order to constantly monitor the impact of the aggregated interference it causes. In this paper, an algorithm is developed for maximizing the CRN throughput (the PC optimization objective) and simultaneously learning how to mitigate PU interference (the optimization problem constraint) by using only the MCC information. Ideal approaches for this problem setting with high convergence rate are the cutting plane methods (CPM). Here, we focus on the analytic center cutting plane method (ACCPM) and the center of gravity cutting plane method (CGCPM) whose effectiveness in the proposed simultaneous PC and interference channel learning algorithm is demonstrated through numerical simulations.


title: Modeling and Scaling of Hysteresis in Magnetic Materials. Frequency, Pick of Induction and Temperature Dependence

abstract: Recently introduced model of magnetic hysteresis was extended into set of the following features: frequency, pick of induction and temperature of specimen. Group theoretical classification of hysteresis loops' sets is presented. An effect analogous to the Zeeman splitting has been revealed in the set of the all hysteresis loops.


title: Vestoids, Part II: The basaltic nature and HED meteorite analogs for eight Vp-type asteroids and their associations with (4) Vesta

abstract: Improving constraints on the abundance of basaltic asteroids in the main asteroid belt is necessary for better understanding the thermal and collisional environment in the early solar system, for more rigorously identifying the genetic family for (4) Vesta, for determining the effectiveness of Yarkovsky/YORP in dispersing asteroid families, and for better quantifying the population of basaltic asteroids in the outer main belt (a greater than 2.5 AU) that are likely unrelated to (4) Vesta. NIR spectral observations in this work were obtained for the Vp-type asteroids (2011) Veteraniya, (5875) Kuga, (8149) Ruff, (9147) Kourakuen, (9553) Colas, (15237) 1988 RL6, (31414) Rotaryusa, and (32940) 1995 UW4 during August and September 2014 utilizing the SpeX spectrograph at the NASA Infrared Telescope Facility (IRTF), Mauna Kea, Hawaii. Spectral band parameter (band centers, Band Area Ratios) and mineralogical analysis (pyroxene chemistry) for each average asteroid NIR reflectance spectrum suggests a howardite-eucrite-diogenite (HED) meteorite analog for each asteroid. (5875) Kuga is most closely associated with the eucrite meteorites, (31414) Rotaryusa is most closely associated with the diogenites, and the remaining other six asteroids are most closely associated with the howardite meteorites. Along with orbital locations in the inner main belt and in the vicinity of (4) Vesta, the existing evidence suggests that these eight Vp-type asteroids are also likely Vestoids.


title: Replica-symmetric approach to the typical eigenvalue fluctuations of Gaussian random matrices

abstract: We discuss an approach to compute the first and second moments of the number of eigenvalues $I_N$ that lie in an arbitrary interval of the real line for $N \times N$ Gaussian random matrices. The method combines the standard replica-symmetric theory with a perturbative expansion of the saddle-point action up to $O(1/N)$ ($N \gg 1$), leading to the correct logarithmic scaling of the variance $\langle I_{N}^{2} \rangle - \langle I_N \rangle^2 = O(\ln N)$ as well as to an analytical expression for the $O(1/N)$ correction to the average $\langle I_N \rangle/N$. Standard results for the number variance at the local scaling regime are recovered in the limit of a vanishing interval. The limitations of the replica-symmetric method are unveiled by comparing our results with those derived through exact methods. The present work represents an important step to study the fluctuations of $I_N$ in non-invariant random matrix ensembles, where the joint distribution of eigenvalues is not known.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP


title: Cross-platform validation and analysis environment for particle physics

abstract: A multi-platform validation and analysis framework for public Monte Carlo simulation for high-energy particle collisions is discussed. The front-end of this framework uses the Python programming language, while the back-end is written in Java, which provides a multi-platform environment that can be run from a web browser and can easily be deployed at the grid sites. The analysis package includes all major software tools used in high-energy physics, such as Lorentz vectors, jet algorithms, histogram packages, graphic canvases, and tools for providing data access. This multi-platform software suite, designed to minimize OS-specific maintenance and deployment time, is used for online validation of Monte Carlo event samples through a web interface.


title: Towards Accurate Modeling of Moving Contact Lines

abstract: A main challenge in numerical simulations of moving contact line problems is that the adherence, or no-slip boundary condition leads to a non-integrable stress singularity at the contact line. In this report we perform the first steps in developing the macroscopic part of an accurate multiscale model for a moving contact line problem in two space dimensions. We assume that a micro model has been used to determine a relation between the contact angle and the contact line velocity. An intermediate region is introduced where an analytical expression for the velocity exists. This expression is used to implement boundary conditions for the moving contact line at a macroscopic scale, along a fictitious boundary located a small distance away from the physical boundary. Model problems where the shape of the interface is constant thought the simulation are introduced. For these problems, experiments show that the errors in the resulting contact line velocities converge with the grid size $h$ at a rate of convergence $p\approx 2$.


title: Bell's Nonlocality Can be Tested through Einstein-Podolsky-Rosen Steering

abstract: Quantum nonlocality has recently been classified into three distinct types: quantum entanglement, Einstein-Podolsky-Rosen (EPR) steering, and Bell's nonlocality. Experimentally Bell's nonlocality is usually tested by quantum violation of the Clause-Horne-Shimony-Holt (CHSH) inequality in the two-qubit system. Bell's nonlocality is the strongest type of nonlocality, also due this reason Bell-test experiments have encountered both the locality loophole and the detection loophole for a very long time. As a weaker nonlocality, EPR steering naturally escapes from the locality loophole and is correspondingly easier to be demonstrated without the detection loophole. In this work, we trigger an extraordinary approach to investigate Bell's nonlocality, which is strongly based on the EPR steering. We present a theorem, showing that for any two-qubit state $\tau$, if its mapped state $\rho$ is EPR steerable, then the state $\tau$ must be Bell nonlocal. The result not only pinpoints a deep connection between EPR steering and Bell's nonlocality, but also sheds a new light to realize a loophole-free Bell-test experiment (without the CHSH inequality) through the violation of steering inequality.


title: Vector-valued spectra of Banach algebra valued continuous functions

abstract: Given a compact space $X$, a commutative Banach algebra $A$, and an $A$-valued function algebra $\mathscr{A}$ on $X$, the notions of vector-valued spectrum of functions $f\in\mathscr{A}$ are discussed. The $A$-valued spectrum $\vec{SP}_A(f)$ of every $f\in\mathscr{A}$ is defined in such a way that $f(X) \subset \vec{SP}_A(f)$. Utilizing the $A$-characters introduced in (M. Abtahi, \textit{Vector-valued characters on vector-valued function algebras}, \texttt{arXiv:1509.09215 [math.FA]}), it is proved that $\vec{SP}_A(f) = \{\Psi(f):\text{$\Psi$ is an $A$-character of $\mathscr{A}$}\}$. For the so-called natural $A$-valued function algebras, such as $C(X,A)$ and $Lip(X,A)$, we see that $\vec{SP}_A(f)=f(X)$. When $A = \mathbb{C}$, Banach $A$-valued function algebras reduce to Banach function algebras, $A$-characters reduce to characters, and $A$-valued spectrums reduce to usual spectrums.


title: Metallicity determination of M dwarfs - High-resolution IR spectroscopy

abstract: Context. Several new techniques to determine the metallicity of M dwarfs with better precision have been developed over the last decades. However, most of these studies were based on empirical methods. In order to enable detailed abundance analysis, standard methods established for warmer solar-like stars, i.e. model-dependent methods using fitting of synthetic spectra, still need to be used. Aims. In this work we continue the reliability confirmation and development of metallicity determinations of M dwarfs using high- resolution infrared spectra. The reliability was confirmed though analysis of M dwarfs in four binary systems with FGK dwarf companions and by comparison with previous optical studies of the FGK dwarfs. Methods. The metallicity determination was based on spectra taken in the J band (1.1-1.4 {\mu}m) with the CRIRES spectrograph. In this part of the infrared, the density of stellar molecular lines is limited, reducing the amount of blends with atomic lines enabling an accurate continuum placement. Lines of several atomic species were used to determine the stellar metallicity. Results. All binaries show excellent agreement between the derived metallicity of the M dwarf and its binary companion. Our results are also in good agreement with values found in the literature. Furthermore, we propose an alternative way to determine the effective temperature of M dwarfs of spectral types later than M2 through synthetic spectral fitting of the FeH lines in our observed spectra. Conclusions. We have confirmed that a reliable metallicity determination of M dwarfs can be achieved using high-resolution infrared spectroscopy. We also note that metallicites obtained with photometric metallicity calibrations available for M dwarfs only partly agree with the results we obtain from high-resolution spectroscopy.


title: A long-period Cepheid variable in the starburst cluster VdBH222

abstract: Galactic starburst clusters play a twin role in astrophysics, serving as laboratories for the study of stellar physics and also delineating the structure and recent star formation history of the Milky Way. In order to exploit these opportunities we have undertaken a multi-epoch spectroscopic survey of the red supergiant dominated young massive clusters thought to be present at both near and far ends of the Galactic Bar. Significant spectroscopic variability suggestive of radial pulsations was found for the yellow supergiant VdBH 222 #505. Follow-up photometric investigations revealed modulation with a period of ~23.325d; both timescale and pulsational profile are consistent with a Cepheid classification. As a consequence #505 may be recognised as one of the longest period Galactic cluster Cepheids identified to date and hence of considerable use in constraining the bright end of the period/luminosity relation at solar metallicities. In conjunction with extant photometry we infer a distance of ~6kpc for VdBH222 and an age of ~20Myr. This results in a moderate reduction in both integrated cluster mass (~2x10^4Msun) and the initial stellar masses of the evolved cluster members (~10Msun). As such, VdBH222 becomes an excellent test-bed for studying the properties of some of the lowest mass stars observed to undergo type-II supernovae. Moreover, the distance is in tension with a location of VdBH 222 at the far end of the Galactic Bar. Instead a birthsite in the near 3kpc arm is suggested; providing compelling evidence of extensive recent star formation in a region of the inner Milky Way which has hitherto been thought to be devoid of such activity.


title: A universal scaling law for the evolution of granular gases

abstract: Dry, freely evolving granular materials in a dilute gaseous state coalesce into dense clusters only due to dissipative interactions. This clustering transition is important for a number of problems ranging from geophysics to cosmology. Here we show that the evolution of a dilute, freely cooling granular gas is determined in a universal way by the ratio of inertial flow and thermal velocities, that is, the Mach number. Theoretical calculations and direct numerical simulations of the granular Navier--Stokes equations show that irrespective of the coefficient of restitution, density or initial velocity distribution, the density fluctuations follow a universal quadratic dependence on the system's Mach number. We find that the clustering exhibits a scale-free dynamics but the clustered state becomes observable when the Mach number is approximately of $\mathcal{O}(1)$. Our results provide a method to determine the age of a granular gas and predict the macroscopic appearance of clusters.


title: Bangs and Meteors from the Quiet Comet 15P/Finlay

abstract: Jupiter-family comet 15P/Finlay has been reportedly quiet in activity for over a century but has harbored two outbursts during its 2014/2015 perihelion passage. Here we present an analysis of these two outbursts using a set of cometary observations. The outbursts took place between 2014 Dec. 15.4--16.0 UT and 2015 Jan. 15.5--16.0 UT as constrained by ground-based and spacecraft observations. We find a characteristic ejection speed of $V_0=300$ to $650 \mathrm{m \cdot s^{-1}}$ for the ejecta of the first outburst and $V_0=550$ to $750 \mathrm{m \cdot s^{-1}}$ for that of the second outburst using a Monte Carlo dust model. The mass of the ejecta is calculated to be $M_\mathrm{d}=2$ to $3\times10^5 \mathrm{kg}$ for the first outburst and $M_\mathrm{d}=4$ to $5\times10^5 \mathrm{kg}$ for the second outburst, corresponds to less than $10^{-7}$ of the nucleus mass. The specific energy of the two outbursts is found to be $0.3$ to $2\times10^5 \mathrm{J \cdot kg^{-1}}$. We also revisit the long-standing puzzle of the non-detection of the hypothetical Finlayid meteor shower by performing a cued search using the 13-year data from the Canadian Meteor Orbit Radar, which does not reveal any positives. The Earth will pass the 2014/2015 outburst ejecta around 2021 Oct. 6 at 22 h UT to Oct. 7 at 1 h UT, with a chance for some significant meteor activity in the radio range, which may provide further clues to the Finlayid puzzle. A southerly radiant in the constellation of Ara will favor the observers in the southern tip of Africa.


title: A 'Gibbs-Newton' Technique for Enhanced Inference of Multivariate Polya Parameters and Topic Models

abstract: Hyper-parameters play a major role in the learning and inference process of latent Dirichlet allocation (LDA). In order to begin the LDA latent variables learning process, these hyper-parameters values need to be pre-determined. We propose an extension for LDA that we call 'Latent Dirichlet allocation Gibbs Newton' (LDA-GN), which places non-informative priors over these hyper-parameters and uses Gibbs sampling to learn appropriate values for them. At the heart of LDA-GN is our proposed 'Gibbs-Newton' algorithm, which is a new technique for learning the parameters of multivariate Polya distributions. We report Gibbs-Newton performance results compared with two prominent existing approaches to the latter task: Minka's fixed-point iteration method and the Moments method. We then evaluate LDA-GN in two ways: (i) by comparing it with standard LDA in terms of the ability of the resulting topic models to generalize to unseen documents; (ii) by comparing it with standard LDA in its performance on a binary classification task.


title: A topological study of gravity waves generated by moving bodies using the method of steepest descents

abstract: The standard analytical approach for studying gravity free-surface waves generated by a moving body often relies upon a linearization of the physical geometry, where the body is considered asymptotically small in one or several of its dimensions. In this paper, a methodology that avoids any such geometrical simplification is presented for the case of flows at low speeds. The approach is made possible through a reduction of the water-wave equations to a complex-valued integral equation that can be studied using the method of steepest descents. The main result is a theory that establishes a correspondence between a given physical flow geometry, with the topology of the Riemann surface formed by the steepest descent paths. Then, when a geometrical feature of the body is modified, a corresponding change to the Riemann surface is observed, and the resultant effects to the water waves can be derived. This visual procedure is demonstrated for the case of two-dimensional free-surface flow past a surface-piercing ship and over an angled step in a channel.


title: Do the $P_c^+$ Pentaquarks Have Strange Siblings?

abstract: The recent LHCb discovery of states $P_c^+(4380)$, $P_c^+(4450)$, believed to be $c\bar c uud$ pentaquark resonances, begs the question of whether equivalent states with $c\bar c \to s\bar s$ exist, and how they might be produced. The precise analogue to the $P_c^+$ discovery channel $\Lambda_b \to J/\psi \, K^- \! p$, namely, $\Lambda_c \to \phi \pi^0 \! p$, is feasible for this study and indeed is less Cabibbo-suppressed, although its limited phase space suggests that evidence of a $s\bar s uud$ resonance $P_s^+$ would be confined to the kinematic endpoint region.