title: D2KE: From Distance to Kernel and Embedding

abstract: For many machine learning problem settings, particularly with structured inputs such as sequences or sets of objects, a distance measure between inputs can be specified more naturally than a feature representation. However, most standard machine models are designed for inputs with a vector feature representation. In this work, we consider the estimation of a function $f:\mathcal{X} \rightarrow \R$ based solely on a dissimilarity measure $d:\mathcal{X}\times\mathcal{X} \rightarrow \R$ between inputs. In particular, we propose a general framework to derive a family of \emph{positive definite kernels} from a given dissimilarity measure, which subsumes the widely-used \emph{representative-set method} as a special case, and relates to the well-known \emph{distance substitution kernel} in a limiting case. We show that functions in the corresponding Reproducing Kernel Hilbert Space (RKHS) are Lipschitz-continuous w.r.t. the given distance metric. We provide a tractable algorithm to estimate a function from this RKHS, and show that it enjoys better generalizability than Nearest-Neighbor estimates. Our approach draws from the literature of Random Features, but instead of deriving feature maps from an existing kernel, we construct novel kernels from a random feature map, that we specify given the distance measure. We conduct classification experiments with such disparate domains as strings, time series, and sets of vectors, where our proposed framework compares favorably to existing distance-based learning methods such as $k$-nearest-neighbors, distance-substitution kernels, pseudo-Euclidean embedding, and the representative-set method.


title: Axisymmetric magnetic modes of neutron stars having mixed poloidal and toroidal magnetic fields

abstract: We calculate axisymmetric magnetic modes of a neutron star possessing a mixed poloidal and toroidal magnetic field, where the toroidal field is assumed to be proportional to a dimensionless parameter $\zeta_0$. Here, we assume an isentropic structure for the neutron star and consider no effects of rotation. Ignoring the equilibrium deformation due to the magnetic field, we employ a polytrope of the index $n=1$ as the background model for our modal analyses. For the mixed poloidal and toroidal magnetic field with $\zeta_0\not=0$, axisymmetric spheroidal and toroidal modes are coupled. We compute axisymmetric spheroidal and toroidal magnetic modes as a function of the parameter $\zeta_0$ from $0$ to $\sim 1$ for the surface field strengths $B_S=10^{14}$G and $10^{15}$G. We find that the frequency $\omega$ of the magnetic modes decreases with increasing $\zeta_0$. We also find that the frequency of the spheroidal magnetic modes is almost exactly proportional to $B_S$ for $\zeta_0\lesssim 1$ but that this proportionality holds only when $\zeta_0\ll 1$ for the toroidal magnetic modes. The wave patterns of the spheroidal magnetic modes and toroidal magnetic modes are not strongly affected by the coupling so long as $\zeta_0\lesssim 1$. We find no unstable modes having $\omega^2<0$.


title: Fluctuations in Arctic Sea Ice Extent: Comparing Observations and Climate Models

abstract: The fluctuation statistics of the observed sea ice extent during the satellite era are compared with model output from CMIP5 models using a multi-fractal time series method. The two robust features of the observations are that on annual to bi-annual time scales the ice extent exhibits white noise structure and there is a decadal scale trend associated with the decay of the ice cover. It is shown that (i) there is a large inter-model variability in the time scales extracted from the models, (ii) none of the models exhibit the decadal time scales found in the satellite observations, (iii) {5} of the 21 models examined exhibit the observed white noise structure, and (iv) the multi-model ensemble mean exhibits neither the observed white noise structure nor the observed decadal trend. It is proposed that the observed fluctuation statistics produced by this method serve as an appropriate test bed for modeling studies.


title: Pathway towards programmable wave anisotropy in cellular metamaterials

abstract: In this work, we provide a proof-of-concept experimental demonstration of the wave control capabilities of cellular metamaterials endowed with populations of tunable electromechanical resonators. Each independently tunable resonator comprises a piezoelectric patch and a resistor-inductor shunt, and its resonant frequency can be seamlessly re-programmed without interfering with the cellular structure's default properties. We show that, by strategically placing the resonators in the lattice domain and by deliberately activating only selected subsets of them, chosen to conform to the directional features of the beamed wave response, it is possible to override the inherent wave anisotropy of the cellular medium. The outcome is the establishment of tunable spatial patterns of energy distillation resulting in a non-symmetric correction of the wavefields.


title: Vertex nomination: The canonical sampling and the extended spectral nomination schemes

abstract: Suppose that one particular block in a stochastic block model is of interest, but block labels are only observed for a few of the vertices in the network. Utilizing a graph realized from the model and the observed block labels, the vertex nomination task is to order the vertices with unobserved block labels into a ranked nomination list with the goal of having an abundance of interesting vertices near the top of the list. There are vertex nomination schemes in the literature, including the optimally precise canonical nomination scheme~$\mathcal{L}^C$ and the consistent spectral partitioning nomination scheme~$\mathcal{L}^P$. While the canonical nomination scheme $\mathcal{L}^C$ is provably optimally precise, it is computationally intractable, being impractical to implement even on modestly sized graphs. With this in mind, an approximation of the canonical scheme---denoted the {\it canonical sampling nomination scheme} $\mathcal{L}^{CS}$---is introduced; $\mathcal{L}^{CS}$ relies on a scalable, Markov chain Monte Carlo-based approximation of $\mathcal{L}^{C}$, and converges to $\mathcal{L}^{C}$ as the amount of sampling goes to infinity. The spectral partitioning nomination scheme is also extended to the {\it extended spectral partitioning nomination scheme}, $\mathcal{L}^{EP}$, which introduces a novel semisupervised clustering framework to improve upon the precision of $\mathcal{L}^P$. Real-data and simulation experiments are employed to illustrate the precision of these vertex nomination schemes, as well as their empirical computational complexity. Keywords: vertex nomination, Markov chain Monte Carlo, spectral partitioning, Mclust MSC[2010]: 60J22, 65C40, 62H30, 62H25


title: A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows VI: Dilute Disperse Gas-Particle Multiphase System

abstract: In this paper, a unified gas kinetic scheme for multiphase dilute gas-particle system is proposed. The UGKS multiphase (UGKS-M) is a finite volume method, which captures flow physics in the regimes from collisionless multispecies transport to the two-fluid hydrodynamic Navier-Stokes (NS) solution with the variation of Knudsen number, and from granular flow regime to dusty gas dynamics with the variation of Stokes number. The main reason for preserving the multiscale nature in UGKS-M is mainly coming from the direct modeling of the flow physics in the scales of discrete cell size and time step, where the ratio of the time step over the particle collision time determines flow behavior in different regimes. For the particle phase, the integral solution of the kinetic equation is used in the construction of the numerical flux, which takes into account the particle transport, collision, and acceleration. The gas phase, which is assumed to be in the continuum flow regime, evolves numerically by the gas kinetic scheme (GKS), which is a subset of the UGKS for the Navier-Stokes (NS) solutions. The interaction between the gas and particle phase is calculated based on a velocity space mapping method, which solves accurately the kinetic acceleration process. The stability of UGKS-M is determined by the CFL condition only. With the inclusion of the material temperature evolution equation of solid particles, the UGKS-M conserves the total mass, momentum, and energy for the whole multiphase system. In the numerical tests, the UGKS-M shows good multiscale property in capturing the particle trajectory crossing (PTC), particle wall reflecting phenomena, and vortex-induced segregation of inertial particles under different Stokes numbers. The scheme is also applied to simulate shock induced fluidization problem and the simulation results agree well with experimental measurement.


title: Disjoint Multi-task Learning between Heterogeneous Human-centric Tasks

abstract: Human behavior understanding is arguably one of the most important mid-level components in artificial intelligence. In order to efficiently make use of data, multi-task learning has been studied in diverse computer vision tasks including human behavior understanding. However, multi-task learning relies on task specific datasets and constructing such datasets can be cumbersome. It requires huge amounts of data, labeling efforts, statistical consideration etc. In this paper, we leverage existing single-task datasets for human action classification and captioning data for efficient human behavior learning. Since the data in each dataset has respective heterogeneous annotations, traditional multi-task learning is not effective in this scenario. To this end, we propose a novel alternating directional optimization method to efficiently learn from the heterogeneous data. We demonstrate the effectiveness of our model and show performance improvements on both classification and sentence retrieval tasks in comparison to the models trained on each of the single-task datasets.


title: Superconvergent recovery of Raviart--Thomas mixed finite elements on triangular grids

abstract: For the second lowest-order Raviart--Thomas mixed method, we prove that the canonical interpolant and finite element solution for the vector variable in elliptic problems are superclose in the $H(\text{div})$-norm on mildly structured meshes, where most pairs of adjacent triangles form approximate parallelograms. We then develop a family of postprocessing operators for Raviart--Thomas mixed elements on triangular grids by using the idea of local least squares fittings. Super-approximation property of the postprocessing operators for the lowest and second lowest order Raviart--Thomas elements is proved under mild conditions. Combining the supercloseness and super-approximation results, we prove that the postprocessed solution superconverges to the exact solution in the $L^2$-norm on mildly structured meshes.


title: Indirect unitarity violation entangled with matter effects in reactor antineutrino oscillations

abstract: If finite but tiny masses of the three active neutrinos are generated via the canonical seesaw mechanism with three heavy sterile neutrinos, the 3\times 3 Pontecorvo-Maki-Nakagawa-Sakata neutrino mixing matrix V will not be exactly unitary. This kind of indirect unitarity violation can be probed in a precision reactor antineutrino oscillation experiment, but it may be entangled with terrestrial matter effects as both of them are very small. We calculate the probability of \overline{\nu}_e \to \overline{\nu}_e oscillations in a good analytical approximation, and find that, besides the zero-distance effect, the effect of unitarity violation is always smaller than matter effects, and their entanglement does not appear until the next-to-leading-order oscillating terms are taken into account. Given a 20-kiloton JUNO-like liquid scintillator detector, we reaffirm that terrestrial matter effects should not be neglected but indirect unitarity violation makes no difference, and demonstrate that the experimental sensitivities to the neutrino mass ordering and a precision measurement of \theta_{12} and \Delta_{21} \equiv m^2_2 - m^2_1 are robust.


