title: Fermi Monitoring of Radio-loud Narrow Line Seyfert 1 Galaxies

abstract: We present detailed analysis of the gamma-ray flux variability and spectral properties of the five radio-loud narrow line Seyfert 1 (RL-NLSy1) galaxies, detected by the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope, namely 1H 0323+342, SBS 0846+513, PMN J0948+0022, PKS 1502+036, and PKS 2004$-$447. The first three sources show significant flux variations including the rapid variability of a few hours by 1H 0323+342. The average gamma-ray spectrum of 1H 0323+342 and PMN J0948+0022 shows deviation from a simple power law (PL) behavior, whereas for other three sources, PL model gives better fit. The spectra of 1H 0323+342, SBS 0846+513, and PMN J0948+0022, respectively in low, flaring, and moderately active states, show significant curvature. Such curvature in the gamma-ray spectrum of 1H 0323+342 and PMN J0948+0022 can be due to emission region located inside the broad line region (BLR) where the primary mechanism of the gamma-ray emission is inverse-Compton (IC) scattering of BLR photons occurring in the Klein-Nishina regime. The gamma-ray emission of SBS 0846+513 is explained by IC scattering of dusty torus photons which puts the emission region to be outside the BLR and thus under the Thomson regime. Therefore, the observed curvature of SBS 0846+513 could be intrinsic to the particle energy distribution. The presence of curvature in the gamma-ray spectrum and flux variability amplitudes of some of the RL-NLSy1 galaxies suggest that these sources could be akin to low/moderate jet power flat spectrum radio quasars.


title: Low-energy neutron-deuteron reactions with N3LO chiral forces

abstract: We solve three-nucleon Faddeev equations with nucleon-nucleon and three-nucleon forces derived consistently in the framework of chiral perturbation theory at next-to-next-to-next-to-leading order in the chiral expansion. In this first investigation we include only matrix elements of the three-nucleon force for partial waves with the total two-nucleon (three-nucleon) angular momenta up to 3 (5/2). Low-energy neutron-deuteron elastic scattering and deuteron breakup reaction are studied. Emphasis is put on Ay puzzle in elastic scattering and cross sections in symmetric-space-star and neutron-neutron quasi-free-scattering breakup configurations, for which large discrepancies between data and theory have been reported.


title: Canonical bases for the quantum linear supergroups

abstract: We give a combinatorial construction for the canonical bases of the $\pm$-parts of the quantum enveloping superalgebra $\bfU(\mathfrak{gl}_{m|n})$ and discuss their relationship with the Kazhdan-Lusztig bases for the quantum Schur superalgebras $\bsS(m|n,r)$ introduced in \cite{DR}. We will also extend this relationship to the induced bases for simple polynomial representations of $\bfU(\mathfrak{gl}_{m|n})$.


title: Impact of point defects on the electronic structure of paramagnetic CrN

abstract: This paper presents first principles calculations of paramagnetic cubic CrN$_x$ with the aim to provide deeper insight into recently published transmission electron microscopy-based study on this material. Among several types of point defects which may result in N-deficient material, N vacancy is found to be energetically preferred to Cr interstitial and anti-sites. Electron Energy Loss Near Edge Structure of N K-edge transition is calculated for various concentrations of N vacancies in CrN$_x$, yielding the same trends as experimentally observed. Analysis of the electronic structure reveals decreased charge transfer from Cr sites with increased N vacancy content, hence increasing the metallic character of the defected material. Finally, the electronic structure is found to be strongly dependent on the local environment (i.e. presence of the N vacancies).


title: cuDNN: Efficient Primitives for Deep Learning

abstract: We present a library of efficient implementations of deep learning primitives. Deep learning workloads are computationally intensive, and optimizing their kernels is difficult and time-consuming. As parallel architectures evolve, kernels must be reoptimized, which makes maintaining codebases difficult over time. Similar issues have long been addressed in the HPC community by libraries such as the Basic Linear Algebra Subroutines (BLAS). However, there is no analogous library for deep learning. Without such a library, researchers implementing deep learning workloads on parallel processors must create and optimize their own implementations of the main computational kernels, and this work must be repeated as new parallel processors emerge. To address this problem, we have created a library similar in intent to BLAS, with optimized routines for deep learning workloads. Our implementation contains routines for GPUs, although similarly to the BLAS library, these routines could be implemented for other platforms. The library is easy to integrate into existing frameworks, and provides optimized performance and memory usage. For example, integrating cuDNN into Caffe, a popular framework for convolutional networks, improves performance by 36% on a standard model while also reducing memory consumption.


title: Scheduling for Multi-Camera Surveillance in LTE Networks

abstract: Wireless surveillance in cellular networks has become increasingly important, while commercial LTE surveillance cameras are also available nowadays. Nevertheless, most scheduling algorithms in the literature are throughput, fairness, or profit-based approaches, which are not suitable for wireless surveillance. In this paper, therefore, we explore the resource allocation problem for a multi-camera surveillance system in 3GPP Long Term Evolution (LTE) uplink (UL) networks. We minimize the number of allocated resource blocks (RBs) while guaranteeing the coverage requirement for surveillance systems in LTE UL networks. Specifically, we formulate the Camera Set Resource Allocation Problem (CSRAP) and prove that the problem is NP-Hard. We then propose an Integer Linear Programming formulation for general cases to find the optimal solution. Moreover, we present a baseline algorithm and devise an approximation algorithm to solve the problem. Simulation results based on a real surveillance map and synthetic datasets manifest that the number of allocated RBs can be effectively reduced compared to the existing approach for LTE networks.


title: Change Point Detection in Correlation Networks

abstract: Many systems of interacting elements can be conceptualized as networks, where network nodes represent the elements and network ties represent interactions between the elements. In systems where the underlying network evolves in time, it is useful to determine the points in time where the network structure changes significantly as these may correspond also to functional change points. We propose a method for detecting these change points in correlation networks that, unlike previous change point detection methods designed for time series data, requires no distributional assumptions. We investigate the difficulty of change point detection near the boundaries of data in correlation networks and demonstrate the power of our method and a competing method through simulation. We also show the generalizable nature of our method by applying it to stock price data as well as fMRI data.


title: Predictability of IL-28B-polymorphism on protease-inhibitor-based triple-therapy in chronic HCV-genotype-1 patients: A meta-analysis

abstract: AIM: To investigate the predictability of interleukin-28B single nucleotide polymorphism rs12979860 with respect to sustained virological response (SVR) in chronically hepatitis C virus (HCV) genotype-1 patients treated with a protease-inhibitor and pegylated interferon-$\alpha$ (Peg-INF-$\alpha$) based triple-therapy. METHODS: We searched PubMed, the Cochrane Library and Web of Knowledge for studies regarding the interleukin 28B (IL-28B)-genotype and protease-inhibitor based triple-therapy. Ten studies with 2707 patients were included into this meta-analysis. We used regression methods in order to investigate determinants of SVR. RESULTS: IL-28B-CC-genotype patients achieved higher SVR rates (odds 5.34, CI: 3.81-7.49) than IL-28B-non-CC-genotype patients (1.88, CI: 1.43-2.48) receiving triple-therapy. The line of therapy (treatment-na\"ive or -experienced for Peg-INF-$\alpha$) did not affect the predictive value of IL-28B (p=0.1). IL-28B-CC-genotype patients treated with protease inhibitor-based triple-therapy consisting of Boceprevir, Simeprevir, Telaprevir or Vaniprevir showed odds of 1.86, 9.77, 4.51 and 0.89, respectively. The odds for CC genotype patients treated with Faldaprevir cannot be quantified, as only a single study with a 100% SVR rate was available. CONCLUSION: IL-28B-SNP predicts the outcome for chronic HCV genotype-1 patients receiving protease inhibitor-based triple-therapy. The predictive value varies between the different protease inhibitors.


title: A new ray-tracing scheme for 3D diffuse radiation transfer on highly parallel architectures

abstract: We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of mesh grids, and is devised to compute the radiation transfer along each light-ray completely in parallel with appropriate grouping of the light-rays. We find that the performance of our scheme scales well with the number of adopted CPU cores and GPUs, and also that our scheme is nicely parallelized on a multi-node system by adopting the multiple wave front scheme, and the performance scales well with the amount of the computational resources. As numerical tests to validate our scheme and to give a physical criterion for the angular resolution of our ray-tracing scheme, we perform several numerical simulations of the photo-ionization of neutral hydrogen gas by ionizing radiation sources without the "on-the-spot" approximation, in which the transfer of diffuse radiation by radiative recombination is incorporated in a self-consistent manner.


title: Vortices at the magnetic equator generated by hybrid Alfven resonant waves

abstract: We performed three-dimensional magnetohydrodynamic simulations of shear Alfven waves in a full field line system with magnetosphere-ionosphere coupling and plasma non-uniformities. Feedback instability of the Alfven resonant modes showed various nonlinear features under the field line cavities: i) a secondary flow shear instability occurs at the magnetic equator, ii) trapping of the ionospheric Alfven resonant modes facilitates deformation of field-aligned current structures, and iii) hybrid Alfven resonant modes grow to cause vortices and magnetic oscillations around the magnetic equator. Essential features in the initial brightening of auroral arc at substorm onsets could be explained by the dynamics of Alfven resonant modes, which are the nature of the field line system responding to a background rapid change.


title: On the Factorization of Rational Discrete-Time Spectral Densities

abstract: In this paper, we consider an arbitrary matrix-valued, rational spectral density $\Phi(z)$. We show with a constructive proof that $\Phi(z)$ admits a factorization of the form $\Phi(z)=W^\top (z^{-1})W(z)$, where $W(z)$ is stochastically minimal. Moreover, $W(z)$ and its right inverse are analytic in regions that may be selected with the only constraint that they satisfy some symplectic-type conditions. By suitably selecting the analyticity regions, this extremely general result particularizes into a corollary that may be viewed as the discrete-time counterpart of the matrix factorization method devised by Youla in his celebrated work (Youla, 1961).


title: On consecutive edge magic total labeling of connected bipartite graphs

