title: RBF-PU Interpolation with Variable Subdomain Sizes and Shape Parameters

abstract: In this paper, we deal with the challenging computational issue of interpolating large data sets, with eventually non-homogeneous densities. To such scope, the Radial Basis Function Partition of Unity (RBF-PU) method has been proved to be a reliable numerical tool. However, there are not available techniques enabling us to efficiently select the sizes of the local PU subdomains which, together with the value of the RBF shape parameter, greatly influence the accuracy of the final fit. Thus here, by minimizing an \emph{a priori} error estimate, we propose a RBF-PU method by suitably selecting variable shape parameters and subdomain sizes. Numerical results and applications show performaces of the interpolation technique.


title: The multiple soliton and peakon solutions of the Dullin-Gottwald-Holm equation

abstract: Explicit multi-soliton and multi-peakon solutions of the Dullin-Gottwald-Holm equation are constructed via Darboux transformation and direct computation, respectively. To this end we first map the Dullin-Gottwald-Holm equation to a negative order KdV equation by a reciprocal transformation. Then we use the Darboux matrix approach to derive multi-soliton solutions of the Dullin-Gottwald-Holm equation from the solutions of the negative order KdV equation. Finally, we find multi-peakon solutions of the Dullin-Gottwald-Holm equation in weak sense. For $\alpha=\frac{1}{2}$ and $\alpha\neq \frac{1}{2}$, several types of two-peakon solutions are discussed in detail. Moreover, the dynamic behaviors of the obtained solutions are illustrated through some figures.


title: Measuring Coherence with Entanglement Concurrence

abstract: Quantum coherence is a fundamental manifestation of the quantum superposition principle. Recently, Baumgratz \emph{et al}. [Phys. Rev. Lett. \textbf{113}, 140401 (2014)] presented a rigorous framework to quantify coherence from the view of theory of physical resource. Here we propose a new valid quantum coherence measure which is a convex roof measure, for a quantum system of arbitrary dimension, essentially using the generalized Gell-Mann matrices. Rigorous proof shows that the proposed coherence measure fulfills all the requirements dictated by the resource theory of quantum coherence measures. Moreover, strong links between the resource frameworks of coherence concurrence and entanglement concurrence is derived, which shows that any degree of coherence with respect to some reference basis can be converted to entanglement via incoherent operations. Our work provides a clear quantitative and operational connection between coherence and entanglement based on two kinds of concurrence. This new coherence measure, coherence concurrence, may be beneficial to the study of quantum coherence.


title: Event reweighting with the NuWro neutrino interaction generator

abstract: Event reweighting has been implemented in the NuWro neutrino event generator for a number of free theory parameters in the interaction model. Event reweighting is a key analysis technique, used to efficiently study the effect of neutrino interaction model uncertainties. This opens up the possibility for NuWro to be used as a primary event generator by experimental analysis groups. A preliminary model tuning to ANL and BNL data of quasi-elastic and single pion production events was performed to validate the reweighting engine.


title: The effect of delay on contact tracing

abstract: We consider a model for an infectious disease in the onset of an outbreak. We introduce contact tracing incorporating a tracing delay. The effect of randomness in the delay and the effect of the length of this delay in comparison to the infectious period of the disease respectively to a latency period on the effect of tracing, given e.g. by the change of the reproduction number, is analyzed. We focus particularly on the effect of randomness in the tracing delay.


title: Computation of Kazhdan-Lusztig polynomials and some applications to finite groups

abstract: We discuss a practical algorithm to compute parabolic Kazhdan-Lusztig polynomials. As an application we compute Kazhdan-Lusztig polynomials which are needed to evaluate a character formula for reductive groups due to Lusztig. Some coefficients of these polynomials have interesting interpretations for certain finite groups. We find examples of finite dimensional modules for finite groups with much higher dimensional first cohomology group than in all previously known cases. Some of these examples lead to the construction of finite groups with many maximal subgroups, contradicting an old conjecture by G.~E.~Wall.


title: UY Puppis - A New Anomalous Z Cam Type Dwarf Nova

abstract: The defining characteristic of Z Cam stars are standstills in their light curves. Some Z Cams exhibit atypical behaviour by going into outburst from a standstill. It has previously been suggested that UY Pup had been a Z Cam star, but it was ruled out due to its long-term light curve. However, in December 2015 UY Pup went into outburst and unexpectedly entered into a short standstill instead of returning to quiescence. Furthermore, UY Pup exhibited additional unusual behaviour with two outbursts detected during its standstill. After this standstill UY Pup made a brief excursion to a quiescence state and slowly rose to a longer and well-defined standstill, where it again went into another outburst. Through comparative analysis, research, and observational data of UY Pup it is evident and thus concluded that it is indeed a Z Cam star, in which renders it to be one of only four known anomalous Z Cam stars.


title: The Chern-Simons invariants for the double of a compression body

abstract: Given a 3-manifold that can be written as the double of a compression body, we compute the Chern-Simons critical values for arbitrary compact connected structure groups. We also show that the moduli space of flat connections is connected when there are no reducibles.


title: A Hecke algebra attached to mod 2 modular forms of level 5

abstract: Let $F$ be the element $\sum_{n\ \mathit{odd},\ n>0}x^{n^{2}}$ of $Z/2[[x]]$. Set $G=F(x^{5})$, $D=F(x)+F(x^{25})$. For $k>0$, $(k,10)=1$, define $D_{k}$ as follows. $D_{1}=D$, $D_{3}=D^{8}/G$, $D_{7}=D^{2}G$, $D_{9}=D^{4}G$; furthermore $D_{k+10}=G^{2}D_{k}$. Using modular forms of level $\Gamma_{0}(5)$ we show that the space $W$ spanned by the $D_{k}$ is stabilized by the formal Hecke operators $T_{p}:Z/2[[x]]\rightarrow Z/2[[x]]$, $p\ne 2$ or $5$. And we determine the structure of the (completed) shallow Hecke algebra attached to $W$. This algebra proves to be a power series ring in $T_{3}$ and $T_{7}$ with an element of square $0$ adjoined. As Hecke module, $W$ identifies with a certain subquotient of the space of mod~2 modular forms of level $\Gamma_{0}(5)$, and our Hecke algebra result parallels findings in level 1 (by J.-L. Nicolas and J.-P. Serre) and in level $\Gamma_{0}(3)$ by us.


title: 30 years, some 700 integrals, and 1 dessert, or: Electroweak two-loop corrections to the Zbb vertex

abstract: The one-loop corrections to the weak mixing angle $\sin^2\theta_{eff}^b$ derived from the $Z{\bar b}b$ vertex, are known since 1985. It took another 30 years to calculate the complete electroweak two-loop corrections to $\sin^2\theta_{eff}^b$. The main obstacle was the calculation of the O(700) bosonic two-loop vertex integrals with up to three mass scales, at $s=M_Z^2$. We did not perform the usual integral reduction and master evaluation, but chose a completely numerical approach, using two different calculational chains. One method relies on publicly available sector decomposition implementations. Further, we derived Mellin-Barnes (MB) representations, exploring the publicly available MB suite. We had to supplement the MB suite by two new packages: AMBRE~3, a Mathematica program, for the efficient treatment of non-planar integrals and MBnumerics for advanced numerics in the Minkowskian space-time. Our preliminary result for LL2016, the "dessert", for the electroweak bosonic two-loop contributions to $\sin^2\theta_{eff}^b$ is: $\Delta \sin^2\theta_{eff}^{b(\alpha^2,\rm bos)} = \sin^2\theta_W ~ \Delta\kappa_b^{(\alpha^2,bos)}$, with $\Delta\kappa_b^{(\alpha^2,bos)} = -1.0276 x 10^{-4}$. This contribution is about a quarter of the corresponding fermionic corrections and of about the same magnitude as several of the known higher-order QCD corrections. The $\sin^2\theta_{eff}^b$ is now predicited in the Standard Model with a relative error of $10^{-4}$ [1].


title: On the nature and correction of the spurious S-wise spiral galaxy winding bias in Galaxy Zoo 1

abstract: The Galaxy Zoo 1 catalog displays a bias towards the S-wise winding direction in spiral galaxies which has yet to be explained. The lack of an explanation confounds our attempts to verify the Cosmological Principle, and has spurred some debate as to whether a bias exists in the real universe. The bias manifests not only in the obvious case of trying to decide if the universe as a whole has a winding bias, but also in the more insidious case of selecting which galaxies to include in a winding direction survey. While the former bias has been accounted for in a previous image-mirroring study, the latter has not. Furthermore, the bias has never been {\em corrected} in the GZ1 catalog, as only a small sample of the GZ1 catalog was re-examined during the mirror study. We show that the existing bias is a human {\em selection} effect rather than a human chirality bias. In effect, the excess S-wise votes are spuriously "stolen" from the elliptical and edge-on-disk categories, not the Z-wise category. Thus, when selecting a set of spiral galaxies by imposing a threshold $T$ so that $\max(P_S,P_Z) > T$ or $P_S+P_Z>T$, we spuriously select more S-wise than Z-wise galaxies. We show that when a provably unbiased machine selects which galaxies are spirals independent of their chirality, the S-wise surplus vanishes, even if humans are still used to determine the chirality. Thus, when viewed across the entire GZ1 sample (and by implication, the Sloan catalog), the winding direction of arms in spiral galaxies as viewed from Earth is consistent with the flip of a fair coin.


title: $C^3$-index: A PageRank based multi-faceted metric for authors' performance measurement

abstract: Ranking scientific authors is an important but challenging task, mostly due to the dynamic nature of the evolving scientific publications. The basic indicators of an author's productivity and impact are still the number of publications and the citation count (leading to the popular metrics such as h-index, g-index etc.). H-index and its popular variants are mostly effective in ranking highly-cited authors, thus fail to resolve ties while ranking medium-cited and low-cited authors who are majority in number. Therefore, these metrics are inefficient to predict the ability of promising young researchers at the beginning of their career. In this paper, we propose $C^3$-index that combines the effect of citations and collaborations of an author in a systematic way using a weighted multi-layered network to rank authors. We conduct our experiments on a massive publication dataset of Computer Science and show that - (i) $C^3$-index is consistent over time, which is one of the fundamental characteristics of a ranking metric, (ii) $C^3$-index is as efficient as h-index and its variants to rank highly-cited authors, (iii) $C^3$-index can act as a conflict resolution metric to break ties in the ranking of medium-cited and low-cited authors, (iv) $C^3$-index can also be used to predict future achievers at the early stage of their career.


