title: ALMA detection of a tentative nearly edge-on rotating disk around the nearby AGB star R~Doradus

abstract: A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature, in combination with previous SPHERE data and theoretical work point towards the presence of a compact differentially rotating disk, orientated nearly edge-on. We model the $^{\rm 28}$SiO ($v=1,~J=8\to7$) emission with a disk model. We estimate the disk mass and angular momentum to be $3 \times 10^{-6}$ Solar masses and $5 \times 10^{40}\ {\rm m^2 kg/s}$. The latter presents an `angular momentum problem' that may be solved by assuming that the disk is the result of wind-companion interactions with a companion of at least 2.5 earth masses, located at 6 AU, the tentatively determined location of the disk's inner rim. An isolated clump of emission is detected to the south-east with a velocity that is high compared to the previously determined terminal velocity of the wind. Its position and mean velocity suggest that it may be associated with a companion planet, located at the disk's inner rim.


title: 3$\times$3 transfer matrix modelling of integrated optical devices

abstract: Unlike common devices based on ring resonators, the structure in Fig. 1.a involves not only 2$\times$2 couplers but also a 3$\times$3 coupler, which means that a 3$\times$3 transfer matrix approach is required to model the system. To the best of our knowledge, no such a model has been developed before. The only model available in the literature is based on a clever recursive 2$\times$2 transfer matrix model [1], which requires lengthy calculations that depend on the chosen boundary conditions and on the particular geometry chosen. The scope of this document is instead to show how to generalize the standard 2$\times$2 transfer matrix approach to cover any system with 3$\times$3 couplers, and calculate the transfer matrix of any complicated system just as a product of simple 3$\times$3 matrices.


title: Highly sensitive NO2 sensors by pulsed laser deposition on graphene

abstract: Graphene as a single-atomic-layer material is fully exposed to environment and has therefore a great potential for creating of sensitive gas sensors. However, in order to realize this potential for different polluting gases, graphene has to be functionalized - adsorption centers of different type and with high affinity to target gases have to be created at its surface. In this present work, modification of graphene by small amounts of laser ablated materials is introduced for this purpose as a versatile and precise tool. The approach was demonstrated with two very different materials chosen for pulsed laser deposition (PLD), a metal (Ag) and a dielectric oxide (ZrO2). It was shown that the gas response and its recovery rate can be significantly enhanced by choosing the PLD target material and deposition conditions. The response to NO2 gas in air was amplified up to 40 times in case of PLD-modified graphene in comparison with pristine graphene and reached 7-8% at 40 ppb of NO2 and 20-30% at 1 ppm of N2. These results were obtained after PLD in gas environment (5 x 10-2 mbar oxygen or nitrogen) and atomic areal densities of deposited materials of were about 10 15 cm-2. The ultimate level of NO2 detection in air, as extrapolated from the experimental data obtained at room temperature under mild UV-excitation, was below 1 ppb.


title: On the dimension datum of a subgroup. II

abstract: This paper studies three aspects around dimension datum: (1), a generalization of the dimension datum, which we call the tau-dimension datum; (2), dimension data of disconnected subgroups; (3), compactness of isospectral sets of normal homogeneous spaces.


title: A Numerical Model for the Construction of Finite Blaschke Products with Preassigned Distinct Critical Points

abstract: We present a numerical model for determining a finite Blaschke product of degree $n+1$ having $n$ preassigned distinct critical points $z_1,\dots,z_n$ in the complex (open) unit disk $\mathbb{D}$. The Blaschke product is uniquely determined up to postcomposition with conformal automorphisms of $\mathbb{D}$. The proposed method is based on the construction of a sparse nonlinear system where the data dependency is isolated to two vectors and on a certain transformation of the critical points. The efficiency and accuracy of the method is illustrated in several examples.


title: Hurwitz Stacks of Groups Extensions and Irreducibility

abstract: We study the irreducible components of special loci of curves whose group of symmetries is given as certain group extension. We introduce some relative Hurwitz data, which we show by using mixed \'etale cohomology theory, identifies some irreducible components for rational and normal non-abelian special loci and Hurwitz spaces. A heuristic, that is supported by three classes of examples, provides an additional context for building further irreducible loci.


title: Smartphone based Driving Style Classification Using Features Made by Discrete Wavelet Transform

abstract: Smartphones consist of different sensors, which provide a platform for data acquisition in many scientific researches such as driving style identification systems. In the present paper, smartphone data are used to evaluate the driving styles based on maneuvers analysis. The data obtained for each maneuver is the speed of the vehicle steering and the vehicle's direct and lateral acceleration. To classify the drivers based on their driving style, machine-learning algorithms can be used on these data. However, these data usually contains more information than it is needed and cause a bad effect on the learning accuracy. In addition, they may transfer some wrong information to the learning algorithm. Thus, we used Haar discrete wavelet transformation to remove noise effects. Then, we get the discrete wavelet transformation with four levels from smartphone sensors data, which include low-to-high frequencies, respectively. The obtained features vector for each maneuver includes the raw signal variance as well as the variance of the wavelet transform components. On these vectors, we use the k-nearest neighbors algorithm for features selection. Then, we use SVM, RBF and MLP neural networks on these features to separate braking and dangerous speed maneuvers from the safe ones as well as dangerous turning, U-turn and lane-changing maneuvers. The results are very interesting.


title: How sure are we? Two approaches to statistical inference

abstract: Suppose you are told that taking a statin will reduce your risk of a heart attack or stroke by 3% in the next ten years, or that women have better emotional intelligence than men. You may wonder how accurate the 3% is, or how confident we should be about the assertion about women's emotional intelligence, bearing in mind that these conclusions are only based on samples of data? My aim here is to present two statistical approaches to questions like these. Approach 1 is often called null hypothesis testing but I prefer the phrase "baseline hypothesis": this is the standard approach in many areas of inquiry but is fraught with problems. Approach 2 can be viewed as a generalisation of the idea of confidence intervals, or as the application of Bayes' theorem. Unlike Approach 1, Approach 2 provides a tentative estimate of the probability of hypotheses of interest. For both approaches, I explain, from first principles, building only on "common sense" statistical concepts like averages and randomness, both how to derive answers, and the rationale behind the answers. This is achieved by using computer simulation methods (resampling and bootstrapping using a spreadsheet available on the web) which avoid the use of probability distributions (t, normal, etc). Such a minimalist, but reasonably rigorous, analysis is particularly useful in a discipline like statistics which is widely used by people who are not specialists. My intended audience includes both statisticians, and users of statistical methods who are not statistical experts.


title: Inverse limits of Macaulay's inverse systems

abstract: Generalizing a result of Masuti and the second author, we describe inverse limits of Macaulay's inverse systems for Cohen-Macaulay factor algebras of formal power series or polynomial rings over an infinite field. On the way we find a strictness result for filtrations defined by regular sequences. It generalizes both a lemma of Uli Walther and the Rees isomorphism.


title: Approximating Dominating Set on Intersection Graphs of Rectangles and L-frames

abstract: We consider the Minimum Dominating Set (MDS) problem on the intersection graphs of geometric objects. Even for simple and widely-used geometric objects such as rectangles, no sub-logarithmic approximation is known for the problem and (perhaps surprisingly) the problem is NP-hard even when all the rectangles are "anchored" at a diagonal line with slope -1 (Pandit, CCCG 2017). In this paper, we first show that for any $\epsilon>0$, there exists a $(2+\epsilon)$-approximation algorithm for the MDS problem on "diagonal-anchored" rectangles, providing the first $O(1)$-approximation for the problem on a non-trivial subclass of rectangles. It is not hard to see that the MDS problem on "diagonal-anchored" rectangles is the same as the MDS problem on "diagonal-anchored" L-frames: the union of a vertical and a horizontal line segment that share an endpoint. As such, we also obtain a $(2+\epsilon)$-approximation for the problem with "diagonal-anchored" L-frames. On the other hand, we show that the problem is APX-hard in case the input L-frames intersect the diagonal, or the horizontal segments of the L-frames intersect a vertical line. However, as we show, the problem is linear-time solvable in case the L-frames intersect a vertical as well as a horizontal line. Finally, we consider the MDS problem in the so-called "edge intersection model" and obtain a number of results, answering two questions posed by Mehrabi (WAOA 2017).


title: Shellability of face posets of electrical networks and the CW poset property

abstract: We prove a conjecture of Thomas Lam that the face posets of stratified spaces of planar resistor networks are shellable. These posets are called uncrossing partial orders. This shellability result combines with Lam's previous result that these same posets are Eulerian to imply that they are CW posets, namely that they are face posets of regular CW complexes. Certain subsets of uncrossing partial orders are shown to be isomorphic to type A Bruhat order intervals; our shelling is shown to coincide on these intervals with a Bruhat order shelling which was constructed by Matthew Dyer using a reflection order. Our shelling for uncrossing posets also yields an explicit shelling for each interval in the face posets of the edge product spaces of phylogenetic trees, namely in the Tuffley posets, by virtue of each interval in a Tuffley poset being isomorphic to an interval in an uncrossing poset. This yields a more explicit proof of the result of Gill, Linusson, Moulton and Steel that the CW decomposition of Moulton and Steel for the edge product space of phylogenetic trees is a regular CW decomposition.


title: Maximal antipodal sets in irreducible compact symmetric spaces

abstract: We give an explicit classification of maximal antipodal sets in any irreducible compact symmetric space except for spin groups and half spin groups, and some quotient symmetric spaces associated to them.


title: Experimentally Generated Randomness Certified by the Impossibility of Superluminal Signals

abstract: From dice to modern complex circuits, there have been many attempts to build increasingly better devices to generate random numbers. Today, randomness is fundamental to security and cryptographic systems, as well as safeguarding privacy. A key challenge with random number generators is that it is hard to ensure that their outputs are unpredictable. For a random number generator based on a physical process, such as a noisy classical system or an elementary quantum measurement, a detailed model describing the underlying physics is required to assert unpredictability. Such a model must make a number of assumptions that may not be valid, thereby compromising the integrity of the device. However, it is possible to exploit the phenomenon of quantum nonlocality with a loophole-free Bell test to build a random number generator that can produce output that is unpredictable to any adversary limited only by general physical principles. With recent technological developments, it is now possible to carry out such a loophole-free Bell test. Here we present certified randomness obtained from a photonic Bell experiment and extract 1024 random bits uniform to within $10^{-12}$. These random bits could not have been predicted within any physical theory that prohibits superluminal signaling and allows one to make independent measurement choices. To certify and quantify the randomness, we describe a new protocol that is optimized for apparatuses characterized by a low per-trial violation of Bell inequalities. We thus enlisted an experimental result that fundamentally challenges the notion of determinism to build a system that can increase trust in random sources. In the future, random number generators based on loophole-free Bell tests may play a role in increasing the security and trust of our cryptographic systems and infrastructure.


title: Collective motion in prolate {\gamma}-rigid nuclei within minimal length concept via a quantum perturbation method

abstract: Based on the minimal length concept, inspired by Heisenberg algebra, a closed analytical formula is derived for the energy spectrum of the prolate {\gamma}-rigid Bohr-Mottelson Hamiltonian of nuclei, within a quantum perturbation method (QPM), by considering a scaled Davidson potential in \b{eta} shape variable. In the resulting solution, called X(3)-D-ML, the ground state and the first \b{eta}-band are all studied as a function of the free parameters. The fact of introducing the minimal length concept with a QPM makes the model very flexible and a powerful approach to describe nuclear collective excitations of a variety of vibrational-like nuclei. The introduction of scaling parameters in the Davidson potential enables us to get a physical minimum of this latter in comparison with previous works. The analysis of the corrected wave function, as well as the probability density distribution, shows that the minimal length parameter has a physical upper bound limit.


title: The Climate of the Khagan. Observations on palaeo-environmental Factors of the History of the Avars (6th-9th century AD)

abstract: Based on palaeoenvironmental, historical and archaeological data, the paper proposes possible climatic impacts on the history of the Avar Khaganate, which comprised the Carpathian Basin between the late 6th and the early 9th century AD. While the establishment of the Avars in East Central Europe took place within a period characterised by cold and dry climatic conditions (recently identified as Late Antique Little Ice Age), more stable climatic parameters may have favoured the stabilisation of Avar rule after a crisis in the aftermath of 626 AD. Data indicates growth of settlement and agricultural activity up to the mid-8th century. These developments did not necessarily strengthen central power, but may have contributed to a greater autonomy of various groups on the basis of increased resources. The Khaganate quickly disintegrated faced by the Carolingian advance of the 790s; the last decades of documented Avar presence were again accompanied by environmental vicissitudes.


title: Quasi-Optimality of an Adaptive Finite Element Method for Cathodic Protection

abstract: In this work, we derive a reliable and efficient residual-typed error estimator for the finite element approximation of a 2d cathodic protection problem governed by a steady-state diffusion equation with a nonlinear boundary condition. We propose a standard adaptive finite element method involving the D\"{o}rfler marking and a minimal refinement without the interior node property. Furthermore, we establish the contraction property of this adaptive algorithm in terms of the sum of the energy error and the scaled estimator. This essentially allows for a quasi-optimal convergence rate in terms of the number of elements over the underlying triangulation. Numerical experiments are provided to confirm this quasi-optimality.


title: Efficient construction of threshold networks of stock markets

abstract: Although the threshold network is one of the most used tools to characterize the underlying structure of a stock market, the identification of the optimal threshold to construct a reliable stock network remains challenging. In this paper, the concept of dynamic consistence between the threshold network and the stock market is proposed. The optimal threshold is estimated by maximizing the consistence function. The application of this procedure to stocks belonging to Standard \& Pool's 500 Index from January 2006 to December 2011 yields the threshold value 0.28. In analyzing topological characteristics of the generated network, three globally financial crises can be distinguished well from the evolutionary perspective.


title: Folding a 3D Euclidean space

abstract: This paper considers an extension of origami geometry to the case of "folding" a three dimensional (3D) space along a plane. First, all possible incidence constraints between given points, lines and planes are analyzed by using the geometry of reflections. Next, a set of 3D elementary fold operations is defined, which satisfy specific combinations of constraints with a finite number of solutions. The set consists of 47 valid fold operations, and solutions to some of them are explored to determine their number and conditions of existence.


title: Graphene functionalised by laser ablated V2O5 as highly sensitive NH3 sensor

abstract: Graphene has been recognized as a promising gas sensing material. The response of graphene-based sensors can be radically improved by introducing defects in graphene using, e. g., metal or metal oxide nanoparticles. We have functionalised CVD grown, single layer graphene by applying pulsed laser deposition (PLD) of V2O5 which resulted in a thin V2O5 layer on graphene with average thickness of ~0.6 nm. According to Raman analysis, PLD process also induced defects in graphene. Compared to unmodified graphene, the obtained chemiresistive sensor showed considerable improvement of sensing ammonia at room temperature. In addition, also the response time, sensitivity and reversibility were essentially enhanced due to graphene functionalisation by laser deposited V2O5. This can be explained by increased surface density of gas adsorption sites introduced by high energy atoms in laser ablation plasma and formation of nanophase boundaries between deposited V2O5 and graphene.


title: Glyph: Symbolic Regression Tools

abstract: We present Glyph - a Python package for genetic programming based symbolic regression. Glyph is designed for usage let by numerical simulations let by real world experiments. For experimentalists, glyph-remote provides a separation of tasks: a ZeroMQ interface splits the genetic programming optimization task from the evaluation of an experimental (or numerical) run. Glyph can be accessed at http://github.com/ambrosys/glyph . Domain experts are be able to employ symbolic regression in their experiments with ease, even if they are not expert programmers. The reuse potential is kept high by a generic interface design. Glyph is available on PyPI and Github.


title: On the standard $L$-function for $GSp_{2n} \times GL_1$ and algebraicity of symmetric fourth $L$-values for $GL_2$

abstract: We prove an explicit integral representation -- involving the pullback of a suitable Siegel Eisenstein series -- for the twisted standard $L$-function associated to a holomorphic vector-valued Siegel cusp form of degree $n$ and arbitrary level. In contrast to all previously proved pullback formulas in this situation, our formula involves only scalar-valued functions despite being applicable to $L$-functions of vector-valued Siegel cusp forms. The key new ingredient in our method is a novel choice of local vectors at the archimedean place which allows us to exactly compute the archimedean local integral. By specializing our integral representation to the case $n=2$ we are able to prove a reciprocity law -- predicted by Deligne's conjecture -- for the critical special values of the twisted standard $L$-function for vector-valued Siegel cusp forms of degree 2 and arbitrary level. This arithmetic application generalizes previously proved critical-value results for the full level case. By specializing further to the case of Siegel cusp forms obtained via the Ramakrishnan--Shahidi lift, we obtain a reciprocity law for the critical special values of the symmetric fourth $L$-function of a classical newform.


title: Six-vertex model and non-linear differential equations II. Continuous symmetries

abstract: This paper is a continuation of our previous work "Six-vertex model and non-linear differential equations I. Spectral problem" in which we have put forward a method for studying the spectrum of the six-vertex model based on non-linear differential equations. Here we intend to elaborate on that approach and also discuss properties of the spectrum unveiled by the aforementioned differential formulation of the transfer matrix's eigenvalue problem. In particular, we intend to demonstrate how this differential approach allows one to study continuous symmetries of the transfer matrix's spectrum through the Lie groups method.


title: Further Consequences of the Colorful Helly Hypothesis

abstract: Let $\mathcal{F}$ be a family of convex sets in ${\mathbb R}^d$, which are colored with $d+1$ colors. We say that $\mathcal{F}$ satisfies the Colorful Helly Property if every rainbow selection of $d+1$ sets, one set from each color class, has a non-empty common intersection. The Colorful Helly Theorem of Lov\'asz states that for any such colorful family $\mathcal{F}$ there is a color class $\mathcal{F}_i\subset \mathcal{F}$, for $1\leq i\leq d+1$, whose sets have a non-empty intersection. We establish further consequences of the Colorful Helly hypothesis. In particular, we show that for each dimension $d\geq 2$ there exist numbers $f(d)$ and $g(d)$ with the following property: either one can find an additional color class whose sets can be pierced by $f(d)$ points, or all the sets in $\mathcal{F}$ can be crossed by $g(d)$ lines.


title: Power Distribution System Synchrophasors with Non-Gaussian Errors: Real-World Measurements and Analysis

abstract: This letter studies the synchrophasor measurement error of electric power distribution systems with on-line and off-line measurements using graphical and numerical tests. It demonstrates that the synchrophasor measurement error follows a non-Gaussian distribution instead of the traditionally-assumed Gaussian distribution. It suggests the need to use non-Gaussian or Gaussian mixture models to represent the synchrophasor measurement error. These models are more realistic to accurately represent the error than the traditional Gaussian model. The measurements and underlying analysis will be helpful for the understanding of distribution system measurement characteristics, and also for the modeling and simulation of distribution system applications.


title: A Programmable CMOS Transceiver for Structural Health Monitoring

abstract: We describe a highly-integrated CMOS transceiver for active structural health monitoring (SHM). The chip actuates piezoelectric transducers and also senses ultrasound waves received by the same or another transducer. The transmitter uses an integer-N frequency synthesizer and pulse-width modulation (PWM) to generate low-distortion, band-limited waveforms up to 12.7 Vpp with center frequency from 0.1-2.75 MHz. The integrated offset-canceling fully-differential receiver has programmable gain and bandwidth, and uses quadrature demodulation to extract both amplitude and phase of the received waveforms for further signal processing. The transceiver was fabricated in a 0.5 um CMOS process and has been validated using (2D) damage localization on an SHM test bed.


title: Free Energies and the Reversed HLS Inequality

abstract: We prove reversed Hardy-Littlewood-Sobolev inequalities by carefully studying the natural associated free energies with direct methods of calculus of variations. Tightness is obtained by a dyadic argument, which quantifies the relative strength of the entropy functional versus the interaction energy. The existence of optimizers is shown in the class of $\prob$. With respect to their regularity, we study conditions for optimizers to be bounded functions. In a related model, we show the condensation phenomena, which suggests that optimizers are not in general regular.


title: Five-dimensional warped product space-time with time-dependent warping and a scalar field in the bulk

abstract: We have considered gravity in a five-dimensional warped product space-time, with a time-dependent warp factor and a time-dependent extra dimension. The five-dimensional field equations are derived for a spatially flat FRW brane and the energy conditions and the nature of bulk geometry have been examined. It is found that the expansion of the four-dimensional universe depends on its location along the extra dimension and will be different at different locations in the bulk spacetime. At low energies, the trapping of fields within the brane implies a specific correlation between the warp factor and the extra-dimensional scale factor. Generally, the bulk is not conformally flat. At high energies, the bulk is assumed to be sourced by a scalar field with self-interaction. The analysis shows that the potential of the scalar field source of gravity at a given position along the fifth dimension, is related to the Hubble parameter on the brane at that position in the bulk.


title: Approximate Method of Variational Bayesian Matrix Factorization/Completion with Sparse Prior

abstract: We derive analytical expression of matrix factorization/completion solution by variational Bayes method, under the assumption that observed matrix is originally the product of low-rank dense and sparse matrices with additive noise. We assume the prior of sparse matrix is Laplace distribution by taking matrix sparsity into consideration. Then we use several approximations for derivation of matrix factorization/completion solution. By our solution, we also numerically evaluate the performance of sparse matrix reconstruction in matrix factorization, and completion of missing matrix element in matrix completion.

categories: eess.SP cond-mat.dis-nn cs.IT cs.LG math.IT


title: Identifying Components from Object-Oriented APIs Based on Dynamic Analysis

abstract: The reuse at the component level is generally more effective than the one at the object-oriented class level. This is due to the granularity level where components expose their functionalities at an abstract level compared to the fine-grained object-oriented classes. Moreover, components clearly define their dependencies through their provided and required interfaces in an explicit way that facilitates the understanding of how to reuse these components. Therefore, several component identification approaches have been proposed to identify components based on the analysis object-oriented software applications. Nevertheless, most of the existing component identification approaches did not consider co-usage dependencies between API classes to identify classes/methods that can be reused to implement a specific scenario. In this paper, we propose an approach to identify reusable software components in object-oriented APIs, based on the interactions between client applications and the targeted API. As we are dealing with actual clients using the API, dynamic analysis allows to better capture the instances of API usage. Approaches using static analysis are usually limited by the difficulty of handling dynamic features such as polymorphism and class loading. We evaluate our approach by applying it to three Java APIs with eight client applications from the DaCapo benchmark. DaCapo provides a set of pre-defined usage scenarios. The results show that our component identification approach has a very high precision.


title: Chemi-net: a graph convolutional network for accurate drug property prediction

abstract: Absorption, distribution, metabolism, and excretion (ADME) studies are critical for drug discovery. Conventionally, these tasks, together with other chemical property predictions, rely on domain-specific feature descriptors, or fingerprints. Following the recent success of neural networks, we developed Chemi-Net, a completely data-driven, domain knowledge-free, deep learning method for ADME property prediction. To compare the relative performance of Chemi-Net with Cubist, one of the popular machine learning programs used by Amgen, a large-scale ADME property prediction study was performed on-site at Amgen. The results showed that our deep neural network method improved current methods by a large margin. We foresee that the significantly increased accuracy of ADME prediction seen with Chemi-Net over Cubist will greatly accelerate drug discovery.


title: Adiabatic absorbers in photonics simulations with the volume integral equation method

abstract: This paper describes the implementation and performance of adiabatic absorbing layers in an FFT-accelerated volume integral equation (VIE) method for simulating truncated nanophotonics structures. At the truncation sites, we place absorbing regions in which the conductivity is increased gradually in order to minimize reflections. In the continuous setting, such adiabatic absorbers have been shown via coupled-mode theory to produce reflections that diminish at a rate related to the smoothness of the absorption profile function. The VIE formulation we employ relies on uniform discretizations of the geometry over which the continuously varying fields and material properties are represented by piecewise constant functions. Such a discretization enables the acceleration of the method via the FFT and, furthermore, the introduction of varying absorption can be performed in a straightforward manner without compromising this speedup. We demonstrate that, in spite of the crude discrete approximation to the smooth absorption profiles, our approach recovers the theoretically predicted reflection behavior of adiabatic absorbers. We thereby show that the FFT-accelerated VIE method is an effective and fast simulation tool for nanophotonics simulations.


title: A Sharp Phase Field Method

abstract: Phase field modelling offers an extremely general framework to predict microstructural evolutions in complex systems. However, its computational implementation requires a discretisation scheme with a grid spacing small enough to preserve the continuous character of the theory. We present here a new formulation, which is intrinsically discrete, in which the interfaces are resolved with essentially one grid point with no pinning on the grid and an accurate rotational invariance, improving drastically the numerical capabilities of the method. We show that interfacial kinetic properties are reproduced with a high accuracy. Finally, we apply the model to a situation where conserved and non-conserved fields are coupled.


title: Trianguloids and Triangulations of Root Polytopes

abstract: Triangulations of a product of two simplices and, more generally, of root polytopes are closely related to Gelfand-Kapranov-Zelevinsky's theory of discriminants, to tropical geometry, tropical oriented matroids, and to generalized permutohedra. We introduce a new approach to these objects, identifying a triangulation of a root polytope with a certain bijection between lattice points of two generalized permutohedra. In order to study such bijections, we define trianguloids as edge-colored graphs satisfying simple local axioms. We prove that trianguloids are in bijection with triangulations of root polytopes.


title: Modelling the bow shock Pulsar Wind Nebulae propagating through a non-uniform ISM

abstract: Many pulsars propagate through the interstellar medium (ISM) with supersonic velocities, and their pulsar winds interact with the interstellar medium (ISM), forming bow shocks and magnetotails (PWN). We model the propagation of pulsars through the inhomogeneous ISM using non-relativistic axisymmetric magneto-hydrodynamic (MHD) simulations. We take into account the wind from the star, and the azimuthal and poloidal components of the magnetic field, and investigate the PWN at different levels of magnetization (the ratio of magnetic to matter energy-densities) in the wind. We consider the interaction of PWN with small-scale and large-scale imhomogeneities in the ISM at different values of magnetization. We conclude that the inhomogeneities in the ISM can change the shapes of the bow shocks and magnetotails at different values of the magnetization.We compare the results of our simulations with the images of the Guitar Nebula and other PWN that show irregularities in the shapes of their bow shocks and magnetotails. We conclude that these irregularities may be caused by the interaction of PWN with the inhomogeneities in the ISM.


title: "Shadowy" modes in Higher-Order Scalar-Tensor theories

abstract: We consider Higher-Order Scalar-Tensor theories which appear degenerate when restricted to the unitary gauge but are not degenerate in an arbitrary gauge. We dub them U-degenerate theories. We provide a full classification of theories that are either DHOST or U-degenerate and that are quadratic in second derivatives of the scalar field, and discuss its extension to cubic and higher order theories. Working with a simple example of U-degenerate theory, we find that, for configurations in which the scalar field gradient is time-like, the apparent extra mode in such a theory can be understood as a generalized instantaneous, or "shadowy" mode, which does not propagate. Appropriate boundary conditions, required by the elliptic nature of part of the equations of motion, lead to the elimination of the apparent instability associated with this extra mode.


title: Maxwell-Higgs vortices with internal structure

abstract: Vortices are considered in relativistic Maxwell-Higgs systems in interaction with a neutral scalar field. The gauge field interacts with the neutral field via the presence of generalized permeability, and the charged and neutral scalar fields interact in a way dictated by the presence of first order differential equations that solve the equations of motion. The neutral field may be seen as the source field of the vortex, and we study some possibilities, which modify the standard Maxwell-Higgs solution and include internal structure to the vortex.


title: Gradients on Sets

abstract: For a locally Lipschitz continuous function $f:X\to\mathbb{R}$ the generalized gradient $\partial f(x)$ of Clarke is used to develop some (set-valued) gradient on a set $A\subset X$. Existence, uniqueness and some approximation are considered for optimal descent directions on set $A$. The results serve as basis for nonsmooth numerical descent algorithms that can be found in subsequent papers.


title: Induced Saturation of Graphs

abstract: A graph $G$ is $H$-saturated for a graph $H$, if $G$ does not contain a copy of $H$ but adding any new edge to $G$ results in such a copy. An $H$-saturated graph on a given number of vertices always exists and the properties of such graphs, for example their highest density, have been studied intensively. A graph $G$ is $H$-induced-saturated if $G$ does not have an induced subgraph isomorphic to $H$, but adding an edge to $G$ from its complement or deleting an edge from $G$ results in an induced copy of $H$. It is not immediate anymore that $H$-induced-saturated graphs exist. In fact, Martin and Smith (2012) showed that there is no $P_4$-induced-saturated graph. Behrens et.al. (2016) proved that if $H$ belongs to a few simple classes of graphs such as a class of odd cycles of length at least $5$, stars of size at least $2$, or matchings of size at least $2$, then there is an $H$-induced-saturated graph. This paper addresses the existence question for $H$-induced-saturated graphs. It is shown that Cartesian products of cliques are $H$-induced-saturated graphs for $H$ in several infinite families, including large families of trees. A complete characterization of all connected graphs $H$ for which a Cartesian product of two cliques is an $H$-induced-saturated graph is given. Finally, several results on induced saturation for prime graphs and families of graphs are provided.


title: Resonances in the asteroid and trans-Neptunian belts: a brief review

abstract: Mean motion resonances play a fundamental role in the dynamics of the small bodies of the Solar System. The last decades of the 20th century gave us a detailed description of the dynamics as well as the process of capture of small bodies in coplanar or small inclination resonant orbits. More recently, semianalytical or numerical methods allowed us to explore the behavior of resonant motions for arbitrary inclination orbits. The emerging dynamics is very rich, including large orbital changes due to secular effects inside mean motion resonances. The process of capture in highly inclined or retrograde resonant orbits was addressed showing that the capture in retrograde resonances is more efficient than in direct ones. A new terminology appeared in order to characterize the properties of the resonances. Numerical explorations in the transneptunian region showed the relevance and the particular dynamics of the exterior resonances with Neptune which can account for some of the known high perihelion orbits in the scattered disk. Moreover, several asteroids evolving in resonance with planets other than Jupiter or Neptune were found and a large number of asteroids in three-body resonances were identified.


title: On the use of the KMR unintegrated parton distribution functions

abstract: We discuss the unintegrated parton distribution functions (UPDFs) introduced by Kimber, Martin and Ryskin (KMR), which are frequently used in phenomenological analyses of hard processes with transverse momenta of partons taken into account. We demonstrate numerically that the commonly used differential definition of the UPDFs leads to erroneous results for large transverse momenta. We identify the reason for that, being the use of the ordinary PDFs instead of the cutoff dependent distribution functions. We show that in phenomenological applications, the integral definition of the UPDFs with the ordinary PDFs can be used.


title: Coordinating users of shared facilities via data-driven predictive assistants and game theory

abstract: We study data-driven assistants that provide congestion forecasts to users of shared facilities (roads, cafeterias, etc.), to support coordination between them, and increase efficiency of such collective systems. Key questions are: (1) when and how much can (accurate) predictions help for coordination, and (2) which assistant algorithms reach optimal predictions? First we lay conceptual ground for this setting where user preferences are a priori unknown and predictions influence outcomes. Addressing (1), we establish conditions under which self-fulfilling prophecies, i.e., "perfect" (probabilistic) predictions of what will happen, solve the coordination problem in the game-theoretic sense of selecting a Bayesian Nash equilibrium (BNE). Next we prove that such prophecies exist even in large-scale settings where only aggregated statistics about users are available. This entails a new (nonatomic) BNE existence result. Addressing (2), we propose two assistant algorithms that sequentially learn from users' reactions, together with optimality/convergence guarantees. We validate one of them in a large real-world experiment.


title: 3D Video Quality Metric for Mobile Applications

abstract: In this paper, we propose a new full-reference quality metric for mobile 3D content. Our method is modeled around the Human Visual System, fusing the information of both left and right channels, considering color components, the cyclopean views of the two videos and disparity. Our method is assessing the quality of 3D videos displayed on a mobile 3DTV, taking into account the effect of resolution, distance from the viewers eyes, and dimensions of the mobile display. Performance evaluations showed that our mobile 3D quality metric monitors the degradation of quality caused by several representative types of distortion with 82 percent correlation with results of subjective tests, an accuracy much better than that of the state of the art mobile 3D quality metric.


title: Link prediction for interdisciplinary collaboration via co-authorship network

abstract: We analyse the Publication and Research (PURE) data set of University of Bristol collected between $2008$ and $2013$. Using the existing co-authorship network and academic information thereof, we propose a new link prediction methodology, with the specific aim of identifying potential interdisciplinary collaboration in a university-wide collaboration network.


title: On the nonlinear wave equation with time periodic potential

abstract: It is known that for some time periodic potentials $q(t, x) \geq 0$ having compact support with respect to $x$ some solutions of the Cauchy problem for the wave equation $\partial_t^2 u - \Delta_x u + q(t,x)u = 0$ have exponentially increasing energy as $t \to \infty$. We show that if one adds a nonlinear defocusing interaction $|u|^ru, 2\leq r < 4,$ then the solution of the nonlinear wave equation exists for all $t \in {\mathbb R}$ and its energy is polynomially bounded as $t \to \infty$ for every choice of $q$. Moreover, we prove that the zero solution of the nonlinear wave equation is instable if the corresponding linear equation has the property mentioned above.


title: Tropical integrable systems and Young tableaux: Shape equivalence and Littlewood-Richardson correspondence

abstract: We present a new characterization of the shape equivalent class and the Littlewood-Richardson correspondence of Young tableaux in terms of tropical (ultradiscrete) integrable systems. As an application, an alternative proof of the "shape change theorem" is given.


title: Joint Recognition of Handwritten Text and Named Entities with a Neural End-to-end Model

abstract: When extracting information from handwritten documents, text transcription and named entity recognition are usually faced as separate subsequent tasks. This has the disadvantage that errors in the first module affect heavily the performance of the second module. In this work we propose to do both tasks jointly, using a single neural network with a common architecture used for plain text recognition. Experimentally, the work has been tested on a collection of historical marriage records. Results of experiments are presented to show the effect on the performance for different configurations: different ways of encoding the information, doing or not transfer learning and processing at text line or multi-line region level. The results are comparable to state of the art reported in the ICDAR 2017 Information Extraction competition, even though the proposed technique does not use any dictionaries, language modeling or post processing.


title: Land cover mapping at very high resolution with rotation equivariant CNNs: towards small yet accurate models

abstract: In remote sensing images, the absolute orientation of objects is arbitrary. Depending on an object's orientation and on a sensor's flight path, objects of the same semantic class can be observed in different orientations in the same image. Equivariance to rotation, in this context understood as responding with a rotated semantic label map when subject to a rotation of the input image, is therefore a very desirable feature, in particular for high capacity models, such as Convolutional Neural Networks (CNNs). If rotation equivariance is encoded in the network, the model is confronted with a simpler task and does not need to learn specific (and redundant) weights to address rotated versions of the same object class. In this work we propose a CNN architecture called Rotation Equivariant Vector Field Network (RotEqNet) to encode rotation equivariance in the network itself. By using rotating convolutions as building blocks and passing only the the values corresponding to the maximally activating orientation throughout the network in the form of orientation encoding vector fields, RotEqNet treats rotated versions of the same object with the same filter bank and therefore achieves state-of-the-art performances even when using very small architectures trained from scratch. We test RotEqNet in two challenging sub-decimeter resolution semantic labeling problems, and show that we can perform better than a standard CNN while requiring one order of magnitude less parameters.


title: Superfluid liquid crystals: pasta phases in neutron star crusts

abstract: The pasta phases predicted to occur near the inner boundary of the crust of a neutron star resemble liquid crystals, a smectic A in the case of sheet-like nuclei (lasagna) and the columnar phase in the case of rod-like nuclei (spaghetti). An important difference compared with usual liquid crystals is that the nucleons are superfluid. We develop the hydrodynamic equations for this system and use them to study collective oscillations. Nucleon superfluidity leads to important qualitative differences in the spectra of these oscillations and also increases their frequencies compared with ordinary liquid crystals. We discuss a number of directions for future work.


title: Identification of SDSSJ141324.27+530527.0 as A New "Changing-Look" Quasar with a "Turn-on" Transition

abstract: We here report an identification of SDSSJ141324+530527.0 (SBS1411+533) at $z=0.456344$ as a new "changing-look" quasar with a "turn-on" spectral type transition from Type-1.9/2 to Type-1 within a rest frame time scale of 1-10 yr by a comparison of our new spectroscopic observation and the Sloan Digital Sky Survey (SDSS) archive data base. The SDSS DR7 spectrum taken in 2003 is dominated by a starlight emission from host galaxies redward of the Balmer limit, and has non-detectable broad H$\beta$ line. The new spectrum taken by us on June 1st, 2017 and SDSS DR14 spectrum taken on May 29, 2017 indicate that the object is of a typical quasar spectrum with a blue continuum and strong Balmer broad emission lines. In addition, an intermediate spectral type can be identified in the SDSS DR13 spectrum taken in 2015. The invariability of the line wing of MgII$\lambda2800$ emission and time scale argument (The invariability of [OIII]$\lambda$5007 line blue asymmetry) suggests that a variation of obscuration (an accelerating outflow) is not a favorable scenario. The time scale argument allows us to believe the type transition is possibly caused by either a viscous radial inflow or a disk instability around a $\sim5-9\times10^{7}M_\odot$ black hole.


title: Study for material analogs of FeSb$_{2}$: material design for thermoelectric materials

abstract: Using the \emph{ab initio} evolutionary algorithm (implemented in USPEX) and electronic structure calculations we investigate the properties of a new thermoelectric material FeSbAs, which is a material analog of the enigmatic thermoelectric FeSb$_{2}$. We utilize the density functional theory and the Gutzwiller method to check the energetics. We find that FeSbAs can be made thermodynamically stable above $\sim$30 GPa. We investigate the electronic structure and thermoelectric properties of FeSbAs based on the density functional theory and compare with those of FeSb$_{2}$. Above 50 K, FeSbAs has higher Seebeck coefficients than FeSb$_{2}$. Upon doping, the figure of merit becomes larger for FeSbAs than for FeSb$_{2}$. Another material analog FeSbP, was also investigated, and found thermodynamically unstable even at very high pressure. Regarding FeSb$_{2}$ as a member of a family of compounds (FeSb$_{2}$, FeSbAs, and FeSbP) we elucidate what are the chemical handles that control the gaps in this series. We also investigate solubility (As or P for Sb in FeSb$_{2}$) we found As to be more soluble. Finally, we study a two-band model for thermoelectric properties and find that the temperature dependent chemical potential and the presence of the ionized impurities are important to explain the extremum in the Seebeck coefficient exhibited in experiments for FeSb$_{2}$.


title: Extreme Khovanov spectra

abstract: We prove that the spectrum constructed by Gonz\'alez-Meneses, Manch\'on and the second author is stably homotopy equivalent to the Khovanov spectrum of Lipshitz and Sarkar at its extreme quantum grading.


title: Online Controlled Experiments for Personalised e-Commerce Strategies: Design, Challenges, and Pitfalls

abstract: Online controlled experiments are the primary tool for measuring the causal impact of product changes in digital businesses. It is increasingly common for digital products and services to interact with customers in a personalised way. Using online controlled experiments to optimise personalised interaction strategies is challenging because the usual assumption of statistically equivalent user groups is violated. Additionally, challenges are introduced by users qualifying for strategies based on dynamic, stochastic attributes. Traditional A/B tests can salvage statistical equivalence by pre-allocating users to control and exposed groups, but this dilutes the experimental metrics and reduces the test power. We present a stacked incrementality test framework that addresses problems with running online experiments for personalised user strategies. We derive bounds that show that our framework is superior to the best simple A/B test given enough users and that this condition is easily met for large scale online experiments. In addition, we provide a test power calculator and describe a selection of pitfalls and lessons learnt from our experience using it.


title: Oncilla robot: a versatile open-source quadruped research robot with compliant pantograph legs

abstract: We present Oncilla robot, a novel mobile, quadruped legged locomotion machine. This large-cat sized, 5.1 robot is one of a kind of a recent, bioinspired legged robot class designed with the capability of model-free locomotion control. Animal legged locomotion in rough terrain is clearly shaped by sensor feedback systems. Results with Oncilla robot show that agile and versatile locomotion is possible without sensory signals to some extend, and tracking becomes robust when feedback control is added (Ajaoolleian 2015). By incorporating mechanical and control blueprints inspired from animals, and by observing the resulting robot locomotion characteristics, we aim to understand the contribution of individual components. Legged robots have a wide mechanical and control design parameter space, and a unique potential as research tools to investigate principles of biomechanics and legged locomotion control. But the hardware and controller design can be a steep initial hurdle for academic research. To facilitate the easy start and development of legged robots, Oncilla-robot's blueprints are available through open-source. [...]


title: Stability, well-posedness and blow-up criterion for the Incompressible Slice Model

abstract: In atmospheric science, slice models are frequently used to study the behaviour of weather, and specifically the formation of atmospheric fronts, whose prediction is fundamental in meteorology. In 2013, Cotter and Holm introduced a new slice model, which they formulated using Hamilton's variational principle, modified for this purpose. In this paper, we show the local existence and uniqueness of strong solutions of the related ISM (Incompressible Slice Model). The ISM is a modified version of the Cotter-Holm Slice Model (CHSM) that we have obtained by adapting the Lagrangian function in Hamilton's principle for CHSM to the Euler-Boussinesq Eady incompressible case. Besides proving local existence and uniqueness, in this paper we also construct a blow-up criterion for the ISM, and study Arnold's stability around a restricted class of equilibrium solutions. These results establish the potential applicability of the ISM equations in physically meaningful situations.


title: Vector-valued higher depth quantum modular forms and higher Mordell integrals

abstract: In this paper, we prove vector-valued higher depth quantum modular properties arising from characters of certain vertex algebras. We then find two-dimensional Mordell integral representations for their errors of modularity.


title: Heterogeneous Doppler Spread-based CSI Estimation Planning for TDD Massive MIMO

abstract: Massive multi-input multi-output (Massive MIMO) has been recognized as a key technology to meet the demand for higher data capacity and massive connectivity. Nevertheless, the number of active users is restricted due to training overhead and the limited coherence time. Current wireless systems assume the same coherence slot duration for all users, regardless of their heterogeneous Doppler spreads. In this paper, we exploit this neglected degree of freedom in addressing the training overhead bottleneck. We propose a new uplink training scheme where the periodicity of pilot transmission differs among users based on their actual channel coherence times. Since the changes in the wireless channel are, primarily, due to movement, uplink training decisions are optimized, over long time periods, while considering the evolution of the users channels and locations. Owing to the different rates of the wireless channel and location evolution, a two time scale control problem is formulated. In the fast time scale, an optimal training policy is derived by choosing which users are requested to send their pilots. In the slow time scale, location estimation decisions are optimized. Simulation results show that the derived training policies provide a considerable improvement of the cumulative average spectral efficiency even with partial location knowledge.


title: A brief guide to reversing and extended symmetries of dynamical systems

abstract: The reversing symmetry group is a well-studied extension of the symmetry group of a dynamical system, the latter being defined by the action of a single homeomorphism on a topological space. While it is traditionally considered in nonlinear dynamics, where the space is simple but the map is complicated, it has an interesting counterpart in symbolic dynamics, where the map is simple but the space is not. Moreover, there is an interesting extension to the case of higher-dimensional shifts, where a similar concept can be introduced via the centraliser and the normaliser of the acting group in the full automorphism group of the shift space. We recall the basic notions and review some of the known results, in a fairly informal manner, to give a first impression of the phenomena that can show up in the extension from the centraliser to the normaliser, with some emphasis on recent developments.


title: Characterization of Strict Positive Definiteness on products of complex spheres

abstract: In this paper we consider Positive Definite functions on products $\Omega_{2q}\times\Omega_{2p}$ of complex spheres, and we obtain a condition, in terms of the coefficients in their disc polynomial expansions, which is necessary and sufficient for the function to be Strictly Positive Definite. The result includes also the more delicate cases in which $p$ and/or $q$ can be $1$ or $\infty$. The condition we obtain states that a suitable set in $\mathbb{Z}^2$, containing the indexes of the strictly positive coefficients in the expansion, must intersect every product of arithmetic progressions.


title: A quantum photonics model for non-classical light generation using integrated nanoplasmonic cavity-emitter systems

abstract: The implementation of non-classical light sources is becoming increasingly important for various quantum applications. A particularly interesting approach is to integrate such functionalities on a single chip as this could pave the way towards fully scalable quantum photonic devices. Several approaches using dielectric systems have been investigated in the past. However, it is still not understood how on-chip nanoplasmonic antennas, interacting with a single quantum emitter, affect the quantum statistics of photons reflected or transmitted in the guided mode of a waveguide. Here we investigate a quantum photonic platform consisting of an evanescently coupled nanoplasmonic cavity-emitter system and discuss the requirements for non-classical light generation. We develop an analytical model that incorporates quenching due to the nanoplasmonic cavity to predict the quantum statistics of the transmitted and reflected guided waveguide light under weak coherent pumping. The analytical predictions match numerical simulations based on a master equation approach. It is moreover shown that for resonant excitation the degree of anti-bunching in transmission is maximized for an optimal cavity modal volume $V_{c}$ and cavity-emitter distance $s$. In reflection, perfectly anti-bunched light can only be obtained for specific $(V_{c},s)$ combinations. Finally, our model also applies to dielectric cavities and as such can guide future efforts in the design and development of on-chip non-classical light sources using dielectric and nanoplasmonic cavity-emitter systems.


title: Compact Binary Mergers and The Event Rate of Fast Radio Bursts

abstract: Fast radio bursts (FRBs) are usually suggested to be associated with mergers of compact binaries consisting of white dwarfs (WDs), neutron stars (NSs), or black holes (BHs). We test these models by fitting the observational distributions in both redshift and isotropic energy of 22 Parkes FRBs, where, as usual, the rates of compact binary mergers (CBMs) are connected with cosmic star formation rates by a power-law distributed time delay. It is found that the observational distributions can well be produced by the CBM model with a characteristic delay time from several ten to several hundred Myr and an energy function index $1.2\lesssim\gamma\lesssim1.7$, where a tentative fixed spectral index $\beta=0.8$ is adopted for all FRBs. Correspondingly, the local event rate of FRBs is constrained to $(3-6)\times10^4f_{\rm b}^{-1}(\mathcal T/270\rm s)^{-1}(\mathcal A/2\pi)^{-1}\rm ~Gpc^{-3}yr^{-1}$ for an adopted minimum FRB energy of $E_{\min}=3\times10^{39}$ erg, where $f_{\rm b}$ is the beaming factor of the radiation, $\mathcal T$ is the duration of each pointing observation, and $\mathcal A $ is the sky area of the survey. This event rate, about an order of magnitude higher than the rates of NS-NS/NS-BH mergers, indicates that the most promising origin of FRBs in the CBM scenario could be mergers of WD-WD binaries. Here a massive WD could be produced since no FRB was found to be associated with a type Ia supernova. Alternatively, if actually all FRBs can repeat on a timescale much longer than the period of current observations, then they could also originate from a young active NS that forms from relatively rare NS-NS mergers and accretion-induced collapses of WD-WD binaries.


title: Consistent sets of lines with no colorful incidence

abstract: We consider incidences among colored sets of lines in $\mathbb{R}^d$ and examine whether the existence of certain concurrences between lines of $k$ colors force the existence of at least one concurrence between lines of $k+1$ colors. This question is relevant for problems in 3D reconstruction in computer vision.


title: On the critical energy required for homogeneous nucleation in bubble chambers employed in dark matter searches

abstract: Two equations for the calculation of the critical energy required for homogeneous nucleation in a superheated liquid, and the related critical radius of the nucleated vapour bubble, are obtained, the former by the direct application of the first law of thermodynamics, the latter by considering that the bubble formation implies the overcoming of a barrier of the free enthalpy potential. Comparisons with the currently used relationships demonstrate that the sensitivity of the bubble chambers employed in dark matter searches can be sometimes notably overestimated.


title: Crystal growth in nano-confinement: Subcritical cavity formation and viscosity effects

abstract: We report on the modeling of the formation of a cavity at the surface of crystals confined by a flat wall during growth in solution. Using a continuum thin film model, we discuss two phenomena that could be observed when decreasing the thickness of the liquid film between the crystal and the wall down to the nanoscale. First, in the presence of an attractive van der Waals contribution to the disjoining pressure, the formation of the cavity becomes sub-critical, i.e., discontinuous. In addition, there is a minimum supersaturation required to form a cavity. Second, when the thickness of the liquid film between the crystal and the substrate reaches the nanoscale, viscosity becomes relevant and hinders the formation of the cavity. We demonstrate that there is a critical value of the viscosity above which no cavity will form. The critical viscosity increases as the square of the thickness of the liquid film. A quantitative discussion of model materials such as Calcite, Sodium Chlorate, Glucose and Sucrose is provided.


title: Dirichlet problems for fully nonlinear equations with \lq subquadratic \lq Hamiltonians

abstract: For a class of fully nonlinear equations having second order operators which may be singular or degenerate when the gradient of the solutions vanishes, and having first order terms with power growth, we prove the existence and uniqueness of suitably defined viscosity solution of Dirichlet problem and we further show that it is a Lipschitz continuous function.


title: Maximal ideals in rings of real measurable functions

abstract: Let $ M (X)$ be the ring of all real measurable functions on a measurable space $(X, \mathscr{A})$. In this article, we show that every ideal of $M(X)$ is a $Z^{\circ}$-ideal. Also, we give several characterizations of maximal ideals of $M(X)$, mostly in terms of certain lattice-theoretic properties of $\mathscr{A}$. The notion of $T$-measurable space is introduced. Next, we show that for every measurable space $(X,\mathscr{A})$ there exists a $T$-measurable space $(Y,\mathscr{A}^{\prime})$ such that $M(X)\cong M(Y)$ as rings. The notion of compact measurable space is introduced. Next, we prove that if $(X, \mathscr{A})$ and $(Y, \mathfrak{M^{\prime}})$ are two compact $T$-measurable spaces, then $X\cong Y$ as measurable spaces if and only if $M(X)\cong M (Y)$ as rings.


title: Graph Partition Neural Networks for Semi-Supervised Classification

abstract: We present graph partition neural networks (GPNN), an extension of graph neural networks (GNNs) able to handle extremely large graphs. GPNNs alternate between locally propagating information between nodes in small subgraphs and globally propagating information between the subgraphs. To efficiently partition graphs, we experiment with several partitioning algorithms and also propose a novel variant for fast processing of large scale graphs. We extensively test our model on a variety of semi-supervised node classification tasks. Experimental results indicate that GPNNs are either superior or comparable to state-of-the-art methods on a wide variety of datasets for graph-based semi-supervised classification. We also show that GPNNs can achieve similar performance as standard GNNs with fewer propagation steps.


title: Donor and Acceptor Characteristics of Native Point Defects in GaN

abstract: The semiconducting behaviour and optoelectronic response of gallium nitride is governed by point defect processes, which, despite many years of research, remain poorly understood. The key difficulty in the description of the dominant charged defects is determining a consistent position of the corresponding defect levels, which is difficult to derive using standard supercell calculations. In a complementary approach, we take advantage of the embedded cluster methodology that provides direct access to a common zero of the electrostatic potential for all point defects in all charge states. Charged defects polarise a host dielectric material with long-range forces that strongly affect the outcome of defect simulations; to account for the polarisation we couple embedding with the hybrid quantum mechanical/molecular mechanical (QM/MM) approach and investigate the structure, formation and ionisation energies, and equilibrium concentrations of native point defects in wurtzite GaN at a chemically accurate hybrid-density-functional-theory level. N vacancies are the most thermodynamically favourable native defects in GaN, which contribute to the n-type character of as-grown GaN but are not the main source, a result that is consistent with experiment. Our calculations show no native point defects can form thermodynamically stable acceptor states. GaN can be easily doped n-type, but, in equilibrium conditions at moderate temperatures acceptor dopants will be compensated by N vacancies and no significant hole concentrations will be observed, indicating non-equilibrium processes must dominate in p-type GaN. We identify spectroscopic signatures of native defects in the infrared, visible and ultraviolet luminescence ranges and complementary spectroscopies. Crucially, we calculate the effective-mass-like-state levels associated with electrons and holes bound in diffuse orbitals...


title: Frequency-modulated SSFP with radial sampling and subspace reconstruction: A time-efficient alternative to phase-cycled bSSFP

abstract: Purpose: A novel subspace-based reconstruction method for frequency-modulated balanced steady-state free precession (fmSSFP) MRI is presented. In this work, suitable data acquisition schemes, subspace sizes, and efficiencies for banding removal are investigated. Theory and Methods: By combining a fmSSFP MRI sequence with a 3D stack-of-stars trajectory, scan efficiency is maximized as spectral information is obtained without intermediate preparation phases. A memory-efficient reconstruction routine is implemented by introducing the low-frequency Fourier transform as a subspace which allows for the formulation of a convex reconstruction problem. The removal of banding artifacts is investigated by comparing the proposed acquisition and reconstruction technique to phase-cycled bSSFP MRI. Aliasing properties of different undersampling schemes are analyzed and water/fat separation is demonstrated by reweighting the reconstructed subspace coefficients to generate virtual spectral responses in a post-processing step. Results: A simple root-of-sum-of-squares combination of the reconstructed subspace coefficients yields high-SNR images with the characteristic bSSFP contrast but without banding artifacts. Compared to Golden-Angle trajectories, turn-based sampling schemes were superior in minimizing aliasing across reconstructed subspace coefficients. Water/fat separated images of the human knee were obtained by reweighting subspace coefficients. Conclusion: The novel subspace-based fmSSFP MRI technique emerges as a time-efficient alternative to phase-cycled bSSFP. The method does not need intermediate preparation phases, offers high SNR and avoids banding artifacts. Reweighting of the reconstructed subspace coefficients allows for generating virtual spectral responses with applications to water/fat separation.


title: Growth monitoring with sub-monolayer sensitivity via real time thermal conductance measurements

abstract: Growth monitoring during the early stages of film formation is of prime importance to understand the growth process, the microstructure and thus the overall layer properties. In this work, we demonstrate that phonons can be used as sensitive probes to monitor real time evolution of film microstructure during growth, from incipient clustering to continuous film formation. For that purpose, a silicon nitride membrane-based sensor has been fabricated to measure in-plane thermal conductivity of thin film samples. Operating with the 3{\omega}-V\"olklein method at low frequencies, the sensor shows an exceptional resolution down to {\Delta}({\kappa}*t)=0.065 nm*W/(m*K), enabling accurate measurements. Validation of the sensor performance is done with organic and metallic thin films. In both cases, at early stages of growth, we observe an initial reduction of the effective thermal conductance of the supporting amorphous membrane, K, related with the surface phonon scattering enhanced by the incipient nanoclusters formation. As clusters develop, K reaches a minimum at the percolation threshold. Subsequent island percolation produces a sharp increase of the conductance and once the surface coverage is completed K increases linearly with thickness The thermal conductivity of the deposited films is obtained from the variation of K with thickness.


title: Two-Layered Falsification of Hybrid Systems guided by Monte Carlo Tree Search

abstract: Few real-world hybrid systems are amenable to formal verification, due to their complexity and black box components. Optimization-based falsification---a methodology of search-based testing that employs stochastic optimization---is attracting attention as an alternative quality assurance method. Inspired by the recent works that advocate coverage and exploration in falsification, we introduce a two-layered optimization framework that uses Monte Carlo tree search (MCTS), a popular machine learning technique with solid mathematical and empirical foundations. MCTS is used in the upper layer of our framework; it guides the lower layer of local hill-climbing optimization, thus balancing exploration and exploitation in a disciplined manner.


title: Extreme Events: Mechanisms and Prediction

abstract: Extreme events, such as rogue waves, earthquakes and stock market crashes, occur spontaneously in many dynamical systems. Because of their usually adverse consequences, quantification, prediction and mitigation of extreme events are highly desirable. Here, we review several aspects of extreme events in phenomena described by high-dimensional, chaotic dynamical systems. We specially focus on two pressing aspects of the problem: (i) Mechanisms underlying the formation of extreme events and (ii) Real-time prediction of extreme events. For each aspect, we explore methods relying on models, data or both. We discuss the strengths and limitations of each approach as well as possible future research directions.

categories: nlin.CD math.DS physics.ao-ph physics.flu-dyn


title: Particle-in-Cell Modeling of Laser Thomson Scattering in Low-Density Plasmas at Elevated Laser Intensities

abstract: Incoherent Thomson scattering is a non-intrusive technique commonly used for measuring local plasma density. Within low-density, low-temperature plasma's and for sufficient laser intensity, the laser may perturb the local electron density via the ponderomotive force, causing the diagnostic to become intrusive and leading to erroneous results. A theoretical model for this effect is validated numerically via kinetic simulations of a quasi-neutral plasma using the Particle-in-Cell technique.


title: On the uniqueness of the steady-state solution of the Lindblad-Gorini-Kossakowski-Sudarshan equation

abstract: The aims of this paper are two. The first is to give a brief review of the most relevant theoretical results concerning the uniqueness of the steady-state solution of the Lindblad-Gorini-Kossakowski-Sudarshan master equation and the criteria which guarantee relaxingness and irreducibility of dynamical semigroups. In particular, we test and discuss their physical meaning by considering their applicability to the characterisation of the simplest open quantum system \emph{i.e.} a two-level system coupled to a bath of harmonic oscillators at zero temperature. The second aim is to provide a set of sufficient conditions which guarantees the uniqueness of the steady-state solution and its attractivity. Starting from simple assumptions, we derive simple criteria that can be efficiently exploited to characterise the behavior of dissipative systems of spins and bosons (with truncated Fock space), and a wide variety of other open quantum systems recently studied.


title: Doping and tilting on optics in noncentrosymmetric multi-Weyl semimetals

abstract: Weyl semimetal (WSM) feature tilted Dirac cones and can be type I or II depending on the magnitude of the tilt parameter ($C$). The boundary between the two types is at $C=1$ where the cones are tipped and there is a Lifshitz transition. The topological charge of a WSM is one. In multi-Weyl it can be two or more depending on the value of the winding number $J$. We calculate the absorptive part of the AC optical conductivity both along the tilt direction ($\sigma_{zz}$) and perpendicular to it ($\sigma_{xx}$) as a function of the tilt ($C$) and chemical potential ($\mu$). For zero tilt there is a discontinuous rise in both $\sigma_{xx}$ and $\sigma_{zz}$ at photon energy $\Omega=2\mu$ followed by the usual linear in $\Omega$ law for $\sigma_{xx}$ at $J=1,2$ and $\sigma_{zz}$ at $J=1$. For $J=2$ and $\sigma_{zz}$ the interband background is constant rather than linear in $\Omega$. For type I there is a readjustment of optical spectral weight as the tilt is increased. The absorption starts from zero at $2\mu/(1+C)$ and then rises in a quasilinear fashion till it merge with the usual undoped untilted interband background at $2\mu/(1-C)$. The discontinuous rise at twice the chemical potential of the untilted case is lost. For type II the interband background of the undoped untilted case is never recovered. For noncentrosymmetric materials the energies of a pair of opposite chirality Weyl nodes become shifted by $\pm Q_{0}$ and this leads to two separate absorption edges corresponding to the effective chemical potential of each of the two nodes at $2(\mu+\chi Q_{0})$ depending on chirality $\chi=\pm$. We provide analytic expressions for the conductivity in this case which depend only on the ratio $Q_{0}/\mu$ and tilt when plotted against $\Omega/\mu$. The signature of finite energy shift $Q_{0}$ is more pronounced for $\sigma_{zz}$ and $J=2$ than for the other cases.


title: Description of 2-local derivations on some Lie rings of skew-adjoint matrices

abstract: In the present paper we prove that every 2-local inner derivation on the Lie ring of skew-symmetric matrices over a commutative ring is an inner derivation. We also apply our technique to various Lie algebras of infinite dimensional skew-adjoint matrix-valued maps on a set and prove that every 2-local spatial derivation on such algebras is a spatial derivation.


title: A mullti- or many- objective evolutionary algorithm with global loop update

abstract: Multi- or many-objective evolutionary algorithm- s(MOEAs), especially the decomposition-based MOEAs have been widely concerned in recent years. The decomposition-based MOEAs emphasize convergence and diversity in a simple model and have made a great success in dealing with theoretical and practical multi- or many-objective optimization problems. In this paper, we focus on update strategies of the decomposition- based MOEAs, and their criteria for comparing solutions. Three disadvantages of the decomposition-based MOEAs with local update strategies and several existing criteria for comparing solutions are analyzed and discussed. And a global loop update strategy and two hybrid criteria are suggested. Subsequently, an evolutionary algorithm with the global loop update is implement- ed and compared to several of the best multi- or many-objective optimization algorithms on two famous unconstraint test suites with up to 15 objectives. Experimental results demonstrate that unlike evolutionary algorithms with local update strategies, the population of our algorithm does not degenerate at any generation of its evolution, which guarantees the diversity of the resulting population. In addition, our algorithm wins in most instances of the two test suites, indicating that it is very compet- itive in terms of convergence and diversity. Running results of our algorithm with different criteria for comparing solutions are also compared. Their differences are very significant, indicating that the performance of our algorithm is affected by the criterion it adopts.


title: On hadron deformation: a model independent extraction of EMR from pion photoproduction data

abstract: The multipole content of pion photoproduction at the $\Delta^+ (1232)$ resonance has been extracted from a data set dominated by recent Mainz Microtron (MAMI) precision measurements. The analysis has been carried out in the Athens Model Independent Analysis Scheme (AMIAS), thus eliminating any model bias. The benchmark quantity for nucleon deformation, $EMR = E2/M1 = E_{1+}^{3/2}/M_{1+}^{3/2}$, was determined to be $-2.5 \pm 0.4_{stat+syst}$, thus reconfirming in a model independent way that the conjecture of baryon deformation is valid. The derived multipole amplitudes provide stringent constraints on QCD simulations and QCD inspired models striving to describe hadronic structure. They are in good agreement with phenomenological models which explicitly incorporate pionic degrees of freedom and with lattice QCD calculations.


title: Oriented cobordism of random manifolds

abstract: We introduce a general framework allowing the systematic study of random manifolds. In order to do so, we will put ourselves in a more general context than usual by allowing the underlying probability space to be non commutative. We introduce in this paper the oriented cobordism groups of random manifolds, which we compute in dimensions 0 and 1, and we prove the surjectivity of the corresponding Thom-Pontryagin homomorphism. Non commutative and usual (commutative) random manifolds are naturally related in this setting since two commutative random manifolds can be cobordant through a non commutative one, even if they are not cobordant as usual random manifolds. The main interest of our approach is that expected characteristic numbers can be generalized to the non commutative case and remain naturally invariant up to cobordism. This includes Hirzebruch signature, Pontryagin numbers and, more generally, the expected index of random elliptic differential operators


title: Visible Light-Activated Photosensitization of Hybridization of Far-red Fluorescent Protein and Silk

abstract: Fluorescent proteins often result in phototoxicity and cytotoxicity, in particular because some red fluorescent proteins (RFP) produce and release reactive oxygen species (ROS). The photogeneration of ROS is considered as a detrimental side effect in cellular imaging or is proactively utilized for ablating cancerous tissue. As ancient textiles or biomaterials, silk produced by silkworms can directly be used as fabrics or be processed into materials and structures to host other functional nanomaterials. We report that transgenic fusion of far-red fluorescent protein (mKate2) with silk provides a photosensitizer hybridization platform for photoinducible control of ROS. Taking advantage of green (visible) light activation, native and regenerated mKate2 silk can produce and release superoxide and singlet oxygen, in a comparable manner of visible light-driven plasmonic photocatalysis. Thus, the genetic expression of mKate2 in silk offers immediately exploitable and scalable photocatalyst-like biomaterials. We further envision that mKate2 silk could potentially rule out hazardous concerns associated with foreign semiconductor photocatalytic nanomaterials.


title: Nonequilibrium phase transitions of sheared colloidal microphases: Results from dynamical density functional theory

abstract: By means of classical density functional theory and its dynamical extension, we consider a colloidal fluid with spherically-symmetric competing interactions, which are well known to exhibit a rich bulk phase behavior. This includes complex three-dimensional periodically ordered cluster phases such as lamellae, two-dimensional hexagonally packed cylinders, gyroid structures or spherical micelles. While the bulk phase behavior has been studied extensively in earlier work, in this paper we focus on such structures confined between planar repulsive walls under shear flow. For sufficiently high shear rates, we observe that microphase separation can become fully suppressed. For lower shear rates, however, we find that e.g. the gyroid structure undergoes a kinetic phase transition to a hexagonally packed cylindrical phase, which is found experimentally and theoretically in amphiphilic block copolymer systems. As such, besides the known similarities between the latter and colloidal systems regarding the equilibrium phase behavior, our work reveals further intriguing non-equilibrium relations between copolymer melts and colloidal fluids with competing interactions.


title: Reduced Basis Kriging for Big Spatial Fields

abstract: In spatial statistics, a common method for prediction over a Gaussian random field (GRF) is maximum likelihood estimation combined with kriging. For massive data sets, kriging is computationally intensive, both in terms of CPU time and memory, and so fixed rank kriging has been proposed as a solution. The method however still involves operations on large matrices, so we develop an alteration to this method by utilizing the approximations made in fixed rank kriging combined with restricted maximum likelihood estimation and sparse matrix methodology. Experiments show that our methodology can provide additional gains in computational efficiency over fixed-rank kriging without loss of accuracy in prediction. The methodology is applied to climate data archived by the United States National Climate Data Center, with very good results.


title: ORGaNICs: A Theory of Working Memory in Brains and Machines

abstract: Working memory is a cognitive process that is responsible for temporarily holding and manipulating information. Most of the empirical neuroscience research on working memory has focused on measuring sustained activity in prefrontal cortex (PFC) and/or parietal cortex during simple delayed-response tasks, and most of the models of working memory have been based on neural integrators. But working memory means much more than just holding a piece of information online. We describe a new theory of working memory, based on a recurrent neural circuit that we call ORGaNICs (Oscillatory Recurrent GAted Neural Integrator Circuits). ORGaNICs are a variety of Long Short Term Memory units (LSTMs), imported from machine learning and artificial intelligence. ORGaNICs can be used to explain the complex dynamics of delay-period activity in prefrontal cortex (PFC) during a working memory task. The theory is analytically tractable so that we can characterize the dynamics, and the theory provides a means for reading out information from the dynamically varying responses at any point in time, in spite of the complex dynamics. ORGaNICs can be implemented with a biophysical (electrical circuit) model of pyramidal cells, combined with shunting inhibition via a thalamocortical loop. Although introduced as a computational theory of working memory, ORGaNICs are also applicable to models of sensory processing, motor preparation and motor control. ORGaNICs offer computational advantages compared to other varieties of LSTMs that are commonly used in AI applications. Consequently, ORGaNICs are a framework for canonical computation in brains and machines.


title: Spinodal superlattices of topological insulators

abstract: Spinodal decomposition is proposed for stabilizing self-assembled interfaces between topological insulators (TIs) by combining layers of iso-structural and iso-valent TlBi$X_2$ ($X$=S, Se, Te) materials. The composition range for gapless states is addressed concurrently to the study of thermodynamically driven boundaries. By tailoring composition, the TlBiS$_2$-TlBiTe$_2$ system might produce both spinodal superlattices and two dimensional eutectic microstructures, either concurrently or separately. The dimensions and topological nature of the metallic channels are determined by following the spatial distribution of the charge density and the spin-texture. The results validate the proof of concept for obtaining spontaneously forming two-dimensional TI-conducting channels embedded into three-dimensional insulating environments without any vacuum interfaces. Since spinodal decomposition is a controllable kinetic phenomenon, its leverage could become the long-sought enabler for effective TI technological deployment.


title: Theory of Single Charge Exchange Heavy Ion Reactions

abstract: The theory of heavy ion single charge exchange reactions is reformulated. In momentum space the reaction amplitude factorizes into a product of projectile and target transition form factors, folded with the nucleon-nucleon isovector interaction and a distortion coefficient which accounts for initial and final state ion-ion elastic interactions. The multipole structure of the transition form factors is studied in detail for Fermi-type non-spin flip and Gamow-Teller-type spin flip transitions, also serving to establish the connection to nuclear beta decay. The reaction kernel is evaluated for central and rank-2 tensor interactions. Initial and final state elastic ion-ion interaction are shown to be dominated by the imaginary part of the optical potential allowing to evaluate the reaction coefficients in the strong absorption limit, realized by the black disk approximation. In that limit the distortion coefficient is evaluated in closed form, revealing the relation to the total reaction cross section and the geometry of the transition form factors. It is shown that at small momentum transfer distortion effects reduce to a simple scaling factor, allowing to define reduced forward-angle cross section which is given by nuclear matrix elements of beta decay-type. The response function formalism is used to describe nuclear charge changing transitions. Spectral distributions obtained by a self-consistent HFB and QRPA approach are discussed for $\tau_\pm$ excitations of $^{18}O$ and $^{40}Ca$, respectively, and compared to spectroscopic data. The interplay of nuclear structure and reaction dynamics is illustrated for the single charge exchange reaction $^{18}O+^{40}Ca \to ^{18}F+^{40}K$ at $T_{lab}=270$ MeV.


title: Secure Wireless Powered and Cooperative Jamming D2D Communications

abstract: This paper investigates a secure wireless-powered device-to-device (D2D) communication network in the presence of multiple eavesdroppers, where a hybrid base station (BS) in a cellular network not only provides power wirelessly for the D2D transmitter to guarantee power efficiency for the D2D network, but also serves as a cooperative jammer (CJ) to interfere with the eavesdroppers. The cellular and D2D networks can belong to different service providers, which means that the D2D transmitter would need to pay for the energy service released by the hybrid BS to guarantee secure D2D communication. In order to exploit the hierarchical interaction between the BS and the D2D transmitter, we first formulate a Stackelberg game based energy trading scheme, where the quadratic energy cost model is considered. Then, a non-energy trading based Stackelberg game is investigated to study the reversed roles of the BS and the D2D users. For comparison, we also formulate and resolve the social welfare optimization problem. We derive the closed-form Stackelberg equilibriums of the formulated games and the optimal solutions for the social welfare optimization problem. Simulation results are provided to validate our proposed schemes to highlight the importance of energy trading interaction between cellular and D2D networks.


title: Search for high-mass resonances in dilepton final states in proton-proton collisions at $\sqrt{s}=$ 13 TeV

abstract: A search is presented for new high-mass resonances decaying into electron or muon pairs. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the CMS experiment at the LHC in 2016, corresponding to an integrated luminosity of 36 fb$^{-1}$. Observations are in agreement with standard model expectations. Upper limits on the product of a new resonance production cross section and branching fraction to dileptons are calculated in a model-independent manner. This permits the interpretation of the limits in models predicting a narrow dielectron or dimuon resonance. A scan of different intrinsic width hypotheses is performed. Limits are set on the masses of various hypothetical particles. For the Z$'_\mathrm{SSM}$ (Z$'_{\psi}$) particle, which arises in the sequential standard model (superstring-inspired model), a lower mass limit of 4.50 (3.90) TeV is set at 95% confidence level. The lightest Kaluza-Klein graviton arising in the Randall-Sundrum model of extra dimensions, with coupling parameters $k/\overline{M}_\mathrm{Pl}$ of 0.01, 0.05, and 0.10, is excluded at 95% confidence level below 2.10, 3.65, and 4.25 TeV, respectively. In a simplified model of dark matter production via a vector or axial vector mediator, limits at 95% confidence level are obtained on the masses of the dark matter particle and its mediator.


title: Measurement of e+e- --> D Dbar Cross Sections at the psi(3770) Resonance

abstract: We report new measurements of the cross sections for the production of D Dbar final states at the psi(3770) resonance. Our data sample consists of an integrated luminosity of 2.93/fb of e+e- annihilation data produced by the BEPCII collider and collected and analyzed with the BESIII detector. We exclusively reconstruct three D0 and six D+ hadronic decay modes and use the ratio of the yield of fully reconstructed D Dbar events ("double tags") to the yield of all reconstructed D or Dbar mesons ("single tags") to determine the number of D0 D0bar and D+D- events, benefiting from the cancellation of many systematic uncertainties. Combining these yields with an independent determination of the integrated luminosity of the data sample, we find the cross sections to be \sigma(e+e- --> D0 D0bar)=(3.615 +- 0.010 +- 0.038) nb and \sigma(e+e- --> D+D-)=(2.830 +- 0.011 +- 0.026) nb, where the uncertainties are statistical and systematic, respectively.


title: SDN for End-Nodes: Scenario Analysis and Architectural Guidelines

abstract: The advent of SDN has brought a plethora of new architectures and controller designs for many use-cases and scenarios. Existing SDN deployments focus on campus, datacenter and WAN networks. However, little research efforts have been devoted to the scenario of effectively controlling a full deployment of end-nodes (e.g. smartphones) that are transient and scattered across the Internet. In this paper, we present a rigorous analysis of the challenges associated with an SDN architecture for end-nodes, show that such challenges are not found in existing SDN scenarios, and provide practical design guidelines to address them. Then, and following these guidelines we present a reference architecture based on a decentralized, distributed and symmetric controller with a connectionless pull-oriented southbound and an intent-driven northbound. Finally, we measure a proof-of-concept deployment to assess the validity of the analysis as well as the architecture.


title: High-dimensional Stochastic Inversion via Adjoint Models and Machine Learning

abstract: Performing stochastic inversion on a computationally expensive forward simulation model with a high-dimensional uncertain parameter space (e.g. a spatial random field) is computationally prohibitive even with gradient information provided. Moreover, the `nonlinear' mapping from parameters to observables generally gives rise to non-Gaussian posteriors even with Gaussian priors, thus hampering the use of efficient inversion algorithms designed for models with Gaussian assumptions. In this paper, we propose a novel Bayesian stochastic inversion methodology, characterized by a tight coupling between a gradient-based Langevin Markov Chain Monte Carlo (LMCMC) method and a kernel principal component analysis (KPCA). This approach addresses the `curse-of-dimensionality' via KPCA to identify a low-dimensional feature space within the high-dimensional and nonlinearly correlated spatial random field. Moreover, non-Gaussian full posterior probability distribution functions are estimated via an efficient LMCMC method on both the projected low-dimensional feature space and the recovered high-dimensional parameter space. We demonstrate this computational framework by integrating and adapting recent developments such as data-driven statistics-on-manifolds constructions and reduction-through-projection techniques to solve inverse problems in linear elasticity.


