title: Bayesian Sparse Global-Local Shrinkage Regression for Selection of Grouped Variables

abstract: Most estimates for penalised linear regression can be viewed as posterior modes for an appropriate choice of prior distribution. Bayesian shrinkage methods, particularly the horseshoe estimator, have recently attracted a great deal of attention in the problem of estimating sparse, high-dimensional linear models. This paper extends these ideas, and presents a Bayesian grouped model with continuous global-local shrinkage priors to handle complex group hierarchies that include overlapping and multilevel group structures. As the posterior mean is never a sparse estimate of the linear model coefficients, we extend the recently proposed decoupled shrinkage and selection (DSS) technique to the problem of selecting groups of variables from posterior samples. To choose a final, sparse model, we also adapt generalised information criteria approaches to the DSS framework. To ensure that sparse groups, in which only a few predictors are active, can be effectively identified, we provide an alternative degrees of freedom estimator for sparse Bayesian linear models that takes into account the effects of shrinkage on the model coefficients. Simulations and real data analysis using our proposed method show promising performance in terms of correct identification of active and inactive groups, and prediction, in comparison with a Bayesian grouped slab-and-spike approach.


title: Realising Interactions Between Dark Matter and Dark Energy Using $k$-essence Cosmology

abstract: In this paper we exploit dynamics of a $k-$essence scalar field to realise interactions between dark components of universe resulting in a evolution consistent with observed features of late time phase of cosmic evolution. Stress energy tensor corresponding to a $k-$essence Lagrangian $L=V( \phi)F(X)$ (where $X=\frac{1}{2}g^{\mu\nu}\nabla_\mu\phi \nabla_\nu\phi$) is shown to be equivalent to an ideal fluid with two components having same equation of state. Stress energy tensor of one of the components may be generated from a constant potential $k-$essence Lagrangian of form $L_1=V_0F(X)$ ($V_0$ constant) and that of other from another Lagrangian of form $L_2=V_1(\phi)F(X)$ with $V=V_0 + V_1(\phi)$. We have shown that, the unified dynamics of dark matter and dark energy described by a single scalar field $\phi$ driven by a $k-$essence Lagrangian $L= V(\phi)F(X)$ may be viewed in terms of diffusive interactions between the two hypothetical fluid components `1' and `2' with stress energy tensors equivalent to that of Lagrangians $L_1$ and $L_2$ respectively. The energy transfer between the fluid components is determined by functions $V(\phi)$, $F(X)$ and their derivatives. Such a realisation is shown to be consistent with the Supernova Ia data with certain constraints on the temporal behaviour of $k-$essence potential $V(\phi)$. We have described a methodology to obtain such constraints.


title: The norm of the harmonic Bergman projection and Besov space

abstract: We estimate the norm of the harmonic Bergman projection in the context of harmonic Besov spaces. We obtain the two-side norm estimates in general Lp-case.


title: Endurance of Quantum Coherence in Born-Markov Open Quantum Systems

abstract: Quantum coherence, the physical property underlying fundamental phenomena such as multi-particle interference and entanglement, has emerged as a valuable resource upon which exotic modern technologies are founded. In general, the most prominent adversary of quantum coherence is noise arising from the interaction of the associated dynamical system with its environment. Under certain conditions, however, the existence of noise may drive quantum and classical systems to endure intriguing nontrivial effects. Along these lines, here we demonstrate, both theoretically and experimentally, that when two indistinguishable particles co-propagate through quantum networks affected by noise, the system always evolves into a steady state in which coherences between certain separable states perpetually prevail. Furthermore, we show that the same steady state with surviving quantum coherences is reached irrespectively of the configuration in which the particles are prepared.


title: Higher-order parabolic equations with VMO assumptions and general boundary conditions with variable leading coefficients

abstract: We prove weighted mixed $L_{p}(L_{q})$-estimates, with $p,q\in(1,\infty)$, for higher-order elliptic and parabolic equations on the half space $\mathbb{R}^{d+1}_{+}$ and on domains with general boundary conditions which satisfy the Lopatinskii--Shapiro condition. We assume that the elliptic operators $A$ have leading coefficients which are in the class of vanishing mean oscillations both in the time and the space variables, and that the boundary conditions have variable leading coefficients. The proofs are based on and generalize the estimates recently obtained by the authors in [DG17]. [DG17] H.Dong and C.Gallarati. Higher order elliptic and parabolic equations with VMO assumptions and general boundary conditions. Accepted for publication in Journal of Functional Analysis. See: arXiv:1702.03254.


title: The hydrogen atom according to relativistic wave mechanics -- amplitude functions and circulating electronic currents

abstract: The solution of Dirac's equation for the hydrogen atom according to relativistic wave mechanics yields for each state a vectorial amplitude function with four components, two large and two small. Each such component has its characteristic surface of constant amplitude, of which we plot several examples. For each state of the hydrogen atom there is both a density of electronic charge surrounding the atomic nucleus and an electronic current circulating about the polar axis; the latter generates a magnetic dipolar moment that agrees precisely with experiment.


title: Quantum, noncommutative and MOND corrections to the entropic law of gravitation

abstract: Quantum and noncommutative corrections to the Newtonian law of inertia are considered in the general setting of Verlinde's entropic force postulate. We demonstrate that the form for the modified Newtonian dynamics (MOND) emerges in a classical setting by seeking appropriate corrections in the entropy.We estimate the correction term by using concrete coherent states in the standard and generalized versions of Heisenberg's uncertainty principle. Using Jackiw's direct and analytic method we compute the explicit wavefunctions for these states producing minimal length as well as minimal products. Subsequently we derive a further selection criterion restricting the free parameters in the model in providing a canonical formulation of the quantum corrected Newtonian law by setting up the Lagrangian and Hamiltonian for the system.


title: Mean square of zeta function, circle problem and divisor problem revisited

abstract: This paper is closely related to the recent work [BW17] of the same authors and our purpose is to elaborate more on some of the results and methods from [BW17]. More specifically our goal is two-fold. Firstly, we will indicate how a simple variant related to Section 4 in [BW17] leads to the following improvements of Theorem 3 in [BW17]


title: Dipolar-glass behaviour of an insulating film containing nanogranular Fe particles

abstract: We perform Metropolis Monte Carlo simulations of the behaviour of a film of insulating material containing a distribution of magnetic nanoparticles. We assume that these particles only interact through dipolar forces and we find that their behaviour at low T shows characteristics of a spin-glass with a freezing T_f at which the linear susceptibility and the specific heat show a maximum. We obtain the spin glass order parameter as a function of temperature and we also calculate the time auto-correlation of the spin at the center of the system. We find that these results are consistent with the temperature dependence of the variance and the mean of the local field at the central spin.


title: Model Selection Confidence Sets by Likelihood Ratio Testing

abstract: The traditional activity of model selection aims at discovering a single model superior to other candidate models. In the presence of pronounced noise, however, multiple models are often found to explain the same data equally well. To resolve this model selection ambiguity, we introduce the general approach of model selection confidence sets (MSCSs) based on likelihood ratio testing. A MSCS is defined as a list of models statistically indistinguishable from the true model at a user-specified level of confidence, which extends the familiar notion of confidence intervals to the model-selection framework. Our approach guarantees asymptotically correct coverage probability of the true model when both sample size and model dimension increase. We derive conditions under which the MSCS contains all the relevant information about the true model structure. In addition, we propose natural statistics based on the MSCS to measure importance of variables in a principled way that accounts for the overall model uncertainty. When the space of feasible models is large, MSCS is implemented by an adaptive stochastic search algorithm which samples MSCS models with high probability. The MSCS methodology is illustrated through numerical experiments on synthetic data and real data examples.


title: End-to-End Audiovisual Fusion with LSTMs

abstract: Several end-to-end deep learning approaches have been recently presented which simultaneously extract visual features from the input images and perform visual speech classification. However, research on jointly extracting audio and visual features and performing classification is very limited. In this work, we present an end-to-end audiovisual model based on Bidirectional Long Short-Term Memory (BLSTM) networks. To the best of our knowledge, this is the first audiovisual fusion model which simultaneously learns to extract features directly from the pixels and spectrograms and perform classification of speech and nonlinguistic vocalisations. The model consists of multiple identical streams, one for each modality, which extract features directly from mouth regions and spectrograms. The temporal dynamics in each stream/modality are modeled by a BLSTM and the fusion of multiple streams/modalities takes place via another BLSTM. An absolute improvement of 1.9% in the mean F1 of 4 nonlingusitic vocalisations over audio-only classification is reported on the AVIC database. At the same time, the proposed end-to-end audiovisual fusion system improves the state-of-the-art performance on the AVIC database leading to a 9.7% absolute increase in the mean F1 measure. We also perform audiovisual speech recognition experiments on the OuluVS2 database using different views of the mouth, frontal to profile. The proposed audiovisual system significantly outperforms the audio-only model for all views when the acoustic noise is high.


title: Flexible Network Binarization with Layer-wise Priority

abstract: How to effectively approximate real-valued parameters with binary codes plays a central role in neural network binarization. In this work, we reveal an important fact that binarizing different layers has a widely-varied effect on the compression ratio of network and the loss of performance. Based on this fact, we propose a novel and flexible neural network binarization method by introducing the concept of layer-wise priority which binarizes parameters in inverse order of their layer depth. In each training step, our method selects a specific network layer, minimizes the discrepancy between the original real-valued weights and its binary approximations, and fine-tunes the whole network accordingly. During the iteration of the above process, it is significant that we can flexibly decide whether to binarize the remaining floating layers or not and explore a trade-off between the loss of performance and the compression ratio of model. The resulting binary network is applied for efficient pedestrian detection. Extensive experimental results on several benchmarks show that under the same compression ratio, our method achieves much lower miss rate and faster detection speed than the state-of-the-art neural network binarization method.