abstract: Since Sedl\'{a}$\breve{\mbox{c}}$ek introduced the notion of magic labeling of a graph in 1963, a variety of magic labelings of a graph have been defined and studied. In this paper, we study consecutive edge magic labelings of a connected bipartite graph. We make a very useful observation that there are only four possible values of $b$ for which a connected bipartite graph has a $b$-edge consecutive magic labeling. On the basis of this fundamental result, we deduce various interesting results on consecutive edge magic labelings of bipartite graphs, especially caterpillars and lobsters, which extends the results given by Sugeng and Miller.


title: First-principles simulation of capacitive charging of graphene and implications for supercapacitor design

abstract: Supercapacitors store energy via the formation of an electric double layer, which generates a strong electric field at the electrode-electrolyte interface. Unlike conventional metallic electrodes, graphene-derived materials suffer from a low electronic density of states (i.e., quantum capacitance), which limits their ability to redistribute charge and efficiently screen this field. To explore these effects, we introduce a first-principles approach based on the effective screening medium framework, which is used to directly simulate the charge storage behavior of single- and multi-layered graphene in a way that more closely approximates operating devices. We demonstrate that the presence of the interfacial field significantly alters the capacitance in electrodes thinner than a few graphene layers, deriving in large part from intrinsic space-charge screening limitations. The capacitance is also found to be highly sensitive to the gap between the electrode and the solvent (contact layer), which offers possibilities for tuning the interfacial capacitance of the electrode by proper engineering of the electrolyte. Our results offer an alternative interpretation of discrepancies between experimental measurements and fixed-band models, and provide specific implications for improving graphene-based devices.


title: Space-Efficient Path-Reporting Approximate Distance Oracles

abstract: We consider approximate {\em path-reporting} distance oracles, distance labeling and labeled routing with extremely low space requirement, for general undirected graphs. For distance oracles, we show how to break the n\log n space bound of Thorup and Zwick if approximate {\em paths} rather than distances need to be reported. For approximate distance labeling and labeled routing, we break the previously best known space bound of O(log n) words per vertex. The cost for such space efficiency is an increased stretch.


title: Origin of the mixed-order transition in multiplex networks: the Ashkin-Teller model

abstract: Recently, diverse phase transition (PT) types have been obtained in multiplex networks, such as discontinuous, continuous, and mixed-order PTs. However, they emerge from individual systems, and there is no theoretical understanding of such PTs in a single framework. Here, we study a spin model called the Ashkin-Teller (AT) model in a mono-layer scale-free network; this can be regarded as a model of two species of Ising spin placed on each layer of a double-layer network. The four-spin interaction in the AT model represents the inter-layer interaction in the multiplex network. Diverse PTs emerge depending on the inter-layer coupling strength and network structure. Especially, we find that mixed-order PTs occur at the critical end points. The origin of such behavior is explained in the framework of Landau-Ginzburg theory.


title: Polynomial-time Algorithms for Weighted Efficient Domination Problems in AT-free Graphs and Dually Chordal Graphs

abstract: An efficient dominating set (or perfect code) in a graph is a set of vertices the closed neighborhoods of which partition the vertex set of the graph. The minimum weight efficient domination problem is the problem of finding an efficient dominating set of minimum weight in a given vertex-weighted graph; the maximum weight efficient domination problem is defined similarly. We develop a framework for solving the weighted efficient domination problems based on a reduction to the maximum weight independent set problem in the square of the input graph. Using this approach, we improve on several previous results from the literature by deriving polynomial-time algorithms for the weighted efficient domination problems in the classes of dually chordal and AT-free graphs. In particular, this answers a question by Lu and Tang regarding the complexity of the minimum weight efficient domination problem in strongly chordal graphs.


title: Stretching of BDT-gold molecular junctions: thiol or thiolate termination?

abstract: It is often assumed that the hydrogen atoms in the thiol groups of a benzene-1,4-dithiol dissociate when Au-benzene-1,4-dithiol-Au junctions are formed. We demonstrate, by stability and transport properties calculations, that this assumption can not be made. We show that the dissociative adsorption of methanethiol and benzene-1,4-dithiol molecules on a flat Au(111) surface is energetically unfavorable and that the activation barrier for this reaction is as high as 1 eV. For the molecule in the junction, our results show, for all electrode geometries studied, that the thiol junctions are energetically more stable than their thiolate counterparts. Due to the fact that density functional theory (DFT) within the local density approximation (LDA) underestimates the energy difference between the lowest unoccupied molecular orbital and the highest occupied molecular orbital by several electron-volts, and that it does not capture the renormalization of the energy levels due to the image charge effect, the conductance of the Au-benzene-1,4-dithiol-Au junctions is overestimated. After taking into account corrections due to image charge effects by means of constrained-DFT calculations and electrostatic classical models, we apply a scissor operator to correct the DFT energy levels positions, and calculate the transport properties of the thiol and thiolate molecular junctions as a function of the electrodes separation.


title: Electromagnetic modelling of superconductors with a smooth current-voltage relation: variational principle and coils from a few turns to large magnets

abstract: Many large-scale applications require electromagnetic modelling with extensive numerical computations, such as magnets or 3-dimensional (3D) objects like transposed conductors or motors and generators. Therefore, it is necessary to develop computationally time-efficient but still accurate numerical methods. This article develops a general variational formalism for any ${\bf E}({\bf J})$ relation and applies it to model coated-conductor coils containing up to thousands of turns, taking magnetization currents fully into account. The variational principle, valid for any 3D situation, restricts the computations to the sample volume, reducing the computation time. However, no additional magnetic materials interacting with the superconductor are taken directly into account. Regarding the coil modelling, we use a power law $E(J)$ relation with magnetic field-dependent critical current density, $J_c$, and power law exponent, $n$. We test the numerical model by comparing the results to analytical formulas for thin strips and experiments for stacks of pancake coils, finding a very good agreement. Afterwards, we model a magnet-size coil of 4000 turns (stack of 20 pancake coils of 200 turns each). We found that the AC loss is mainly due to magnetization currents. We also found that for an $n$ exponent of 20, the magnetization currents are greatly suppressed after 1 hour relaxation. In addition, in coated conductor coils magnetization currents have an important impact on the generated magnetic field; which should be taken into account for magnet design. In conclusion, the presented numerical method fulfills the requirements for electromagnetic design of coated conductor windings.


title: Comparing Apples and Oranges: Query Tradeoff in Submodular Maximization

abstract: Fast algorithms for submodular maximization problems have a vast potential use in applicative settings, such as machine learning, social networks, and economics. Though fast algorithms were known for some special cases, only recently Badanidiyuru and Vondr\'{a}k (2014) were the first to explicitly look for such algorithms in the general case of maximizing a monotone submodular function subject to a matroid independence constraint. The algorithm of Badanidiyuru and Vondr\'{a}k matches the best possible approximation guarantee, while trying to reduce the number of value oracle queries the algorithm performs. Our main result is a new algorithm for this general case which establishes a surprising tradeoff between two seemingly unrelated quantities: the number of value oracle queries and the number of matroid independence queries performed by the algorithm. Specifically, one can decrease the former by increasing the latter and vice versa, while maintaining the best possible approximation guarantee. Such a tradeoff is very useful since various applications might incur significantly different costs in querying the value and matroid independence oracles. Furthermore, in case the rank of the matroid is $O(n^c)$, where $n$ is the size of the ground set and $c$ is an absolute constant smaller than $1$, the total number of oracle queries our algorithm uses can be made to have a smaller magnitude compared to that needed by Badanidiyuru and Vondr\'{a}k. We also provide even faster algorithms for the well studied special cases of a cardinality constraint and a partition matroid independence constraint, both of which capture many real-world applications and have been widely studied both theorically and in practice.


title: Coherent forward scattering in 2D disordered systems

abstract: We present a detailed numerical and theoretical analysis of the recently discovered phenomenon of coherent forward scattering. This effect manifests itself as a macroscopic interference peak in the forward direction of the momentum distribution of a matter wave launched with finite velocity in a random potential. Focusing on the two-dimensional case, we show that coherent forward scattering generally arises due the confinement of the wave in a finite region of space, and explain under which conditions it can be seen as a genuine signature of Anderson localization.


title: A new formulation of higher parallel transport in higher gauge theory

abstract: In this technical paper, we present a new formulation of higher parallel transport in strict higher gauge theory required for the rigorous construction of Wilson lines and surfaces. Our approach is based on an original notion of Lie crossed module cocycle and cocycle 1- and 2-gauge transformation with a non standard double category theoretic interpretation. We show its equivalence to earlier formulations.


title: Almost-toric Hypersurfaces

abstract: An almost-toric hypersurface is parameterized by monomials multiplied by polynomials in one extra variable. We determine the Newton polytope of such a hypersurface, and apply this to give an algorithm for computing the implicit equation.


title: Degenerate affine Grassmannians and loop quivers

abstract: We study the connection between the affine degenerate Grassmannians in type $A$, quiver Grassmannians for one vertex loop quivers and affine Schubert varieties. We give an explicit description of the degenerate affine Grassmannian of type $GL_n$ and identify it with semi-infinite orbit closure of type $A_{2n-1}$. We show that principal quiver Grassmannians for the one vertex loop quiver provide finite-dimensional approximations of the degenerate affine Grassmannian. Finally, we give an explicit description of the degenerate affine Grassmannian of type $A_1^{(1)}$, propose a conjectural description in the symplectic case and discuss the generalization to the case of the affine degenerate flag varieties.


title: A characterization of barycentrically preassociative functions

abstract: We provide a characterization of the variadic functions which are barycentrically preassociative as compositions of length-preserving associative string functions with one-to-one unary maps. We also discuss some consequences of this characterization.


title: Solar neutrino with Borexino: results and perspectives

abstract: Borexino is a unique detector able to perform measurement of solar neutrinos fluxes in the energy region around 1 MeV or below due to its low level of radioactive background. It was constructed at the LNGS underground laboratory with a goal of solar $^{7}$Be neutrino flux measurement with 5\% precision. The goal has been successfully achieved marking the end of the first stage of the experiment. A number of other important measurements of solar neutrino fluxes have been performed during the first stage. Recently the collaboration conducted successful liquid scintillator repurification campaign aiming to reduce main contaminants in the sub-MeV energy range. With the new levels of radiopurity Borexino can improve existing and challenge a number of new measurements including: improvement of the results on the Solar and terrestrial neutrino fluxes measurements; measurement of pp and CNO solar neutrino fluxes; search for non-standard interactions of neutrino; study of the neutrino oscillations on the short baseline with an artificial neutrino source (search for sterile neutrino) in context of SOX project.


