title: Real-Time Bidding by Reinforcement Learning in Display Advertising

abstract: The majority of online display ads are served through real-time bidding (RTB) --- each ad display impression is auctioned off in real-time when it is just being generated from a user visit. To place an ad automatically and optimally, it is critical for advertisers to devise a learning algorithm to cleverly bid an ad impression in real-time. Most previous works consider the bid decision as a static optimization problem of either treating the value of each impression independently or setting a bid price to each segment of ad volume. However, the bidding for a given ad campaign would repeatedly happen during its life span before the budget runs out. As such, each bid is strategically correlated by the constrained budget and the overall effectiveness of the campaign (e.g., the rewards from generated clicks), which is only observed after the campaign has completed. Thus, it is of great interest to devise an optimal bidding strategy sequentially so that the campaign budget can be dynamically allocated across all the available impressions on the basis of both the immediate and future rewards. In this paper, we formulate the bid decision process as a reinforcement learning problem, where the state space is represented by the auction information and the campaign's real-time parameters, while an action is the bid price to set. By modeling the state transition via auction competition, we build a Markov Decision Process framework for learning the optimal bidding policy to optimize the advertising performance in the dynamic real-time bidding environment. Furthermore, the scalability problem from the large real-world auction volume and campaign budget is well handled by state value approximation using neural networks.


title: Axial and polar gravitational wave equations in a de Sitter expanding universe by Laplace transform

abstract: In this paper we study the propagation in a de Sitter universe of gravitational waves generated by perturbating some unspecified spherical astrophysical object in the frequencies domain. We obtain the axial and polar perturbation equations in a cosmological de Sitter universe in the usual comoving coordinates, the coordinates we occupy in our galaxy. We write down the relevant equations in terms of Laplace transform with respect to the comoving time $t$ instead of the usual Fourier one that is no longer available in a cosmological context. Both axial and polar perturbation equations are expressed in terms of a non trivial mixture of retarded-advanced metric coefficients with respect to the Laplace parameter $s$ (complex translation). The axial case is studied in more detail. In particular, the axial perturbations can be reduced to a master linear second-order differential equation in terms of the Regge-Wheeler function $Z$ where a coupling with a retarded $Z$ with respect to the cosmological time $t$ is present. It is shown that a de Sitter expanding universe can change the frequency $\omega$ of a gravitational wave as perceived by a comoving observer. The polar equations are much more involved. Nevertheless, we show that also the polar perturbations can be expressed in terms of four independent integrable differential equations.


title: Recurrent fast radio bursts from collisions of neutron stars in the evolved stellar clusters

abstract: We propose the model describing the observed multiple fast radio bursts due to the close encounters and collisions of neutron stars in the central clusters of the evolved galactic nuclei. The subsystem of neutron star cluster may originate in the dense galactic nucleus evolutionary in the combined processes of stellar and dynamical evolution. The neutron stars in the compact cluster can produce the short living binaries with the highly eccentric orbits, and finally collide after several orbital revolutions. Fast radio bursts may be produced during the close periastron approach and at the process of the final binary merging. In the sufficiently dense star cluster the neutron stars collisions can be very frequent. Therefore, this model can explain in principle the observed recurrent (multiple) fast radio bursts, analogous to the observed ones from the source FRB 121102. Among the possible observational signatures of the proposed model may be the registration of the gravitational wave bursts by the laser interferometers LIGO/VIRGO or by the next generation of gravitational wave detectors.


title: Fresnel coefficients and Fabry-Perot formula for spatially dispersive metallic layers

abstract: The repulsion between free electrons inside a metal makes its optical response spatially dispersive, so that it is not described by Drude's model but by a hydrodynamic model. We give here fully analytic results for a metallic slab in this framework, thanks to a two-modes cavity formalism leading to a Fabry-Perot formula, and show that a simplification can be made that preserves the accuracy of the results while allowing much simpler analytic expressions. For metallic layers thicker than 2.7 nm modified Fresnel coefficients can actually be used to accurately predict the response of any multilayer with spatially dispersive metals (for reflection, transmission or the guided modes). Finally, this explains why adding a small dielectric layer[Y. Luo et al., Phys. Rev. Lett. 111, 093901 (2013)] allows to reproduce the effects of nonlocality in many cases, and especially for multilayers.


title: Dynamic Complexity under Definable Changes

abstract: This paper studies dynamic complexity under definable change operations in the DynFO framework by Patnaik and Immerman. It is shown that for changes definable by parameter-free first-order formulas, all (uniform) $AC^1$ queries can be maintained by first-order dynamic programs. Furthermore, many maintenance results for single-tuple changes are extended to more powerful change operations: (1) The reachability query for undirected graphs is first-order maintainable under single tuple changes and first-order defined insertions, likewise the reachability query for directed acyclic graphs under quantifier-free insertions. (2) Context-free languages are first-order maintainable under $\Sigma_1$-defined changes. These results are complemented by several inexpressibility results, for example, that the reachability query cannot be maintained by quantifier-free programs under definable, quantifier-free deletions.


title: ATK-ForceField: A New Generation Molecular Dynamics Software Package

abstract: ATK-ForceField is a software package for atomistic simulations using classical interatomic potentials. It is implemented as a part of the Atomistix ToolKit (ATK), which is a Python programming environment that makes it easy to create and analyze both standard and highly customized simulations. This paper will focus on the atomic interaction potentials, molecular dynamics, and geometry optimization features of the software, however, many more advanced modeling features are available. The implementation details of these algorithms and their computational performance will be shown. We present three illustrative examples of the types of calculations that are possible with ATK-ForceField: modeling thermal transport properties in a silicon germanium crystal, vapor deposition of selenium molecules on a selenium surface, and a simulation of creep in a copper polycrystal.


title: Network Topology Modulation for Energy and Data Transmission in Internet of Magneto-Inductive Things

abstract: Internet-of-things (IoT) architectures connecting a massive number of heterogeneous devices need energy efficient, low hardware complexity, low cost, simple and secure mechanisms to realize communication among devices. One of the emerging schemes is to realize simultaneous wireless information and power transfer (SWIPT) in an energy harvesting network. Radio frequency (RF) solutions require special hardware and modulation methods for RF to direct current (DC) conversion and optimized operation to achieve SWIPT which are currently in an immature phase. On the other hand, magneto-inductive (MI) communication transceivers are intrinsically energy harvesting with potential for SWIPT in an efficient manner. In this article, novel modulation and demodulation mechanisms are presented in a combined framework with multiple-access channel (MAC) communication and wireless power transmission. The network topology of power transmitting active coils in a transceiver composed of a grid of coils is changed as a novel method to transmit information. Practical demodulation schemes are formulated and numerically simulated for two-user MAC topology of small size coils. The transceivers are suitable to attach to everyday objects to realize reliable local area network (LAN) communication performances with tens of meters communication ranges. The designed scheme is promising for future IoT applications requiring SWIPT with energy efficient, low cost, low power and low hardware complexity solutions.


title: Size dependence of the surface tension of a free surface of an isotropic fluid

abstract: We report on the size dependence of the surface tension of a free surface of an isotropic fluid. The size dependence of the surface tension is evaluated based on the Gibbs-Tolman-Koenig-Buff equation for positive and negative values of curvatures and the Tolman lengths. For all combinations of positive and negative signs of curvature and the Tolman length, we succeed to have a continuous function, avoiding the existing discontinuity at zero curvature (flat interfaces). As an example, a water droplet in the thermodynamical equilibrium with the vapor is analyzed in detail. The size dependence of the surface tension and the Tolman length are evaluated with the use of experimental data of the International Association for the Properties of Water and Steam. The evaluated Tolman length of our approach is in good agreement with molecular dynamics and experimental data


title: Ru-doping on iron based pnictides: the "unfolded" dominant role of structural effects for superconductivity

abstract: We present an ab-initio study of Ru substitution in two different compounds, BaFe2As2 and LaFeAsO, pure and F-doped. Despite the many similarities among them, Ru substitution has very different effects on these compounds. By means of an unfolding technique, which allows us to trace back the electronic states into the primitive cell of the pure compounds, we are able to disentangle the effects brought by the local structural deformations and by the impurity potential to the states at the Fermi level. Our results are compared with available experiments and show: i) satisfying agreement of the calculated electronic properties with experiments, confirming the presence of a magnetic order on a short range scale; ii) Fermi surfaces strongly dependent on the internal structural parameters, more than on the impurity potential. These results enter a widely discussed field in the literature and provide a better understanding of the role of Ru in iron pnictides: although isovalent to Fe, the Ru-Fe substitution leads to changes in the band structure at the Fermi level mainly related to local structural modifications.


title: High performance graphene oxide-based humidity sensor integrated on a photonic crystal cavity

abstract: We report a high performance relative humidity (RH) microsensor by coating a few-layer graphene oxide (GO) flake over a photonic crystal (PC) cavity. Since the GO layer has a high water-activity and interacts with the evanescent cavity mode strongly, the exposure of GO-PC cavity in varied humidity levels results in significant resonant wavelength shifts, showing a slope of 0.68 nm/%RH in the range of 60%-85% RH. By interrogating the power variation of the cavity reflection, the microsensor presents an ultrahigh sensitivity exceeding 3.9 dB/%RH. Relying on the unimpeded permeation of water molecules through the GO interlayers and microscale distribution of the cavity mode, the integrated sensor has a response time less than 100 ms, which promises successful measurements of human breathing. Combining with the ease of fabrication, this high performance RH sensor has potential applications requiring optical access, device compactness, and fast dynamic response.


title: Orbital-driven two-dome superconducting phases in iron-based superconductors

