title: On Distributed Linear Estimation With Observation Model Uncertainties

abstract: We consider distributed estimation of a Gaussian source in a heterogenous bandwidth constrained sensor network, where the source is corrupted by independent multiplicative and additive observation noises, with incomplete statistical knowledge of the multiplicative noise. For multi-bit quantizers, we derive the closed-form mean-square-error (MSE) expression for the linear minimum MSE (LMMSE) estimator at the FC. For both error-free and erroneous communication channels, we propose several rate allocation methods named as longest root to leaf path, greedy and integer relaxation to (i) minimize the MSE given a network bandwidth constraint, and (ii) minimize the required network bandwidth given a target MSE. We also derive the Bayesian Cramer-Rao lower bound (CRLB) and compare the MSE performance of our proposed methods against the CRLB. Our results corroborate that, for low power multiplicative observation noises and adequate network bandwidth, the gaps between the MSE of our proposed methods and the CRLB are negligible, while the performance of other methods like individual rate allocation and uniform is not satisfactory.


title: Irrational points on random hyperelliptic curves

abstract: We consider genus $g$ hyperelliptic curves over $\mathbb{Q}$ with a rational Weierstrass point, ordered by height. If $d < g$ is odd, we prove, under an assumption, that there exists $B_d$ such that a positive proportion of these curves have at most $B_d$ points of degree $d$. If $d < g$ is even, we conditionally bound degree $d$ points not pulled back from points of degree $d/2$ on the projective line. Furthermore, we show one may take $B_2=24$ and B_3=114$. Our proofs proceed by refining recent work of Park, which applied tropical geometry to symmetric power Chabauty, and then applying results of Bhargava and Gross on average ranks of Jacobians of hyperelliptic curves.


title: Kinetics of the simplest Criegee intermediate reaction with ozone studied by mid-infrared quantum cascade laser spectrometer

abstract: The kinetics of CH$_2$OO reaction with ozone has been studied by monitoring CH$_2$OO with time-resolved infrared (IR) absorption spectroscopy, which utilized fast chirped IR pulse train from a quantum cascade laser [J. Chem. Phys. 146, 244302 (2017)]. CH$_2$OO was prepared by photolyzing a gas mixture of CH$_2$I$_2$/O$_2$/O$_3$ at 352 nm; the photolysis wavelength was chosen to minimize the photodissociation of O$_3$. The measured rate coefficient at 298 K and 30 Torr is (6.7$\pm$0.5)x10$^{-14}$ cm$^3$sec$^{-1}$, independent of pressure from 30 to 100 Torr. The result indicates that previous ab initio calculations either underestimated or overestimated this reaction rate by one order of magnitude or more. The result also implies that in laboratory studies of ozonolysis of alkenes, the reaction of Criegee intermediate with ozone may play a role. However, this reaction would not compete with other CH$_2$OO sinks in the atmosphere.


title: The Mating Rituals of Deep Neural Networks: Learning Compact Feature Representations through Sexual Evolutionary Synthesis

abstract: Evolutionary deep intelligence was recently proposed as a method for achieving highly efficient deep neural network architectures over successive generations. Drawing inspiration from nature, we propose the incorporation of sexual evolutionary synthesis. Rather than the current asexual synthesis of networks, we aim to produce more compact feature representations by synthesizing more diverse and generalizable offspring networks in subsequent generations via the combination of two parent networks. Experimental results were obtained using the MNIST and CIFAR-10 datasets, and showed improved architectural efficiency and comparable testing accuracy relative to the baseline asexual evolutionary neural networks. In particular, the network synthesized via sexual evolutionary synthesis for MNIST had approximately double the architectural efficiency (cluster efficiency of 34.29X and synaptic efficiency of 258.37X) in comparison to the network synthesized via asexual evolutionary synthesis, with both networks achieving a testing accuracy of ~97%.


title: Transient behavior of the solutions to the second order difference equations by the renormalization method based on Newton-Maclaurin expansion

abstract: The renormalization method based on the Newton-Maclaurin expansion is applied to study the transient behavior of the solutions to the difference equations as they tend to the steady-states. The key and also natural step is to make the renormalization equations to be continuous such that the elementary functions can be used to describe the transient behavior of the solutions to difference equations. As the concrete examples, we deal with the important second order nonlinear difference equations with a small parameter. The result shows that the method is more natural than the multi-scale method.

categories: math-ph math.AP math.CA math.DS math.MP physics.comp-ph


title: Optimal integrability threshold for Gibbs measures associated with focusing NLS on the torus

abstract: We study an optimal mass threshold for normalizability of the Gibbs measures associated with the focusing mass-critical nonlinear Schr\"odinger equation on the one-dimensional torus. In an influential paper, Lebowitz, Rose, and Speer (1988) proposed a critical mass threshold given by the mass of the ground state on the real line. We provide a proof for the optimality of this critical mass threshold. The proof also applies to the two-dimensional radial problem posed on the unit disc. In this case, we answer a question posed by Bourgain and Bulut (2014) on the optimal mass threshold. Furthermore, in the one-dimensional case, we show that the Gibbs measure is indeed normalizable at the optimal mass threshold, thus answering an open question posed by Lebowitz, Rose, and Speer (1988). This normalizability at the optimal mass threshold is rather striking in view of the minimal mass blowup solution for the focusing quintic nonlinear Schr\"odinger equation on the one-dimensional torus.


title: Properties of Kinetic Transition Networks for Atomic Clusters and Glassy Solids

abstract: A database of minima and transition states corresponds to a network where the minima represent nodes and the transition states correspond to edges between the pairs of minima they connect via steepest-descent paths. Here we construct networks for small clusters bound by the Morse potential for a selection of physically relevant parameters, in two and three dimensions. The properties of these unweighted and undirected networks are analysed to examine two features: whether they are small-world, where the shortest path between nodes involves only a small number or edges; and whether they are scale-free, having a degree distribution that follows a power law. Small-world character is present, but statistical tests show that a power law is not a good fit, so the networks are not scale-free. These results for clusters are compared with the corresponding properties for the molecular and atomic structural glass formers ortho-terphenyl and binary Lennard-Jones. These glassy systems do not show small-world properties, suggesting that such behaviour is linked to the structure-seeking landscapes of the Morse clusters.

categories: cond-mat.dis-nn cond-mat.stat-mech physics.comp-ph stat.AP


title: Spectral theory of multiplication operators on Hardy-Sobolev spaces

abstract: For a pointwise multiplier $\varphi$ of the Hardy-Sobolev space $H^2_\beta$ on the open unit ball $\bn$ in $\cn$, we study spectral properties of the multiplication operator $M_\varphi: H^2_\beta\to H^2_\beta$. In particular, we compute the spectrum and essential spectrum of $M_\varphi$ and develop the Fredholm theory for these operators.


title: Finite energy solutions to inhomogeneous nonlinear elliptic equations with sub-natural growth terms

abstract: We obtain necessary and sufficient conditions for the existence of a positive finite energy solution to the inhomogeneous quasilinear elliptic equation \[ -\Delta_{p} u = \sigma u^{q} + \mu \quad \text{on} \;\; \mathbb{R}^n \] in the sub-natural growth case $0<q<p-1$, where $\Delta_{p}$ ($1<p<\infty$) is the $p$-Laplacian, and $\sigma$, $\mu$ are positive Borel measures on $\mathbb{R}^n$. Uniqueness of such a solution is established as well. Similar inhomogeneous problems in the sublinear case $0<q<1$ are treated for the fractional Laplace operator $(-\Delta)^{\alpha}$ in place of $-\Delta_{p}$, on $\mathbb{R}^n$ for $0<\alpha<\frac{n}{2}$, and on an arbitrary domain $\Omega \subset \mathbb{R}^n$ with positive Green's function in the classical case $\alpha = 1$.


title: 2D Bayesian automated tilted-ring fitting of disk galaxies in large HI galaxy surveys: 2DBAT

abstract: We present a novel algorithm based on a Bayesian method for 2D tilted-ring analysis of disk galaxy velocity fields. Compared to the conventional algorithms based on a chi-squared minimisation procedure, this new Bayesian-based algorithm suffers less from local minima of the model parameters even with highly multi-modal posterior distributions. Moreover, the Bayesian analysis, implemented via Markov Chain Monte Carlo (MCMC) sampling, only requires broad ranges of posterior distributions of the parameters, which makes the fitting procedure fully automated. This feature will be essential when performing kinematic analysis on the large number of resolved galaxies expected to be detected in neutral hydrogen (HI) surveys with the Square Kilometre Array (SKA) and its pathfinders. The so-called '2D Bayesian Automated Tilted-ring fitter' (2DBAT) implements Bayesian fits of 2D tilted-ring models in order to derive rotation curves of galaxies. We explore 2DBAT performance on (a) artificial HI data cubes built based on representative rotation curves of intermediate-mass and massive spiral galaxies, and (b) Australia Telescope Compact Array (ATCA) HI data from the Local Volume HI Survey (LVHIS). We find that 2DBAT works best for well-resolved galaxies with intermediate inclinations (20 deg < i < 70 deg), complementing three-dimensional techniques better suited to modelling inclined galaxies.


title: Geometry of Information Integration

abstract: Information geometry is used to quantify the amount of information integration within multiple terminals of a causal dynamical system. Integrated information quantifies how much information is lost when a system is split into parts and information transmission between the parts is removed. Multiple measures have been proposed as a measure of integrated information. Here, we analyze four of the previously proposed measures and elucidate their relations from a viewpoint of information geometry. Two of them use dually flat manifolds and the other two use curved manifolds to define a split model. We show that there are hierarchical structures among the measures. We provide explicit expressions of these measures.


title: Habitable Climate Scenarios for Proxima Centauri b With a Dynamic Ocean

