title: Magnetic reconnection at the earliest stage of solar flux emergence

abstract: On 2016 September 20, the Interface Region Imaging Spectrograph observed an active region during its earliest emerging phase for almost 7 hours. The Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory observed continuous emergence of small-scale magnetic bipoles with a rate of $\sim$10$^{16}$ Mx~s$^{-1}$. The emergence of magnetic fluxes and interactions between different polarities lead to frequent occurrence of ultraviolet (UV) bursts, which exhibit as intense transient brightenings in the 1400 \AA{} images. In the meantime, discrete small patches with the same magnetic polarity tend to move together and merge, leading to enhancement of the magnetic fields and thus formation of pores (small sunspots) at some locations. The spectra of these UV bursts are characterized by the superposition of several chromospheric absorption lines on the greatly broadened profiles of some emission lines formed at typical transition region temperatures, suggesting heating of the local materials to a few tens of thousands of kelvin in the lower atmosphere by magnetic reconnection. Some bursts reveal blue and red shifts of $\sim$100~km~s$^{-1}$ at neighboring pixels, indicating the spatially resolved bidirectional reconnection outflows. Many such bursts appear to be associated with the cancellation of magnetic fluxes with a rate of the order of $\sim$10$^{15}$ Mx~s$^{-1}$. We also investigate the three-dimensional magnetic field topology through a magneto-hydrostatic model and find that a small fraction of the bursts are associated with bald patches (magnetic dips). Finally, we find that almost all bursts are located in regions of large squashing factor at the height of $\sim$1 Mm, reinforcing our conclusion that these bursts are produced through reconnection in the lower atmosphere.


title: Weak convergence of branched conformal immersions with uniformly bounded areas and Willmore energies

abstract: In this paper, we firstly extend Theorem 5.1.1 in \cite {Helein} due to H\'elein to a rescaled branched conformal immersed sequence(c.f. Theorem 1.5). By virtue of this local convergence theorem, we study the blowup behavior of a sequence of branched conformal immersions of closed Riemannian surface in $\mathbb{R}^{n}$ with uniformly bounded areas and Willmore energies. Furthermore, we prove that the integral identity of Gauss curvature is true.


title: Yamabe Equation on Some Complete Noncompact Manifolds

abstract: In this paper, we consider the Yamabe equation on a complete noncompact Riemannian manifold and find some geometric conditions on the manifold such that the Yamabe problem admits a bounded positive solution.


title: Cosmological Theories Of The Extra Terms

abstract: This paper purports to have: Introduced a new formulation of Quantum Mechanics, explained the apparent disconnect between Quantum Mechanics and General Relativity, explained the observed far field expansion of the Universe (Dark Energy), supplied an argument which goes towards explaining away Dark Matter (there are modelling difficulties) and not explained, on the basis of gravitational theory, the Voyager Anomaly.


title: Beautiful mathematics for beauty-full and other multi-heavy hadronic systems

abstract: In most non-perturbative methods in hadron physics the calculations are started with a correlation function in terms of some interpolating and transition currents in $ x $-space. For simplicity, the calculations are then transformed to the momentum space by a Fourier transformation. To suppress the contributions of the higher states and continuum; and enhance the ground state contribution, Borel transformation as well as continuum subtraction are applied by the help of quark-hadron duality assumption. In the present study we work out the mathematics required for these processes in the case of light and multi-heavy hadrons. We address a well-known problem in subtraction of the effects of the higher states and continuum and discuss how we find finite results without any divergence by using an appropriate representation of the modified Bessel functions, appearing in the heavy quark propagator, and successive applications of the Borel transformations, which lead to more suppression of the higher states and continuum contributions. The results obtained can be used in determination of the spectroscopic and decay properties of the multi-heavy standard and non-conventional (exotic) systems in many non-perturbative methods, specially the QCD sum rules.


title: Decoupled Learning for Conditional Adversarial Networks

abstract: Incorporating encoding-decoding nets with adversarial nets has been widely adopted in image generation tasks. We observe that the state-of-the-art achievements were obtained by carefully balancing the reconstruction loss and adversarial loss, and such balance shifts with different network structures, datasets, and training strategies. Empirical studies have demonstrated that an inappropriate weight between the two losses may cause instability, and it is tricky to search for the optimal setting, especially when lacking prior knowledge on the data and network. This paper gives the first attempt to relax the need of manual balancing by proposing the concept of \textit{decoupled learning}, where a novel network structure is designed that explicitly disentangles the backpropagation paths of the two losses. Experimental results demonstrate the effectiveness, robustness, and generality of the proposed method. The other contribution of the paper is the design of a new evaluation metric to measure the image quality of generative models. We propose the so-called \textit{normalized relative discriminative score} (NRDS), which introduces the idea of relative comparison, rather than providing absolute estimates like existing metrics.


title: Tunable multi-color coherent light generation in single MgO: PPLN bulk crystal

abstract: We report a theoretical design analysis of domain-engineered periodically poled lithium niobate (PPLN) for wavelength conversion of near-infrared sources to generate coherent light in the visible spectral range. Our analysis on the spectral outputs show that with a proper design of the quasi phase matching (QPM) periods, tunable, multiple nonlinear optical processes can be simultaneously phase matched in a single segmented crystal. We show that a three-segment single PPLN crystal has potential to generate violet (432 nm), blue (490 nm) and orange (600 nm) wavelengths by simultaneous sum frequency and second harmonic generation processes. Such a design scheme has promising potential for a compact, robust and tunable multi-colored visible light source which can find various applications such as in biomedicine, high-density optical data storage and laser based color displays.


title: Attentive Recurrent Tensor Model for Community Question Answering

abstract: A major challenge to the problem of community question answering is the lexical and semantic gap between the sentence representations. Some solutions to minimize this gap includes the introduction of extra parameters to deep models or augmenting the external handcrafted features. In this paper, we propose a novel attentive recurrent tensor network for solving the lexical and semantic gap in community question answering. We introduce token-level and phrase-level attention strategy that maps input sequences to the output using trainable parameters. Further, we use the tensor parameters to introduce a 3-way interaction between question, answer and external features in vector space. We introduce simplified tensor matrices with L2 regularization that results in smooth optimization during training. The proposed model achieves state-of-the-art performance on the task of answer sentence selection (TrecQA and WikiQA datasets) while outperforming the current state-of-the-art on the tasks of best answer selection (Yahoo! L4) and answer triggering task (WikiQA).


title: NOOP: A Domain-Theoretic Model of Nominally-Typed OOP

abstract: The majority of industrial-strength object-oriented (OO) software is written using nominally-typed OO programming languages. Extant domain-theoretic models of OOP developed to analyze OO type systems miss, however, a crucial feature of these mainstream OO languages: nominality. This paper presents the construction of NOOP as the first domain-theoretic model of OOP that includes full class/type names information found in nominally-typed OOP. Inclusion of nominal information in objects of NOOP and asserting that type inheritance in statically-typed OO programming languages is an inherently nominal notion allow readily proving that type inheritance and subtyping are completely identified in these languages. This conclusion is in full agreement with intuitions of developers and language designers of these OO languages, and contrary to the belief that "inheritance is not subtyping," which came from assuming non-nominal (a.k.a., structural) models of OOP. To motivate the construction of NOOP, this paper briefly presents the benefits of nominal-typing to mainstream OO developers and OO language designers, as compared to structural-typing. After presenting NOOP, the paper further briefly compares NOOP to the most widely known domain-theoretic models of OOP. Leveraging the development of NOOP, the comparisons presented in this paper provide clear, brief and precise technical and mathematical accounts for the relation between nominal and structural OO type systems. NOOP, thus, provides a firmer semantic foundation for analyzing and progressing nominally-typed OO programming languages.


title: A Rate-Optimal Construction of Codes with Sequential Recovery with Low Block Length

abstract: An erasure code is said to be a code with sequential recovery with parameters $r$ and $t$, if for any $s \leq t$ erased code symbols, there is an $s$-step recovery process in which at each step we recover exactly one erased code symbol by contacting at most $r$ other code symbols. In earlier work by the same authors, presented at ISIT 2017, we had given a construction for binary codes with sequential recovery from $t$ erasures, with locality parameter $r$, which were optimal in terms of code rate for given $r,t$, but where the block length was large, on the order of $r^{c^t}$, for some constant $c >1$. In the present paper, we present an alternative construction of a rate-optimal code for any value of $t$ and any $r\geq3$, where the block length is significantly smaller, on the order of $r^{\frac{5t}{4}+\frac{7}{4}}$ (in some instances of order $r^{\frac{3t}{2}+2}$). Our construction is based on the construction of certain kind of tree-like graphs with girth $t+1$. We construct these graphs and hence the codes recursively.


title: On the variety of triangles for a hyper-Kaehler fourfold constructed by Debarre and Voisin

abstract: We study the similarities between the Fano varieties of lines on a cubic fourfold, a hyper-Kaehler fourfold studied by Beauville and Donagi, and the hyper-Kaehler fourfold constructed by Debarre and Voisin. We exhibit an analog of the notion of "triangle" for these varieties and prove that the 6-dimensional variety of "triangles" is a Lagrangian subvariety in the cube of the constructed hyper-Kaehler fourfold.


title: The isoperimetric quotient of a convex body decreases monotonically under the Eikonal abrasion model

abstract: We show that under the Eikonal abrasion model, prescribing uniform normal speed in the direction of the inward surface normal, the isoperimetric quotient of a convex shape is decreasing monotonically.


title: Depth CNNs for RGB-D scene recognition: learning from scratch better than transferring from RGB-CNNs

abstract: Scene recognition with RGB images has been extensively studied and has reached very remarkable recognition levels, thanks to convolutional neural networks (CNN) and large scene datasets. In contrast, current RGB-D scene data is much more limited, so often leverages RGB large datasets, by transferring pretrained RGB CNN models and fine-tuning with the target RGB-D dataset. However, we show that this approach has the limitation of hardly reaching bottom layers, which is key to learn modality-specific features. In contrast, we focus on the bottom layers, and propose an alternative strategy to learn depth features combining local weakly supervised training from patches followed by global fine tuning with images. This strategy is capable of learning very discriminative depth-specific features with limited depth images, without resorting to Places-CNN. In addition we propose a modified CNN architecture to further match the complexity of the model and the amount of data available. For RGB-D scene recognition, depth and RGB features are combined by projecting them in a common space and further leaning a multilayer classifier, which is jointly optimized in an end-to-end network. Our framework achieves state-of-the-art accuracy on NYU2 and SUN RGB-D in both depth only and combined RGB-D data.