title: Structural-elastic determination of the mechanical lifetime of biomolecules

abstract: The lifetime of protein domains and ligand-receptor complexes under force is crucial for mechanosensitive functions, while many aspects of how force affects the lifetime still remain poorly understood. Here, we report a new analytical expression of the force-dependent molecular lifetime to understand transitions overcoming a single barrier. Unlike previous models derived in the framework of Kramers theory that requires a presumed one-dimensional free energy landscape, our model is derived based on the structural-elastic properties of molecules which is not restricted by the shape and dimensionality of the underlying free energy landscape. Importantly, the parameters of this model provide direct information of the structural-elastic features of the molecules between the transition and the native states. We demonstrate the applications of this model by applying it to explain complex force-dependent lifetime data for several molecules reported in recent experiments, and predict the structural-elastic properties of the transition states of these molecules.


title: Destination Choice Game: A Spatial Interaction Theory on Human Mobility

abstract: With remarkable significance in migration prediction, global disease mitigation, urban planning and many others, an arresting challenge is to predict human mobility fluxes between any two locations. A number of methods have been proposed against the above challenge, including the gravity model, the intervening opportunity model, the radiation model, the population-weighted opportunity model, and so on. Despite their theoretical elegance, all models ignored an intuitive and important ingredient in individual decision about where to go, that is, the possible congestion on the way and the possible crowding in the destination. Here we propose a microscopic mechanism underlying mobility decisions, named destination choice game (DCG), which takes into account the crowding effects resulted from spatial interactions among individuals. In comparison with the state-of-the-art models, the present one shows more accurate prediction on mobility fluxes across wide scales from intracity trips to intercity travels, and further to internal migrations. The well-known gravity model is proved to be the equilibrium solution of a degenerated DCG neglecting the crowding effects in the destinations.


title: DESlib: A Dynamic ensemble selection library in Python

abstract: DESlib is an open-source python library providing the implementation of several dynamic selection techniques. The library is divided into three modules: (i) \emph{dcs}, containing the implementation of dynamic classifier selection methods (DCS); (ii) \emph{des}, containing the implementation of dynamic ensemble selection methods (DES); (iii) \emph{static}, with the implementation of static ensemble techniques. The library is fully documented (documentation available online on Read the Docs), has a high test coverage (codecov.io) and is part of the scikit-learn-contrib supported projects. Documentation, code and examples can be found on its GitHub page: https://github.com/scikit-learn-contrib/DESlib.


title: Median Shapes

abstract: We introduce and begin to explore the mean and median of finite sets of shapes represented as integral currents. The median can be computed efficiently in practice, and we focus most of our theoretical and computational attention on medians. We consider questions on the existence and regularity of medians. While the median might not exist in all cases, we show that a mass-regularized median is guaranteed to exist. When the input shapes are modeled by integral currents with shared boundaries in codimension $1$, we show that the median is guaranteed to exist, and is contained in the \emph{envelope} of the input currents. On the other hand, we show that medians can be \emph{wild} in this setting, and smooth inputs can generate non-smooth medians. For higher codimensions, we show that \emph{books} are minimizing for a finite set of $1$-currents in $\Bbb{R}^3$ with shared boundaries. As part of this proof, we present a new result in graph theory---that \emph{cozy} graphs are \emph{comfortable}---which should be of independent interest. Further, we show that regular points on the median have book-like tangent cones in this case. From the point of view of computation, we study the median shape in the settings of a finite simplicial complex. When the input shapes are represented by chains of the simplicial complex, we show that the problem of finding the median shape can be formulated as an integer linear program. This optimization problem can be solved as a linear program in practice, thus allowing one to compute median shapes efficiently. We provide open source code implementing our methods, which could also be used by anyone to experiment with ideas of their own. The software could be accessed at https://github.com/tbtraltaa/medianshape.


title: Anti-screening of the Galileon force around a disk center hole

abstract: The Vainshtein mechanism is known as an efficient way of screening the fifth force around a matter source in modified gravity. This has been verified mainly in highly symmetric matter configurations. To study how the Vainshtein mechanism works in a less symmetric setup, we numerically solve the scalar field equation around a disk with a hole at its center in the cubic Galileon theory. We find, surprisingly, that the Galileon force is enhanced, rather than suppressed, in the vicinity of the hole. This anti-screening effect is larger for a thinner, less massive disk with a smaller hole. At this stage our setup is only of academic interest and its astrophysical consequences are unclear, but this result implies that the Vainshtein screening mechanism around less symmetric matter configurations is quite nontrivial.


title: Abstract Family-based Model Checking using Modal Featured Transition Systems: Preservation of CTL* (Extended Version)

abstract: Variational systems allow effective building of many custom variants by using features (configuration options) to mark the variable functionality. In many of the applications, their quality assurance and formal verification are of paramount importance. Family-based model checking allows simultaneous verification of all variants of a variational system in a single run by exploiting the commonalities between the variants. Yet, its computational cost still greatly depends on the number of variants (often huge). In this work, we show how to achieve efficient family-based model checking of CTL* temporal properties using variability abstractions and off-the-shelf (single-system) tools. We use variability abstractions for deriving abstract family-based model checking, where the variability model of a variational system is replaced with an abstract (smaller) version of it, called modal featured transition system, which preserves the satisfaction of both universal and existential temporal properties, as expressible in CTL*. Modal featured transition systems contain two kinds of transitions, termed may and must transitions, which are defined by the conservative (over-approximating) abstractions and their dual (under-approximating) abstractions, respectively. The variability abstractions can be combined with different partitionings of the set of variants to infer suitable divide-and-conquer verification plans for the variational system. We illustrate the practicality of this approach for several variational systems.


title: Hadronic Paschen-Back effect

abstract: We find a novel phenomenon induced by the interplay between a strong magnetic field and finite orbital angular momenta in hadronic systems, which is analogous to the Paschen-Back effect observed in the field of atomic physics. This effect allows the wave functions to drastically deform. We discuss anisotropic decay from the deformation as a possibility to measure the strength of the magnetic field in heavy-ion collision at LHC, RHIC and SPS, which has not experimentally been measured. As an example we investigate charmonia with a finite orbital angular momentum in a strong magnetic field. We calculate the mass spectra and mixing rates. To obtain anisotropic wave functions, we apply the cylindrical Gaussian expansion method, where the Gaussian bases to expand the wave functions have different widths along transverse and longitudinal directions in the cylindrical coordinate.

categories: hep-ph hep-ex nucl-ex nucl-th physics.atom-ph


title: A velocity-space adaptive unified gas kinetic scheme for continuum and rarefied flows

abstract: In this paper, a unified gas kinetic scheme with adaptive velocity space (AUGKS) for multiscale flow transport will be developed. In near-equilibrium flow regions, particle distribution function is close to the Chapman-Enskog expansion and can be formulated with a continuous velocity space, where only macroscopic conservative variables are updated. With the emerging of non-equilibrium effects, the AUGKS automatically switches to a discrete velocity space to follow the evolution of particle distribution function. Based on the Chapman-Enskog expansion, a criterion is proposed in this paper to quantify the intensity of non-equilibrium effects and is used for the continuous-discrete velocity space transformation. Following the scale-dependent local evolution solution, the AUGKS presents the discretized gas dynamic equations directly on the cell size and time step scales, i.e., the so-called direct modeling method. As a result, the scheme is able to capture the cross-scale flow physics from particle transport to hydrodynamic wave propagation, and provides a continuous variation of solutions from the Boltzmann to the Navier-Stokes. Under the unified framework, different from conventional DSMC-NS hybrid method, the AUGKS does not need a buffer zone to match up kinetic and hydrodynamic solutions. Instead, a continuous and discrete particle velocity space is naturally connected, which is feasible for the numerical simulations with unsteadiness or complex geometries. Compared with the asymptotic preserving (AP) methods which solves kinetic equations uniformly over the entire flow field with discretized velocity space, the current velocity-space adaptive unified scheme speeds up the computation and reduces the memory requirement in multiscale flow problems, and maintains the equivalent accuracy. The AUGKS provides an effective tool for non-equilibrium flow simulations.


title: Coproducts in brane topology

abstract: We extend the loop product and the loop coproduct to the mapping space from the $k$-dimensional sphere, or more generally from any $k$-manifold, to a $k$-connected space with finite dimensional rational homotopy group, $k\geq 1$. The key to extending the loop coproduct is the fact that the embedding $M\rightarrow M^{S^{k-1}}$ is of "finite codimension" in a sense of Gorenstein spaces. Moreover, we prove the associativity, commutativity, and Frobenius compatibility of them.


title: A Framework for Input-Output Analysis of Wall-Bounded Shear Flows

abstract: We propose a framework to understand input-output amplification properties of non- linear partial differential equation (PDE) models of wall-bounded shear flows, which are spatially invariant in one coordinate (e.g., streamwise-constant plane Couette flow). Our methodology is based on the notion of dissipation inequalities in control theory. In particular, we consider flows with body and other forcings, for which we study the input- to-output properties, including energy growth, worst-case disturbance amplification, and stability to persistent disturbances. The proposed method can be applied to a large class of flow configurations as long as the base flow is described by a polynomial. This includes many examples in both channel flows and pipe flows, e.g., plane Couette flow, and Hagen-Poiseuille flow. The methodology we use is numerically implemented as the solution of a (convex) optimization problem. We use the framework to study input-output amplification mechanisms in rotating Couette flow, plane Couette flow, plane Poiseuille flow, and Hagen-Poiseuille flow. In addition to showing that the application of the proposed framework leads to results that are consistent with theoretical and experimental amplification scalings obtained in the literature through linearization around the base flow, we demonstrate that the stability bounds to persistent forcings can be used as a means to predict transition to turbulence in wall-bounded shear flows.