abstract: The nearby exoplanet Proxima Centauri b will be a prime future target for characterization, despite questions about its retention of water. Climate models with static oceans suggest that an Earth-like Proxima b could harbor a small dayside region of surface liquid water at fairly warm temperatures despite its weak instellation. We present the first 3-dimensional climate simulations of Proxima b with a dynamic ocean. We find that an ocean-covered Proxima b could have a much broader area of surface liquid water but at much colder temperatures than previously suggested, due to ocean heat transport and depression of the freezing point by salinity. Elevated greenhouse gas concentrations do not necessarily produce more open ocean area because of possible dynamic regime transitions. For an evolutionary path leading to a highly saline present ocean, Proxima b could conceivably be an inhabited, mostly open ocean planet dominated by halophilic life. For an ocean planet in 3:2 spin-orbit resonance, a permanent tropical waterbelt exists for moderate eccentricity. Simulations of Proxima Centauri b may also be a model for the habitability of planets receiving similar instellation from slightly cooler or warmer stars, e.g., in the TRAPPIST-1, LHS 1140, GJ 273, and GJ 3293 systems.


title: How do stars gain their mass? A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions

abstract: Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.


title: Conditions for Coronal Observation at the Lijiang Observatory in 2011

abstract: The sky brightness is a critical parameter for estimating the coronal observation conditions for solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Tele- scope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths accounts for 40.41 % of the total observing time in a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.


title: Focusing Attention: Towards Accurate Text Recognition in Natural Images

abstract: Scene text recognition has been a hot research topic in computer vision due to its various applications. The state of the art is the attention-based encoder-decoder framework that learns the mapping between input images and output sequences in a purely data-driven way. However, we observe that existing attention-based methods perform poorly on complicated and/or low-quality images. One major reason is that existing methods cannot get accurate alignments between feature areas and targets for such images. We call this phenomenon "attention drift". To tackle this problem, in this paper we propose the FAN (the abbreviation of Focusing Attention Network) method that employs a focusing attention mechanism to automatically draw back the drifted attention. FAN consists of two major components: an attention network (AN) that is responsible for recognizing character targets as in the existing methods, and a focusing network (FN) that is responsible for adjusting attention by evaluating whether AN pays attention properly on the target areas in the images. Furthermore, different from the existing methods, we adopt a ResNet-based network to enrich deep representations of scene text images. Extensive experiments on various benchmarks, including the IIIT5k, SVT and ICDAR datasets, show that the FAN method substantially outperforms the existing methods.


title: Event-Triggered Stabilization of Nonlinear Systems with Time-Varying Sensing and Actuation Delay

abstract: This paper studies the problem of stabilization of a nonlinear system with time-varying delays in both sensing and actuation using event-triggered control. Our proposed strategy seeks to opportunistically minimize the number of control updates while guaranteeing stabilization and builds on predictor feedback to compensate for arbitrarily large known time-varying delays. We establish, using a Lyapunov approach, the global asymptotic stability of the closed-loop system as long as the open-loop system is globally input-to-state stabilizable in the absence of time delays and sampling. We further prove that the proposed event-triggered law has inter-event times that are uniformly lower bounded and hence does not exhibit Zeno behavior. For the particular case of a stabilizable linear system, we show global exponential stability of the closed-loop system and analyze the trade-off between the rate of exponential convergence and a bound on the sampling frequency. We illustrate these results in simulation and also examine the properties of the proposed event-triggered strategy beyond the class of systems for which stabilization can be guaranteed.


title: High cumulants of conserved charges and their statistical uncertainties

abstract: We study the influence of measured high cumulants of conserved charges on their associated statistical uncertainties in relativistic heavy-ion collisions. With a given number of events, the measured cumulants randomly fluctuate with an approximately normal distribution, while the estimated statistical uncertainties are found to be correlated with corresponding values of the obtained cumulants. Generally, with a given number of events, the larger the cumulants we measure, the larger the statistical uncertainties that are estimated. The error-weighted averaged cumulants are dependent on statistics. Despite this effect, however, it is found that the three sigma rule of thumb is still applicable when the statistics are above one million.


title: Topological phase transitions in the photonic spin Hall effect

abstract: The recent synthesis of two-dimensional staggered materials opens up burgeoning opportunities to study optical spin-orbit interactions in semiconducting Dirac-like systems. We unveil topological phase transitions in the photonic spin Hall effect in the graphene family materials. It is shown that an external static electric field and a high frequency circularly polarized laser allow for active on-demand manipulation of electromagnetic beam shifts. The spin Hall effect of light presents a rich dependence with radiation degrees of freedom, material properties, and features non-trivial topological properties. We discover that photonic Hall shifts are sensitive to spin and valley properties of the charge carries, providing a unprecedented pathway to investigate spintronics and valleytronics in staggered 2D semiconductors.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.other physics.optics


title: The Detection Rates of Merging Binary Black Holes Originating from Star Clusters and Their Mass Function

abstract: Advanced LIGO achieved the first detection of the gravitational wave, which was from a merging binary black hole (BBH). In the near future, more merger events will be observed, and the mass distribution of them will become available. The mass distribution of merger events reflects the evolutionary path of BBHs: dynamical formation in dense star clusters or common envelope evolution from primordial binaries. In this paper, we estimate the detection rate of merging BBHs which dynamically formed in dense star clusters by combining the results of $N$-body simulations, modeling of globular clusters, and cosmic star-cluster formation history. We estimate that the merger rate density in the local universe within the redshift of 0.1 is 13--57 Gpc$^{-3}$ yr$^{-1}$. We find that the detection rate is 0.23--4.6 per year for the current sensitivity limit and that it would increase to 5.1--99 per year for the designed sensitivity which will be achieved in 2019. The distribution of merger rate density in the local universe as a function of redshifted chirp mass has a peak close to the low-mass end.The chirp mass function of the detected mergers, on the other hand, has a peak at the high-mass end, but is almost flat. This difference is simply because the detection range is larger for more massive BBHs.


title: Underlying one-step methods and nonautonomous stability of general linear methods

abstract: We generalize the theory of underlying one-step methods to strictly stable general linear methods (GLMs) solving nonautonomous ordinary differential equations (ODEs) that satisfy a global Lipschitz condition. We combine this theory with the Lyapunov and Sacker-Sell spectral stability theory for one-step methods developed in [34,35,36] to analyze the stability of a strictly stable GLM solving a nonautonomous linear ODE. These results are applied to develop a stability diagnostic for the solution of nonautonomous linear ODEs by strictly stable GLMs.


title: A philosophical argument against time machines

abstract: Time machines are predictions of Einstein's theory of general relativity and provide a myriad of unsolved paradoxes. Convincing and general arguments against time machines and their paradoxes are missing in physics and philosophy so far. In this article, a philosophical argument against time machines is given. When thought of as a process, individuation refuses the idea of time machines, in particular travels into the past. With the aid of Nietzsche-Heraclitus' philosophy of becoming and Simondon's notion of process of individuation, I propose that time machines are modern fables, created by the man of $\textit{ressentiment}$. In the $\textit{amor fati}$ formula of Nietzsche, I suggest the antipode to time machines.


title: Vogan classes in type $B_n$

abstract: Consider a weighted Coxeter system $(W,S,\mathscr{L})$. Via its associated Iwahori-Hecke algebra, we may determine the partition of $W$ into Kazhdan-Lusztig cells. In this paper, we use the theory of Vogan classes introduced by Bonnaf\'e--Geck to obtain a combinatorial description of the left cells of type $B_n$ when the ratio of the weights of the first to second generator is $n-1$. We further give information on the left cells when this ratio lies in the interval $(n-2,n-1)$.


title: Lattice-based designs possessing quasi-optimal separation distance on all projections

abstract: Experimental designs that spread out points apart from each other on projections are important for computer experiments when not necessarily all factors have substantial influence on the response. We provide a theoretical framework to generate designs that possess quasi-optimal separation distance on all of the projections and quasi-optimal fill distance on univariate margins. The key is to use special techniques to rotate certain lattices. One such type of design is densest packing-based maximum projection designs, which outperform existing types of space-filling designs in many scenarios. Computer code to generate these designs is provided in R package LatticeDesign.


title: Ranking ideas for diversity and quality

abstract: When selecting ideas or trying to find inspiration, designers often must sift through hundreds or thousands of ideas. This paper provides an algorithm to rank design ideas such that the ranked list simultaneously maximizes the quality and diversity of recommended designs. To do so, we first define and compare two diversity measures using Determinantal Point Processes (DPP) and additive sub-modular functions. We show that DPPs are more suitable for items expressed as text and that a greedy algorithm diversifies rankings with both theoretical guarantees and empirical performance on what is otherwise an NP-Hard problem. To produce such rankings, this paper contributes a novel way to extend quality and diversity metrics from sets to permutations of ranked lists. These rank metrics open up the use of multi-objective optimization to describe trade-offs between diversity and quality in ranked lists. We use such trade-off fronts to help designers select rankings using indifference curves. However, we also show that rankings on trade-off front share a number of top-ranked items; this means reviewing items (for a given depth like the top 10) from across the entire diversity-to-quality front incurs only a marginal increase in the number of designs considered. While the proposed techniques are general purpose enough to be used across domains, we demonstrate concrete performance on selecting items in an online design community (OpenIDEO), where our approach reduces the time required to review diverse, high-quality ideas from around 25 hours to 90 minutes. This makes evaluation of crowd-generated ideas tractable for a single designer. Our code is publicly accessible for further research.


title: Mid-infrared Plasmonic Circular Dichroism Generated by Graphene Nanodisk Assemblies