title: Recovering lost 21 cm radial modes via cosmic tidal reconstruction

abstract: 21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe, at redshifts $z$ from 1 to 10. Unfortunately, many of the key cross-correlations with the cosmic microwave background and photo-$z$ galaxies have been thought to be impossible due to the foreground contamination for radial modes with small wave numbers. In this paper, we apply tidal reconstruction to the simulated 21 cm fields and recover the lost large-scale radial modes successfully. We estimate the detectability of the cross-correlation signals and find they can be detected at high significance with current 21 cm experiments. The tidal field reconstruction method opens up a new set of possibilities to probe the Universe and is extremely valuable not only for 21 cm surveys but also for cosmic microwave background and photometric-redshift observations.


title: Extended skyrmion lattice scattering and long-time memory in the chiral magnet Fe$_{1-x}$Co$_x$Si

abstract: Small angle neutron scattering measurements on a bulk single crystal of the doped chiral magnet Fe$_{1-x}$Co$_x$Si with $x$=0.3 reveal a pronounced effect of the magnetic history and cooling rates on the magnetic phase diagram. The extracted phase diagrams are qualitatively different for zero and field cooling and reveal a metastable skyrmion lattice phase outside the A-phase for the latter case. These thermodynamically metastable skyrmion lattice correlations coexist with the conical phase and can be enhanced by increasing the cooling rate. They appear in a wide region of the phase diagram at temperatures below the $A$-phase but also at fields considerably smaller or higher than the fields required to stabilize the A-phase.


title: Hierarchical DSSC structures based on single walled TiO2 nanotube arrays reach back-side illumination solar light conversion efficiency of 8%

abstract: In the present work we introduce a path to the controlled construction of DSSCs based on hierarchically structured single walled, self-organized TiO2 layers. In a first step we describe a simple approach to selectively remove the inner detrimental shell of anodic TiO2 nanotubes (NTs). This then allows controlled well-defined layer-by-layer decoration of these TiO2-NT walls with TiO2 nanoparticles (this in contrast to conventional TiO2 nanotubes). We show that such defined multiple layered decoration can be optimized to build dye sensitized solar cells that (under back-side illumination conditions) can yield solar light conversion efficiencies in the range of 8 %. The beneficial effects observed can be ascribed to a combination of three factors : 1) improved electronic properties of the single walled tubes themselves, 2) a further improvement of the electronic properties by the defined TiCl4 treatment, and 3) a higher specific dye loading that becomes possible for the layer-by-layer decorated single walled tubes.


title: On derivatives of Siegel-Eisenstein series over global function fields

abstract: The aim of this article is to study the derivative of "incoherent" Siegel-Eisenstein series on symplectic groups over function fields. By the Siegel-Weil formula for "coherent" Siegel-Eisenstein series, we can relate the non-singular Fourier coefficients of the derivative in question to the arithmetic of quadratic forms. Restricting to the special case when the incoherent quadratic space has dimension 2, we explicitly compute all the Fourier coefficients, and connect the derivative with the special cycles on the coarse moduli schemes of rank 2 Drinfeld modules with "complex multiplication."


title: DSValidator: An Automated Counterexample Reproducibility Tool for Digital Systems (Tool Demonstration)

abstract: An automated counterexample reproducibility tool based on MATLAB is presented, called DSValidator, with the goal of reproducing counterexamples that refute specific properties related to digital systems. We exploit counterexamples generated by the Digital System Verifier (DSVerifier), which is a model checking tool based on satisfiability modulo theories for digital systems. DSValidator reproduces the execution of a digital system, relating its input with the counterexample, in order to establish trust in a verification result. We show that DSValidator can validate a set of intricate counterexamples for digital controllers used in a real quadrotor attitude system within seconds and also expose incorrect verification results in DSVerifier. The resulting toolbox leverages the potential of combining different verification tools for validating digital systems via an exchangeable counterexample format.


title: The planar $L_p$-Minkowski problem for $0< p<1$

abstract: Necessary and sufficient conditions for the existence of solutions to the asymmetric $L_p$ Minkowski problem in $\mathbb{R}^2$ are established for $0 < p < 1$.


title: Nonlinear Sturm Oscillation: from the interval to a star

abstract: The Sturm oscillation property, i.e. that the $n$-th eigenfunction of a Sturm-Liouville operator on an interval has $n -1$ zeros (nodes), has been well studied. This result is known to hold when the interval is replaced by a metric (quantum) tree graph. We prove that the solutions of the real stationary nonlinear Schr\"odinger equation on an interval satisfy a nonlinear version of the Sturm oscillation property. However, we show that unlike the linear theory, the nonlinear version of the Sturm oscillation breaks down already for a star graph. We point out conditions under which this violation can be assured.


title: Hawking radiation of scalar and vector particles from 5D Myers-Perry black holes

abstract: In the present paper we explore the Hawking radiation as a quantum tunneling effect from a rotating 5 dimensional Myers-Perry black hole (5D-MPBH) with two independent angular momentum components. First, we investigate the Hawking temperature by considering the tunneling of massive scalar particles and spin-1 vector particles from the 5D-MPBH in the Painlev\'{e} coordinates and then in the corotating frames. More specifically, we solve the Klein-Gordon and Proca equations by applying the WKB method and Hamilton-Jacobi equation in both cases. Finally, we recover the Hawking temperature and show that coordinates systems do not affect the Hawking temperature.


title: Creep and fluidization in thermal amorphous solids

abstract: When submitted to a constant mechanical load, many amorphous solids display power law creep followed by fluidization. A fundamental understanding of these processes is still far from being achieved. Here, we characterize creep and fluidization on the basis of a mesoscopic viscoplastic model that includes thermally activated yielding events and a broad distribution of energy barriers, which may be lowered under the effect of a local deformation. We relate the creep exponent observed before fluidization to the width of barrier distribution and to the specific form of stress redistribution following yielding events. We show that Andrade creep is accompanied by local strain-hardening driven by stress redistribution and find that the fluidization depends exponentially on the applied stress. The simulation results are interpreted in the light of a mean-field analysis, and should help in rationalizing the creep phenomenology of amorphous solids.


title: Fabrication and Analysis of Three-Layer All-Silicon Interference Optical Filter with Sub-Wavelength Structure toward High Performance Terahertz Optics

abstract: We propose an all-silicon multi-layer interference filter composed solely of silicon with sub-wavelength structure (SWS) in order to realize high performance optical filters operating in the THz frequency region with robustness against cryogenic thermal cycling and mechanical damage. We demonstrate fabrication of a three-layer prototype using well-established common micro-electro-mechanical systems (MEMS) technologies as a first step toward developing practical filters. The measured transmittance of the three-layer filter agrees well with the theoretical transmittances calculated by a simple thin-film calculation with effective refractive indices as well as a rigorous coupled-wave analysis simulation. We experimentally show that SWS layers can work as homogeneous thin-film interference layers with effective refractive indices even if there are multiple SWS layers in a filter.


title: Complexity Issues in Computing Spectra, Pseudospectra and Resolvents

abstract: We display methods that allow for computations of spectra, pseudospectra and resolvents of linear operators on Hilbert spaces and also elements in unital Banach algebras. The paper considers two different approaches, namely, pseudospectral techniques and polynomial numerical hull theory. The former is used for Hilbert space operators whereas the latter can handle the general case of elements in a Banach algebra. This approach leads to multicentric holomorphic calculus. We also discuss some new types of pseudospectra and the recently defined Solvability Complexity Index


title: Imbalanced magnetohydrodynamic turbulence modified by velocity shear in the solar wind

abstract: We study incompressible imbalanced magnetohydrodynamic turbulence in the presence of background velocity shears. Using scaling arguments, we show that the turbulent cascade is significantly accelerated when the background velocity shear is stronger than the velocity shears in the subdominant Alfv% \'{e}n waves at the injection scale. The spectral transport is then controlled by the background shear rather than the turbulent shears and the Tchen spectrum with spectral index $-1$ is formed. This spectrum extends from the injection scale to the scale of the spectral break where the subdominant wave shear becomes equal to the background shear. The estimated spectral breaks and power spectra are in good agreement with those observed in the fast solar wind. The proposed mechanism can contribute to enhanced turbulent cascades and modified $-1$ spectra observed in the fast solar wind with strong velocity shears. This mechanism can also operate in many other astrophysical environments where turbulence develops on top of non-uniform plasma flows.


title: Numerical and analytical results for the two disks capacitor problem

abstract: In this paper we study the two disks capacitor, for equal and different radii. The new results obtained allow a complete characterization of capacity coefficients and forces at short distances. An extensive numerical calculation confirms the theoretical results. The study shows the existence of a hierarchy in the divergent behavior of the capacitance coefficients and this implies some unusual behavior of the forces, strictly related to the dimensionality of the near-contact zone between electrodes.


title: Equivalence of Weighted Anchored and ANOVA Spaces of Functions with Mixed Smoothness of Order one in $L_p$

abstract: We consider $\gamma$-weighted anchored and ANOVA spaces of functions with mixed first order partial derivatives bounded in a weighted $L_p$ norm with $1 \leq p \leq \infty$. The domain of the functions is $D^d$, where $D \subseteq \mathbb{R}$ is a bounded or unbounded interval. We provide conditions on the weights $\gamma$ that guarantee that anchored and ANOVA spaces are equal (as sets of functions) and have equivalent norms with equivalence constants uniformly or polynomially bounded in $d$. Moreover, we discuss applications of these results to integration and approximation of functions on $D^d$.


title: On one-dimensional compressible Navier-Stokes equations for a reacting mixture in unbounded domains

abstract: In this paper we consider the one-dimensional Navier-Stokes system for a heat-conducting, compressible reacting mixture which describes the dynamic combustion of fluids of mixed kinds on unbounded domains. This model has been discussed on bounded domains by Chen (SIAM Jour. Math. Anal., 23 (1992), 609--634) and Chen-Hoff-Trivisa (Arch. Rat. Mech. Anal. 166 (2003), 321--358) among others, in which the reaction rate function is a discontinuous function obeying the Arrhenius Law. We prove the global existence of weak solutions to this model on one-dimensional unbounded domains with large initial data in $H^1$. Moreover, the large-time behaviour of the weak solution is identified and proved. In particular, the uniform-in-time bounds for the temperature and specific volume have been established via energy estimates. For this purpose we utilise techniques developed by Kazhikhov and coauthors ({\it cf.} Siber. Math. Jour. 23 (1982), 44--49; Jour. Appl. Math. Mech., 41 (1977), 273--282), as well as a crucial estimate in the recent work by Li-Liang (Arch. Rat. Mech. Anal. 220 (2016), 1195--1208). Several new estimates are also established, in order to treat the unbounded domain and the reacting terms.