title: High-Dimensional Entanglement in States with Positive Partial Transposition

abstract: Genuine high-dimensional entanglement, i.e. the property of having a high Schmidt number, constitutes a resource in quantum communication, overcoming limitations of low-dimensional systems. States with a positive partial transpose (PPT), on the other hand, are generally considered weakly entangled, as they can never be distilled into pure entangled states. This naturally raises the question, whether high Schmidt numbers are possible for PPT states. Volume estimates suggest that optimal, i.e. linear, scaling in local dimension should be possible, albeit without providing an insight into the possible slope. We provide the first explicit construction of a family of PPT states that achieves linear scaling in local dimension and we prove that random PPT states typically share this feature. Our construction also allows us to answer a recent question by Chen et al. on the existence of PPT states whose Schmidt number increases by an arbitrarily large amount upon partial transposition. Finally, we link the Schmidt number to entangled sub-block matrices of a quantum state. We use this connection to prove that quantum states invariant under partial transposition on the smaller of their two subsystems cannot have maximal Schmidt number. This generalizes a well-known result by Kraus et al. We also show that the Schmidt number of absolutely PPT states cannot be maximal, contributing to an open problem in entanglement theory.


title: Dihedral Group, 4-Torsion on an Elliptic Curve, and a Peculiar Eigenform Modulo 4

abstract: We work out a non-trivial example of lifting a so-called weak eigenform to a true, characteristic 0 eigenform. The weak eigenform is closely related to Ramanujan's tau function whereas the characteristic 0 eigenform is attached to an elliptic curve defined over ${\mathbb Q}$. We produce the lift by showing that the coefficients of the initial, weak eigenform (almost all) occur as traces of Frobenii in the Galois representation on the 4-torsion of the elliptic curve. The example is remarkable as the initial form is known not to be liftable to any characteristic 0 eigenform of level 1. We use this example as illustrating certain questions that have arisen lately in the theory of modular forms modulo prime powers. We give a brief survey of those questions.


title: Paraphrasing Complex Network: Network Compression via Factor Transfer

abstract: Many researchers have sought ways of model compression to reduce the size of a deep neural network (DNN) with minimal performance degradation in order to use DNNs in embedded systems. Among the model compression methods, a method called knowledge transfer is to train a student network with a stronger teacher network. In this paper, we propose a novel knowledge transfer method which uses convolutional operations to paraphrase teacher's knowledge and to translate it for the student. This is done by two convolutional modules, which are called a paraphraser and a translator. The paraphraser is trained in an unsupervised manner to extract the teacher factors which are defined as paraphrased information of the teacher network. The translator located at the student network extracts the student factors and helps to translate the teacher factors by mimicking them. We observed that our student network trained with the proposed factor transfer method outperforms the ones trained with conventional knowledge transfer methods.


title: Singularly perturbed forward-backward stochastic differential equations: application to the optimal control of bilinear systems

abstract: We study linear-quadratic stochastic optimal control problems with bilinear state dependence for which the underlying stochastic differential equation (SDE) consists of slow and fast degrees of freedom. We show that, in the same way in which the underlying dynamics can be well approximated by a reduced order effective dynamics in the time scale limit (using classical homogenziation results), the associated optimal expected cost converges in the time scale limit to an effective optimal cost. This entails that we can well approximate the stochastic optimal control for the whole system by the reduced order stochastic optimal control, which is clearly easier to solve because of lower dimensionality. The approach uses an equivalent formulation of the Hamilton-Jacobi-Bellman (HJB) equation, in terms of forward-backward SDEs (FBSDEs). We exploit the efficient solvability of FBSDEs via a least squares Monte Carlo algorithm and show its applicability by a suitable numerical example.


title: M4CD: A Robust Change Detection Method for Intelligent Visual Surveillance

abstract: In this paper, we propose a robust change detection method for intelligent visual surveillance. This method, named M4CD, includes three major steps. Firstly, a sample-based background model that integrates color and texture cues is built and updated over time. Secondly, multiple heterogeneous features (including brightness variation, chromaticity variation, and texture variation) are extracted by comparing the input frame with the background model, and a multi-source learning strategy is designed to online estimate the probability distributions for both foreground and background. The three features are approximately conditionally independent, making multi-source learning feasible. Pixel-wise foreground posteriors are then estimated with Bayes rule. Finally, the Markov random field (MRF) optimization and heuristic post-processing techniques are used sequentially to improve accuracy. In particular, a two-layer MRF model is constructed to represent pixel-based and superpixel-based contextual constraints compactly. Experimental results on the CDnet dataset indicate that M4CD is robust under complex environments and ranks among the top methods.


title: How to Simulate Patchy Particles

abstract: Patchy particles is the name given to a large class of systems of mesoscopic particles characterized by a repulsive core and a discrete number of short-range and highly directional interaction sites. Numerical simulations have contributed significantly to our understanding of the behaviour of patchy particles, but, although simple in principle, advanced simulation techniques are often required to sample the low temperatures and long time scales associated with their self-assembly behaviour. In this work we review the most popular simulation techniques that have been used to study patchy particles, with a special focus on Monte Carlo methods. We cover many of the tools required to simulate patchy systems, from interaction potentials to biased moves, cluster moves, and free energy methods. The review is complemented by an educationally-oriented Monte Carlo computer code that implements all the techniques described in the text to simulate a well-known tetrahedral patchy particle model.


title: Adaptive importance sampling with forward-backward stochastic differential equations

abstract: We describe an adaptive importance sampling algorithm for rare events that is based on a dual stochastic control formulation of a path sampling problem. Specifically, we focus on path functionals that have the form of cumulate generating functions, which appear relevant in the context of, e.g.~molecular dynamics, and we discuss the construction of an optimal (i.e. minimum variance) change of measure by solving a stochastic control problem. We show that the associated semi-linear dynamic programming equations admit an equivalent formulation as a system of uncoupled forward-backward stochastic differential equations that can be solved efficiently by a least squares Monte Carlo algorithm. We illustrate the approach with a suitable numerical example and discuss the extension of the algorithm to high-dimensional systems.


title: Craig Interpolation and Access Interpolation with Clausal First-Order Tableaux

abstract: We develop foundations for computing Craig interpolants and similar intermediates of two given formulas with first-order theorem provers that construct clausal tableaux. Provers that can be understood in this way include efficient machine-oriented systems based on calculi of two families: goal-oriented like model elimination and the connection method, and bottom-up like the hyper tableau calculus. The presented method for Craig-Lyndon interpolation involves a lifting step where terms are replaced by quantified variables, similar as known for resolution-based interpolation, but applied to a differently characterized ground formula and proven correct more abstractly on the basis of Herbrand's theorem, independently of a particular calculus. Access interpolation is a recent form of interpolation for database query reformulation that applies to first-order formulas with relativized quantifiers and constrains the quantification patterns of predicate occurrences. It has been previously investigated in the framework of Smullyan's non-clausal tableaux. Here, in essence, we simulate these with the more machine-oriented clausal tableaux through structural constraints that can be ensured either directly by bottom-up tableau construction methods or, for closed clausal tableaux constructed with arbitrary calculi, by postprocessing with restructuring transformations.


title: Quantum Fluctuation and Dissipation in Holographic Theories: A Unifying Study Scheme

abstract: Motivated by the wide range of applicability of the fluctuation and dissipation phenomena in non-equilibrium systems, we provide a universal study scheme for the dissipation of the energy and the corresponding Brownian motion analysis of massive particles due to quantum and thermal fluctuations in a wide class of strongly coupled quantum field theories. The underlying reason for the existence of such unified study scheme, is that our analytic methods turn out to heavily depend on the order of the Bessel functions $\nu$, describing the string fluctuations attached to the particle. Different values of the order are associated to different theories. The two-point function of the fluctuations exhibits two different late time behaviors, depending purely on the value of the order of Bessel functions. We then find that the coefficients and observables associated with the stochastic motion at zero and finite temperature, depend on the scales of the theory through powers of the order $\nu$. Moreover, the fluctuation-dissipation theorem is verified from the bulk perspective to be universally satisfied for the whole class of theories. Finally, we show that the analysis of certain types of Dp-brane fluctuations can be mapped one-to-one to the string fluctuations and therefore the stochastic brane observables can be read from the string ones. In the closing remarks we demonstrate how our analysis accommodates known results as special cases and provide more applications.


title: On ranks of polynomials

abstract: Let $V$ be a vector space over a field $k, P:V\to k, d\geq 3$. We show the existence of a function $C(r,d)$ such that $rank (P)\leq C(r,d)$ for any field $k,char (k)>d$, a finite-dimensional $k$-vector space $V$ and a polynomial $P:V\to k$ of degree $d$ such that $rank(\partial P/\partial t)\leq r$ for all $t\in V-0$. Our proof of this theorem is based on the application of results on Gowers norms for finite fields $k$. We don't know a direct proof in the case when $k=\mathbb C$.


title: A class of fully nonlinear equations

abstract: In this paper we consider a class of fully nonlinear equations which cover the equation introduced by S. Donaldson a decade ago and the equation introduced by Gursky-Streets recently. We solve the equation with uniform weak $C^2$ estimates, which hold for degenerate case.


title: Convolutional Neural Networks over Control Flow Graphs for Software Defect Prediction

abstract: Existing defects in software components is unavoidable and leads to not only a waste of time and money but also many serious consequences. To build predictive models, previous studies focus on manually extracting features or using tree representations of programs, and exploiting different machine learning algorithms. However, the performance of the models is not high since the existing features and tree structures often fail to capture the semantics of programs. To explore deeply programs' semantics, this paper proposes to leverage precise graphs representing program execution flows, and deep neural networks for automatically learning defect features. Firstly, control flow graphs are constructed from the assembly instructions obtained by compiling source code; we thereafter apply multi-view multi-layer directed graph-based convolutional neural networks (DGCNNs) to learn semantic features. The experiments on four real-world datasets show that our method significantly outperforms the baselines including several other deep learning approaches.