abstract: It is very interesting to bring plasmonic circular dichroism spectroscopy to the mid-infrared spectral interval, and there are two reasons for this. This spectral interval is very important for thermal bio-imaging and, simultaneously, this spectral range includes vibrational lines of many chiral biomolecules. Here we demonstrate that graphene plasmons indeed offer such opportunity. In particular, we show that chiral graphene assemblies consisting of a few graphene nanodisks can generate strong circular dichroism (CD) in the mid-infrared interval. The CD signal is generated due to the plasmon-plasmon coupling between adjacent nanodisks in the specially designed chiral graphene assemblies. Because of the large dimension mismatch between the thickness of a graphene layer and the incoming light's wavelength, three-dimensional configurations with a total height of a few hundred nanometers are necessary to obtain a strong CD signal in the mid-infrared range. The mid-infrared CD strength is mainly governed by the total dimensions (total height and helix scaffold radius) of the graphene nanodisk assembly, and by the plasmon-plasmon interaction strength between its constitutive nanodisks. Both positive and negative CD bands can be observed in the graphene assembly array. The frequency interval of the plasmonic CD spectra overlaps with the vibrational modes of some important biomolecules, such as DNA and many different peptides, giving rise to the possibility of enhancing the vibrational optical activity of these molecular species by attaching them to the graphene assemblies. Simultaneously the spectral range of chiral mid-infrared plasmons in our structures appears near the typical wavelength of the human-body thermal radiation and, therefore, our chiral metastructures can be potentially utilized as optical components in thermal imaging devices.


title: Nil Clean Ideal

abstract: Motivated by the concept of clean ideals, we introduce the notion of nil clean ideals of a ring. We define an ideal $I$ of a ring $R$ to be nil clean ideal if every element of $I$ can be written as a sum of an idempotent and a nilpotent element of $R$. In this article we discuss various properties of nil clean ideals.


title: Formulation of Deep Reinforcement Learning Architecture Toward Autonomous Driving for On-Ramp Merge

abstract: Multiple automakers have in development or in production automated driving systems (ADS) that offer freeway-pilot functions. This type of ADS is typically limited to restricted-access freeways only, that is, the transition from manual to automated modes takes place only after the ramp merging process is completed manually. One major challenge to extend the automation to ramp merging is that the automated vehicle needs to incorporate and optimize long-term objectives (e.g. successful and smooth merge) when near-term actions must be safely executed. Moreover, the merging process involves interactions with other vehicles whose behaviors are sometimes hard to predict but may influence the merging vehicle optimal actions. To tackle such a complicated control problem, we propose to apply Deep Reinforcement Learning (DRL) techniques for finding an optimal driving policy by maximizing the long-term reward in an interactive environment. Specifically, we apply a Long Short-Term Memory (LSTM) architecture to model the interactive environment, from which an internal state containing historical driving information is conveyed to a Deep Q-Network (DQN). The DQN is used to approximate the Q-function, which takes the internal state as input and generates Q-values as output for action selection. With this DRL architecture, the historical impact of interactive environment on the long-term reward can be captured and taken into account for deciding the optimal control policy. The proposed architecture has the potential to be extended and applied to other autonomous driving scenarios such as driving through a complex intersection or changing lanes under varying traffic flow conditions.


title: Counting the number of Feynman Graphs in QCD

abstract: Information about the number of Feynman graphs for a given physical process in a given field theory is especially useful for confirming the result of a Feynman graph generator used in an automatic system of perturbative calculations. A method of counting the number of Feynman graphs with weight of symmetry factor was established based on zero-dimensional field theory, and was used in scalar theories and QED. In this article this method is generalized to more complicated models by direct calculation of generating functions on a symbolic calculating system. This method is applied to QCD with and without counter terms, where many higher order are being calculated automatically.


title: Observation of the isospin breaking decay $\Upsilon(10860)\to\Upsilon(1S)f_0(980)\to\Upsilon(1S)\eta\pi^0$ with the Belle II detector

abstract: Using the Belle data on the decay $\Upsilon(10860)\to \Upsilon(1S)f_0(980)\to\Upsilon(1S)\pi^+\pi^-$, we estimate the fraction of the process $\Upsilon(10860)\to\Upsilon(1S)f_0(980)\to \Upsilon(1S)\eta\pi^0$, caused by the mixing of $a^0_0(980)$ and $f_0(980)$ resonances that breaks the isotopic invariance due to the $K^+$ and $K^0$ meson mass difference. With an instantaneous luminosity of $8\times10^{35}$ cm$^2\cdot$\,s$^{-1}$ at the SuperKEKB collider, one can collect with the Belle II detector about a hundred $\Upsilon(10860)\to\Upsilon(1S)\eta\pi^0$ events in the narrow region of the $\eta\pi^0$ invariant mass near the $K\bar K$ thresholds in a time roughly spent on the present Belle experiment.


title: An Alternative Approach to Functional Linear Partial Quantile Regression

abstract: Functional data such as curves and surfaces have become more and more common with modern technological advancements. The use of functional predictors remains challenging due to its inherent infinite-dimensionality. The common practice is to project functional data into a finite dimensional space. The popular partial least square (PLS) method has been well studied for the functional linear model [1]. As an alternative, quantile regression provides a robust and more comprehensive picture of the conditional distribution of a response when it is non-normal, heavy-tailed, or contaminated by outliers. While partial quantile regression (PQR) was proposed in [2], no theoretical guarantees were provided due to the iterative nature of the algorithm and the non-smoothness of quantile loss function. To address these issues, we propose an alternative PQR (APQR) formulation with guaranteed convergence. This novel formulation motivates new theories and allows us to establish asymptotic properties. Numerical studies on a benchmark dataset show the superiority of our new approach. We also apply our novel method to a functional magnetic resonance imaging (fMRI) data to predict attention deficit hyperactivity disorder (ADHD) and a diffusion tensor imaging (DTI) dataset to predict Alzheimer's disease (AD).


title: Jastrow form of the Ground State Wave Functions for Fractional Quantum Hall States

abstract: The topological morphology--order of zeros at the positions of electrons with respect to a specific electron--of Laughlin state at filling fractions $1/m$ ($m$ odd) is homogeneous as every electron feels zeros of order $m$ at the positions of other electrons. Although fairly accurate ground state wave functions for most of the other quantum Hall states in the lowest Landau level are quite well-known, it had been an open problem in expressing the ground state wave functions in terms of flux-attachment to particles, {\em a la}, this morphology of Laughlin state. With a very general consideration of flux-particle relations only, in spherical geometry, we here report a novel method for determining morphologies of these states. Based on these, we construct almost exact ground state wave-functions for the Coulomb interaction. Although the form of interaction may change the ground state wave-function, the same morphology constructs the latter irrespective of the nature of the interaction between electrons.


title: Lozenge tilings of a halved hexagon with an array of triangles removed from the boundary, part II

abstract: Proctor's work on staircase plane partitions yields an enumeration of lozenge tilings of a halved hexagon on the triangular lattice. Rohatgi later extended this tiling enumeration to a halved hexagon with a triangle cut off from the boundary. In the previous paper, the author proved a common generalization of Proctor's and Rohatgi's results by enumerating lozenge tilings of a halved hexagon in the case an array of an arbitrary number of triangles has been removed from a non-staircase side. In this paper, we consider the other case when the array of triangles has been removed from the staircase side of the halved hexagon. Our result also implies an explicit formula for the number of tilings of a hexagon with an array of triangles missing on the symmetry axis.


title: Cosmic Magnetism in Centimeter and Meter Wavelength Radio Astronomy

abstract: Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography.


title: Deep Embedding Convolutional Neural Network for Synthesizing CT Image from T1-Weighted MR Image

abstract: Recently, more and more attention is drawn to the field of medical image synthesis across modalities. Among them, the synthesis of computed tomography (CT) image from T1-weighted magnetic resonance (MR) image is of great importance, although the mapping between them is highly complex due to large gaps of appearances of the two modalities. In this work, we aim to tackle this MR-to-CT synthesis by a novel deep embedding convolutional neural network (DECNN). Specifically, we generate the feature maps from MR images, and then transform these feature maps forward through convolutional layers in the network. We can further compute a tentative CT synthesis from the midway of the flow of feature maps, and then embed this tentative CT synthesis back to the feature maps. This embedding operation results in better feature maps, which are further transformed forward in DECNN. After repeat-ing this embedding procedure for several times in the network, we can eventually synthesize a final CT image in the end of the DECNN. We have validated our proposed method on both brain and prostate datasets, by also compar-ing with the state-of-the-art methods. Experimental results suggest that our DECNN (with repeated embedding op-erations) demonstrates its superior performances, in terms of both the perceptive quality of the synthesized CT image and the run-time cost for synthesizing a CT image.


title: Non-stationary Markovian Replication Process causing Diverse Diffusions

abstract: We introduce a single generative mechanism with which it is able to describe diverse non-stationary diffusions. A non-stationary Markovian replication process for steps is considered, for which we analytically derive time-evolution of the probability distribution of the walker's displacement and the generalized telegrapher equation with time-varying coefficients, and find that diffusivity can be determined by temporal changes of replication of a immediate step. By controlling the replications, we realize the diverse diffusions such as alternating diffusions, superdiffusions, subdiffusions, and marginal diffusions which are originated from oscillating, increasing,decreasing, and slowly increasing or decreasing replications with time, respectively.


title: Laughlin wave function, Berry Phase and Quantization

abstract: In this paper, we show that the Laughlin wave function is a Hamiltonian and its associated Berry connection as the Schr\"odinger equation by transforming the Schr\"odinger equation into the Kirchhoff equation which describes the evolution of $n$ point vortices in Hydrodynamics. This helps us to view the Berry connection associated with Laughlin wave function or Schr\"odinger equation is not Hermitian, therefore we propose a self adjoint model Hamiltonian for the fractional quantum Hall effect, from the study of Kirchhoff equation, using superymmetric quantum mechanics whose solutions are complex Hermite polynomials. The Schr\"odinger equation as Berry connection allows us to formulate the problem of quantisation in terms of topology, as the quantum numbers are topological invariant which arise due to singularities in the Kirchhoff equation. Quantisation arises as it continuously connects the topologically inequivalent Hamiltonians in the Hilbert space.