title: The cycle structure of a Markoff automorphism over finite fields

abstract: We begin an investigation of the action of pseudo-Anosov elements of $\mathrm{Out}(\mathbf{F}_{2})$ on the Markoff-type varieties \[ \mathbb{X}_{\kappa}:\:x^{2}+y^{2}+z^{2}=xyz+2+\kappa \] over finite fields $\mathbb{F}_{p}$ with $p$ prime. We first make a precise conjecture about the permutation group generated by $\mathrm{Out}(\mathbf{F}_{2})$ on $\mathbb{X}_{-2}(\mathbb{F}_{p})$ that shows there is no obstruction at the level of the permutation group to a pseudo-Anosov acting `generically'. We prove that this conjecture is sharp. We show that for a fixed pseudo-Anosov $g\in\mathrm{Out}(\mathbf{F}_{2})$, there is always an orbit of $g$ of length $\geq C\log p+O(1)$ on $\mathbb{X}_{\kappa}(\mathbb{F}_{p})$ where $C>0$ is given in terms of the eigenvalues of $g$ viewed as an element of $\mathrm{GL}_{2}(\mathbf{Z})$. This improves on a result of Silverman (2007) that applies to general morphisms of quasi-projective varieties. We have discovered that the asymptotic $(p\to\infty)$ behavior of the longest orbit of a fixed pseudo-Anosov $g$ acting on $\mathbb{X}_{-2}(\mathbb{F}_{p})$ is dictated by a dichotomy that we describe both in combinatorial terms and in algebraic terms related to Gauss's ambiguous binary quadratic forms, following Sarnak. This dichotomy is illustrated with numerics, based on which we formulate a precise conjecture.


title: Polar Coordinates and Noncommutative Phase Space

abstract: The so-called Weyl transform is a linear map from a commutative algebra of functions to a noncommutative algebra of linear operators, characterized by an action on Cartesian coordinate functions of the form $(x, y) \mapsto (X, Y)$ such that $XY -YX = i\epsilon I$, i.e. the defining relation for the Heisenberg Lie algebra. Study of this transform has been expansive. We summarize many important results from the literature. The primary goal of this work is to prove the final result: the realization of the polar transformation $(\rho, e^{i\theta}) \mapsto (R, e^{i\Theta})$ in terms of explicit orthogonal function expansions, while starting from elementary principles and utilizing minimal machinery. Our results are not strictly original but their presentation here is intended to simplify introduction to these subjects in a novel manner.


title: On the topology of the Lorenz system

abstract: We present a new paradigm for three dimensional chaos, and specifically for the Lorenz equations. The main difficulty in these equations and for a generic flow in dimension three is the existence of singularities. We show how to use knot theory as a way to remove the singularities. Specifically, we claim: (1) For certain parameters, the Lorenz system has an invariant one dimensional curve, which is a trefoil knot. The knot is a union of invariant manifolds of the singular points. (2) The flow is topologically equivalent to an Anosov flow on the complement of this curve, and even to a geodesic flow. (3) When varying the parameters, the system exhibits topological phase transitions, i.e. for special parameter values, it will be topologically equivalent to an Anosov flow on a knot complement, and different knots appear for different parameter values. The steps of a mathematical proof of these statements are at different stages. Some have been proven, for some we present numerical evidence and some are still conjectural.


title: An Alternating Automaton for First-Order Linear Temporal Logic--Tech Report

abstract: In this paper we give automata-based representation of LTL-FO$^+$ properties. LTL-FO$^+$ is an extension of LTL that includes first-order quantification over bounded variable, thus greatly increasing the expressivity of the language. An automata representation of this formalism allows greater ease in writing and understanding properties, as well as in performing manipulations, such as negation or emptiness checking. The automata representation of an LTL-FO$^+$ formula has finite size regardless of the domain of quantified variables, and the number of states that is linear in the size of the property.


title: The limits of weak selection and large population size in evolutionary game theory

abstract: Evolutionary game theory is a mathematical approach to studying how social behaviors evolve. In many recent works, evolutionary competition between strategies is modeled as a stochastic process in a finite population. In this context, two limits are both mathematically convenient and biologically relevant: weak selection and large population size. These limits can be combined in different ways, leading to potentially different results. We consider two orderings: the wN limit, in which weak selection is applied before the large population limit, and the Nw limit, in which the order is reversed. Formal mathematical definitions of the Nw and wN limits are provided. Applying these definitions to the Moran process of evolutionary game theory, we obtain asymptotic expressions for fixation probability and conditions for success in these limits. We find that the asymptotic expressions for fixation probability, and the conditions for a strategy to be favored over a neutral mutation, are different in the Nw and wN limits. However, the ordering of limits does not affect the conditions for one strategy to be favored over another.


title: Structure of parabolically induced modules for affine Kac-Moody algebras

abstract: The main result of the paper establishes the irreducibility of a large family of nonzero central charge induced modules over Affine Lie algebras for any non standard parabolic subalgebra. It generalizes all previously known partial results and provides a a construction of many new irreducible modules.


title: Ab initio prediction of stable nanotwin double layers and 4O structure in Ni$_{2}$MnGa

abstract: The ab initio electronic structure calculations of the Ni$_{2}$MnGa The alloy indicate that the orthorhombic 4O structure exhibits the lowest energy compared to all known martensitic structures. The 4O structure is formed by nanotwin double layers, i.e., oppositely oriented nanotwins consisting of two (101) lattice planes of nonmodulated martensitic structure. It exhibits the lowest occupation of density of states at the Fermi level. The total energy 1.98 meV/atom below the energy of nonmodulated martensite is achieved within structural relaxation by shifting Mn and Ga atoms at the nanotwin boundaries. The same atomic shift can also be found in other martensitic nanotwinned or modulated structures such as 10M and 14M, which indicates the importance of the nanotwin double layer for the stability of these structures. Our discovery shows that the nanotwinning or modulation is a natural property of low-temperature martensitic phases in Ni-Mn-Ga alloys.


title: Fluctuations of Functions of Wigner Matrices

abstract: We show that matrix elements of functions of $N\times N$ Wigner matrices fluctuate on a scale of order $N^{-1/2}$ and we identify the limiting fluctuation. Our result holds for any function $f$ of the matrix that has bounded variation and thus considerably relaxes the regularity requirement imposed in [7,11].


title: Radiation of non-relativistic particle on a conducting sphere and a string of spheres

abstract: The radiation arising under uniform motion of non-relativistic charged particle by (or through) perfectly conducting sphere is considered. The rigorous results are obtained using the method of images known from electrostatics.


title: Deep image mining for diabetic retinopathy screening

abstract: Deep learning is quickly becoming the leading methodology for medical image analysis. Given a large medical archive, where each image is associated with a diagnosis, efficient pathology detectors or classifiers can be trained with virtually no expert knowledge about the target pathologies. However, deep learning algorithms, including the popular ConvNets, are black boxes: little is known about the local patterns analyzed by ConvNets to make a decision at the image level. A solution is proposed in this paper to create heatmaps showing which pixels in images play a role in the image-level predictions. In other words, a ConvNet trained for image-level classification can be used to detect lesions as well. A generalization of the backpropagation method is proposed in order to train ConvNets that produce high-quality heatmaps. The proposed solution is applied to diabetic retinopathy (DR) screening in a dataset of almost 90,000 fundus photographs from the 2015 Kaggle Diabetic Retinopathy competition and a private dataset of almost 110,000 photographs (e-ophtha). For the task of detecting referable DR, very good detection performance was achieved: $A_z = 0.954$ in Kaggle's dataset and $A_z = 0.949$ in e-ophtha. Performance was also evaluated at the image level and at the lesion level in the DiaretDB1 dataset, where four types of lesions are manually segmented: microaneurysms, hemorrhages, exudates and cotton-wool spots. The proposed detector outperforms recent algorithms trained to detect those lesions specifically, as well as competing heatmap generation algorithms for ConvNets. This detector is part of the Messidor system for mobile eye pathology screening. Because it does not rely on expert knowledge or manual segmentation for detecting relevant patterns, the proposed solution is a promising image mining tool, which has the potential to discover new biomarkers in images.


title: Higher Commutator Theory for Congruence Modular Varieties

abstract: We develop the basic properties of the higher commutator for congruence modular varieties.


title: Differential-free characterisation of smooth mappings with controlled growth

abstract: In this paper we give some generalizations and improvements of the Pavlovi\'c result on the Holland-Walsh type characterization of the Bloch space of continuously differentiable (smooth) functions in the unit ball in $\mathbf{R}^m$.


title: Reinforcement Learning in Conflicting Environments for Autonomous Vehicles

abstract: In this work, we investigate the application of Reinforcement Learning to two well known decision dilemmas, namely Newcomb's Problem and Prisoner's Dilemma. These problems are exemplary for dilemmas that autonomous agents are faced with when interacting with humans. Furthermore, we argue that a Newcomb-like formulation is more adequate in the human-machine interaction case and demonstrate empirically that the unmodified Reinforcement Learning algorithms end up with the well known maximum expected utility solution.


title: STEPS: Predicting place attributes via spatio-temporal analysis

abstract: In recent years, a vast amount of research has been conducted on learning people's interests from their actions. Yet their collective actions also allow us to learn something about the world, in particular, infer attributes of places people visit or interact with. Imagine classifying whether a hotel has a gym or a swimming pool, or whether a restaurant has a romantic atmosphere without ever asking its patrons. Algorithms we present can do just that. Many web applications rely on knowing attributes of places, for instance, whether a particular restaurant has WiFi or offers outdoor seating. Such data can be used to support a range of user experiences, from explicit query-driven search to personalized place recommendations. However, obtaining these attributes is generally difficult, with existing approaches relying on crowdsourcing or parsing online reviews, both of which are noisy, biased, and have limited coverage. Here we present a novel approach to classifying place attributes, which learns from patrons' visit patterns based on anonymous observational data. Our method, STEPS, learns from aggregated sequences of place visits. For example, if many people visit the restaurant on a Saturday evening, coming from a luxury hotel or theater, and stay for a long time, then this restaurant is more likely to have a romantic atmosphere. On the other hand, if most people visit the restaurant on weekdays, coming from work or a grocery store, then the restaurant is less likely to be romantic. We show that such transition features are highly predictive of place attributes. In an extensive empirical evaluation, STEPS nearly doubled the coverage of a state of the art approach thanks to learning from observational location data, which allowed our method to reason about many more places.


title: Automatic Identification of Sarcasm Target: An Introductory Approach

abstract: Past work in computational sarcasm deals primarily with sarcasm detection. In this paper, we introduce a novel, related problem: sarcasm target identification i.e., extracting the target of ridicule in a sarcastic sentence). We present an introductory approach for sarcasm target identification. Our approach employs two types of extractors: one based on rules, and another consisting of a statistical classifier. To compare our approach, we use two baselines: a na\"ive baseline and another baseline based on work in sentiment target identification. We perform our experiments on book snippets and tweets, and show that our hybrid approach performs better than the two baselines and also, in comparison with using the two extractors individually. Our introductory approach establishes the viability of sarcasm target identification, and will serve as a baseline for future work.


title: Bounds in Cohen's idempotent theorem

abstract: We show that if $G$ is a finite Abelian group and $f$ is an integer-valued map on $G$ with algebra norm at most $M$ then there is some $L < \exp(M^{4+o(1)})$, cosets of (possibly different) subgroups $W_1,...,W_L$, and $s_1,...,s_L \in \{-1,1\}$ such that $f=\sum_i{s_i1_{W_i}}$.


title: Equivalence between contextuality and negativity of the Wigner function for qudits

abstract: Contextuality and negativity of the Wigner function are two notions of non-classicality for quantum systems. Howard, Wallman, Veitch and Emerson proved recently that these two notions coincide for qudits in odd prime dimension. This equivalence is particularly important since it promotes contextuality as a ressource that magic states must possess in order to allow for a quantum speed-up. We propose a simple proof of the equivalence between contextuality and negativity of the Wigner function based on character theory. This simplified approach allows us to generalize this equivalence to multiple qudits and to any qudit system of odd local dimension.


title: Inferring parameters of prey switching in a plankton 1 predator--2 prey system with a linear preference tradeoff

abstract: We construct two ordinary-differential-equation models of a predator feeding adaptively on two prey types, and we evaluate the models' ability to fit data on freshwater plankton. We model the predator's switch from one prey to the other in two different ways: (1) smooth switching using a hyperbolic tangent function; and (2) by incorporating a parameter that changes abruptly across the switching boundary as a system variable that is coupled to the population dynamics. We conduct linear stability analyses, use approximate Bayesian computation (ABC) combined with a population Monte Carlo (PMC) method to fit model parameters, and compare model predictions quantitatively to data for ciliate predators and their two algal prey groups collected from Lake Constance on the German--Swiss--Austrian border. We show that the two models fit the data well when the smooth transition is steep, supporting the simplifying assumption of a discontinuous prey switching behavior for this scenario. We thus conclude that prey switching is a possible mechanistic explanation for the observed ciliate--algae dynamics in Lake Constance in spring, but that these data cannot distinguish between the details of prey switching that are encoded in these different models.


title: Observation of Upsilon(1S) pair production in proton-proton collisions at sqrt(s) = 8 TeV

abstract: Pair production of Upsilon(1S) mesons is observed at the LHC in proton-proton collisions at sqrt(s) = 8 TeV by the CMS experiment in a data sample corresponding to an integrated luminosity of 20.7 inverse femtobarns. Both Upsilon(1S) candidates are fully reconstructed via their decays to mu+ mu-. The fiducial acceptance region is defined by an absolute Upsilon(1S) rapidity smaller than 2.0. The fiducial cross section for the production of Upsilon(1S) pairs, assuming that both mesons decay isotropically, is measured to be 68.8 +/- 12.7 (stat) +/- 7.4 (syst) +/- 2.8 (B) pb, where the third uncertainty comes from the uncertainty in the branching fraction of Upsilon(1S) decays to mu+ mu-. Assuming instead that the Upsilon(1S) mesons are produced with different polarizations leads to variations in the measured cross section in the range from -38% to +36%.


title: A statistical approach to covering lemmas

abstract: We discuss a statistical variant of Ruzsa's covering lemma and use it to show that if G is an Abelian group of bounded exponent and A in G has |A+A| < K|A| then the subgroup generated by A has size at most exp(O(K log^22K))|A|, where the constant in the big-O depends on the exponent of the group only.


title: Analysis of minimal representations of SL(n,R)

abstract: Some minimal representations of SL(n,R) can be realized on a Hilbert space of holomorphic functions. This is the analogue of the Brylinski-Kostant model. They can also be realized on a Hilbert space of homogeneous functions on ${\bboard R}^n$. This is the analogue of the Kobayashi-Orsted model. We will describe the two realizations and a transformation which maps one model to the other. It can be seen as an analogue of the classical Bargmann transform.


title: A boundary integral method for the general conjugation problem in multiply connected circle domains

abstract: We present a boundary integral method for solving a certain class of Riemann-Hilbert problems known as the general conjugation problem. The method is based on a uniquely solvable boundary integral equation with the generalized Neumann kernel. We present also an alternative proof for the existence and uniqueness of the solution of the general conjugation problem.


title: Excess area dependent scaling behavior of nano-sized membrane tethers

abstract: Thermal fluctuations in cell membranes manifest as an excess area (${\cal A}_{\rm ex}$) which governs a multitude of physical process at the sub-micron scale. We present a theoretical framework, based on an in silico tether pulling method, which may be used to reliably estimate ${\cal A}_{\rm ex}$ in live cells. The tether forces estimated from our simulations compare well with our experimental measurements for tethers extracted from ruptured GUVs and HeLa cells. We demonstrate the significance and validity of our method by showing that all our calculations along with experiments of tether extraction in 15 different cell types collapse onto two unified scaling relationships mapping tether force, tether radius, bending stiffness $\kappa$, and membrane tension $\sigma$. We show that $R_{\rm bead}$, the size of the wetting region, is an important determinant of the radius of the extracted tether, which is equal to $\xi=\sqrt{\kappa/2\sigma}$ (a characteristic length scale of the membrane) for $R_{\rm bead}<\xi$, and is equal to $R_{\rm bead}$ for $R_{\rm bead}>\xi$. We also find that the estimated excess area follows a linear scaling behavior that only depends on the true value of ${\cal A}_{\rm ex}$ for the membrane, based on which we propose a self-consistent technique to estimate the range of excess membrane areas in a cell.

categories: physics.bio-ph cond-mat.mes-hall cond-mat.soft cond-mat.stat-mech physics.comp-ph


title: Local Maxima and Improved Exact Algorithm for MAX-2-SAT

abstract: Given a MAX-2-SAT instance, we define a local maximum to be an assignment such that changing any single variable reduces the number of satisfied clauses. We consider the question of the number of local maxima that an instance of MAX-2-SAT can have. We give upper bounds in both the sparse and nonsparse case, where the sparse case means that there is a bound $d$ on the average number of clauses involving any given variable. The bounds in the nonsparse case are tight up to polylogarithmic factors, while in the sparse case the bounds are tight up to a multiplicative factor in $d$ for large $d$. Additionally, we generalize to the question of assignments which are maxima up to changing $k> 1$ variables simultaneously; in this case, we give explicit constructions with large (in a sense explained below) numbers of such maxima in the sparse case. The basic idea of the upper bound proof is to consider a random assignment to some subset of the variables and determine the probability that some fraction of the remaining variables can be fixed without considering interactions between them. The bounded results hold in the case of weighted MAX-2-SAT as well. Using this technique and combining with ideas from Ref. 6, we find an algorithm for weighted MAX-2-SAT which is faster for large $d$ than previous algorithms which use polynomial space; this algorithm does require an additional bounds on maximum weights and degree.


title: General Central Limit Theorems for Associated Sequences

abstract: In this paper, we provide general central limit theorems (CLT's) for associated random variables (rv's) following the approaches used by Newman (1980) and Olivera et al.(2012). Given some assumptions, a Lyapounov-Feller-Levy type theorem is stated. We next specify different particular CLT versions of associated sequences based on moment conditions. A comparison study with available CTL's is performed. As a by-product, we complete an important available theorem where an assumption was missing.


title: Depinning asymptotics in ergodic media

abstract: We study speeds of fronts in bistable, spatially inhomogeneous media at parameter regimes where speeds approach zero. We provide a set of conceptual assumptions under which we can prove power-law asymptotics for the speed, with exponent depending a local dimension of the ergodic measure near extremal values. We also show that our conceptual assumptions are satisfied in a context of weak inhomogeneity of the medium and almost balanced kinetics, and compare asymptotics with numerical simulations.


title: Surface tension and a self-consistent theory of soft composite solids with elastic inclusions

abstract: The importance of surface tension effects is being recognized in the context of soft composite solids, where it is found to significantly affect the mechanical properties, such as the elastic response to an external stress. It has recently been discovered that Eshelby's inclusion theory breaks down when the inclusion size approaches the elastocapillary length $L\equiv\gamma/E$, where $\gamma$ is the inclusion/host surface tension and $E$ is the host Young's modulus. Extending our recent results for liquid inclusions, here we model the elastic behavior of a non-dilute distribution of isotropic elastic spherical inclusions in a soft isotropic elastic matrix, subject to a prescribed infinitesimal far-field loading. Within our framework, the composite stiffness is uniquely determined by the elastocapillary length $L$, the spherical inclusion radius $R$, and the stiffness contrast parameter ${C}$, which is the ratio of the inclusion to the matrix stiffness. We compare the results with those from the case of liquid inclusions, and we derive an analytical expression for elastic cloaking of the composite by the inclusions. Remarkably, we find that the composite stiffness is influenced significantly by surface tension even for inclusions two orders of magnitude more stiff than the host matrix. Finally, we show how to simultaneously determine the surface tension and the inclusion stiffness using two independent constraints provided by global and local measurements.


title: Independent Component Analysis by Entropy Maximization with Kernels

abstract: Independent component analysis (ICA) is the most popular method for blind source separation (BSS) with a diverse set of applications, such as biomedical signal processing, video and image analysis, and communications. Maximum likelihood (ML), an optimal theoretical framework for ICA, requires knowledge of the true underlying probability density function (PDF) of the latent sources, which, in many applications, is unknown. ICA algorithms cast in the ML framework often deviate from its theoretical optimality properties due to poor estimation of the source PDF. Therefore, accurate estimation of source PDFs is critical in order to avoid model mismatch and poor ICA performance. In this paper, we propose a new and efficient ICA algorithm based on entropy maximization with kernels, (ICA-EMK), which uses both global and local measuring functions as constraints to dynamically estimate the PDF of the sources with reasonable complexity. In addition, the new algorithm performs optimization with respect to each of the cost function gradient directions separately, enabling parallel implementations on multi-core computers. We demonstrate the superior performance of ICA-EMK over competing ICA algorithms using simulated as well as real-world data.


title: Broken sublattice symmetry states in Bernal stacked multilayer graphene

abstract: We analyze the ordered phases of Bernal stacked multilayer graphene in the presence of interaction induced band gaps due to sublattice symmetry breaking potentials, whose solutions can be analyzed in terms of light-mass and heavy-mass pseudospin doublets which have the same Chern numbers but opposite charge polarization directions. The application of a perpendicular external electric field reveals an effective Hund's rule for the ordering of the sublattice pseudospin doublets in a tetralayer, while a similar but more complex phase diagram develops with increasing layer number.


title: Ultrafast pulse generation in a mode-locked Erbium chip waveguide laser

abstract: We report mode-locked ~1550 nm output of transform-limited ~180 fs pulses from a large mode-area (diameter ~ 50 {\mu}m) guided-wave erbium fluorozirconate glass laser. The passively mode-locked oscillator generates pulses with 25 nm bandwidth at 156 MHz repetition rate and peak-power of 260 W. Scalability to higher repetition rate is demonstrated by transform-limited 410 fs pulse output at 1.3 GHz. To understand the origins of the broad spectral output, the laser cavity is simulated by using a numerical solution to the Ginzburg-Landau equation. This paper reports the widest bandwidth and shortest pulses achieved from an ultra-fast laser inscribed waveguide laser.


title: Efficient quantum circuits for continuous-time quantum walks on composite graphs

abstract: In this paper, we investigate the simulation of continuous-time quantum walks on specific classes of graphs, for which it is possible to fast-forward the time-evolution operator to achieve constant-time simulation complexity and to perform the simulation exactly, while maintaining $\mbox{poly}(\mbox{log}(n))$ efficiency. In particular, we discuss two classes of composite graphs, commuting graphs and Cartesian product of graphs, that contain classes of graphs which can be simulated in this fashion. This allows us to identify new families of graphs that we can efficiently simulate in a quantum circuit framework, providing practical and explicit means to explore quantum-walk based algorithms in laboratories.


title: Fast and Reliable Parameter Estimation from Nonlinear Observations

abstract: In this paper we study the problem of recovering a structured but unknown parameter ${\bf{\theta}}^*$ from $n$ nonlinear observations of the form $y_i=f(\langle {\bf{x}}_i,{\bf{\theta}}^*\rangle)$ for $i=1,2,\ldots,n$. We develop a framework for characterizing time-data tradeoffs for a variety of parameter estimation algorithms when the nonlinear function $f$ is unknown. This framework includes many popular heuristics such as projected/proximal gradient descent and stochastic schemes. For example, we show that a projected gradient descent scheme converges at a linear rate to a reliable solution with a near minimal number of samples. We provide a sharp characterization of the convergence rate of such algorithms as a function of sample size, amount of a-prior knowledge available about the parameter and a measure of the nonlinearity of the function $f$. These results provide a precise understanding of the various tradeoffs involved between statistical and computational resources as well as a-prior side information available for such nonlinear parameter estimation problems.


title: Symplectic Geometric Algorithm for Quaternion Kinematical Differential Equation

abstract: Solving quaternion kinematical differential equations is one of the most significant problems in the automation, navigation, aerospace and aeronautics literatures. Most existing approaches for this problem neither preserve the norm of quaternions nor avoid errors accumulated in the sense of long term time. We present symplectic geometric algorithms to deal with the quaternion kinematical differential equation by modeling its time-invariant and time-varying versions with Hamiltonian systems by adopting a three-step strategy. Firstly, a generalized Euler's formula for the autonomous quaternion kinematical differential equation are proved and used to construct symplectic single-step transition operators via the centered implicit Euler scheme for autonomous Hamiltonian system. Secondly, the symplecitiy, orthogonality and invertibility of the symplectic transition operators are proved rigorously. Finally, the main results obtained are generalized to design symplectic geometric algorithm for the time-varying quaternion kinematical differential equation which is a non-autonomous and nonlinear Hamiltonian system essentially. Our novel algorithms have simple algorithmic structures and low time complexity of computation, which are easy to be implemented with real-time techniques. The correctness and efficiencies of the proposed algorithms are verified and validated via numerical simulations.


title: Tsallis thermostatics as a statistical physics of random chains

abstract: In this paper we point out that the generalized statistics of Tsallis-Havrda-Charv\'at can be conveniently used as a conceptual framework for statistical treatment of random chains. In particular, we use the path-integral approach to show that the ensuing partition function can be identified with the partition function of a fluctuating oriented random loop of arbitrary length and shape in a background scalar potential. To put some meat on the bare bones, we illustrate this with two statistical systems; Schultz-Zimm polymer and relativistic particle. Further salient issues such as the $PSL(2;R)$ transformation properties of Tsallis' inverse-temperature parameter and a grandcanonical ensemble of fluctuating random loops related to the Tsallis-Havrda-Charv\'at statistics are also briefly discussed.

categories: cond-mat.stat-mech math-ph math.MP quant-ph


title: Exact solution of a generalized two-sites Bose-Hubbard model

abstract: I introduce a new parametrization of a bosonic Lax operator for the algebraic Bethe ansatz method with the $gl(2)$-invariant $R$-matrix and use it to present the exact solution of a generalized two-sites Bose-Hubbard model with asymmetric tunnelling. In the no interaction limit I show that the Bethe ansatz equations can be written as a $S^{N-1}$ sphere, where $N$ is the total number of atoms in the condensate.

categories: cond-mat.quant-gas hep-th math-ph math.MP nlin.SI quant-ph


title: Nonlinear $E$-mode clustering in Lagrangian space

abstract: We study the nonlinear $E$-mode clustering in Lagrangian space by using large scale structure $N$-body simulations and use the displacement field information in Lagrangian space to recover the primordial linear density field. We find that, compared to Eulerian nonlinear density fields, the $E$-mode displacement fields in Lagrangian space improves the cross-correlation scale $k$ with initial density field by a factor of 6-7, containing two orders of magnitude more primordial information. This illustrates ability of potential density reconstruction algorithms, to improve the baryonic acoustic oscillation measurements from current and future large scale structure surveys.


title: Unusual UUDD magnetic chain structure of the spin-1/2 tetragonally distorted spinel GeCu2O4

abstract: GeCu2O4 exhibits a tetragonal spinel structure due to the strong Jahn-Teller distortion associated with Cu2+ ions. We show that its magnetic structure can be described as slabs composed of a pair of layers with orthogonally oriented spin 1/2 Cu chains in the basal ab plane. The spins between the two layers within a slab are collinearly aligned while the spin directions of neighboring slabs are perpendicular to each other. Interestingly, we find that spins along each chain form an unusual up-up-down-down (UUDD) pattern, suggesting a non-negligible nearest-neighbor biquadratic exchange interaction in the effective classical spin Hamiltonian. We hypothesize that spin-orbit coupling and orbital mixing of Cu2+ ions in this system is non-negligible, which calls for future calculations using perturbation theory with extended Hilbert (spin and orbital) space and calculations based on density functional theory including spin-orbit coupling and looking at the global stability of the UUDD state.


title: Thermal cascaded lattice Boltzmann method

abstract: In this paper, a thermal cascaded lattice Boltzmann method (TCLBM) is developed in combination with the double-distribution-function (DDF) approach. A density distribution function relaxed by the cascaded scheme is employed to solve the flow field, and a total energy distribution function relaxed by the BGK scheme is used to solve temperature field, where two distribution functions are coupled naturally. The forcing terms are incorporated by means of central moments, which is consistent with the previous force scheme [Premnath \emph{et al.}, Phys. Rev. E \textbf{80}, 036702 (2009)] but the derivation is more intelligible and the evolution process is simpler. In the method, the viscous heat dissipation and compression work are taken into account, the Prandtl number and specific-heat ratio are adjustable, the external force is considered directly without the Boussinesq assumption, and the low-Mach number compressible flows can also be simulated. The forcing scheme is tested by simulating a steady Taylor-Green flow, and the new TCLBM is then tested by numerical experiments of the thermal Couette flow, low-Mach number shock tube problem, and Rayleigh-B\'{e}nard convection. The simulation results agree well with the analytical solutions and/or results given by previous researchers.


title: A Simple Picosecond Pulse Generator Based on a Pair of Step Recovery Diodes

abstract: A picosecond pulse generator based on a pair of step recovery diodes (SRD), leveraging the transient response of the SRD PN junction and controlling the pulse width by a resistor, is proposed. We first explain the operation principle of the device, decomposing the pulse generation into different phases, and then demonstrate an experimental prototype with two different resistance, and hence pulse width, values.


title: Online Classification with Complex Metrics

abstract: We present a framework and analysis of consistent binary classification for complex and non-decomposable performance metrics such as the F-measure and the Jaccard measure. The proposed framework is general, as it applies to both batch and online learning, and to both linear and non-linear models. Our work follows recent results showing that the Bayes optimal classifier for many complex metrics is given by a thresholding of the conditional probability of the positive class. This manuscript extends this thresholding characterization -- showing that the utility is strictly locally quasi-concave with respect to the threshold for a wide range of models and performance metrics. This, in turn, motivates simple normalized gradient ascent updates for threshold estimation. We present a finite-sample regret analysis for the resulting procedure. In particular, the risk for the batch case converges to the Bayes risk at the same rate as that of the underlying conditional probability estimation, and the risk of proposed online algorithm converges at a rate that depends on the conditional probability estimation risk. For instance, in the special case where the conditional probability model is logistic regression, our procedure achieves $O(\frac{1}{\sqrt{n}})$ sample complexity, both for batch and online training. Empirical evaluation shows that the proposed algorithms out-perform alternatives in practice, with comparable or better prediction performance and reduced run time for various metrics and datasets.


title: Role of a triangle singularity in the $\gamma p\rightarrow K^+ \Lambda(1405)$ reaction

abstract: We show the effects of a triangle singularity mechanism for the $\gamma p \to K^+\Lambda(1405)$ reaction. The mechanism has a $N^*$ resonance around 2030 MeV, which decays into $K^* \Sigma$. The $K^*$ decays to $K^+ \pi$, and the $\pi \Sigma$ merge to form the $\Lambda(1405)$. This mechanism produces a peak around $\sqrt{s} = 2110$ MeV, and has its largest contribution around cos$\theta=0$. The addition of this mechanism to other conventional ones, leads to a good reproduction of ${\rm d}\sigma/{\rm dcos}\theta$ and the integrated cross section around this energy, providing a solution to a problem encountered in previous theoretical models.


title: Verified Parallel String Matching in Haskell

abstract: In this paper, we prove correctness of parallelizing a string matcher using Haskell as a theorem prover. We use refinement types to specify correctness properties, Haskell terms to express proofs and Liquid Haskell to check correctness of proofs. First, we specify and prove that a class of monoid morphisms can be parallelized via parallel monoid concatenation. Then, we encode string matching as a morphism to get a provably correct parallel transformation. Our 1839LoC prototype proof shows that Liquid Haskell can be used as a fully expressive theorem prover on realistic Haskell implementations.


title: Cross Device Matching for Online Advertising with Neural Feature Ensembles : First Place Solution at CIKM Cup 2016

abstract: We describe the 1st place winning approach for the CIKM Cup 2016 Challenge. In this paper, we provide an approach to reasonably identify same users across multiple devices based on browsing logs. Our approach regards a candidate ranking problem as pairwise classification and utilizes an unsupervised neural feature ensemble approach to learn latent features of users. Combined with traditional hand crafted features, each user pair feature is fed into a supervised classifier in order to perform pairwise classification. Lastly, we propose supervised and unsupervised inference techniques.


title: Iterative coupling of flow, geomechanics and adaptive phase-field fracture including level-set crack width approaches

abstract: In this work, we present numerical studies of fixed-stress iterative coupling for solving flow and geomechanics with propagating fractures in a porous medium. Specifically, fracture propagations are described by employing a phase-field approach. The extension to fixed-stress splitting to propagating phase-field fractures and systematic investigation of its properties are important enhancements to existing studies. Moreover, we provide an accurate computation of the fracture opening using level-set approaches and a subsequent finite element interpolation of the width. The latter enters as fracture permeability into the pressure diffraction problem which is crucial for fluid filled fractures. Our developments are substantiated with several numerical tests that include comparisons of computational cost for iterative coupling and nonlinear and linear iterations as well as convergence studies in space and time.


title: Adaptive enriched Galerkin methods for miscible displacement problems with entropy residual stabilization

abstract: We present a novel approach to the simulation of miscible displacement by employing adaptive enriched Galerkin finite element methods (EG) coupled with entropy residual stabilization for transport. In particular, numerical simulations of viscous fingering instabilities in heterogeneous porous media and Hele-Shaw cells are illustrated. EG is formulated by enriching the conforming continuous Galerkin finite element method (CG) with piecewise constant functions. The method provides locally and globally conservative fluxes, which is crucial for coupled flow and transport problems. Moreover, EG has fewer degrees of freedom in comparison with discontinuous Galerkin (DG) and an efficient flow solver has been derived which allows for higher order schemes. Dynamic adaptive mesh refinement is applied in order to save computational cost for large-scale three dimensional applications. In addition, entropy residual based stabilization for high order EG transport systems prevents any spurious oscillations. Numerical tests are presented to show the capabilities of EG applied to flow and transport.


title: Tangle Floer homology and cobordisms between tangles

abstract: We introduce a generalization of oriented tangles, which are still called tangles, so that they are in one-to-one correspondence with the sutured manifolds. We define cobordisms between sutured manifolds (tangles) by generalizing cobordisms between oriented tangles. For every commutative algebra A over Z/2Z, we define A-Tangles to be the category consisting of A-tangles, which are balanced tangles with A-colorings of the tangle strands and fixed SpinC structures, and A-cobordisms as morphisms. An A-cobordism is a cobordism with a compatible A-coloring and an affine set of SpinC structures. Associated with every A-module M we construct a functor $HF^M$ from A-Tangles to A-Modules, called the tangle Floer homology functor, where A-Modules denotes the the category of A-modules and A-homomorphisms between them. Moreover, for any A-tangle T the A-module $HF^M(T)$ is the extension of sutured Floer homology defined in an earlier work of the authors. In particular, this construction generalizes the 4-manifold invariants of Ozsvath and Szabo. Moreover, applying the above machinery to decorated cobordisms between links, we get functorial maps on link Floer homology.


title: Analysis of Count Data by Transmuted Geometric Distribution

abstract: Transmuted geometric distribution (TGD) was recently introduced and investigated by Chakraborty and Bhati (2016). This is a flexible extension of geometric distribution having an additional parameter that determines its zero inflation as well as the tail length. In the present article we further study this distribution for some of its reliability, stochastic ordering and parameter estimation properties. In parameter estimation among others we discuss an EM algorithm and the performance of estimators is evaluated through extensive simulation. For assessing the statistical significance of additional parameter, Likelihood ratio test, the Rao's score tests and the Wald's test are developed and its empirical power via simulation were compared. We have demonstrate two applications of (TGD) in modeling real life count data.


title: An Alternative Method to Implement Contact Angle Boundary Condition on Immersed Surfaces for Phase-Field Simulations

abstract: In this paper, we propose an alternative approach to implement the contact angle boundary condition on immersed surfaces for phase-field simulations of two-phase flows using the Cahn-Hilliard equation on a Cartesian mesh. This simple and effective method was inspired by previous works on the geometric formulation of the wetting boundary condition. In two dimensions, by making full use of the hyperbolic tangent profile of the order parameter, we were able to obtain its unknown value at a ghost point from the information at only one point in the fluid. This is in contrast with previous approaches using interpolations involving several points. The special feature allows this method to be easily implemented on immersed surfaces (including curved ones) that cut through the grid lines. It is verified through the study of two examples: (1) the shape of a drop on a circular cylinder with different contact angles; (2) the spreading of a drop on an embedded inclined wall with a given contact angle.


title: On the hyperplane conjecture for periods of Calabi-Yau hypersurfaces in $\mathbb P^n$

abstract: In [HLY1], Hosono, Lian, and Yau posed a conjecture characterizing the set of solutions to certain Gelfand-Kapranov-Zelevinsky hypergeometric equations which are realized as periods of Calabi-Yau hypersurfaces in a Gorenstein Fano toric variety $X$. We prove this conjecture in the case where $X$ is a complex projective space.


title: On Unifying Multi-View Self-Representations for Clustering by Tensor Multi-Rank Minimization

abstract: In this paper, we address the multi-view subspace clustering problem. Our method utilizes the circulant algebra for tensor, which is constructed by stacking the subspace representation matrices of different views and then rotating, to capture the low rank tensor subspace so that the refinement of the view-specific subspaces can be achieved, as well as the high order correlations underlying multi-view data can be explored.} By introducing a recently proposed tensor factorization, namely tensor-Singular Value Decomposition (t-SVD) \cite{kilmer13}, we can impose a new type of low-rank tensor constraint on the rotated tensor to capture the complementary information from multiple views. Different from traditional unfolding based tensor norm, this low-rank tensor constraint has optimality properties similar to that of matrix rank derived from SVD, so the complementary information among views can be explored more efficiently and thoroughly. The established model, called t-SVD based Multi-view Subspace Clustering (t-SVD-MSC), falls into the applicable scope of augmented Lagrangian method, and its minimization problem can be efficiently solved with theoretical convergence guarantee and relatively low computational complexity. Extensive experimental testing on eight challenging image dataset shows that the proposed method has achieved highly competent objective performance compared to several state-of-the-art multi-view clustering methods.


title: The quadratic regulator problem and the Riccati equation for a process governed by a linear Volterra integrodifferential equations

abstract: In this paper we study the quadratic regulator problem for a process governed by a Volterra integral equation in ${\mathbb R}^n$. Our main goal is the proof that it is possible to associate a Riccati differential equation to this quadratic control problem, which leads to the feedback form of the optimal control. This is in contrast with previous papers on the subject, which confine themselves to study the Fredholm integral equation which is solved by the optimal control.


title: Linear optical quantum metrology with single photons --- Experimental errors, resource counting, and quantum Cram\'er-Rao bounds

abstract: Quantum number-path entanglement is a resource for super-sensitive quantum metrology and in particular provides for sub-shotnoise or even Heisenberg-limited sensitivity. However, such number-path entanglement has thought to have been resource intensive to create in the first place --- typically requiring either very strong nonlinearities, or nondeterministic preparation schemes with feed-forward, which are difficult to implement. Recently in [Phys. Rev. Lett. 114, 170802 (2015)] we showed that number-path entanglement from a BosonSampling inspired interferometer can be used to beat the shot-noise limit. In this manuscript we compare and contrast different interferometric schemes, discuss resource counting, calculate exact quantum Cram\'er-Rao bounds, and study details of experimental errors.


title: Developing and Assessing MATLAB Exercises for Active Concept Learning

abstract: New technologies, such as MOOCs, provide innovative methods to tackle new challenges in teaching and learning, such as globalization and changing contemporary culture and to remove the limits of conventional classrooms. However, they also bring challenges in course delivery and assessment, due to factors such as less direct student-instructor interaction. These challenges are especially severe in engineering education, which relies heavily on experiential learning, such as computer simulations and laboratory exercises, to assist students in understanding concepts. As a result, effective design of experiential learning components is extremely critical for engineering MOOCs. In this paper, we will share our experience gained through developing and offering a MOOC on communication systems, with special focus on the development and assessment of MATLAB exercises for active concept learning. Our approach introduced students to concepts using learning components commonly provided by many MOOC platforms (e.g., online lectures and quizzes), and augmented the student experience with MATLAB based computer simulations and exercises to enable more concrete and detailed understanding of the material. We describe here a systematic approach to MATLAB problem design and assessment, based on our experience with the MATLAB server provided by MathWorks and integrated with the edX MOOC platform. We discuss the effectiveness of the instructional methods as evaluated through students' learning performance. We analyze the impact of the course design tools from both the instructor and the student perspective.


title: Quantum Trajectories: Dirac, Moyal and Bohm

abstract: We recall Dirac's early proposals to develop a description of quantum phenomena in terms of a non-commutative algebra in which he suggested a way to construct what he called `quantum trajectories'. Generalising these ideas, we show how they are related to weak values and explore their use in the experimental construction of quantum trajectories. We discuss covering spaces which play an essential role in accounting for the `wave' properties of quantum particles. We briefly point out how new mathematical techniques take us beyond Hilbert space and into a deeper structure which connects with the algebras originally introduced by Born, Heisenberg and Jordan. This enables us to bring out the geometric aspects of quantum phenomena.