title: Entanglement between total intensity and polarization for pairs of coherent states

abstract: We examine entanglement between number and polarization, or between number and relative phase, for pairs of coherent states and two-mode squeezed vacuum via linear entropy and covariance criteria. We consider the embedding of the two-mode Hilbert space in a larger space to get a well-defined factorization of the number-phase variables. This can be regarded as a kind of proto-entanglement that can be extracted and converted into real particle entanglement via feasible experimental procedures. In particular this reveals interesting entanglement properties of pairs of coherent states.


title: Universal Origin for Environment-Assisted Quantum Transport in Exciton Transfer Networks

abstract: Environment-assisted quantum transport (ENAQT) is the possibility of an external environment to enhance transport efficiency of quantum particles. This idea has generated much excitement over recent years, especially due to the experimentally-motivated possibility of ENAQT in photo-synthetic exciton transfer complexes. Many theoretical calculations have shown ENAQT, but the explanations for its origin differ, and a universal explanation has been elusive. Here we demonstrate a universal origin for ENAQT in quantum networks with a dephasing environment, based on a relation between exciton current and occupation within a Markovian open quantum system approach. We show that ENAQT appears due to two competing processes, namely the tendency of dephasing to make the exciton population uniform, and the formation of an exciton density gradient, defined by the source and the sink. Furthermore, we find a geometric condition on the network for the appearance of ENAQT, relevant to natural and artificial systems.


title: N=4 Super-Schwarzian Theory on the Coadoint Orbit and PSU(1,1|2)

abstract: An N=4 super-Schwarzian theory is formulated by the coadjoint orbit method. It is discovered that the action has symmetry under PSU(1,1|2).


title: Curvature-based Comparison of Two Neural Networks

abstract: In this paper we show the similarities and differences of two deep neural networks by comparing the manifolds composed of activation vectors in each fully connected layer of them. The main contribution of this paper includes 1) a new data generating algorithm which is crucial for determining the dimension of manifolds; 2) a systematic strategy to compare manifolds. Especially, we take Riemann curvature and sectional curvature as part of criterion, which can reflect the intrinsic geometric properties of manifolds. Some interesting results and phenomenon are given, which help in specifying the similarities and differences between the features extracted by two networks and demystifying the intrinsic mechanism of deep neural networks.


title: Frequently Occurring Reconnection Jets from Sunspot Light Bridges

abstract: Solid evidence of magnetic reconnection is rarely reported within sunspots, the darkest regions with the strongest magnetic fields and lowest temperatures in the solar atmosphere. Using the world's largest solar telescope, the 1.6-meter Goode Solar Telescope, we detect prevalent reconnection through frequently occurring fine-scale jets in the H${\alpha}$ line wings at light bridges, the bright lanes that may divide the dark sunspot core into multiple parts. Many jets have an inverted Y-shape, shown by models to be typical of reconnection in a unipolar field environment. Simultaneous spectral imaging data from the Interface Region Imaging Spectrograph show that the reconnection drives bidirectional flows up to 200~km~s$^{-1}$, and that the weakly ionized plasma is heated by at least an order of magnitude up to $\sim$80,000 K. Such highly dynamic reconnection jets and efficient heating should be properly accounted for in future modeling efforts of sunspots. Our observations also reveal that the surge-like activity previously reported above light bridges in some chromospheric passbands such as the H${\alpha}$ core has two components: the ever-present short surges likely to be related to the upward leakage of magnetoacoustic waves from the photosphere, and the occasionally occurring long and fast surges that are obviously caused by the intermittent reconnection jets.


title: Nonlinear operations on a class of modulation spaces

abstract: We discuss when the nonlinear operation $f\mapsto F(f)$ maps the modulation space $M^{p,q}_s(\mathbb{R}^n)$ ($1 \leq p,q \leq \infty$) to the same space again. It is known that $M^{p,q}_s(\mathbb{R}^n)$ is a multiplication algebra when $s > n-n/q$, hence it is true for this space if $F$ is entire. We claim that it is still true for non-analytic $F$ when $q\geq4/3$.


title: On Taylor coefficients of smooth functions

abstract: We study the Borel map, which maps infinitely differentiable functions on an interval to the jets of their Taylor coefficients at a given point in the interval. Our main results include a complete description of the image of the Borel map for Beurling classes of smooth functions and a moment-type summation method which allows one to recover a function from its Taylor jet. A surprising feature of this description is an unexpected threshold at the logarithmic class. Another interesting finding is a "duality" between non-quasianalytic and quasianalytic classes, which reduces the description of the image of the Borel map for non-quasianalytic classes to the one for the corresponding quasianalytic classes, and complements classical results of Carleson and Ehrenpreis.


title: Decoupled Learning for Factorial Marked Temporal Point Processes

abstract: This paper introduces the factorial marked temporal point process model and presents efficient learning methods. In conventional (multi-dimensional) marked temporal point process models, event is often encoded by a single discrete variable i.e. a marker. In this paper, we describe the factorial marked point processes whereby time-stamped event is factored into multiple markers. Accordingly the size of the infectivity matrix modeling the effect between pairwise markers is in power order w.r.t. the number of the discrete marker space. We propose a decoupled learning method with two learning procedures: i) directly solving the model based on two techniques: Alternating Direction Method of Multipliers and Fast Iterative Shrinkage-Thresholding Algorithm; ii) involving a reformulation that transforms the original problem into a Logistic Regression model for more efficient learning. Moreover, a sparse group regularizer is added to identify the key profile features and event labels. Empirical results on real world datasets demonstrate the efficiency of our decoupled and reformulated method. The source code is available online.


title: NNLL resummation for the associated production of a top pair with a heavy boson at the LHC

abstract: In this presentation we review recent results on the resummation of soft gluon emission corrections for the associated production of a top-quark pair with a heavy boson (Higgs/W/Z) at the Large Hadron Collider (LHC). We develop a parton level Monte Carlo based on a soft-gluon resummation formula valid up to next-to-next-to-leading logarithmic (NNLL) accuracy. With this tool we study the impact of the newly computed corrections to the total cross sections and some important differential distributions.


title: Embedding Learning Through Multilingual Concept Induction

abstract: We present a new method for estimating vector space representations of words: embedding learning by concept induction. We test this method on a highly parallel corpus and learn semantic representations of words in 1259 different languages in a single common space. An extensive experimental evaluation on crosslingual word similarity and sentiment analysis indicates that concept-based multilingual embedding learning performs better than previous approaches.


title: Dual-wavelength vortex beam with high stability in diode-pumped Yb:CaGdAlO4 Laser

abstract: We present a stable dual-wavelength vortex beam carrying orbital angular momentum (OAM) with two spectral peaks separated by a few terahertz in diode-pumped Yb:CaGdAlO4 (CALGO) laser. The dual-wavelength spectrum is controlled by the pump power and off-axis loss in laser resonator, arising from the broad emission bandwidth of Yb:CALGO. The OAM beam is obtained by a pair of cylindrical lens that serves as an {\pi}/2 convertor for the high-order Hermite-Gaussian modes. The stability is verifed that the 1 h OAM beam with two spectral peaks at 1046.1 nm and 1057.2 nm (3.01 THz interval) can steadily operate for more than three hours. It has great potential for scaling the application for OAM beams in Terahertz spectroscopy, high-resolution interferometry, and so on.


title: A new nonparametric test for two sample multivariate location problem with application to astronomy