title: PlayeRank: data-driven performance evaluation and player ranking in soccer via a machine learning approach

abstract: The problem of evaluating the performance of soccer players is attracting the interest of many companies and the scientific community, thanks to the availability of massive data capturing all the events generated during a match (e.g., tackles, passes, shots, etc.). Unfortunately, there is no consolidated and widely accepted metric for measuring performance quality in all of its facets. In this paper, we design and implement PlayeRank, a data-driven framework that offers a principled multi-dimensional and role-aware evaluation of the performance of soccer players. We build our framework by deploying a massive dataset of soccer-logs and consisting of millions of match events pertaining to four seasons of 18 prominent soccer competitions. By comparing PlayeRank to known algorithms for performance evaluation in soccer, and by exploiting a dataset of players' evaluations made by professional soccer scouts, we show that PlayeRank significantly outperforms the competitors. We also explore the ratings produced by {\sf PlayeRank} and discover interesting patterns about the nature of excellent performances and what distinguishes the top players from the others. At the end, we explore some applications of PlayeRank -- i.e. searching players and player versatility --- showing its flexibility and efficiency, which makes it worth to be used in the design of a scalable platform for soccer analytics.


title: The Depoissonisation quintet: Rice-Poisson-Mellin-Newton-Laplace

abstract: This paper is devoted to the Depoissonnisation process, which is central in various analyses of the AofA domain. We first recall the two possible paths that may be used in this process. The first path, called here the Depoissonisation path, is better studied and is proven to apply in any practical situation; however, it often uses technical tools, that are not so easy to deal with. Moreover, the various results are scattered in the litterature, and the most recent results are not well known within the AofA domain. The present paper gathers in Section 2 all these results in a survey style. The second path, called here the Rice-Mellin path, is less often used within the AofA domain. It is often very easy to apply, but it needs a tameness condition, which appears {\em a priori} to be quite restrictive, and is not deeply studiedin the litterature. In Section 3, the paper precisely describes the Rice-Mellin path, together with its tameness condition, in a survey style, too. Finally, in Section 4, the paper presents original results for the Rice-Mellin path: it exhibits a framework, of practical use, where the tameness condition is proven to hold. It then proves that the Rice-Mellin path is both of easy and practical use : even though (much?) less general than the Depoissonisation path, it is easier to apply.


title: A second order scheme for a Robin boundary condition in random walk algorithms

abstract: Random Walk (RW) is a common numerical tool for modeling the Advection-Diffusion equation. In this work, we develop a second order scheme for incorporating a heterogeneous reaction (i.e., a Robin boundary condition) in the RW model. In addition, we apply the approach in two test cases. We compare the second order scheme with the first order one as well as with analytical and other numerical solution. We show that the new scheme can reduce the computational error significantly, relative to the first order scheme. This reduction comes at no additional computational cost.


title: American Options in the Hobson-Rogers Model

abstract: In this article, we consider a risky asset $X$ for which evolution follows a model proposed by D.G. Hobson and L.C.G. Rogers\cite{HR98}. We assume that the volatility of $X$ depends on the ratio of the present value and the exponentially weighted average of the past value. Using the Markovian modelling of the enlarged two-dimensional process, we show that, for the American put option with $X$ as the underlying asset, the continuation region and the stopped region are separated a striking curve . This striking curve lies between the two striking curves from the basic BSM model, yet is {\em not} monotone.


title: Regularity of entropy, geodesic currents and entropy at infinity

abstract: In this work, we introduce the notion of entropy at infinity, and define a wide class of noncompact manifolds with negative curvature, those which admit a critical gap between entropy at infinity and topological entropy. We call them strongly positively recurrent manifolds (SPR), and provide many examples. We show that dynamically, they behave as compact manifolds. In particular, they admit a finite measure of maximal entropy. Using the point of view of currents at infinity, we show that on these SPR manifolds the topological entropy of the geodesic flow varies in a C 1 -way along (uniformly) C 1 -perturbations of the metric. This result generalizes former work of Katok (1982) and Katok-Knieper-Weiss (1991) in the compact case.


title: The Soul Conjecture in Alexandrov Geometry in dimension 4

abstract: In this paper, we prove the Soul Conjecture in Alexandrov geometry in dimension $4$, i.e. if $X$ is a complete non-compact $4$-dimensional Alexandrov space of non-negative curvature and positive curvature around one point, then a soul of $X$ is a point.


title: Heat Coulomb Blockade of One Ballistic Channel

abstract: Quantum mechanics and Coulomb interaction dictate the behavior of small circuits. The thermal implications cover fundamental topics from quantum control of heat to quantum thermodynamics, with prospects of novel thermal machines and an ineluctably growing influence on nanocircuit engineering. Experimentally, the rare observations thus far include the universal thermal conductance quantum and heat interferometry. However, evidences for many-body thermal effects paving the way to markedly different heat and electrical behaviors in quantum circuits remain wanting. Here we report on the observation of the Coulomb blockade of electronic heat flow from a small metallic circuit node, beyond the widespread Wiedemann-Franz law paradigm. We demonstrate this thermal many-body phenomenon for perfect (ballistic) conduction channels to the node, where it amounts to the universal suppression of precisely one quantum of conductance for the transport of heat, but none for electricity. The inter-channel correlations that give rise to such selective heat current reduction emerge from local charge conservation, in the floating node over the full thermal frequency range ($\lesssim$temperature$\times k_\mathrm{B}/h$). This observation establishes the different nature of the quantum laws for thermal transport in nanocircuits.

categories: cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el quant-ph


title: Idempotent geometry in generic algebras

abstract: Using the syzygy method, established in our earlier paper, we characterize the combinatorial stratification of the variety of two-dimensional real generic algebras. We show that there exist exactly three different homotopic types of such algebras and relate this result to potential applications and known facts from qualitative theory of quadratic ODEs. The genericity condition is crucial. For example, the idempotent geometry in Clifford algebras or Jordan algebras of Clifford type is very different: such algebras always contain nontrivial submanifolds of idempotents.


title: Breeze: Sharing Biofeedback Through Wearable Technologies

abstract: Digitally presenting physiological signals as biofeedback to users raises awareness of both body and mind. This paper describes the effectiveness of conveying a physiological signal often overlooked for communication: breathing. We present the design and development of digital breathing patterns and their evaluation along three output modalities: visual, audio, and haptic. We also present Breeze, a wearable pendant placed around the neck that measures breathing and sends biofeedback in real-time. We evaluated how the breathing patterns were interpreted in a fixed environment and gathered qualitative data on the wearable device's design. We found that participants intentionally modified their own breathing to match the biofeedback, as a technique for understanding the underlying emotion. Our results describe how the features of the breathing patterns and the feedback modalities influenced participants' perception. We include guidelines and suggested use cases, such as Breeze being used by loved ones to increase connectedness and empathy.


title: The algebraic de Rham realization of the elliptic polylogarithm via the Poincar\'e bundle

abstract: In this paper, we describe the algebraic de Rham realization of the elliptic polylogarithm for arbitrary families of elliptic curves in terms of the Poincar\'e bundle. Our work builds on previous work of Scheider and generalizes results of Bannai-Kobayashi-Tsuji and Scheider. As an application, we compute the de Rham Eisenstein classes explicitly in terms of certain algebraic Eisenstein series.


title: Quantum phase transitions of a two-leg bosonic ladder in an artificial gauge field

abstract: We consider a two leg bosonic ladder in a $U(1)$ gauge field with both interleg hopping and interleg repulsion. As a function of the flux, the interleg interaction converts the commensurate-incommensurate transition from the Meissner to a Vortex phase, into an Ising-type of transition towards a density wave phase. A disorder point is also found after which the correlation functions develop a damped sinusoid behavior signaling a melting of the vortex phase. We discuss the differences on the phase diagram for attractive and repulsive interleg interaction. In particular, we show how repulsion favors the Meissner phase at low-flux and a phase with a second incommensuration in the correlation functions for intermediate flux, leading to a richer phase diagram than in the case of interleg attraction. The effect of the temperature on the chiral current is also discussed.


title: Interacting hadron resonance gas model in K-matrix formalism

abstract: An extension of Hadron Resonance Gas (HRG) model is constructed to include interactions using relativistic virial expansion of partition function. The non-interacting part of the expansion contains all the stable baryons and mesons and the interacting part contains all the higher mass resonances which decay into two stable hadrons. The virial coefficients are related to the phase shifts which are calculated using K-matrix formalism in the present work. We have calculated various thermodynamics quantities like pressure, energy density, and entropy density of the system. A comparison of thermodynamic quantities with non interacting HRG model, calculated using the same number of hadrons, shows that the results of above formalism are larger. A good agreement between equation of state calculated in K-matrix formalism and lattice QCD simulations is observed. Specifically the lattice QCD calculated interaction measure is well described in our formalism. We have also calculated second order fluctuations and correlations of conserved charges in K-matrix formalism. We observe a good agreement of second order fluctuations and baryon-strangeness correlation with lattice data below the cross-over temperature.


title: The syntomic realization of the elliptic polylogarithm via the Poincar\'e bundle

abstract: We give an explicit description of the syntomic elliptic polylogarithm on the universal elliptic curve over the ordinary locus of the modular curve in terms of certain $p$-adic analytic moment functions associated to Katz' two-variable $p$-adic Eisenstein measure. The present work generalizes previous results of Bannai-Kobayashi-Tsuji and Bannai-Kings on the syntomic Eisenstein classes.


title: High magnetocaloric efficiency of a NiFe/NiCu/CoFe/MnIr multilayer in a small magnetic field