title: Asymptotic of Non-Crossings probability of Additive Wiener Fields

abstract: Let $W_i=\{W_i(t_i), t_i\in \R_+\}, i=1,2,\ldots,d$ are independent Wiener processes. $W=\{W(\mathbf{t}),t\in \R_+^d\}$ be the additive Wiener field define as the sum of $W_i$. For any trend $f$ in $\kHC$ (the reproducing kernel Hilbert Space of $W$), we derive upper and lower bounds for the boundary non-crossing probability $$P_f=P\{\sum_{i=1}^{d}W_i(t_i) +f(\mathbf{t})\leq u(\mathbf{t}), \mathbf{t}\in\R_+^d\},$$ where $u: \R_+^d\rightarrow \R_+$ is a measurable function. Furthermore, for large trend functions $\gamma f>0$, we show that the asymptotically relation $\ln P_{\gamma f}\sim \ln P_{\gamma \underline{f}}$ as $\gamma \to \IF$, where $\underline{f}$ is the projection of $f$ on some closed convex subset of $\kHC$.


title: Emergence of topological semimetals in gap closing in semiconductors without inversion symmetry

abstract: A band gap for electronic states in crystals governs various properties of solids, such as transport, optical and magnetic properties. Its estimation and control have been an important issue in solid state physics. The band gap can be controlled externally by various parameters, such as pressure, atomic compositions and external field. Sometimes, the gap even collapses by tuning some parameter. In the field of topological insulators, such closing of the gap at a time-reversal invariant momentum indicates a band inversion, i.e. it leads to a topological phase transition from a normal insulator to a topological insulator. Here we show that the gap losing in inversion-asymmetric crystals is universal, in the sense that the gap closing always leads either to a Weyl semimetal or a nodal-line semimetal, from an exhaustive study on possible space groups. We here consider three-dimensional spinful systems with time-reversal symmetry. The space group of the system and the wavevector at the gap closing uniquely determine which possibility occurs and where the gap-closing points or lines lie in the wavevector space after closing of the gap. In particular, we show that an insulator-to-insulator transition never happens, which is in sharp contrast with inversion-symmetric systems.