abstract: This paper provides a nonparametric test for the identity of two multivariate continuous distribution functions (d.f.'s) when they differ in locations. The test uses Wilcoxon rank-sum statistics on distances between observations for each of the components and is unaffected by outliers. It is numerically compared with two existing procedures in terms of power. The simulation study shows that its power is strictly increasing in the sample sizes and/or in the number of components. The applicability of this test is demonstrated by use of two astronomical data sets on early-type galaxies.


title: Guidelines for Systematic Mapping Studies in Security Engineering

abstract: Security engineering in the software lifecycle aims at protecting information and systems to guarantee confidentiality, integrity, and availability. As security engineering matures and the number of research papers grows, there is an increasing need for papers that summarize results and provide an overview of the area. A systematic mapping study "maps" a research area by classifying papers to identify which topics are well-studied and which need additional study. Therefore, systematic mapping studies are becoming increasingly important in security engineering. This chapter provides methodological support for systematic mapping studies in security engineering based on examples from published security engineering papers. Because security engineering is similar to software engineering in that it bridges research and practice, researchers can use the same basic systematic mapping process, as follows: (1) study planning, (2) searching for studies, (3) study selection, (4) study quality assessment, (5) data extraction, (6) data classification, (7) data analysis, and (8) reporting of results. We use published mapping studies to describe the tailoring of this process for security engineering. In addition to guidance on how to perform systematic mapping studies in security engineering, this chapter should increase awareness in the security engineering community of the need for additional mapping studies.


title: Accelerating AdS black holes as the holographic heat engines in a benchmarking scheme

abstract: We investigate the properties of holographic heat engines with an uncharged accelerating non-rotating AdS black hole as the working substance in a benchmarking scheme. We find that the efficiencies of the black hole heat engines can be influenced by both the size of the benchmark circular cycle and the cosmic string tension as a thermodynamic variable. In general, the efficiency can be increased by enlarging the cycle, but is still constrained by a universal bound $2\pi/(\pi+4)$ as expected. A cross-comparison of the efficiencies of the accelerating black hole heat engines and Schwarzschild-AdS black hole heat engines suggests that the acceleration also increases the efficiency although the amount of increase is not remarkable.


title: Recent Results on Classifying Risk-Based Testing Approaches

abstract: In order to optimize the usage of testing efforts and to assess risks of software-based systems, risk-based testing uses risk (re-)assessments to steer all phases in a test process. Several risk-based testing approaches have been proposed in academia and/or applied in industry, so that the determination of principal concepts and methods in risk-based testing is needed to enable a comparison of the weaknesses and strengths of different risk-based testing approaches. In this chapter we provide an (updated) taxonomy of risk-based testing aligned with risk considerations in all phases of a test process. It consists of three top-level classes, i.e., contextual setup, risk assessment, and risk-based test strategy. This taxonomy provides a framework to understand, categorize, assess and compare risk-based testing approaches to support their selection and tailoring for specific purposes. Furthermore, we position four recent risk-based testing approaches into the taxonomy in order to demonstrate its application and alignment with available risk-based testing approaches.


title: Lower bounds on the growth of Sobolev norms in some linear time dependent Schr\"odinger equations

abstract: In this paper we consider linear, time dependent Schr\"odinger equations of the form $i \partial_t \psi = K_0 \psi + V(t) \psi $, where $K_0$ is a positive self-adjoint operator with discrete spectrum and whose spectral gaps are asymptotically constant. We give a strategy to construct bounded perturbations $V(t)$ such that the Hamiltonian $K_0 + V(t)$ generates unbounded orbits. We apply our abstract construction to three cases: (i) the Harmonic oscillator on $\mathbb R$, (ii) the half-wave equation on $\mathbb T$ and (iii) the Dirac-Schr\"odinger equation on the sphere. In each case, $V(t)$ is a smooth and periodic in time pseudodifferential operator and the Schr\"odinger equation has solutions fulfilling $\| \psi(t) \|_r \gtrsim |t|^{r }$ as $|t| \gg 1$.


title: Curriculum Guidelines for Undergraduate Programs in Data Science

abstract: The Park City Math Institute (PCMI) 2016 Summer Undergraduate Faculty Program met for the purpose of composing guidelines for undergraduate programs in Data Science. The group consisted of 25 undergraduate faculty from a variety of institutions in the U.S., primarily from the disciplines of mathematics, statistics and computer science. These guidelines are meant to provide some structure for institutions planning for or revising a major in Data Science.


title: Combinatorial proofs and generalizations of conjectures related to Euler's partition theorem

abstract: In a recent work, Andrews gave analytic proofs of two conjectures concerning some variations of two combinatorial identities between partitions of a positive integer into odd parts and partitions into distinct parts discovered by Beck. Subsequently, using the same method as Andrews, Chern presented the analytic proof of another Beck's conjecture relating the gap-free partitions and distinct partitions with odd length. However, the combinatorial interpretations of these conjectures are still unclear and required. In this paper, motivated by Glaisher's bijection, we give the combinatorial proofs of these three conjectures directly or by proving more generalized results.


title: Preferential Attachment Graphs with Planted Communities

abstract: A variation of the preferential attachment random graph model of Barab\'asi and Albert is defined that incorporates planted communities. The graph is built progressively, with new vertices attaching to the existing ones one-by-one. At every step, the incoming vertex is randomly assigned a label, which represents a community it belongs to. This vertex then chooses certain vertices as its neighbors, with the choice of each vertex being proportional to the degree of the vertex multiplied by an affinity depending on the labels of the new vertex and a potential neighbor. It is shown that the fraction of half-edges attached to vertices with a given label converges almost surely for some classes of affinity matrices. In addition, the empirical degree distribution for the set of vertices with a given label converges to a heavy tailed distribution, such that the tail decay parameter can be different for different communities. Our proof method may be of independent interest, both for the classical Barab\'asi -Albert model and for other possible extensions.


title: Thermally activated vapor bubble nucleation: the Landau-Lifshitz/Van der Waals approach

abstract: Vapor bubbles are formed in liquids by two mechanisms: evaporation (temperature above the boiling threshold) and cavitation (pressure below the vapor pressure). The liquid resists in these metastable (overheating and tensile, respectively) states for a long time since bubble nucleation is an activated process that needs to surmount the free energy barrier separating the liquid and the vapor states. The bubble nucleation rate is difficult to assess and, typically, only for extremely small systems treated at atomistic level of detail. In this work a powerful approach, based on a continuum diffuse interface modeling of the two-phase fluid embedded with thermal fluctuations (Fluctuating Hydrodynamics) is exploited to study the nucleation process in homogeneous conditions, evaluating the bubble nucleation rates and following the long term dynamics of the metastable system, up to the bubble coalescence and expansion stages. In comparison with more classical approaches, this methodology allows on the one hand to deal with much larger systems observed for a much longer times than possible with even the most advanced atomistic models. On the other it extends contin- uum formulations to thermally activated processes, impossible to deal with in a purely determinist setting.


title: Community Recovery in a Preferential Attachment Graph

abstract: A message passing algorithm is derived for recovering communities within a graph generated by a variation of the Barab\'{a}si-Albert preferential attachment model. The estimator is assumed to know the arrival times, or order of attachment, of the vertices. The derivation of the algorithm is based on belief propagation under an independence assumption. Two precursors to the message passing algorithm are analyzed: the first is a degree thresholding (DT) algorithm and the second is an algorithm based on the arrival times of the children (C) of a given vertex, where the children of a given vertex are the vertices that attached to it. Comparison of the performance of the algorithms shows it is beneficial to know the arrival times, not just the number, of the children. The probability of correct classification of a vertex is asymptotically determined by the fraction of vertices arriving before it. Two extensions of Algorithm C are given: the first is based on joint likelihood of the children of a fixed set of vertices; it can sometimes be used to seed the message passing algorithm. The second is the message passing algorithm. Simulation results are given.


title: A Next-Best-Smell Approach for Remote Gas Detection with a Mobile Robot

abstract: The problem of gas detection is relevant to many real-world applications, such as leak detection in industrial settings and landfill monitoring. Using mobile robots for gas detection has several advantages and can reduce danger for humans. In our work, we address the problem of planning a path for a mobile robotic platform equipped with a remote gas sensor, which minimizes the time to detect all gas sources in a given environment. We cast this problem as a coverage planning problem by defining a basic sensing operation -- a scan with the remote gas sensor -- as the field of "view" of the sensor. Given the computing effort required by previously proposed offline approaches, in this paper we suggest a online coverage algorithm, called Next-Best-Smell, adapted from the Next-Best-View class of exploration algorithms. Our algorithm evaluates candidate locations with a global utility function, which combines utility values for travel distance, information gain, and sensing time, using Multi-Criteria Decision Making. In our experiments, conducted both in simulation and with a real robot, we found the performance of the Next-Best-Smell approach to be comparable with that of the state-of-the-art offline algorithm, at much lower computational cost.


title: Fluctuations, renormalizations, and a convective instability in driven wet active matter

abstract: A scalar model of wet active matter in the presence of an imposed temperature gradient, or chemical potential gradient, is considered. It is shown that there is a convective instability driven by a (negative) activity parameter. In this non-equilibrium steady state the generic long-ranged correlations are computed and compared and contrasted with the analogous results in a passive fluid. In addition, the non-equilibrium Casimir pressure or force is computed. Singularities in various physical quantities as the instability is approached are determined. Finally, we give the generalized Lorenz equations characterizing the fluid behavior above the instability and contrast these equations to the Lorenz equations for the Rayleigh-Bernard instability in a passive fluid.


title: Charge transfer in H + He$^+$ and H$^+$ + He collisions in excited states

abstract: We present a theoretical study of charge transfer in collisions of excited ($n=2,3$) hydrogen atoms with He$^+$ and in collisions of excited ($n=2,3$) helium atoms with H$^+$, extending the results of Phys. Rev. A 82 012708 (2010). A combination of quantum-mechanical and semi-classical approaches is employed to calculate the charge-exchange cross sections at collision energies from 0.1 eV/u up to 1 keV/u. These methods are based on accurate ab initio potential energy curves and non-adiabatic couplings for the molecular ion HeH$^+$. Charge transfer can occur either in singlet or in triplet states, and the differences between the singlet and triplet spin manifolds are discussed. The dependence of the cross section on the quantum numbers $n$ and $l$ of the initial state is demonstrated. The isotope effect on the charge transfer cross sections, arising at low collision energy when H is substituted by D or T, is investigated. Finally, the impact of the present calculations on models of laboratory plasmas is discussed.


title: ERIM: Secure, Efficient In-process Isolation with Memory Protection Keys (MPK)

abstract: Isolating sensitive state and data can increase the security and robustness of many applications. Examples include protecting cryptographic keys against exploits like OpenSSL's Heartbleed bug or protecting a language runtime from native libraries written in unsafe languages. When runtime references across isolation boundaries occur relatively infrequently, then conventional page-based hardware isolation can be used, because the cost of kernel- or hypervisor-mediated domain switching is tolerable. However, some applications, such as the isolation of cryptographic session keys in network-facing services, require very frequent domain switching. In such applications, the overhead of kernel- or hypervisor-mediated domain switching is prohibitive. In this paper, we present ERIM, a novel technique that provides hardware-enforced isolation with low overhead on x86 CPUs, even at high switching rates (ERIM's measured overhead is less than 1% for 100,000 switches per second). The key idea is to combine protection keys (MPKs), a feature recently added to x86 that allows protection domain switches in userspace, with binary inspection to prevent circumvention. We show that ERIM can be applied with little effort to new and existing applications, doesn't require compiler changes, can run on a stock Linux kernel, and has low runtime overhead even at high domain switching rates.


title: Charm production in interactions of antiproton with proton and nuclei at ${\bar P}ANDA$ energies

abstract: We study the production of charmed baryons in the antiproton-proton and antiproton-nucleus interactions within a fully covariant model that is based on an effective Lagrangian approach. The baryon production proceeds via the $t$-channel $D^0$ and $D^{*0}$ meson-exchange diagrams. We have also explored the production of the charm-baryon hypernucleus $^{16}_{\Lambda_c^+}$O in the antiproton - $^{16}$O collisions. For antiproton beam momenta of interest to the ${\bar P}ANDA$ experiment, the 0$^\circ$ differential cross sections for the formation of $^{16}_{\Lambda_c^+}$O hypernuclear states with simple particle-hole configurations, have magnitudes in the range of a few $\mu$b/sr.