title: Monotonically controlled integrals

abstract: The monotonically controlled integral defined by Bendov\'a and Mal\'y, which is equivalent to the Denjoy-Perron integral, admits a natural parameter $\alpha>0$ thereby leading to the whole scale of integrals called $\alpha$-monotonically controlled integrals. While the power of these integrals is easily seen to increase with increasing $\alpha$, our main results show that their exact dependence on $\alpha$ is rather curious. For $\alpha<1$ they do not even contain the Lebesgue integral, for $1\le\alpha\le 2$ they coincide with the Denjoy-Perron integral, and for $\alpha>2$ they are mutually different and not even contained in the Denjoy-Khintchine integral.


title: A non-Standard Standard Model

abstract: We examine the Standard Model under the electroweak symmetry group $U_{EW}(2)$ subject to the Lie algebra condition $\mathfrak{u}_{EW}(2)\not\cong \mathfrak{su}_{I}(2)\oplus \mathfrak{u}_{Y}(1)$. Physically, the condition ensures that all electroweak gauge bosons interact with each other prior to symmetry breaking. This represents a crucial shift in the identification of physical gauge bosons: Unlike the Standard Model which posits a change of Lie algebra basis induced by spontaneous symmetry breaking, here the basis is unaltered and $A,\,Z^0,\,W^\pm$ represent the physical bosons both before and after spontaneous symmetry breaking. Our choice of $\mathfrak{u}_{EW}(2)$ requires some modification of the matter field representation of the Standard Model. For $U_{EW}(2)$, there are two pertinent representations ${\mathbf{2}}$ and its $U(2)$-conjugate ${\mathbf{2^c}}$ related by a global gauge transformation that squares to minus the identity. The product group structure calls for strong-electroweak degrees of freedom in the $(\mathbf{3},\mathbf{2})$ and the $(\mathbf{3},{\mathbf{2^c}})$ of $SU_C(3)\times U_{EW}(2)$ that possess integer electric charge just like leptons. These degrees of freedom play the role of quarks, and they lead to a modified Lagrangian that nevertheless reproduces transition rates and cross sections equivalent to the Standard Model. The close resemblance between quark and lepton electroweak doublets suggests a mechanism for a speculative phase transition between quarks and leptons that stems from the product structure of the symmetry group. Our hypothesis is that the strong and electroweak bosons see each other as a source of decoherence. In effect, lepton representations get identified with the $SU(3)$-trace-reduced quark representations. This mechanism allows for possible extensions of the Standard Model that don't require large inclusive multiplets of matter fields.


title: Zoom Out-and-In Network with Map Attention Decision for Region Proposal and Object Detection

abstract: In this paper, we propose a zoom-out-and-in network for generating object proposals. A key observation is that it is difficult to classify anchors of different sizes with the same set of features. Anchors of different sizes should be placed accordingly based on different depth within a network: smaller boxes on high-resolution layers with a smaller stride while larger boxes on low-resolution counterparts with a larger stride. Inspired by the conv/deconv structure, we fully leverage the low-level local details and high-level regional semantics from two feature map streams, which are complimentary to each other, to identify the objectness in an image. A map attention decision (MAD) unit is further proposed to aggressively search for neuron activations among two streams and attend the most contributive ones on the feature learning of the final loss. The unit serves as a decisionmaker to adaptively activate maps along certain channels with the solely purpose of optimizing the overall training loss. One advantage of MAD is that the learned weights enforced on each feature channel is predicted on-the-fly based on the input context, which is more suitable than the fixed enforcement of a convolutional kernel. Experimental results on three datasets, including PASCAL VOC 2007, ImageNet DET, MS COCO, demonstrate the effectiveness of our proposed algorithm over other state-of-the-arts, in terms of average recall (AR) for region proposal and average precision (AP) for object detection.


title: Natural Language Inference over Interaction Space

abstract: Natural Language Inference (NLI) task requires an agent to determine the logical relationship between a natural language premise and a natural language hypothesis. We introduce Interactive Inference Network (IIN), a novel class of neural network architectures that is able to achieve high-level understanding of the sentence pair by hierarchically extracting semantic features from interaction space. We show that an interaction tensor (attention weight) contains semantic information to solve natural language inference, and a denser interaction tensor contains richer semantic information. One instance of such architecture, Densely Interactive Inference Network (DIIN), demonstrates the state-of-the-art performance on large scale NLI copora and large-scale NLI alike corpus. It's noteworthy that DIIN achieve a greater than 20% error reduction on the challenging Multi-Genre NLI (MultiNLI) dataset with respect to the strongest published system.


title: Superconductivity in Engineered Two-Dimensional Electron Gases

abstract: We consider Kohn-Luttinger mechanism for superconductivity in a two-dimensional electron gas confined to a narrow well between two metallic planes with two occupied subbands with Fermi momenta $k_{FL} > k_{FS}$. On the basis of a perturbative analysis, we conclude that non-s-wave superconductivity emerges even when the bands are parabolic. We analyze the conditions that maximize $T_c$ as a function of the distance to the metallic planes, the ratio $k_{FL}/k_{FS}$, and $r_s$, which measures the strength of Coulomb correlations. The largest attraction is in p-wave and d-wave channels, of which p-wave is typically the strongest. For $r_s = O(1)$ we estimate that the dimensionless coupling $\lambda \approx 10^{-1}$, but it likely continues increasing for larger $r_s$ (where we lose theoretical control).


title: Towards microscopic studies of survival probabilities of compound superheavy nuclei

abstract: The microscopic approach of fission rates and neutron emission rates in compound nuclei have been applied to $^{258}$No and $^{286}$Cn. The microscopic framework is based on the finite-temperature Skyrme-Hartree-Fock+BCS calculations, in which the fission barriers and mass parameters are self-consistently temperature dependent. The fission rates from low to high temperatures can be obtained based on the imaginary free energy method. The neutron emission rates are obtained with neutron gases at surfaces. Finally the survival probabilities of superheavy nuclei can be calculated microscopically. The microscopic approach has been compared with the widely used statistical models. Generally, there are still large uncertainties in descriptions of fission rates.


title: A posteriori error analysis for random scalar conservation laws using the Stochastic Galerkin method

abstract: In this article we present an a posteriori error estimator for the spatial-stochastic error of a Galerkin-type discretisation of an initial value problem for a random hyperbolic conservation law. For the stochastic discretisation we use the Stochastic Galerkin method and for the spatial-temporal discretisation of the Stochastic Galerkin system a Runge-Kutta Discontinuous Galerkin method. The estimator is obtained using smooth reconstructions of the discrete solution. Combined with the relative entropy stability framework of Dafermos \cite{dafermos2005hyperbolic}, this leads to computable error bounds for the space-stochastic discretisation error. Moreover, it turns out that the error estimator admits a splitting into one part representing the spatial error, and a remaining term, which can be interpreted as the stochastic error. This decomposition allows us to balance the errors arising from spatial and stochastic discretisation. We conclude with some numerical examples confirming the theoretical findings.


title: Comments on Gravitoelectromagnetism of Ummarino and Gallerati in "Superconductor in a weak static gravitational field" vs Other Versions

abstract: Recently reported [Eur. Phys. J. C., ${\bf 77}$, 549 (2017)] gravitoelectromagnetic equations of Ummarino and Gallerati (UG) in their linearized version of General Relativity (GR) are shown to match with (a) our previously reported special relativistic Maxwellian Gravity equations in the non-relativistic limit and with (b) the non-relativistic equations derived here, when the speed of gravity $c_g$ (an undetermined parameter of the theory here) is set equal to $c$ (the speed of light in vacuum). Seen in the light of our new results, the UG equations satisfy the Correspondence Principle (cp), while many other versions of linearized GR equations that are being (or may be) used to interpret the experimental data defy the cp. Such new findings assume significance and relevance in the contexts of recent detection of gravitational waves and the gravitomagnetic field of the spinning earth and their interpretations. Being well-founded and self-consistent, the equations may be of interest and useful to researchers exploring the phenomenology of gravitomagnetism, gravitational waves and the novel interplay of gravity with different states of matter in flat space-time like UG's interesting work on superconductors in weak gravitational fields.


title: Stochasticity in the 21cm power spectrum at the epoch of reionization and cosmic dawn

abstract: The 21cm neutral hydrogen line is likely to be a key probe for studying the epoch of reionization and comic dawn in the forthcoming decades. This prospect stimulates the development of the theoretical basis for simulating the power spectrum of this line. Because of the beam size of the upcoming radio telescopes at high redshifts, most of the theoretical models are focused on the inhomogeneities on scales above few comoving megaparsecs. Therefore, smaller scales are often neglected and modeled with approximated sub-grid models. In this study we explore whether the perturbations on small scales ($\lesssim 1h^{-1}\mathrm{Mpc}$) can affect the 21cm signal on larger scales. Two distinct mechanism are discussed. First, we show that during the cosmic dawn small scale perturbations regulate the formation time of the first Lyman-alpha Emitters (LAE), and consequently the coupling timing of spin and kinetic temperatures. Due to the low number density of LAE, the inhomogeneity of coupling includes the shot noise and manifests itself in the observed 21cm power spectrum. Second mechanism works during the reionization when the ionized bubbles actively grow and overlap. Small scales perturbations affect the galactic properties and merger histories, and consequently the number of ionizing photons produced by each galaxy. The ionizing photons bring the perturbations from the galactic scales to the scales of ionizing fronts, affecting the 21cm power spectrum. We conclude that these two effects introduce stochasticity in the potentially observed 21cm power spectrum and, moreover, might give another perspective into the physics of the first galaxies.