title: On the equivalence of proportional-integral and proportional-resonant controllers with anti-windup

abstract: It is shown that proportional-integral (PI) control in the synchronously rotating (d, q)-reference frame and proportional-resonant (PR) control in the stationary ({\alpha}, \b{eta})-reference frame, both with anti-windup, are equivalent if and only if their implementation is done correctly in state space and the controller parameters and the initial values are identical. It is shown that an equivalence in the frequency domain does only hold if simplifying assumptions are satisfied. As consequence of the equivalence, both closed- loop control performances are identical with respect to closed-loop dynamics and steady-state accuracy. The control performance will only differ if their implementation is not done correctly or the time delay induced by the voltage source inverter becomes significant. To the best knowledge of the author, equivalence of PR and PI controllers with anti-windup has not been shown before (in particular not in state space).


title: Kinematic Space and Wormholes

abstract: The kinematic space could play a key role in constructing the bulk geometry from dual CFT. In this paper, we study the kinematic space from geometric points of view, without resorting to differential entropy. We find that the kinematic space could be intrinsically defined in the embedding space. For each oriented geodesic in the Poincar\'e disk, there is a corresponding point in the kinematic space. This point is the tip of the causal diamond of the disk whose intersection with the Poincar\'e disk determines the geodesic. In this geometric construction, the causal structure in the kinematic space can be seen clearly. Moreover, we find that every transformation in the $SL(2,\mathbb{R})$ leads to a geodesic in the kinematic space. In particular, for a hyperbolic transformation defining a BTZ black hole, it is a timelike geodesic in the kinematic space. We show that the horizon length of the static BTZ black hole could be computed by the geodesic length of corresponding points in the kinematic space. Furthermore, we discuss the fundamental regions in the kinematic space for the BTZ blackhole and multi-boundary wormholes.