title: On generalized choice and coloring numbers

abstract: A well-known result of Alon shows that the coloring number of a graph is bounded by a function of its choosability. We explore this relationship in a more general setting with relaxed assumptions on color classes, encoded by a graph parameter.


title: Composite Behavioral Modeling for Identity Theft Detection in Online Social Networks

abstract: In this work, we aim at building a bridge from poor behavioral data to an effective, quick-response, and robust behavior model for online identity theft detection. We concentrate on this issue in online social networks (OSNs) where users usually have composite behavioral records, consisting of multi-dimensional low-quality data, e.g., offline check-ins and online user generated content (UGC). As an insightful result, we find that there is a complementary effect among different dimensions of records for modeling users' behavioral patterns. To deeply exploit such a complementary effect, we propose a joint model to capture both online and offline features of a user's composite behavior. We evaluate the proposed joint model by comparing with some typical models on two real-world datasets: Foursquare and Yelp. In the widely-used setting of theft simulation (simulating thefts via behavioral replacement), the experimental results show that our model outperforms the existing ones, with the AUC values $0.956$ in Foursquare and $0.947$ in Yelp, respectively. Particularly, the recall (True Positive Rate) can reach up to $65.3\%$ in Foursquare and $72.2\%$ in Yelp with the corresponding disturbance rate (False Positive Rate) below $1\%$. It is worth mentioning that these performances can be achieved by examining only one composite behavior (visiting a place and posting a tip online simultaneously) per authentication, which guarantees the low response latency of our method. This study would give the cybersecurity community new insights into whether and how a real-time online identity authentication can be improved via modeling users' composite behavioral patterns.


title: Direction of cascades in a magnetofluid model with electron skin depth and ion sound Larmor radius scales

abstract: The direction of cascades in a two-dimensional model that takes electron inertia and ion sound Larmor radius into account is studied, resulting in analytical expressions for the absolute equilibrium states of the energy and helicities. It is found that typically both the energy and magnetic helicity at scales shorter than electron skin depth have direct cascade, while at large scales the helicity has an inverse cascade as established earlier for reduced magnetohydrodynamics (MHD). It is also found that the introduction of gyro-effects allows for the existence of negative temperature (conjugate to energy) states and the condensation of energy to the large scales. Comparisons between two- and three-dimensional extended MHD models (MHD with two-fluid effects) show qualitative agreement between the two.


title: Artificial Impostors for Location Privacy Preservation

abstract: The progress of location-based services has led to serious concerns on location privacy leakage. For effective and efficient location privacy preservation (LPP), existing methods are still not fully competent. They are often vulnerable under the identification attack with side information, or hard to be implemented due to the high computational complexity. In this paper, we pursue the high protection efficacy and low computational complexity simultaneously. We propose a scalable LPP method based on the paradigm of counterfeiting locations. To make fake locations extremely plausible, we forge them through synthesizing artificial impostors (AIs). The AIs refer to the synthesized traces which have similar semantic features to the actual traces, and do not contain any target location. Two dedicated techniques are devised: the sampling-based synthesis method and population-level semantic model. They play significant roles in two critical steps of synthesizing AIs. We conduct experiments on real datasets in two cities (Shanghai, China and Asturias, Spain) to validate the high efficacy and scalability of the proposed method. In these two datasets, the experimental results show that our method achieves the preservation efficacy of $97.65\%$ and $96.12\%$, and its run time of building the generators is only $230.47$ and $215.92$ seconds, respectively. This study would give the research community new insights into improving the practicality of the state-of-the-art LPP paradigm via counterfeiting locations.


title: Channel Input Adaptation via Natural Type Selection

abstract: For the model of communication through a discrete memoryless channel using i.i.d. random block codes, where the channel is changing slowly from block to block, we propose a stochastic algorithm for adaptation of the generating distribution of the code in the process of continuous reliable communication. The purpose of the algorithm is to match the generating distribution $Q(x)$ to the changing channel $P(y\,|\,x)$, so that reliable communication is maintained at some constant rate $R$. This is achieved by a feedback of one bit per transmitted block. The feedback bit is determined by the joint type of the last transmitted codeword and the received block, a constant threshold $T>R$, and some conditional distribution $\Phi(x\,|\,y)$. Depending on the value of the feedback bit, the system parameters $Q(x)$ and $\Phi(x\,|\,y)$ are both updated according to the joint type of the last transmitted and received blocks, or remain unchanged. We show that, under certain technical conditions, the iterations of the algorithm lead to a distribution $Q(x)$, which guarantees reliable communication for all rates below the threshold $T$, provided that the discrete memoryless channel capacity of $P(y\,|\,x)$ stays above $T$.


title: Doppler magnification in flux-limited galaxy number counts with finite redshift bin width

abstract: In this paper, I investigated the Doppler magnification effect in the flux-limited galaxy number counts with finite redshift bin width. In contradiction with our intuition, the number counts correction formula for different redshift bin are different. For the redshift window as delta function and constant function, the correction formulas have been derived analytically in this paper. These two windows correspond to two extreme cases(i.e. extremely narrow bin and extremely broad bin). An simulation have been implemented to test these formulas. The simulation results indicate the perturbation changes smoothly from one extreme case to another while extending the redshift bin width from half bin width $0.05$ to $0.4$. As a result, Doppler magnification caused number counts perturbation of finite redshift bin width can NOT be compute via redshift integration of the perturbation formula derived with delta function redshift window. It is only a good approximation when the bin width is small. These results are important for galaxy redshift survey kinematic dipole estimation, and meaningful for the relativistic galaxy number counts angular power spectrum estimation, such as CLASSgal and CAMB source.


title: Neural Multi-task Learning in Automated Assessment

abstract: Grammatical error detection and automated essay scoring are two tasks in the area of automated assessment. Traditionally these tasks have been treated independently with different machine learning models and features used for each task. In this paper, we develop a multi-task neural network model that jointly optimises for both tasks, and in particular we show that neural automated essay scoring can be significantly improved. We show that while the essay score provides little evidence to inform grammatical error detection, the essay score is highly influenced by error detection.


title: Dense Recurrent Neural Networks for Scene Labeling

abstract: Recently recurrent neural networks (RNNs) have demonstrated the ability to improve scene labeling through capturing long-range dependencies among image units. In this paper, we propose dense RNNs for scene labeling by exploring various long-range semantic dependencies among image units. In comparison with existing RNN based approaches, our dense RNNs are able to capture richer contextual dependencies for each image unit via dense connections between each pair of image units, which significantly enhances their discriminative power. Besides, to select relevant and meanwhile restrain irrelevant dependencies for each unit from dense connections, we introduce an attention model into dense RNNs. The attention model enables automatically assigning more importance to helpful dependencies while less weight to unconcerned dependencies. Integrating with convolutional neural networks (CNNs), our method achieves state-of-the-art performances on the PASCAL Context, MIT ADE20K and SiftFlow benchmarks.


title: The longitudinal leading-twist distribution amplitude of $J/\psi$ meson within background field theory

abstract: We make a detailed study on the $J/\psi$ meson longitudinal leading-twist distribution amplitude $\phi_{2;J/\psi}^\|$ by using the QCD sum rules within the background field theory. By keeping all the non-perturbative condensates up to dimension-six, we obtain accurate QCD sum rules for the moments $\langle\xi_{n;J/\psi}^\|\rangle$. The first three ones are $\langle\xi_{2;J/\psi}^\|\rangle=0.083(12)$, $\langle\xi_{4;J/\psi}^\|\rangle=0.015(5)$ and $\langle\xi_{6;J/\psi}^\|\rangle=0.003(2)$, leading to a single peaked behavior for $\phi_{2;J/\psi}^\|$ which is sharper than the previous ones around the region of $x\sim0.5$. As an application, we adopt the QCD light-cone sum rules to calculate the $B_c$ meson semileptonic decay $B_c^+ \to J/\psi \ell^+ \nu_\ell$. We obtain $\Gamma(B_c^+ \to J/\psi \ell^+ \nu_\ell) = (89.67^{+24.76}_{-19.06}) \times 10^{-15}~{\rm GeV}$ and $\Re(J/\psi \ell^+ \nu_\ell) = 0.217^{+0.069}_{-0.057}$, which agree with the next-to-leading order pQCD prediction and the new CDF measurement within errors.


title: Hurricane's Maximum Potential Intensity and the Gravitational Power of Precipitation

abstract: The concept of Maximum Potential Intensity (MPI) is widely used in tropical cyclone research to estimate the minimum central pressure and the maximum velocity of tropical storms from environmental parameters. The MPI pressure derives from consideration of an idealized thermodynamic cycle, while the MPI velocity additionally requires information about real-time power and heat flows within the storm. Recently MPI velocity was proposed to be a substantial overestimate (by 10-30 percent) presumably neglecting the power needed to lift precipitating water (the gravitational power of precipitation). This conclusion did not involve a theoretical analysis of the MPI concept but was based on observed hurricane rainfall to estimate gravitational power of precipitation. However, since the MPI pressure estimate does explicitly account for lifting water, and the MPI velocity derives from this pressure, the question arises whether a correction to MPI velocity is required. Here we represent and examine the MPI derivations in their most general form and show that although a correction to MPI velocity is justified, it is an order of magnitude or so smaller than originally proposed. We show that the neglect of gravitational power of precipitation in the MPI velocity estimate was caused by the incomplete formulation of the general relationship between pressure work and dissipation of kinetic energy, taken per unit time and integrated over the storm. We highlight the importance of an internally consistent framework to estimate both storm energy and power and provide some perspectives for further investigating the role of moisture.


title: $R^4 $ terms in supergravities via T-duality constraint