title: Small ball estimates for quasi-norms

abstract: This note contains two types of small ball estimates for random vectors in finite dimensional spaces equipped with a quasi-norm. In the first part, we obtain bounds for the small ball probability of random vectors under some smoothness assumptions on their density function. In the second part, we obtain Littlewood-Offord type estimates for quasi-norms. This generalizes a result which was previously obtained by Friedland and Sodin and by Rudelson and Vershynin.


title: SimNets: A Generalization of Convolutional Networks

abstract: We present a deep layered architecture that generalizes classical convolutional neural networks (ConvNets). The architecture, called SimNets, is driven by two operators, one being a similarity function whose family contains the convolution operator used in ConvNets, and the other is a new soft max-min-mean operator called MEX that realizes classical operators like ReLU and max pooling, but has additional capabilities that make SimNets a powerful generalization of ConvNets. Three interesting properties emerge from the architecture: (i) the basic input to hidden layer to output machinery contains as special cases kernel machines with the Exponential and Generalized Gaussian kernels, the output units being "neurons in feature space" (ii) in its general form, the basic machinery has a higher abstraction level than kernel machines, and (iii) initializing networks using unsupervised learning is natural. Experiments demonstrate the capability of achieving state of the art accuracy with networks that are an order of magnitude smaller than comparable ConvNets.


title: Multi-resource fairness: Objectives, algorithms and performance

abstract: Designing efficient and fair algorithms for sharing multiple resources between heterogeneous demands is becoming increasingly important. Applications include compute clusters shared by multi-task jobs and routers equipped with middleboxes shared by flows of different types. We show that the currently preferred objective of Dominant Resource Fairness has a significantly less favorable efficiency-fairness tradeoff than alternatives like Proportional Fairness and our proposal, Bottleneck Max Fairness. In addition to other desirable properties, these objectives are equally strategyproof in any realistic scenario with dynamic demand.


title: Lee-Yang zero distribution of high temperature QCD and Roberge-Weiss phase transition

abstract: Canonical partition functions and Lee-Yang zeros of QCD at finite density and high temperature are studied. Recent lattice simulations have confirmed that the free energy of QCD is a quartic function of quark chemical potential at temperature slightly above pseudo-critical temperature $T_c$, as in the case with a gas of free massless fermions. We present analytic derivation of the canonical partition functions and Lee-Yang zeros for this type of free energy using the saddle point approximation. We also perform lattice QCD simulation in a canonical approach using the fugacity expansion of the fermion determinant, and carefully examine its reliability. By comparing the analytic and numerical results, we conclude that the canonical partition functions follow the Gaussian distribution of the baryon number, and the accumulation of Lee-Yang zeros of these canonical partition functions exhibit the first-order Roberge-Weiss phase transition. We discuss the validity and applicable range of the result and its implications both for theoretical and experimental studies.


title: The Very Fast Evolution of the VLTP Object V4334 Sgr