title: Mass-invariance of the iron enrichment in the hot haloes of massive ellipticals, groups, and clusters of galaxies

abstract: X-ray measurements find systematically lower Fe abundances in the X-ray emitting haloes pervading groups ($kT\lesssim1.7$ keV) than in clusters of galaxies. These results have been difficult to reconcile with theoretical predictions. However, models using incomplete atomic data or the assumption of isothermal plasmas may have biased the best fit Fe abundance in groups and giant elliptical galaxies low. In this work, we take advantage of a major update of the atomic code in the spectral fitting package SPEX to re-evaluate the Fe abundance in 43 clusters, groups, and elliptical galaxies (the CHEERS sample) in a self-consistent analysis and within a common radius of 0.1$r_{500}$. For the first time, we report a remarkably similar average Fe enrichment in all these systems. Unlike previous results, this strongly suggests that metals are synthesised and transported in these haloes with the same average efficiency across two orders of magnitude in total mass. We show that the previous metallicity measurements in low temperature systems were biased low due to incomplete atomic data in the spectral fitting codes. The reasons for such a code-related Fe bias, also implying previously unconsidered biases in the emission measure and temperature structure, are discussed.


title: Mass-Metallicity Relation from Cosmological Hydrodynamical Simulations and X-ray Observations of Galaxy Groups and Clusters