abstract: It has been speculated in the literature that the effective actions of string theories at any order of $\alpha'$ should be invariant under the Buscher rules plus their higher covariant derivative corrections. This may be used as a constraint to find effective actions at any order of $\alpha'$, in particular, the metric, the B-field and the dilaton couplings in supergravities at order $\alpha'^3$ up to an overall factor. For the simple case of zero B-field and diagonal metric in which we have done the calculations explicitly, we have found that the constraint fixes almost all the seven independent Riemann curvature couplings. There is only one term which is not fixed, because when metric is diagonal, the reduction of two $R^4$ terms become identical. The Riemann curvature couplings that the T-duality constraint produces for both type II and Heterotic theories are fully consistent with the existing couplings in the literature which have been found by the S-matrix and by the sigma-model approaches.


title: Distributed Multi-User Wireless Charging Power Allocation

abstract: Wireless power charging enables portable devices to be permanently unplugged. Due to its low transmission power and low transmission efficiency, it requires much longer time slot to charge users compared with that for data transmission in wireless communication networks. Besides, each user's demand urgency needs to be taken into consideration for power allocation. Therefore, new algorithms are essential for wireless power allocation in multi-user wireless charging networks. In this paper, this problem is formulated as a static noncooperative game. It is shown that there exists a unique Nash equilibrium, which is the static state of the wireless power charging network. A distributed power allocation algorithm is proposed to compute the Nash equilibrium of the game. The main result of the paper consists of rigorous analysis of the distributed algorithm for power allocation. The algorithm is shown to converge to Nash equilibrium of the game with exponentially convergence rate for arbitrary initial value with synchronous scheduling. Moreover, the distributed algorithm is also convergence guaranteed with asynchronous scheduling under communication delay and packet drops. Numerical simulations prove the correctness of the analysis and demonstrate the fast convergence of the algorithm and the robustness to synchronous scheduling.


title: Peak-locking centroid bias in Shack-Hartmann wavefront sensing

abstract: Shack-Hartmann wavefront sensing relies on accurate spot centre measurement. Several algorithms were developed with this aim, mostly focused on precision, i.e. minimizing random errors. In the solar and extended scene community, the importance of the accuracy (bias error due to peak-locking, quantisation or sampling) of the centroid determination was identified and solutions proposed. But these solutions only allow partial bias corrections. To date, no systematic study of the bias error was conducted. This article bridges the gap by quantifying the bias error for different correlation peak-finding algorithms and types of sub-aperture images and by proposing a practical solution to minimize its effects. Four classes of sub-aperture images (point source, elongated laser guide star, crowded field and solar extended scene) together with five types of peak-finding algorithms (1D parabola, the centre of gravity, Gaussian, 2D quadratic polynomial and pyramid) are considered, in a variety of signal-to-noise conditions. The best performing peak-finding algorithm depends on the sub-aperture image type, but none is satisfactory to both bias and random errors. A practical solution is proposed that relies on the anti-symmetric response of the bias to the sub-pixel position of the true centre. The solution decreases the bias by a factor of ~7 to values of < 0.02 pix. The computational cost is typically twice of current cross-correlation algorithms.


title: System-of-Systems Viewpoint for System Architecture Documentation

abstract: Context: The systems comprising a system of systems (SoS) are inde- pendently acquired, operated, and managed. Frequently, the architecture documentation of these existing systems addresses only a stand-alone perspective, and must be augmented to address concerns that arise in the integrated SoS. Objective: We evaluated an architecture documentation viewpoint to address the concerns of a SoS architect about a constituent system, to support SoS design and analysis involving that constituent system. Method: We performed an expert review of documentation produced by applying the viewpoint to a system, using the active review method. Results: The expert panel was able to used a view constructed using the baseline version of the viewpoint to answer questions related to all SoS architect concerns about a constituent system, except for questions concerning the interaction of the constituent system with the platform and network infrastructure. Conclusions: We found that the expert panel was unable to answer certain questions because the baseline version of the viewpoint had a gap in coverage related to relationship of software units of execution (e.g., processes or services) to computers and networks. The viewpoint was revised to add a Deployment Model to address these concerns, and is included in an appendix.


title: Global and Concrete Quantizations on General Type I Groups

abstract: In recent papers and books, a global quantization has been developed for unimodular groups of type I. It involves operator-valued symbols defined on the product between the group $\mathsf{G}$ and its unitary dual $\widehat{\mathsf{G}}$, composed of equivalence classes of irreducible representations. For compact or for graded Lie groups, this has already been developed into a powerful pseudo-differential calculus. In the present article we extend the formalism to arbitrary locally compact groups of type I, making use of the Fourier theory of non-unimodular second countable groups. The unitary dual and its Plancherel measure being quite abstract in general, we put into evidence situations in which concrete forms are available. Kirillov theory and parametrizations of large parts of $\widehat{\mathsf{G}}$ allow rewriting the basic formulae in a manageable form. Some examples of completely solvable groups are worked out.


title: Direct observation of atom-ion non-equilibrium sympathetic cooling

abstract: Sympathetic cooling is the process of energy exchange between a system and a colder bath. We investigate this fundamental process in an atom-ion experiment where the system is composed of a single ion, trapped in a radio-frequency Paul trap, and prepared in a coherent state of ~200 K and the bath is an ultracold cloud of atoms at {\mu}K temperature. We directly observe the sympathetic cooling dynamics with single-shot energy measurements during one, to several, collisions in two distinct regimes. In one, collisions predominantly cool the system with very efficient momentum transfer leading to cooling in only a few collisions. In the other, collisions can both cool and heat the system due to the non-equilibrium dynamics of the atom-ion collisions in the presence of the ion-trap's oscillating electric fields. While the bulk of our observations agree well with a molecular dynamics simulation of hard-sphere (Langevin) collisions, a measurement of the scattering angle distribution reveals forward-scattering (glancing) collisions which are beyond the Langevin model. This work paves the way for further non-equilibrium and collision dynamics studies using the well-controlled atom-ion system.


title: Information loss in effective field theory: entanglement and thermal entropies

abstract: Integrating out high energy degrees of freedom to yield a low energy effective field theory leads to a loss of information with a concomitant increase in entropy. We obtain the effective field theory of a light scalar field interacting with heavy fields after tracing out the heavy degrees of freedom from the time evolved density matrix. The initial density matrix describes the light field in its ground state and the heavy fields in equilibrium at a common temperature $T$. For $T=0$, we obtain the reduced density matrix in a perturbative expansion, it reveals an emergent mixed state as a consequence of the entanglement between light and heavy fields. We obtain the effective action that determines the time evolution of the \emph{reduced} density matrix for the light field in a non-perturbative Dyson resummation of one-loop correlations of the heavy fields. The Von-Neumann \emph{entanglement entropy} associated with the reduced density matrix is obtained for the non-resonant and resonant cases in the asymptotic long time limit. In the non-resonant case the reduced density matrix displays an \emph{incipient} thermalization albeit with a wave-vector, time and coupling dependent \emph{effective temperature} as a consequence of memory of initial conditions. The entanglement entropy is time independent and is the \emph{thermal entropy} for this effective, non-equilibrium temperature. In the resonant case the light field fully \emph{thermalizes} with the heavy fields, the reduced density matrix looses memory of the initial conditions and the entanglement entropy becomes the \emph{thermal entropy} of the light field. We discuss the relation between the entanglement entropy ultraviolet divergences and renormalization.

categories: hep-th cond-mat.stat-mech hep-ph quant-ph


title: UAV-Enabled Cooperative Jamming for Improving Secrecy of Ground Wiretap Channel

abstract: This letter proposes a novel UAV-enabled mobile jamming scheme to improve the secrecy rate of ground wiretap channel. Specifically, a UAV is employed to transmit jamming signals to combat against eavesdropping. Such a mobile jamming scheme is particularly appealing since the UAV-enabled jammer can fly close to the eavesdropper and opportunistically jam it by leveraging the UAV's mobility. We aim to maximize the average secrecy rate by jointly optimizing the UAV's trajectory and jamming power over a given flight period. To make the problem more tractable, we drive a closed-form lower bound for the achievable secrecy rate, based on which the UAV's trajectory and transmit power are optimized alternately by an efficient iterative algorithm applying the block coordinate descent and successive convex optimization techniques. Simulation results demonstrate that the proposed joint design can significantly enhance the secrecy rate of the considered wiretap system as compared to benchmark schemes.


title: Shubnikov-de Haas-like Quantum Oscillations in Artificial One-Dimensional LaAlO3/SrTiO3 Electron Channels

abstract: The widely reported magnetoresistance oscillations in LaAlO3/SrTiO3 heterostructures have invariably been attributed to the Shubnikov-de Haas (SdH) effect, despite a pronounced inconsistency with low-field Hall resistance measurements. Here we report SdH-like resistance oscillations in quasi-1D electron waveguides created at the LaAlO3/SrTiO3 interface by conductive atomic force microscopy lithography. These oscillations can be directly attributed to magnetic depopulation of magnetoelectric subbands. Our results suggest that the SdH oscillations in 2D SrTiO3-based systems may originate from naturally forming quasi-1D channels.


title: On the modified method of simplest equation and the nonlinear Schr\"odinger equation

abstract: We consider an extension of the methodology of the modified method of simplest equation to the case of use of two simplest equations. The extended methodology is applied for obtaining exact solutions of model nonlinear partial differential equations for deep water waves: the nonlinear Schr\"odinger equation. It is shown that the methodology works also for other equations of the nonlinear Schr\"odinger kind.


title: Dynein catch bond as a mediator of codependent bidirectional cellular transport

abstract: Intracellular bidirectional transport of cargo on microtubule filaments is achieved by the collective action of oppositely directed dynein and kinesin motors. Experiments have found that in certain cases, inhibiting the activity of one type of motor results in an overall decline in the motility of the cellular cargo in both directions. This counter-intuitive observation, referred to as {\em paradox of codependence} is inconsistent with the existing paradigm of a mechanistic tug-of-war between oppositely directed motors. Unlike kinesin, dynein motors exhibit catchbonding, wherein the unbinding rates of these motors decrease with increasing force on them. Incorporating this catchbonding behavior of dynein in a theoretical model, we show that the functional divergence of the two motors species manifests itself as an internal regulatory mechanism, and leads to codependent transport behaviour in biologically relevant regimes. Using analytical methods and stochastic simulations, we analyse the processivity characteristics and probability distribution of run times and pause times of transported cellular cargoes. We show that catchbonding can drastically alter the transport characteristics and also provide a plausible resolution of the paradox of codependence.


