title: Black is to Criminal as Caucasian is to Police: Detecting and Removing Multiclass Bias in Word Embeddings

abstract: Online texts -- across genres, registers, domains, and styles -- are riddled with human stereotypes, expressed in overt or subtle ways. Word embeddings, trained on these texts, perpetuate and amplify these stereotypes, and propagate biases to machine learning models that use word embeddings as features. In this work, we propose a method to debias word embeddings in multiclass settings such as race and religion, extending the work of (Bolukbasi et al., 2016) from the binary setting, such as binary gender. Next, we propose a novel methodology for the evaluation of multiclass debiasing. We demonstrate that our multiclass debiasing is robust and maintains the efficacy in standard NLP tasks.


title: Explicit Stencil Computation Schemes Generated by Poisson's Formula for the 2D Wave Equation

abstract: A new approach to building explicit time-marching stencil computation schemes for the transient 2D acoustic wave equation is implemented. It is based on using Poisson's formula and its three time level modification combined with polynomial stencil interpolation of the solution at each time-step and exact integration. The time-stepping algorithm consists of two explicit stencil computation procedures: a first time-step procedure incorporating the initial conditions and a two-step scheme for the second and next time-steps. Three particular explicit stencil schemes (with five, nine and 13 space points) are constructed using this approach. Their stability regions are presented. Accuracy advantages of the new schemes in comparison with conventional finite-difference schemes are demonstrated by simulation using an exact benchmark solution.


title: Simple Question Answering with Subgraph Ranking and Joint-Scoring

abstract: Knowledge graph based simple question answering (KBSQA) is a major area of research within question answering. Although only dealing with simple questions, i.e., questions that can be answered through a single knowledge base (KB) fact, this task is neither simple nor close to being solved. Targeting on the two main steps, subgraph selection and fact selection, the research community has developed sophisticated approaches. However, the importance of subgraph ranking and leveraging the subject--relation dependency of a KB fact have not been sufficiently explored. Motivated by this, we present a unified framework to describe and analyze existing approaches. Using this framework as a starting point, we focus on two aspects: improving subgraph selection through a novel ranking method and leveraging the subject--relation dependency by proposing a joint scoring CNN model with a novel loss function that enforces the well-order of scores. Our methods achieve a new state of the art (85.44% in accuracy) on the SimpleQuestions dataset.


title: Inclusive scattering matrix and scattering of quasiparticles

abstract: The quantum theory can be formulated in the language of positive functionals on Weyl or Clifford algebra ($L$-functionals). It is shown that this language gives simple understanding of diagrams of Keldysh formalism (that coincide in our case with the diagrams of thermo-field dynamics). The matrix elements of the scattering matrix in the formalism of $L$-functionals are related to inclusive cross-sections, therefore we suggest the name "inclusive scattering matrix" for this notion. The inclusive scattering matrix can be expressed in terms of on-shell values of generalized Green functions. This notion is necessary if we want to analyze collisions of quasiparticles.

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


title: Algebraic Quantum Field Theory -- an introduction

abstract: We give a pedagogical introduction to algebraic quantum field theory (AQFT), with the aim of explaining its key structures and features. Topics covered include: algebraic formulations of quantum theory and the GNS representation theorem, the appearance of unitarily inequivalent representations in QFT (exemplified by the van Hove model), the main assumptions of AQFT and simple models thereof, the spectrum condition, Reeh--Schlieder theorem, split property, the universal type of local algebras, and the theory of superselection sectors. The abstract discussion is illustrated by concrete examples. One of our concerns is to explain various ways in which quantum field theory differs from quantum mechanics, not just in terms of technical detail, but in terms of physical content. The text is supplemented by exercises and appendices that enlarge on some of the relevant mathematical background. These notes are based on lectures given by CJF for the International Max Planck Research School at the Albert Einstein Institute, Golm (October, 2018) and by KR at the Raman Research Institute, Bangalore (January, 2019).


title: Separating effect from significance in Markov chain tests

abstract: We give qualitative and quantitative improvements to theorems which enable significance testing in Markov Chains, with a particular eye toward the goal of enabling strong, interpretable, and statistically rigorous claims of political gerrymandering. Our results can be used to demonstrate at a desired significance level that a given Markov Chain state (e.g., a districting) is extremely unusual (rather than just atypical) with respect to the fragility of its characteristics in the chain. We also provide theorems specialized to leverage quantitative improvements when there is a product structure in the underlying probability space, as can occur due to geographical constraints on districtings.


title: Noncommutative scalar field in the non-extremal Reissner-Nordstr\"om background: QNM spectrum

abstract: In our previous work [18] we constructed a model of a noncommutative, charged and massive scalar field based on the angular twist. Then we used this model to analyze the motion of the scalar field in the Reissner-Nordstr\"om black hole background. In particular, we determined the QNM spectrum analytically in the near-extremal limit. To broaden our analysis, in this paper we apply a well defined numerical method, the continued fraction method and calculate the QNM spectrum for a non-extremal Reissner-Nordstr\"om black hole. To check the validity of our analytic calculations, we compare results of the continued fraction method in the near extremal limit with the analytic results obtained in the previous paper. We find that the results are in good agreement. For completeness, we also study the QNM spectrum in the WKB approximation.


title: M/M/$c$ Queues and the Poisson Clumping Heuristic

abstract: In continuous time, customers arrive at random. Each waits until one of $c$ servers is available; each thereafter departs at random. The distribution of maximum line length of idle customers was studied over 25 years ago. We revisit two good approximations of this, employing a discrete Gumbel formulation and detailed graphics to describe simulation outcomes.


title: Evaluating KGR10 Polish word embeddings in the recognition of temporal expressions using BiLSTM-CRF

abstract: The article introduces a new set of Polish word embeddings, built using KGR10 corpus, which contains more than 4 billion words. These embeddings are evaluated in the problem of recognition of temporal expressions (timexes) for the Polish language. We described the process of KGR10 corpus creation and a new approach to the recognition problem using Bidirectional Long-Short Term Memory (BiLSTM) network with additional CRF layer, where specific embeddings are essential. We presented experiments and conclusions drawn from them.


title: Wide-range Prandtl/Schmidt number power spectrum of optical turbulence and its application to oceanic light propagation

abstract: Light influenced by the turbulent ocean can be fully characterized with the help of the power spectrum of the water's refractive index fluctuations, resulting from the combined effect of two scalars, temperature and salinity concentration advected by the velocity field. The Nikishovs' model [ Fluid Mech. Res. 27, 8298 (2000)] frequently used in the analysis of light evolution through the turbulent ocean channels is the linear combination of the temperature spectrum, the salinity spectrum and their co-spectrum, each being described by an approximate expression developed by Hill [ J. Fluid Mech. 88, 541562 (1978)] in the first of his four suggested models. The fourth of the Hill's models provides much more precise power spectrum than the first one expressed via a non-linear differential equation that does not have a closed-form solution. We develop an accurate analytic approximation to the fourth Hill's model valid for Prandtl/Schmidt numbers in the interval [3, 3000] and use it for the development of a more precise oceanic power spectrum. To illustrate the advantage of our model, we include numerical examples relating to the spherical wave scintillation index evolving in the underwater turbulent channels with different average temperatures, and, hence, different Prandtl numbers for temperature and different Schmidt numbers for salinity. Since our model is valid for a large range of Prandtl number (or/and Schmidt number), it can be readily adjusted to oceanic waters with seasonal or extreme average temperature and/or salinity or any other turbulent fluid with one or several advected quantities.


title: Task Oriented Channel State Information Quantization

abstract: In this paper, we propose a new perspective for quantizing a signal and more specifically the channel state information (CSI). The proposed point of view is fully relevant for a receiver which has to send a quantized version of the channel state to the transmitter. Roughly, the key idea is that the receiver sends the right amount of information to the transmitter so that the latter be able to take its (resource allocation) decision. More formally, the decision task of the transmitter is to maximize an utility function u(x;g) with respect to x (e.g., a power allocation vector) given the knowledge of a quantized version of the function parameters g. We exhibit a special case of an energy-efficient power control (PC) problem for which the optimal task oriented CSI quantizer (TOCQ) can be found analytically. For more general utility functions, we propose to use neural networks (NN) based learning. Simulations show that the compression rate obtained by adapting the feedback information rate to the function to be optimized may be significantly increased.


title: A comparative study of physics-informed neural network models for learning unknown dynamics and constitutive relations

abstract: We investigate the use of discrete and continuous versions of physics-informed neural network methods for learning unknown dynamics or constitutive relations of a dynamical system. For the case of unknown dynamics, we represent all the dynamics with a deep neural network (DNN). When the dynamics of the system are known up to the specification of constitutive relations (that can depend on the state of the system), we represent these constitutive relations with a DNN. The discrete versions combine classical multistep discretization methods for dynamical systems with neural network based machine learning methods. On the other hand, the continuous versions utilize deep neural networks to minimize the residual function for the continuous governing equations. We use the case of a fedbatch bioreactor system to study the effectiveness of these approaches and discuss conditions for their applicability. Our results indicate that the accuracy of the trained neural network models is much higher for the cases where we only have to learn a constitutive relation instead of the whole dynamics. This finding corroborates the well-known fact from scientific computing that building as much structural information is available into an algorithm can enhance its efficiency and/or accuracy.

categories: cs.LG cs.NA math.DS math.NA physics.comp-ph


title: Minimal kernels and compact analytic objects in complex surfaces