title: A new earthquake location method based on the waveform inversion

abstract: In this paper, a new earthquake location method based on the waveform inversion is proposed. As is known to all, the waveform misfit function is very sensitive to the phase shift between the synthetic waveform signal and the real waveform signal. Thus, the convergence domain of the conventional waveform based earthquake location methods is very small. In present study, by introducing and solving a simple sub-optimization problem, we greatly expand the convergence domain of the waveform based earthquake location method. According to a large number of numerical experiments, the new method expands the range of convergence by several tens of times. This allows us to locate the earthquake accurately even from some relatively bad initial values.


title: The first Cheeger constant of a simplex

abstract: The coboundary expansion generalizes the classical graph expansion to the case of the general simplicial complexes, and allows the definition of the higher-dimensional Cheeger constants $h_k(X)$ for an arbitrary simplicial complex $X$, and any $k\geq 0$. In this paper we investigate the value of $h_1(\Delta^{[n]})$ - the first Cheeger constant of a simplex with $n$ vertices. It is known, due to the pioneering work of Meshulam and Wallach, that \[\lceil n/3\rceil\geq h_1(\Delta^{[n]})\geq n/3, \textrm{ for all } n,\] and that the equality $h_1(\Delta^{[n]})=n/3$ is achieved when $n$ is divisible by $3$. Here we expand on these results. First, we show that \[h_1(\Delta^{[n]})=n/3, \textrm{ whenever }n\textrm{ is not a power of }2.\] So the sharp equality holds on a set whose density goes to $1$. Second, we show that \[h_1(\Delta^{[n]})=n/3+O(1/n),\textrm{ when }n\textrm{ is a power of }2.\] In other words, as $n$ goes to infinity, the value $h_1(\Delta^{[n]})-n/3$ is either $0$ or goes to $0$ very rapidly. Our methods include recasting the original question in purely graph-theoretic language, followed by a detailed investigation of a specific graph family, the so-called {\it staircase graphs}. These are defined by associating a graph to every partition, and appear to be especially suited to gain information about the first Cheeger constant of a simplex.