title: Time series kernel similarities for predicting Paroxysmal Atrial Fibrillation from ECGs

abstract: We tackle the problem of classifying Electrocardiography (ECG) signals with the aim of predicting the onset of Paroxysmal Atrial Fibrillation (PAF). Atrial fibrillation is the most common type of arrhythmia, but in many cases PAF episodes are asymptomatic. Therefore, in order to help diagnosing PAF, it is important to design procedures for detecting and, more importantly, predicting PAF episodes. We propose a method for predicting PAF events whose first step consists of a feature extraction procedure that represents each ECG as a multi-variate time series. Successively, we design a classification framework based on kernel similarities for multi-variate time series, capable of handling missing data. We consider different approaches to perform classification in the original space of the multi-variate time series and in an embedding space, defined by the kernel similarity measure. We achieve a classification accuracy comparable with state of the art methods, with the additional advantage of detecting the PAF onset up to 15 minutes in advance.


title: Linking and Cutting Spanning Trees

abstract: We consider the problem of uniformly generating a spanning tree, of a connected undirected graph. This process is useful to compute statistics, namely for phylogenetic trees. We describe a Markov chain for producing these trees. For cycle graphs we prove that this approach significantly outperforms existing algorithms. For general graphs we obtain no analytical bounds, but experimental results show that the chain still converges quickly. This yields an efficient algorithm, also due to the use of proper fast data structures. To bound the mixing time of the chain we describe a coupling, which we analyse for cycle graphs and simulate for other graphs.


title: Monocular Imaging-based Autonomous Tracking for Low-cost Quad-rotor Design - TraQuad

abstract: TraQuad is an autonomous tracking quadcopter capable of tracking any moving (or static) object like cars, humans, other drones or any other object on-the-go. This article describes the applications and advantages of TraQuad and the reduction in cost (to about 250$) that has been achieved so far using the hardware and software capabilities and our custom algorithms wherever needed. This description is backed by strong data and the research analyses which have been drawn out of extant information or conducted on own when necessary. This also describes the development of completely autonomous (even GPS is optional) low-cost drone which can act as a major platform for further developments in automation, transportation, reconnaissance and more. We describe our ROS Gazebo simulator and our STATUS algorithms which form the core of our development of our object tracking drone for generic purposes.


title: Multi-soliton states under triangular spatial modulation of the quadratic nonlinearity

abstract: We introduce multi-soliton sets in the two-dimensional medium with the second-harmonic-generating nonlinearity subject to spatial modulation in the form of a triangle of singular peaks. Various families of symmetric and asymmetric sets are constructed, and their stability is investigated. Stable symmetric patterns may be built of 1, 4, or 7 individual solitons, while stable asymmetric ones contain 1, 2, or 3 solitons. Symmetric and asymmetric patterns may demonstrate mutual bistability. The shift of the asymmetric single-soliton state from the central position is accurately predicted analytically. Vortex rings composed of three solitons are produced too.


title: Multiplicity of singularities is not a bi-Lipschitz invariant

abstract: It was conjectured that multiplicity of a singularity is bi-Lipschitz invariant. We disprove this conjecture, constructing examples of bi-Lipschitz equivalent complex algebraic singularities with different values of multiplicity.


title: Mean field dynamics of a quantum tracer particle interacting with a boson gas

abstract: We consider the dynamics of a heavy quantum tracer particle coupled to a non-relativistic boson field in ${\mathbb R}^3$. The pair interactions of the bosons are of mean-field type, with coupling strength proportional to $\frac1N$ where $N$ is the expected particle number. Assuming that the mass of the tracer particle is proportional to $N$, we derive generalized Hartree equations in the limit $N\rightarrow\infty$. Moreover, we prove the global well-posedness of the associated Cauchy problem for sufficiently weak interaction potentials.


title: Next-to-leading order QCD and electroweak corrections to Higgs-strahlung processes at the LHC

abstract: In this paper we calculate the total and fiducial cross sections as well as differential distributions for the Higgs-strahlung or VH process $p p \to VH \to l\nu_l/l^- l^+ + H$, (V = W or Z, l=e,$\mu$) including QCD and electro-weak corrections up to next-to-leading order before and after reweighting photon PDFs of NNPDF2.3qed, NNPDF3.0qed, MRST2004qed, CT14QEDinc, and LUXqed at the LHC with 13 TeV and Higgs boson mass $\ M_{H}=125$ GeV. The predictions from the various photon PDFs before and after reweighting against each other are in good agreement. The photon PDF uncertainties of the photon-induced cross sections decrease significantly with the reweighting PDFs.


title: Interpolation of functional by integral continued C-fractions

abstract: The functional interpolation problem on a continual set of nodes by an integral continued C-fraction is studied. The necessary and sufficient conditions for its solvability are found. As a particular case, the considered integral continued fraction contains a standard interpolation continued C-fraction which is used to approximate the functions of one variable.


title: Basic Model of Purposeful Kinesis

abstract: The notions of taxis and kinesis are introduced and used to describe two types of behavior of an organism in non-uniform conditions: (i) Taxis means the guided movement to more favorable conditions; (ii) Kinesis is the non-directional change in space motion in response to the change of conditions. Migration and dispersal of animals has evolved under control of natural selection. In a simple formalisation, the strategy of dispersal should increase Darwinian fitness. We introduce new models of purposeful kinesis with diffusion coefficient dependent on fitness. The local and instant evaluation of Darwinian fitness is used, the reproduction coefficient. New models include one additional parameter, intensity of kinesis, and may be considered as the {\em minimal models of purposeful kinesis}. The properties of models are explored by a series of numerical experiments. It is demonstrated how kinesis could be beneficial for assimilation of patches of food or of periodic fluctuations. Kinesis based on local and instant estimations of fitness is not always beneficial: for species with the Allee effect it can delay invasion and spreading. It is proven that kinesis cannot modify stability of positive homogeneous steady states.


title: Outage Analysis of Energy Harvesting based Opportunistic Cooperative Communication Systems

abstract: Wireless energy harvesting constitutes an effective way to prolong the lifetime of wireless networks. In this paper, an opportunistic decode-and-forward cooperative communication system is investigated, where the energy-constrained relays harvest energy from the received information signal and the co-channel interference (CCI) signals and then use that harvested energy to forward the correctly decoded signal to the destination. Different from conventional relaying system with constant energy supplier, the transmission power of the energy-constrained relay depends on the available energy that harvested, which is a random process. Best relay selection that takes into account all the impacting ingredients on the received signal quality at the destination is deployed. The exact closed-form expression of the outage probability is derived, and the optimal value of the energy harvesting ratio for achieving minimum outage is numerically investigated. In addition, the impacts of the CCI signals on the system's outage and diversity performances are analyzed. It is shown that the proposed relaying scheme can achieve full diversity order equal to the number of relays without the need of fixed power supplier.


title: Combined Constraints on the Equation of State of Dense Neutron-Rich Matter from Terrestrial Experiments and Observations of Neutron Stars

abstract: Within the parameter space of equation of state (EOS) of dense neutron-rich matter limited by existing constraints mainly from terrestrial nuclear experiments, we investigate how the neutron star maximum mass $M_{\rm{max}}>2.01\pm0.04$ M$_\odot$, radius $10.62<R_{\rm{1.4}}< 12.83$ km and tidal deformability $\Lambda_{1.4}\leq800$ of canonical neutron stars all together constrain the EOS of dense neutron-rich nucleonic matter. While the 3-D parameter space of $K_{\rm{sym}}$ (curvature of nuclear symmetry energy), $J_{\rm{sym}}$ and $J_0$ (skewness of the symmetry energy and EOS of symmetric nuclear matter, respectively) are narrowed down significantly by the observational constraints, more data are needed to pin down the individual values of $K_{\rm{sym}}$, $J_{\rm{sym}}$ and $J_0$ with quantified uncertainties. The $J_0$ largely controls the maximum mass of neutron stars. While the EOS with $J_0=0$ is sufficiently stiff to support neutron stars as massive as 2.37 M$_{\odot}$, to support the hyperthetical ones as massive as 2.74 M$_{\odot}$ (composite mass of GW170817) requires $J_0$ to be larger than its currently known maximum value of about 400 MeV and beyond the causality limit. The upper limit on the tidal deformability of $\Lambda_{1.4}=800$ from the recent observation of GW170817 is found to provide upper limits on some EOS parameters consistent with but far less restrictive than the existing constraints of other observables studied.

categories: nucl-th astro-ph.HE astro-ph.SR hep-ph nucl-ex


title: Time-Delay Origins of Fundamental Tradeoffs Between Risk of Large Fluctuations and Network Connectivity

abstract: For the class of noisy time-delay linear consensus networks, we obtain explicit formulas for risk of large fluctuations of a scalar observable as a function of Laplacian spectrum and its eigenvectors. It is shown that there is an intrinsic tradeoff between risk and effective resistance of the underlying coupling graph of the network. The main implication is that increasing network connectivity, increases the risk of large fluctuations. For vector-valued observables, we obtain computationally tractable lower and upper bounds for joint risk measures. Then, we study behavior of risk measures for networks with specific graph topologies and show how risk scales with network size.


title: Bulk Heating Effects as Tests for Collapse Models

abstract: We discuss limits on the noise strength parameter in mass-proportional-coupled wave function collapse models implied by bulk heating effects, and examine the role of the noise power spectrum in comparing experiments of different types. This comparison utilizes a calculation of the rate of heating through phonon excitation implied by a general noise power spectrum $\lambda(\omega)$. We find that in the standard heating formula, the reduction rate $\lambda$ is replaced by $\lambda_{\rm eff}=\frac{2}{3 \pi^{3/2}} \int d^3w e^{-\vec w^2} \vec w^2 \lambda(\omega_L(\vec w/r_c))$, with $\omega_L(\vec q)$ the longitudinal acoustic phonon frequency as a function of wave number $\vec q$, and with $r_C$ the noise correlation length. Hence if the noise power spectrum is cut off below $\omega_L(|\vec q| \sim r_c^{-1})$, the bulk heating rate is sharply reduced, allowing compatibility of current experimental results. We suggest possible new bulk heating experiments that can be performed subject to limits placed by natural heating from radioactivity and cosmic rays. The proposed experiments exploit the vanishing of thermal transport in the low temperature limit.