abstract: V4334 Sgr (Sakurai's object) is an enigmatic evolved star that underwent a very late thermal pulse a few years before its discovery in 1996. It ejected a new, hydrogen-deficient nebula in the process. Emission lines from the newly ejected gas were first discovered in 1998 (He I 1083 nm) and 2001 (optical). We have monitored the optical emission spectrum since. From 2001 through 2007 the optical spectrum showed an exponential decline in flux, consistent with a shock that occurred around 1998 and started cooling soon after that. In this paper we show that since 2008 the line fluxes have been continuously rising again. Our preliminary interpretation is that this emission comes from a region close to the central star, and is excited by a second shock. This shock may have been induced by an increase in the stellar mass loss and wind velocity associated with a rise in the stellar temperature.


title: A Posteriori Error Bounds for Two Point Boundary Value Problems: A Green's Function Approach

abstract: We present a computer assisted method for generating existence proofs and a posteriori error bounds for solutions to two point boundary value problems (BVPs). All truncation errors are accounted for and, if combined with interval arithmetic to bound the rounding errors, the computer generated results are mathematically rigorous. The method is formulated for $n$-dimensional systems and does not require any special form for the vector field of the differential equation. It utilizes a numerically generated approximation to the BVP fundamental solution and Green's function and thus can be applied to stable BVPs whose initial value problem is unstable. The utility of the method is demonstrated on a pair of singularly perturbed model BVPs and by using it to rigorously show the existence of a periodic orbit in the Lorenz system.


title: Malliavin and flow regularity of SDEs. Application to the study of densities and the stochastic transport equation

abstract: In this work we present a condition for the regularity, in both space and Malliavin sense, of strong solutions to SDEs driven by Brownian motion. We conjecture that this condition is optimal. As a consequence, we are able to improve the regularity of densities of such solutions. We also apply these results to construct a classical solution to the stochastic transport equation when the drift is Lipschitz.


title: Weak Value and the Wave-Particle Duality

abstract: The weak value, introduced by Aharonov et al. to extend the conventional scope of physical observables in quantum mechanics, is an intriguing concept which sheds new light on quantum foundations and is also useful for precision measurement, but it poses serious questions on its physical meaning due to the unconventional features including the complexity of its value. In this paper we point out that the weak value has a direct connection with the wave-particle duality, in the sense that the wave nature manifests itself in the imaginary part while the particle nature in the real part. This is illustrated by the double slit experiment, where we argue, with no conflict with complementarity, that the trajectory of the particle can be inferred based on the weak value without destroying the interference.


title: The thermodynamics of a black hole in equilibrium implies the breakdown of Einstein equations on a macroscopic near-horizon shell

abstract: We study a black hole of mass $M$, enclosed within a spherical box, in equilibrium with its Hawking radiation. We show that the spacetime geometry inside the box is described by the Oppenheimer-Volkoff equations for radiation, except for a {\em thin shell} around the horizon. We use the maximum entropy principle to show that the invariant width of the shell is of order $\sqrt{M}$, its entropy is of order $M$ and its temperature of order $1/\sqrt{M}$ (in Planck units). Thus, the width of the shell is {\em much larger} than the Planck length. Our approach is to insist on thermodynamic consistency when classical general relativity coexists with the Hawking temperature in the description of a gravitating system. No assumptions about an underlying theory are made and no restrictions are placed on the origins of the new physics near the horizon. We only employ classical general relativity and the principles of thermodynamics. Our result is strengthened by an analysis of the trace anomaly associated to the geometry inside the box, i.e., the regime where quantum field effects become significant corresponds to the shells of maximum entropy around the horizon.


title: 3alpha clustering in the excited states of 16C

abstract: The alpha cluster states of 16C are investigated by using the antisymmetrized molecular dynamics. It is shown that two different types of alpha cluster states exist: triangular and linear-chain states. The former has an approximate isosceles triangular configuration of alpha particles surrounded by four valence neutrons occupying sd-shell, while the latter has the linearly aligned alpha particles with two sd-shell neutrons and two pf-shell neutrons. It is found that the structure of the linear-chain state is qualitatively understood in terms of the 3/2 pi- and 1/2 sigma- molecular orbit as predicted by molecular-orbital model, but there exists non-negligible Be+alpha+2n correlation. The band-head energies of the triangular and linear-chain rotational bands are 8.0 and 15.5 MeV, and the latter is close to the He+Be threshold energy. It is also shown that the linear-chain state becomes the yrast sstate at J=10 with excitation energy 27.8 MeV owing to its very large moment-of-inertia comparable with hyperdeformation.


title: The Murchison Widefield Array Commissioning Survey: A Low-Frequency Catalogue of 14,110 Compact Radio Sources over 6,100 Square Degrees

abstract: We present the results of an approximately 6,100 square degree 104--196MHz radio sky survey performed with the Murchison Widefield Array during instrument commissioning between 2012 September and 2012 December: the Murchison Widefield Array Commissioning Survey (MWACS). The data were taken as meridian drift scans with two different 32-antenna sub-arrays that were available during the commissioning period. The survey covers approximately 20.5 h < Right Ascension (RA) < 8.5 h, -58 deg < Declination (Dec) < -14 deg over three frequency bands centred on 119, 150 and 180 MHz, with image resolutions of 6--3 arcmin. The catalogue has 3-arcmin angular resolution and a typical noise level of 40 mJy/beam, with reduced sensitivity near the field boundaries and bright sources. We describe the data reduction strategy, based upon mosaiced snapshots, flux density calibration and source-finding method. We present a catalogue of flux density and spectral index measurements for 14,110 sources, extracted from the mosaic, 1,247 of which are sub-components of complexes of sources.


title: Resolving the Clumpy Structure of the Outflow Winds in the Gravitationally Lensed Quasar SDSS J1029+2623

abstract: We study the geometry and the internal structure of the outflowing wind from the accretion disk of a quasar by observing multiple sightlines with the aid of strong gravitational lensing. Using Subaru/HDS, we performed high-resolution ($R$ $\sim$ 36,000) spectroscopic observations of images A and B of the gravitationally lensed quasar SDSS J1029+2623 (at $z_{em}$ $\sim$ 2.197) whose image separation angle, $\theta$ $\sim$ 22$^{\prime\prime}\!\!$.5, is the largest among those discovered so far. We confirm that the difference in absorption profiles in the images A and B discovered by Misawa et al. (2013) remains unchanged since 2010, implying the difference is not due to time variability of the absorption profiles over the delay between the images, $\Delta t$ $\sim$ 744 days, but rather due to differences along the sightlines. We also discovered time variation of C IV absorption strength in both images A and B, due to change of ionization condition. If a typical absorber's size is smaller than its distance from the flux source by more than five orders of magnitude, it should be possible to detect sightline variations among images of other smaller separation, galaxy-scale gravitationally lensed quasars.


title: A Hilbert Space setting for higher spin interactions which replaces Gauge Theory

abstract: The recently discovered Hilbert space description of renormalizable interactions of higher spin (equal or bigger than 1) fields requires to replace the pointlocal s=1 vectorpotentials of indefinite metric (Krein space) BRST gauge theory by their stringlike counterpart in Hilbert space. It is shown that the Hilbert space positivity leads to new properties outside the conceptual range of the gauge theoretic description: topological aspects of Wilson loops, induced normalization terms (in particular Mexican hat type potentials for massive vectormesons coupled to Hermitian scalar field) and a possible role of string-localization in confinerment.


title: On the maximization of likelihoods belonging to the exponential family using ideas related to the Levenberg-Marquardt approach

abstract: The Levenberg-Marquardt algorithm is a flexible iterative procedure used to solve non-linear least squares problems. In this work we study how a class of possible adaptations of this procedure can be used to solve maximum likelihood problems when the underlying distributions are in the exponential family. We formally demonstrate a local convergence property and we discuss a possible implementation of the penalization involved in this class of algorithms. Applications to real and simulated compositional data show the stability and efficiency of this approach.


title: Maximization of surface-enhanced transversal magneto-optic Kerr effect in Au/Co/Au thin films

abstract: In order to maximize the transversal magneto-optic Kerr effect (T-MOKE) of a Au/Co/Au structure, we propose a method to obtain the optimum thickness values. A criteria based on preserving good plasmonic properties has been included as part of this method. Using the theoretical prediction, we grew Au/Co/Au trilayers and perform optical andMO characterization using the Kretschmann configuration. The results admit very easy interpretation in terms of the interaction between the magneto-optical and plasmonic properties dictating the optimal thicknesses of the structure. Moreover, we have grown and characterized the optimized structure finding good agreement with theory reaching, for a 532nm green laser, amaximal surface magneto-optic (MO) signal enhancement of close to nine folds with respect to the signal without plasmonic excitation.


title: Squarefree monomial modules and extremal Betti numbers

abstract: Let K be a field and S a polynomial ring in a finite number of variables over K. Let F be a finitely generated graded free S-module. We examine some classes of squarefree monomial submodules of F. Hence, we focalize our attention on the Betti table of such classes in order to analyze the behavior of their extremal Betti numbers.


title: Nodal discontinuous Galerkin methods for fractional diffusion equations on 2D domain with triangular meshes

abstract: This paper, as the sequel to previous work, develops numerical schemes for fractional diffusion equations on a two-dimensional finite domain with triangular meshes. We adopt the nodal discontinuous Galerkin methods for the full spatial discretization by the use of high-order nodal basis, employing multivariate Lagrange polynomials defined on the triangles. Stability analysis and error estimates are provided, which shows that if polynomials of degree $N$ are used, the methods are (N+1)-th order accurate for general triangulations. Finally, the performed numerical experiments confirm the optimal order of convergence.


title: Relaxation of regularity for the Westervelt equation by nonlinear damping with application in acoustic-acoustic and elastic-acoustic coupling

abstract: In this paper we show local (and partially global) in time existence for the Westervelt equation with several versions of nonlinear damping. This enables us to prove well-posedness with spatially varying $L_\infty$-coefficients, which includes the situation of interface coupling between linear and nonlinear acoustics as well as between linear elasticity and nonlinear acoustics, as relevant, e.g., in high intensity focused ultrasound (HIFU) applications.


title: Sharp magnetization jump at the first-order superconducting transition in Sr2RuO4

abstract: The magnetization and magnetic torque of a high-quality single crystal of Sr$_2$RuO$_4$ have been measured down to 0.1 K under a precise control of the magnetic-field orientation. When the magnetic field is applied exactly parallel to the $ab$ plane, a sharp magnetization jump $4\pi\delta M$ of $(0.74 \pm 0.15)$ G at the upper critical field $H_{{\rm c2},{ab}} \sim 15$ kOe with a field hysteresis of 100 Oe is observed at low temperatures, evidencing a first-order superconducting-normal transition. A strong magnetic torque appearing when $H$ is slightly tilted away from the $ab$ plane confirms an intrinsic anisotropy $\varGamma=\xi_a/\xi_c$ of as large as 60 even at 100 mK, in contrast with the observed $H_{{\rm c2}}$ anisotropy of $\sim 20$. The present results raise fundamental issues in both the existing spin-triplet and spin-singlet scenarios, providing, in turn, crucial hints toward the resolution of the superconducting nature of Sr$_2$RuO$_4$.


title: Optimization of radiation hardness and charge collection of edgeless silicon pixel sensors for photon science

abstract: Recent progress in active-edge technology of silicon sensors enables the development of large-area tiled silicon pixel detectors with small dead space between modules by utilizing edgeless sensors. Such technology has been proven in successful productions of ATLAS and Medipix-based silicon pixel sensors by a few foundries. However, the drawbacks of edgeless sensors are poor radiation hardness for ionizing radiation and non-uniform charge collection by edge pixels. In this work, the radiation hardness of edgeless sensors with different polarities has been investigated using Synopsys TCAD with X-ray radiation-damage parameters implemented. Results show that if no conventional guard ring is present, none of the current designs are able to achieve a high breakdown voltage (typically < 30 V) after irradiation to a dose of ~10 MGy. In addition, a charge-collection model has been developed and was used to calculate the charges collected by the edge pixels of edgeless sensors when illuminated with X-rays. The model takes into account the electric field distribution inside the pixel sensor, the absorption of X-rays, drift and diffusion of electrons and holes, charge sharing effect, and threshold settings in ASICs. It is found that the non-uniform charge collection of edge pixels is caused by the strong bending of electric field and the non-uniformity depends on bias voltage, sensor thickness and distance from active edge to the last pixel ("edge space"). In particular, the last few pixels close to the active edge of the sensor are not sensitive to low-energy X-rays (< 10 keV) especially for sensors with thicker Si and smaller edge space. The results from the model calculation have been compared to measurements and good agreement was obtained. The model has been used to optimize the edge design.


title: Siklos waves with torsion in 3D

abstract: Starting from the Siklos waves in general relativity with a cosmological constant, interpreted as gravitational waves on the anti-de Sitter background, a new class of exact torsion waves is constructed in the framework of three-dimensional gravity with propagating torsion. In the asymptotic limit, the geometry of torsion waves takes the anti-de Sitter form. In the sector with massless torsion modes, we found a set of asymptotic conditions that leads to the conformal asymptotic symmetry.


title: Critical behavior and continuum scaling of $3D$ $Z(N)$ lattice gauge theories

abstract: Three-dimensional $Z(N)$ lattice gauge theories are studied numerically at finite temperature for $N$ = 5, 6, 8, 12, 13, 20 and for $N_t$=2,4,8. For each model the location of phase transitions and its critical indices are determined. The scaling of critical points with $N$ is proposed. The data obtained enable us to verify the scaling near the continuum limit for the $Z(N)$ models at finite temperatures.


title: Noether analysis of Scalar-Tensor Cosmology

abstract: A scalar--tensor theory of gravity, containing an arbitrary coupling function $F(\phi)$ and a general potential $V(\phi)$, is considered in the context of a spatially flat FLRW model. The use of reparametrization invariance enables a particular lapse parametrization in which the mini--superspace metric completely specifies the dynamics of the system. A requirement of existence of the maximal possible number of autonomous integrals of motion is imposed. This leads to a flat mini--superspace metric realized by a particular relation between the coupling function and the potential. The space of solutions is completely described in terms of the three autonomous integrals of motion constructed by the Killing fields of the mini--supermetric and an additional rheonomous emanating from the homothetic field. The solutions contain the arbitrary function which remains after the imposition of the relation between $F(\phi)$ and $V(\phi)$. To exemplify the use of the general results, we select some particular cases and study their physical implications through an effective energy--momentum tensor, which tends out to be that of a perfect fluid.


title: Pressure pulsations in reciprocating pump piping systems Part 1: Modelling

abstract: A distributed parameter model of pipeline transmission line behaviour is presented, based on a Galerkin method incorporating frequency-dependent friction. This is readily interfaced to an existing model of the pumping dynamics of a plunger pump to allow time-domain simulations of pipeline pressure pulsations in both suction and delivery lines. A new model for the pump inlet manifold is also proposed.


title: Multi-step Uniformization with Steady-State Detection in Nonstationary M/M/s Queuing Systems

abstract: A new approach to the steady state detection in the uniformization method of solving continuous time Markov chains is introduced. The method is particularly useful in solving inhomogenous CTMC's in multiple steps, where the desired error bound of the whole solution can be distributed not proportionally to the lengths of the respective intervals, but rather in a way, that maximizes the chances of detecting a steady state. Additionally, the convergence properties of the underlying DTMC are used to further enhance the computational savings due to the steady state detection. The method is applied to the problem of modeling a Call Center using inhomogenous CTMC model of a M(t)/M(t)/s(t) queuing system.


title: CuFeAs: A New Member in the 111-Family of Iron-Pnictides

abstract: We have synthesized CuFeAs, a new iron-pnictide compound with a layered tetragonal Cu2Sb type structure (space group P4/nmm: a = b = 3.7442(2) {\AA} and c = 5.8925(4) {\AA}) that is identical to that of 111-type iron-based superconductors. Our measurements suggest that in low applied magnetic field it undergoes an antiferromagnetic transition below TN ~ 9 K. When compared with the ground state of CuFeSb, recently reported 111-type ferromagnetic material (TC ~ 375 K), it has important implication with regard to the nature of Fe-Fe magnetic interaction in Fe-As materials. CuFeAs does not exhibit superconductivity down to 2 K.


title: A random pointwise ergodic theorem with Hardy field weights

abstract: Let $a_n$ be the random increasing sequence of natural numbers which takes each value independently with probability $n^{-a}$, $0 < a < 1/2$, and let $p(n) = n^{1+\epsilon}$, $0 < \epsilon < 1$. We prove that, almost surely, for every measure-preserving system $(X,T)$ and every $f \in L^1(X)$ the modulated, random averages \[ \frac{1}{N} \sum_{n = 1}^N e(p(n)) T^{a_n(\omega)} f\] converge to $0$ pointwise almost everywhere.


title: Functional limit theorems for the Bouchaud trap model with slowly varying traps

abstract: We consider the Bouchaud trap model on the integers in the case that the trap distribution has a slowly varying tail at infinity. Our main result is a functional limit theorem for the model under the annealed law, analogous to the functional limit theorems previously established in the literature in the case of integrable or regularly varying trap distribution. Reflecting the fact that the clock process is dominated in the limit by the contribution from the deepest-visited trap, the limit process for the model is a spatially-subordinated Brownian motion whose associated clock process is an extremal process.


title: Explicit Construction of Equivalence Bimodules between Noncommutative Solenoids

abstract: Let $p\in \mathbb{N}$ be prime, and let $\theta$ be irrational. The authors have previously shown that the noncommutative $p$-solenoid corresponding to the multiplier of the group $\left(\mathbb{Z}\left[\frac{1}{p}\right]\right)^2$ parametrized by $\alpha=(\theta+1, (\theta+1)/p,\cdots, (\theta+1)/p^j,\cdots )$ is strongly Morita equivalent to the noncommutative solenoid on $\left(\mathbb{Z}\left[\frac{1}{p}\right]\right)^2$ coming from the multiplier $\beta= (1-\frac{\theta+1}{\theta},1-\frac{\theta+1}{p\theta}, \cdots, 1-\frac{\theta+1}{p^j\theta}, \cdots )$ . The method used a construction of Rieffel referred to as the "Heisenberg bimodule" in which the two noncommutative solenoid corresponds to two different twisted group algebras associated to dual lattices in $(\mathbb Q_p\times \mathbb R)^2.$ In this paper, we make three additional observations: first, that at each stage, the subalgebra given by the irrational rotation algebra corresponding to $\alpha_{2j}=(\theta+1)/p^{2j}$ is strongly Morita equivalent to the irrational rotation algebra corresponding to the irrational rotation algebra corresponding to $\beta_{2j}= 1-\frac{\theta+1}{p^{2j}\theta}$ by a different construction of Rieffel, secondly, that that Rieffel's Heisenberg module relating the two non commutative solenoids can be constructed as the closure of a nested sequence of function spaces associated to a multiresolution analysis for a $p$-adic wavelet, and finally, at each stage, the equivalence bimodule between $A_{\alpha_{2j}}$ and $A_{\beta_{2j}}$ can be identified with the subequivalence bimodules arising from the $p$-adic MRA. Aside from its instrinsic interest, we believe this construction will guide us in our efforts to show that certain necessary conditions for two noncommutative solenoids to be strongly Morita equivalent are also sufficient.


title: Do the constants of nature couple to strong gravitational fields?

abstract: Recently, white dwarf stars have found a new use in the fundamental physics community. Many prospective theories of the fundamental interactions of Nature allow traditional constants, like the fine structure constant $\alpha$, to vary in some way. A study by Berengut et al. (2013) used the Fe/Ni V line measurements made by Preval et al. (2013) from the hot DA white dwarf G191-B2B, in an attempt to detect any variation in $\alpha$. It was found that the Fe V lines indicated an increasing alpha, whereas the Ni V lines indicated a decreasing alpha. Possible explanations for this could be misidentification of the lines, inaccurate atomic data, or wavelength dependent distortion in the spectrum. We examine the first two cases by using a high S/N reference spectrum from the hot sdO BD+28$^{\circ}$4211 to calibrate the Fe/Ni V atomic data. With this new data, we re-evaluate the work of Berengut et al. (2013) to derive a new constraint on the variation of alpha in a gravitational field.


title: New Candidate Eruptive Young Stars in Lynds 1340

abstract: We report on the discovery of three candidate eruptive young stars, found during our comprehensive multi-wavelength study of the young stellar population of the dark cloud L1340. These stars are as follows. (1) IRAS 02224+7227 (2MASS 02270555+7241167, HH 487 S) exhibited FUor-like spectrum in our low-resolution optical spectra. The available photometric data restrict its luminosity to 23 L_sun < L_bol < 59 L_sun. (2) 2MASS 02263797+7304575, identified as a classical T Tauri star during our H alpha survey, exhibited an EXor type brightening in 2005 November, at the time of the SDSS observations of the region. (3) 2MASS 02325605+7246055, a low-mass embedded young star, associated with a fan-shaped infrared nebula, underwent an outburst between the DSS1 and DSS2 surveys, leading to the appearance of a faint optical nebula. Our [SII] and H alpha images, as well as the Spitzer IRAC 4.5 micron images revealed Herbig-Haro objects associated with this star. Our results suggest that amplitudes and time scales of outbursts do not necessarily correlate with the evolutionary stage of the stars.


title: A possible solution to the Lyman/Balmer line problem in hot DA white dwarfs

abstract: Arguably, the best method for determining the effective temperature ($T_{\mathrm{eff}}$) and surface gravity (log $g$) of a DA white dwarf is by fitting the Hydrogen Lyman and Balmer absorption features. However, as has been shown for white dwarfs with $T_{\mathrm{eff}}$>50,000K, the calculated value from the Lyman and Balmer lines are discrepant, which worsens with increasing temperature. Many different solutions have been suggested, ranging from the input physics used to calculate the models, to interstellar reddening. We will focus on the former, and consider three variables. The first is the atomic data used, namely the number of transitions included in line blanketing treatments and the photoionization cross sections. The second is the stark broadening treatment used to synthesise the Lyman and Balmer line profiles, namely the calculations performed by Lemke (1997) and Tremblay & Bergeron (2009). Finally, the third is the atmospheric content. The model grids are calculated with a pure H composition, and a metal polluted composition using the abundances of Preval et al. (2013). We present the preliminary results of our analysis, whereby we have determined the $T_{\mathrm{eff}}$ for a small selection of white dwarfs. We plan to extend our analysis by allowing metallicity to vary in future model grids.


title: The upgrade of the LHCb Vertex Locator

abstract: The LHCb experiment is set for a significant upgrade, which will be ready for Run 3 of the LHC in 2020. This upgrade will allow LHCb to run at a significantly higher instantaneous luminosity and collect an integrated luminosity of $50\,\text{fb}^{-1}$ by the end of Run 4. In this process the Vertex locator (VELO) detector will be upgraded to a pixel-based silicon detector. The upgraded VELO will improve upon the current detector by being closer to the beam and having lower material modules with microchannel cooling and a thinner RF-foil. Simulations have shown that it will maintain its excellent performance, even after the radiation damage caused by collecting an integrated luminosity of $50\,\text{fb}^{-1}$.


title: Gaussian Tree Constraints Applied to Acoustic Linguistic Functional Data

abstract: Evolutionary models of languages are usually considered to take the form of trees. With the development of so-called tree constraints the plausibility of the tree model assumptions can be addressed by checking whether the moments of observed variables lie within regions consistent with trees. In our linguistic application, the data set comprises acoustic samples (audio recordings) from speakers of five Romance languages or dialects. We wish to assess these functional data for compatibility with a hereditary tree model at the language level. A novel combination of canonical function analysis (CFA) with a separable covariance structure provides a method for generating a representative basis for the data. This resulting basis is formed of components which emphasize language differences whilst maintaining the integrity of the observational language-groupings. A previously unexploited Gaussian tree constraint is then applied to component-by-component projections of the data to investigate adherence to an evolutionary tree. The results indicate that while a tree model is unlikely to be suitable for modeling all aspects of the acoustic linguistic data, certain features of the spoken Romance languages highlighted by the separable-CFA basis may indeed be suitably modeled as a tree.


title: Why has spacetime torsion such negligible effect on our universe?

abstract: We attempt an answer to the question as to why the evolution of four-dimensional universe is governed by spacetime curvature but not torsion. An answer is found if there is an additional compact spacelike dimension with a warped geometry, with torsion caused by a Kalb-Ramond (KR) antisymmetric tensor field in the bulk. Starting from a Randall-Sundrum type of warped extra dimension, and including the inevitable back reaction ensuing from the radius stabilization mechanism, we show that there is always an extra exponential suppression of the KR field on the four-dimensional projection that constitutes our visible universe. The back reaction is found to facilitate the process of such suppression.


title: Invariants of third-order ordinary differential equations $y'''=f(x,y,y',y'')$ via fiber preserving transformations

abstract: Bagderina \cite{Bagderina2008} solved the equivalence problem for scalar third-order ordinary differential equations (ODEs), quadratic in the second-order derivative, via point transformations. However, the question is open for the general class $y'''=f(x,y,y',y'')$ which is not quadratic in the second-order derivative. We utilize Lie's infinitesimal method to study the differential invariants of this general class under pseudo-group of fiber preserving equivalence transformations $\bar{x}=\phi(x), \bar{y}=\psi(x,y)$. As a result, all third-order differential invariants of this group and the invariant differentiation operators are determined. This leads to simple necessary explicit conditions for a third-order ODE to be equivalent to the respective canonical form under the considered group of transformations. Applications motivated by the literature are presented.


title: Comment on 'Enhancement in the production of nucleating clusters due to dimethylamine and large uncertainties in the thermochemistry of amine-enhanced nucleation' by Nadykto et al., Chem. Phys. Lett., 609, 42-49 (2014)

abstract: We comment on a study by Nadykto et al. recently published in this journal. Earlier work from our group has been misrepresented in this study, and we feel that the claims made need to be amended. Also the analysis of Nadykto et al. concerning the implications of their own density functional calculations is incomplete. We present cluster formation simulations allowing more conclusions to be drawn from their data, and also compare them to recent experimental results not cited in their work.


title: Second order statistics of robust estimators of scatter. Application to GLRT detection for elliptical signals

abstract: A central limit theorem for bilinear forms of the type $a^*\hat{C}_N(\rho)^{-1}b$, where $a,b\in{\mathbb C}^N$ are unit norm deterministic vectors and $\hat{C}_N(\rho)$ a robust-shrinkage estimator of scatter parametrized by $\rho$ and built upon $n$ independent elliptical vector observations, is presented. The fluctuations of $a^*\hat{C}_N(\rho)^{-1}b$ are found to be of order $N^{-\frac12}$ and to be the same as those of $a^*\hat{S}_N(\rho)^{-1}b$ for $\hat{S}_N(\rho)$ a matrix of a theoretical tractable form. This result is exploited in a classical signal detection problem to provide an improved detector which is both robust to elliptical data observations (e.g., impulsive noise) and optimized across the shrinkage parameter $\rho$.


title: Feature Learning from Incomplete EEG with Denoising Autoencoder

abstract: An alternative pathway for the human brain to communicate with the outside world is by means of a brain computer interface (BCI). A BCI can decode electroencephalogram (EEG) signals of brain activities, and then send a command or an intent to an external interactive device, such as a wheelchair. The effectiveness of the BCI depends on the performance in decoding the EEG. Usually, the EEG is contaminated by different kinds of artefacts (e.g., electromyogram (EMG), background activity), which leads to a low decoding performance. A number of filtering methods can be utilized to remove or weaken the effects of artefacts, but they generally fail when the EEG contains extreme artefacts. In such cases, the most common approach is to discard the whole data segment containing extreme artefacts. This causes the fatal drawback that the BCI cannot output decoding results during that time. In order to solve this problem, we employ the Lomb-Scargle periodogram to estimate the spectral power from incomplete EEG (after removing only parts contaminated by artefacts), and Denoising Autoencoder (DAE) for learning. The proposed method is evaluated with motor imagery EEG data. The results show that our method can successfully decode incomplete EEG to good effect.


title: The Formation of Kappa-Distribution Accelerated Electron Populations in Solar Flares

abstract: Driven by recent RHESSI observations of confined loop-top hard X-ray sources in solar flares, we consider stochastic acceleration of electrons in the presence of Coulomb collisions. If electron escape from the acceleration region can be neglected, the electron distribution function is determined by a balance between diffusive acceleration and collisions. Such a scenario admits a stationary solution for the electron distribution function that takes the form of a kappa distribution. We show that the evolution toward this kappa distribution involves a "wave front" propagating forwards in velocity space, so that electrons of higher energy are accelerated later; the acceleration time scales with energy according to $\tau_{\rm acc} \sim E^{3/2}$. At sufficiently high energies escape from the finite-length acceleration region will eventually dominate. For such energies, the electron velocity distribution function is obtained by solving a time-dependent Fokker-Planck equation in the "leaky-box" approximation. Solutions are obtained in the limit of a small escape rate from an acceleration region that can effectively be considered a thick target.


title: $B_s$ mesons: semileptonic and nonleptonic decays

abstract: In this contribution we compute some nonleptonic and semileptonic decay widths of $B_s$ mesons, working in the context of constituent quark models \cite{Albertus:2014gba, Albertus:2014bfa}. For the case of semileptonic decays we consider reactions leading to kaons or different $J^\pi$ $D_s$ mesons. The study of nonleptonic decays has been done in the factorisation approximation and includes the final states enclosed in Table 2.


title: Topological nature and the multiple Dirac cones hidden in Bismuth high-Tc superconductors

abstract: Recent theoretical studies employing density-functional theory have predicted BaBiO$_{3}$ (when doped with electrons) and YBiO$_{3}$ to become a topological insulator (TI) with a large topological gap (~ 0.7 eV). This, together with the natural stability against surface oxidation, makes the Bismuth-Oxide family of special interest for possible applications in quantum information and spintronics. The central question, we study here, is whether the hole-doped Bismuth Oxides, i.e. Ba$_{1-x}$K$_{x}$BiO$_{3}$ and BaPb$_{1-x}$Bi$_{x}$O$_{3}$, which are "high-Tc" bulk superconducting near 30 K, additionally display in the further vicinity of their Fermi energy $E_{F}$ a topological gap with a Dirac-type of topological surface state. Our electronic structure calculations predict the K-doped family to emerge as a TI, with a topological gap above $E_{F}$. Thus, these compounds can become superconductors with hole-doping and potential TIs with additional electron doping. Furthermore, we predict the Bismuth-Oxide family to contain an additional Dirac cone below $E_{F}$ for further hole doping, which manifests these systems to be candidates for both electron- and hole-doped topological insulators.


title: Reentrance of Berezinskii-Kosterlitz-Thouless-like transitions in three-state Potts antiferromagnetic thin film

abstract: Using Monte Carlo simulations and finite-size scaling, we study three-state Potts antiferromagnet on layered square lattice with two and four layers $L_z=2$ and $4$. As temperature decreases, the system develops quasi-long-range order via a Berezinskii-Kosterlitz-Thouless transition at finite temperature $T_{c1}$. For $L_z=4$, as temperature is further lowered, a long-range order breaking the $Z_6$ symmetry develops at a second transition at $T_{c2} < T_{c1}$. The transition at $T_{c2}$ is also Berezinskii-Kosterlitz-Thouless-like, but has magnetic critical exponent $\eta=1/9$ instead of the conventional value $\eta = 1/4$. The emergent $U(1)$ symmetry is clearly demonstrated in the quasi-long-range ordered region $T_{c2} \leq T \leq T_{c1}$.


title: Geometric mean method for judgement matrices: Formulas for errors

abstract: The geometric mean method (GMM) and the eigenvector method (EM) are well-known approaches to deriving information from pairwise comparison matrices in decision making processes. However, the original algorithms of these methods are logically incomplete and have significant drawbacks: their actual numerical errors are unknown and their reliability is doubted by different rank reversal phenomena including the GMM-EM ones. Recently (Tomashevskii 2015: Eur J Oper Res 240:774-780), the actual numerical errors were found for EM, and it was shown that all EM rank reversal phenomena have the same cause and are eliminated when the errors are taken into account. In this paper the similar ap-proach is used for GMM: we associate GMM with some standard measuring procedure, ana-lyze this procedure from the viewpoint of measurement theory, and find the actual GMM errors. We show that the GMM and the GMM-EM rank reversal phenomena are eliminated when the GMM and the EM errors are taken into account. The GMM decision support tool, which has all components of a standard measuring tool, is composed of pairwise comparisons as an initial measuring procedure, GMM as a data processor, and the obtained formulas for GMM errors as an error indicator. This tool is analogous to the EM decision support tool received earlier in the above paper. It is shown that the EM and the GMM versions of the tool are equally suitable to measure and rank any comparable elements with positive numerical values. We also analyze the Geometric Consistency Index usually used in the original GMM to measure of the inconsistency and to accept or reject an inconsistent pairwise comparison matrix, and show that this index is not an acceptable GMM error indicator.


title: The Sequential Empirical Process of a Random Walk in Random Scenery

abstract: A random walk in random scenery $(Y_n)_{n\in\mathbb{N}}$ is given by $Y_n=\xi_{S_n}$ for a random walk $(S_n)_{n\in\mathbb{N}}$ and iid random variables $(\xi_n)_{n\in\mathbb{Z}}$. In this paper, we will show the weak convergence of the sequential empirical process, i.e. the centered and rescaled empirical distribution function. The limit process shows a new type of behavior, combining properties of the limit in the independent case (roughness of the paths) and in the long range dependent case (self-similarity).


title: Bifurcation to Chaos in the complex Ginzburg-Landau equation with large third-order dispersion

abstract: We give an analytic proof of the existence of Shilnikov chaos in complex Ginzburg-Landau equation subject to a large third-order dispersion perturbation.


title: Analytical Evaluation of Saturation Throughput of a Cognitive 802.11-based WLAN Overlaid on a WiMAX-TDD Network

abstract: This paper analyzes the saturation throughput of a cognitive single hop WLAN overlaid on a primary IEEE 802.16e TDD WiMAX network. After the contention among the secondary nodes, the winner node transmits its data packet in the empty slots of downlink subframes of WiMAX. Regarding the OFDMA structure as well as time-scheduled resources in WiMAX, the time duration of opportunities for the secondary network does not follow simple exponential on-off pattern. To model the dynamic behavior of opportunities for secondary nodes as well as contentions to exploit the opportunities, we propose an analytical model comprised of a discrete-time Markov chain and two inter-related open multi-class queueing networks. The effects of random number of empty slots at different frames as the result of random amount of download data, random packet transmission time at WLAN due to random opportunities in different frames, the dependency of the number of empty slots at consecutive WiMAX frames, and the details of 802.11 MAC protocol are included in our analytical approach. We compare the effect of two resource allocations, i.e., horizontal and vertical striping on the saturation throughput of the cognitive WLAN. Simulation results confirm the accuracy of our analytical approach in different conditions.


title: Multiscale Bernstein polynomials for densities

abstract: Our focus is on constructing a multiscale nonparametric prior for densities. The Bayes density estimation literature is dominated by single scale methods, with the exception of Polya trees, which favor overly-spiky densities even when the truth is smooth. We propose a multiscale Bernstein polynomial family of priors, which produce smooth realizations that do not rely on hard partitioning of the support. At each level in an infinitely-deep binary tree, we place a beta dictionary density; within a scale the densities are equivalent to Bernstein polynomials. Using a stick-breaking characterization, stochastically decreasing weights are allocated to the finer scale dictionary elements. A slice sampler is used for posterior computation, and properties are described. The method characterizes densities with locally-varying smoothness, and can produce a sequence of coarse to fine density estimates. An extension for Bayesian testing of group differences is introduced and applied to DNA methylation array data.


title: Dynamics of Users Activity on Web-Blogs

abstract: Activity of users on Internet discussion forums is analyzed. The rank of users is shown to be approximated better by stretched-exponential function than by Zipfs law. Cumulative distribution function is found as an excellent tool in analysis of the dynamics of the collective social phenomena. We are able to approximate the number of blog comments with time by simple functions that resemble Fermi-Dirac distribution function: probability of posting a comment is given by $P(t)=P_0(t) / (1+P_0(t)$, where $P_0(t)=exp(a \cdot log(t/t_0))$. It is argued that dynamics of blog entries by a specific user ought to be related to personality of each user.


title: Entropy of foliations with leafwise Finsler structure

abstract: We extend the notion of the geometric entropy of foliation to foliated manifolds equipped with leafwise Finsler structure. We study the relation between the geometric entropy and the topological entropy of the holonomy pseudogroup. The case of foliated manifold with leafwise Randers structure. In this case the estimates for one dimensional foliation defined by a vector field in term of topological entropy of a flow are presented.


title: Enhancement of the thermoelectric power by electronic correlations in bad metals: a study of the Kelvin formula

abstract: In many strongly correlated electron metals the thermoelectric power has a non-monotonic temperature dependence and values that are orders of magnitude larger than for elemental metals. Kelvin proposed a particularly simple expression for the thermopower in terms of the temperature dependence of the chemical potential. We consider a Hubbard model on an anisotropic triangular lattice at half filling, a minimal effective Hamiltonian for several classes of organic charge transfer salts. The finite temperature Lanczos method is used to calculate the temperature dependence of the thermopower using the Kelvin formula. We find that electronic correlations significantly enhance the magnitude of the thermopower and lead to a non-monotonic temperature dependence. The latter reflects a crossover with increasing temperature from a Fermi liquid to a bad metal. Although, the Kelvin formula gives a semi-quantitative description of some experimental results it cannot describe the directional dependence of the sign of the thermopower in some materials.


title: Effect of diffusion in simple discontinuous absorbing transition models

abstract: Discontinuous transitions into absorbing states require an effective mechanism that prevents the stabilization of low density states. They can be found in different systems, such as lattice models or stochastic differential equations (e.g. Langevin equations). Recent results for the latter approach have shown that the inclusion of limited diffusion suppresses discontinuous transitions, whereas they are maintained for larger diffusion strengths. Here we give a further step by addressing the effect of diffusion in two simple lattice models originally presenting discontinuous absorbing transitions. They have been studied via mean-field theory (MFT) and distinct sort of numerical simulations. For both cases, results suggest that the diffusion does not change the order of the transition, regardless its strength and thus, in partial contrast with results obtained from Langevin approach. Also, all transitions present a common finite size scaling behavior that is similar to discontinuous absorbing transitions studied in Phys. Rev. E {\bf 89}, 022104 (2014).


title: Unexpected impact of magnetic disorder on multiband superconductivity

abstract: We analyze how the magnetic disorder affects the properties of the two-band $s_\pm$ and $s_{++}$ models, which are subject of hot discussions regarding iron-based superconductors and other multiband systems like MgB$_2$. We show that there are several cases when the transition temperature $T_c$ is not fully suppressed by magnetic impurities in contrast to the Abrikosov-Gor'kov theory, but a saturation of $T_c$ takes place in the regime of strong disorder. These cases are: (1) the purely interband impurity scattering, (2) the unitary scattering limit. We show that in the former case the $s_\pm$ gap is preserved, while the $s_{++}$ state transforms into the $s_\pm$ state with increasing magnetic disorder. For the case (2), the gap structure remains intact.


title: Improved Algorithms for Parity and Streett objectives

abstract: The computation of the winning set for parity objectives and for Streett objectives in graphs as well as in game graphs are central problems in computer-aided verification, with application to the verification of closed systems with strong fairness conditions, the verification of open systems, checking interface compatibility, well-formedness of specifications, and the synthesis of reactive systems. We show how to compute the winning set on $n$ vertices for (1) parity-3 (aka one-pair Streett) objectives in game graphs in time $O(n^{5/2})$ and for (2) k-pair Streett objectives in graphs in time $O(n^2 + nk \log n)$. For both problems this gives faster algorithms for dense graphs and represents the first improvement in asymptotic running time in 15 years.


title: Re-interpretation of "Bar slowdown and the distribution of dark matter in barred galaxies" by Athanassoula

abstract: Athanassoula (2014) has claimed that measurements of the ratio of corotation radius to bar length in galaxies do not place a constraint on the disk to halo mass ratio. Her conclusion was based on a series of simulations published by Athanassoula et al. (2013). Here we show that these results are, in fact, consistent with previous work on the slow down of bars due to dynamical friction because gas inflow rearranges the disk mass before the bar settles. It therefore remains true that a bar rotating sufficiently fast that corotation is not far beyond the bar end requires a near maximum disk.


title: Subgroup depth and twisted coefficients

abstract: Danz computes the depth of certain twisted group algebra extensions in Comm. Alg. (2011), which are less than the values of the depths of the corresponding untwisted group algebra extensions in Burciu et al, I.E.J.A. (2011). In this paper, we show that the closely related h-depth of any group crossed product algebra extension is less than or equal to the h-depth of the corresponding (finite rank) group algebra extension. A convenient theoretical underpinning to do so is provided by the entwining structure of a right $H$-comodule algebra A and a right H-module coalgebra C for a Hopf algebra H. Then A tensor C is an A-coring, where corings have a notion of depth extending h-depth. This coring is Galois in certain cases where C is the quotient module Q of a coideal subalgebra R < H. We note that this applies for the group crossed product algebra extension, so that the depth of this Galois coring is less than the h-depth of H in G. Along the way, we show that subgroup depth behaves exactly like combinatorial depth with respect to the core of a subgroup, and extend results in Kadison J.Pure Appl.Alg. (2014) to coideal subalgebras of finite dimension.


title: Hydrogen and muonic-Hydrogen Atomic Spectra in Non-commutative Space-Time

abstract: Comparing electronic Hydrogen with muonic Hydrogen shows that the discrepancy in measurement of the Lamb shift in the both systems are relatively of order of $(\frac{m_\mu}{m_e})^{4-5}$. We explore the spectrum of Hydrogen atom in noncommutative $QED$ to compare the noncommutative effects on the both bound states. We show that in the Lorentz violating noncommutative QED the ratio of NC-corrections is $(\frac{m_\mu}{m_e})^3$ while in the Lorentz conserving NCQED is $(\frac{m_\mu}{m_e})^5$. An uncertainty about $1 \,Hz\ll 3\,kHz$ in the Lamb shift of Hydrogen atom leads to an NC correction about $10 \,MHz$ in the Lorentz violating noncommutative QED and about $400 \,GHz$ in the Lorentz conserving noncommutative QED.


title: Asymptotic Representations of Quantum Affine Superalgebras

abstract: We study representations of the quantum affine superalgebra associated with a general linear Lie superalgebra. In the spirit of Hernandez-Jimbo, we construct inductive systems of Kirillov-Reshetikhin modules based on a cyclicity result that we established previously on tensor products of these modules, and realize their inductive limits as modules over its Borel subalgebra, the so-called $q$-Yangian. A new generic asymptotic limit of the same inductive systems is proposed, resulting in modules over the full quantum affine superalgebra. We derive generalized Baxter's relations in the sense of Frenkel-Hernandez for representations of the full quantum group.


title: Effect of charge carrier relaxation during hopping process on electroluminescence in organic solids

abstract: Energetic disorder in disordered organic solids has been found to alter their physical parameters. Here, we have demonstrated, by means of Monte-Carlo simulation and experiments, that the electroluminescence (EL) spectrum is dependent on energetic disorder. This dependence has been attributed to the charge carrier relaxation during hopping process. The dependence of EL spectrum on energetic disorder makes it temperature dependent and temperature dependence has been found to vary with energetic disorder in a variety of materials. The simulation has been performed by taking the relaxation of charge carriers via transport energy in the Gaussian density of states. An analytical equation was established for spectral shift as a function of transport energy.


title: New results on torus cube packings and tilings

abstract: We consider sequential random packing of integral translate of cubes $[0,N]^n$ into the torus $Z^n / 2NZ^n$. Two special cases are of special interest: (i) The case $N=2$ which corresponds to a discrete case of tilings (considered in \cite{cubetiling,book}) (ii) The case $N=\infty$ corresponds to a case of continuous tilings (considered in \cite{combincubepack,book}) Both cases correspond to some special combinatorial structure and we describe here new developments.


title: Multirhythmicity in an optoelectronic oscillator with large delay

abstract: An optoelectronic oscillator exhibiting a large delay in its feedback loop is studied both experimentally and theoretically. We show that multiple square-wave oscillations may coexist for the same values of the parameters (multirhythmicity). Depending on the sign of the phase shift, these regimes admit either periods close to an integer fraction of the delay or periods close to an odd integer fraction of twice the delay. These periodic solutions emerge from successive Hopf bifurcation points and stabilize at a finite amplitude following a scenario similar to Eckhaus instability in spatially extended systems. We find quantitative agreements between experiments and numerical simulations. The linear stability of the square-waves is substantiated analytically by determining stable fixed points of a map.


title: MIARMA: An information preserving method for filling gaps in time series. Application to CoRoT light curves

abstract: The method here presented intends to minimize the effect of the gaps in the power spectra by gap-filling preserving the original information, that is, in the case of asteroseismology, the stellar oscillation frequency content. We make use of a forward-backward predictor based on autoregressive moving average modelling (ARMA) in the time domain. The method MIARMA is particularly suitable for replacing invalid data such as those present in the light curves of the CoRoT satellite due to the pass through the South Atlantic Anomaly, and eventually for the data gathered by the NASA planet hunter Kepler. We select a sample of stars from the ultra-precise photometry collected by the asteroseismic camera on board the CoRoT satellite: the {\delta} Scuti star HD 174966, showing periodic variations of the same order as the CoRoT observational window, the Be star HD 51193, showing longer time variations, and the solar-like HD 49933, with rapid time variations. We showed that in some cases linear interpolations are less reliable to what was believed. In particular: the power spectrum of HD 174966 is clearly aliased when this interpolation is used for filling the gaps; the light curve of HD 51193 presents a much more aliased spectrum than expected for a low frequency harmonic signal; and finally, although the linear interpolation does not affect noticeably the power spectrum of the CoRoT light curve of the solar-like star HD 49933, the ARMA interpolation showed rapid variations previously unidentified that ARMA interprets as a signal. In any case, the ARMA interpolation method provides a cleaner power spectrum, that is, less contaminated by spurious frequencies. In conclusion, MIARMA appears to be a suitable method for filling gaps in the light curves of pulsating stars observed by CoRoT since the method preserves their frequency content, which is a necessary condition for asteroseismic studies.


title: The Gaia-ESO Survey: Extracting diffuse interstellar bands from cool star spectra, DIB-based interstellar medium line-of-sight structures at the kpc scale

abstract: We study how diffuse interstellar bands (DIBs) measured toward distance-distributed target stars can be used to locate dense interstellar (IS) clouds in the Galaxy and probe a line-of-sight (LOS) kinematical structure, a potential useful tool when gaseous absorption lines are saturated or not available in the spectral range. Cool target stars are numerous enough for this purpose. We have devised automated DIB fitting methods appropriate to cool star spectra and multiple IS components. The data is fitted with a combination of a synthetic stellar spectrum, a synthetic telluric transmission, and empirical DIB profiles. In parallel, stellar distances and extinctions are estimated self-consistently by means of a 2D Bayesian method, from spectroscopically-derived stellar parameters and photometric data. We have analyzed Gaia-ESO Survey (GES) and previously recorded spectra that probe between $\sim$ 2 and 10 kpc long LOS in five different regions of the Milky Way. Depending on the observed spectral intervals, we extracted one or more of the following DIBs: $\lambda\lambda$ 6283.8, 6613.6 and 8620.4. For each field, we compared the DIB strengths with the Bayesian distances and extinctions, and the DIB Doppler velocities with the HI emission spectra. For all fields, the DIB strength and the target extinction are well correlated. In case of targets widely distributed in distance, marked steps in DIBs and extinction radial distance profiles match with each other and broadly correspond to the expected locations of spiral arms. For all fields, the DIB velocity structure agrees with HI emission spectra and all detected DIBs correspond to strong NaI lines. This illustrates how DIBs can be used to locate the Galactic interstellar gas and to study its kinematics at the kpc scale.


title: Critical and subcritical fractional problems with vanishing potentials

abstract: We investigate a class of nonlinear nonautonomous scalar field equations with fractional diffusion, critical power nonlinearity and a subcritical term. The involved potentials are allowed for vanishing behavior at infinity. The problem is set on the whole space and compactness issues have to be tackled.


title: Enhanced brightness and photostability of cyanine dyes by supramolecular containment

abstract: Ultrasensitive detection and real-time monitoring of biological processes can benefit significantly from the improved brightness and photostability of the popular organic dyes such as cyanines. Here, using a model cyanine dye, Cy3, we demonstrate that brightness and photostability of the dye is significantly altered when trapped in a molecular container, e.g. cucurbit[n]urils (CB[n]) and cyclodextrins (CD).Through computational modeling, we predicted that Cy3 forms a stable inclusion complex with three different hosts, CB[7], beta-CD, and methyl-beta-CD, which was further confirmed by single-molecule diffusion measurements using fluorescence correlation spectroscopy. The effect of supramolecular encapsulation on Cy3 photophysical properties was found to be highly host-specific. Up to three-fold increase in brightness of Cy3 was observed when the dye was trapped in methyl-beta-CD, due to an increase in both dye absorption and quantum yield. Steady-state and time-resolved spectroscopy of the various complexes revealed that host polarizability and restricted mobility of the dye in the host both contribute to the observed increase in molecular brightness. Furthermore, entrapment of the dye molecule in CDs resulted in a marked increase in dye photostability, whereas the dye degraded more rapidly in CB[7]. These results suggest that the changes in photophysical properties of the dye afforded by supramolecular encapsulation are highly dependent on the host molecule. The reported improvement in brightness and photostability together with the excellent biocompatibility of cyclodextrins makes supramolecular encapsulation a viable strategy for routine dye enhancement.


title: Spacetime condensation in (2+1)-dimensional CDT from a Horava-Lifshitz minisuperspace model

abstract: A spacetime condensation phenomenon underlies the emergence of a macroscopic universe in causal dynamical triangulations, where the time extension of the condensate is strictly smaller than the total time. It has been known for some time that the volumes of spatial slices in the bulk of the macroscopic universe follow a time evolution which resembles that of a sphere, and their effective dynamics is well described by a minisuperspace reduction of the general relativistic action. More recently, it has been suggested that the same minusuperspace model can also provide an understanding of the condensation phenomenon itself, thus explaining the presence of an extended droplet of spacetime connected to a stalk of minimal spatial extension. We show here that a minisuperspace model based on the general relativistic action fails in that respect for the (2+1)-dimensional case, while a successful condensation is obtained from a minisuperspace model of Horava-Lifshitz gravity.


title: An introduction to Docker for reproducible research, with examples from the R environment

abstract: As computational work becomes more and more integral to many aspects of scientific research, computational reproducibility has become an issue of increasing importance to computer systems researchers and domain scientists alike. Though computational reproducibility seems more straight forward than replicating physical experiments, the complex and rapidly changing nature of computer environments makes being able to reproduce and extend such work a serious challenge. In this paper, I explore common reasons that code developed for one research project cannot be successfully executed or extended by subsequent researchers. I review current approaches to these issues, including virtual machines and workflow systems, and their limitations. I then examine how the popular emerging technology Docker combines several areas from systems research - such as operating system virtualization, cross-platform portability, modular re-usable elements, versioning, and a `DevOps' philosophy, to address these challenges. I illustrate this with several examples of Docker use with a focus on the R statistical environment.


title: The low Sr/Ba ratio on some extremely metal-poor stars

abstract: It has been noted that, in classical extremely metal-poor (EMP) stars, the abundance ratio of Sr and Ba, is always higher than [Sr/Ba] = -0.5, the value of the solar r-only process; however, a handful of EMP stars have recently been found with a very low Sr/Ba ratio. We try to understand the origin of this anomaly by comparing the abundance pattern of the elements in these stars and in the classical EMP stars. Four stars with very low Sr/Ba ratios were observed and analyzed within LTE approximation through 1D (hydrostatic) model atmosphere, providing homogeneous abundances of nine neutron-capture elements. In CS 22950-173, the only turnoff star of the sample, the Sr/Ba ratio is, in fact, found to be higher than the r-only solar ratio, so the star is discarded. The remaining stars (CS 29493-090, CS 30322-023, HE 305-4520) are cool evolved giants. They do not present a clear carbon enrichment. The abundance patterns of the neutron-capture elements in the three stars are strikingly similar to a theoretical s-process pattern. This pattern could at first be attributed to pollution by a nearby AGB, but none of the stars presents a clear variation in the radial velocity indicating the presence of a companion. The stellar parameters seem to exclude any internal pollution in a TP-AGB phase for at least two of these stars. The possibility that the stars are early-AGB stars polluted during the core He flash does not seem compatible with the theory.


title: Determination of normalized magnetic eigenfields in microwave cavities

abstract: The magnetic field integral equation for axially symmetric cavities with perfectly conducting surfaces is discretized according to a high-order convergent Fourier--Nystr\"om scheme. The resulting solver is used to determine eigenwavenumbers and normalized magnetic eigenfields to very high accuracy in the entire computational domain.


title: Braidlab: A Software Package for Braids and Loops

abstract: Braidlab is a Matlab package for analyzing data using braids. It was designed to be fast, so it can be used on relatively large problems. It uses the object-oriented features of Matlab to provide a class for braids on punctured disks and a class for equivalence classes of simple closed loops. The growth of loops under iterated action by braids is used to compute the topological entropy of braids, as well as for determining the equality of braids. This guide is a survey of the main capabilities of braidlab, with many examples; the help messages of the various commands provide more details. Some of the examples contain novel observations, such as the existence of cycles of the linear effective action for arbitrary braids.


title: Recurrence plots from altimetry data of some lakes in Africa

abstract: The paper shows recurrence plots obtained from time series of the level variations of four lakes in Africa (Nasser, Tana, Chad and Kainji). The data, coming from remote sensing, are provided by the United States Department of Agriculture. The recurrence plots allow a good visual comparison of the behaviours of local drainage basins.


title: Effects of the symmetry energy on properties of neutron star crusts near the neutron drip density

abstract: We study the effects of the symmetry energy on the neutron drip density and properties of nuclei in neutron star crusts. The nonuniform matter around the neutron drip point is calculated by using the Thomas--Fermi approximation with the relativistic mean-field model. The neutron drip density and the composition of the crust are found to be correlated with the symmetry energy and its slope. We compare the self-consistent Thomas--Fermi approximation with other treatments of surface and Coulomb energies and find that these finite-size effects play an essential role in determining the equilibrium state at low density.


title: A General Duality Relation with Applications in Quantitative Risk Management

abstract: A fundamental problem in risk management is the robust aggregation of different sources of risk in a situation where little or no data are available to infer information about their dependencies. A popular approach to solving this problem is to formulate an optimization problem under which one maximizes a risk measure over all multivariate distributions that are consistent with the available data. In several special cases of such models, there exist dual problems that are easier to solve or approximate, yielding robust bounds on the aggregated risk. In this chapter we formulate a general optimization problem, which can be seen as a doubly infinite linear programming problem, and we show that the associated dual generalizes several well known special cases and extends to new risk management models we propose.


title: Production, characterization and operation of $^{76}$Ge enriched BEGe detectors in GERDA

abstract: The GERmanium Detector Array (GERDA) at the Gran Sasso Underground Laboratory (LNGS) searches for the neutrinoless double beta decay (0{\nu}{\beta}{\beta}) of $^{76}$Ge. Germanium detectors made of material with an enriched $^{76}$Ge fraction act simultaneously as sources and detectors for this decay. During Phase I of the experiment mainly refurbished semi-coaxial Ge detectors from former experiments were used. For the upcoming Phase II, 30 new $^{76}$Ge enriched detectors of broad energy germanium (BEGe)-type were produced. A subgroup of these detectors has already been deployed in GERDA during Phase I. The present paper reviews the complete production chain of these BEGe detectors including isotopic enrichment, purification, crystal growth and diode production. The efforts in optimizing the mass yield and in minimizing the exposure of the $^{76}$Ge enriched germanium to cosmic radiation during processing are described. Furthermore, characterization measurements in vacuum cryostats of the first subgroup of seven BEGe detectors and their long-term behavior in liquid argon are discussed. The detector performance fulfills the requirements needed for the physics goals of GERDA Phase~II.


title: On the Possible Detection of Low Frequency Periodic Signals in Gravitational Wave Interferometers

abstract: We carried out a computer simulation of a large gravitational wave (GW) interferometer using the specifications of the LIGO instruments. We find that if in addition to the carrier, a single sideband offset from the carrier by the fsr frequency (the free spectral range of the arm cavities) is injected, it is equally sensitive to GW signals as is the carrier. The amplitude of the fsr sideband signal in the DC region is generally much less subject to noise than the carrier, and this makes possible the detection of periodic signals with frequencies well below the so-called seismic wall.