title: Bosonic-Seesaw Portal Dark Matter

abstract: We discuss a new type of Higgs-portal dark matter (DM)-production mechanism, called bosonic-seesaw portal (BSP) scenario. The BS provides the dynamical origin of the electroweak symmetry breaking, triggered by mixing between the elementary Higgs and a composite Higgs generated by a new-color strong dynamics, hypercolor (HC). At the HC strong coupling scale, the classical-scale invariance assumed in the model is dynamically broken as well as the "chiral" symmetry present in the HC sector. In addition to the composite Higgs, HC baryons emerge to potentially be stable because of the unbroken HC baryon number symmetry. Hence the lightest HC baryon can be a DM candidate. Of interest in the present scenario is that HC pions can be as heavy as the HC baryon due to the possibly enhanced-explicit "chiral"-breaking effect triggered after the BS mechanism, so the HC baryon pair cannot annihilate into HC pions. As in the standard setup of freeze-in scenario, it is assumed that the DM was never in the thermal equilibrium, which ends up with no thermal abundance. It is then the non-thermal BSP process that crucially comes into the game below the HC scale: the HC baryon significantly couples to the standard-model Higgs via the BS mechanism, and can non-thermally be produced from the thermal plasma below the HC scale, which turns out to allow the TeV mass scale for the composite baryonic DM, much smaller than the generic bound placed in the conventional thermal freeze-out scenario, to account for the observed relic abundance. Thus the DM can closely be related to the mechanism of the electroweak symmetry breaking.


title: Higher moments of multiplicity fluctuations in a hadron-resonance gas with exact conservation laws

abstract: Higher moments of multiplicity fluctuations of hadrons produced in central nucleus-nucleus collisions are studied within the hadron-resonance gas model in the canonical ensemble. Exact conservation of three charges, baryon number, electric charge, and strangeness, is enforced in the large volume limit. Moments up to the forth order of various particles are calculated at SPS, RHIC and LHC energies. The asymptotic fluctuations within a simplified model with only one conserved charge in the canonical ensemble are discussed where simple analytical expressions for moments of multiplicity distributions can be obtained. Moments products of net-proton, net-kaon, and net-charge distributions in Au + Au collisions at RHIC energies are calculated. The pseudo-rapidity coverage dependence of net-charge fluctuation is discussed.


title: A Monte Carlo Study of Flux Ratios of Raman Scattered O~VI Features at 6825 \AA\ and 7082 \AA\ in Symbiotic Stars

abstract: Symbiotic stars are regarded as wide binary systems consisting of a hot white dwarf and a mass losing giant. They exhibit unique spectral features at 6825 \AA\ and 7082 \AA, which are formed via Raman scattering of \ion{O}{6}$\lambda\lambda$ 1032 and 1038 with atomic hydrogen. We adopt a Monte Carlo technique to generate the same number of \ion{O}{6}$\lambda$1032 and $\lambda$1038 line photons and compute the flux ratio $F(6825)/F(7082)$ of these Raman scattered \ion{O}{6} features formed in neutral regions with a simple geometric shape as a function of \ion{H}{1} column density $N_{HI}$. In cylindrical and spherical neutral regions with the \ion{O}{6} source embedded inside, the flux ratio $F(6825)/F(7082)$ shows an overall decrease from 3 to 1 as $N_{HI}$ increases in the range $10^{22-24}{\rm\ cm^{-2}}$. In the cases of a slab geometry and other geometries with the \ion{O}{6} source outside the \ion{H}{1} region, Rayleigh escape operates to lower the flux ratio considerably. For moderate values of $N_{HI}\sim 10^{23}{\rm\ cm^{-2}}$ the flux ratio behaves in a complicated way to exhibit a broad bump with a peak value of 3.5 in the case of a sphere geometry. We find that the ratio of Raman conversion efficiencies of \ion{O}{6}$\lambda\lambda$1032, 1038 ranges from 0.8 to 3.5. Our high resolution spectra of 'D' type HM~Sge and 'S' type AG~Dra obtained with the {\it Canada-France-Hawaii-Telescope} show that the flux ratio $F(6825)/F(7082)$ of AG~Dra is significantly smaller than that of HM~Sge, implying that 'S' type symbiotics are characterized by higher $N_{HI}$ than 'D' type symbiotics.


title: Finite Temperature Corrections to Tachyon Mass in Intersecting D-Branes