abstract: The isothermal magnetic entropy changes are studied in Ni80Fe20/Ni67Cu33/Co90Cu10/Mn80Ir20 stacks at temperatures near the Curie point of the Ni67Cu33 spacer by applying magnetic fields in a few tens of Oersted. Such low fields were sufficient for toggling magnetic moments in the soft ferromagnetic (FM) layer (Ni80Fe20). It is found out that this switching provides the magnetic entropy change, which is up to 20 times larger than that achievable in a single Ni67Cu33 film subjected to such low fields. Our finding holds promise to be utilized in the magnetocaloric devices that would be based on FM/PM/FM heterostructures and would operate with moderate magnetic fields.


title: Angular clustering of point sources at 150 MHz in the TGSS survey

abstract: We study the angular clustering of point sources in The GMRT (Giant Meter Wave Telescope) Sky Survey (TGSS). The survey at 150 MHz with delta > -53.5 degrees has a sky coverage of 3.6 pi steradians, i.e., 90% of the whole sky. We created subsamples by applying different total flux thresholds limit (S >> 5 sigma) for good completeness and measured the angular correlation function omega(theta) of point sources at large scales ( >= 1 degree). We find that the amplitude of angular clustering is higher for brighter subsamples, this indicates that higher threshold flux samples are hosted by massive halos and cluster strongly: this conclusions is based on the assumption that the redshift distribution of sources does not change with flux and this is supported by models of radio sources. We compare our results with other low-frequency studies of clustering of point sources and verify that the amplitude of clustering varies with the flux limit. We quantify this variation as a power law dependence of the amplitude of correlation function with the flux limit. This dependence can be used to estimate foreground contamination due to clustering of point sources for low frequency HI intensity mapping surveys for studying the epoch of reionisation.


title: Algebraic torus actions on contact manifolds

abstract: We prove the LeBrun-Salamon Conjecture in low dimensions. More precisely, we show that a contact Fano manifold X of dimension 2n+1 that has reductive automorphism group of rank at least n-2 is necessarily homogeneous. This implies that any positive quaternion-Kahler manifold of real dimension at most 16 is necessarily a symmetric space, one of the Wolf spaces. A similar result about contact Fano manifolds of dimension at most 9 with reductive automorphism group also holds. The main difficulty in approaching the conjecture is how to recognize a homogeneous space in an abstract variety. We contribute to such problem in general, by studying the action of algebraic torus on varieties and exploiting Bialynicki-Birula decomposition and equivariant Riemann-Roch theorems. From the point of view of T-varieties (that is, varieties with a torus action), our result is about high complexity T-manifolds. The complexity here is at most 1/2(dim X+5) with dim X arbitrarily high, but we require this special (contact) structure of X. Previous methods for studying T-varieties in general usually only apply for complexity at most 2 or 3.


title: Evidence of Intermediate-Scale Energy Spectrum Anisotropy of Cosmic Rays E$\geq$10$^{19.2}$ eV with the Telescope Array Surface Detector

abstract: An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).


title: Zero-Knowledge Password Policy Check from Lattices

abstract: Passwords are ubiquitous and most commonly used to authenticate users when logging into online services. Using high entropy passwords is critical to prevent unauthorized access and password policies emerged to enforce this requirement on passwords. However, with current methods of password storage, poor practices and server breaches have leaked many passwords to the public. To protect one's sensitive information in case of such events, passwords should be hidden from servers. Verifier-based password authenticated key exchange, proposed by Bellovin and Merrit (IEEE S\&P, 1992), allows authenticated secure channels to be established with a hash of a password (verifier). Unfortunately, this restricts password policies as passwords cannot be checked from their verifier. To address this issue, Kiefer and Manulis (ESORICS 2014) proposed zero-knowledge password policy check (ZKPPC). A ZKPPC protocol allows users to prove in zero knowledge that a hash of the user's password satisfies the password policy required by the server. Unfortunately, their proposal is not quantum resistant with the use of discrete logarithm-based cryptographic tools and there are currently no other viable alternatives. In this work, we construct the first post-quantum ZKPPC using lattice-based tools. To this end, we introduce a new randomised password hashing scheme for ASCII-based passwords and design an accompanying zero-knowledge protocol for policy compliance. Interestingly, our proposal does not follow the framework established by Kiefer and Manulis and offers an alternate construction without homomorphic commitments. Although our protocol is not ready to be used in practice, we think it is an important first step towards a quantum-resistant privacy-preserving password-based authentication and key exchange system.


title: Using Longitudinal Targeted Maximum Likelihood Estimation in Complex Settings with Dynamic Interventions