abstract: Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($\beta\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.


title: Spectroscopic characterisation of the stellar content of ultra diffuse galaxies

abstract: Understanding the peculiar properties of Ultra Diffuse Galaxies (UDGs) via spectroscopic analysis is a challenging task requiring very deep observations and exquisite data reduction. In this work we perform one of the most complete characterisations of the stellar component of UDGs to date using deep optical spectroscopic data from OSIRIS at GTC. We measure radial and rotation velocities, star formation histories (SFH) and mean population parameters, such as ages and metallicities, for a sample of five UDG candidates in the Coma cluster. From the radial velocities, we confirm the Coma membership of these galaxies. We find that their rotation properties, if detected at all, are compatible with dwarf-like galaxies. The SFHs of the UDG are dominated by old (~ 7 Gyr), metal-poor ([M/H] ~ -1.1) and alpha-enhanced ([Mg/Fe] ~ 0.4) populations followed by a smooth or episodic decline which halted ~ 2 Gyr ago, possibly a sign of cluster-induced quenching. We find no obvious correlation between individual SFH shapes and any UDG morphological properties. The recovered stellar properties for UDGs are similar to those found for DDO44, a local UDG analogue resolved into stars. We conclude that the UDGs in our sample are extended dwarfs whose properties are likely the outcome of both internal processes, such as bursty SFHs and/or high-spin haloes, as well as environmental effects within the Coma cluster.


title: On the existence of a scalar pressure field in the Br\"odinger problem

abstract: This work deals with the entropic regularization of the Brenier problem for perfect incompressible fluids introduced by Arnaudon, Cruzeiro, L\'eonard and Zambrini. We show that as in the original setting, there exists a scalar pressure field which is interpreted as the Lagrange multiplier associated to the incompressibility constraint. The proof goes through a reformulation of the problem in PDE terms.