title: Composition by Conversation

abstract: Most musical programming languages are developed purely for coding virtual instruments or algorithmic compositions. Although there has been some work in the domain of musical query languages for music information retrieval, there has been little attempt to unify the principles of musical programming and query languages with cognitive and natural language processing models that would facilitate the activity of composition by conversation. We present a prototype framework, called MusECI, that merges these domains, permitting score-level algorithmic composition in a text editor while also supporting connectivity to existing natural language processing frameworks.


title: Circumstellar Light Echo as a Possible Origin of the Polarization of Type IIP Supernovae

abstract: Type IIP supernovae (SNe IIP) are the most common class of core-collapse SNe. They often show rapid increase of polarization degree in the late phase. This time evolution is generally believed to originate from the emergence of an inner aspherical core, while an effect of polarized-scattered echoes by circumstellar (CS) dust around the SN may also substantially contribute to this polarization feature. In this study, we examine the effects of the scatted echoes on the SN polarization through radiative transfer simulations for various geometry and amount of CS dust. It has been found that asymmetrically-distributed CS dust, which is generally inferred for red supergiants, could reproduce the observed polarization features. We have applied our results to SNe 2004dj and 2006ov, deriving the geometry and amount of CS dust to explain their observed polarization features in this scenario. For both SNe, the blob-like or bipolar distribution of CS dust rather than the disk-like distribution is favored. The derived dust mass $M_{\mathrm{dust}}$ in the blob model (the bipolar CS dust model) for SNe 2004dj and 2006ov are $\sim 7.5 \times 10^{-4}$ M$_{\odot}$ ($\sim 8.5 \times10^{-4}$ M$_{\odot}$) and $\sim 5.2 \times 10^{-4}$ M$_{\odot}$ ($\sim 1.3 \times10^{-3}$ M$_{\odot}$), respectively. Even in the case where this process would not play a dominant role in the observed polarization signals, this effect should in principle contribute to it, the strength of which depends on the nature of CS dust. Therefore, this effect must be taken into account in discussing multi-dimensional structure of an SN explosion through polarimetric observations.


title: Density of a minimal submanifold and total curvature of its boundary

abstract: Given a piecewise smooth submanifold $\Gamma^{n-1} \subset \R^m$ and $p \in \R^m$, we define the {\em vision angle} $\Pi_p(\Gamma)$ to be the $(n-1)$-dimensional volume of the radial projection of $\Gamma$ to the unit sphere centered at $p$. If $p$ is a point on a stationary $n$-rectifiable set $\Sigma \subset \R^m$ with boundary $\Gamma$, then we show the density of $\Sigma$ at $p$ is $\leq$ the density at its vertex $p$ of the cone over $\Gamma$. It follows that if $\Pi_p(\Gamma)$ is less than twice the volume of $S^{n-1}$, for all $p \in \Gamma$, then $\Sigma$ is an embedded submanifold. As a consequence, we prove that given two $n$-planes $R^n_1, R^n_2$ in $\R^m$ and two compact convex hypersurfaces $\Gamma_i$ of $R^n_i, i=1,2$, a nonflat minimal submanifold spanned by $\Gamma:=\Gamma_1\cup\Gamma_2$ is embedded.


title: A Non-commutative Cryptosystem Based on Quaternion Algebras

abstract: We propose BQTRU, a non-commutative NTRU-like cryptosystem over quaternion algebras. This cryptosystem uses bivariate polynomials as the underling ring. The multiplication operation in our cryptosystem can be performed with high speed using quaternions algebras over finite rings. As a consequence, the key generation and encryption process of our cryptosystem is faster than NTRU in comparable parameters. Typically using Strassen's method, the key generation and encryption process is approximately $16/7$ times faster than NTRU for an equivalent parameter set. Moreover, the BQTRU lattice has a hybrid structure that makes inefficient standard lattice attacks on the private key. This entails a higher computational complexity for attackers providing the opportunity of having smaller key sizes. Consequently, in this sense, BQTRU is more resistant than NTRU against known attacks at an equivalent parameter set. Moreover, message protection is feasible through larger polynomials and this allows us to obtain the same security level as other NTRU-like cryptosystems but using lower dimensions.


title: Cooperative Curve Tracking in Two Dimensions Without Explicit Estimation of the Field Gradient

abstract: We design a control law for two agents to successfully track a level curve in the plane without explicitly estimating the field gradient. The velocity of each agent is decomposed along two mutually perpendicular directions, and separate control laws are designed along each direction. We prove that the formation center will converge to the neighborhood of the level curve with the desired level value. The algorithm is tested on some test functions used in optimization problems in the presence of noise. Our results indicate that in spite of the control law being simple and gradient-free, we are able to successfully track noisy planar level curves fast and with a high degree of accuracy.


title: An unsupervised long short-term memory neural network for event detection in cell videos

abstract: We propose an automatic unsupervised cell event detection and classification method, which expands convolutional Long Short-Term Memory (LSTM) neural networks, for cellular events in cell video sequences. Cells in images that are captured from various biomedical applications usually have different shapes and motility, which pose difficulties for the automated event detection in cell videos. Current methods to detect cellular events are based on supervised machine learning and rely on tedious manual annotation from investigators with specific expertise. So that our LSTM network could be trained in an unsupervised manner, we designed it with a branched structure where one branch learns the frequent, regular appearance and movements of objects and the second learns the stochastic events, which occur rarely and without warning in a cell video sequence. We tested our network on a publicly available dataset of densely packed stem cell phase-contrast microscopy images undergoing cell division. This dataset is considered to be more challenging that a dataset with sparse cells. We compared our method to several published supervised methods evaluated on the same dataset and to a supervised LSTM method with a similar design and configuration to our unsupervised method. We used an F1-score, which is a balanced measure for both precision and recall. Our results show that our unsupervised method has a higher or similar F1-score when compared to two fully supervised methods that are based on Hidden Conditional Random Fields (HCRF), and has comparable accuracy with the current best supervised HCRF-based method. Our method was generalizable as after being trained on one video it could be applied to videos where the cells were in different conditions. The accuracy of our unsupervised method approached that of its supervised counterpart.


title: A deep generative model for gene expression profiles from single-cell RNA sequencing

abstract: We propose a probabilistic model for interpreting gene expression levels that are observed through single-cell RNA sequencing. In the model, each cell has a low-dimensional latent representation. Additional latent variables account for technical effects that may erroneously set some observations of gene expression levels to zero. Conditional distributions are specified by neural networks, giving the proposed model enough flexibility to fit the data well. We use variational inference and stochastic optimization to approximate the posterior distribution. The inference procedure scales to over one million cells, whereas competing algorithms do not. Even for smaller datasets, for several tasks, the proposed procedure outperforms state-of-the-art methods like ZIFA and ZINB-WaVE. We also extend our framework to account for batch effects and other confounding factors, and propose a Bayesian hypothesis test for differential expression that outperforms DESeq2.


title: Radiating fluid sphere immersed in an anisotropic atmosphere

abstract: We model a radiating star undergoing dissipative gravitational collapse in the form of radial heat flux. The exterior of the collapsing star is described by the generalised Vaidya solution representing a mixture of null radiation and strings. Our model generalises previously known results of constant string density atmosphere to include inhomogeneities in the exterior spacetime. By utilising a causal heat transport equation of the Maxwell-Cattaneo form we show that relaxational effects are enhanced in the presence of inhomogeneities due to the string density.


title: Field-Induced Instability of a Gapless Spin Liquid with a Spinon Fermi Surface

abstract: The ground state of the quantum kagome antiferromagnet Zn-brochantite, ZnCu$_3$(OH)$_6$SO$_4$, which is one of only a few known spin-liquid (SL) realizations in two or three dimensions, has been described as a gapless SL with a spinon Fermi surface. Employing nuclear magnetic resonance in a broad magnetic-field range down to millikelvin temperatures, we show that in applied magnetic fields this enigmatic state is intrinsically unstable against a SL with a full or a partial gap. A similar instability of the gapless Fermi-surface SL was previously encountered in an organic triangular-lattice antiferromagnet, suggesting a common destabilization mechanism that most likely arises from spinon pairing. A salient property of this instability is that an infinitesimal field suffices to induce it, as predicted theoretically for some other types of gapless SL's.


title: Short Wavelength Geodesic Acoustic Mode Excitation by Energetic Particles

abstract: Taking the collisionless damping of geodesic acoustic mode (GAM) as an example, the physics processes underlying wave particle resonances in the short wavelength limit are clarified. As illus- trative application, GAM excitation by energetic particles in short wavelength limit is investigated assuming a single pitch angle slowing-down fast ion equilibrium distribution function. Conditions for this energetic particle-induced GAM (EGAM) to be unstable are discussed.


title: Sampling for approximating $R$-limited functions

abstract: $R$-limited functions are multivariate generalization of band-limited functions whose Fourier transforms are supported within a compact region $R\subset\mathbb{R}^{n}$. In this work, we generalize sampling and interpolation theorems for band-limited functions to $R$-limited functions. More precisely, we investigated the following question: "For a function compactly supported within a region similar to $R$, does there exist an $R$-limited function that agrees with the function over its support for a desired accuracy?". Starting with the Fourier domain definition of an $R$-limited function, we write the equivalent convolution and a discrete Fourier transform representations for $R$-limited functions through approximation of the convolution kernel using a discrete subset of Fourier basis. The accuracy of the approximation of the convolution kernel determines the accuracy of the discrete Fourier representation. Construction of the discretization can be achieved using the tools from approximation theory as demonstrated in the appendices. The main contribution of this work is proving the equivalence between the discretization of the Fourier and convolution representations of $R$-limited functions. Here discrete convolution representation is restricted to shifts over a compactly supported region similar to $R$. We show that discrete shifts for the convolution representation are equivalent to the spectral parameters used in discretization of the Fourier representation of the convolution kernel. This result is a generalization of the cardinal theorem of interpolation of band-limited functions. The error corresponding to discrete convolution representation is also bounded by the approximation of the convolution kernel using discretized Fourier basis.