abstract: Longitudinal targeted maximum likelihood estimation (LTMLE) has very rarely been used to estimate dynamic treatment effects in the context of time-dependent confounding affected by prior treatment when faced with long follow-up times, multiple time-varying confounders, and complex associational relationships simultaneously. Reasons for this include the potential computational burden, technical challenges, restricted modeling options for long follow-up times, and limited practical guidance in the literature. However, LTMLE has desirable asymptotic properties, i.e. it is doubly robust, and can yield valid inference when used in conjunction with machine learning. We use a topical and sophisticated question from HIV treatment research to show that LTMLE can be used successfully in complex realistic settings and compare results to competing estimators. Our example illustrates the following practical challenges common to many epidemiological studies 1) long follow-up time (30 months), 2) gradually declining sample size 3) limited support for some intervention rules of interest 4) a high-dimensional set of potential adjustment variables, increasing both the need and the challenge of integrating appropriate machine learning methods 5) consideration of collider bias. Our analyses, as well as simulations, shed new light on the application of LTMLE in complex and realistic settings: we show that (i) LTMLE can yield stable and good estimates, even when confronted with small samples and limited modeling options; (ii) machine learning utilized with a small set of simple learners (if more complex ones can't be fitted) can outperform a single, complex model, which is tailored to incorporate prior clinical knowledge; (iii) performance can vary considerably depending on interventions and their support in the data, and therefore critical quality checks should accompany every LTMLE analysis.


title: Longitudinal phase space reconstruction simulation studies using a novel X-band transverse deflecting structure at the SINBAD facility at DESY

abstract: A transverse deflecting structure (TDS) is a well-known device for characterizing the longitudinal properties of an electron bunch in a linear accelerator. The standard use of such a cavity involves streaking the bunch along a transverse axis and analysing the image on a screen downstream to find the bunch length and slice properties along the other transverse axis. A novel X-band deflecting structure, which will allow the polarization of the deflecting field to be adjusted, is currently being designed as part of a collaboration between CERN, DESY and PSI. This new design will allow bunches to be streaked at any transverse angle within the cavity, which will open up the possibility of new measurement techniques, which could be combined to characterize the 6D phase space distribution of bunches. In this paper, a method is presented for reconstructing the longitudinal phase space distribution of bunches by using the TDS in combination with a dipole. Simulations of this technique for the SINBAD-ARES beamline are presented and the key limitations related to temporal resolution and induced energy spread are discussed.


title: The step Sidorenko property and non-norming edge-transitive graphs

abstract: Sidorenko's Conjecture asserts that every bipartite graph H has the Sidorenko property, i.e., a quasirandom graph minimizes the density of H among all graphs with the same edge density. We study a stronger property, which requires that a quasirandom multipartite graph minimizes the density of H among all graphs with the same edge densities between its parts; this property is called the step Sidorenko property. We show that many bipartite graphs fail to have the step Sidorenko property and use our results to show the existence of a bipartite edge-transitive graph that is not weakly norming; this answers a question of Hatami [Israel J. Math. 175 (2010), 125-150].


title: TIFR Near Infrared Imaging Camera-II on the 3.6-m Devasthal Optical Telescope

abstract: TIFR Near Infrared Imaging Camera-II is a closed-cycle Helium cryo-cooled imaging camera equipped with a Raytheon 512 x 512 pixels InSb Aladdin III Quadrant focal plane array having sensitivity to photons in the 1-5 microns wavelength band. In this paper, we present the performance of the camera on the newly installed 3.6-m Devasthal Optical Telescope (DOT) based on the calibration observations carried out during 2017 May 11-14 and 2017 October 7-31. After the preliminary characterization, the camera has been released to the Indian and Belgian astronomical community for science observations since 2017 May. The camera offers a field-of-view of ~86.5 arcsec x 86.5 arcsec on the DOT with a pixel scale of 0.169 arcsec. The seeing at the telescope site in the near-infrared bands is typically sub-arcsecond with the best seeing of ~0.45 arcsec realized in the near-infrared K-band on 2017 October 16. The camera is found to be capable of deep observations in the J, H and K bands comparable to other 4-m class telescopes available world-wide. Another highlight of this camera is the observational capability for sources up to Wide-field Infrared Survey Explorer (WISE) W1-band (3.4 microns) magnitudes of 9.2 in the narrow L-band (nbL; lambda_{cen} ~3.59 microns). Hence, the camera could be a good complementary instrument to observe the bright nbL-band sources that are saturated in the Spitzer-Infrared Array Camera ([3.6] <= 7.92 mag) and the WISE W1-band ([3.4] <= 8.1 mag). Sources with strong polycyclic aromatic hydrocarbon (PAH) emission at 3.3 microns are also detected. Details of the observations and estimated parameters are presented in this paper.


title: Three-body correlations in direct reactions: Example of $^{6}$Be populated in $(p,n)$ reaction

abstract: The $^{6}$Be continuum states were populated in the charge-exchange reaction $^1$H($^{6}$Li,$^{6}$Be)$n$ collecting very high statistics data ($\sim 5 \times 10^6$ events) on the three-body $\alpha$+$p$+$p$ correlations. The $^{6}$Be excitation energy region below $\sim 3$ MeV is considered, where the data are dominated by contributions from the $0^+$ and $2^+$ states. It is demonstrated how the high-statistics few-body correlation data can be used to extract detailed information on the reaction mechanism. Such a derivation is based on the fact that highly spin-aligned states are typically populated in the direct reactions.


title: Cycles of Singularities appearing in the Resolution Problem in positive Characteristic

abstract: We present a hypersurface singularity in positive characteristic which is defined by a purely inseparable power series, and a sequence of point blowups so that, after applying the blowups to the singularity, the same type of singularity reappears after the last blowup, with just certain exponents of the defining power series shifted upwards. The construction hence yields a cycle. Iterating this cycle leads to an infinite increase of the residual order of the defining power series. This disproves a theorem claimed by Moh about the stability of the residual order under sequences of blowups. It is not a counter-example to the resolution in positive characteristic since larger centers are also permissible and prevent the phenomenon from happening.


title: SIR epidemics and vaccination on random graphs with clustering

abstract: In this paper we consider SIR epidemics on random graphs with clustering. To incorporate group structure of the underlying social network, we use a generalized version of the configuration model in which each node is a member of a specified number of triangles. SIR epidemics on this type of graph have earlier been investigated under the assumption of homogeneous infectivity and also under the assumption of Poisson transmission and recovery rates. We extend known results from literature by relaxing the assumption of homogeneous infectivity. An important special case of the epidemic model analyzed in this paper is epidemics in continuous time with arbitrary infectious period distribution. We use branching process approximations of the spread of the disease to provide expressions for the basic reproduction number R0, the probability of a major outbreak and the expected final size. In addition, the impact of random vaccination with a perfect vaccine on the final outcome of the epidemic is investigated. We find that, for this particular model, R0 equals the perfect vaccine-associated reproduction number. Generalizations to groups larger than three are discussed briefly.


title: 'Getting out of the closet': Scientific authorship of literary fiction and knowledge transfer

abstract: Some scientists write literary fiction books in their spare time. If these books contain scientific knowledge, literary fiction becomes a mechanism of knowledge transfer. In this case, we could conceptualize literary fiction as non-formal knowledge transfer. We model knowledge transfer via literary fiction as a function of the type of scientist (academic or non-academic) and his/her scientific field. Academic scientists are those employed in academia and public research organizations whereas non-academic scientists are those with a scientific background employed in other sectors. We also distinguish between direct knowledge transfer (the book includes the scientist's research topics), indirect knowledge transfer (scientific authors talk about their research with cultural agents) and reverse knowledge transfer (cultural agents give scientists ideas for future research). Through mixed-methods research and a sample from Spain, we find that scientific authorship accounts for a considerable percentage of all literary fiction authorship. Academic scientists do not transfer knowledge directly so often as non-academic scientists, but the former engage into indirect and reverse transfer knowledge more often than the latter. Scientists from History stand out in direct knowledge transfer. We draw propositions about the role of the academic logic and scientific field on knowledge transfer via literary fiction. We advance some tentative conclusions regarding the consideration of scientific authorship of literary fiction as a valuable knowledge transfer mechanism.


title: Nuclear spin-lattice relaxation in p-type GaAs

abstract: Spin-lattice relaxation of the nuclear spin system in p-type GaAs is studied using a three-stage experimental protocol including optical pumping and measuring the difference of the nuclear spin polarization before and after a dark interval of variable length. This method allows us to measure the spin-lattice relaxation time $T_1$ of optically pumped nuclei "in the dark", that is, in the absence of illumination. The measured $T_1$ values fall into the sub-second time range, being three orders of magnitude shorter than in earlier studied n-type GaAs. The drastic difference is further emphasized by magnetic-field and temperature dependences of $T_1$ in p-GaAs, showing no similarity to those in n-GaAs. This unexpected behavior is explained within a developed theoretical model involving quadrupole relaxation of nuclear spins, which is induced by electric fields within closely spaced donor-acceptor pairs.


title: Distributional Term Set Expansion

abstract: This paper is a short empirical study of the performance of centrality and classification based iterative term set expansion methods for distributional semantic models. Iterative term set expansion is an interactive process using distributional semantics models where a user labels terms as belonging to some sought after term set, and a system uses this labeling to supply the user with new, candidate, terms to label, trying to maximize the number of positive examples found. While centrality based methods have a long history in term set expansion, we compare them to classification methods based on the the Simple Margin method, an Active Learning approach to classification using Support Vector Machines. Examining the performance of various centrality and classification based methods for a variety of distributional models over five different term sets, we can show that active learning based methods consistently outperform centrality based methods.


title: Parameter estimation for discretely-observed linear birth-and-death processes

abstract: Birth-and-death processes are widely used to model the development of biological populations. Although they are relatively simple models, their parameters can be challenging to estimate, because the likelihood can become numerically unstable when data arise from the most common sampling schemes, such as annual population censuses. Simple estimators may be based on an embedded Galton-Watson process, but this presupposes that the observation times are equi-spaced. We estimate the birth, death, and growth rates of a linear birth-and-death process whose population size is periodically observed via an embedded Galton-Watson process, and by maximizing a saddlepoint approximation to the likelihood. We show that a Gaussian approximation to the saddlepoint-based likelihood connects the two approaches, we establish consistency and asymptotic normality of quasi-likelihood estimators, compare our estimators on some numerical examples, and apply our results to census data for two endangered bird populations and the H1N1 influenza pandemic.


title: Multilevel nested simulation for efficient risk estimation

abstract: We investigate the problem of computing a nested expectation of the form $\mathbb{P}[\mathbb{E}[X|Y] \!\geq\!0]\!=\!\mathbb{E}[\textrm{H}(\mathbb{E}[X|Y])]$ where $\textrm{H}$ is the Heaviside function. This nested expectation appears, for example, when estimating the probability of a large loss from a financial portfolio. We present a method that combines the idea of using Multilevel Monte Carlo (MLMC) for nested expectations with the idea of adaptively selecting the number of samples in the approximation of the inner expectation, as proposed by (Broadie et al., 2011). We propose and analyse an algorithm that adaptively selects the number of inner samples on each MLMC level and prove that the resulting MLMC method with adaptive sampling has an $\mathcal{O}\left( \varepsilon^{-2}|\log\varepsilon|^2 \right)$ complexity to achieve a root mean-squared error $\varepsilon$. The theoretical analysis is verified by numerical experiments on a simple model problem. We also present a stochastic root-finding algorithm that, combined with our adaptive methods, can be used to compute other risk measures such as Value-at-Risk (VaR) and Conditional Value-at-Risk (CVaR), with the latter being achieved with $\mathcal{O}\left(\varepsilon^{-2}\right)$ complexity.


title: Insisting on the role of experimental data: the pseudoscalar-pole piece to the $(g_\mu-2)$ and the $|V_{ub}|$ from $B \to \pi \ell \nu_{\ell}$ differential branching ratio

abstract: We employ a mathematical framework based on rational approximants in order to calculate meson form factors. The method profits from unitary, is systematic and data based, and is able to ascribe a systematic uncertainty which provides for the desired model independence. Two examples are discussed: the transition form factor entering the pseudoscalar-pole piece of the hadronic light-by-light contribution to the anomalous magnetic moment of the muon, and the $B \to \pi$ form factor participating the $B\to\pi\ell\nu_{\ell}$ differential branching ratios which allows to determine the $|V_{ub}|$ CKM parameter.


title: Formation of recurring slope lineae on Mars by rarefied gas-triggered granular flows

abstract: Active dark flows known as recurring slope lineae have been observed on the warmest slopes of equatorial Mars. The morphology, composition and seasonality of the lineae suggest a role of liquid water in their formation. However, internal and atmospheric sources of water appear to be insufficient to sustain the observed slope activity. Experimental evidence suggests that under the low atmospheric pressure at the surface of Mars, gas can flow upwards through porous Martian soil due to thermal creep under surface regions heated by the Sun, and disturb small particles. Here we present numerical simulations to demonstrate that such a dry process involving the pumping of rarefied gas in the Martian soil due to temperature contrasts can explain the formation of the recurring slope lineae. In our simulations, solar irradiation followed by shadow significantly reduces the angle of repose due to the resulting temporary temperature gradients over shaded terrain, and leads to flow at intermediate slope angles. The simulated flow locations are consistent with observed recurring slope lineae that initiate in rough and bouldered terrains with local shadows over the soil. We suggest that this dry avalanche process can explain the formation of the recurring slope lineae on Mars without requiring liquid water or CO2 frost activity.


title: Discovery of a strain-stabilised smectic electronic order in LiFeAs

abstract: In many high temperature superconductors, small orthorhombic distortions of the lattice structure result in surprisingly large symmetry breaking of the electronic states and macroscopic properties, an effect often referred to as nematicity. To directly study the impact of symmetry-breaking lattice distortions on the electronic states, using low-temperature scanning tunnelling microscopy we image at the atomic scale the influence of strain-tuned lattice distortions on the correlated electronic states in the iron-based superconductor LiFeAs, a material which in its ground state is tetragonal, with four-fold ($C_4$) symmetry. Our experiments uncover a new strain-stabilised modulated phase which exhibits a smectic order in LiFeAs, an electronic state which not only breaks rotational symmetry but also reduces translational symmetry. We follow the evolution of the superconducting gap from the unstrained material with $C_4$ symmetry through the new nematic phase with two-fold ($C_2$) symmetry and charge-density-wave order to a state where superconductivity is completely suppressed.


title: The $\varepsilon$-form of the differential equations for Feynman integrals in the elliptic case

abstract: Feynman integrals are easily solved if their system of differential equations is in $\varepsilon$-form. In this letter we show by the explicit example of the kite integral family that an $\varepsilon$-form can even be achieved, if the Feynman integrals do not evaluate to multiple polylogarithms. The $\varepsilon$-form is obtained by a (non-algebraic) change of basis for the master integrals.


title: Rayleigh-Taylor turbulence with singular nonuniform initial conditions

abstract: We perform direct numerical simulations of three dimensional Rayleigh-Taylor turbulence with a nonuniform singular initial temperature background. In such conditions, the mixing layer evolves under the driving of a varying effective Atwood number; the long-time growth is still self-similar, but not anymore proportional to $t^2$ and depends on the singularity exponent $c$ of the initial profile $\Delta T \propto z^c$. We show that the universality is recovered when looking at the efficiency, defined as the ratio of the variation rates of the kinetic energy over the heat flux. A closure model is proposed that is able to reproduce analytically the time evolution of the mean temperature profiles, in excellent agreement with the numerical results. Finally, we reinterpret our findings in the light of spontaneous stochasticity where the growth of the mixing layer is mapped into the propagation of a wave of turbulent fluctuations on a rough background.


title: PyFml - a Textual Language For Feature Modeling

abstract: The Feature model is a typical approach to capture variability in a software product line design and implementation. For that, most works automate feature model using a limited graphical notation represented by propositional logic and implemented by Prolog or Java programming languages. These works do not properly combine the extensions of classical feature models and do not provide scalability to implement large size problem issues. In this work, we propose a textual feature modeling language based on Python programming language (PyFML), that generalizes the classical feature models with instance feature cardinalities and attributes which be extended with highlight of replication and complex logical and mathematical cross-tree constraints. textX Meta-language is used for building PyFML to describe and organize feature model dependencies, and PyConstraint Problem Solver is used to implement feature model variability and its constraints validation. The work provides a textual human-readable language to represent feature model and maps the feature model descriptions directly into the object-oriented representation to be used by Constraint Problem Solver for computation. Furthermore, the proposed PyFML makes the notation of feature modeling more expressive to deal with complex software product line representations and using PyConstraint Problem Solver


title: Recursive Chaining of Reversible Image-to-image Translators For Face Aging

abstract: This paper addresses the modeling and simulation of progressive changes over time, such as human face aging. By treating the age phases as a sequence of image domains, we construct a chain of transformers that map images from one age domain to the next. Leveraging recent adversarial image translation methods, our approach requires no training samples of the same individual at different ages. Here, the model must be flexible enough to translate a child face to a young adult, and all the way through the adulthood to old age. We find that some transformers in the chain can be recursively applied on their own output to cover multiple phases, compressing the chain. The structure of the chain also unearths information about the underlying physical process. We demonstrate the performance of our method with precise and intuitive metrics, and visually match with the face aging state-of-the-art.


title: Totally non congruence Veech groups

abstract: Veech groups are discrete subgroups of SL(2, R) which play an important role in the theory of translation surfaces. For a special class of translation surfaces called origamis or square-tiled surfaces their Veech groups are subgroups of finite index of SL(2, Z). We show that each stratum of the space of translation surfaces contains infinitely many origamis whose Veech group is a totally non congruence group, i.e. it surjects to SL(2, Z/nZ) for any n.


title: Parity-Time Symmetry meets Photonics: A New Twist in non-Hermitian Optics

abstract: In the past decade, the concept of parity-time ($\mathcal{PT}$) symmetry, originally introduced in non-Hermitian extensions of quantum mechanical theories, has come into thinking of photonics, providing a fertile ground for studying, observing, and utilizing some of the peculiar aspects of $\mathcal{PT}$ symmetry in optics. Together with related concepts of non-Hermitian physics of open quantum systems, such as non-Hermitian degeneracies (exceptional points) and spectral singularities, $\mathcal{PT}$ symmetry represents one among the most fruitful ideas introduced in optics in the past few years. Judicious tailoring of optical gain and loss in integrated photonic structures has emerged as a new paradigm in shaping the flow of light in unprecedented ways, with major applications encompassing laser science and technology, optical sensing, and optical material engineering. In this perspective, I review some of the main achievements and emerging areas of $\mathcal{PT}$-symmetric and non-Hermtian photonics, and provide an outline of challenges and directions for future research in one of the fastest growing research area of photonics.


title: Directed cycles have the edge-Erd\H os-P\'osa property

abstract: In this short note we prove that for every $k\in \mathbb{N}$ there is a $t_k\in\mathbb{N}$ such that for every digraph $G$ there are either $k$ edge-disjoint directed cycles in $G$ or a set $X$ of at most $t_k$ edges such that $G-X$ contains no directed cycle.


title: Not to Cry Wolf: Distantly Supervised Multitask Learning in Critical Care

abstract: Patients in the intensive care unit (ICU) require constant and close supervision. To assist clinical staff in this task, hospitals use monitoring systems that trigger audiovisual alarms if their algorithms indicate that a patient's condition may be worsening. However, current monitoring systems are extremely sensitive to movement artefacts and technical errors. As a result, they typically trigger hundreds to thousands of false alarms per patient per day - drowning the important alarms in noise and adding to the exhaustion of clinical staff. In this setting, data is abundantly available, but obtaining trustworthy annotations by experts is laborious and expensive. We frame the problem of false alarm reduction from multivariate time series as a machine-learning task and address it with a novel multitask network architecture that utilises distant supervision through multiple related auxiliary tasks in order to reduce the number of expensive labels required for training. We show that our approach leads to significant improvements over several state-of-the-art baselines on real-world ICU data and provide new insights on the importance of task selection and architectural choices in distantly supervised multitask learning.


title: Inversions of synthetic umbral flashes: effects of the scanning time on the inferred atmospheres

abstract: The use of instruments that record narrow band images at selected wavelengths is a common approach in solar observations. They allow the scanning of a spectral line by sampling the Stokes profiles with 2D images at each line position, but require a compromise between spectral resolution and temporal cadence. We evaluate the impact of the time-dependent acquisition of different wavelengths on the inversion of spectropolarimetric profiles from chromospheric lines during umbral flashes. Simulations of non-linear wave propagation in a sunspot were performed with the code MANCHA. Synthetic Stokes parameters in the Ca II 8542 A line in NLTE were computed for an umbral flash using the code NICOLE. Artificial profiles with the same wavelength coverage and temporal cadence from reported observations were constructed and inverted. The inferred atmospheric stratifications were compared with the original models. The inferred atmospheres provide a reasonable characterization of the thermodynamic properties of the atmosphere during most of the phases of the umbral flash. Only at the early stages of the flash, when the shock wave reaches the formation height of the line, the Stokes profiles present apparent wavelength shifts and other spurious deformations. These features are misinterpreted by the inversion code, which can return unrealistic atmospheric models from a good fit of the Stokes profiles. The misguided results include flashed atmospheres with strong downflows, even though the simulation exhibits upflows during the umbral flash, and large variations in the magnetic field strength. Our analyses validate the inversion of Stokes profiles acquired by sequentially scanning certain selected wavelengths of a line profile, even in the case of rapidly-changing events such as umbral flashes. However, the inversions are unreliable during a short period at the development phase of the flash.


title: A co-located partitions strategy for parallel CFD-DEM couplings

abstract: In this work, a new partition-collocation strategy for the parallel execution of CFD--DEM couplings is investigated. Having a good parallel performance is a key issue for an Eulerian-Lagrangian software that aims to be applied to solve industrially significant problems, as the computational cost of these couplings is one of their main drawback. The approach presented here consists in co-locating the overlapping parts of the simulation domain of each software on the same MPI process, in order to reduce the cost of the data exchanges. It is shown how this strategy allows reducing memory consumption and inter-process communication between CFD and DEM to a minimum and therefore to overcome an important parallelization bottleneck identified in the literature. Three benchmarks are proposed to assess the consistency and scalability of this approach. A coupled execution on 280 cores shows that less than 0.1% of the time is used to perform inter-physics data exchange.


title: Facebook Use of Sensitive Data for Advertising in Europe

abstract: The upcoming European General Data Protection Regulation (GDPR) prohibits the processing and exploitation of some categories of personal data (health, political orientation, sexual preferences, religious beliefs, ethnic origin, etc.) due to the obvious privacy risks that may be derived from a malicious use of such type of information. These categories are referred to as sensitive personal data. Facebook has been recently fined EUR 1.2M in Spain for collecting, storing and processing sensitive personal data for advertising purposes. This paper quantifies the portion of Facebook users in the European Union (EU) who are labeled with interests linked to sensitive personal data. The results of our study reveal that Facebook labels 73% EU users with sensitive interests. This corresponds to 40% of the overall EU population. We also estimate that a malicious third-party could unveil the identity of Facebook users that have been assigned a sensitive interest at a cost as low as EUR 0.015 per user. Finally, we propose and implement a web browser extension to inform Facebook users of the sensitive interests Facebook has assigned them.


title: Tackling Multilabel Imbalance through Label Decoupling and Data Resampling Hybridization

abstract: The learning from imbalanced data is a deeply studied problem in standard classification and, in recent times, also in multilabel classification. A handful of multilabel resampling methods have been proposed in late years, aiming to balance the labels distribution. However these methods have to face a new obstacle, specific for multilabel data, as is the joint appearance of minority and majority labels in the same data patterns. We proposed recently a new algorithm designed to decouple imbalanced labels concurring in the same instance, called REMEDIAL (\textit{REsampling MultilabEl datasets by Decoupling highly ImbAlanced Labels}). The goal of this work is to propose a procedure to hybridize this method with some of the best resampling algorithms available in the literature, including random oversampling, heuristic undersampling and synthetic sample generation techniques. These hybrid methods are then empirically analyzed, determining how their behavior is influenced by the label decoupling process. As a result, a noteworthy set of guidelines on the combined use of these techniques can be drawn from the conducted experimentation.


title: The Interaction of Two Dyons in The Near Field Limit

abstract: We study the interaction of two dyons in the region of their cores where they are non-linear and non-Abelian. We assume the superposition of two dyons as a solution of the equation of motion. The terms due to the non-linearity of the strength tensor are considered as the perturbation terms which deforms the profile function of two individual dyons. As a result, the profile function of dyons are obtained to be dependent on the polar angle and the spherical symmetry is lost.


title: Dealing with Difficult Minority Labels in Imbalanced Mutilabel Data Sets

abstract: Multilabel classification is an emergent data mining task with a broad range of real world applications. Learning from imbalanced multilabel data is being deeply studied latterly, and several resampling methods have been proposed in the literature. The unequal label distribution in most multilabel datasets, with disparate imbalance levels, could be a handicap while learning new classifiers. In addition, this characteristic challenges many of the existent preprocessing algorithms. Furthermore, the concurrence between imbalanced labels can make harder the learning from certain labels. These are what we call \textit{difficult} labels. In this work, the problem of difficult labels is deeply analyzed, its influence in multilabel classifiers is studied, and a novel way to solve this problem is proposed. Specific metrics to assess this trait in multilabel datasets, called \textit{SCUMBLE} (\textit{Score of ConcUrrence among iMBalanced LabEls}) and \textit{SCUMBLELbl}, are presented along with REMEDIAL (\textit{REsampling MultilabEl datasets by Decoupling highly ImbAlanced Labels}), a new algorithm aimed to relax label concurrence. How to deal with this problem using the R mldr package is also outlined.


title: Cometary impactors on the TRAPPIST-1 planets can destroy all planetary atmospheres and rebuild secondary atmospheres on planets f, g, h

abstract: The TRAPPIST-1 system is unique in that it has a chain of seven terrestrial Earth-like planets located close to or in its habitable zone. In this paper, we study the effect of potential cometary impacts on the TRAPPIST-1 planets and how they would affect the primordial atmospheres of these planets. We consider both atmospheric mass loss and volatile delivery with a view to assessing whether any sort of life has a chance to develop. We ran N-body simulations to investigate the orbital evolution of potential impacting comets, to determine which planets are more likely to be impacted and the distributions of impact velocities. We consider three scenarios that could potentially throw comets into the inner region (i.e within 0.1au where the seven planets are located) from an (as yet undetected) outer belt similar to the Kuiper belt or an Oort cloud: Planet scattering, the Kozai-Lidov mechanism and Galactic tides. For the different scenarios, we quantify, for each planet, how much atmospheric mass is lost and what mass of volatiles can be delivered over the age of the system depending on the mass scattered out of the outer belt. We find that the resulting high velocity impacts can easily destroy the primordial atmospheres of all seven planets, even if the mass scattered from the outer belt is as low as that of the Kuiper belt. However, we find that the atmospheres of the outermost planets f, g and h can also easily be replenished with cometary volatiles (e.g. $\sim$ an Earth ocean mass of water could be delivered). These scenarios would thus imply that the atmospheres of these outermost planets could be more massive than those of the innermost planets, and have volatiles-enriched composition.


title: Nonnegative PARAFAC2: a flexible coupling approach

abstract: Modeling variability in tensor decomposition methods is one of the challenges of source separation. One possible solution to account for variations from one data set to another, jointly analysed, is to resort to the PARAFAC2 model. However, so far imposing constraints on the mode with variability has not been possible. In the following manuscript, a relaxation of the PARAFAC2 model is introduced, that allows for imposing nonnegativity constraints on the varying mode. An algorithm to compute the proposed flexible PARAFAC2 model is derived, and its performance is studied on both synthetic and chemometrics data.


title: Robust Continuous Co-Clustering

abstract: Clustering consists of grouping together samples giving their similar properties. The problem of modeling simultaneously groups of samples and features is known as Co-Clustering. This paper introduces ROCCO - a Robust Continuous Co-Clustering algorithm. ROCCO is a scalable, hyperparameter-free, easy and ready to use algorithm to address Co-Clustering problems in practice over massive cross-domain datasets. It operates by learning a graph-based two-sided representation of the input matrix. The underlying proposed optimization problem is non-convex, which assures a flexible pool of solutions. Moreover, we prove that it can be solved with a near linear time complexity on the input size. An exhaustive large-scale experimental testbed conducted with both synthetic and real-world datasets demonstrates ROCCO's properties in practice: (i) State-of-the-art performance in cross-domain real-world problems including Biomedicine and Text Mining; (ii) very low sensitivity to hyperparameter settings; (iii) robustness to noise and (iv) a linear empirical scalability in practice. These results highlight ROCCO as a powerful general-purpose co-clustering algorithm for cross-domain practitioners, regardless of their technical background.


title: Semiprojectivity and semiinjectivity in different categories

abstract: Projectivity and injectivity are fundamental notions in category theory. We consider natural weakenings termed semiprojectivity and semiinjectivity, and study these concepts in different categories. For example, in the category of metric spaces, (semi)injective objects are precisely the absolute (neighborhood) retracts. We show that the trivial group is the only semiinjective group, while every free product of a finitely presented group and a free group is semiprojective. To a compact, metric space $X$ we associate the commutative C*-algebra $C(X)$. This association is contravariant, whence semiinjectivity of $X$ is related to semiprojectivity of $C(X)$. Together with Adam S{\o}rensen, we showed that $C(X)$ is semiprojective in the category of all C*-algebras if and only if $X$ is an absolute neighborhood retract with dimension at most one.


title: Motion of interfaces for a damped hyperbolic Allen-Cahn equation

abstract: Consider the Allen-Cahn equation $u_t=\varepsilon^2\Delta u-F'(u)$, where $F$ is a double well potential with wells of equal depth, located at $\pm1$. There are a lot of papers devoted to the study of the limiting behavior of the solutions as the diffusion coefficient $\varepsilon\to0^+$, and it is well known that, if the initial datum $u(\cdot,0)$ takes the values $+1$ and $-1$ in the regions $\Omega_+$ and $\Omega_-$, then the "interface" connecting $\Omega_+$ and $\Omega_-$ moves with normal velocity equal to the sum of its principal curvatures, i.e. the interface moves by mean curvature flow. This paper concerns with the motion of the inteface for a damped hyperbolic Allen-Cahn equation, in a bounded domain of $\mathbb{R}^n$, for $n=2$ or $n=3$. In particular, we focus the attention on radially simmetric solutions, studying in detail the differences with the classic parabolic case, and we prove that, under appropriate assumptions on the initial data $u(\cdot,0)$ and $u_t(\cdot,0)$, the interface moves by mean curvature as $\varepsilon\to0^+$ also in the hyperbolic framework.


title: Super-blockers and the effect of network structure on information cascades

abstract: Modelling information cascades over online social networks is important in fields from marketing to civil unrest prediction, however the underlying network structure strongly affects the probability and nature of such cascades. Even with simple cascade dynamics the probability of large cascades are almost entirely dictated by network properties, with well-known networks such as Erdos-Renyi and Barabasi-Albert producing wildly different cascades from the same model. Indeed, the notion of 'superspreaders' has arisen to describe highly influential nodes promoting global cascades in a social network. Here we use a simple model of global cascades to show that the presence of locality in the network increases the probability of a global cascade due to the increased vulnerability of connecting nodes. Rather than 'super-spreaders', we find that the presence of these highly connected 'super-blockers' in heavy-tailed networks in fact reduces the probability of global cascades, while promoting information spread when targeted as the initial spreader.


title: Virtual depth by active background suppression: Revisiting the cosmic muon induced background of GERDA Phase II

abstract: In-situ production of long-lived isotopes by cosmic muon interactions may generate a non-negligible background for deep underground rare event searches. Previous Monte Carlo studies for the GERDA experiment at LNGS identified the delayed decays of $^{77}$Ge and its metastable state $^{77m}$Ge as dominant cosmogenic background in the search for neutrinoless double beta decay of $^{76}$Ge. This might limit the sensitivity of next generation experiments aiming for increased $^{76}$Ge mass at background-free conditions and thereby define a minimum depth requirement. A re-evaluation of the $^{77(m)}$Ge background for the GERDA experiment has been carried out by a set of Monte Carlo simulations. The obtained $^{77(m)}$Ge production rate is (0.21$\pm$0.01) nuclei/(kg$\cdot$yr). After application of state-of-the-art active background suppression techniques and simple delayed coincidence cuts this corresponds to a background contribution of (2.7$\pm$0.3)$\cdot10^{-6}$ cts/(keV$\cdot$kg$\cdot$yr). The suppression achieved by this strategy equals an effective muon flux reduction of more than one order of magnitude. This virtual depth increase opens the way for next generation rare event searches.


title: An extended Lagrangian formalism

abstract: A simple formal procedure makes the main properties of the lagrangian binomial extendable to functions depending to any kind of order of the time--derivatives of the lagrangian coordinates. Such a broadly formulated binomial can provide the lagrangian components, in the classical sense of the Newton's law, for a quite general class of forces. At the same time, the generalized equations of motions recover some of the classical alternative formulations of the Lagrangian equations.


title: Wall effects on spatial correlations of non-affine strain in a 3D model glass

abstract: Effects of hard planar walls with a particle scale roughness on the spatial correlations of non-affine strain in amorphous solids are investigated via molecular dynamics simulations. When determined within layers parallel to the wall plane, normalized non-affine strain correlations are enhanced within layers closer to the wall. The amplitude of these correlations, on the other hand, is found to be suppressed by the wall. While the former is connected to the effects of a hard boundary on the continuum mechanics scale, the latter is attributed to molecular scale wall effects on the size of the region (nearest neighbor cage), explored by particles on intermediate times scales.


title: Projection methods based on spline quasi-interpolation for Urysohn integral equations

abstract: In this paper we propose projection methods based on spline quasi-interpolating projectors of degree $d$ and class $C^{d-1}$ on a bounded interval for the numerical solution of nonlinear integral equations. We prove that they have high order of convergence $2d+2$ if $d$ is odd and $2d+3$ if $d$ is even. We also present the implementation details of the above methods. Finally, we provide numerical tests, that confirm the theoretical results. Moreover, we compare the theoretical and numerical results with those obtained by using a collocation method based on the same spline quasi-interpolating projectors.