abstract: In this paper, we want to study the link between the presence of compact objects with some analytic structure and the global geometry of a weakly complete surface. We begin with a brief survey of some now classic results on the local geometry around a (complex) curve, which depends on the sign of its self-intersection and, in the flat case, on some more refined invariants (see the works of Grauert, Suzuki, Ueda). Then, we recall some results about the propagation of compact curves and the existence of holomorphic functions (from the works of Nishino and Ohsawa). With such considerations in mind, we give an overview of the classification results for weakly complete surfaces that we obtained in two joint papers with Slodkowski (see [MST18], [MST17] and we present some new results which stem from this somehow more local (or less global) viewpoint (see Sections 3.2, 3.3 and 4).


title: Difference of source regions between fast and slow coronal mass ejections

abstract: Coronal mass ejections (CMEs) are tightly related to filament eruptions and usually are their continuation in the upper solar corona. It is common practice to divide all observed CMEs into fast and slow ones. Fast CMEs usually follow eruptive events in active regions near big sunspot groups and associated with major solar flares. Slow CMEs are more related to eruptions of quiescent prominences located far from active regions. We analyze ten eruptive events with particular attention to the events on 2013 September 29 and on 2016 January 26, one of which was associated with a fast CME, while another was followed by a slow CME. We estimated the initial store of free magnetic energy in the two regions and show the resemblance of pre-eruptive situations. The difference of late behaviour of the two eruptive prominences is a consequence of the different structure of magnetic field above the filaments. We estimated this structure on the basis of potential magnetic field calculations. Analysis of other eight events confirmed that all fast CMEs originate in regions with rapidly changing with height value and direction of coronal magnetic field.


title: Transferring Knowledge Fragments for Learning Distance Metric from A Heterogeneous Domain

abstract: The goal of transfer learning is to improve the performance of target learning task by leveraging information (or transferring knowledge) from other related tasks. In this paper, we examine the problem of transfer distance metric learning (DML), which usually aims to mitigate the label information deficiency issue in the target DML. Most of the current Transfer DML (TDML) methods are not applicable to the scenario where data are drawn from heterogeneous domains. Some existing heterogeneous transfer learning (HTL) approaches can learn target distance metric by usually transforming the samples of source and target domain into a common subspace. However, these approaches lack flexibility in real-world applications, and the learned transformations are often restricted to be linear. This motivates us to develop a general flexible heterogeneous TDML (HTDML) framework. In particular, any (linear/nonlinear) DML algorithms can be employed to learn the source metric beforehand. Then the pre-learned source metric is represented as a set of knowledge fragments to help target metric learning. We show how generalization error in the target domain could be reduced using the proposed transfer strategy, and develop novel algorithm to learn either linear or nonlinear target metric. Extensive experiments on various applications demonstrate the effectiveness of the proposed method.


title: Quality-Aware Broadcasting Strategies for Position Estimation in VANETs

abstract: The dissemination of vehicle position data all over the network is a fundamental task in Vehicular Ad Hoc Network (VANET) operations, as applications often need to know the position of other vehicles over a large area. In such cases, inter-vehicular communications should be exploited to satisfy application requirements, although congestion control mechanisms are required to minimize the packet collision probability. In this work, we face the issue of achieving accurate vehicle position estimation and prediction in a VANET scenario. State of the art solutions to the problem try to broadcast the positioning information periodically, so that vehicles can ensure that the information their neighbors have about them is never older than the inter-transmission period. However, the rate of decay of the information is not deterministic in complex urban scenarios: the movements and maneuvers of vehicles can often be erratic and unpredictable, making old positioning information inaccurate or downright misleading. To address this problem, we propose to use the Quality of Information (QoI) as the decision factor for broadcasting. We implement a threshold-based strategy to distribute position information whenever the positioning error passes a reference value, thereby shifting the objective of the network to limiting the actual positioning error and guaranteeing quality across the VANET. The threshold-based strategy can reduce the network load by avoiding the transmission of redundant messages, as well as improving the overall positioning accuracy by more than 20% in realistic urban scenarios.


title: Analyzing and Interpreting Neural Networks for NLP: A Report on the First BlackboxNLP Workshop

abstract: The EMNLP 2018 workshop BlackboxNLP was dedicated to resources and techniques specifically developed for analyzing and understanding the inner-workings and representations acquired by neural models of language. Approaches included: systematic manipulation of input to neural networks and investigating the impact on their performance, testing whether interpretable knowledge can be decoded from intermediate representations acquired by neural networks, proposing modifications to neural network architectures to make their knowledge state or generated output more explainable, and examining the performance of networks on simplified or formal languages. Here we review a number of representative studies in each category.


title: New type Pythagorean fuzzy soft set and decision-making application

abstract: We define the Pythagorean fuzzy parameterized soft set and investigate some properties of the new set. Further, we propose to the solution of decision-making application for the Pythagorean fuzzy parameterized soft set and other related concepts.


title: On the Reduciblity of a Certain Type of Rank 3 Uniform Oriented Matroid by a Point

abstract: For a positive integer $n\geq 3$, the sides and diagonals of a convex $n$-gon divide the interior of the convex $n$-gon into finitely (polynomial in $n$) many regions bounded by them. In this article, we associate to every region a unique $n$-cycle in the symmetric group $S_n$ of a certain type (defined as $2$-standard consecutive cycle) by studying point arrangements in the plane. Then we find that there are more (exponential in $n$) number of such cycles leading to the conclusion that not every region labelled by a cycle appears in every convex $n$-gon. In fact most of them do not occur in any given single convex $n$-gon. Later in the main theorem of this article we characterize combinatorially those cycles (defined as definite cycles) whose corresponding regions occur in every convex $n$-gon and those cycles (defined as indefinite cycles) whose corresponding regions do not occur in every convex $n$-gon. As a consequence we characterize those one point extensions of a uniform rank $3$ convex oriented matroid for which the one point extension is reducible by the, one point, when it lies inside the convex hull.


title: Plasma living in a curved surface at some special temperature

abstract: The simplest statistical mechanics model of a Coulomb plasma in two spatial dimensions admits an exact analytic solution at some special temperature in several (curved) surfaces. We present in a unifying perspective these solutions for the (non-quantum) plasma, made of point particles carrying an absolute charge $e$, in thermal equilibrium at a temperature $T=e^2/2k_B$, with $k_B$ Boltzmann's constant, discussing the importance of having an exact solution, the role of the curvature of the surface, and the densities of the plasma.

categories: cond-mat.stat-mech cond-mat.mtrl-sci physics.class-ph physics.plasm-ph


title: Extracting analytical one-loop amplitudes from numerical evaluations

abstract: In this article we present a method to generate analytic expressions for the integral coefficients of loop amplitudes using numerical evaluations only. We use high-precision arithmetic to explore the singularity structure of the coefficients and decompose them into parts of manageable complexity. To illustrate the usability of our method we provide analytical expressions for all helicity configurations of the colour-ordered six-point gluon amplitudes at one loop with a gluon in the loop.


title: Determining the recurrence timescale of long-lasting YSO outbursts

abstract: We have determined the rate of large accretion events in class I and II young stellar objects (YSOs) by comparing the all-sky digitised photographic plate surveys provided by SuperCOSMOS with the latest data release from Gaia (DR2). The long mean baseline of 55 years along with a large sample of class II YSOs ($\simeq$15,000) allows us to study approximately 1 million YSO-years. We find 139 objects with $\Delta R\geq1$~mag, most of which are found at amplitudes between 1 and 3 mag. The majority of YSOs in this group show irregular variability or long-lasting fading events, which is best explained as hot spots due to accretion or by variable extinction. There is a tail of YSOs at $\Delta R\geq3$~mag and they seem to represent a different population. Surprisingly many objects in this group show high-amplitude irregular variability over timescales shorter than 10 years, in contrast with the view that high-amplitude objects always have long outbursts. However, we find 6 objects that are consistent with undergoing large, long lasting accretion events, 3 of them previously unknown. This yields an outburst recurrence timescale of 112 kyr, with a 68\% confidence interval [74 to 180] kyr. This represents the first robust determination of the outburst rate in class II YSOs and shows that YSOs in their planet-forming stage do in fact undergo large accretion events, and with timescales of $\simeq$100,000 years. In addition, we find that outbursts in the class II stage are $\simeq$10 times less frequent than during the class I stage.


title: A geometric approach to K-homology for Lie manifolds

abstract: We show that the computation of the Fredholm index of a fully elliptic pseudodifferential operator on an integrated Lie manifold can be reduced to the computation of the index of a Dirac operator, perturbed by a smoothing operator, canonically associated, via the so-called clutching map. To this end we adapt to our framework ideas coming from Baum-Douglas geometric $K$-homology and in particular we introduce a notion of geometric cycles, that can be categorized as a variant of the famous geometric $K$-homology groups, for the specific situation here. We also define a comparison map between this geometric $K$-homology theory and a relative $K$-theory group, directly associated to a fully elliptic pseudodifferential operator.


title: A genus formula for the positive \'{e}tale wild kernel

abstract: Let $F$ be a number field and let $i\geq 2$ be an integer. In this paper, we study the positive \'{e}tale wild kernel $\mathrm{WK}^{\mbox{\'{e}t},+}_{2i-2}F$, which is the twisted analogue of the $2$-primary part of the narrow class group. If $E/F$ is a Galois extension of number fields with Galois group $G$, we prove a genus formula relating the order of the groups $ (\mathrm{WK}^{\mbox{\'{e}t},+}_{2i-2}E)_{G}$ and $\mathrm{WK}^{\mbox{\'{e}t},+}_{2i-2}F$.


title: Infinite families of 3-designs from APN functions

abstract: Combinatorial $t$-designs have nice applications in coding theory, finite geometries and several engineering areas. The objective of this paper is to study how to obtain $3$-designs with $2$-transitive permutation groups. The incidence structure formed by the orbits of a base block under the action of the general affine groups, which are $2$-transitive, is considered. A characterization of such incidence structure to be a $3$-design is presented, and a sufficient condition for the stabilizer of a base block to be trivial is given. With these general results, infinite families of $3$-designs are constructed by employing APN functions. Some $3$-designs presented in this paper give rise to self-dual binary codes or linear codes with optimal or best parameters known. Several conjectures on $3$-designs and binary codes are also presented.


title: The Complexity of the Ideal Membership Problem for Constrained Problems Over the Boolean Domain

abstract: Given an ideal $I$ and a polynomial $f$ the Ideal Membership Problem is to test if $f\in I$. This problem is a fundamental algorithmic problem with important applications and notoriously intractable. We study the complexity of the Ideal Membership Problem for combinatorial ideals that arise from constrained problems over the Boolean domain. As our main result, we identify the borderline of tractability. By using Gr\"{o}bner bases techniques, we extend Schaefer's dichotomy theorem [STOC, 1978] which classifies all Constraint Satisfaction Problems over the Boolean domain to be either in P or NP-hard. Moreover, our result implies necessary and sufficient conditions for the efficient computation of Theta Body SDP relaxations, identifying therefore the borderline of tractability for constraint language problems. This paper is motivated by the pursuit of understanding the recently raised issue of bit complexity of Sum-of-Squares proofs [O'Donnell, ITCS, 2017]. Raghavendra and Weitz [ICALP, 2017] show how the Ideal Membership Problem tractability for combinatorial ideals implies bounded coefficients in Sum-of-Squares proofs.


title: Abusive Language Detection with Graph Convolutional Networks

abstract: Abuse on the Internet represents a significant societal problem of our time. Previous research on automated abusive language detection in Twitter has shown that community-based profiling of users is a promising technique for this task. However, existing approaches only capture shallow properties of online communities by modeling follower-following relationships. In contrast, working with graph convolutional networks (GCNs), we present the first approach that captures not only the structure of online communities but also the linguistic behavior of the users within them. We show that such a heterogeneous graph-structured modeling of communities significantly advances the current state of the art in abusive language detection.


title: Alternative approach to gravity and MOND

abstract: The classical gravitational two-body problem is generalized in order to be applicable also to weak gravitational fields. The equation of motion holds both for terrestrial and large cosmic scales, the Newtonian gravitational law represents a mathematical limit of the generalized form. Motivation comes from observational results on rotation curves of galaxies. Existence of a dark matter is not assumed. The crucial laws of physics hold and also the potential energy of the system is symmetric with respect to masses of the two bodies. Shortcomings of the results published for decades, including MOND theories and false-yet-familiar approaches, are overcome. The impact on searching for a fundamental physical theory is stressed. Some of the conventional ideas of the past centuries do not hold for the zone of small accelerations, e.g., the principle of least action using the Lagrangian density of potentials and fields does not work. We may look forward to great changes in our understanding of the evolution of the Universe.


title: Bargmann-Wigner Equations, Fermion-Boson Correspondence and Superradiant Problem in Curved Spacetime

abstract: Bargmann-Wigner equations and their solutions are studied in (3+1)-dimensional curved spacetime. Fermion-Boson correspondence for bi-spinor case is studied through the Bargmann-Wigner equations and solutions over curved spacetime. As an application to scattering phenomena of massive Fermions and Bosons on the rotating black holes, the superradiance with negative energy $(\omega<0)$ and positive effective energy in co-rotating coordinate system near horizon $(\omega-m\Omega_{H}>0)$ is possible to occur as stable physical states in Kerr spacetime.


title: Adiabatic limits, Theta functions, and geometric quantization

abstract: Let $\pi\colon (M,\omega)\to B$ be a non-singular Lagrangian torus fibration on a complete base $B$ with prequantum line bundle $(L,\nabla^L)\to (M,\omega)$. Compactness on $M$ is not assumed. For a positive integer $N$ and a compatible almost complex structure $J$ on $(M,\omega)$ invariant along the fiber of $\pi$, let $D$ be the associated Spin${}^c$ Dirac operator with coefficients in $L^{\otimes N}$. First, in the case where $J$ is integrable, under certain technical condition on $J$, we give a complete orthogonal system $\{ \vartheta_b\}_{b\in B_{BS}}$ of the space of holomorphic $L^2$-sections of $L^{\otimes N}$ indexed by the Bohr-Sommerfeld points $B_{BS}$ such that each $\vartheta_b$ converges to a delta-function section supported on the corresponding Bohr-Sommerfeld fiber $\pi^{-1}(b)$ by the adiabatic(-type) limit. We also explain the relation of $\vartheta_b$ with Jacobi's theta functions when $(M,\omega)$ is $T^{2n}$. Next, in the case where $J$ is not integrable, we give an orthogonal family $\{ {\tilde \vartheta}_b\}_{b\in B_{BS}}$ of $L^2$-sections of $L^{\otimes N}$ indexed by $B_{BS}$ which has the same property as above, and show that each $D{\tilde \vartheta}_b$ converges to $0$ by the adiabatic(-type) limit with respect to the $L^2$-norm.


title: On the number of non-G-equivalent minimal abelian codes

abstract: Let $G$ be a finite abelian group. Ferraz, Guerreiro and Polcino Milies prove that the number of $G$-equivalence classes of minimal abelian codes is equal to the number of $G$-isomorphism classes of subgroups for which corresponding quotients are cyclic. In this article, we prove that the notion of $G$-isomorphism is equivalent to the notion of isomorphism on the set of all subgroups $H$ of $G$ with the property that $G/H$ is cyclic. As an application, we calculate the number of non-$G$-equivalent minimal abelian codes for some specific family of abelian groups. We also prove that the number of non-$G$-equivalent minimal abelian codes is equal to number of divisors of the exponent of $G$ if and only if for each prime $p$ dividing the order of $G$, the Sylow $p$-subgroups of $G$ are homocyclic.


title: Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library

abstract: In this paper we employ two implementations of the fictitious domain (FD) method to simulate water-entry and water-exit problems and demonstrate their ability to simulate practical marine engineering problems. In FD methods, the fluid momentum equation is extended within the solid domain using an additional body force that constrains the structure velocity to be that of a rigid body. Using this formulation, a single set of equations is solved over the entire computational domain. The constraint force is calculated in two distinct ways: one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method and another using a fully-Eulerian approach of the Brinkman penalization (BP) method. Both FSI strategies use the same multiphase flow algorithm that solves the discrete incompressible Navier-Stokes system in conservative form. A consistent transport scheme is employed to advect mass and momentum in the domain, which ensures numerical stability of high density ratio multiphase flows involved in practical marine engineering applications. Example cases of a free falling wedge (straight and inclined) and cylinder are simulated, and the numerical results are compared against benchmark cases in literature.


title: Differentiable Sampling with Flexible Reference Word Order for Neural Machine Translation

abstract: Despite some empirical success at correcting exposure bias in machine translation, scheduled sampling algorithms suffer from a major drawback: they incorrectly assume that words in the reference translations and in sampled sequences are aligned at each time step. Our new differentiable sampling algorithm addresses this issue by optimizing the probability that the reference can be aligned with the sampled output, based on a soft alignment predicted by the model itself. As a result, the output distribution at each time step is evaluated with respect to the whole predicted sequence. Experiments on IWSLT translation tasks show that our approach improves BLEU compared to maximum likelihood and scheduled sampling baselines. In addition, our approach is simpler to train with no need for sampling schedule and yields models that achieve larger improvements with smaller beam sizes.


title: Multicolour chain avoidance in the boolean lattice

abstract: Given a collection of colored chain posets, we estimate the number of colored subsets of the boolean lattice which avoid all chains in the collection.


title: Heterogeneous Multi-task Metric Learning across Multiple Domains

abstract: Distance metric learning (DML) plays a crucial role in diverse machine learning algorithms and applications. When the labeled information in target domain is limited, transfer metric learning (TML) helps to learn the metric by leveraging the sufficient information from other related domains. Multi-task metric learning (MTML), which can be regarded as a special case of TML, performs transfer across all related domains. Current TML tools usually assume that the same feature representation is exploited for different domains. However, in real-world applications, data may be drawn from heterogeneous domains. Heterogeneous transfer learning approaches can be adopted to remedy this drawback by deriving a metric from the learned transformation across different domains. But they are often limited in that only two domains can be handled. To appropriately handle multiple domains, we develop a novel heterogeneous multi-task metric learning (HMTML) framework. In HMTML, the metrics of all different domains are learned together. The transformations derived from the metrics are utilized to induce a common subspace, and the high-order covariance among the predictive structures of these domains is maximized in this subspace. There do exist a few heterogeneous transfer learning approaches that deal with multiple domains, but the high-order statistics (correlation information), which can only be exploited by simultaneously examining all domains, is ignored in these approaches. Compared with them, the proposed HMTML can effectively explore such high-order information, thus obtaining more reliable feature transformations and metrics. Effectiveness of our method is validated by the extensive and intensive experiments on text categorization, scene classification, and social image annotation.


title: Symmetry and monotonicity properties of singular solutions to some cooperative semilinear elliptic systems involving critical nonlinearities

abstract: We investigate qualitative properties of positive singular solutions of some elliptic systems in bounded and unbounded domains. We deduce symmetry and monotonicity properties via the moving plane procedure. Moreover, in the unbounded case, we study some cooperative elliptic systems involving critical nonlinearities in $\mathbb{R}^n$.


title: Convolutive Blind Source Separation on Surface EMG Signals for Respiratory Diagnostics and Medical Ventilation Control

abstract: The electromyogram (EMG) is an important tool for assessing the activity of a muscle and thus also a valuable measure for the diagnosis and control of respiratory support. In this article we propose convolutive blind source separation (BSS) as an effective tool to pre-process surface electromyogram (sEMG) data of the human respiratory muscles. Specifically, the problem of discriminating between inspiratory, expiratory and cardiac muscle activity is addressed, which currently poses a major obstacle for the clinical use of sEMG for adaptive ventilation control. It is shown that using the investigated broadband algorithm, a clear separation of these components can be achieved. The algorithm is based on a generic framework for BSS that utilizes multiple statistical signal characteristics. Apart from a four-channel FIR structure, there are no further restrictive assumptions on the demixing system.


title: ContextDesc: Local Descriptor Augmentation with Cross-Modality Context

abstract: Most existing studies on learning local features focus on the patch-based descriptions of individual keypoints, whereas neglecting the spatial relations established from their keypoint locations. In this paper, we go beyond the local detail representation by introducing context awareness to augment off-the-shelf local feature descriptors. Specifically, we propose a unified learning framework that leverages and aggregates the cross-modality contextual information, including (i) visual context from high-level image representation, and (ii) geometric context from 2D keypoint distribution. Moreover, we propose an effective N-pair loss that eschews the empirical hyper-parameter search and improves the convergence. The proposed augmentation scheme is lightweight compared with the raw local feature description, meanwhile improves remarkably on several large-scale benchmarks with diversified scenes, which demonstrates both strong practicality and generalization ability in geometric matching applications.


title: An initial overview of the extent and structure of recent star formation within the Serpens Molecular Cloud using Gaia Data Release 2

abstract: The dense clusters within the Serpens Molecular Cloud are among the most active regions of nearby star formation. In this paper, we use Gaia DR2 parallaxes and proper motions to statistically measure 1167 kinematic members of Serpens, few of which were previously identified, to evaluate the star formation history of the complex. The optical members of Serpens are concentrated in three distinct groups located at 380 to 480 pc; the densest clusters are still highly obscured by optically-thick dust and have few optical members. The total population of young stars and protostars in Serpens is at least 2000 stars, including past surveys that were most sensitive to protostars and disks, and may be far higher. Distances to dark clouds measured from deficits in star counts are consistent with the distances to the optical star clusters. The Serpens Molecular Cloud is seen in the foreground of the Aquila Rift, dark clouds located at 600 to 700 pc, and behind patchy extinction, here called the Serpens Cirrus, located at ~250 pc. Based on the lack of a distributed population of older stars, the star formation rate throughout the Serpens Molecular Cloud increased by at least a factor of 20 within the past ~5 Myr. The optically bright stars in Serpens Northeast are visible because their natal molecular cloud has been eroded and not because they were flung outwards from a central factory of star formation. The separation between subclusters of 20 to 100 pc and the absence of an older population leads to speculation that an external forcing was needed to trigger the active star formation.


title: Reexamination of the warm inflation curvature perturbations spectrum

abstract: The two approaches to compute perturbations in warm inflation are examined. It is shown that both approaches lead to different expressions for the amplitude of the primordial spectrum, with a difference between them of $\Upsilon/(4H)$ at leading order, where $\Upsilon$ is the dissipation coefficient. In terms of observables, this discrepancy can lead to the spectral index differing by up to order $\mathcal{O}(10^{-3})$, which is within precision demands for current CMB data. Thus, it is important to resolve this ambiguity to have reliable predictions from warm inflation. For this we prove the extent of this discrepancies by deriving a formula for the spectral index and the tensor-to-scalar ratio in each approach. In doing so, we find disparities to be more noticeable in a regime where dissipation is comparable with the expansion rate, which is a very important regime from a phenomenological point of view. To determine the extent of the discrepancy, several cases are examined, including quadratic, quartic and hybrid potentials with quadratic and $T-$dependent dissipative coefficients. The origin of the discrepancy is found to be due to the approximation performed in one of the methods, which underestimates the variation of the momentum perturbation with expansion. Once this is corrected, both approaches are then in agreement.


title: Describing phase transitions in field theory by self-similar approximants

abstract: Self-similar approximation theory is shown to be a powerful tool for describing phase transitions in quantum field theory. Self-similar approximants present the extrapolation of asymptotic series in powers of small variables to the arbitrary values of the latter, including the variables tending to infinity. The approach is illustrated by considering three problems: (i) The influence of the coupling parameter strength on the critical temperature of the O(N)-symmetric multicomponent field theory. (ii) The calculation of critical exponents for the phase transition in the O(N)-symmetric field theory. (iii) The evaluation of deconfinement temperature in quantum chromodynamics. The results are in good agreement with the available numerical calculations, such as Monte Carlo simulations, Pade-Borel summation, and lattice data.


title: Large Margin Multi-modal Multi-task Feature Extraction for Image Classification

abstract: The features used in many image analysis-based applications are frequently of very high dimension. Feature extraction offers several advantages in high-dimensional cases, and many recent studies have used multi-task feature extraction approaches, which often outperform single-task feature extraction approaches. However, most of these methods are limited in that they only consider data represented by a single type of feature, even though features usually represent images from multiple modalities. We therefore propose a novel large margin multi-modal multi-task feature extraction (LM3FE) framework for handling multi-modal features for image classification. In particular, LM3FE simultaneously learns the feature extraction matrix for each modality and the modality combination coefficients. In this way, LM3FE not only handles correlated and noisy features, but also utilizes the complementarity of different modalities to further help reduce feature redundancy in each modality. The large margin principle employed also helps to extract strongly predictive features so that they are more suitable for prediction (e.g., classification). An alternating algorithm is developed for problem optimization and each sub-problem can be efficiently solved. Experiments on two challenging real-world image datasets demonstrate the effectiveness and superiority of the proposed method.


title: Atomistic-Object Kinetic Monte Carlo simulations of irradiation damage in tungsten

abstract: We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-the-fly. Elastic interactions are handled within objects with empirical potentials at short distances, and between spatially distinct objects using the dipole tensor formalism. The model is shown to evolve mobile interstitial clusters in tungsten faster than an equivalent molecular dynamics simulation, even at elevated temperatures. We apply the model to the evolution of complex defects generated using molecular dynamics simulations of primary radiation damage in tungsten. We show that we can evolve defect structures formed in cascade simulations to experimentally observable timescales of seconds while retaining atomistic detail. We conclude that the first few nanoseconds of simulation following cascade initiation would be better performed using molecular dynamics, as this will capture some of the near-temperature-independent evolution of small highly-mobile interstitial clusters. We also conclude that, for the 20keV PKA cascades annealing simulations considered here, internal relaxations of sessile objects difficult to capture using conventional object KMC with idealised object geometries establish the conditions for long timescale evolution.


title: On Functions Weakly Computable by Pushdown Petri Nets and Related Systems

abstract: We consider numerical functions weakly computable by grammar-controlled vector addition systems (GVASes, a variant of pushdown Petri nets). GVASes can weakly compute all fast growing functions $F_\alpha$ for $\alpha<\omega^\omega$, hence they are computationally more powerful than standard vector addition systems. On the other hand they cannot weakly compute the inverses $F_\alpha^{-1}$ or indeed any sublinear function. The proof relies on a pumping lemma for runs of GVASes that is of independent interest.


title: Geostatistical Modeling of Positive Definite Matrices: An Application to Diffusion Tensor Imaging

abstract: Geostatistical modeling for continuous point-referenced data has been extensively applied to neuroimaging because it produces efficient and valid statistical inference. However, diffusion tensor imaging (DTI), a neuroimaging characterizing the brain structure produces a positive definite (p.d.) matrix for each voxel. Current geostatistical modeling has not been extended to p.d. matrices because introducing spatial dependence among positive definite matrices properly is challenging. In this paper, we use the spatial Wishart process, a spatial stochastic process (random field) where each p.d. matrix-variate marginally follows a Wishart distribution, and spatial dependence between random matrices is induced by latent Gaussian processes. This process is valid on an uncountable collection of spatial locations and is almost surely continuous, leading to a reasonable means of modeling spatial dependence. Motivated by a DTI dataset of cocaine users, we propose a spatial matrix-variate regression model based on the spatial Wishart process. A problematic issue is that the spatial Wishart process has no closed-form density function. Hence, we propose approximation methods to obtain a feasible working model. A local likelihood approximation method is also applied to achieve fast computation. The simulation studies and real data analysis demonstrate that the working model produces reliable inference and improved performance compared to other methods.


title: Leveraging the Invariant Side of Generative Zero-Shot Learning

abstract: Conventional zero-shot learning (ZSL) methods generally learn an embedding, e.g., visual-semantic mapping, to handle the unseen visual samples via an indirect manner. In this paper, we take the advantage of generative adversarial networks (GANs) and propose a novel method, named leveraging invariant side GAN (LisGAN), which can directly generate the unseen features from random noises which are conditioned by the semantic descriptions. Specifically, we train a conditional Wasserstein GANs in which the generator synthesizes fake unseen features from noises and the discriminator distinguishes the fake from real via a minimax game. Considering that one semantic description can correspond to various synthesized visual samples, and the semantic description, figuratively, is the soul of the generated features, we introduce soul samples as the invariant side of generative zero-shot learning in this paper. A soul sample is the meta-representation of one class. It visualizes the most semantically-meaningful aspects of each sample in the same category. We regularize that each generated sample (the varying side of generative ZSL) should be close to at least one soul sample (the invariant side) which has the same class label with it. At the zero-shot recognition stage, we propose to use two classifiers, which are deployed in a cascade way, to achieve a coarse-to-fine result. Experiments on five popular benchmarks verify that our proposed approach can outperform state-of-the-art methods with significant improvements.


title: Gaussian belief propagation solvers for nonsymmetric systems of linear equations

abstract: In this paper, we argue for the utility of deterministic inference in the classical problem of numerical linear algebra, that of solving a linear system. We show how the Gaussian belief propagation solver, known to work for symmetric matrices can be modified to handle nonsymmetric matrices. Furthermore, we introduce a new algorithm for matrix inversion that corresponds to the generalized belief propagation derived from the cluster variation method (or Kikuchi approximation). We relate these algorithms to LU and block LU decompositions and provide certain guarantees based on theorems from the theory of belief propagation. All proposed algorithms are compared with classical solvers (e.g., Gauss-Seidel, BiCGSTAB) with application to linear elliptic equations. We also show how the Gaussian belief propagation can be used as multigrid smoother, resulting in a substantially more robust solver than the one based on the Gauss-Seidel iterative method.


title: Weighted Point Cloud Augmentation for Neural Network Training Data Class-Imbalance

abstract: Recent developments in the field of deep learning for 3D data have demonstrated promising potential for end-to-end learning directly from point clouds. However, many real-world point clouds contain a large class im-balance due to the natural class im-balance observed in nature. For example, a 3D scan of an urban environment will consist mostly of road and facade, whereas other objects such as poles will be under-represented. In this paper we address this issue by employing a weighted augmentation to increase classes that contain fewer points. By mitigating the class im-balance present in the data we demonstrate that a standard PointNet++ deep neural network can achieve higher performance at inference on validation data. This was observed as an increase of F1 score of 19% and 25% on two test benchmark datasets; ScanNet and Semantic3D respectively where no class im-balance pre-processing had been performed. Our networks performed better on both highly-represented and under-represented classes, which indicates that the network is learning more robust and meaningful features when the loss function is not overly exposed to only a few classes.


title: Barium & related stars and their white-dwarf companions II. Main-sequence and subgiant stars

abstract: Barium (Ba) dwarfs and CH subgiants are the less-evolved analogues of Ba and CH giants. They are F- to G-type main-sequence stars polluted with heavy elements by a binary companion when the latter was on the Asymptotic Giant Branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or CH subgiants. Combining radial-velocity measurements from HERMES and SALT high-resolution spectra with radial-velocity data from CORAVEL and CORALIE, we determine the orbital parameters of 27 systems. We also derive their masses by comparing their location in the Hertzsprung-Russell diagram with evolutionary models. We confirm that Ba dwarfs and CH subgiants are not at different evolutionary stages and have similar metallicities, despite their different names. Additionally, Ba giants appear significantly more massive than their main-sequence analogues. This is likely due to observational biases against the detection of hotter main-sequence post-mass-transfer objects. Combining our spectroscopic orbits with the Hipparcos astrometric data, we derive the orbital inclinations and the mass of the WD companion for four systems. Since this cannot be done for all systems in our sample yet (but should be with upcoming Gaia data releases), we also analyse the mass-function distribution of our binaries. We can model this distribution with very narrow mass distributions for the two components and random orbital orientation on the sky. Finally, based on BINSTAR evolutionary models, we suggest that the orbital evolution of low-mass Ba systems can be affected by a second phase of interaction along the Red Giant Branch of the Ba star, impacting on the eccentricities and periods of the giants.


title: Deep Learning Sentiment Analysis of Amazon.com Reviews and Ratings

abstract: Our study employs sentiment analysis to evaluate the compatibility of Amazon.com reviews with their corresponding ratings. Sentiment analysis is the task of identifying and classifying the sentiment expressed in a piece of text as being positive or negative. On e-commerce websites such as Amazon.com, consumers can submit their reviews along with a specific polarity rating. In some instances, there is a mismatch between the review and the rating. To identify the reviews with mismatched ratings we performed sentiment analysis using deep learning on Amazon.com product review data. Product reviews were converted to vectors using paragraph vector, which then was used to train a recurrent neural network with gated recurrent unit. Our model incorporated both semantic relationship of review text and product information. We also developed a web service application that predicts the rating score for a submitted review using the trained model and if there is a mismatch between predicted rating score and submitted rating score, it provides feedback to the reviewer.


title: A General Framework for the Semantics of Type Theory

abstract: We propose an abstract notion of a type theory to unify the semantics of various type theories including Martin-L\"{o}f type theory, two-level type theory and cubical type theory. We establish basic results in the semantics of type theory: every type theory has a bi-initial model; every model of a type theory has its internal language; the category of theories over a type theory is bi-equivalent to a full sub-2-category of the 2-category of models of the type theory.


title: A Survey of Distributed Consensus Protocols for Blockchain Networks

abstract: Since the inception of Bitcoin, cryptocurrencies and the underlying blockchain technology have attracted an increasing interest from both academia and industry. Among various core components, consensus protocol is the defining technology behind the security and performance of blockchain. From incremental modifications of Nakamoto consensus protocol to innovative alternative consensus mechanisms, many consensus protocols have been proposed to improve the performance of the blockchain network itself or to accommodate other specific application needs. In this survey, we present a comprehensive review and analysis on the state-of-the-art blockchain consensus protocols. To facilitate the discussion of our analysis, we first introduce the key definitions and relevant results in the classic theory of fault tolerance which help to lay the foundation for further discussion. We identify five core components of a blockchain consensus protocol, namely, block proposal, block validation, information propagation, block finalization, and incentive mechanism. A wide spectrum of blockchain consensus protocols are then carefully reviewed accompanied by algorithmic abstractions and vulnerability analyses. The surveyed consensus protocols are analyzed using the five-component framework and compared with respect to different performance metrics. These analyses and comparisons provide us new insights in the fundamental differences of various proposals in terms of their suitable application scenarios, key assumptions, expected fault tolerance, scalability, drawbacks and trade-offs. We believe this survey will provide blockchain developers and researchers a comprehensive view on the state-of-the-art consensus protocols and facilitate the process of designing future protocols.


title: Extension of the sub-grid-scale gradient model for compressible magnetohydrodynamics turbulent instabilities

abstract: Performing accurate large eddy simulations in compressible, turbulent magnetohydrodynamics is more challenging than in non-magnetized fluids due to the complex interplay between kinetic, magnetic and internal energy at different scales. Here we extend the sub-grid-scale gradient model, so far used in the momentum and induction equations, to account also for the unresolved scales in the energy evolution equation of a compressible ideal MHD fluid with a generic equation of state. We assess the model by considering box simulations of the turbulence triggered across a shear layer by the Kelvin-Helmholtz instability, testing cases where the small-scale dynamics cannot be fully captured by the resolution considered, such that the efficiency of the simulated dynamo effect depends on the resolution employed. This lack of numerical convergence is actually a currently common issue in several astrophysical problems, where the integral and fastest-growing-instability scales are too far apart to be fully covered numerically. We perform a-priori and a-posteriori tests of the extended gradient model. In the former, we find that, for many different initial conditions and resolutions, the gradient model outperforms other commonly used models in terms of correlation with the residuals coming from the filtering of a high-resolution run. In the second test, we show how a low-resolution run with the gradient model is able to quantitatively reproduce the evolution of the magnetic energy (the integrated value and the spectral distribution) coming from higher-resolution runs. This extension is the first step towards the implementation in relativistic magnetohydrodynamics.


title: Completely Unsupervised Speech Recognition By A Generative Adversarial Network Harmonized With Iteratively Refined Hidden Markov Models

abstract: Producing a large annotated speech corpus for training ASR systems remains difficult for more than 95% of languages all over the world which are low-resourced, but collecting a relatively big unlabeled data set for such languages is more achievable. This is why some initial effort have been reported on completely unsupervised speech recognition learned from unlabeled data only, although with relatively high error rates. In this paper, we develop a Generative Adversarial Network (GAN) to achieve this purpose, in which a Generator and a Discriminator learn from each other iteratively to improve the performance. We further use a set of Hidden Markov Models (HMMs) iteratively refined from the machine generated labels to work in harmony with the GAN. The initial experiments on TIMIT data set achieve an phone error rate of 33.1%, which is 8.5% lower than the previous state-of-the-art.


title: Thin-shell wormholes from Kiselev black holes

abstract: This paper discusses the theoretical construction of thin-shell wormholes from Kiselev black holes. We assume a barotropic equation of state for the exotic matter on the shell. While most of these wormholes are unstable to linearized radial perturbations, a limit argument is used to show that under certain conditions, stable solutions can be found.


title: Asymmetric Metal-dielectric Meta-cylinders and their Potential Applications from Engineering Scattering Patterns to Spatial Optical Signal Processing

abstract: We propose a novel type of bi-anisotropic hybrid metal-dielectric structure comprising dielectric and metallic cylindrical wedges wherein the composite meta-cylinder enables advanced control of electric, magnetic and magnetoelectric resonances. We establish a theoretical framework in which the electromagnetic response of this meta-atom is described through the electric and magneticmultipole moments. The complete dynamic polarizability tensor, expressed in a compact form, is derived as a function of the Mie scattering coefficients. Flexibility in the design makes the proposed meta-cylinder a viable candidate for various applications in the microscopic (single meta-atom) and macroscopic (metasurface) levels. We show that the highly versatile bi-anisotropic meta-atom is amenable to being designed for the desired electromagnetic response, such as electric dipole-free and zero/near-zero (backward and forward) scattering at the microscopic level. In addition, we show that the azimuthal asymmetry gives rise to normal polarizability components which are vital elements in synthesizing asymmetric Optical Transfer Function (OTF) at the macroscopic level.We conduct a precise inspection, from the microscopic to the macroscopic level, of the metasurface synthesis for emphasizing on the role of normal polarizability components for spatial optical signal processing. It is shown that this simple two-dimensional asymmetric meta-atom can perform first-order differentiation and edge detection at normal illumination.The results reported herein contribute toward improving the physical understanding of wave interaction with artificial materials composed of asymmetric elongated metal-dielectric inclusions and open the potential of its application in spatial signal and image processing.


title: Evaporating black-holes, wormholes, and vacuum polarisation: must they always conserve charge?

abstract: A careful examination of the fundamentals of electromagnetic theory shows that due to the underlying mathematical assumptions required for Stokes' Theorem, global charge conservation cannot be guaranteed in topologically non-trivial spacetimes. However, in order to break the charge conservation mechanism we must also allow the electromagnetic excitation fields D, H to possess a gauge freedom, just as the electromagnetic scalar and vector potentials phi and A do. This has implications for the treatment of electromagnetism in spacetimes where black holes both form and then evaporate, as well as extending the possibilities for treating vacuum polarisation. Using this gauge freedom of D, H we also propose an alternative to the accepted notion that a charge passing through a wormhole necessarily leads to an additional (effective) charge on the wormhole's mouth.


title: Software Engineering in Civic Tech: A Case Study about Code for Ireland

abstract: Civic grassroots have proven their ability to create useful and scalable software that addresses pressing social needs. Although software engineering plays a fundamental role in the process of creating civic technology, academic literature that analyses the software development processes of civic tech grassroots is scarce. This paper aims to advance the understanding of how civic grassroots tackle the different activities in their software development processes. In this study, we followed the formation of two projects in a civic tech group (Code for Ireland) seeking to understand how their development processes evolved over time, and how the group carried out their work in creating new technology. Our preliminary findings show that such groups are capable of setting up systematic software engineering processes that address software specification, development, validation, and evolution. While they were able to deliver software according to self-specified quality standards, the group has challenges in requirements specification, stakeholder engagement, and reorienting from development to product delivery. Software engineering methods and tools can effectively support the future of civic technologies and potentially improve their management, quality, and durability.


title: Augmenting adaptive immunity: progress and challenges in the quantitative engineering and analysis of adaptive immune receptor repertoires

abstract: The adaptive immune system is a natural diagnostic and therapeutic. It recognizes threats earlier than clinical symptoms manifest and neutralizes antigen with exquisite specificity. Recognition specificity and broad reactivity is enabled via adaptive B- and T-cell receptors: the immune receptor repertoire. The human immune system, however, is not omnipotent. Our natural defense system sometimes loses the battle to parasites and microbes and even turns against us in the case of cancer, autoimmune and inflammatory disease. A long-standing dream of immunoengineers has been, therefore, to mechanistically understand how the immune system sees, reacts and remembers antigens. Only very recently, experimental and computational methods have achieved sufficient quantitative resolution to start querying and engineering adaptive immunity with great precision. In specific, these innovations have been applied with the greatest fervency and success in immunotherapy, autoimmunity and vaccine design. The work here highlights advances, challenges and future directions of quantitative approaches which seek to advance the fundamental understanding of immunological phenomena, and reverse engineer the immune system to produce auspicious biopharmaceutical drugs and immunodiagnostics. Our review indicates that the merger of fundamental immunology, computational immunology and digital-biotechnology minimizes black box engineering, thereby advancing both immunological knowledge and as well immunoengineering methodologies.


title: Improved supervised learning methods for EoR parameters reconstruction

abstract: Within the next few years, the Square Kilometer Array (SKA) or one of its pathfinders will hopefully provide a detection of the 21-cm signal fluctuations from the Epoch of Reionization (EoR). Then, the main goal will be to accurately constrain the underlying astrophysical parameters. Currently, this is mainly done with Bayesian inference using Markov Chain Monte Carlo sampling. Recently, studies using neural networks trained to performed inverse modelling have shown interesting results. We build on these by improving the accuracy of the predictions using neural network and exploring other supervised learning methods: the kernel and ridge regressions. Based on a large training set of 21-cm power spectra, we compare the performances of these supervised learning methods. When using an un-noised signal as input, we improve on previous neural network accuracy by one order of magnitude and, using local ridge kernel regression, we gain another factor of a few. We then reach a rms prediction error of a few percents of the 1-sigma confidence level due to SKA thermal noise (as estimated with Bayesian inference). This last performance level requires optimizing the hyper-parameters of the method: how to do that perfectly in the case of an unknown signal remains an open question. For an input signal altered by a SKA-type thermal noise, our neural network recovers the astrophysical parameter values with an error within half of the 1$\sigma$ confidence level due to the SKA thermal noise. This accuracy improves to 10$\%$ of the 1$\sigma$ level when using the local ridge kernel regression (with optimized hyper-parameters). We are thus reaching a performance level where supervised learning methods are a viable alternative to determine the best-fit parameters values.


title: Enumeration degrees and non-metrizable topology

abstract: The enumeration degrees of sets of natural numbers can be identified with the degrees of difficulty of enumerating neighborhood bases of points in a universal second-countable $T_0$-space (e.g. the $\omega$-power of the Sierpi\'nski space). Hence, every represented second-countable $T_0$-space determines a collection of enumeration degrees. For instance, Cantor space captures the total degrees, and the Hilbert cube captures the continuous degrees by definition. Based on these observations, we utilize general topology (particularly non-metrizable topology) to establish a classification theory of enumeration degrees of sets of natural numbers.


title: A systematic TMRT observational study of Galactic $^{12}$C/$^{13}$C ratios from Formaldehyde

abstract: We present observations of the C-band $1_{10}-1_{11}$ (4.8 GHz) and Ku-band $2_{11}-2_{12}$ (14.5 GHz) K-doublet lines of H$_2$CO and the C-band $1_{10}-1_{11}$ (4.6 GHz) line of H$_2$$^{13}$CO toward a large sample of Galactic molecular clouds, through the Shanghai Tianma 65-m radio telescope (TMRT). Our sample with 112 sources includes strong H$_2$CO sources from the TMRT molecular line survey at C-band and other known H$_2$CO sources. All three lines are detected toward 38 objects (43 radial velocity components) yielding a detection rate of 34\%. Complementary observations of their continuum emission at both C- and Ku-bands were performed. Combining spectral line parameters and continuum data, we calculate the column densities, the optical depths and the isotope ratio H$_2$$^{12}$CO/H$_2$$^{13}$CO for each source. To evaluate photon trapping caused by sometimes significant opacities in the main isotopologue's rotational mm-wave lines connecting our measured K-doublets, and to obtain $^{12}$C/$^{13}$C abundance ratios, we used the RADEX non-LTE model accounting for radiative transfer effects. This implied the use of the new collision rates from \citet{Wiesenfeld2013}. Also implementing distance values from trigonometric parallax measurements for our sources, we obtain a linear fit of $^{12}$C/$^{13}$C = (5.08$\pm$1.10)D$_{GC}$ + (11.86$\pm$6.60), with a correlation coefficient of 0.58. D$_{GC}$ refers to Galactocentric distances. Our $^{12}$C/$^{13}$C ratios agree very well with the ones deduced from CN and C$^{18}$O but are lower than those previously reported on the basis of H$_2$CO, tending to suggest that the bulk of the H$_2$CO in our sources was formed on dust grain mantles and not in the gas phase.


title: Vector SIMP dark matter with approximate custodial symmetry

abstract: We consider a novel scenario for Vector Strongly Interacting Massive Particle (VSIMP) dark matter with local $SU(2)_X\times U(1)_{Z'}$ symmetry in the dark sector. Similarly to the Standard Model, after the dark symmetry is broken spontaneously by the VEVs of dark Higgs fields, the approximate custodial symmetry determines comparable but split masses for $SU(2)_X$ gauge bosons. In this model, we show that $U(1)_{Z'}$-charged gauge boson of $SU(2)_X$ ($X_\pm$) becomes a natural candidate for SIMP dark matter, annihilating through $3\rightarrow 2$ or forbidden $2\rightarrow 2$ annihilations due to gauge self-interactions. On the other hand, the $U(1)_{Z'}$-neutral gauge boson of $SU(2)_X$ achieves the kinetic equilibrium of dark matter through a gauge kinetic mixing between $U(1)_{Z'}$ and Standard Model. We present the parameter space for the correct relic density in our model and discuss in detail the current constraints and projections from colliders and direct detection experiments.


title: Polarized positron beams via intense two-color laser pulses

abstract: Generation of ultrarelativistic polarized positrons during interaction of an ultrarelativistic electron beam with a counterpropagating two-color petawatt laser pulse is investigated theoretically. Our Monte Carlo simulation based on a semi-classical model, incorporates photon emissions and pair productions, using spin-resolved quantum probabilities in the local constant field approximation, and describes the polarization of electrons and positrons for the pair production and photon emission processes, as well as the classical spin precession in-between. The main reason of the polarization is shown to be the spin-asymmetry of the pair production process in strong external fields, combined with the asymmetry of the two-color laser field. Employing a feasible scenario, we show that highly polarized positron beams, with a polarization degree of $\zeta\approx 60\%$, can be produced in a femtosecond time scale, with a small angular divergence, $\sim 74$ mrad, and high density $\sim 10^{14}$ cm$^{-3}$. The laser-driven positron source, along with laser wakefield acceleration, may pave the way to small scale facilities for high energy physics studies.


title: Uncertainty relation under information measurement and feedback control

abstract: Here, we investigate the uncertainty of dynamical observables in classical systems manipulated by repeated measurements and feedback control; the precision should be enhanced in the presence of an external controller but limited by the amount of information obtained from the measurements. We prove that the entropy production and the information quantity constrain from below the fluctuation of arbitrary observables that are antisymmetric under time reversal. The information term is the sum of the mutual entropy production and the Kullback--Leibler divergence, which characterises the irreversibility of the measurement outcomes. The result holds for finite observation times and for both continuous- and discrete-time systems. We apply the derived relation to study the precision of a flashing Brownian ratchet.


title: Convex Sobolev inequalities related to unbalanced optimal transport

abstract: We study the behaviour of various Lyapunov functionals (relative entropies) along the solutions of a family of nonlinear drift-diffusion-reaction equations coming from statistical mechanics and population dynamics. These equations can be viewed as gradient flows over the space of Radon measures equipped with the Hellinger-Kantorovich distance. The driving functionals of the gradient flows are not assumed to be geodesically convex or semi-convex. We prove new isoperimetric-type functional inequalities, allowing us to control the relative entropies by their productions, which yields the exponential decay of the relative entropies.


title: Thermodynamic coupling rule for far-from-equilibrium systems

abstract: The energy conversion efficiency of far-from-equilibrium systems is generally limited by irreversible thermodynamic fluxes that make contact with different heat baths. For complex systems, the states of the maximum efficiency and the minimum entropy production are usually not equivalent. Here we show that the proper adjustments of the interaction between the energy and matter currents offer some important criteria for the performance characterizations of thermal agents, regardless of the system types and transition protocols. The universal thermodynamic coupling rule plays a critical role in irreversible processes. A double quantum dot system is applied to demonstrate that the performances of heat engines or refrigrators can be enhanced by suitably adjusting the coupling strength between thermodynamic fluxes.


title: Quark and gluon condensates in QCD reevaluated

abstract: We compute the quark and gluon condensates in $QCD$ with $N$ colors and $N_f$ flavors based on the renormalization group equations and on the knowledge of a single scale $\Lambda_{QCD}$. For $N=3$ and $N_f=3$ our findings are in the good range for $0.2\leq\Lambda_{QCD}\leq 0.36$ and in excellent agreement with the results in the literature from sum rules for a value $\Lambda_{QCD}=0.28$.


title: ADMM for Block Circulant Model Predictive Control

abstract: This paper deals with model predictive control problems for large scale dynamical systems with cyclic symmetry. Based on the properties of block circulant matrices, we introduce a complex-valued coordinate transformation that block diagonalizes and truncates the original finite-horizon optimal control problem. Using this coordinate transformation, we develop a modified alternating direction method of multipliers (ADMM) algorithm for general constrained quadratic programs with block circulant blocks. We test our modified algorithm in two different simulated examples and show that our coordinate transformation significantly increases the computation speed.


title: Revisiting Adversarial Autoencoder for Unsupervised Word Translation with Cycle Consistency and Improved Training

abstract: Adversarial training has shown impressive success in learning bilingual dictionary without any parallel data by mapping monolingual embeddings to a shared space. However, recent work has shown superior performance for non-adversarial methods in more challenging language pairs. In this work, we revisit adversarial autoencoder for unsupervised word translation and propose two novel extensions to it that yield more stable training and improved results. Our method includes regularization terms to enforce cycle consistency and input reconstruction, and puts the target encoders as an adversary against the corresponding discriminator. Extensive experimentations with European, non-European and low-resource languages show that our method is more robust and achieves better performance than recently proposed adversarial and non-adversarial approaches.


title: The experiment paradox in physics

abstract: Modern physics is founded on two mainstays: mathematical modelling and empirical verification. These two assumptions are prerequisite for the objectivity of scientific discourse. Here we show, however, that they are contradictory, leading to the `experiment paradox'. We reveal that any experiment performed on a physical system is - by necessity - invasive and thus establishes inevitable limits to the accuracy of any mathematical model. We track its manifestations in both classical and quantum physics and show how it is overcome `in practice' via the concept of environment. We argue that the scientific pragmatism ordains two methodological principles of compressibility and stability.

categories: physics.hist-ph gr-qc hep-th physics.class-ph quant-ph


title: Three-dimensional local anisotropy of velocity fluctuations in the solar wind

abstract: We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations ({\delta}B) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations ({\delta}V) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of {\delta}B in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes.


title: The Innovative Design of the PANDA Barrel DIRC

abstract: The Barrel DIRC of the PANDA experiment at FAIR will cleanly separate pions from kaons for the physics program of PANDA. Innovative solutions for key components of the detector sitting in the strong magnetic field of the compact PANDA target spectrometer as well as two reconstruction methods were developed in an extensive prototype program. The technical design and present results from the test beam campaigns at the CERN PS in 2017 and 2018 are discussed.


title: Temperature dependence of the probability of "small heating" and total losses of ucns on the surface of fomblin oils of different molecular mass

abstract: We measured the temperature dependence of the probability of small heating and total losses of UCNs on the PFPE Fomblin Y surface with various molecular masses Mw=2800, 3300, 6500 amu in the temperature range of 100-300 K. The probability of small heating sharply decreases with increasing Mw and decreasing temperature. The probability of total loss weakly decreases with decreasing temperature and takes the minimum value at Mw=3300 amu. As this oil provides a homogeneous surface with minimal probabilities of small heating and total losses of UCNs, it is the preferred candidate for experiments on measuring the neutron lifetime.


title: Quark masses, CKM angles and Lepton Flavour Universality violation

abstract: A properly defined and suitably broken $U(2)$ flavour symmetry leads to successful quantitative relations between quark mass ratios and CKM angles. At the same time the intrinsic distinction introduced by $U(2)$ between the third and the first two families of quarks and leptons may support anomalies in charged and neutral current semi-leptonic $B$-decays of the kind tentatively observed in current flavour experiments. We show how this is possible by the exchange of the $(3,1)_{2/3}$ vector leptoquark in two $U(2)$-models with significantly different values of Lepton Flavour Universality violation, observable in foreseen experiments.


title: Ab initio and nuclear inelastic scattering studies of Fe$_3$Si/GaAs heterostructures

abstract: The structure and dynamical properties of the Fe$_3$Si/GaAs(001) interface are investigated by density functional theory and nuclear inelastic scattering measurements. The stability of four different atomic configurations of the Fe$_3$Si/GaAs multilayers is analyzed by calculating the formation energies and phonon dispersion curves. The differences in charge density, magnetization, and electronic density of states between the configurations are examined. Our calculations unveil that magnetic moments of the Fe atoms tend to align in a plane parallel to the interface, along the [110] direction of the Fe$_3$Si crystallographic unit cell. In some configurations, the spin polarization of interface layers is larger than that of bulk Fe$_3$Si. The effect of the interface on element-specific and layer-resolved phonon density of states is discussed. The Fe-partial phonon density of states measured for the Fe$_3$Si layer thickness of three monolayers is compared with theoretical results obtained for each interface atomic configuration. The best agreement is found for one of the configurations with a mixed Fe-Si interface layer, which reproduces the anomalous enhancement of the phonon density of states below 10 meV


title: ASAP: Architecture Search, Anneal and Prune

abstract: Automatic methods for Neural Architecture Search (NAS) have been shown to produce state-of-the-art network models. Yet, their main drawback is the computational complexity of the search process. As some primal methods optimized over a discrete search space, thousands of days of GPU were required for convergence. A recent approach is based on constructing a differentiable search space that enables gradient-based optimization, which reduces the search time to a few days. While successful, it still includes some noncontinuous steps, e.g., the pruning of many weak connections at once. In this paper, we propose a differentiable search space that allows the annealing of architecture weights, while gradually pruning inferior operations. In this way, the search converges to a single output network in a continuous manner. Experiments on several vision datasets demonstrate the effectiveness of our method with respect to the search cost and accuracy of the achieved model. Specifically, with $0.2$ GPU search days we achieve an error rate of $1.68\%$ on CIFAR-10.


title: From top-hat masking to smooth transitions: P-filter and its application to polarized microwave sky maps

abstract: In CMB science, the simplest idea to remove a contaminated sky region is to multiply the sky map with a mask that is 0 for the contaminated region and 1 elsewhere, which is also called a top-hat masking. Although it is easy to use, such top-hat masking is known to suffer from various leakage problems. Therefore, we want to extend the top-hat masking to a series of semi-analytic functions called the P-filters. Most importantly, the P-filters can seamlessly realize the core idea of masking in CMB science, and, meanwhile, guarantee continuity up to the first derivative everywhere. The P-filters can significantly reduce many leakage problems without additional cost, including the leakages due to low-, high-, and band-pass filtering, and the E-to-E, B-to-B, B-to-E, and E-to-B leakages. The workings of the P-filter are illustrated by using the WMAP and Planck polarization sky maps. By comparison to the corresponding WMAP/Planck masks, we show that the P-filter performs much better than top-hat masking, and meanwhile, has the potential to supersede the principal idea of masking in CMB science. Compared to mask apodization, the P-filter is ``outward'', that tends to make proper use of the region that was marked as 0; whereas apodization is ``inward'', that always kills more signal in the region marked as 1.


title: Local study of the insulating quantum kagome antiferromagnets YCu3(OH)6OxCl3-x (x=0,1/3)

abstract: The quantum kagome antiferromagnets YCu3(OH)6OxCl3-x (x=0,1/3) are produced using a unified solid state synthesis route for polycrystalline samples. From structural refinements based on neutron diffraction data, we clarify the structure of the Y3Cu9(OH)18OCl8 (x=1/3) compound and provide a revised chemical formula. We use muon spin relaxation, as a local probe of magnetism, to investigate the exotic low temperature properties in the two compounds. In agreement with the low temperature neutron diffraction data, we find no evidence for long range ordering in both materials but they exhibit distinct ground states: while disordered static magnetism develops in the x=0 compound, we conclude on the stabilization of a quantum spin liquid in the x=1/3 one, since the local fields remain fully dynamical. Our findings are in contrast to previous reports based on thermodynamical measurements only. We then discuss their origin on the basis of structural details and specific heat measurements. In particular, the x=1/3 compound appears to realize an original spatially anisotropic kagome model.


title: Weighted Reservoir Sampling from Distributed Streams

abstract: We consider message-efficient continuous random sampling from a distributed stream, where the probability of inclusion of an item in the sample is proportional to a weight associated with the item. The unweighted version, where all weights are equal, is well studied, and admits tight upper and lower bounds on message complexity. For weighted sampling with replacement, there is a simple reduction to unweighted sampling with replacement. However, in many applications the stream has only a few heavy items which may dominate a random sample when chosen with replacement. Weighted sampling \textit{without replacement} (weighted SWOR) eludes this issue, since such heavy items can be sampled at most once. In this work, we present the first message-optimal algorithm for weighted SWOR from a distributed stream. Our algorithm also has optimal space and time complexity. As an application of our algorithm for weighted SWOR, we derive the first distributed streaming algorithms for tracking \textit{heavy hitters with residual error}. Here the goal is to identify stream items that contribute significantly to the residual stream, once the heaviest items are removed. Residual heavy hitters generalize the notion of $\ell_1$ heavy hitters and are important in streams that have a skewed distribution of weights. In addition to the upper bound, we also provide a lower bound on the message complexity that is nearly tight up to a $\log(1/\epsilon)$ factor. Finally, we use our weighted sampling algorithm to improve the message complexity of distributed $L_1$ tracking, also known as count tracking, which is a widely studied problem in distributed streaming. We also derive a tight message lower bound, which closes the message complexity of this fundamental problem.


title: Dynamics of a class A nonlinear mirror mode-locked laser

abstract: Using a delay differential equation model we study theoretically the dynamics of a unidirectional class-A ring laser with a nonlinear amplifying loop mirror. We perform linear stability analysis of the CW regimes in the large delay limit and demonstrate that these regimes can be destabilized via modulational and Turing-type instabilities, as well as by an instability leading to the appearance of square-waves. We investigate the formation of square-waves and mode-locked pulses in the system. We show that mode-locked pulses are asymmetric with exponential decay of the trailing edge in positive time and faster-than-exponential (super-exponential) decay of the leading edge in negative time. We discuss asymmetric interaction of these pulses leading to a formation of harmonic mode-locked regimes.


title: A multiple criteria nominal classification method in a web-based platform: Demonstration in a case of recruitment for the Portuguese Army

abstract: CAT-SD (CATegorization by Similarity-Dissimilarity) is a recently developed method for handling nominal classification problems in the context of Multiple Criteria Decision Aiding (MCDA). This paper describes the design and implementation of this method, as well as an application dealing with a recruitment process in the Special Forces of the Portuguese Army. In addition, it proposes interaction protocols to elicit the preference parameters of the method to facilitate the construction of a decision model when the analyst guides the decision maker. CAT-SD has been implemented in DecSpace, a user-friendly on-line platform for supporting decision aiding processes using one or more MCDA methods. The study related to the Portuguese Army Special Forces recruitment presents and demonstrates how these protocols and a tool like DecSpace can facilitate the process of applying the method in real-world scenarios.


title: A note on Cunningham's algorithm for matroid intersection

abstract: In the matroid intersection problem, we are given two matroids of rank $r$ on a common ground set $E$ of $n$ elements and the goal is to find the maximum set that is independent in both matroids. In this note, we show that Cunningham's algorithm for matroid intersection can be implemented to use $O(nr\log^2(r))$ independent oracle calls.


title: Results on heavy ion physics at LHCb

abstract: In the last years, the \lhcb experiment established itself as an important contributor to heavy ion physics by exploiting some of its specific features. Production of particles, notably heavy flavour states, can be studied in p-p, p-Pb and Pb-Pb collisions at LHC energies in the forward rapidity region (pseudorapidity between 2 and 5), providing measurements which are highly complementary to the other LHC experiments. Moreover, owing to its forward geometry, the detector is also well suited to study fixed-target collisions, obtained by impinging the LHC beams on gas targets with different mass numbers. In this configuration, p-A collisions can be studied at the relatively unexplored scale of sqrt(sNN) ~ 100 GeV, also providing valuable inputs to cosmic ray physics. An overview of the measurements obtained so far by the LHCb ion program is presented.


title: Computing the Mostar index in networks with applications to molecular graphs

abstract: Recently, a bond-additive topological descriptor, named as the Mostar index, has been introduced as a measure of peripherality in networks. For a connected graph $G$, the Mostar index is defined as $Mo(G) = \sum_{e=uv \in E(G)} |n_u(e) - n_v(e)|$, where for an edge $e=uv$ we denote by $n_u(e)$ the number of vertices of $G$ that are closer to $u$ than to $v$ and by $n_v(e)$ the number of vertices of $G$ that are closer to $v$ than to $u$. In this paper, we generalize the definition of the Mostar index to weighted graphs and prove that the Mostar index of a weighted graph can be computed in terms of Mostar indices of weighted quotient graphs. As a consequence, we show that the Mostar index of a benzenoid system can be computed in sub-linear time with respect to the number of vertices. Finally, our method is applied to some benzenoid systems and to a fullerene patch.


title: Neutrino-Dark Matter Portals

abstract: The nature of dark matter is one of the open problems of the Standard Model of particle physics. Despite the great experimental efforts, we have not yet found a positive signal of its interactions with ordinary matter. One possible explanation would be that the dark matter particle is primarily coupled to another elusive particle, neutrinos. In this work we study this possibility with several realisations.


title: The influence of building interactions on seismic and elastic surface waves

abstract: We outline some recent research advances on the control of elastic waves in thin and thick plates, that have occurred since the large scale experiment [Phys. Rev. Lett. 112, 133901, 2014] that demonstrated significant interaction of surface seismic waves with holes structuring sedimentary soils at the meter scale. We further investigate the seismic wave trajectories in soils structured with buildings. A significant substitution of soils by inclusions, acting as foundations, raises the question of the effective dynamic properties of these structured soils. Buildings, in the case of perfect elastic conditions for both soil and buildings, are shown to interact and strongly influence elastic surface waves; such site-city seismic interactions were pointed out in [Bulletin of Seismological Society of America 92, 794-811, 2002], and we investigate a variety of scenarios to illustrate the variety of behaviours possible.


title: Bakry-\'Emery Ricci curvature of doubly warped product of weighted spaces

abstract: We introduce a notion of doubly warped product of weighted graphs that is consistent with the doubly warped product in the Riemannian setting. We establish various discrete Bakry-\'Emery Ricci curvature-dimension bounds for such warped products in terms of the curvature of the constituent graphs. This requires deliberate analysis of the quadratic forms involved, prompting the introduction of some crucial notions such as curvature saturation at a vertex. In the spirit of being thorough and to provide a frame of reference, we also introduce the $\left(R_1,R_2\right)$-doubly warped products of smooth measure spaces and establish $\N$-Bakry-\'Emery Ricci curvature (lower) bounds thereof in terms of those of the factors. At the end of these notes, we present examples and demonstrate applications of warped products with some toy models.


title: Thermal counting statistics in an atomic two-mode squeezed vacuum state

abstract: We measure the population distribution in one of the atomic twin beams generated by four-wave mixing in an optical lattice. Although the produced two-mode squeezed vacuum state is pure, each individual mode is described as a statistical mixture. We confirm the prediction that the particle number follows an exponential distribution when only one spatio-temporal mode is selected. We also show that this distribution accounts well for the contrast of an atomic Hong--Ou--Mandel experiment. These experiments constitute an important validation of our twin beam source in view of a future test of a Bell inequalities.


title: The $b_1$ resonance in coupled $\pi\omega$, $\pi\phi$ scattering from lattice QCD

abstract: We present the first lattice QCD calculation of coupled $\pi\omega$ and $\pi\phi$ scattering, incorporating coupled $S$ and $D$-wave $\pi\omega$ in $J^P=1^+$. Finite-volume spectra in three volumes are determined via a variational analysis of matrices of two-point correlation functions, computed using large bases of operators resembling single-meson, two-meson and three-meson structures, with the light-quark mass corresponding to a pion mass of $m_\pi \approx 391$ MeV. Utilizing the relationship between the discrete spectrum of finite-volume energies and infinite-volume scattering amplitudes, we find a narrow axial-vector resonance ($J^{PC}=1^{+-}$), the analogue of the $b_1$ meson, with mass $m_{R}\approx1380$ MeV and width $\Gamma_{R}\approx 91$ MeV. The resonance is found to couple dominantly to $S$-wave $\pi\omega$, with a much-suppressed coupling to $D$-wave $\pi\omega$, and a negligible coupling to $\pi\phi$ consistent with the `OZI rule'. No resonant behavior is observed in $\pi\phi$, indicating the absence of a putative low-mass $Z_s$ analogue of the $Z_c$ claimed in $\pi J/\psi$. In order to minimally present the contents of a unitary three-channel scattering matrix, we introduce an $n$-channel generalization of the traditional two-channel Stapp parameterization.


title: Diabetes Mellitus Forecasting Using Population Health Data in Ontario, Canada

abstract: Leveraging health administrative data (HAD) datasets for predicting the risk of chronic diseases including diabetes has gained a lot of attention in the machine learning community recently. In this paper, we use the largest health records datasets of patients in Ontario,Canada. Provided by the Institute of Clinical Evaluative Sciences (ICES), this database is age, gender and ethnicity-diverse. The datasets include demographics, lab measurements,drug benefits, healthcare system interactions, ambulatory and hospitalizations records. We perform one of the first large-scale machine learning studies with this data to study the task of predicting diabetes in a range of 1-10 years ahead, which requires no additional screening of individuals.In the best setup, we reach a test AUC of 80.3 with a single-model trained on an observation window of 5 years with a one-year buffer using all datasets. A subset of top 15 features alone (out of a total of 963) could provide a test AUC of 79.1. In this paper, we provide extensive machine learning model performance and feature contribution analysis, which enables us to narrow down to the most important features useful for diabetes forecasting. Examples include chronic conditions such as asthma and hypertension, lab results, diagnostic codes in insurance claims, age and geographical information.