title: Casimir squared correction to the standard rotator Hamiltonian for the O($n$) sigma-model in the delta-regime

abstract: In a previous paper we found that the isospin susceptibility of the O($n$) sigma-model calculated in the standard rotator approximation differs from the next-to-next to leading order chiral perturbation theory result in terms vanishing like $1/\ell\,,$ for $\ell=L_t/L\to\infty$ and further showed that this deviation could be described by a correction to the rotator spectrum proportional to the square of the quadratic Casimir invariant. Here we confront this expectation with analytic nonperturbative results on the spectrum in 2 dimensions, by Balog and Heged\"us for $n=3,4$ and by Gromov, Kazakov and Vieira for $n=4$. We also consider the case of 3 dimensions.


title: Magnetic Imprisonment of Dusty Pinballs by a Supernova Remnant

abstract: Motivated by recent measurements of deposits of $^{60}$Fe on the ocean floor and the lunar surface, we model the transport of dust grains containing $^{60}$Fe from a near-Earth (i.e., within 100 pc) supernova (SN). We inject dust grains into the environment of an SN remnant (SNR) and trace their trajectories by applying a 1D hydrodynamic description assuming spherical symmetry to describe the plasma dynamics, and include a rudimentary, 3D magnetic field description to examine its influence on charged dust grains. We assume the interstellar medium (ISM) magnetic fields are turbulent, and are amplified by the SNR shock, while the SN wind and ejecta fields are negligible. We examine the various influences on the dust grains within the SNR to determine when/if the dust decouples from the plasma, how much it is sputtered, and where within the SNR the dust grains are located. We find that Rayleigh-Taylor instabilities are important for dust survival, as they influence the location of the SN's reverse shock. We find that the presence of a magnetic field within the shocked ISM material limits the passage of SN dust grains, with the field either reflecting or trapping the grains within the heart of the SNR. These results have important implications for {\it in situ} $^{60}$Fe measurements, and for dust evolution in SNRs generally.


title: On Utility Maximisation Under Model Uncertainty in Discrete-Time Markets

abstract: We study the problem of maximising terminal utility for an agent facing model uncertainty, in a frictionless discrete-time market with one safe asset and finitely many risky assets. We show that an optimal investment strategy exists if the utility function, defined either over the positive real line or over the whole real line, is bounded from above. We further find that the boundedness assumption can be dropped provided that we impose suitable integrability conditions, related to some strengthened form of no-arbitrage. These results are obtained in an alternative framework for model uncertainty, where all possible dynamics of the stock prices are represented by a collection of stochastic processes on the same filtered probability space, rather than by a family of probability measures.


title: Geolocating social media posts for emergency mapping

abstract: The demo will illustrate the features of a webGIS interface to support the rapid mapping activities after a natural disaster, with the goal of providing additional information from social media to the mapping operators. This demo shows the first results of the E2mC H2020 European project, where the goal is to extract precisely located information from available social media sources, providing accurate geolocating functionalities and, starting from posts searched in Twitter, extending the social media exploration to Flickr, YouTube, and Instagram.


title: Testing the Number of Regimes in Markov Regime Switching Models

abstract: Markov regime switching models have been used in numerous empirical studies in economics and finance. However, the asymptotic distribution of the likelihood ratio test statistic for testing the number of regimes in Markov regime switching models has been an unresolved problem. This paper derives the asymptotic distribution of the likelihood ratio test statistic for testing the null hypothesis of $M_0$ regimes against the alternative hypothesis of $M_0 + 1$ regimes for any $M_0 \geq 1$ both under the null hypothesis and under local alternatives. We show that the contiguous alternatives converge to the null hypothesis at a rate of $n^{-1/8}$ in regime switching models with normal density. The asymptotic validity of the parametric bootstrap is also established.


title: Unpacking the Social Media Bot: A Typology to Guide Research and Policy

abstract: Amidst widespread reports of digital influence operations during major elections, policymakers, scholars, and journalists have become increasingly interested in the political impact of social media 'bots.' Most recently, platform companies like Facebook and Twitter have been summoned to testify about bots as part of investigations into digitally-enabled foreign manipulation during the 2016 US Presidential election. Facing mounting pressure from both the public and from legislators, these companies have been instructed to crack down on apparently malicious bot accounts. But as this article demonstrates, since the earliest writings on bots in the 1990s, there has been substantial confusion as to exactly what a 'bot' is and what exactly a bot does. We argue that multiple forms of ambiguity are responsible for much of the complexity underlying contemporary bot-related policy, and that before successful policy interventions can be formulated, a more comprehensive understanding of bots --- especially how they are defined and measured --- will be needed. In this article, we provide a history and typology of different types of bots, provide clear guidelines to better categorize political automation and unpack the impact that it can have on contemporary technology policy, and outline the main challenges and ambiguities that will face both researchers and legislators concerned with bots in the future.


title: A Convex Optimization Framework for Constrained Concurrent Motion Control of a Hybrid Redundant Surgical System

abstract: We present a constrained motion control framework for a redundant surgical system designed for minimally invasive treatment of pelvic osteolysis. The framework comprises a kinematics model of a six Degrees-of-Freedom (DoF) robotic arm integrated with a one DoF continuum manipulator as well as a novel convex optimization redundancy resolution controller. To resolve the redundancy resolution problem, formulated as a constrained l2-regularized quadratic minimization, we study and evaluate the potential use of an optimally tuned alternating direction method of multipliers (ADMM) algorithm. To this end, we prove global convergence of the algorithm at linear rate and propose expressions for the involved parameters resulting in a fast convergence. Simulations on the robotic system verified our analytical derivations and showed the capability and robustness of the ADMM algorithm in constrained motion control of our redundant surgical system.


title: Interpolation between H\" older and Lebesgue spaces with applications

abstract: Classical interpolation inequality of the type $\|u\|_{X}\leq C\|u\|_{Y}^{\theta}\|u\|_{Z}^{1-\theta}$ is well known in the case when $X$, $Y$, $Z$ are Lebesgue spaces. In this paper we show that this result may be extended by replacing norms $\|\cdot\|_{Y}$ or $\|\cdot\|_{X}$ by suitable H\" older semi-norm. We shall even prove sharper version involving weak Lorentz norm. We apply this result to prove the Gagliardo--Nirenberg inequality for a wider scale of parameters.


title: Sector-Based Radio Resource Allocation (SBRRA) Algorithm for Better Quality of Service and Experience in Device-to-Device (D2D) Communication

abstract: The mounting content sharing among users has resulted in a considerable rise in wireless data traffic, pressurizing the cellular networks to undergo a suitable upheaval. A competent technology of the fifth-generation networks (5G) for efficiently supporting proximity-based applications is Device-to-Device (D2D) communication, underlaying cellular networks. Significant advances have been made till date, for allocating resources to D2D users in cellular networks, such that sharing of spectral resources between cellular and D2D users is carried out in a coordinated manner. In this paper, a sector-based radio resource allocation (SBRRA) algorithm for resource block allocation to D2D pairs has been proposed, where the number of resource blocks (RBs) is allocated to each D2D pair in an adaptive manner, based on the demanded application by each pair. Different applications demand a varying number of RBs, in accordance with their priority. This algorithm focusses on the use of sectored antennas at the base station, for a better performance and low complexity. Extensive simulations are carried out, considering real-time scenario, for ensuring satisfactory Quality of Service (QoS) and Quality of Experience (QoE) by the users. The efficiency of the proposed scheme is proved by comparing it with the RB allocation using Hidden Markov Model (HMM).


title: Scene recognition with CNNs: objects, scales and dataset bias

abstract: Since scenes are composed in part of objects, accurate recognition of scenes requires knowledge about both scenes and objects. In this paper we address two related problems: 1) scale induced dataset bias in multi-scale convolutional neural network (CNN) architectures, and 2) how to combine effectively scene-centric and object-centric knowledge (i.e. Places and ImageNet) in CNNs. An earlier attempt, Hybrid-CNN, showed that incorporating ImageNet did not help much. Here we propose an alternative method taking the scale into account, resulting in significant recognition gains. By analyzing the response of ImageNet-CNNs and Places-CNNs at different scales we find that both operate in different scale ranges, so using the same network for all the scales induces dataset bias resulting in limited performance. Thus, adapting the feature extractor to each particular scale (i.e. scale-specific CNNs) is crucial to improve recognition, since the objects in the scenes have their specific range of scales. Experimental results show that the recognition accuracy highly depends on the scale, and that simple yet carefully chosen multi-scale combinations of ImageNet-CNNs and Places-CNNs, can push the state-of-the-art recognition accuracy in SUN397 up to 66.26% (and even 70.17% with deeper architectures, comparable to human performance).


title: Light-controlled Assembly of Active Colloidal Molecules

abstract: Thanks to a constant energy input, active matter can self-assemble into phases with complex architectures and functionalities such as living clusters that dynamically form, reshape and break-up, which are forbidden in equilibrium materials by the entropy maximization (or free energy minimization) principle. The challenge to control this active self-assembly has evoked widespread efforts typically hinging on engineering of the properties of individual motile constituents. Here, we provide a different route, where activity occurs as an emergent phenomenon only when individual building blocks bind together in a way that we control by laser light. Using experiments and simulations of two species of immotile microspheres, we exemplify this route by creating active molecules featuring a complex array of behaviors, becoming migrators, spinners and rotators. The possibility to control the dynamics of active self-assembly via light-controllable nonreciprocal interactions will inspire new approaches to understand living matter and to design active materials.


title: Counter-propagating waves in a system of transport-reaction equations

abstract: Hyperbolic transport-reaction equations are abundant in the description of movement of motile organisms. Here, we focus on system of four coupled transport-reaction equations that arises from an age-structuring of a species of turning individuals. The highlight consists of the explicit construction and characterization of counter-propagating traveling waves, patterns which have been observed in bacterial colonies. Stability analysis reveals conditions for the wave formation as well as pulsating-in-time spatially constant solutions.