title: Observing the Galactic Plane with the Cherenkov Telescope Array

abstract: The Cherenkov Telescope Array is a next generation ground-based gamma-ray observatory de- signed to detect photons in the 20 GeV to 300 TeV energy range. With a sensitivity improvement of up to one order of magnitude on the entire energy range with respect to currently operating facilities, coupled with significantly better angular resolution, the array will be used to address many open questions in high-energy astrophysics. In addition, CTA will explore the ultra-high energy (E >50 TeV) window with great sensitivity for the first time. CTA is expected to reveal a detailed picture of the Galactic plane at the highest energies, and to discover around one hundred new supernova remnants and many hundreds of pulsar wind nebulae, according to current population estimates. The ability of the observatory to resolve such a large number of Galactic sources is one of the challenges to be faced. In this paper, we will present the first simulated scan of the Galactic plane with a realistic observation strategy, with particular attention to the potential source confusion. We will also present prospects for morphological studies of extended sources, such as the young SNR RX J1713.7-39.


title: Gaussian multiplicative chaos through the lens of the 2D Gaussian free field

abstract: The aim of this review-style paper is to provide a concise, self-contained and unified presentation of the construction and main properties of Gaussian multiplicative chaos (GMC) measures for log-correlated fields in 2D in the subcritical regime. By considering the case of the 2D Gaussian free field, we review convergence, uniqueness and characterisations of the measures; revisit Kahane's convexity inequalities and existence and scaling of moments; discuss the measurability of the underlying field with respect to the GMC measure and present a KPZ relation for scaling exponents.


title: Mechanics of ELM control coil induced alpha particle transport

abstract: Using the orbit-following code ASCOT, we model alpha particle transport in ITER under the influence of ELM control coils (ECCs), toroidal field ripple, and test blanket modules, with emphasis on how the plasma response (PR) modifies the transport mechanisms and fast ion loads on the divertor. We found that while PR shields the plasma by healing broken flux surfaces, it also opens a new loss channel for marginally trapped particles: PR causes strong toroidal variation of the poloidal field near the X-point which leads to de-localisation of banana tips and collisionless transport. The reduction in passing particle losses and the increase in marginally trapped particle losses shift divertor loads from targets to the dome and under-the-dome structures. The plasma response was calculated by both MARS-F and JOREK codes. The new transport mechanism was stronger for PR calculated by JOREK which, unlike MARS-F, explicitly includes the X-point.


title: Zeros of the deformed exponential function

abstract: Let $f(x)=\sum_{n=0}^{\infty}\frac{1}{n!}q^{n(n-1)/2}x^n$ ($0<q<1$) be the deformed exponential function. It is known that the zeros of $f(x)$ are real and form a negative decreasing sequence $(x_k)$ ($k\ge 1$). We investigate the complete asymptotic expansion for $x_{k}$ and prove that for any $n\ge1$, as $k\to \infty$, \begin{align*} x_k=-kq^{1-k}\Big(1+\sum_{i=1}^{n}C_i(q)k^{-1-i}+o(k^{-1-n})\Big), \end{align*} where $C_i(q)$ are some $q$ series which can be determined recursively. We show that each $C_{i}(q)\in \mathbb{Q}[A_0,A_1,A_2]$, where $A_{i}=\sum_{m=1}^{\infty}m^i\sigma(m)q^m$ and $\sigma(m)$ denotes the sum of positive divisors of $m$. When writing $C_{i}$ as a polynomial in $A_0, A_1$ and $A_2$, we find explicit formulas for the coefficients of the linear terms by using Bernoulli numbers. Moreover, we also prove that $C_{i}(q)\in \mathbb{Q}[E_2,E_4,E_6]$, where $E_2$, $E_4$ and $E_6$ are the classical Eisenstein series of weight 2, 4 and 6, respectively.


title: An additional structure over integer rings $\mathbb{Z}_{p^r}^n$

abstract: We present an algebraic structure in modules over integer rings with cardinality prime powers, which allows to define bases. With such structure, we prove a similar version for the basis extension theorem of linear algebra over fields. Moreover, we exhibit results involving the modules and their duals.


title: Linguistic Features of Genre and Method Variation in Translation: A Computational Perspective

abstract: In this paper we describe the use of text classification methods to investigate genre and method variation in an English - German translation corpus. For this purpose we use linguistically motivated features representing texts using a combination of part-of-speech tags arranged in bigrams, trigrams, and 4-grams. The classification method used in this paper is a Bayesian classifier with Laplace smoothing. We use the output of the classifiers to carry out an extensive feature analysis on the main difference between genres and methods of translation.


title: An overview of the escape dynamics in the Henon-Heiles Hamiltonian system

abstract: The aim of this work is to revise but also explore even further the escape dynamics in the H\'{e}non-Heiles system. In particular, we conduct a thorough and systematic numerical investigation distinguishing between trapped (ordered and chaotic) and escaping orbits, considering only unbounded motion for several energy levels. It is of particular interest, to locate the basins of escape towards the different escape channels and relate them with the corresponding escape periods of the orbits. In order to elucidate the escape process we conduct a thorough investigation in several types of two-dimensional planes and also in a three-dimensional subspace of the entire four-dimensional phase space. We classify extensive samples of orbits by integrating numerically the equations of motion as well as the variational equations. In an attempt to determine the regular or chaotic nature of trapped motion, we apply the SALI method, as an accurate chaos detector. It was found, that in all studied cases regions of non-escaping orbits coexist with several basins of escape. Most of the current outcomes have been compared with previous related work.


title: Distributed Scheduling in Time Dependent Environments: Algorithms and Analysis

abstract: Consider the problem of a multiple access channel in a time dependent environment with a large number of users. In such a system, mostly due to practical constraints (e.g., decoding complexity), not all users can be scheduled together, and usually only one user may transmit at any given time. Assuming a distributed, opportunistic scheduling algorithm, we analyse the system's properties, such as delay, QoS and capacity scaling laws. Specifically, we start with analyzing the performance while \emph{assuming the users are not necessarily fully backlogged}, focusing on the queueing problem and, especially, on the \emph{strong dependence between the queues}. We first extend a known queueing model by Ephremides and Zhu, to give new results on the convergence of the probability of collision to its average value (as the number of users grows), and hence for the ensuing system performance metrics, such as throughput and delay. This model, however, is limited in the number of users one can analyze. We thus suggest a new model, which is much simpler yet can accurately describes the system behaviour when the number of users is large. We then proceed to the analysis of this system under the assumption of time dependent channels. Specifically, we assume each user experiences a different channel state sequence, expressing different channel fluctuations (specifically, the Gilbert-Elliott model). The system performance under this setting is analysed, along with the channel capacity scaling laws.


title: Data-Driven Filtered Reduced Order Modeling Of Fluid Flows

abstract: We propose a data-driven filtered reduced order model (DDF-ROM) framework for the numerical simulation of fluid flows. The novel DDF-ROM framework consists of two steps: (i) In the first step, we use explicit ROM spatial filtering of the nonlinear PDE to construct a filtered ROM. This filtered ROM is low-dimensional, but is not closed (because of the nonlinearity in the given PDE). (ii) In the second step, we use data-driven modeling to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved modes. To this end, we use a quadratic ansatz to model this interaction and close the filtered ROM. To find the new coefficients in the closed filtered ROM, we solve an optimization problem that minimizes the difference between the full order model data and our ansatz. We emphasize that the new DDF-ROM is built on general ideas of spatial filtering and optimization and is independent of (restrictive) phenomenological arguments. We investigate the DDF-ROM in the numerical simulation of a 2D channel flow past a circular cylinder at Reynolds number $Re=100$. The DDF-ROM is significantly more accurate than the standard projection ROM. Furthermore, the computational costs of the DDF-ROM and the standard projection ROM are similar, both costs being orders of magnitude lower than the computational cost of the full order model. We also compare the new DDF-ROM with modern ROM closure models in the numerical simulation of the 1D Burgers equation. The DDF-ROM is more accurate and significantly more efficient than these ROM closure models.


title: New results on Initial State and Quarkonia with ALICE

abstract: The study of quarkonia in heavy-ion collisions has been the subject of intense experimental and theoretical effort, ever since their production was predicted to be sensitive to the formation of a deconfined state of strongly-interacting matter, known as the Quark-Gluon Plasma (QGP). In p-Pb collisions, Cold Nuclear Matter (CNM) effects, such as nuclear shadowing or partonic energy loss, are expected to influence quarkonium production. The study of such system is therefore crucial to shed light on the mechanisms taking place at the initial-state of quarkonium production, and to disentangle the cold and hot nuclear effects envisioned in Pb-Pb collisions. The ALICE experiment at the LHC, is capable of reconstructing J/$\psi$, $\psi$(2S) and $\Upsilon$ states at forward rapidity through their $\mu^{\rm{+}}\mu^{\rm{-}}$ decay channel, as well as J/$\psi$ at central rapidity through their $e^{\rm{+}}e^{\rm{-}}$ decay channel, down to zero transverse momentum. A review of the main ALICE findings from the measurements of the inclusive quarkonium yields in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV, collected during the LHC Run I period, as well as more recent results from J/$\psi$ measurements in p-Pb at $\sqrt{s_{\rm{NN}}}$ = 8.16 TeV, from LHC Run II period, will be presented in this paper.


title: Laser-induced electronic bridge for characterization of the $~^{229\rm{m}}\rm{Th} \rightarrow ~^{229\rm{g}}\rm{Th}$ nuclear transition with a tunable optical laser