title: Sharp Bounds for Generalized Uniformity Testing

abstract: We study the problem of generalized uniformity testing \cite{BC17} of a discrete probability distribution: Given samples from a probability distribution $p$ over an {\em unknown} discrete domain $\mathbf{\Omega}$, we want to distinguish, with probability at least $2/3$, between the case that $p$ is uniform on some {\em subset} of $\mathbf{\Omega}$ versus $\epsilon$-far, in total variation distance, from any such uniform distribution. We establish tight bounds on the sample complexity of generalized uniformity testing. In more detail, we present a computationally efficient tester whose sample complexity is optimal, up to constant factors, and a matching information-theoretic lower bound. Specifically, we show that the sample complexity of generalized uniformity testing is $\Theta\left(1/(\epsilon^{4/3}\|p\|_3) + 1/(\epsilon^{2} \|p\|_2) \right)$.

categories: cs.DS cs.IT cs.LG math.IT math.ST stat.TH


title: On the torsion subgroups of the modular Jacobians

abstract: For any positive integer $N$, we prove that the rational torsion subgroup of $J_0(N)$ agrees with its rational cuspidal subgroups up to a factor of $6N\prod_{p\mid N}(p^2-1)$. Moreover, for modular Jacobians of the form $J_0(DC)$ with $D$ a positive square-free integer and $C$ any positive divisor of $D$, we prove that the $\psi$-part of the torsion subgroup of $J_0(DC)$ agrees with the $\psi$-part of its cuspidal subgroup up to a factor of $6D\prod_{p\mid D}(p^2-1)$, where $\psi$ is any quadratic character of conductor dividing $C$.


title: Approximating fractional derivative of the Gaussian function and Dawson's integral

abstract: A new method for approximating fractional derivatives of the Gaussian function and Dawson's integral are presented. Unlike previous approaches, which are dominantly based on some discretization of Riemann-Liouville integral using polynomial or discrete Fourier basis, we take an alternative approach which is based on expressing the Riemann-Liouville definition of the fractional integral for the semi-infinite axis in terms of a moment problem. As a result, fractional derivatives of the Gaussian function and Dawson's integral are expressed as a weighted sum of complex scaled Gaussian and Dawson's integral. Error bounds for the approximation are provided. Another distinct feature of the proposed method compared to the previous approaches, it can be extended to approximate partial derivative with respect to the order of the fractional derivative which may be used in PDE constraint optimization problems.


title: Relation of $a^\dagger a$ terms to higher-order terms in the adiabatic expansion for large-amplitude collective motion

abstract: We investigate the relation of $a^\dagger a$ terms in the collective operator to the higher-order terms in the adiabatic self-consistent collective coordinate (ASCC) method. In the ASCC method, a state vector is written as $e^{i\hat G(q,p,n)}|\phi(q)\rangle$ with $\hat G(q,p,n)$ which is a function of collective coordinate $q$, its conjugate momentum $p$ and the particle number $n$. According to the generalized Thouless theorem, $\hat G$ can be written as a linear combination of two-quasiparticle creation and annihilation operators $a^\dagger_\mu a^\dagger_\nu$ and $a_\nu a_\mu$. We show that, if $a^\dagger a$ terms are included in $\hat G(q,p,n)$, it corresponds to the higher-order terms in the adiabatic expansion of $\hat G$. This relation serves as a prescription to determine the higher-order collective operators from the $a^\dagger a$ part of the collective operator, once it is given without solving the higher-order equations of motion.


title: Asynchronous COMID: the theoretic basis for transmitted data sparsification tricks on Parameter Server

abstract: Asynchronous FTRL-proximal and L2 norm done at server are two widely used tricks in Parameters Server which is an implement of delayed SGD. Their commonness is leaving parts of updating computation on server which reduces the burden of network via making transmitted data sparse. But above tricks' convergences are not well-proved. In this paper, based on above commonness, we propose a more general algorithm named as asynchronous COMID and prove its convergence. We prove that asynchronous FTRL-proximal and L2 norm done at server are applications of asynchronous COMID, which demonstrates the convergences of above two tricks. Then, we conduct experiments to verify theoretical results. Experimental results show that compared with delayed SGD on Parameters Server, asynchronous COMID reduces the burden of the network without any harm on the mathematical convergence speed and final output.


title: Probabilistic Analysis Based On Symbolic Game Semantics and Model Counting

abstract: Probabilistic program analysis aims to quantify the probability that a given program satisfies a required property. It has many potential applications, from program understanding and debugging to computing program reliability, compiler optimizations and quantitative information flow analysis for security. In these situations, it is usually more relevant to quantify the probability of satisfying/violating a given property than to just assess the possibility of such events to occur. In this work, we introduce an approach for probabilistic analysis of open programs (i.e. programs with undefined identifiers) based on game semantics and model counting. We use a symbolic representation of algorithmic game semantics to collect the symbolic constraints on the input data (context) that lead to the occurrence of the target events (e.g. satisfaction/violation of a given property). The constraints are then analyzed to quantify how likely is an input to satisfy them. We use model counting techniques to count the number of solutions (from a bounded integer domain) that satisfy given constraints. These counts are then used to assign probabilities to program executions and to assess the probability for the target event to occur at the desired level of confidence. Finally, we present the results of applying our approach to several interesting examples and illustrate the benefits they may offer.


title: Parametric Markov Chains: PCTL Complexity and Fraction-free Gaussian Elimination

abstract: Parametric Markov chains have been introduced as a model for families of stochastic systems that rely on the same graph structure, but differ in the concrete transition probabilities. The latter are specified by polynomial constraints for the parameters. Among the tasks typically addressed in the analysis of parametric Markov chains are (1) the computation of closed-form solutions for reachabilty probabilities and other quantitative measures and (2) finding symbolic representations of the set of parameter valuations for which a given temporal logical formula holds as well as (3) the decision variant of (2) that asks whether there exists a parameter valuation where a temporal logical formula holds. Our contribution to (1) is to show that existing implementations for computing rational functions for reachability probabilities or expected costs in parametric Markov chains can be improved by using fraction-free Gaussian elimination, a long-known technique for linear equation systems with parametric coefficients. Our contribution to (2) and (3) is a complexity-theoretic discussion of the model checking problem for parametric Markov chains and probabilistic computation tree logic (PCTL) formulas. We present an exponential-time algorithm for (2) and a PSPACE upper bound for (3). Moreover, we identify fragments of PCTL and subclasses of parametric Markov chains where (1) and (3) are solvable in polynomial time and establish NP-hardness for other PCTL fragments.


title: On the Complexity of Model Checking for Syntactically Maximal Fragments of the Interval Temporal Logic HS with Regular Expressions

abstract: In this paper, we investigate the model checking (MC) problem for Halpern and Shoham's interval temporal logic HS. In the last years, interval temporal logic MC has received an increasing attention as a viable alternative to the traditional (point-based) temporal logic MC, which can be recovered as a special case. Most results have been obtained under the homogeneity assumption, that constrains a proposition letter to hold over an interval if and only if it holds over each component state. Recently, Lomuscio and Michaliszyn proposed a way to relax such an assumption by exploiting regular expressions to define the behaviour of proposition letters over intervals in terms of their component states. When homogeneity is assumed, the exact complexity of MC is a difficult open question for full HS and for its two syntactically maximal fragments AA'BB'E' and AA'EB'E'. In this paper, we provide an asymptotically optimal bound to the complexity of these two fragments under the more expressive semantic variant based on regular expressions by showing that their MC problem is AEXP_pol-complete, where AEXP_pol denotes the complexity class of problems decided by exponential-time bounded alternating Turing Machines making a polynomially bounded number of alternations.


title: Thermodynamics and relativistic kinetic theory for q-generalized Bose-Einstein and Fermi-Dirac systems

abstract: The thermodynamics and covariant kinetic theory have been elaborately investigated in a non-extensive environment considering the non-extensive generalization of Bose-Einstein (BE) and Fermi-Dirac (FD) statistics. Starting with Tsallis' entropy formula, the fundamental principles of thermostatistics have been established for a grand canonical system having q-generalized BE/FD degrees of freedom. The many particle kinetic theory has been set up in terms of the relativistic transport equation with q-generalized Uehling-Uhlenbeck collision term. The conservation laws have been realized in terms of appropriate moments of the transport equation. The thermodynamic quantities have been obtained in weak non-extensive environment for a massive pion-nucleon and a massless quark-gluon system with non-zero baryon chemical potential. In order to get an estimate of the impact of non-extensivity on the system dynamics, the q-modified Debye mass and hence the q-modified effective coupling have been estimated for a quark-gluon system.


title: An Existence Theorem of Nash Equilibrium in Coq and Isabelle

abstract: Nash equilibrium (NE) is a central concept in game theory. Here we prove formally a published theorem on existence of an NE in two proof assistants, Coq and Isabelle: starting from a game with finitely many outcomes, one may derive a game by rewriting each of these outcomes with either of two basic outcomes, namely that Player 1 wins or that Player 2 wins. If all ways of deriving such a win/lose game lead to a game where one player has a winning strategy, the original game also has a Nash equilibrium. This article makes three other contributions: first, while the original proof invoked linear extension of strict partial orders, here we avoid it by generalizing the relevant definition. Second, we notice that the theorem also implies the existence of a secure equilibrium, a stronger version of NE that was introduced for model checking. Third, we also notice that the constructive proof of the theorem computes secure equilibria for non-zero-sum priority games (generalizing parity games) in quasi-polynomial time.


title: The Satisfiability Problem for Boolean Set Theory with a Choice Correspondence