abstract: Recent several experiments revealed that novel bipartite magnetic/superconducting phases widely exist in iron pnictides and chalcogenides. Nevertheless, the origin of the two-dome superconducting phases in iron-based compounds still remains unclear. Here we theoretically investigated the electronic structures, magnetic and superconducting properties of three representative iron-based systems, i.e. LaFeAsO1-xHx, LaFeAs1-xPxO and KFe2As2. We found that in addition to the degenerate in-plane anisotropic xz/yz orbitals, the quasi-degenerate in-plane isotropic orbitals drive these systems entering into the second parent phase. Moreover, the second superconducting phase is contributed by the isotropic orbitals rather than the anisotropic ones in the first superconducting phase, indicating an orbital-selective pairing state. These results imply an orbital-driven mechanism and shed light on the understanding of the two-dome magnetic/superconducting phases in iron-based compounds.


title: Stellarator bootstrap current and plasma flow velocity at low collisionality

abstract: The bootstrap current and flow velocity of a low-collisionality stellarator plasma are calculated. As far as possible, the analysis is carried out in a uniform way across all low-collisionality regimes in general stellarator geometry, assuming only that the confinement is good enough that the plasma is approximately in local thermodynamic equilibrium. It is found that conventional expressions for the ion flow speed and bootstrap current in the low-collisionality limit are accurate only in the $1/\nu$-collisionality regime and need to be modified in the $\sqrt{\nu}$-regime. The correction due to finite collisionality is also discussed and is found to scale as $\nu^{2/5}$.


title: Untwisting two-way transducers in elementary time