abstract: We continue with the analysis of finite temperature corrections to the Tachyon mass in intersecting branes which was initiated in arxiv:1403.0389. In this paper we extend the computation to the case of intersecting D3-branes by considering a setup of two intersecting branes in flat-space background. A holographic model dual to BCS superconductor consisting of intersecting D8-branes in D4-brane background was proposed in arxiv:1104.2843. The background considered here is a simplified configuration of this dual model. We compute the one-loop Tachyon amplitude in the Yang-Mills approximation and show that the result is finite. Analyzing the amplitudes further we numerically compute the transition temperature at which the Tachyon becomes massless. The analytic expressions for the one-loop amplitudes obtained here reduce to those for intersecting D1-branes obtained in arxiv:1403.0389 as well as those for intersecting D2-branes.


title: An Efficient Web Traffic Defence Against Timing-Analysis Attacks

abstract: We introduce a new class of lower overhead tunnel that is resistant to traffic analysis. The tunnel opportunistically reduces the number of dummy packets transmitted during busy times when many flows are simultaneously active while maintaining well-defined privacy properties. We find that the dummy packet overhead is typically less than 20% on lightly loaded links and falls to zero as the traffic load increases i.e. the tunnel is capacity-achieving. The additional latency incurred is less than 100ms. We build an experimental prototype of the tunnel and carry out an extensive performance evaluation that demonstrates its effectiveness under a range of network conditions and real web page fetches.


title: Weyl Ferroelectric Semimetal

abstract: The recent discoveries of ferroelectric metal and Weyl semimetal (WSM) have stimulated a natural question: whether these two exotic states of matter can coexist in a single material or not. These two discoveries ensure us that physically it is possible since both of them share the same necessary condition, the broken inversion symmetry. Here, by using first-principles calculations, we demonstrate that the experimentally synthesized nonmagnetic HgPbO$_3$ represents a unique example of such hybrid "\emph{Weyl ferroelectric semimetal}". Its centrosymmetric $R\bar{3}c$ phase will undergo a ferroelectric phase transition to the ferroelectric $R3c$ structure. Both phases are metallic and the ferroelectric phase owns a spontaneous polarization of 33 $\mu$C/cm$^2$. Most importantly, it also harbors six pairs of chiral Weyl nodes around the Fermi level to be an oxide WSM. The structural symmetry broken phase transition induces a topological phase transition. The coexistence of ferroelectricity and Weyl nodes in HgPbO$_3$ is an ideal platform for exploring multiphase interaction and mutual control. The Weyl nodes can be tuned by external pulse electric field, which is promising for potential applications of integrated topotronic and ferroelectric devices.


title: Physics potential of searching for $0\nu\beta\beta$ decays in JUNO

abstract: In the past few decades, numerous searches have been made for the neutrinoless double-beta decay (0$\nu\beta\beta$) process, aiming to establish whether neutrinos are their own antiparticles (Majorana neutrinos), but no 0$\nu\beta\beta$ decay signal has yet been observed. A number of new experiments are proposed but they ultimately suffer from a common problem: the sensitivity may not increase indefinitely with the target mass. We have performed a detailed analysis of the physics potential by using the Jiangmen Underground Neutrino Observatory (JUNO) to improve the sensitivity to 0$\nu\beta\beta$ up to a few meV, a major step forward with respect to the experiments currently being planned. JUNO is a 20 kton low-background liquid scintillator (LS) detector with 3\%/$\sqrt{E \text{(MeV)}}$ energy resolution, now under construction. It is feasible to build a balloon filled with enriched xenon gas (with $^{136}$Xe up to 80\%) dissolved in LS, inserted into the central region of the JUNO LS. The energy resolution is $\sim$1.9\% at the $Q$-value of $^{136}$Xe 0$\nu\beta\beta$ decay. Ultra-low background is the key for 0$\nu\beta\beta$ decay searches. Detailed studies of background rates from intrinsic 2$\nu\beta\beta$ and $^{8}$B solar neutrinos, natural radioactivity, and cosmogenic radionuclides (including light isotopes and $^{137}$Xe) were performed and several muon veto schemes were developed. We find that JUNO has the potential to reach a sensitivity (at 90\% C. L.) to $T^{0\nu\beta\beta}_{1/2}$ of $1.8\times10^{28}$ yr ($5.6\times10^{27}$ yr) with $\sim$50 tons (5 tons) of fiducial $^{136}$Xe and 5 years exposure, while in the 50-ton case the corresponding sensitivity to the effective neutrino mass, $m_{\beta\beta}$, could reach (5--12) meV, covering completely the allowed region of inverted neutrino mass ordering.


title: Bounded embeddings of graphs in the plane

abstract: A drawing in the plane ($\mathbb{R}^2$) of a graph $G=(V,E)$ equipped with a function $\gamma: V \rightarrow \mathbb{N}$ is \emph{$x$-bounded} if (i) $x(u) <x(v)$ whenever $\gamma(u)<\gamma(v)$ and (ii) $\gamma(u)\leq\gamma(w)\leq \gamma(v)$, where $uv\in E$ and $\gamma(u)\leq \gamma(v)$, whenever $x(w)\in x(uv)$, where $x(.)$ denotes the projection to the $x$-axis. We prove a characterization of isotopy classes of graph embeddings in the plane containing an $x$-bounded embedding. Then we present an efficient algorithm, that relies on our result, for testing the existence of an $x$-bounded embedding if the given graph is a tree or generalized $\Theta$-graph. This partially answers a question raised recently by Angelini et al. and Chang et al., and proves that c-planarity testing of flat clustered graphs with three clusters is tractable if each connected component of the underlying abstract graph is a tree.


title: Sequential decision problems, dependent types and generic solutions

abstract: We present a computer-checked generic implementation for solving finite-horizon sequential decision problems. This is a wide class of problems, including inter-temporal optimizations, knapsack, optimal bracketing, scheduling, etc. The implementation can handle time-step dependent control and state spaces, and monadic representations of uncertainty (such as stochastic, non-deterministic, fuzzy, or combinations thereof). This level of genericity is achievable in a programming language with dependent types (we have used both Idris and Agda). Dependent types are also the means that allow us to obtain a formalization and computer-checked proof of the central component of our implementation: Bellman's principle of optimality and the associated backwards induction algorithm. The formalization clarifies certain aspects of backwards induction and, by making explicit notions such as viability and reachability, can serve as a starting point for a theory of controllability of monadic dynamical systems, commonly encountered in, e.g., climate impact research.


title: Note on the butterfly effect in holographic superconductor models

abstract: In this note we remark that the butterfly effect can be used to diagnose the phase transition of superconductivity in a holographic framework. Specifically, we compute the butterfly velocity in a charged black hole background as well as anisotropic backgrounds with Q-lattice structure. In both cases we find its derivative to the temperature is discontinuous at critical points. We also propose that the butterfly velocity can signalize the occurrence of thermal phase transition in general holographic models.


title: Another characterization of homogeneous Poisson processes

abstract: For a general renewal process $N$ (allowing delay, defect and multiple simultaneous arrivals) the independence of the first renewal epochs of the marked processes got from $N$ by Bernoulli $0$/$1$ thinning is characterized. This independence is well-known to hold true in the case of homogeneous Poisson processes; by way of corollary one obtains the interesting observation that, when coupled with some minimal extra conditions, it in fact already identifies them. The proof is analytic in character.


title: Revisiting optimal eavesdropping in quantum cryptography: Optimal interaction is unique up to rotation of the underlying basis

abstract: A general framework of optimal eavesdropping on BB84 protocol was provided by Fuchs et al. [Phys. Rev. A, 1997]. An upper bound on mutual information was derived, which could be achieved by a specific type of interaction and the corresponding measurement. However, uniqueness of optimal interaction was posed as an unsolved problem there and it has remained open for almost two decades now. In this paper, we solve this open problem and establish the uniqueness of optimal interaction up to rotation. The specific choice of optimal interaction by Fuchs et al. is shown to be a special case of the form derived in our work.


title: Two are Better than One: An Ensemble of Retrieval- and Generation-Based Dialog Systems

abstract: Open-domain human-computer conversation has attracted much attention in the field of NLP. Contrary to rule- or template-based domain-specific dialog systems, open-domain conversation usually requires data-driven approaches, which can be roughly divided into two categories: retrieval-based and generation-based systems. Retrieval systems search a user-issued utterance (called a query) in a large database, and return a reply that best matches the query. Generative approaches, typically based on recurrent neural networks (RNNs), can synthesize new replies, but they suffer from the problem of generating short, meaningless utterances. In this paper, we propose a novel ensemble of retrieval-based and generation-based dialog systems in the open domain. In our approach, the retrieved candidate, in addition to the original query, is fed to an RNN-based reply generator, so that the neural model is aware of more information. The generated reply is then fed back as a new candidate for post-reranking. Experimental results show that such ensemble outperforms each single part of it by a large margin.


title: Holomorphic Gromov's partial order

abstract: As in [5], we study holomorphic maps of positive degree between compact complex manifolds, and prove that any holomorphic map of degree one from a compact complex manifold to itself is biholomorphic. This conclusion confirms that under a mild restriction the holomorphic Gromov relation ">_" is indeed a partial order.


title: In silico modeling of the rheological properties of covalently crosslinked collagen triple helices

abstract: Biomimetic hydrogels based on natural polymers are a promising class of biomaterial, mimicking the natural extra-cellular matrix of biological tissues and providing cues for cell attachment, proliferation and differentiation. With a view to providing an upstream method to guide subsequent experimental design, the aim of this study was to introduce a mathematical model that described the rheological properties of a hydrogel system based on covalently crosslinked collagen triple helices. In light of their organization, such gels exhibit limited collagen bundling that cannot be described by existing fibril network models. The model presented here treats collagen triple helices as discrete semi-flexible polymers, permits full access to metrics for network microstructure, and should provide a comprehensive understanding of the parameter space associated with the development of such multi-functional materials. Triple helical hydrogel networks were experimentally obtained via reaction of type I collagen with both aromatic and aliphatic diacids. The complex modulus G* was found from rheological testing in linear shear and quantitatively compared to model predictions. In silico data from the computational model successfully described the experimental trends in hydrogel storage modulus with either (i) the concentration of collagen triple helices during crosslinking reaction or (ii) the type of crosslinking segment introduced in resulting hydrogel networks. This approach may pave the way to a step change in the rational design of biomimetic triple helical collagen systems with controlled multi-functionality.