abstract: Given a set U of alternatives, a choice (correspondence) on U is a contractive map c defined on a family Omega of nonempty subsets of U. Semantically, a choice c associates to each menu A in Omega a nonempty subset c(A) of A comprising all elements of A that are deemed selectable by an agent. A choice on U is total if its domain is the powerset of U minus the empty set, and partial otherwise. According to the theory of revealed preferences, a choice is rationalizable if it can be retrieved from a binary relation on U by taking all maximal elements of each menu. It is well-known that rationalizable choices are characterized by the satisfaction of suitable axioms of consistency, which codify logical rules of selection within menus. For instance, WARP (Weak Axiom of Revealed Preference) characterizes choices rationalizable by a transitive relation. Here we study the satisfiability problem for unquantified formulae of an elementary fragment of set theory involving a choice function symbol c, the Boolean set operators and the singleton, the equality and inclusion predicates, and the propositional connectives. In particular, we consider the cases in which the interpretation of c satisfies any combination of two specific axioms of consistency, whose conjunction is equivalent to WARP. In two cases we prove that the related satisfiability problem is NP-complete, whereas in the remaining cases we obtain NP-completeness under the additional assumption that the number of choice terms is constant.


title: MK-fuzzy Automata and MSO Logics

abstract: We introduce MK-fuzzy automata over a bimonoid K which is related to the fuzzification of the McCarthy-Kleene logic. Our automata are inspired by, and intend to contribute to, practical applications being in development in a project on runtime network monitoring based on predicate logic. We investigate closure properties of the class of recognizable MK-fuzzy languages accepted by MK-fuzzy automata as well as of deterministically recognizable MK-fuzzy languages accepted by their deterministic counterparts. Moreover, we establish a Nivat-like result for recognizable MK-fuzzy languages. We introduce an MK-fuzzy MSO logic and show the expressive equivalence of a fragment of this logic with MK-fuzzy automata, i.e., a B\"uchi type theorem.


title: Robust Exponential Worst Cases for Divide-et-Impera Algorithms for Parity Games

abstract: The McNaughton-Zielonka divide et impera algorithm is the simplest and most flexible approach available in the literature for determining the winner in a parity game. Despite its theoretical worst-case complexity and the negative reputation as a poorly effective algorithm in practice, it has been shown to rank among the best techniques for the solution of such games. Also, it proved to be resistant to a lower bound attack, even more than the strategy improvements approaches, and only recently a family of games on which the algorithm requires exponential time has been provided by Friedmann. An easy analysis of this family shows that a simple memoization technique can help the algorithm solve the family in polynomial time. The same result can also be achieved by exploiting an approach based on the dominion-decomposition techniques proposed in the literature. These observations raise the question whether a suitable combination of dynamic programming and game-decomposition techniques can improve on the exponential worst case of the original algorithm. In this paper we answer this question negatively, by providing a robustly exponential worst case, showing that no intertwining of the above mentioned techniques can help mitigating the exponential nature of the divide et impera approaches.


title: Dynamics and Coalitions in Sequential Games