abstract: Functional transductions realized by two-way transducers (equivalently, by streaming transducers and by MSO transductions) are the natural and standard notion of "regular" mappings from words to words. It was shown recently (LICS'13) that it is decidable if such a transduction can be implemented by some one-way transducer, but the given algorithm has non-elementary complexity. We provide an algorithm of different flavor solving the above question, that has double exponential space complexity. We further apply our technique to decide whether the transduction realized by a two-way transducer can be implemented by a sweeping transducer, with either known or unknown number of passes.


title: Size-dependent frequency bands in the ferromagnetic resonance of a Fe-nanocube

abstract: Using full micromagnetic simulations we calculate the spectra of ferromagnetic resonance (FMR) for an iron (core-shell) nanocube and show that the FMR characteristics are strongly size dependent. For instance, for a $40~$nm it is found that, in contrast to a macrospin picture, the spectrum of the iron nanocube possesses two bands centered around $0.4~$T and $\approx 0.1~$T. The peaks originate from the surface anisotropy induced by the strong demagnetizing fields (DMFs) of iron. Further simulations reveal that for $\approx 20~$nm nanocubes the macrospin model becomes viable. Above $40~$nm we find a broad band for FMR absorption. Our results point to possible interpretations of existing FMR experimental observations for the system studied here.


title: Splashing of impacting drops

abstract: We investigate the impact velocity beyond which the ejection of smaller droplets from the main droplet (splashing) occurs for droplets impacting a smooth surface. We examine its dependence on the surface wetting properties and droplet surface tension. We show that the splashing velocity is independent of the wetting properties of the surface, but increases roughly linearly with increasing surface tension of the liquid. A preexisting splashing model is considered that predicts the splashing velocity by incorporating the air viscosity. The model is consistent with our data, but the calculation is complex. To address this issue, we propose a simplification, assuming atmospheric conditions and low viscosities, and show that the simplification gives an equally good prediction of the splashing velocity from a simple analytical formula.


title: Essential surfaces in graph pairs

abstract: A well known question of Gromov asks whether every one-ended hyperbolic group $\Gamma$ has a surface subgroup. We give a positive answer when $\Gamma$ is the fundamental group of a graph of free groups with cyclic edge groups. As a result, Gromov's question is reduced (modulo a technical assumption on 2-torsion) to the case when $\Gamma$ is rigid. We also find surface subgroups in limit groups. It follows that a limit group with the same profinite completion as a free group must in fact be free, which answers a question of Remeslennikov in this case.


title: Using red clump stars to correct the Gaia DR1 parallaxes

abstract: Recent results have suggested that there is tension between the Gaia DR1 TGAS distances and the distances obtained using luminosities determined by eclipsing binaries or asteroseismology on red giant stars. We use the $K_s$-band luminosities of red clump stars, identified and characterized by asteroseismology, to make independent distance estimates. Our results suggest that Gaia TGAS distances contain a systematic error that decreases with increasing distance. We propose a correction to mitigate this offset as a function of parallax that is valid for the Kepler field and values of parallax that are less than ${\sim} 1.6 \rm \, mas$. For parallaxes greater than this we find agreement with previously published values. We note that the TGAS distances to the red clump stars of the open cluster M67 show a high level of disagreement that is difficult to correct for.


title: Time-domain Studies of M31

abstract: M31, our closest neighboring galaxy, is a stepping stone to studies of stellar evolution, star formation, galaxy evolution, and cosmology. However, due to the difficulties of performing photometry in such crowded fields and the lack of wide cameras to encompass the entire galaxy, there has not been a complete census of the stellar contents of M31. The advent of wide-field camera provides us a unique opportunity to have a complete view of our neighboring galaxy, enabling an inventory of its variable content. We present a review of recent progresses of wide-field, high cadence surveys of M31, covering different population of variables and transients. We also outline future studies enabled by on-going and up-coming facilities.


title: Evolution of piecewise polynomial wave functions

abstract: For a non-relativistic particle, we consider the evolution of wave functions that consist of polynomial segments, usually joined smoothly together. These spline wave functions are compact (that is, they are initially zero outside a finite region), but they immediately extend over all available space as they evolve. The simplest splines are the square and triangular wave functions in one dimension, but very complicated splines have been used in physics. In general the evolution of such spline wave functions can be expressed in terms of antiderivatives of the propagator; in the case of a free particle or an oscillator, all the evolutions are expressed exactly in terms of Fresnel integrals. Some extensions of these methods to two and three dimensions are discussed.


title: Tangle-tree duality in abstract separation systems

abstract: We prove a general width duality theorem for combinatorial structures with well-defined notions of cohesion and separation. These might be graphs and matroids, but can be much more general or quite different. The theorem asserts a duality between the existence of high cohesiveness somewhere local and a global overall tree structure. We describe cohesive substructures in a unified way in the format of tangles: as orientations of low-order separations satisfying certain consistency axioms. These axioms can be expressed without reference to the underlying structure, such as a graph or matroid, but just in terms of the poset of the separations themselves. This makes it possible to identify tangles, and apply our tangle-tree duality theorem, in very diverse settings. Our result implies all the classical duality theorems for width parameters in graph minor theory, such as path-width, tree-width, branch-width or rank-width. It yields new, tangle-type, duality theorems for tree-width and path-width. It implies the existence of width parameters dual to cohesive substructures such as $k$-blocks, edge-tangles, or given subsets of tangles, for which no width duality theorems were previously known. Abstract separation systems can be found also in structures quite unlike graphs and matroids. For example, our theorem can be applied to image analysis by capturing the regions of an image as tangles of separations defined as natural partitions of its set of pixels. It can be applied in big data contexts by capturing clusters as tangles. It can be applied in the social sciences, e.g. by capturing as tangles the few typical mindsets of individuals found by a survey. It could also be applied in pure mathematics, e.g. to separations of compact manifolds.


title: Noise-enhanced chaos in a weakly coupled GaAs/(Al,Ga)As superlattice

abstract: Noise-enhanced chaos in a doped, weakly coupled GaAs/Al_{0.45}Ga_{0.55}As superlattice has been observed at room temperature in experiments as well as in the results of the simulation of nonlinear transport based on a discrete tunneling model. When external noise is added, both the measured and simulated current-versus-time traces contain irregularly spaced spikes for particular applied voltages, which separate a regime of periodic current oscillations from a region of no current oscillations at all. In the voltage region without current oscillations, the electric field profile consist of a low-field domain near the emitter contact separated by a domain wall consisting of a charge accumulation layer from a high-field regime closer to the collector contact. With increasing noise amplitude, spontaneous chaotic current oscillations appear over a wider bias voltage range. For these bias voltages, the domain boundary between the two electric-field domains becomes unstable, and very small current or voltage fluctuations can trigger the domain boundary to move toward the collector and induce chaotic current spikes. The experimentally observed features are qualitatively very well reproduced by the simulations. Increased noise can consequently enhance chaotic current oscillations in semiconductor superlattices.


title: Heterogeneous domain adaptation: An unsupervised approach

abstract: Domain adaptation leverages the knowledge in one domain - the source domain - to improve learning efficiency in another domain - the target domain. Existing heterogeneous domain adaptation research is relatively well-progressed, but only in situations where the target domain contains at least a few labeled instances. In contrast, heterogeneous domain adaptation with an unlabeled target domain has not been well-studied. To contribute to the research in this emerging field, this paper presents: (1) an unsupervised knowledge transfer theorem that guarantees the correctness of transferring knowledge; and (2) a principal angle-based metric to measure the distance between two pairs of domains: one pair comprises the original source and target domains and the other pair comprises two homogeneous representations of two domains. The theorem and the metric have been implemented in an innovative transfer model, called a Grassmann-Linear monotonic maps-geodesic flow kernel (GLG), that is specifically designed for heterogeneous unsupervised domain adaptation (HeUDA). The linear monotonic maps meet the conditions of the theorem and are used to construct homogeneous representations of the heterogeneous domains. The metric shows the extent to which the homogeneous representations have preserved the information in the original source and target domains. By minimizing the proposed metric, the GLG model learns the homogeneous representations of heterogeneous domains and transfers knowledge through these learned representations via a geodesic flow kernel. To evaluate the model, five public datasets were reorganized into ten HeUDA tasks across three applications: cancer detection, credit assessment, and text classification. The experiments demonstrate that the proposed model delivers superior performance over the existing baselines.


title: An RKHS model for variable selection in functional regression

abstract: A mathematical model for variable selection in functional regression models with scalar response is proposed. By "variable selection" we mean a procedure to replace the whole trajectories of the functional explanatory variables with their values at a finite number of carefully selected instants (or "impact points"). The basic idea of our approach is to use the Reproducing Kernel Hilbert Space (RKHS) associated with the underlying process, instead of the more usual L2[0,1] space, in the definition of the linear model. This turns out to be especially suitable for variable selection purposes, since the finite-dimensional linear model based on the selected "impact points" can be seen as a particular case of the RKHS-based linear functional model. In this framework, we address the consistent estimation of the optimal design of impact points and we check, via simulations and real data examples, the performance of the proposed method.


title: pNRQCD determination of E1 radiative transitions

abstract: This contribution contains the first numerical computation of the complete set of relativistic corrections of relative order $v^{2}$ for electric dipole (E1) transitions in heavy quarkonium; in particular, for the processes $\chi_{bJ}(1P) \to \Upsilon(1S) + \gamma$ with $J=0,\,1,\,2$. We assume that the momentum transfer of the heavy mesons involved in the reactions lies in the weak-coupling regime of the low-energy effective field theory potential non-relativistic QCD (pNRQCD) and thus a full perturbative calculation can be performed.


title: On Centroidal Dynamics and Integrability of Average Angular Velocity

abstract: In the literature on robotics and multibody dynamics, the concept of average angular velocity has received considerable attention in recent years. We address the question of whether the average angular velocity defines an orientation framethat depends only on the current robot configuration and provide a simple algebraic condition to check whether this holds. In the language of geometric mechanics, this condition corresponds to requiring the flatness of the mechanical connection associated to the robotic system. Here, however, we provide both a reinterpretation and a proof of this result accessible to readers with a background in rigid body kinematics and multibody dynamics but not necessarily acquainted with differential geometry, still providing precise links to the geometric mechanics literature. This should help spreading the algebraic condition beyond the scope of geometric mechanics,contributing to a proper utilization and understanding of the concept of average angular velocity.


title: Restriction of representations of metaplectic $GL_{2}(F)$ to tori

abstract: Let $F$ be a non-Archimedean local field. We study the restriction of an irreducible admissible genuine representations of the two fold metaplectic cover $\widetilde{GL}_{2}(F)$ of $GL_{2}(F)$ to the inverse image in $\widetilde{GL}_{2}(F)$ of a maximal torus in $GL_{2}(F)$.


title: A comment of the combinatorics of the vertex operator $\Gamma_{(t|X)}$

abstract: The Jacobi--Trudi identity associates a symmetric function to any integer sequence. Let $\Gamma_{(t|X)}$ be the vertex operator defined by $\Gamma_{(t|X)} s_\alpha =\sum_{n \in \mathbb{Z}} s_{(n,\alpha)} [X] t^n$. We provide a combinatorial proof for the identity $\Gamma_{(t|X)} s_\alpha = \sigma[tX] s_{\alpha}\big[x-1/t\big] $ due to Thibon et al. We include an overview of all the combinatorial ideas behind this beautiful identity, including a combinatorial description for the expansion of $s_{(n,\alpha)} [X] $ in the Schur basis, for any integer value of $n$.


title: Study of $ B_s^0\rightarrow l^+ l^-\gamma $ decays in the nonuniversal Z' model

abstract: We investigate $ B_s^0\rightarrow l^+ l^-\gamma $ decays in a nonuniversal Z' model derived from extension of the standard model (SM). Considering the Z'-mediated flavor changing neutral current (FCNC) effects we calculate the branching ratio and forwardbackward asymmetry for $ B_s^0\rightarrow l^+ l^-\gamma $ decay processes. We compare the obtained results with predictions of the SM and discuss the sensitivity of the observables to Z' boson coupling parameters. We find the branching ratios are enhanced by one order from SM predictions in Z' model scenario. We also observe the variation of forward-backward asymmetry with the Z' boson coupling parameters portrait discrimination between NP effects and SM results.


title: Cluster algebras and symmetrizable matrices

abstract: In this paper, we study combinatorial properties of quasi-Cartan companions defined by the c-vectors of acyclic skew-symmetrizable cluster algebras. In particular, we show that the diagram of any skew-symmetrizable matrix associated with an acyclic cluster algebra has an admissible cut of edges.


title: Gravity Amplitudes as Generalized Double Copies

abstract: Whenever the integrand of a gauge-theory loop amplitude can be arranged into a form where the BCJ duality between color and kinematics is manifest, a corresponding gravity integrand can be obtained simply via the double-copy procedure. However, finding such gauge-theory representations can be challenging, especially at high loop orders. Here we show that we can instead start from generic gauge-theory integrands, where the duality is not manifest, and apply a modified double-copy procedure to obtain gravity integrands that include contact terms generated by violations of dual Jacobi identities. We illustrate this with three-, four- and five-loop examples in N=8 supergravity.


title: Mean field approach to first and second order phase transitions in ionic liquids

abstract: Ionic liquids are solvent-free electrolytes, some of which possess an intriguing self-assembly property. Using a mean-field framework (based on Onsager's relations) we show that bulk nano-structures arise via type-I and II phase transitions (PT), which directly affect the electrical double layer (EDL) structure. Ginzburg-Landau equation is derived and PT are related to temperature, potential and interactions. The type-I PT occurs for geometrically dissimilar anion/cation ratio and, surprisingly, is induced by perturbations on order of thermal fluctuations. Finally, we compare the insights with the decaying charge layers within the EDL, as widely observed in experiments.


title: Accurate formation energies of charged defects in solids: a systematic approach

abstract: Defects on surfaces of semiconductors have a strong effect on their reactivity and catalytic properties. The concentration of different charge states of defects is determined by their formation energies. First-principles calculations are an important tool for computing defect formation energies and for studying the microscopic environment of the defect. The main problem associated with the widely used supercell method in these calculations is the error in the electrostatic energy, which is especially pronounced in calculations that involve surface slabs and 2D materials. We present an internally consistent approach for calculating defect formation energies in inhomogeneous and anisotropic dielectric environments, and demonstrate its applicability to the cases of the positively charged Cl vacancy on the NaCl (100) surface and the negatively charged S vacancy in monolayer MoS2.


title: Magnus expansions and pseudospectra of Master Equations

abstract: New directions in research on master equations are showcased by example. Magnus expansions, time-varying rates, and pseudospectra are highlighted. Exact eigenvalues are found and contrasted with the large errors produced by standard numerical methods in some cases. Isomerisation provides a running example and an illustrative application to chemical kinetics. We also give a brief example of the totally asymmetric exclusion process.


title: Maps on positive definite operators preserving the quantum $\chi_\alpha^2$-divergence

abstract: We describe the structure of all bijective maps on the cone of positive definite operators acting on a finite and at least two-dimensional complex Hilbert space which preserve the quantum $\chi_\alpha^2$-divergence for some $\alpha \in [0,1]$. We prove that any such transformation is necessarily implemented by either a unitary or an antiunitary operator. Similar results concerning maps on the cone of positive semidefinite operators as well as on the set of all density operators are also derived.


title: A Chunk Caching Location and Searching Scheme in Content Centric Networking

abstract: Content Centric Networking (CCN) is a new network infrastructure around content dissemination and retrieval, shift from host addresses to named data. Each CCN router has a cache to store the chunks passed by it. Therefore the caching strategy about chunk placement can greatly affect the whole CCN performance. This paper proposes an implicit coordinate chunk caching location and searching scheme (CLS) in CCN hierarchical infrastructure. In CLS, there is at most one copy of a chunk cached on the path between a server and a leaf router. This copy is pulled down one level towards the leaf router by a request or pushed up one level towards the server by the cache eviction. Thus, it is possible to store more diverse contents in the whole CCN and improve the network performance. Plus, in order to reduce the server workload and file download time, a caching trail of chunk is created to direct the following request where to find the chunk. Extensive test-bed experiments have been performed to evaluate the proposed scheme in terms of a wide range of performance metrics. The results show that the proposed scheme outperforms existing algorithms.


title: A reduced fast construction of polynomial lattice point sets with low weighted star discrepancy

abstract: The weighted star discrepancy is a quantitative measure for the performance of point sets in quasi-Monte Carlo algorithms for numerical integration. We consider polynomial lattice point sets, whose generating vectors can be obtained by a component-by-component construction to ensure a small weighted star discre-pancy. Our aim is to significantly reduce the construction cost of such generating vectors by restricting the size of the set of polynomials from which we select the components of the vectors. To gain this reduction we exploit the fact that the weights of the spaces we consider decay very fast.


title: Modeling and Reasoning About Wireless Networks: A Graph-based Calculus Approach

abstract: We propose a graph-based process calculus for modeling and reasoning about wireless networks with local broadcasts. Graphs are used at syntactical level to describe the topological structures of networks. This calculus is equipped with a reduction semantics and a labelled transition semantics. The former is used to define weak barbed congruence. The latter is used to define a parameterized weak bisimulation emphasizing locations and local broadcasts. We prove that weak bisimilarity implies weak barbed congruence. The potential applications are illustrated by some examples and two case studies.


title: The heavy path approach to Galton-Watson trees with an application to Apollonian networks

abstract: We study the heavy path decomposition of conditional Galton-Watson trees. In a standard Galton-Watson tree conditional on its size $n$, we order all children by their subtree sizes, from large (heavy) to small. A node is marked if it is among the $k$ heaviest nodes among its siblings. Unmarked nodes and their subtrees are removed, leaving only a tree of marked nodes, which we call the $k$-heavy tree. We study various properties of these trees, including their size and the maximal distance from any original node to the $k$-heavy tree. In particular, under some moment condition, the $2$-heavy tree is with high probability larger than $cn$ for some constant $c > 0$, and the maximal distance from the $k$-heavy tree is $O(n^{1/(k+1)})$ in probability. As a consequence, for uniformly random Apollonian networks of size $n$, the expected size of the longest simple path is $\Omega(n)$.


title: Why It Takes So Long to Connect to a WiFi Access Point

abstract: Today's WiFi networks deliver a large fraction of traffic. However, the performance and quality of WiFi networks are still far from satisfactory. Among many popular quality metrics (throughput, latency), the probability of successfully connecting to WiFi APs and the time cost of the WiFi connection set-up process are the two of the most critical metrics that affect WiFi users' experience. To understand the WiFi connection set-up process in real-world settings, we carry out measurement studies on $5$ million mobile users from $4$ representative cities associating with $7$ million APs in $0.4$ billion WiFi sessions, collected from a mobile "WiFi Manager" App that tops the Android/iOS App market. To the best of our knowledge, we are the first to do such large scale study on: how large the WiFi connection set-up time cost is, what factors affect the WiFi connection set-up process, and what can be done to reduce the WiFi connection set-up time cost. Based on the measurement analysis, we develop a machine learning based AP selection strategy that can significantly improve WiFi connection set-up performance, against the conventional strategy purely based on signal strength, by reducing the connection set-up failures from $33\%$ to $3.6\%$ and reducing $80\%$ time costs of the connection set-up processes by more than $10$ times.


title: Focused time-lapse inversion of radio and audio magnetotelluric data

abstract: Geoelectrical techniques are widely used to monitor groundwater processes, while surprisingly few studies have considered audio (AMT) and radio (RMT) magnetotellurics for such purposes. In this numerical investigation, we analyze to what extent inversion results based on AMT and RMT monitoring data can be improved by (1) time-lapse difference inversion; (2) incorporation of statistical information about the expected model update (i.e., the model regularization is based on a geostatistical model); (3) using alternative model norms to quantify temporal changes (i.e., approximations of l1 and Cauchy norms using iteratively reweighted least-squares), (4) constraining model updates to predefined ranges (i.e., using Lagrange Multipliers to only allow either increases or decreases of electrical resistivity with respect to background conditions). To do so, we consider a simple illustrative model and a more realistic test case related to seawater intrusion. The results are encouraging and show significant improvements when using time-lapse difference inversion with non l2 model norms. Artifacts that may arise when imposing compactness of regions with temporal changes can be suppressed through inequality constraints to yield models without oscillations outside the true region of temporal changes. Based on these results, we recommend approximate l1-norm solutions as they can resolve both sharp and smooth interfaces within the same model.


title: Effect of loss on the topological features of dimer chain described by the extended Aubry-Andr\'e-Harper model

abstract: By introducing loss to one sublattice of a dimer chain described by the extended Aubry-Andr\'e or Harper (AAH) model, we study the topological features including the edge states, spectrum and winding number of the chain. We find that the parameter region for the system to have real band-gap-closing is increased due to the loss, and the average displacement of the single excitation can still witness the topological features of the chain in the presence of loss. The robustness of the zero energy eigenstate against four kinds of disorders is also examined. A feasible experiment setup based on coupled waveguides to observe the prediction of this paper is proposed.


title: Dynamics of subarcsecond bright dots in the transition region above sunspot and their relation to penumbral micro-jets

abstract: Recent high-resolution observations reveal that subarcsecond bright dots (BDs) with sub-minute lifetimes appears ubiquitously in the transition region (TR) above sunspot penumbra. The presence of penumbral micro-jets (PMJs) in the chromosphere have also been reported earlier. It was proposed that both the PMJs and BDs are formed due to magnetic reconnection process and may play an important role in heating of the penumbra. Using simultaneous observation of the chromosphere from the Solar Optical Telescope (SOT) aboard Hinode and the TR from the Interface Region Imaging Spectrograph (IRIS), we study the dynamics of BDs and their relation with PMJs. We find two types of BDs, one which is related to PMJs and the others which do not show any visible dynamics in the SOT Ca II H images. From a statistical analysis we show that these two types have different properties. The BDs which are related to PMJs always appear at the top of the PMJs, the vast majority of which show inward motion and originate before the generation of the PMJs. These results may indicate that the reconnection occurs at the lower coronal/TR height and initiates PMJs at the chromosphere. This formation mechanism is in contrast with the currently believed formation of PMJs by reconnection in the (upper) photosphere between differently inclined fields.


title: Transmission in graphene-topological insulator heterostructures

abstract: We investigate scattering of the topological surface state of a three-dimensional time-reversal invariant topological insulator when graphene is deposited on the topological-insulator surface. Specifically, we consider the (111) surface of a Bi$_2$Se$_3$-like topological insulator. We present a low-energy model for the bulk graphene-topological insulator heterostructure and we calculate the transmission probability at zigzag and armchair edges of the deposited graphene, and the conductance through graphene nanoribbon barriers and show that its features can be understood from antiresonances in the transmission probability.


title: Visualisation of Multi-mission Astronomical Data with ESASky

abstract: ESASky is a science-driven discovery portal to explore the multi-wavelength sky and visualise and access multiple astronomical archive holdings. The tool is a web application that requires no prior knowledge of any of the missions involved and gives users world-wide simplified access to the highest-level science data products from multiple astronomical space-based astronomy missions plus a number of ESA source catalogues. The first public release of ESASky features interfaces for the visualisation of the sky in multiple wavelengths, the visualisation of query results summaries, and the visualisation of observations and catalogue sources for single and multiple targets. This paper describes these features within ESASky, developed to address use cases from the scientific community. The decisions regarding the visualisation of large amounts of data and the technologies used were made in order to maximise the responsiveness of the application and to keep the tool as useful and intuitive as possible.


title: Visible spectra of (474640) 2004 VN112-2013 RF98 with OSIRIS at the 10.4 m GTC: evidence for binary dissociation near aphelion among the extreme trans-Neptunian objects

abstract: The existence of significant anisotropies in the distributions of the directions of perihelia and orbital poles of the known extreme trans-Neptunian objects (ETNOs) has been used to claim that trans-Plutonian planets may exist. Among the known ETNOs, the pair (474640) 2004 VN112-2013 RF98 stands out. Their orbital poles and the directions of their perihelia and their velocities at perihelion/aphelion are separated by a few degrees, but orbital similarity does not necessarily imply common physical origin. In an attempt to unravel their physical nature, visible spectroscopy of both targets was obtained using the OSIRIS camera-spectrograph at the 10.4 m Gran Telescopio Canarias (GTC). From the spectral analysis, we find that 474640-2013 RF98 have similar spectral slopes (12 vs. 15 %/0.1um), very different from Sedna's but compatible with those of (148209) 2000 CR105 and 2012 VP113. These five ETNOs belong to the group of seven linked to the Planet Nine hypothesis. A dynamical pathway consistent with these findings is dissociation of a binary asteroid during a close encounter with a planet and we confirm its plausibility using N-body simulations. We thus conclude that both the dynamical and spectroscopic properties of 474640-2013 RF98 favour a genetic link and their current orbits suggest that the pair was kicked by a perturber near aphelion.


title: Transient growth of perturbations on scales beyond the accretion disc thickness

abstract: Turbulent state of spectrally stable shear flows may be developed and sustained according to the bypass scenario of transition. If it works in non-magnetised boundless and homogeneous quasi-Keplerian flow, transiently growing shearing vortices should supply turbulence with energy. Employing the large shearing box approximation, as well as a set of global disc models, we study the optimal growth of the shearing vortices in such a flow in the whole range of azimuthal length-scales, $\lambda_y$, as compared to the flow scale-height, $H$. It is shown that with the account of the viscosity the highest possible amplification of shearing vortices, $G_{max}$, attains maximum at $\lambda_y\lesssim H$ and declines towards both the large scales $\lambda_y\gg H$ and the small scales $\lambda_y\ll H$ in a good agreement with analytical estimations based on balanced solutions. We pay main attention to the large-scale vortices $\lambda_y\gg H$, which produce $G_{max}\propto (\Omega/\kappa)^4$, where $\Omega$ and $\kappa$ denote local rotational and epicyclic frequencies, respectively. It is demonstrated that the large-scale vortices acquire high density perturbation as they approach the instant of swing. At the same time, their growth is not affected by bulk viscosity. We check that $G_{max}$ obtained globally is comparable to its local counterpart and the shape and localisation of global optimal vortices can be explained in terms of the local approach. The obtained results allow us to suggest that the critical Reynolds number of subcritical transition to turbulence in quasi-Keplerian flow, as well as the corresponding turbulent effective azimuthal stress should substantially depend on shear rate.

categories: astro-ph.HE astro-ph.EP astro-ph.SR physics.flu-dyn


title: Odin observations of ammonia in the Sgr A +50 km/s Cloud and Circumnuclear Disk

abstract: Context. The Odin satellite is now into its sixteenth year of operation, much surpassing its design life of two years. One of the sources which Odin has observed in great detail is the Sgr A Complex in the centre of the Milky Way. Aims. To study the presence of NH3 in the Galactic Centre and spiral arms. Methods. Recently, Odin has made complementary observations of the 572 GHz NH3 line towards the Sgr A +50 km/s Cloud and Circumnuclear Disk (CND). Results. Significant NH3 emission has been observed in both the +50 km/s Cloud and the CND. Clear NH3 absorption has also been detected in many of the spiral arm features along the line of sight from the Sun to the core of our Galaxy. Conclusions. The very large velocity width (80 km/s) of the NH3 emission associated with the shock region in the southwestern part of the CND may suggest a formation/desorption scenario similar to that of gas-phase H2O in shocks/outflows.


title: V773 Cas, QS Aql, and BR Ind: Eclipsing Binaries as Parts of Multiple Systems

abstract: Eclipsing binaries remain crucial objects for our understanding of the universe. In particular, those that are components of multiple systems can help us solve the problem of the formation of these systems. Analysis of the radial velocities together with the light curve produced for the first time precise physical parameters of the components of the multiple systems V773 Cas, QS Aql, and BR Ind. Their visual orbits were also analyzed, which resulted in slightly improved orbital elements. What is typical for all these systems is that their most dominant source is the third distant component. The system V773 Cas consists of two similar G1-2V stars revolving in a circular orbit and a more distant component of the A3V type. Additionally, the improved value of parallax was calculated to be 17.6 mas. Analysis of QS Aql resulted in the following: the inner eclipsing pair is composed of B6V and F1V stars, and the third component is of about the B6 spectral type. The outer orbit has high eccentricity of about 0.95, and observations near its upcoming periastron passage between the years 2038 and 2040 are of high importance. Also, the parallax of the system was derived to be about 2.89 mas, moving the star much closer to the Sun than originally assumed. The system BR Ind was found to be a quadruple star consisting of two eclipsing K dwarfs orbiting each other with a period of 1.786 days; the distant component is a single-lined spectroscopic binary with an orbital period of about 6 days. Both pairs are moving around each other on their 148 year orbit.


title: Vandermonde Matrices with Nodes in the Unit Disk and the Large Sieve

abstract: We derive bounds on the extremal singular values and the condition number of NxK, with N>=K, Vandermonde matrices with nodes in the unit disk. The mathematical techniques we develop to prove our main results are inspired by a link---first established by by Selberg [1] and later extended by Moitra [2]---between the extremal singular values of Vandermonde matrices with nodes on the unit circle and large sieve inequalities. Our main conceptual contribution lies in establishing a connection between the extremal singular values of Vandermonde matrices with nodes in the unit disk and a novel large sieve inequality involving polynomials in z \in C with |z|<=1. Compared to Baz\'an's upper bound on the condition number [3], which, to the best of our knowledge, constitutes the only analytical result---available in the literature---on the condition number of Vandermonde matrices with nodes in the unit disk, our bound not only takes a much simpler form, but is also sharper for certain node configurations. Moreover, the bound we obtain can be evaluated consistently in a numerically stable fashion, whereas the evaluation of Baz\'an's bound requires the solution of a linear system of equations which has the same condition number as the Vandermonde matrix under consideration and can therefore lead to numerical instability in practice. As a byproduct, our result---when particularized to the case of nodes on the unit circle---slightly improves upon the Selberg-Moitra bound.


title: A new method to suppress the bias in polarized intensity

abstract: Computing polarised intensities from noisy data in Stokes U and Q suffers from a positive bias that should be suppressed. To develop a correction method that, when applied to maps, should provide a distribution of polarised intensity that closely follows the signal from the source. We propose a new method to suppress the bias by estimating the polarisation angle of the source signal in a noisy environment with help of a modified median filter. We then determine the polarised intensity, including the noise, by projection of the observed values of Stokes U and Q onto the direction of this polarisation angle. We show that our new method represents the true signal very well. If the noise distribution in the maps of U and Q is Gaussian, then in the corrected map of polarised intensity it is also Gaussian. Smoothing to larger Gaussian beamsizes, to improve the signal-to-noise ratio, can be done directly with our method in the map of the polarised intensity. Our method also works in case of non-Gaussian noise distributions. The maps of the corrected polarised intensities and polarisation angles are reliable even in regions with weak signals and provide integrated flux densities and degrees of polarisation without the cumulative effect of the bias, which especially affects faint sources. Features at low intensity levels like 'depolarisation canals' are smoother than in the maps using the previous methods, which has broader implications, for example on the interpretation of interstellar turbulence.


title: Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

abstract: Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a two-dimensional (2D) pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are a necessity to stabilize the MCMC inversion results of a highly-discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.


title: A Schur decomposition reveals the richness of structure in homogeneous, isotropic turbulence as a consequence of localised shear

abstract: An improved understanding of turbulence is essential for the effective modelling and control of industrial and geophysical processes. Homogeneous, isotropic turbulence (HIT) is the archetypal field for developing turbulence physics theory. Based on the Schur transform, we introduce an additive decomposition of the velocity gradient tensor into a normal part (containing the eigenvalues) and a non-normal or shear-related tensor. We re-interrogate some key properties of HIT and show that the the tendency of the flow to form disc-like structures is not a property of the normal tensor; it emerges from an interaction with the non-normality. Also, the alignment between the vorticity vector and the second eigenvector of the strain tensor is another consequence of local shear processes.


title: Analysis of eight binaries in Lyncis constellation: RV Lyn, AA Lyn, AH Lyn, CD Lyn, CF Lyn, DR Lyn, EK Lyn, and FS Lyn

abstract: The available photometry mainly from the WASP database was used for the first light curve analysis of eight eclipsing binary systems located in the Lyncis constellation: RV Lyn, AA Lyn, AH Lyn, CD Lyn, CF Lyn, DR Lyn, EK Lyn, and FS Lyn. Most of these eclipsing stars are detached ones, having the orbital periods from 0.54 to 2.3 days. For the systems AA Lyn and CF Lyn a non-negligible third light was detected during the light curve solution. Moreover, 284 new times of minima for these binaries were derived, trying to identify the period variations. For the system CD Lyn a hypothetical third body was detected with the period of about 59 yr.


title: Predicting Citywide Crowd Flows Using Deep Spatio-Temporal Residual Networks

abstract: Forecasting the flow of crowds is of great importance to traffic management and public safety, and very challenging as it is affected by many complex factors, including spatial dependencies (nearby and distant), temporal dependencies (closeness, period, trend), and external conditions (e.g., weather and events). We propose a deep-learning-based approach, called ST-ResNet, to collectively forecast two types of crowd flows (i.e. inflow and outflow) in each and every region of a city. We design an end-to-end structure of ST-ResNet based on unique properties of spatio-temporal data. More specifically, we employ the residual neural network framework to model the temporal closeness, period, and trend properties of crowd traffic. For each property, we design a branch of residual convolutional units, each of which models the spatial properties of crowd traffic. ST-ResNet learns to dynamically aggregate the output of the three residual neural networks based on data, assigning different weights to different branches and regions. The aggregation is further combined with external factors, such as weather and day of the week, to predict the final traffic of crowds in each and every region. We have developed a real-time system based on Microsoft Azure Cloud, called UrbanFlow, providing the crowd flow monitoring and forecasting in Guiyang City of China. In addition, we present an extensive experimental evaluation using two types of crowd flows in Beijing and New York City (NYC), where ST-ResNet outperforms nine well-known baselines.


title: Inflationary Primordial Black Holes as All Dark Matter

abstract: Following a new microlensing constraint on primordial black holes (PBHs) with $\sim10^{20}$--$10^{28}\,\mathrm{g}$~[1], we revisit the idea of PBH as all Dark Matter (DM). We have shown that the updated observational constraints suggest the viable mass function for PBHs as all DM to have a peak at $\simeq 10^{20}\,\mathrm{g}$ with a small width $\sigma \lesssim 0.1$, by imposing observational constraints on an extended mass function in a proper way. We have also provided an inflation model that successfully generates PBHs as all DM fulfilling this requirement.


title: IoFClime: The fuzzy logic and the Internet of Things to control indoor temperature regarding the outdoor ambient conditions

abstract: The Internet of Things is arriving to our homes or cities through fields already known like Smart Homes, Smart Cities, or Smart Towns. The monitoring of environmental conditions of cities can help to adapt the indoor locations of the cities in order to be more comfortable for people who stay there. A way to improve the indoor conditions is an efficient temperature control, however, it depends on many factors like the different combinations of outdoor temperature and humidity. Therefore, adjusting the indoor temperature is not setting a value according to other value. There are many more factors to take into consideration, hence the traditional logic based in binary states cannot be used. Many problems cannot be solved with a set of binary solutions and we need a new way of development. Fuzzy logic is able to interpret many states, more than two states, giving to computers the capacity to react in a similar way to people. In this paper we will propose a new approach to control the temperature using the Internet of Things together its platforms and fuzzy logic regarding not only the indoor temperature but also the outdoor temperature and humidity in order to save energy and to set a more comfortable environment for their users. Finally, we will conclude that the fuzzy approach allows us to achieve an energy saving around 40% and thus, save money.


title: MDPs with Energy-Parity Objectives

abstract: Energy-parity objectives combine $\omega$-regular with quantitative objectives of reward MDPs. The controller needs to avoid to run out of energy while satisfying a parity objective. We refute the common belief that, if an energy-parity objective holds almost-surely, then this can be realised by some finite memory strategy. We provide a surprisingly simple counterexample that only uses coB\"uchi conditions. We introduce the new class of bounded (energy) storage objectives that, when combined with parity objectives, preserve the finite memory property. Based on these, we show that almost-sure and limit-sure energy-parity objectives, as well as almost-sure and limit-sure storage parity objectives, are in $\mathit{NP}\cap \mathit{coNP}$ and can be solved in pseudo-polynomial time for energy-parity MDPs.


title: A Convenient Category for Higher-Order Probability Theory

abstract: Higher-order probabilistic programming languages allow programmers to write sophisticated models in machine learning and statistics in a succinct and structured way, but step outside the standard measure-theoretic formalization of probability theory. Programs may use both higher-order functions and continuous distributions, or even define a probability distribution on functions. But standard probability theory does not handle higher-order functions well: the category of measurable spaces is not cartesian closed. Here we introduce quasi-Borel spaces. We show that these spaces: form a new formalization of probability theory replacing measurable spaces; form a cartesian closed category and so support higher-order functions; form a well-pointed category and so support good proof principles for equational reasoning; and support continuous probability distributions. We demonstrate the use of quasi-Borel spaces for higher-order functions and probability by: showing that a well-known construction of probability theory involving random functions gains a cleaner expression; and generalizing de Finetti's theorem, that is a crucial theorem in probability theory, to quasi-Borel spaces.


title: A uniform stability principle for dual lattices

abstract: We prove a highly uniform stability or "almost-near" theorem for dual lattices of lattices $L \subseteq \Bbb R^n$. More precisely, we show that, for a vector $x$ from the linear span of a lattice $L \subseteq \Bbb R^n$, subject to $\lambda_1(L) \ge \lambda > 0$, to be $\varepsilon$-close to some vector from the dual lattice $L'$ of $L$, it is enough that the inner products $u\,x$ are $\delta$-close (with $\delta < 1/3$) to some integers for all vectors $u \in L$ satisfying $\| u \| \le r$, where $r > 0$ depends on $n$, $\lambda$, $\delta$ and $\varepsilon$, only. This generalizes an earlier analogous result proved for integral vector lattices by M. Ma\v{c}aj and the second author. The proof is nonconstructive, using the ultraproduct construction and a slight portion of nonstandard analysis.


title: Geometric Algebra and Information Geometry for Quantum Computational Software

abstract: The art of quantum algorithm design is highly nontrivial. Grover's search algorithm constitutes a masterpiece of quantum computational software. In this article, we use methods of geometric algebra (GA) and information geometry (IG) to enhance the algebraic efficiency and the geometrical significance of the digital and analog representations of Grover's algorithm, respectively. Specifically, GA is used to describe the Grover iterate and the discretized iterative procedure that exploits quantum interference to amplify the probability amplitude of the target-state before measuring the query register. The transition from digital to analog descriptions occurs via Stone's theorem which relates the (unitary) Grover iterate to a suitable (Hermitian) Hamiltonian that controls Schrodinger's quantum mechanical evolution of a quantum state towards the target state. Once the discrete-to-continuos transition is completed, IG is used to interpret Grover's iterative procedure as a geodesic path on the manifold of the parametric density operators of pure quantum states constructed from the continuous approximation of the parametric quantum output state in Grover's algorithm. Finally, we discuss the dissipationless nature of quantum computing, recover the quadratic speedup relation, and identify the superfluity of the Walsh-Hadamard operation from an IG perspective with emphasis on statistical mechanical considerations.


title: Thermal impact of magmatism in subduction zones

abstract: Magmatism in subduction zones builds continental crust and causes most of Earth's subaerial volcanism. The production rate and composition of magmas are controlled by the thermal structure of subduction zones. A range of geochemical and heat flow evidence has recently converged to indicate that subduction zones are hotter at lithospheric depths beneath the arc than predicted by canonical thermomechanical models, which neglect magmatism. We show that this discrepancy can be resolved by consideration of the heat transported by magma. In our one- and two-dimensional numerical models and scaling analysis, magmatic transport of sensible and latent heat locally alters the thermal structure of canonical models by $\sim$300 K, increasing predicted surface heat flow and mid-lithospheric temperatures to observed values. We find the advection of sensible heat to be larger than the deposition of latent heat. Based on these results we conclude that thermal transport by magma migration affects the chemistry and the location of arc volcanoes.


title: A character of Siegel modular group of level 2 from theta constants

abstract: Given a characteristic, we define a character of the Siegel modular group of level 2, the computations of their values are also obtained. By using our theorems, some key theorems of Igusa [1] can be recovered.


title: Is single-particle interference spooky?

abstract: It is said about quantum interference that "In reality, it contains the only mystery". Indeed, together with non-locality it is often considered as the characteristic feature of quantum theory which can not be explained in any classical way. In this work we are concerned with a restricted setting of a single particle propagating in multi-path interferometric circuits, that is physical realisation of a qudit. It is shown that this framework, including collapse of the wave function, can be simulated with classical resources without violating the locality principle. We present a local ontological model whose predictions are indistinguishable from the quantum case. 'Non-locality' in the model appears merely as an epistemic effect arising on the level of description by agents whose knowledge is incomplete. This result suggests that the real quantum mystery should be sought in the multi-particle behaviour, since single-particle interferometric phenomena are explicable in a classical manner.


title: Lorentz invariance violation in the neutrino sector: a joint analysis from big bang nucleosynthesis and the cosmic microwave background

abstract: We investigate constraints on Lorentz invariance violation in the neutrino sector from a joint analysis of big bang nucleosynthesis and the cosmic microwave background. The effect of Lorentz invariance violation during the epoch of big bang nucleosynthesis changes the predicted helium-4 abundance, which influences the power spectrum of the cosmic microwave background at the recombination epoch. In combination with the latest measurement of the primordial helium-4 abundance, the Planck 2015 data of the cosmic microwave background anisotropies give a strong constraint on the deformation parameter since adding the primordial helium measurement breaks the degeneracy between the deformation parameter and the physical dark matter density.


title: XMM-Newton observation of the nearby pulsar B1133+16

abstract: We constrain the X-ray properties of the nearby $(360\,{\rm pc})$, old ($5\,{\rm Myr}$) pulsar B1133+16 with $\sim 100\,{\rm ks}$ effective exposure time by {\it XMM-Newton}. The observed pulsar flux in the 0.2-3 keV energy range is $\sim 10^{-14} \, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$, which results in the recording of $\sim 600$ source counts with the EPIC pn and MOS detectors. The X-ray radiation is dominated by nonthermal radiation and is well described by both a single power-law model (PL) and a sum of blackbody and power-law emission (BB+PL). The BB+PL model results in a spectral photon index $\Gamma=2.4^{+0.4}_{-0.3}$ and a nonthermal flux in the 0.2-3 keV energy range of $(7\pm 2) \times 10^{-15}\, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$. The thermal emission is consistent with the blackbody emission from a small hot spot with a radius of $R_{\rm pc} \approx 14^{+7}_{-5} \, {\rm m}$ and a temperature of $T_{\rm s} = 2.9^{+0.6}_{-0.4} \, {\rm MK}$. Assuming that the hot spot corresponds to the polar cap of the pulsar, we can use the magnetic flux conservation law to estimate the magnetic field at the surface $B_{\rm s} \approx 3.9 \times 10^{14} \, {\rm G}$. The observations are in good agreement with the predictions of the partially screened gap model, which assumes the existence of small-scale surface magnetic field structures in the polar cap region.


title: The ANTS problem

abstract: We introduce the Ants Nearby Treasure Search (ANTS) problem, which models natural cooperative foraging behavior such as that performed by ants around their nest. In this problem, k probabilistic agents, initially placed at a central location, collectively search for a treasure on the two-dimensional grid. The treasure is placed at a target location by an adversary and the agents' goal is to find it as fast as possible as a function of both k and D, where D is the (unknown) distance between the central location and the target. We concentrate on the case in which agents cannot communicate while searching. It is straightforward to see that the time until at least one agent finds the target is at least $\Omega$(D + D 2 /k), even for very sophisticated agents, with unrestricted memory. Our algorithmic analysis aims at establishing connections between the time complexity and the initial knowledge held by agents (e.g., regarding their total number k), as they commence the search. We provide a range of both upper and lower bounds for the initial knowledge required for obtaining fast running time. For example, we prove that log log k + $\Theta$(1) bits of initial information are both necessary and sufficient to obtain asymptotically optimal running time, i.e., O(D +D 2 /k). We also we prove that for every 0 \textless{} \textless{} 1, running in time O(log 1-- k $\times$(D +D 2 /k)) requires that agents have the capacity for storing $\Omega$(log k) different states as they leave the nest to start the search. To the best of our knowledge, the lower bounds presented in this paper provide the first non-trivial lower bounds on the memory complexity of probabilistic agents in the context of search problems. We view this paper as a "proof of concept" for a new type of interdisciplinary methodology. To fully demonstrate this methodology, the theoretical tradeoff presented here (or a similar one) should be combined with measurements of the time performance of searching ants.


title: Radiation Memory, Boosted Schwarzschild Spacetimes and Supertranslations

abstract: We investigate gravitational radiation memory and its corresponding effect on the asymptotic symmetries of a body whose exterior is a boosted Schwarzschild spacetime. First, in the context of linearized theory, we consider such a Schwarzschild body which is initially at rest, then goes through a radiative stage and finally emerges as a boosted Schwarzschild body. We show that the proper retarded solution of the exterior Schwarzschild spacetime for this process can be described in terms of the ingoing Kerr-Schild form of the Schwarzschild metric for both the initial and final states. An outgoing Kerr-Schild or time symmetric metric does not give the proper solution. The special property of Kerr-Schild metrics that their linearized and nonlinear forms are identical allows us to extend this result to processes in the nonlinear regime. We then discuss how the nonlinear memory effect, and its associated supertranslation, affect angular momentum conservation. Our approach provides a new framework for studying nonlinear aspects of the memory effect.


title: Wiedemann-Franz law demonstration in a student practicum

abstract: Wiedemann-Franz law is a prediction of electronic theory of electric and thermal conductivity in metals, which states that a Lorenz ratio $L=\kappa/(\sigma T)$, where $\kappa$ is a thermal conductivity, $\sigma$ --- electric conductivity and $T$ --- absolute temperature, is a universal constant in certain cases. We present here a simple experimental setup to verify this prediction in a teaching experiment.


title: Majorana Zero Modes Protected by Hopf Invariant in Topologically Trivial Superconductors

abstract: Majorana zero modes are usually attributed to topological superconductors. We study a class of two-dimensional topologically trivial superconductors without chiral edge modes, which nevertheless host robust Majorana zero modes in topological defects. The construction of this minimal single-band model is facilitated by the Hopf map and the Hopf invariant. This work will stimulate investigations of Majorana zero modes in superconductors in the topologically trivial regime.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el quant-ph


title: T2K near detector constraints for oscillation results

abstract: The T2K experiment is a long-baseline accelerator neutrino experiment using a near detector complex ND280 and a far detector, Super-Kamiokande. Neutrino interactions are detected by Cherenkov light in Super-Kamiokande in order to measure neutrino oscillation parameters. The accuracy of the oscillation parameter measurements depends on our knowledge of neutrino interactions, the neutrino flux and the detector response. ND280 is composed of a tracker and several sub-detectors designed to characterise the neutrino beam before oscillation, which allow us to constrain the uncertainties on the neutrino interaction and the accelerator flux models. In this talk we present the result of a fit on the ND280 data, as well as a study of the robustness of the fit to the choice of neutrino interaction model. The constraint obtained on the oscillation parameters is also discussed.


title: Time Complexity Analysis of a Distributed Stochastic Optimization in a Non-Stationary Environment

abstract: In this paper, we consider a distributed stochastic optimization problem where the goal is to minimize the time average of a cost function subject to a set of constraints on the time averages of related stochastic processes called penalties. We assume that the state of the system is evolving in an independent and non-stationary fashion and the "common information" available at each node is distributed and delayed. Such stochastic optimization is an integral part of many important problems in wireless networks such as scheduling, routing, resource allocation and crowd sensing. We propose an approximate distributed Drift- Plus-Penalty (DPP) algorithm, and show that it achieves a time average cost (and penalties) that is within epsilon > 0 of the optimal cost (and constraints) with high probability. Also, we provide a condition on the convergence time t for this result to hold. In particular, for any delay D >= 0 in the common information, we use a coupling argument to prove that the proposed algorithm converges almost surely to the optimal solution. We use an application from wireless sensor network to corroborate our theoretical findings through simulation results.


title: Effect of strong magnetic field on competing order parameters in two-flavor dense quark matter

abstract: We study the effect of strong magnetic field on competing chiral and diquark order parameters in a regime of moderately dense quark matter. The inter-dependence of the chiral and diquark condensates through nonperturbative quark mass and strong coupling effects is analyzed in a two-flavor Nambu-Jona-Lasinio (NJL) model. In the weak magnetic field limit, our results agree qualitatively with earlier zero-field studies in the literature that find a critical coupling ratio $G_D/G_S\sim 1.1$ below which chiral or superconducting order parameters appear almost exclusively. Above the critical ratio, there exists a significant mixed broken phase region where both gaps are non-zero. However, a strong magnetic field $B\gtrsim 10^{18}$ G disrupts this mixed broken phase region and changes a smooth crossover found in the weak-field case to a first-order transition for both gaps at almost the same critical density. Our results suggest that in the two-flavor approximation to moderately dense quark matter, strong magnetic field enhances the possibility of a mixed phase at high density, with implications for the structure, energetics and vibrational spectrum of neutron stars.


title: Dynamical compensation in biological systems as a particular case of structural non-identifiability

abstract: Dynamical compensation (DC) has been recently defined as the ability of a biological system to keep its output dynamics unchanged in the face of varying parameters. This concept is purported to describe a design principle that provides robustness to physiological circuits. Here we note the similitude between DC and Structural Identifiability (SI), and we argue that the former can be explained in terms of (lack of) the latter. We propose to exploit this fact by using currently existing tools for SI analysis to perform DC analysis. We demonstrate the feasibility of this approach with four physiological circuits, for which we confirm the correspondence between DC and lack of SI. We also warn that care should we taken when using an unidentifiable model to extract biological insight, since lack of SI can be the result of an inappropriate choice of model structure and therefore not necessarily a sign of biological robustness.


title: A stochastic Pontryagin maximum principle on the Sierpinski gasket

abstract: In this paper, we consider stochastic control problems on the Sierpinski gasket. An order comparison lemma is derived using heat kernel estimate for Brownian motion on the gasket. Using the order comparison lemma and techniques of BSDEs, we establish a Pontryagin stochastic maximum principle for these control problems. It turns out that the stochastic maximum principle on the Sierpinski gasket involves two necessity equations in contrast to its counterpart on Euclidean spaces. This effect is due to singularity between the Hausdorff measure and the energy dominant measure on the gasket, which is a common feature shared by many fractal spaces. The linear regulator problems on the gasket is also considered as an example.


title: Free multivariate w*-semicrossed products: reflexivity and the bicommutant property

abstract: We study w*-semicrossed products over actions of the free semigroup and the free abelian semigroup on (possibly non-selfadjoint) w*-closed algebras. We show that they are reflexive when the dynamics are implemented by uniformly bounded families of invertible row operators. Combining with results of Helmer we derive that w*-semicrossed products over factors (on a separable Hilbert space) are reflexive. Furthermore we show that w*-semicrossed products of automorphic actions on maximal abelian selfadjoint algebras are reflexive. In all cases we prove that the w*-semicrossed products have the bicommutant property if and only if so does the ambient algebra of the dynamics.


title: The Energetic Cost of Building a Virus

abstract: Viruses are incapable of autonomous energy production. Although many experimental studies make it clear that viruses are parasitic entities that hijack the host's molecular resources, a detailed estimate for the energetic cost of viral synthesis is largely lacking. To quantify the energetic cost of viruses to their hosts, we enumerated the costs associated with two very distinct but representative DNA and RNA viruses, namely T4 and influenza. We found that for these viruses, translation of viral proteins is the most energetically expensive process. Interestingly, the cost of building a T4 phage and a single influenza virus are nearly the same. Due to influenza's higher burst size, however, the overall cost of a T4 phage infection is only 2-3% of the cost of an influenza infection. The costs of these infections relative to their host's estimated energy budget during the infection reveal that a T4 infection consumes about a third of its host's energy budget, where as an influenza infection consumes only 1%. Building on our estimates for T4, we show how the energetic costs of double-stranded DNA viruses scale with virus size, revealing that the dominant cost of building a virus can switch from translation to genome replication above a critical virus size. Lastly, using our predictions for the energetic cost of viruses, we provide estimates for the strengths of selection and genetic drift acting on newly incorporated genetic elements in viral genomes, under conditions of energy limitation.


title: Development of high vorticity in incompressible 3D Euler equations: influence of initial conditions

abstract: The incompressible three-dimensional ideal flows develop very thin pancake-like regions of increasing vorticity. These regions evolve with the scaling $\omega_{\max}(t)\propto\ell(t)^{-2/3}$ between the vorticity maximum and pancake thickness, and provide the leading contribution to the energy spectrum, where the gradual formation of the Kolmogorov interval $E_{k}\propto k^{-5/3}$ is observed for some initial flows [Agafontsev et. al, Phys. Fluids 27, 085102 (2015)]. With the massive numerical simulations, in the present paper we study the influence of initial conditions on the processes of pancake formation and the Kolmogorov energy spectrum development.


title: On the effectiveness of the thermoelectric energy filtering mechanism in low-dimensional superlattices and nano-composites

abstract: Electron energy filtering has been suggested as a promising way to improve the power factor and enhance the ZT figure of merit of thermoelectric materials. In this work we explore the effect that reduced dimensionality has on the success of the energy-filtering mechanism for power factor enhancement. We use the quantum mechanical non-equilibrium Green's function (NEGF) method for electron transport including electron-phonon scattering to explore 1D and 2D superlattice/nanocomposite systems. We find that, given identical material parameters, 1D channels utilize energy filtering more effectively than 2D as they: i) allow one to achieve maximal power factor for smaller well sizes / smaller grains (which is needed to maximize phonon scattering), ii) take better advantage of a lower thermal conductivity in the barrier/boundary materials compared to the well/grain materials in both: enhancing the Seebeck coefficient; and in producing a system which is robust against detrimental random deviations from optimal barrier design. In certain cases we find that the relative advantage can be as high as a factor of 3. We determine that energy-filtering is most effective when the average energy of carrier flow varies the most in the wells and the barriers along the channel, an event which appears when the energy of the carrier flow in the host material is low and when the energy relaxation mean-free-path of carriers is short. Although the ultimate reason these aspects, which cause a 1D system to see greater relative improvement than a 2D, is the 1D system's van Hove singularity in the density-of-states, the insights obtained are general and inform energy-filtering design beyond dimensional considerations.


title: Saturation of resistivity and Kohler's rule in Ni-doped La$_{1.85}$Sr$_{0.15}$CuO$_{4}$ cuprate

abstract: We present the results of electrical transport measurements of La$_{1.85}$Sr$_{0.15}$Cu$_{1-y}$Ni$_{y}$O$_{4}$ thin single-crystal films at magnetic fields up to 9 T. Adding Ni impurity with strong Coulomb scattering potential to slightly underdoped cuprate makes the signs of resistivity saturation at $\rho_{sat}$ visible in the measurement temperature window up to 350 K. Employing the parallel-resistor formalism reveals that $\rho_{sat}$ is consistent with classical Ioffe-Regel-Mott limit and changes with carrier concentration $n$ as $\rho_{sat}\propto 1/\sqrt{n}$. Thermopower measurements show that Ni tends to localize mobile carriers, decreasing their effective concentration as $n\!\cong0.15\!-\!y$. The classical unmodified Kohler's rule is fulfilled for magnetoresistance in the nonsuperconducting part of the phase diagram when applied to the ideal branch in the parallel-resistor model.


title: Geometric tuning of self-propulsion for Janus catalytic particles

abstract: Catalytic swimmers have attracted much attention as alternatives to biological systems for examining collective microscopic dynamics and the response to physico-chemical signals. Yet, understanding and predicting even the most fundamental characteristics of their individual propulsion still raises important challenges. While chemical asymmetry is widely recognized as the cornerstone of catalytic propulsion, different experimental studies have reported that particles with identical chemical properties may propel in opposite directions. Here, we show that, beyond its chemical properties, the detailed shape of a catalytic swimmer plays an essential role in determining its direction of motion, demonstrating the compatibility of the classical theoretical framework with experimental observations.

categories: cond-mat.soft physics.bio-ph physics.chem-ph physics.flu-dyn


title: Asymptotic expansions of holonomy

abstract: Given a principal bundle with a connection, we look for an asymptotic expansion of the holonomy of a loop in terms of its length. This length is defined relative to some Riemannian or sub-Riemannian structure. We are able to give an asymptotic formula that is independent of choice of gauge.


title: Stack Semantics of Type Theory

abstract: We give a model of dependent type theory with one univalent universe and propositional truncation interpreting a type as a stack, generalising the groupoid model of type theory. As an application, we show that countable choice cannot be proved in dependent type theory with one univalent universe and propositional truncation.


title: Loose social organisation of AB strain zebrafish groups in a two-patch environment

abstract: We explore the collective behaviours of 7 group sizes: 1, 2, 3, 5, 7, 10 and 20 AB zebrafish (Danio rerio) in a constraint environment composed of two identical squared rooms connected by a corridor. This simple set-up is similar to a natural patchy environment. We track the positions and the identities of the fish and compute the metrics at the group and at the individual levels. First, we show that the size of the population affects the behaviour of each individual in a group, the cohesion of the groups, the preferential interactions and the transition dynamics between the two rooms. Second, during collective departures, we show that the rankings of exit correspond to the topological organisations of the fish prior to their collective departure with no leadership. This spatial organisation emerge in the group a few seconds before a collective departure. These results provide new evidences on the spatial organisation of the groups and the effect of the population size on individual and collective behaviours in a patchy environment.


title: Relative singular locus and Balmer spectrum of matrix factorizations

abstract: For a separated Noetherian scheme $X$ with an ample family of line bundles and a non-zero-divisor $W\in\Gamma(X,L)$ of a line bundle $L$ on $X$, we classify certain thick subcategories of the derived matrix factorization category ${\rm DMF}(X,L,W)$ of the Landau-Ginzburg model $(X,L,W)$. Furthermore, by using the classification result and the theory of Balmer's tensor triangular geometry, we show that the spectrum of the tensor triangulated category $({\rm DMF}(X,L,W), \otimes^{\frac{1}{2}})$ is homeomorphic to the relative singular locus ${\rm Sing}(X_0/X)$, introduced in this paper, of the zero scheme $X_0\subset X$ of $W$.


title: $G_2$-structures and quantization of non-geometric M-theory backgrounds

abstract: We describe the quantization of a four-dimensional locally non-geometric M-theory background dual to a twisted three-torus by deriving a phase space star product for deformation quantization of quasi-Poisson brackets related to the nonassociative algebra of octonions. The construction is based on a choice of $G_2$-structure which defines a nonassociative deformation of the addition law on the seven-dimensional vector space of Fourier momenta. We demonstrate explicitly that this star product reduces to that of the three-dimensional parabolic constant $R$-flux model in the contraction of M-theory to string theory, and use it to derive quantum phase space uncertainty relations as well as triproducts for the nonassociative geometry of the four-dimensional configuration space. By extending the $G_2$-structure to a $Spin(7)$-structure, we propose a 3-algebra structure on the full eight-dimensional M2-brane phase space which reduces to the quasi-Poisson algebra after imposing a particular gauge constraint, and whose deformation quantisation simultaneously encompasses both the phase space star products and the configuration space triproducts. We demonstrate how these structures naturally fit in with previous occurences of 3-algebras in M-theory.


title: Local Okounkov bodies and limits in prime characteristic

abstract: This article is concerned with the asymptotic behavior of certain sequences of ideals in rings of prime characteristic. These sequences, which we call $p$-families of ideals, are ubiquitous in prime characteristic commutative algebra (e.g., they occur naturally in the theories of tight closure, Hilbert-Kunz multiplicity, and $F$-signature). We associate to each $p$-family of ideals an object in Euclidean space that is analogous to the Newton-Okounkov body of a graded family of ideals, which we call a $p$-body. Generalizing the methods used to establish volume formulas for the Hilbert-Kunz multiplicity and $F$-signature of semigroup rings, we relate the volume of a $p$-body to a certain asymptotic invariant determined by the corresponding $p$-family of ideals. We apply these methods to obtain new existence results for limits in positive characteristic, an analogue of the Brunn-Minkowski theorem for Hilbert-Kunz multiplicity, and a uniformity result concerning the positivity of a $p$-family.


title: Singularity formation and global existence of classical solutions for one dimensional rotating shallow water system

abstract: We study classical solutions of one dimensional rotating shallow water system which plays an important role in geophysical fluid dynamics. The main results contain two contrasting aspects. First, when the solution crosses certain threshold, we prove finite-time singularity formation for the classical solutions by studying the weighted gradients of Riemann invariants and utilizing conservation of physical energy. In fact, the singularity formation will take place for a large class of ${C}^1$ initial data whose gradients and physical energy can be arbitrarily small and higher order derivatives should be large. Second, when the initial data have constant potential vorticity, global existence of small classical solutions is established via studying an equivalent form of a quasilinear Klein-Gordon equation satisfying certain null conditions. In this global existence result, the smallness condition is in terms of the higher order Sobolev norms of the initial data.


title: Formulation, Implementation and Validation of the Horizontal Coupling Method for 1D/2D Shallow Water Flow Models

abstract: One dimensional (1D) simulations of the flow and flooding of open channels are known to be inaccurate as the flow is multi-dimensional in nature, especially at the flooded regions. However, multi-dimensional simulations, even in two dimensions (2D), are computationally expensive, hence the problem of efficiently coupling 2D and 1D simulations for the flow and flooding of open channels has been the subject of much research and is investigated in this paper. We adopt a 1D model with coupling term for the channel flow and the 2D shallow water flow model for the floodplain. The 1D model with coupling term is derived by integrating the 3D Free Surface Euler equations but without imposing any restriction on the channel width variations. %leading to a new coupling term. Finite volume methods are formulated for both the 2D and 1D models including a discrete coupling term in closed form. Coupling is achieved through the discrete coupling term in the 1D model and the lateral numerical fluxes in the 2D model. Since the lateral discharge in the channel cannot be guaranteed to be zero during flooding, we aim to recover the lateral variation by computing two lateral discharges over each cross section and propose to use an ad-hoc model based on the $y$-discharge equation in the 2D model for this purpose. We then propose the numerical scheme for this ad-hoc model following the hydrostatic reconstruction philosophy. Then, we show that the resulting method, named Horizontal Coupling Method (HCM), is well-balanced; we introduce the no-numerical flooding property and also show that the method satisfies the property. Three numerical test cases are used to verify the performance of the method. The results show that the method performs well in both accuracy and efficiency and also approximates the channel lateral discharges with very good accuracy and little computational overhead.


title: Multiprocessor Approximate Message Passing with Column-Wise Partitioning

abstract: Solving a large-scale regularized linear inverse problem using multiple processors is important in various real-world applications due to the limitations of individual processors and constraints on data sharing policies. This paper focuses on the setting where the matrix is partitioned column-wise. We extend the algorithmic framework and the theoretical analysis of approximate message passing (AMP), an iterative algorithm for solving linear inverse problems, whose asymptotic dynamics are characterized by state evolution (SE). In particular, we show that column-wise multiprocessor AMP (C-MP-AMP) obeys an SE under the same assumptions when the SE for AMP holds. The SE results imply that (i) the SE of C-MP-AMP converges to a state that is no worse than that of AMP and (ii) the asymptotic dynamics of C-MP-AMP and AMP can be identical. Moreover, for a setting that is not covered by SE, numerical results show that damping can improve the convergence performance of C-MP-AMP.


title: Difficulty of distinguishing product states locally

abstract: Non-locality without entanglement is a rather counter-intuitive phenomenon in which information may be encoded entirely in product (unentangled) states of composite quantum systems in such a way that local measurement of the subsystems is not enough for optimal decoding. For simple examples of pure product states, the gap in performance is known to be rather small when arbitrary local strategies are allowed. Here we restrict to local strategies readily achievable with current technology; those requiring neither a quantum memory nor joint operations. We show that, even for measurements on pure product states there can be a large gap between such strategies and theoretically optimal performance. Thus even in the absence of entanglement physically realizable local strategies can be far from optimal for extracting quantum information.


title: Local properties of the random Delaunay triangulation model and topological 2D gravity

abstract: Delaunay triangulations provide a bijection between a set of $N+3$ points in the complex plane, and the set of triangulations with given circumcircle intersection angles. The uniform Lebesgue measure on these angles translates into a K\"ahler measure for Delaunay triangulations, or equivalently on the moduli space $\mathcal M_{0,N+3}$ of genus zero Riemann surfaces with $N+3$ marked points. We study the properties of this measure. First we relate it to the topological Weil-Petersson symplectic form on the moduli space $\mathcal M_{0,N+3}$. Then we show that this measure, properly extended to the space of all triangulations on the plane, has maximality properties for Delaunay triangulations. Finally we show, using new local inequalities on the measures, that the volume $\mathcal{V}_N$ on triangulations with $N+3$ points is monotonically increasing when a point is added, $N\to N+1$. We expect that this can be a step towards seeing that the large $N$ limit of random triangulations can tend to the Liouville conformal field theory.


title: Lattice QCD with $N_f = 2+1+1 $ domain-wall quarks

abstract: We perform hybrid Monte Carlo simulation of (2+1+1)-flavors lattice QCD with the optimal domain-wall fermion (which has the effective 4D Dirac operator exactly equal to the Zolotarev optimal rational approximation of the overlap Dirac operator). The gauge ensemble is generated on the $32^3 \times 64 $ lattice with the extent $ N_s = 16 $ in the fifth dimension, and with the plaquette gauge action at $ \beta = 6/g^2 = 6.20 $. The lattice spacing ($ a \simeq 0.063 $ fm) is determined by the Wilson flow, using the value $ \sqrt{t_0} = 0.1416(8) $ fm obtained by the MILC Collaboration for the $(2+1+1)$-flavors QCD. The masses of $s$ and $c$ quarks are fixed by the masses of the vector mesons $ \phi(1020) $ and $ J/\psi(3097) $ respectively; while the mass of the $u/d$ quarks is heavier than their physical values, with the unitary pion mass $ M_\pi \simeq 280$ MeV (and $ M_\pi L \simeq 3 $). We compute the point-to-point quark propagators, and measure the time-correlation functions of meson and baryon interpolators. Our results of the mass spectra of the lowest-lying hadrons containing $ s $ and $ c $ quarks are in good agreement with the high energy experimental values, together with the predictions of the charmed baryons which have not been observed in experiments.