title: Raman Detection Threshold Measurements for Acetic Acid in Martian Regolith Simulant JSC-1 in the Presence of Hydrated Metallic Sulfates

abstract: Several measurements, including the data from the laser spectrometer in the Sample Analysis at Mars instrument suite on the Curiosity Rover in Gale crater, have measured seasonal variations in atmospheric methane at the 0.7 ppbv (parts per billion volume) level. As a result, models have been proposed to understand the methane production including novel chemical, geological, meteorological, and biological mechanisms. Biological models often rely on methanogenic extremophile archaea which might be similar to some Earth based organisms which can be studied exhaustively in a laboratory setting. Such organisms might thrive in a subsurface ecosystem involving water and methane and as such could leave a unique biosignature in the Martian regolith that could be preserved over extended periods of time. The resulting mixture of carboxyl, acetyl and hydroxyl groups blended with the metallic sulfates and regolith constituents could produce a Raman signature detectable on the next rover mission if the concentration is high enough. Here we measure the 532 nm and the 780 nm Raman spectra of a variety of molar concentrations of several hydrated metallic sulfates mixed with JSC-1 Martian regolith simulant to identify the signature concentration thresholds for the isolated C=O Raman peak from acetic acid. We find a Raman peak removed from most fragmented molecular Raman bands near 1608 cm-1 using the 532 nm laser for the Fe and Mg sulfates to yield a threshold in the 120 - 160 ppmv range for a cutoff S/N of 2.


title: On large $q$ expansion in the Sachdev-Ye-Kitaev model

abstract: We consider the Sachdev-Ye-Kitaev (SYK) model where interaction involves $q$ fermions at a time. We find the next order correction to the thermal two-point function in the large $q$ expansion. Using this result we find the next order correction to the SYK free energy.


title: Existence of ground state solution and concentration of maxima for a class of indefinite variational problems

abstract: In this paper we study the existence of ground state solution and concentration of maxima for a class of strongly indefinite problem like $$ \left\{\begin{array}{l} -\Delta u+V(x)u=A(\epsilon x)f(u) \quad \mbox{in} \quad \R^{N}, \\ u\in H^{1}(\R^{N}), \end{array}\right. \eqno{(P)_{\epsilon}} $$ where $N \geq 1$, $\epsilon$ is a positive parameter, $f: \mathbb{R} \to \mathbb{R}$ is a continuous function with subcritical growth and $V,A: \mathbb{R}^{N} \to \mathbb{R}$ are continuous functions verifying some technical conditions. Here $V$ is a $\mathbb{Z}^N$-periodic function, $0 \not\in \sigma(-\Delta + V)$, the spectrum of $-\Delta +V$, and $$ 0 < \inf_{x \in \R^{N}}A(x)\leq \displaystyle\lim_{|x|\rightarrow+\infty}A(x)<\sup_{x \in \R^{N}}A(x). $$


title: Partially Coherent Sources with Inhomogeneous Azimuth

abstract: A new model of physically realizable electromagnetic source is proposed. The source is partially coherent and non-uniformly totally polarized. The coherence and polarization characteristics of this new source are analyzed. The spatial coherence area of the source can be easily modified at will. The state of polarization is linear across the transverse plane of the source with an azimuth that varies from point to point in a different way depending on the selected values of the parameters that define the source.


title: Optimal Co-Optimizing of Distributed Generation and ESS Facilities Considering Voltage Stability

abstract: The paper deals with the optimal sizing and allocation of dispersed generation, and distributed storage systems. The optimization aims at minimizing the sum of the costs sustained by the distributor for the power losses, for network upgrading, for carrying out the reactive power service and the costs of storage and capacitor installation, over a planning period of several years. Also, analysis of the effect of distributed generation capacity and location on voltage stability enhancement of distribution networks are presented in this paper. The analysis is performed using a steady state voltage stability index which can be evaluated at each node of the distribution system. Different optimal capacities and locations are used to check this effect. The location of distributed generation has the main effect on the system voltage stability. Voltage stability should be taken into account as one of the objectives of distributed generation optimal allocation techniques. The analysis is conducted on a 69-node distribution system in order to show the feasibility of the proposed procedure.


title: The jet-ISM interactions in IC 5063

abstract: The interstellar medium of the radio galaxy IC 5063 is highly perturbed by an AGN jet expanding in the gaseous disc of the galaxy. We model this interaction with relativistic hydrodynamic simulations and multiphase initial conditions for the interstellar medium and compare the results with recent observations. As the jets flood through the inter-cloud channels of the disc, they ablate, accelerate, and disperse clouds to velocities exceeding $400 \mbox{km s}^{-1}$. Clouds are also destroyed or displaced in bulk from the central regions of the galaxy. Our models with jet powers of $10^{44} \mbox{erg s}^{-1}$ and $10^{45} \mbox{erg s}^{-1}$ are capable of reproducing many of the features seen in the position-velocity diagram published in Morganti et al. (2015) and confirm the notion that the jet is responsible for the strongly perturbed gas dynamics seen in the ionized, neutral, and molecular gas phases. In our simulations, we also see strong venting of the jet plasma perpendicular to the disc, which entrains clumps and diffuse filaments into the halo of the galaxy. Our simulations are the first 3D hydrodynamic simulations of the jet and ISM of IC 5063.


title: Static and Dynamic Verification of Relational Properties on Self-Composed C Code

abstract: Function contracts are a well-established way of formally specifying the intended behavior of a function. However, they usually only describe what should happen during a single call. Relational properties, on the other hand, link several function calls. They include such properties as non-interference, continuity and monotonicity. Other examples relate sequences of function calls, for instance, to show that decrypting an encrypted message with the appropriate key gives back the original message. Such properties cannot be expressed directly in the traditional setting of modular deductive verification, but are amenable to verification through self-composition. This paper presents a verification technique dedicated to relational properties in C programs and its implementation in the form of a FRAMA-C plugin called RPP and based on self-composition. It supports functions with side effects and recursive functions. The proposed approach makes it possible to prove a relational property, to check it at runtime, to generate a counterexample using testing and to use it as a hypothesis in the subsequent verification. Our initial experiments on existing benchmarks confirm that the proposed technique is helpful for static and dynamic analysis of relational properties.


title: Limiting Distributions of Spectral Radii for Product of Matrices from the Spherical Ensemble

abstract: Consider the product of $m$ independent $n\times n$ random matrices from the spherical ensemble for $m\ge 1$. The spectral radius is defined as the maximum absolute value of the $n$ eigenvalues of the product matrix. When $m=1$, the limiting distribution for the spectral radii has been obtained by Jiang and Qi (2017). In this paper, we investigate the limiting distributions for the spectral radii in general. When $m$ is a fixed integer, we show that the spectral radii converge weakly to distributions of functions of independent Gamma random variables. When $m=m_n$ tends to infinity as $n$ goes to infinity, we show that the logarithmic spectral radii have a normal limit.


title: Density profiles of a self-gravitating lattice gas in one, two, and three dimensions

abstract: We consider a lattice gas in spaces of dimensionality $\mathcal{D}=1,2,3$. The particles are subject to a hardcore exclusion interaction and an attractive pair interaction that satisfies Gauss' law as do Newtonian gravity in $\mathcal{D}=3$, a logarithmic potential in $\mathcal{D}=2$, and a distance-independent force in $\mathcal{D}=1$. Under mild additional assumptions regarding symmetry and fluctuations we investigate equilibrium states of self-gravitating material clusters, in particular radial density profiles for closed and open systems. We present exact analytic results in several instances and high-precision numerical data in others. The density profile of a cluster with finite mass is found to exhibit exponential decay in $\mathcal{D}=1$ and power-law decay in $\mathcal{D}=2$ with temperature-dependent exponents in both cases. In $\mathcal{D}=2$ the gas evaporates in a continuous transition at a nonzero critical temperature. We describe clusters of infinite mass in $\mathcal{D}=3$ with a density profile consisting of three layers (core, shell, halo) and an algebraic large-distance asymptotic decay. In $\mathcal{D}=3$ a cluster of finite mass can be stabilized at $T>0$ via confinement to a sphere of finite radius. In some parameter regime, the gas thus enclosed undergoes a discontinuous transition between distinct density profiles. For the free energy needed to identify the equilibrium state we introduce a construction of gravitational self-energy that works in all $\mathcal{D}$ for the lattice gas. The decay rate of the density profile of an open cluster is shown to transform via a stretched exponential for $1<\mathcal{D}<2$ whereas it crosses over from one power-law at intermediate distances to a different power-law at larger distances for $2<\mathcal{D}<3$.


title: Bayesian Deep Convolutional Encoder-Decoder Networks for Surrogate Modeling and Uncertainty Quantification

abstract: We are interested in the development of surrogate models for uncertainty quantification and propagation in problems governed by stochastic PDEs using a deep convolutional encoder-decoder network in a similar fashion to approaches considered in deep learning for image-to-image regression tasks. Since normal neural networks are data intensive and cannot provide predictive uncertainty, we propose a Bayesian approach to convolutional neural nets. A recently introduced variational gradient descent algorithm based on Stein's method is scaled to deep convolutional networks to perform approximate Bayesian inference on millions of uncertain network parameters. This approach achieves state of the art performance in terms of predictive accuracy and uncertainty quantification in comparison to other approaches in Bayesian neural networks as well as techniques that include Gaussian processes and ensemble methods even when the training data size is relatively small. To evaluate the performance of this approach, we consider standard uncertainty quantification benchmark problems including flow in heterogeneous media defined in terms of limited data-driven permeability realizations. The performance of the surrogate model developed is very good even though there is no underlying structure shared between the input (permeability) and output (flow/pressure) fields as is often the case in the image-to-image regression models used in computer vision problems. Studies are performed with an underlying stochastic input dimensionality up to $4,225$ where most other uncertainty quantification methods fail. Uncertainty propagation tasks are considered and the predictive output Bayesian statistics are compared to those obtained with Monte Carlo estimates.