abstract: We consider N-player non-zero sum games played on finite trees (i.e., sequential games), in which the players have the right to repeatedly update their respective strategies (for instance, to improve the outcome wrt to the current strategy profile). This generates a dynamics in the game which may eventually stabilise to a Nash Equilibrium (as with Kukushkin's lazy improvement), and we argue that it is interesting to study the conditions that guarantee such a dynamics to terminate. We build on the works of Le Roux and Pauly who have studied extensively one such dynamics, namely the Lazy Improvement Dynamics. We extend these works by first defining a turn-based dynamics, proving that it terminates on subgame perfect equilibria, and showing that several variants do not terminate. Second, we define a variant of Kukushkin's lazy improvement where the players may now form coalitions to change strategies. We show how properties of the players' preferences on the outcomes affect the termination of this dynamics, and we thereby characterise classes of games where it always terminates (in particular two-player games).


title: A Parallel Linear Temporal Logic Tableau

abstract: For many applications, we are unable to take full advantage of the potential massive parallelisation offered by supercomputers or cloud computing because it is too hard to work out how to divide up the computation task between processors in such a way to minimise the need for communication. However, a recently developed branch-independent tableaux for the common LTL temporal logic should intuitively be easy to parallelise as each branch can be developed independently. Here we describe a simple technique for partitioning such a tableau such that each partition can be processed independently without need for interprocess communication. We investigate the extent to which this technique improves the performance of the LTL tableau on standard benchmarks and random formulas.


title: LTL to Deterministic Emerson-Lei Automata

abstract: We introduce a new translation from linear temporal logic (LTL) to deterministic Emerson-Lei automata, which are omega-automata with a Muller acceptance condition symbolically expressed as a Boolean formula. The richer acceptance condition structure allows the shift of complexity from the state space to the acceptance condition. Conceptually the construction is an enhanced product construction that exploits knowledge of its components to reduce the number of states. We identify two fragments of LTL, for which one can easily construct deterministic automata and show how knowledge of these components can reduce the number of states. We extend this idea to a general LTL framework, where we can use arbitrary LTL to deterministic automata translators for parts of formulas outside the mentioned fragments. Further, we show succinctness of the translation compared to existing construction. The construction is implemented in the tool Delag, which we evaluate on several benchmarks of LTL formulas and probabilistic model checking case studies.


title: Linear-time Temporal Logic with Event Freezing Functions

abstract: Formal properties represent a cornerstone of the system-correctness proofs based on formal verification techniques such as model checking. Formalizing requirements into temporal properties may be very complex and error prone, due not only to the ambiguity of the textual requirements but also to the complexity of the formal language. Finding a property specification language that balances simplicity, expressiveness, and tool support remains an open problem in many real-world contexts. In this paper, we propose a new temporal logic, which extends First-Order Linear-time Temporal Logic with Past adding two operators "at next" and "at last", which take in input a term and a formula and represent the value of the term at the next state in the future or last state in the past in which the formula holds. We consider different models of time (including discrete, dense, and super-dense time) and Satisfiability Modulo Theories (SMT) of the first-order formulas. The "at next" and "at last" functions can be seen as a generalization of Event-Clock operators and can encode some Metric Temporal operators also with counting. They are useful to formalize properties of component-based models because they allow to express constraints on the data exchanged with messages at different instants of time. We provide a simple encoding into equisatisfiable formulas without the extra functional symbols. We implement a prototype tool support based on SMT-based model checking.


title: Beyond $\omega$BS-regular Languages: $\omega$T-regular Expressions and Counter-Check Automata

abstract: In the last years, various extensions of {\omega}-regular languages have been proposed in the literature, including {\omega}B-regular ({\omega}-regular languages extended with boundedness), {\omega}S-regular ({\omega}-regular languages extended with strict unboundedness), and {\omega}BS-regular languages (the combination of {\omega}B- and {\omega}S-regular ones). While the first two classes satisfy a generalized closure property, namely, the complement of an {\omega}B-regular (resp., {\omega}S-regular) language is an {\omega}S-regular (resp., {\omega}B-regular) one, the last class is not closed under complementation. The existence of non-{\omega}BS-regular languages that are the complements of some {\omega}BS-regular ones and express fairly natural properties of reactive systems motivates the search for other well-behaved classes of extended {\omega}-regular languages. In this paper, we introduce the class of {\omega}T-regular languages, that includes meaningful languages which are not {\omega}BS-regular. We first define it in terms of {\omega}T-regular expressions. Then, we introduce a new class of automata (counter-check automata) and we prove that (i) their emptiness problem is decidable in PTIME and (ii) they are expressive enough to capture {\omega}T-regular languages (whether or not {\omega}T-regular languages are expressively complete with respect to counter-check automata is still an open problem). Finally, we provide an encoding of {\omega}T-regular expressions into S1S+U.


title: Model Checking Social Network Models

abstract: A social network service is a platform to build social relations among people sharing similar interests and activities. The underlying structure of a social networks service is the social graph, where nodes represent users and the arcs represent the users' social links and other kind of connections. One important concern in social networks is privacy: what others are (not) allowed to know about us. The "logic of knowledge" (epistemic logic) is thus a good formalism to define, and reason about, privacy policies. In this paper we consider the problem of verifying knowledge properties over social network models (SNMs), that is social graphs enriched with knowledge bases containing the information that the users know. More concretely, our contributions are: i) We prove that the model checking problem for epistemic properties over SNMs is decidable; ii) We prove that a number of properties of knowledge that are sound w.r.t. Kripke models are also sound w.r.t. SNMs; iii) We give a satisfaction-preserving encoding of SNMs into canonical Kripke models, and we also characterise which Kripke models may be translated into SNMs; iv) We show that, for SNMs, the model checking problem is cheaper than the one based on standard Kripke models. Finally, we have developed a proof-of-concept implementation of the model-checking algorithm for SNMs.


title: A Backward-traversal-based Approach for Symbolic Model Checking of Uniform Strategies for Constrained Reachability

abstract: Since the introduction of Alternating-time Temporal Logic (ATL), many logics have been proposed to reason about different strategic capabilities of the agents of a system. In particular, some logics have been designed to reason about the uniform memoryless strategies of such agents. These strategies are the ones the agents can effectively play by only looking at what they observe from the current state. ATL_ir can be seen as the core logic to reason about such uniform strategies. Nevertheless, its model-checking problem is difficult (it requires a polynomial number of calls to an NP oracle), and practical algorithms to solve it appeared only recently. This paper proposes a technique for model checking uniform memoryless strategies. Existing techniques build the strategies from the states of interest, such as the initial states, through a forward traversal of the system. On the other hand, the proposed approach builds the winning strategies from the target states through a backward traversal, making sure that only uniform strategies are explored. Nevertheless, building the strategies from the ground up limits its applicability to constrained reachability objectives only. This paper describes the approach in details and compares it experimentally with existing approaches implemented into a BDD-based framework. These experiments show that the technique is competitive on the cases it can handle.


title: On the Complexity of ATL and ATL* Module Checking

abstract: Module checking has been introduced in late 1990s to verify open systems, i.e., systems whose behavior depends on the continuous interaction with the environment. Classically, module checking has been investigated with respect to specifications given as CTL and CTL* formulas. Recently, it has been shown that CTL (resp., CTL*) module checking offers a distinctly different perspective from the better-known problem of ATL (resp., ATL*) model checking. In particular, ATL (resp., ATL*) module checking strictly enhances the expressiveness of both CTL (resp., CTL*) module checking and ATL (resp. ATL*) model checking. In this paper, we provide asymptotically optimal bounds on the computational cost of module checking against ATL and ATL*, whose upper bounds are based on an automata-theoretic approach. We show that module-checking for ATL is EXPTIME-complete, which is the same complexity of module checking against CTL. On the other hand, ATL* module checking turns out to be 3EXPTIME-complete, hence exponentially harder than CTL* module checking.


title: ParaPlan: A Tool for Parallel Reachability Analysis of Planar Polygonal Differential Inclusion Systems

abstract: We present the ParaPlan tool which provides the reachability analysis of planar hybrid systems defined by differential inclusions (SPDI). It uses the parallelized and optimized version of the algorithm underlying the SPeeDI tool. The performance comparison demonstrates the speed-up of up to 83 times with respect to the sequential implementation on various benchmarks. Some of the benchmarks we used are randomly generated with the novel approach based on the partitioning of the plane with Voronoi diagrams.


title: Machine learning for analysis of plasma driven Ion source

abstract: Recently, neural networks have found many applications in different fields including Genetics, Pharmacy, Astrophysics and High Energy Physics [1-3]. In the field of accelerator physics it has been used for control systems [4]. In this paper we present the results based on machine learning techniques motivated to predict the behaviour of ion source in terms of composition of the ion beam while using hydrogen gas to produce $H^+$ ions. In the framework of the stellarator type Figure-8 Storage Ring (F8SR) project, a volume type ion source was designed for the low energy ion beam transport experiments. In a first step the functioning of this ion source was studied and the results were published, but only small number of measurements were analysed as the main requirement for the on going experiment was fulfilled. Though at a later stage, more number of measurements were recorded with larger parameter space to investigate the properties of extracted ion beams from this source further. With recent interests and improved techniques in the applications of machine learning algorithms data analysis using neural network has been applied to study the ion beams from this ion source.


title: Bosonization of Cooper pairs and novel Bose-liquid superconductivity and superfluidity in high-$T_c$ cuprates and other systems

abstract: The criteria for bosonization of Cooper pairs and novel Bose-liquid superconductivity and superfluidity in high-$T_c$ cuprates and other pseudogap matters are formulated by using the uncertainty principle and the boson mean field theory. We argue that the BCS-type s- or d-wave superconductivity occurring in the fermionic limit of Cooper pairs (which exist in ordinary metals and heavily overdoped cuprates with large Fermi energies) is not characteristic of underdoped to overdoped cuprates with low Fermi energies. The superconducting order parameter in high-$T_c$ cuprates and other pseudogap matters cannot be determined as the BCS-like (s- or d-wave) gap. We show that the unconventional superconductivity/superfluidity occurring in the bosonic limit of Cooper pairs exists in low Fermi energy systems where the bosonic Cooper pairs are formed at a pseudogap temperature T* above the uperconducting/superfluid transition temperature $T_c$ and then such bosons condenses into a Bose superfluid at $T_c$. The pair condensation of attracting bosons occurs at $T_c$ and then their single particle condensation sets in at $T_c^*$ lower than $T_c$ (in three dimensions (3D)) or at T=0 (in two dimensions (2D)). The coherent single particle and pair condensates of bosons exist as two distinct superfluid phases below $T_c$. By solving the mean field equations for high-$T_c$ cuprates, the novel superconducting states (i.e., a vortex-like state existing below the temperature $T_v=T^{2D}_c$ lower than T* but higher than $T_c=T^{3D}_c$ and two superconducting phases below $T_c$) and their properties characterized by the boson superfluid stiffness are determined and compared with the key experimental findings. The novel Bose-liquid superconductivity and superfluidity exists also in other exotic superconductors/superfluids, including quantum liquids, atomic Fermi gases and low-density nuclear matter.


title: Tensile elasticity of semiflexible polymers with hinge defects

abstract: It has become clear in recent years that the simple uniform wormlike chain model needs to be modified in order to account for more complex behavior which has been observed experimentally in some important biopolymers. For example, the large flexibility of short ds-DNA has been attributed to kink or hinge defects. In this paper, we calculate analytically, within the weak bending approximation, the force-extension relation of a wormlike chain with a permanent hinge defect along its contour. The defect is characterized by its bending energy (which can be zero, in the completely flexible case) and its position along the polymer contour. Besides the bending rigidity of the chain, these are the only parameters which describe our model. We show that a hinge defect causes a significant increase in the differential tensile compliance of a pre-stressed chain. In the small force limit, a hinge defect significantly increases the entropic elasticity. Our results apply to any pair of semiflexible segments connected by a hinge. As such, they may also be relevant to cytoskeletal filaments (F-actin, microtubules), where one may treat the cross-link connecting two filaments as a hinge defect.


title: A bathtub model for the star-forming interstellar medium

abstract: The bathtub model of the star forming interstellar medium is based on the powerful constraint that mass has to be conserved when gas flows through its various thermal and density phases, ending up eventually in a young star or being blown away by stellar feedback. It predicts that the star formation rate of a molecular cloud is not determined by the cloud's mass or its internal collapse timescale, but rather by the accretion rate of new gas. For the most simple case of a constant accretion flow an equilibrium state is reached quickly where the star formation rate equals the accretion rate and where the dense gas mass is constant and independent of time. The mass of the young star cluster, on the other hand, increases linearly with time. The stellar mass fraction therefore represents a sensitive clock to measure the age of the star-forming region. The bathtub model predicts that the efficiency of star formation is small, of order 1%, even in the dense filamentary phases of molecular clouds. It provides a simple explanation for the dense gas fraction of order 10% in molecular clouds and for the large gas depletion timescales of star-forming galaxies of order $5 \times 10^8 - 10^9$ yrs.


title: Truncation in Average and Worst Case Settings for Special Classes of $\infty$-Variate Functions

abstract: The paper considers truncation errors for functions of the form $f(x_1,x_2,\dots)=g(\sum_{j=1}^\infty x_j\,\xi_j)$, i.e., errors of approximating $f$ by $f_k(x_1,\dots,x_k)=g(\sum_{j=1}^k x_j\,\xi_j)$, where the numbers $\xi_j$ converge to zero sufficiently fast and $x_j$'s are i.i.d. random variables. As explained in the introduction, functions $f$ of the form above appear in a number of important applications. To have positive results for possibly large classes of such functions, the paper provides sharp bounds on truncation errors in both the average and worst case settings. In the former case, the functions $g$ are from a Hilbert space $G$ endowed with a zero mean probability measure with a given covariance kernel. In the latter case, the functions $g$ are from a reproducing kernel Hilbert space, or a space of functions satisfying a H\"older condition.


title: Osmotic pressure between arbitrarily charged surfaces: a revisited approach

abstract: The properties of ionic solutions between charged surfaces are often studied within the Poisson-Boltzmann framework, by finding the electrostatic potential profile. For example, the osmotic pressure between two charged planar surfaces can be evaluated by solving coupled equations for the electrostatic potential and osmotic pressure. Such a solution relies on symmetry arguments and is restricted to either equally or oppositely charged surfaces. Here, we provide a different and more efficient scheme to derive the osmotic pressure straight-forwardly, without the need to find the electrostatic potential profile. We derive analytical expressions for the osmotic pressure in terms of the inter-surface separation, salt concentration, and arbitrary boundary conditions. Such results should be useful in force measurement setups, where the force is measured between two differently prepared surfaces, or between two surfaces held at a fixed potential difference. The proposed method can be systematically used for generalized Poisson-Boltzmann theories in planar geometries, as is demonstrated for the sterically modified Poisson-Boltzmann theory.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.chem-ph


title: Smoothness of Flow and Path-by-Path Uniqueness in Stochastic Differential Equations

abstract: We consider the stochastic differential equation $$ X_t = x_0 + \int_0^t f(X_s)ds + \int_0^t\sigma(X_s)dB^{H}_s,$$ with $x_0 \in \mathbb{R}^d$, $d \geq 1$, $f: \mathbb{R}^d \rightarrow \mathbb{R}^d$ is bounded continuous, $\sigma: \mathbb{R}^d \rightarrow \mathbb{R}^{d\times d}$ is a uniformly elliptic, bounded, twice continuously differentiable conservative vector field and $B^H$ is fractional Brownian motion with $H \in (\frac{1}{3}, \frac{1}{2}]$. When $d=1$, $H= \frac{1}{2}$, and $f$ is H\"older continuous, in the spirit of Davie [D07], we establish the existence of a null set $\mathcal{N}$ depending only on $f, \sigma$ such that for all $x_0\in \mathbb{R}$ and $\omega \in \Omega\setminus \mathcal{N}$, the above equation admits a path-by-path unique solution. Our proof is based on establishing the uniform continuous differentiability of the flow associated with the equation. We also establish the path-by-path uniqueness for $d \geq 1$ and $H \in (\frac{1}{3}, \frac{1}{2}]$, but the null set may depend on $x_0$, thus extending a result of Catellier-Gubinelli [CG12].


title: Unreported links between trial registrations and published articles were identified using document similarity measures in a cross-sectional analysis of ClinicalTrials.gov

abstract: Objectives: Trial registries can be used to measure reporting biases and support systematic reviews but 45% of registrations do not provide a link to the article reporting on the trial. We evaluated the use of document similarity methods to identify unreported links between ClinicalTrials.gov and PubMed. Study Design and Setting: We extracted terms and concepts from a dataset of 72,469 ClinicalTrials.gov registrations and 276,307 PubMed articles, and tested methods for ranking articles across 16,005 reported links and 90 manually-identified unreported links. Performance was measured by the median rank of matching articles, and the proportion of unreported links that could be found by screening ranked candidate articles in order. Results: The best performing concept-based representation produced a median rank of 3 (IQR 1-21) for reported links and 3 (IQR 1-19) for the manually-identified unreported links, and term-based representations produced a median rank of 2 (1-20) for reported links and 2 (IQR 1-12) in unreported links. The matching article was ranked first for 40% of registrations, and screening 50 candidate articles per registration identified 86% of the unreported links. Conclusions: Leveraging the growth in the corpus of reported links between ClinicalTrials.gov and PubMed, we found that document similarity methods can assist in the identification of unreported links between trial registrations and corresponding articles.


title: Metric methods for heteroclinic connections in infinite dimensional spaces

abstract: We consider the minimal action problem min \int\_R 1/2 |$\gamma$'|^2 + W($\gamma$) dt among curves lying in a non-locally-compact metric space and connecting two given zeros of W $\ge$ 0. For this problem, the optimal curves are usually called heteroclinic connections. We reduce it, following a standard method, to a geodesic problem of the form min \int\_0^1 K($\gamma$)|$\gamma$'| dt with K = (2W)^(1/2). We then prove existence of curves minimizing this new action under some suitable compactness assumptions on K, which are minimal. The method allows to solve some PDE problems in unbounded domains, in particular in two variables x, y, when y = t and when the metric space is an L^2 space in the first variable x, and the potential W includes a Dirichlet energy in the same variable. We then apply this technique to the problem of connecting, in a functional space, two different heteroclinic connections between two points of the Euclidean space, as it was previously studied by Alama-Bronsard-Gui and by Schatzman more than fifteen years ago. With a very different technique, we are able to recover the same results, and to weaken some assumptions.


title: Detecting a hidden obstacle via the time domain enclosure method. A scalar wave case

abstract: The characterization problem of the existence of an unknown obstacle behind a known obstacle is considered by using a singe observed wave at a place where the wave is generated. The unknown obstacle is invisible from the place by using visible ray. A mathematical formulation of the problem using the classical wave equation is given. The main result consists of two parts: (i) one can make a decision whether the unknown obstacle exists or not behind a known impenetrable obstacle by using a single wave over a finite time interval under some a-priori information on the position of the unknown obstacle; (ii) one can obtain a lower bound of the Euclidean distance of the unknown obstacle to the center point of the support of the initial data of the wave. The proof is based on the idea of the time domain enclosure method and employs some previous results on the Gaussian lower/upper estimates for the heat kernels and domination of semigroups.


title: Crystal field effects on spin pumping

abstract: "Spin pumping" is the injection of spin angular momentum by a time-dependent magnetization into an adjacent normal metal proportional to the spin mixing conductance. We study the role of electrostatic interactions in the form of crystal fields on the pumped spin currents generated by insulators with exchange-coupled local moments at the interface to a metal. The crystal field is shown to render the spin currents anisotropic, which implies that the spin mixing conductance of insulator|normal metal bilayers depends on crystal cut and orientation. We interpret the interface "effective field" (imaginary part of the spin mixing conductance) in terms of the coherent motion of the equilibrium spin density induced by proximity in the normal metal.


title: Updated three-body model of $^6$He $\beta$ decay into the $\alpha$ + $d$ continuum

abstract: The $\beta$-decay process of the $^6$He halo nucleus into the $\alpha+d$ continuum is studied in an updated three-body model. The $^6$He nucleus is described as an $\alpha+n+n$ system in hyperspherical coordinates on a Lagrange-mesh. The shape and absolute values of the transition probability per time and energy units of new experiments are reproduced with a modified $\alpha+d$ potential. The obtained total transition probabilities are $2.48 \times 10^{-6}$ s$^{-1}$ for the full energy region and $2.40 \times 10^{-6}$ s$^{-1}$ for the cut-off $E>150$ keV. The strong cancellation between the internal and halo parts of the $\beta$ decay matrix element is a challenge for future {\it ab initio} calculations.


title: The Wigner branching random walk: Efficient implementation and performance evaluation

abstract: To implement the Wigner branching random walk, the particle carrying a signed weight, either $-1$ or $+1$, is more friendly to data storage and arithmetic manipulations than that taking a real-valued weight continuously from $-1$ to $+1$. The former is called a signed particle and the latter a weighted particle. In this paper, we propose two efficient strategies to realize the signed-particle implementation. One is to interpret the multiplicative functional as the probability to generate pairs of particles instead of the incremental weight, and the other is to utilize a bootstrap filter to adjust the skewness of particle weights. Performance evaluations on the Gaussian barrier scattering (2D) and a Helium-like system (4D) demonstrate the feasibility of both strategies and the variance reduction property of the second approach. We provide an improvement of the first signed-particle implementation that partially alleviates the restriction on the time step and perform a thorough theoretical and numerical comparison among all the existing signed-particle implementations. Details on implementing the importance sampling according to the quasi-probability density and an efficient resampling or particle reduction are also provided.


title: Charting the Parameter Space of the 21-cm Power Spectrum

abstract: The high-redshift 21-cm signal of neutral hydrogen is expected to be observed within the next decade and will reveal epochs of cosmic evolution that have been previously inaccessible. Due to the lack of observations, many of the astrophysical processes that took place at early times are poorly constrained. In recent work we explored the astrophysical parameter space and the resulting large variety of possible global (sky-averaged) 21-cm signals. Here we extend our analysis to the fluctuations in the 21-cm signal, accounting for those introduced by density and velocity, Ly$\alpha$ radiation, X-ray heating, and ionization. While the radiation sources are usually highlighted, we find that in many cases the density fluctuations play a significant role at intermediate redshifts. Using both the power spectrum and its slope, we show that properties of high-redshift sources can be extracted from the observable features of the fluctuation pattern. For instance, the peak amplitude of ionization fluctuations can be used to estimate whether heating occurred early or late and, in the early case, to also deduce the cosmic mean ionized fraction at that time. The slope of the power spectrum has a more universal redshift evolution than the power spectrum itself and can thus be used more easily as a tracer of high-redshift astrophysics. Its peaks can be used, for example, to estimate the redshift of the Ly$\alpha$ coupling transition and the redshift of the heating transition (and the mean gas temperature at that time). We also show that a tight correlation is predicted between features of the power spectrum and of the global signal, potentially yielding important consistency checks.


title: Integrating Specialized Classifiers Based on Continuous Time Markov Chain

abstract: Specialized classifiers, namely those dedicated to a subset of classes, are often adopted in real-world recognition systems. However, integrating such classifiers is nontrivial. Existing methods, e.g. weighted average, usually implicitly assume that all constituents of an ensemble cover the same set of classes. Such methods can produce misleading predictions when used to combine specialized classifiers. This work explores a novel approach. Instead of combining predictions from individual classifiers directly, it first decomposes the predictions into sets of pairwise preferences, treating them as transition channels between classes, and thereon constructs a continuous-time Markov chain, and use the equilibrium distribution of this chain as the final prediction. This way allows us to form a coherent picture over all specialized predictions. On large public datasets, the proposed method obtains considerable improvement compared to mainstream ensemble methods, especially when the classifier coverage is highly unbalanced.


title: Sessile drop on oscillating incline

abstract: Natural or industrial flows of a fluid often involve droplets or bubbles of another fluid, pinned by physical or chemical impurities or by the roughness of the bounding walls. Here we study numerically one drop pinned on a circular hydrophilic patch, on an oscillating incline whose angle is proportional to $\sin(\omega\,t)$. The resulting deformation of the drop is measured by the displacement of its center of mass, which behaves like a driven over-damped linear oscillator with amplitude $A(\omega)$ and phase lag $\varphi(\omega)$. The phase lag is ${\cal O}(\omega)$ at small $\omega$ like a linear oscillator, but the amplitude is ${\cal O}(\omega^{-1})$ at large $\omega$ instead of ${\cal O}(\omega^{-2})$ for a linear oscillator. A heuristic explanation is given for this behaviour. The simulations were performed with the software Comsol in mode Laminar Two-Phase Flow, Level Set, with fluid 1 as engine oil and fluid 2 as water.


title: Beyond 16GB: Out-of-Core Stencil Computations

abstract: Stencil computations are a key class of applications, widely used in the scientific computing community, and a class that has particularly benefited from performance improvements on architectures with high memory bandwidth. Unfortunately, such architectures come with a limited amount of fast memory, which is limiting the size of the problems that can be efficiently solved. In this paper, we address this challenge by applying the well-known cache-blocking tiling technique to large scale stencil codes implemented using the OPS domain specific language, such as CloverLeaf 2D, CloverLeaf 3D, and OpenSBLI. We introduce a number of techniques and optimisations to help manage data resident in fast memory, and minimise data movement. Evaluating our work on Intel's Knights Landing Platform as well as NVIDIA P100 GPUs, we demonstrate that it is possible to solve 3 times larger problems than the on-chip memory size with at most 15\% loss in efficiency


title: Proceedings First Workshop on Formal Verification of Autonomous Vehicles

abstract: These are the proceedings of the workshop on Formal Verification of Autonomous Vehicles, held on September 19th, 2017 in Turin, Italy, as an affiliated workshop of the International Conference on integrated Formal Methods (iFM 2017). The workshop aim is to bring together researchers from the formal verification community that are developing formal methods for autonomous vehicles as well as researchers working, e.g., in the area of control theory or robotics, interested in applying verification techniques for designing and developing of autonomous vehicles.