abstract: An alternative method to determine the excitation energy of the $~^{229\rm{m}}\rm{Th}$ isomer via the laser-induced electronic bridge is investigated theoretically. In the presence of an optical or ultra-violet laser at energies that fulfill a two-photon resonance condition, the excited nuclear state can decay by transfering its energy to the electronic shell. A bound electron is then promoted to an excited state by absorption of a laser photon and simultaneous de-excitation of the nucleus. We present calculated rates for the laser-induced electronic bridge process and discuss the experimental requirements for the corresponding setup. Our results show that depending on the actual value of the nuclear transition energy, the rate can be very high, with an enhancement factor compared to the radiative nuclear decay of up to $10^8$.


title: Stability Analysis of Stellar Radiating Filaments

abstract: The aim of this paper is to perform stability analysis of anisotropic dissipative cylindrical collapsing model in $f(R,T,R_{\mu\nu} T^{\mu\nu})$ gravity. In this context, the modified version of hydrodynamical equation is explored by means of dynamical equations and radial perturbation scheme. We examined the role of adiabatic index, dissipation as well as the particular cosmological model on the onset of dynamical instability of the evolving cylindrical system that was initially in hydrostatic equilibrium with Newtonian and post Newtonian approximations. It is pointed out that extra curvature terms of $f(R,T,R_{\mu\nu}T^{\mu\mu\nu})$ gravity tends to increase the stability, while that heat radiations push the system to enter into unstable window. Further, our results reveal the significance of adiabatic index in the stability analysis of cylindrical celestial model.


title: Above-threshold ionization (ATI) in multicenter molecules: the role of the initial state

abstract: A possible route to extract electronic and nuclear dynamics from molecular targets with attosecond temporal and nanometer spatial resolution is to employ recolliding electrons as `probes'. The recollision process in molecules is, however, very challenging to treat using {\it ab initio} approaches. Even for the simplest diatomic systems, such as H$_2$, today's computational capabilities are not enough to give a complete description of the electron and nuclear dynamics initiated by a strong laser field. As a consequence, approximate qualitative descriptions are called to play an important role. In this contribution we extend the work presented in N. Su\'arez {\it et al.}, Phys.~Rev. A {\bf 95}, 033415 (2017), to three-center molecular targets. Additionally, we incorporate a more accurate description of the molecular ground state, employing information extracted from quantum chemistry software packages. This step forward allows us to include, in a detailed way, both the molecular symmetries and nodes present in the high-occupied molecular orbital. We are able to, on the one hand, keep our formulation as analytical as in the case of diatomics, and, on the other hand, to still give a complete description of the underlying physics behind the above-threshold ionization process. The application of our approach to complex multicenter - with more than 3 centers, targets appears to be straightforward.


title: Role of Topologies in the Study of Dynamical Aspects of Stable and Generic Properties of Space-times in General Relativity

abstract: In [1], the authors have studied stability of certain causal properties of space-times in general relativity. As a continuation of this work, in the present paper, we review and discuss, some more aspects of stability which occur in various situations in the dynamics of general relativity. We argue that not only choice of appropriate topology, but also the nature of perturbation, like metric, matter or initial data, are key factors in deciding whether a property is stable or not. We also study certain properties of space-times which are generic in an appropriate mathematical sense. In particular we argue that Linearization stability of a space-time is a generic property.


title: Tight-binding dispersion of the prismatic pentagonal lattice

abstract: Tight-binding Hamiltonian on the prismatic pentagonal lattice is exactly solved to obtain the analytic expressions of dispersion relations and eigenvectors. This lattice is made of prismatic pentagon which is different from Cairo pentagon. Six different dispersion relations and total density of states are obtained. Dispersion relations are symmetric about the zero energy at a particular point in the parameter space. Although a large gap is found for the Cairo pentagonal lattice, no gap as well as no Dirac cone is found to appear in the tight-binding band structure for this prismatic pentagonal lattice. Instead, a pair of van Hove singularities has been identified at two different energy values in the band structure.


title: Decomposition method related to saturated hyperball packings

abstract: In this paper we study the problem of hyperball (hypersphere) packings in $3$-dimensional hyperbolic space. We introduce a new definition of the non-compact saturated ball packings and describe to each saturated hyperball packing, a new procedure to get a decomposition of 3-dimensional hyperbolic space $\HYP$ into truncated tetrahedra. Therefore, in order to get a density upper bound for hyperball packings, it is sufficient to determine the density upper bound of hyperball packings in truncated simplices.


title: Distillation of scalar exchange by coherent hypernucleus production in antiproton-nucleus collisions

abstract: The total and angular differential cross sections of the coherent process $\bar p + {}^AZ \to {}^A_\Lambda(Z-1) + \bar{\Lambda}$ are evaluated at the beam momenta $1.5\div20$ GeV/c within the meson exchange model with bound proton and $\Lambda$-hyperon wave functions. It is shown that the shape of the beam momentum dependence of the hypernucleus production cross sections with various discrete $\Lambda$ states is strongly sensitive to the presence of the scalar $\kappa$-meson exchange in the $\bar p p \to \bar\Lambda \Lambda$ amplitude. This can be used as a clean test of the exchange by scalar $\pi K$ correlation in coherent $\bar p A$ reactions.


title: Exploring High-order three dimensional Virtual Elements: bases and stabilizations

abstract: We present numerical tests of the Virtual Element Method (VEM) tailored for the discretization of a three dimensional Poisson problem with high-order "polynomial" degree (up to $p=10$). Besides, we discuss possible reasons for which the method could return suboptimal-wrong error convergence curves. Among these motivations, we highlight ill-conditioning of the stiffness matrix and not particularly "clever" choices of the stabilizations. We propose variants of the definition of face/bulk degrees of freedom, as well as of stabilizations, which lead to methods that are much more robust in terms of numerical performances.


title: Negative Differential Mobility and Trapping in Active Matter Systems

abstract: Using simulations, we examine the average velocity as a function of applied drift force for active matter particles moving through a random obstacle array. We find that for low drift force, there is an initial flow regime where the mobility increases linearly with drive, while for higher drift forces a regime of negative differential mobility appears in which the velocity decreases with increasing drive due to the trapping of active particles behind obstacles. A fully clogged regime exists at very high drift forces when all the particles are permanently trapped behind obstacles. We find for increasing activity that the overall mobility is nonmonotonic, with an enhancement of the mobility for small levels of activity and a decrease in mobility for large activity levels. We show how these effects evolve as a function of disk and obstacle density, active run length, drift force, and motor force.


title: Families of invariant tori in KAM theory: interplay of integer characteristics

abstract: The purpose of this brief note is twofold. First, we summarize in a very concise form the principal information on Whitney smooth families of quasi-periodic invariant tori in various contexts of KAM theory. Our second goal is to attract (via an informal discussion and a simple example) the experts' attention to the peculiarities of the so-called excitation of elliptic normal modes in the reversible context 2.


title: Impact of User Height on the Coverage of 3D Beamforming-Enabled Massive MIMO Systems

abstract: In this paper, we perform a coverage analysis of a cellular massive multiple-input multiple-output (MaMIMO) network which adopts 3D beamforming. In contrast to the previous works on 3D beamforming which assume that all users are placed on the ground (i.e, in a 2D environment), we consider a more practical scenario where the users of the network are dropped in a 3D environment with different heights. In this scenario, by adopting a stochastic geometry framework, first we calculate the coverage probability of the network as a function of the users height. Then considering this coverage probability as the objective function, the optimum tilt angle of the base station antenna pattern is found and the effect of users' heights on this tilt angle is investigated. Our numerical results show that by taking the users' heights into account and using the proposed method, a considerable improvement is achieved in the performance of the network compared to other similar 3D beamforming methods that ignore the users'height distribution.


title: The role of tumbling frequency and persistence in optimal run-and-tumble chemotaxis

abstract: One of simplest examples of navigation found in nature is run-and-tumble chemotaxis. Tumbles reorient cells randomly, and cells can drift toward attractants or away from repellents by biasing the frequency of these events. The post-tumble swimming directions are typically correlated with those prior, as measured by the variance of the reorientation angle distribution. This variance can range from large, in the case of bacteria, to so small that tumble events are imperceptible, as observed in choanoflagellates. This raises the question of optimality: why is such a range of persistence observed in nature? Here we study persistent run-and-tumble dynamics, focusing first on the optimisation of the linearised chemotactic response within the two-dimensional parameter space of tumble frequency and angular persistence. Although an optimal persistence does exist for a given tumble frequency, in the full parameter space there is a continuum of optimal solutions. Introducing finite tumble times that depend on the persistence can change this picture, illuminating one possible method for selecting tumble persistence based on species-specific reorientation dynamics. Moving beyond linear theory we find that optimal chemotactic strengths exist, and that these maximise reaction when swimming in a wrong direction, but have little or no reaction when swimming with even the slightest projection along the chemoattractant gradient.


title: Lasserre hierarchy for large scale polynomial optimization in real and complex variables

abstract: We propose general notions to deal with large scale polynomial optimization problems and demonstrate their efficiency on a key industrial problem of the twenty first century, namely the optimal power flow problem. These notions enable us to find global minimizers on instances with up to 4,500 variables and 14,500 constraints. First, we generalize the Lasserre hierarchy from real to complex to numbers in order to enhance its tractability when dealing with complex polynomial optimization. Complex numbers are typically used to represent oscillatory phenomena, which are omnipresent in physical systems. Using the notion of hyponormality in operator theory, we provide a finite convergence criterion which generalizes the Curto-Fialkow conditions of the real Lasserre hierarchy. Second, we introduce the multi-ordered Lasserre hierarchy in order to exploit sparsity in polynomial optimization problems (in real or complex variables) while preserving global convergence. It is based on two ideas: 1) to use a different relaxation order for each constraint, and 2) to iteratively seek a closest measure to the truncated moment data until a measure matches the truncated data. Third and last, we exhibit a block diagonal structure of the Lasserre hierarchy in the presence of commonly encountered symmetries. To the best of our knowledge, the Lasserre hierarchy was previously limited to small scale problems, while we solve a large scale industrial problem with thousands of variables and constraints to global optimality.


title: ProSLAM: Graph SLAM from a Programmer's Perspective

abstract: In this paper we present ProSLAM, a lightweight stereo visual SLAM system designed with simplicity in mind. Our work stems from the experience gathered by the authors while teaching SLAM to students and aims at providing a highly modular system that can be easily implemented and understood. Rather than focusing on the well known mathematical aspects of Stereo Visual SLAM, in this work we highlight the data structures and the algorithmic aspects that one needs to tackle during the design of such a system. We implemented ProSLAM using the C++ programming language in combination with a minimal set of well known used external libraries. In addition to an open source implementation, we provide several code snippets that address the core aspects of our approach directly in this paper. The results of a thorough validation performed on standard benchmark datasets show that our approach achieves accuracy comparable to state of the art methods, while requiring substantially less computational resources.


title: Random cover times using the Poisson cylinder process

abstract: In this paper we deal with the classical problem of random cover times. We investigate the distribution of the time it takes for a Poisson process of cylinders to cover a set $A \subset \mathbb{R}^d.$ This Poisson process of cylinders is invariant under rotations, reflections and translations, and in addition we add a time component so that cylinders are "raining from the sky" at unit rate. Our main results concerns the asymptotic of this cover time as the set $A$ grows. If the set $A$ is discrete and well separated, we show convergence of the cover time to a Gumbel distribution. If instead $A$ has positive box dimension (and satisfies a weak additional assumption), we find the correct rate of convergence.


title: Layered Space-Time Index Coding

abstract: Multicasting $K$ independent messages via multiple-input multiple-output (MIMO) channels to multiple users where each user already has a subset of messages as side information is studied. A general framework of constructing layered space-time index coding (LSTIC) from a large class of space-time block codes (STBC), including perfect STBC, is proposed. We analyze the proposed LSTIC and show that it provides minimum determinant gains that are exponential with the amount of information contained in the side information for any possible side information. When constructed over a perfect STBC, the proposed LSTIC is itself a perfect STBC and hence many desired properties are preserved. To illustrate, we construct LSTIC over the following well-known STBCs: Golden code; $3\times 3$, $4\times 4$, $6\times 6$ perfect STBCs; and Alamouti code. Simulation results show that the obtained side information gain can be well predicted by our analysis.


title: Neural Network Based Nonlinear Weighted Finite Automata

abstract: Weighted finite automata (WFA) can expressively model functions defined over strings but are inherently linear models. Given the recent successes of nonlinear models in machine learning, it is natural to wonder whether ex-tending WFA to the nonlinear setting would be beneficial. In this paper, we propose a novel model of neural network based nonlinearWFA model (NL-WFA) along with a learning algorithm. Our learning algorithm is inspired by the spectral learning algorithm for WFAand relies on a nonlinear decomposition of the so-called Hankel matrix, by means of an auto-encoder network. The expressive power of NL-WFA and the proposed learning algorithm are assessed on both synthetic and real-world data, showing that NL-WFA can lead to smaller model sizes and infer complex grammatical structures from data.


title: Search for the SM Higgs boson produced in association with top quarks at $\sqrt{s}$ = 13 TeV with the ATLAS detector at the LHC

abstract: A review of the searches for the SM Higgs boson produced in association with a pair of top quarks, ${t\bar{t}H}$, using up to 13.3 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector at the LHC is presented. Searches in the diphoton, multilepton and $b\bar{b}$ decay channels are summarised. The combination of these searches yields a measured (expected) significance for the observation of the ${t\bar{t}H}$ production process of $2.8 \sigma$ ($1.8 \sigma$).


title: On the decidability of the existence of polyhedral invariants in transition systems

abstract: Automated program verification often proceeds by exhibiting inductive invariants entailing the desired properties.For numerical properties, a classical class of invariants is convex polyhedra: solution sets of system of linear (in)equalities.Forty years of research on convex polyhedral invariants have focused, on the one hand, on identifying "easier" subclasses, on the other hand on heuristics for finding general convex polyhedra.These heuristics are however not guaranteed to find polyhedral inductive invariants when they exist.To our best knowledge, the existence of polyhedral inductive invariants has never been proved to be undecidable.In this article, we show that the existence of convex polyhedral invariants is undecidable, even if there is only one control state in addition to the "bad" one.The question is still open if one is not allowed any nonlinear constraint.


title: Solving nonlinear optimal control problems with state and control delays by shooting methods combined with numerical continuation on the delays

abstract: In this paper we introduce a new procedure to solve nonlinear optimal control problems with delays which exploits indirect methods combined with numerical homotopy procedures. It is known that solving this kind of problems via indirect methods (which arise from the Pontrya-gin Maximum Principle) is complex and computationally demanding because their implementation is faced to two main difficulties: the extremal equations involve forward and backward terms, and besides, the related shooting method has to be carefully initialized. Here, starting from the solution of the non-delayed version of the optimal control problem, delays are introduced by a numerical continuation. This creates a sequence of optimal delayed solutions that converges to the desired solution. We establish a convergence theorem ensuring the continuous dependence w.r.t. the delay of the optimal state, of the optimal control (in a weak sense) and of the corresponding adjoint vector. The convergence of the adjoint vector represents the most challenging step to prove and it is crucial for the well-posedness of the proposed homotopy procedure. Two numerical examples are proposed and analyzed to show the efficiency of this approach.


title: Generating Music Medleys via Playing Music Puzzle Games

abstract: Generating music medleys is about finding an optimal permutation of a given set of music clips. Toward this goal, we propose a self-supervised learning task, called the music puzzle game, to train neural network models to learn the sequential patterns in music. In essence, such a game requires machines to correctly sort a few multisecond music fragments. In the training stage, we learn the model by sampling multiple non-overlapping fragment pairs from the same songs and seeking to predict whether a given pair is consecutive and is in the correct chronological order. For testing, we design a number of puzzle games with different difficulty levels, the most difficult one being music medley, which requiring sorting fragments from different songs. On the basis of state-of-the-art Siamese convolutional network, we propose an improved architecture that learns to embed frame-level similarity scores computed from the input fragment pairs to a common space, where fragment pairs in the correct order can be more easily identified. Our result shows that the resulting model, dubbed as the similarity embedding network (SEN), performs better than competing models across different games, including music jigsaw puzzle, music sequencing, and music medley. Example results can be found at our project website, https://remyhuang.github.io/DJnet.


title: The Power of Synchronisation: Formal Analysis of Power Consumption in Networks of Pulse-Coupled Oscillators

abstract: We assess the power consumption of network synchronisation protocols, particularly the energy required to synchronise all nodes across a network. We use the widely adopted approach of bio-inspired, pulse-coupled oscillators to achieve network-wide synchronisation and provide an extended formal model of just such a protocol, enhanced with structures for recording energy usage. Exhaustive analysis is then carried out through formal verification, utilising the PRISM model checker to calculate the resources consumed on each possible system execution. This allows us to assess a range of parameter instantiations and to explore trade-offs between power consumption and time to synchronise. This provides a principled basis for the formal analysis of a much broader range of large-scale network protocols.


title: Radio emission from satellite-Jupiter interactions (especially Ganymede)

abstract: Analyzing a database of 26 years of observations of Jupiter from the Nan\c{c}ay Decameter Array, we study the occurrence of Io-independent emissions as a function of the orbital phase of the other Galilean satellites and Amalthea. We identify unambiguously the emissions induced by Ganymede and characterize their intervals of occurrence in CML and Ganymede phase and longitude. We also find hints of emissions induced by Europa and, surprisingly, by Amalthea. The signature of Callisto-induced emissions is more tenuous.


title: Welfare effects of information and rationality in portfolio decisions under parameter uncertainty

abstract: We analyze and quantify, in a financial market with parameter uncertainty and for a Constant Relative Risk Aversion investor, the utility effects of two different boundedly rational (i.e., sub-optimal) investment strategies (namely, myopic and unconditional strategies) and compare them between each other and with the utility effect of full information. We show that effects are mainly caused by full information and predictability, being the effect of learning marginal. We also investigate the saver's decision of whether to manage her/his portfolio personally (DIY investor) or hire, against the payment of a management fee, a professional investor and find that delegation is mainly motivated by the belief that professional advisors are, depending on investment horizon and risk aversion, either better informed ("insiders") or more capable of gathering and processing information rather than their ability of learning from financial data. In particular, for very short investment horizons, delegation is primarily, if not exclusively, motivated by the beliefs that professional investors are better informed.


title: Modified Dark Matter: Relating Dark Energy, Dark Matter and Baryonic Matter

abstract: Modified dark matter (MDM) is a phenomenological model of dark matter, inspired by gravitational thermodynamics. For an accelerating Universe with positive cosmological constant ($\Lambda$), such phenomenological considerations lead to the emergence of a critical acceleration parameter related to $\Lambda$. Such a critical acceleration is an effective phenomenological manifestation of MDM, and it is found in correlations between dark matter and baryonic matter in galaxy rotation curves. The resulting MDM mass profiles, which are sensitive to $\Lambda$, are consistent with observational data at both the galactic and cluster scales. In particular, the same critical acceleration appears both in the galactic and cluster data fits based on MDM. Furthermore, using some robust qualitative arguments, MDM appears to work well on cosmological scales, even though quantitative studies are still lacking. Finally, we comment on certain non-local aspects of the quanta of modified dark matter, which may lead to novel non-particle phenomenology and which may explain why, so far, dark matter detection experiments have failed to detect dark matter particles.


title: Scalable and Efficient Statistical Inference with Estimating Functions in the MapReduce Paradigm for Big Data

abstract: The theory of statistical inference along with the strategy of divide-and-conquer for large- scale data analysis has recently attracted considerable interest due to great popularity of the MapReduce programming paradigm in the Apache Hadoop software framework. The central analytic task in the development of statistical inference in the MapReduce paradigm pertains to the method of combining results yielded from separately mapped data batches. One seminal solution based on the confidence distribution has recently been established in the setting of maximum likelihood estimation in the literature. This paper concerns a more general inferential methodology based on estimating functions, termed as the Rao-type confidence distribution, of which the maximum likelihood is a special case. This generalization provides a unified framework of statistical inference that allows regression analyses of massive data sets of important types in a parallel and scalable fashion via a distributed file system, including longitudinal data analysis, survival data analysis, and quantile regression, which cannot be handled using the maximum likelihood method. This paper investigates four important properties of the proposed method: computational scalability, statistical optimality, methodological generality, and operational robustness. In particular, the proposed method is shown to be closely connected to Hansen's generalized method of moments (GMM) and Crowder's optimality. An interesting theoretical finding is that the asymptotic efficiency of the proposed Rao-type confidence distribution estimator is always greater or equal to the estimator obtained by processing the full data once. All these properties of the proposed method are illustrated via numerical examples in both simulation studies and real-world data analyses.


title: FESTUNG: A MATLAB / GNU Octave toolbox for the discontinuous Galerkin method. Part III: Hybridized discontinuous Galerkin (HDG) formulation

abstract: The third paper in our series on open source MATLAB / GNU Octave implementation of the discontinuous Galerkin (DG) method(s) focuses on a hybridized formulation. The main aim of this ongoing work is to develop rapid prototyping techniques covering a range of standard DG methodologies and suitable for small to medium sized applications. Our FESTUNG package relies on fully vectorized matrix / vector operations throughout, and all details of the implementation are fully documented. Once again, great care was taken to maintain a direct mapping between discretization terms and code routines as well as to ensure full compatibility to GNU Octave. The current work formulates a hybridized DG scheme for linear advection problem, describes hybrid approximation spaces on the mesh skeleton, and compares the performance of this discretization to the standard (element-based) DG method for different polynomial orders.


title: Microhexcavity Plasma Panel Detectors

abstract: Plasma panel detectors are a variant of micropattern detectors that are sensitive to ionizing radiation. They are motivated by the design and operation of plasma display panels. The detectors consist of arrays of electrically and optically isolated pixels defined by metallized cavities embedded in a dielectric substrate. These are hermetically sealed gaseous detectors that use exclusively non-hydrocarbon gas mixtures. The newest variant of these closed-architecture detectors is known as the Microhexcavity plasma panel detector ($\mu$Hex) consisting of 2 mm wide, regular close-packed hexagonal pixels each with a circular thick-film anode. The fabrication, staging, and operation of these detectors is described. Initial tests with the $\mu$Hex detectors operated in Geiger mode yield Volt-level signals in the presence of ionizing radiation. The spontaneous discharge rate in the absence of a source is roughly 3-4 orders of magnitude lower compared to the rates measured using low energy betas.


title: Black holes in vector-tensor theories and their thermodynamics

abstract: In this paper, we study Einstein gravity either minimally or non-minimally coupled to a vector field which breaks the gauge symmetry explicitly in general dimensions. We first consider a minimal theory which is simply the Einstein-Proca theory extended with a quartic self-interaction term for the vector field. We obtain its general static maximally symmetric black hole solution and study the thermodynamics using Wald formalism. The aspects of the solution are much like a Reissner-Nordstr{\o}m black hole in spite of that a global charge cannot be defined for the vector. For non-minimal theories, we obtain a lot of exact black hole solutions, depending on the parameters of the theories. In particular, many of the solutions are general static and have maximal symmetry. However, there are some subtleties and ambiguities in the derivation of the first laws because the existence of an algebraic degree of freedom of the vector in general invalids the Wald entropy formula. The thermodynamics of these solutions deserves further studies.


title: An Efficient Evolutionary Based Method For Image Segmentation

abstract: The goal of this paper is to present a new efficient image segmentation method based on evolutionary computation which is a model inspired from human behavior. Based on this model, a four layer process for image segmentation is proposed using the split/merge approach. In the first layer, an image is split into numerous regions using the watershed algorithm. In the second layer, a co-evolutionary process is applied to form centers of finals segments by merging similar primary regions. In the third layer, a meta-heuristic process uses two operators to connect the residual regions to their corresponding determined centers. In the final layer, an evolutionary algorithm is used to combine the resulted similar and neighbor regions. Different layers of the algorithm are totally independent, therefore for certain applications a specific layer can be changed without constraint of changing other layers. Some properties of this algorithm like the flexibility of its method, the ability to use different feature vectors for segmentation (grayscale, color, texture, etc), the ability to control uniformity and the number of final segments using free parameters and also maintaining small regions, makes it possible to apply the algorithm to different applications. Moreover, the independence of each region from other regions in the second layer, and the independence of centers in the third layer, makes parallel implementation possible. As a result the algorithm speed will increase. The presented algorithm was tested on a standard dataset (BSDS 300) of images, and the region boundaries were compared with different people segmentation contours. Results show the efficiency of the algorithm and its improvement to similar methods. As an instance, in 70% of tested images, results are better than ACT algorithm, besides in 100% of tested images, we had better results in comparison with VSP algorithm.


title: Monte Carlo methods for massively parallel computers

abstract: Applications that require substantial computational resources today cannot avoid the use of heavily parallel machines. Embracing the opportunities of parallel computing and especially the possibilities provided by a new generation of massively parallel accelerator devices such as GPUs, Intel's Xeon Phi or even FPGAs enables applications and studies that are inaccessible to serial programs. Here we outline the opportunities and challenges of massively parallel computing for Monte Carlo simulations in statistical physics, with a focus on the simulation of systems exhibiting phase transitions and critical phenomena. This covers a range of canonical ensemble Markov chain techniques as well as generalized ensembles such as multicanonical simulations and population annealing. While the examples discussed are for simulations of spin systems, many of the methods are more general and moderate modifications allow them to be applied to other lattice and off-lattice problems including polymers and particle systems. We discuss important algorithmic requirements for such highly parallel simulations, such as the challenges of random-number generation for such cases, and outline a number of general design principles for parallel Monte Carlo codes to perform well.

categories: physics.comp-ph cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech


title: Tight Semi-Nonnegative Matrix Factorization

abstract: The nonnegative matrix factorization is a widely used, flexible matrix decomposition, finding applications in biology, image and signal processing and information retrieval, among other areas. Here we present a related matrix factorization. A multi-objective optimization problem finds conical combinations of templates that approximate a given data matrix. The templates are chosen so that as far as possible only the initial data set can be represented this way. However, the templates are not required to be nonnegative nor convex combinations of the original data.


title: A Tutorial on Deep Learning for Music Information Retrieval

abstract: Following their success in Computer Vision and other areas, deep learning techniques have recently become widely adopted in Music Information Retrieval (MIR) research. However, the majority of works aim to adopt and assess methods that have been shown to be effective in other domains, while there is still a great need for more original research focusing on music primarily and utilising musical knowledge and insight. The goal of this paper is to boost the interest of beginners by providing a comprehensive tutorial and reducing the barriers to entry into deep learning for MIR. We lay out the basic principles and review prominent works in this hard to navigate the field. We then outline the network structures that have been successful in MIR problems and facilitate the selection of building blocks for the problems at hand. Finally, guidelines for new tasks and some advanced topics in deep learning are discussed to stimulate new research in this fascinating field.


title: Quantum Hall States and Conformal Field Theory on a Singular Surface

abstract: In [Can et al. 2016], quantum Hall states on singular surfaces were shown to possess an emergent conformal symmetry. In this paper, we develop this idea further and flesh out details on the emergent conformal symmetry in holomorphic adiabatic states, which we define in the paper. We highlight the connection between the universal features of geometric transport of quantum Hall states and holomorphic dimension of primary fields in conformal field theory. In parallel we compute the universal finite-size corrections to the free energy of a critical system on a hyperbolic sphere with conical and cusp singularities, thus extending the result of Cardy and Peschel for critical systems on a flat cone [Cardy, Peschel 1988], and the known results for critical systems on polyhedra and flat branched Riemann surfaces.


title: Identifiability of dynamical networks: which nodes need be measured?

abstract: Much recent research has dealt with the identifiability of a dynamical network in which the node signals are connected by causal linear time-invariant transfer functions and are possibly excited by known external excitation signals and/or unknown noise signals. So far all results on the identifiability of the whole network have assumed that all node signals are measured. Under this assumption, it has been shown that such networks are identifiable only if some prior knowledge is available about the structure of the network, in particular the structure of the excitation. In this paper we present the first results for the situation where not all node signals are measurable, under the assumptions that the topology of the network is known, that each node is excited by a known signal and that the nodes are noise-free. Using graph theoretical properties, we show that the transfer functions that can be identified depend essentially on the topology of the paths linking the corresponding vertices to the measured nodes. An important outcome of our research is that, under those assumptions, a network can often be identified using only a small subset of node measurements.


title: Elastic bending energy: a variational approach

abstract: Geometric continuum models for fluid lipid membranes are considered using classical field theory, within a covariant variational approach. The approach is cast as a higher-derivative Lagrangian formulation of continuum classical field theory, and it can be seen as a covariant version of the field theoretical variational approach that uses the height representation. This novel Lagrangian formulation is presented first for a generic reparametrization invariant geometric model, deriving its equilibrium condition equation, or shape equation, and its linear and angular stress tensors, using the classical Canham-Helfrich elastic bending energy for illustration. The robustness of the formulation is established by extending it to the presence of external forces, and to the case of heterogenous lipid membranes, breaking reparametrization invariance. In addition, a useful and compact general expression for the second variation of the free energy is obtained within the Lagrangian formulation, as a first step towards the study of the stability of membrane configurations. The simple structure of the expressions derived for the basic entities that appear in the mechanics of a lipid membrane is a direct consequence of the well established power of a Lagrangian variational approach. The paper is self-contained, and it is meant to provide, besides a new framework, also a convenient introduction to the mechanics of lipid membranes.


title: Improved thermal lattice Boltzmann model for simulation of liquid-vapor phase change

abstract: In this paper, an improved thermal lattice Boltzmann (LB) model is proposed for simulating liquid-vapor phase change, which is aimed at improving an existing thermal LB model for liquid-vapor phase change [S. Gong and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012)]. First, we emphasize that the replacement of \[{\left( {\rho {c_V}} \right)^{ - 1}}\nabla \cdot \left( {\lambda \nabla T} \right)\] with \[\nabla \cdot \left( {\chi \nabla T} \right)\] is an inappropriate treatment for diffuse interface modeling of liquid-vapor phase change. Furthermore, the error terms ${\partial_{t0}}\left( {Tv} \right) + \nabla \cdot \left( {Tvv} \right)$, which exist in the macroscopic temperature equation recovered from the standard thermal LB equation, are eliminated in the present model through a way that is consistent with the philosophy of the LB method. In addition, the discrete effect of the source term is also eliminated in the present model. Numerical simulations are performed for droplet evaporation and bubble nucleation to validate the capability of the improved model for simulating liquid-vapor phase change. Numerical comparisons show that the aforementioned replacement leads to significant numerical errors and the error terms in the recovered macroscopic temperature equation also result in considerable errors.


title: Life-span of blowup solutions to semilinear wave equation with space-dependent critical damping

abstract: This paper is concerned with the blowup phenomena for initial value problem of semilinear wave equation with critical space-dependent damping term (DW:$V$). The main result of the present paper is to give a solution of the problem and to provide a sharp estimate for lifespan for such a solution when $\frac{N}{N-1}<p\leq p_S(N+V_0)$, where $p_S(N)$ is the Strauss exponent for (DW:$0$). The main idea of the proof is due to the technique of test functions for (DW:$0$) originated by Zhou--Han (2014, MR3169791). Moreover, we find a new threshold value $V_0=\frac{(N-1)^2}{N+1}$ for the coefficient of critical and singular damping $|x|^{-1}$.


title: On Early-stage Debunking Rumors on Twitter: Leveraging the Wisdom of Weak Learners

abstract: Recently a lot of progress has been made in rumor modeling and rumor detection for micro-blogging streams. However, existing automated methods do not perform very well for early rumor detection, which is crucial in many settings, e.g., in crisis situations. One reason for this is that aggregated rumor features such as propagation features, which work well on the long run, are - due to their accumulating characteristic - not very helpful in the early phase of a rumor. In this work, we present an approach for early rumor detection, which leverages Convolutional Neural Networks for learning the hidden representations of individual rumor-related tweets to gain insights on the credibility of each tweets. We then aggregate the predictions from the very beginning of a rumor to obtain the overall event credits (so-called wisdom), and finally combine it with a time series based rumor classification model. Our extensive experiments show a clearly improved classification performance within the critical very first hours of a rumor. For a better understanding, we also conduct an extensive feature evaluation that emphasized on the early stage and shows that the low-level credibility has best predictability at all phases of the rumor lifetime.


title: Explicit Commutativity Conditions for Second-order Linear Time-Varying Systems with Non-Zero Initial Conditions

abstract: Although the explicit commutativitiy conditions for second-order linear time-varying systems have been appeared in some literature, these are all for initially relaxed systems. This paper presents explicit necessary and sufficient commutativity conditions for commutativity of second-order linear time-varying systems with non-zero initial conditions. It has appeared interesting that the second requirement for the commutativity of non-relaxed systems plays an important role on the commutativity conditions when non-zero initial conditions exist. Another highlight is that the commutativity of switched systems is considered and spoiling of commutativity at the switching instants is illustrated for the first time. The simulation results support the theory developed in the paper.


title: Maximum matchings and minimum dominating sets in Apollonian networks and extended Tower of Hanoi graphs

abstract: The Apollonian networks display the remarkable power-law and small-world properties as observed in most realistic networked systems. Their dual graphs are extended Tower of Hanoi graphs, which are obtained from the Tower of Hanoi graphs by adding a special vertex linked to all its three extreme vertices. In this paper, we study analytically maximum matchings and minimum dominating sets in Apollonian networks and their dual graph- s, both of which have found vast applications in various fields, e.g. structural controllability of complex networks. For both networks, we determine their matching number, domination number, the number of maximum matchings, as well as the number of minimum dominating sets.


title: Commutativity of Systems with their Feedback Conjugates

abstract: After introducing commutativity concept and summarizing the relevant literature, this work is focused on the commutativity of feedback conjugates. It is already known that a linear time-varying differential system describing a single input-single output dynamical system is always commutative with its constant gain feedback pairs. In this article, it is proven that among the time-varying feedback conjugates of a linear time-varying system, constant feedback conjugates are the only commutative feedback pairs and any of the time-varying feedback conjugates cannot constitutes a commutative pair of a linear time-varying system.


title: Life-span of solutions to semilinear wave equation with time-dependent critical damping for specially localized initial data

abstract: This paper is concerned with the blowup phenomena for initial value problem of semilinear wave equation with critical time-dependent damping term (DW). The result is the sharp upper bound of lifespan of solution with respect to the small parameter $\ep$ when $p_F(N)\leq p\leq p_0(N+\mu)$, where $p_F(N)$ denotes the Fujita exponent for the nonlinear heat equations and $p_0(n)$ denotes the Strauss exponent for nonlinear wave equation in $n$-dimension with $\mu=0$. Consequently, by connecting the result of D'Abbicco--Lucente--Reissig 2015, our result clarifies the threshold exponent $p_0(N+\mu)$ for dividing blowup phenomena and global existence of small solutions when $N=3$. The crucial idea is to construct suitable test functions satisfying the conjugate linear equation $\pa_t^2\Phi-\Delta \Phi-\pa_t(\frac{\mu}{1+t}\Phi)=0$ of (DW) including the Gauss hypergeometric functions; note that the construction of test functions is different from Zhou--Han in 2014.


title: An Inversion-Based Learning Approach for Improving Impromptu Trajectory Tracking of Robots with Non-Minimum Phase Dynamics

abstract: This paper presents a learning-based approach for impromptu trajectory tracking for non-minimum phase systems, i.e., systems with unstable inverse dynamics. Inversion-based feedforward approaches are commonly used for improving tracking performance; however, these approaches are not directly applicable to non-minimum phase systems due to their inherent instability. In order to resolve the instability issue, existing methods have assumed that the system model is known and used pre-actuation or inverse approximation techniques. In this work, we propose an approach for learning a stable, approximate inverse of a non-minimum phase baseline system directly from its input-output data. Through theoretical discussions, simulations, and experiments on two different platforms, we show the stability of our proposed approach and its effectiveness for high-accuracy, impromptu tracking. Our approach also shows that including more information in the training, as is commonly assumed to be useful, does not lead to better performance but may trigger instability and impact the effectiveness of the overall approach.


title: High-resolution reconstruction of cellular traction-force distributions: the role of physically motivated constraints and compressive regularization

abstract: We develop a method to reconstruct, from measured displacements of an underlying elastic substrate, the spatially dependent forces that cells or tissues impart on it. Given newly available high-resolution images of substrate displacements, it is desirable to be able to reconstruct small scale, compactly supported focal adhesions which are often localized and exist only within the footprint of a cell. In addition to the standard quadratic data mismatch terms that define least-squares fitting, we motivate a regularization term in the objective function that penalizes vectorial invariants of the reconstructed surface stress while preserving boundaries. We solve this inverse problem by providing a numerical method for setting up a discretized inverse problem that is solvable by standard convex optimization techniques. By minimizing the objective function subject to a number of important physically motivated constraints, we are able to efficiently reconstruct stress fields with localized structure from simulated and experimental substrate displacements. Our method incorporates the exact solution for the stress tensor accurate to first-order finite-differences and motivates the use of distance-based cut-offs for data inclusion and problem sparsification.


title: A Review of Evaluation Techniques for Social Dialogue Systems

abstract: In contrast with goal-oriented dialogue, social dialogue has no clear measure of task success. Consequently, evaluation of these systems is notoriously hard. In this paper, we review current evaluation methods, focusing on automatic metrics. We conclude that turn-based metrics often ignore the context and do not account for the fact that several replies are valid, while end-of-dialogue rewards are mainly hand-crafted. Both lack grounding in human perceptions.


title: Commutativity of First-order Discrete-Time Linear Time-varying Systems

abstract: After introducing the concept of commutativity for continuous-time linear time-varying systems, the related literature and the results obtained so far are presented. For a simple introduction of the commutativity of discrete-time linear time-varying systems, the problem is formulated for first-order systems. Finally, explicit necessary and sufficient conditions for the commutativity of first-order discrete-time linear time-varying systems are derived, and their advantageous use in digital system design is illustrated; which are the main objectives of the paper. The results are verified by examples which include an application in amplitude modulation for digital telecommunication.


title: Exploiting skeletal structure in computer vision annotation with Benders decomposition

abstract: Many annotation problems in computer vision can be phrased as integer linear programs (ILPs). The use of standard industrial solvers does not to exploit the underlying structure of such problems eg, the skeleton in pose estimation. The leveraging of the underlying structure in conjunction with industrial solvers promises increases in both speed and accuracy. Such structure can be exploited using Bender's decomposition, a technique from operations research, that solves complex ILPs or mixed integer linear programs by decomposing them into sub-problems that communicate via a master problem. The intuition is that conditioned on a small subset of the variables the solution to the remaining variables can be computed easily by taking advantage of properties of the ILP constraint matrix such as block structure. In this paper we apply Benders decomposition to a typical problem in computer vision where we have many sub-ILPs (eg, partitioning of detections, body-parts) coupled to a master ILP (eg, constructing skeletons). Dividing inference problems into a master problem and sub-problems motivates the development of a plethora of novel models, and inference approaches for the field of computer vision.


title: The Merging Path Plot: adaptive fusing of k-groups with likelihood-based model selection

abstract: There are many statistical tests that verify the null hypothesis: the variable of interest has the same distribution among k-groups. But once the null hypothesis is rejected, how to present the structure of dissimilarity between groups? In this article, we introduce The Merging Path Plot - a methodology, and factorMerger - an R package, for exploration and visualization of k-group dissimilarities. Comparison of k-groups is one of the most important issues in exploratory analyses and it has zillions of applications. The classical solution is to test a~null hypothesis that observations from all groups come from the same distribution. If the global null hypothesis is rejected, a~more detailed analysis of differences among pairs of groups is performed. The traditional approach is to use pairwise post hoc tests in order to verify which groups differ significantly. However, this approach fails with a large number of groups in both interpretation and visualization layer. The~Merging Path Plot methodology solves this problem by using an easy-to-understand description of dissimilarity among groups based on Likelihood Ratio Test (LRT) statistic.


title: Minor Contribution of Quasars to Ionizing Photon Budget at z~6: Update on Quasar Luminosity Function at the Faint-end with Subaru/Suprime-Cam

abstract: We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.


title: Controllability and lack of controllability with smooth controls in viscoelasticity via moment methods

abstract: In this paper we study controllability of a linear equation with persistent memory when the control belongs to $ H^k_0(0,T;L^2(\ZOMq)) $. In the case the memory is zero, our equation is reduced to the wave equation and a result due to Everdoza and Zuazua informally states that smoother targets can be reached by using smoother controls. In this paper we prove that this result can be partially extended to systems with memory, but that the memory is an obstruction to a complete extensions.


title: Risk-Aware Multi-Armed Bandit Problem with Application to Portfolio Selection

abstract: Sequential portfolio selection has attracted increasing interests in the machine learning and quantitative finance communities in recent years. As a mathematical framework for reinforcement learning policies, the stochastic multi-armed bandit problem addresses the primary difficulty in sequential decision making under uncertainty, namely the exploration versus exploitation dilemma, and therefore provides a natural connection to portfolio selection. In this paper, we incorporate risk-awareness into the classic multi-armed bandit setting and introduce an algorithm to construct portfolio. Through filtering assets based on the topological structure of financial market and combining the optimal multi-armed bandit policy with the minimization of a coherent risk measure, we achieve a balance between risk and return.


title: Registering the evolutionary history in individual-based models of speciation

abstract: Understanding the emergence of biodiversity patterns in nature is a central problem in biology. Theoretical models of speciation have addressed this question in the macroecological scale, but little has been investigated in the macroevolutionary context. Knowledge of the evolutionary history allows the study of patterns underlying the processes considered in these models, revealing their signatures and the role of speciation and extinction in shaping macroevolutionary patterns. In this paper we introduce two algorithms to record the evolutionary history of populations in individual-based models of speciation, from which genealogies and phylogenies can be constructed. The first algorithm relies on saving ancestral-descendant relationships, generating a matrix that contains the times to the most recent common ancestor between all pairs of individuals at every generation (the Most Recent Common Ancestor Time matrix, MRCAT). The second algorithm directly records all speciation and extinction events throughout the evolutionary process, generating a matrix with the true phylogeny of species (the Sequential Speciation and Extinction Events, SSEE). We illustrate the use of these algorithms in a spatially explicit individual-based model of speciation. We compare the trees generated via MRCAT and SSEE algorithms with trees inferred by methods that use only genetic distance among extant species, commonly used in empirical studies and applied here to simulated genetic data. Comparisons between tress are performed with metrics describing the overall topology, branch length distribution and imbalance of trees. We observe that both MRCAT and distance-based trees differ from the true phylogeny, with the first being closer to the true tree than the second.


title: Investigating Storage as a Service Cloud Platform: pCloud as a Case Study

abstract: Due to the flexibility, affordability and portability of cloud storage, individuals and companies envisage the cloud storage as one of the preferred storage media nowadays. This attracts the eyes of cyber criminals, since much valuable informa- tion such as user credentials, and private customer records are stored in the cloud. There are many ways for criminals to compromise cloud services; ranging from non-technical attack methods, such as social engineering, to deploying advanced malwares. Therefore, it is vital for cyber forensics examiners to be equipped and informed about best methods for investigation of different cloud platforms. In this chapter, using pCloud (an extensively used online cloud storage service) as a case study, and we elaborate on different kinds of artefacts retrievable during a forensics examination. We carried out our experiments on four different virtual machines running four popular operating systems: a 64 bit Windows 8, Ubuntu 14.04.1 LTS, Android 4.4.2, and iOS 8.1. Moreover, we examined cloud remnants of two different web browsers: Internet Explorer and Google Chrome on Windows. We believe that our study would promote awareness among digital forensic examiners on how to conduct cloud storage forensics examination.


title: Power in High-Dimensional Testing Problems

abstract: Fan et al. (2015) recently introduced a remarkable method for increasing asymptotic power of tests in high-dimensional testing problems. If applicable to a given test, their power enhancement principle leads to an improved test that has the same asymptotic size, uniformly non-inferior asymptotic power, and is consistent against a strictly broader range of alternatives than the initially given test. We study under which conditions this method can be applied and show the following: In asymptotic regimes where the dimensionality of the parameter space is fixed as sample size increases, there often exist tests that can not be further improved with the power enhancement principle. However, when the dimensionality of the parameter space increases sufficiently slowly with sample size and a marginal local asymptotic normality (LAN) condition is satisfied, every test with asymptotic size smaller than one can be improved with the power enhancement principle. While the marginal LAN condition alone does not allow one to extend the latter statement to all rates at which the dimensionality increases with sample size, we give sufficient conditions under which this is the case.


title: Likelihood Approximation With Hierarchical Matrices For Large Spatial Datasets

abstract: We use available measurements to estimate the unknown parameters (variance, smoothness parameter, and covariance length) of a covariance function by maximizing the joint Gaussian log-likelihood function. To overcome cubic complexity in the linear algebra, we approximate the discretized covariance function in the hierarchical (H-) matrix format. The H-matrix format has a log-linear computational cost and storage $O(kn \log n)$, where the rank $k$ is a small integer and $n$ is the number of locations. The H-matrix technique allows us to work with general covariance matrices in an efficient way, since H-matrices can approximate inhomogeneous covariance functions, with a fairly general mesh that is not necessarily axes-parallel, and neither the covariance matrix itself nor its inverse have to be sparse. We demonstrate our method with Monte Carlo simulations and an application to soil moisture data. The C, C++ codes and data are freely available.


title: Multielectronic processes in particle and antiparticle collisions with rare gases

abstract: In this chapter we analyze the multiple ionization by impact of |Z|=1 projectiles: electrons, positrons, protons and antiprotons. Differences and similarities among the cross sections by these four projectiles allows us to have an insight on the physics involved. Mass and charge effects, energy thresholds, and relative importance of collisional and post-collisional processes are discussed. For this purpose, we performed a detailed theoretical-experimental comparison for single up to quintuple ionization of Ne, Ar, Kr and Xe by particles and antiparticles. We include an extensive compilation of the available data for the sixteen collisional systems, and the theoretical cross sections by means of the continuum distorted wave eikonal initial state approximation. We underline here that post-collisional ionization is decisive to describe multiple ionization by light projectiles, covering almost the whole energy range, from threshold to high energies. The normalization of positron and antiproton measurements to electron impact ones, the lack of data in certain cases, and the future prospects are presented and discussed.


title: Liveness-Driven Random Program Generation

abstract: Randomly generated programs are popular for testing compilers and program analysis tools, with hundreds of bugs in real-world C compilers found by random testing. However, existing random program generators may generate large amounts of dead code (computations whose result is never used). This leaves relatively little code to exercise a target compiler's more complex optimizations. To address this shortcoming, we introduce liveness-driven random program generation. In this approach the random program is constructed bottom-up, guided by a simultaneous structural data-flow analysis to ensure that the generator never generates dead code. The algorithm is implemented as a plugin for the Frama-C framework. We evaluate it in comparison to Csmith, the standard random C program generator. Our tool generates programs that compile to more machine code with a more complex instruction mix.


title: NLO exclusive diffractive processes with saturation

abstract: We present two NLO exclusive impact factors computed in the QCD shock wave approach. These are the very first steps towards precision studies of a wide range of high energy exclusive processes with saturation effects in $ep$, $eA$, $pp$ and $pA$ collisions.