title: Statistical Inference for Model Parameters in Stochastic Gradient Descent

abstract: The stochastic gradient descent (SGD) algorithm has been widely used in statistical estimation for large-scale data due to its computational and memory efficiency. While most existing works focus on the convergence of the objective function or the error of the obtained solution, we investigate the problem of statistical inference of true model parameters based on SGD when the population loss function is strongly convex and satisfies certain smoothness conditions. Our main contributions are two-fold. First, in the fixed dimension setup, we propose two consistent estimators of the asymptotic covariance of the average iterate from SGD: (1) a plug-in estimator, and (2) a batch-means estimator, which is computationally more efficient and only uses the iterates from SGD. Both proposed estimators allow us to construct asymptotically exact confidence intervals and hypothesis tests. Second, for high-dimensional linear regression, using a variant of the SGD algorithm, we construct a debiased estimator of each regression coefficient that is asymptotically normal. This gives a one-pass algorithm for computing both the sparse regression coefficients and confidence intervals, which is computationally attractive and applicable to online data.


title: Path to Low Cost Microfluidics

abstract: The paper describes a novel concept for a low cost microfluidic platform utilizing materials and processes used in low cost thermal inkjet printing. The concept re-purposes the jetting elements to create pumps, mixers, and valves all necessary components for the transport of fluids in a broad range of microfluidic applications.


title: A Simplified Model for Fast Optimization of Free-Electron Laser Oscillator

abstract: A simplified theoretical model for free-electron laser oscillator (FELO) simulation which reserves the main physics is proposed. In stead of using traditional macro particles sampling method, the theoretical model takes advantages of low gain theory to calculate the optical power single-pass gain in the undulator analytically, and some reasonable approximations are made to simplify the calculation of power growth in the cavity. The theoretical analysis of single-pass gain, power growth, time-dependent laser profile evolution and cavity desynchronism are accomplished more efficiently. We present the results of infrared wavelength FELO and X-ray FELO with the new model. The results is checked by simulation with GENESIS and OPC which demonstrates the validity of the theoretical model.


title: Anomaly Detection with the Voronoi Diagram Evolutionary Algorithm

abstract: This paper presents the Voronoi diagram-based evolutionary algorithm (VorEAl). VorEAl partitions input space in abnormal/normal subsets using Voronoi diagrams. Diagrams are evolved using a multi-objective bio-inspired approach in order to conjointly optimize classification metrics while also being able to represent areas of the data space that are not present in the training dataset. As part of the paper VorEAl is experimentally validated and contrasted with similar approaches.


title: Upper critical field in the model with finite range interaction between electrons

abstract: We develop a theory of the upper critical field in a BCS superconductor with a nonlocal interaction between electrons. We have shown that the nonlocal interaction is characterized by the parameter $k_F\rho_0$ where $k_F$ is the Fermi momentum and $\rho_0$ is the radius of electron-electron interaction. The presence of the external magnetic field leads to the generation of additional components of the order parameter with different angular momentum. This effect leads to the enhancement of the upper critical field above the orbital limiting field. In addition the upward curvature in the temperature dependence of $H_{c2}(T)$ in the clean limit is predicted. The impurity scattering suppresses the effect in the dirty limit.


title: Patterning Superconductivity in a Topological Insulator

abstract: While topological superconductors are predicted to provide building blocks for fault-tolerant quantum computing, one of the remaining challenges is to find a convenient experimental platform that would allow patterning of circuits. We find that superconductivity can be patterned directly into Bi$_2$Se$_3$ nanostructures by selective doping with palladium (Pd). Superconducting regions are defined by depositing Pd on top of the nanostructures using electron beam lithography, followed by in-situ annealing. Electrical transport measurements at low temperatures show either partial or full superconducting transition, depending on the doping conditions. Structural characterization techniques indicate that Pd remains localized in the targeted areas, making it possible to pattern superconducting circuits of arbitrary shapes in this topological material.


title: Corrected knife-edge-based reconstruction of tightly focused higher order beams

abstract: The knife-edge method is an established technique for profiling of even tightly focused light beams. However the straightforward implementation of this method fails if the materials and geometry of the knife-edges are not chosen carefully or in particular if knife-edges are used that are made of pure materials. In these cases artifacts are introduced in the shape and position of the reconstructed beam profile due to the interaction of the light beam under study with the knife. Hence, corrections to the standard knife-edge evaluation method are required. Here we investigate the knife-edge method for highly focused radially and azimuthally polarized beams and their linearly polarized constituents. We introduce relative shifts for those constituents and report on the consistency with the case of a linearly polarized Gaussian beam. An adapted knife-edge reconstruction technique is presented and proof-of-concept tests demonstrating the reconstruction of beam profiles are shown.


title: Utility Maximization and Indifference Value under Risk and Information Constraints for a Market with a Change Point

abstract: In this article we consider an optimization problem of expected utility maximization of continuous-time trading in a financial market. This trading is constrained by a benchmark for a utility-based shortfall risk measure. The market consists of one asset whose price process is modeled by a Geometric Brownian motion where the market parameters change at a random time. The information flow is modeled by initially and progressively enlarged filtrations which represent the knowledge about the price process, the Brownian motion and the random time. We solve the maximization problem and give the optimal terminal wealth depending on these different filtrations for general utility functions by using martingale representation results for the corresponding filtration. Moreover, for a special utility function and risk measure we calculate the utility indifference value which measures the gain of further information for the investor.


title: Three-dimensional anomalous twisted gauge theories with global symmetry: Implications for quantum spin liquids

abstract: Topological spin liquids can be described by topological gauge theories with global symmetry. Due to the presence of both nontrivial bulk deconfined gauge fluxes and global symmetry, topological spin liquids are examples of the so-called "Symmetry Enriched Topological phases" (SET). In this paper, we find that, in some twisted versions of topological gauge theories (with discrete Abelian gauge group $G_g$), implementing a global symmetry (denoted by $G_s$) is anomalous although symmetry charge carried by topological point-like excitations is normally fractionalized and classified by the second cohomology group. To demonstrate the anomaly, we fully gauge the global symmetry, rendering a new gauge theory that is not gauge invariant. Therefore, the SET order of the ground state is anomalous, which cannot exist in 3D system alone. Such anomalous state construction generalizes the "2D surface topological order" to 3D. A concrete example with $G_g=\mathbb{Z}_2\times\mathbb{Z}_4$ and $G_s=\mathbb{Z}_2$ is calculated.


title: On the solvability of the Brinkman-Forchheimer-extended Darcy equation

abstract: The nonlinear Brinkman-Forchheimer-extended Darcy equation is used to model some porous medium flow in chemical reactors of packed bed type. The results concerning the existence and uniqueness of a weak solution are presented for nonlinear convective flows in medium with nonconstant porosity and for small data. Furthermore, the finite element approximations to the flow profiles in the fixed bed reactor are presented for several Reynolds numbers at the non-Darcy's range.


title: Spectral-Galerkin Approximation and Optimal Error Estimate for Stokes Eigenvalue Problems in Polar Geometries

abstract: In this paper we propose and analyze spectral-Galerkin methods for the Stokes eigenvalue problem based on the stream function formulation in polar geometries. We first analyze the stream function} formulated fourth-order equation under the polar coordinates, then we derive the pole condition and reduce the problem on a circular disk to a sequence of equivalent one-dimensional eigenvalue problems that can be solved in parallel. The novelty of our approach lies in the construction} of suitably weighted Sobolev spaces according to the pole conditions, based on which, the optimal error estimate for approximated eigenvalue of each one dimensional problem can be obtained. Further, we extend our method to the non-separable Stokes eigenvalue problem in an elliptic domain and establish the optimal error bounds. Finally, we provide some numerical experiments to validate our theoretical results and algorithms.


title: Optimality certificates for convex minimization and Helly numbers

abstract: We consider the problem of minimizing a convex function over a subset of R^n that is not necessarily convex (minimization of a convex function over the integer points in a polytope is a special case). We define a family of duals for this problem and show that, under some natural conditions, strong duality holds for a dual problem in this family that is more restrictive than previously considered duals.


title: Minkowski Endomorphisms

abstract: Several open problems concerning Minkowski endomorphisms and Minkowski valuations are solved. More precisely, it is proved that all Minkowski endomorphisms are uniformly continuous but that there exist Minkowski endomorphisms that are not weakly-monotone. This answers questions posed repeatedly by various authors. Furthermore, a recent representation result for Minkowski valuations by Schuster and Wannerer is improved under additional homogeneity assumptions. Also a question related to the structure of Minkowski endomorphisms by the same authors is answered. Finally, it is shown that there exists no McMullen decomposition in the class of continuous, even, $SO(n)$-equivariant and translation invariant Minkowski valuations extending a result by Parapatits and Wannerer.


title: Electrostatics of quadrupolarizable media

abstract: The classical macroscopic Maxwell equations are approximated. They are a corollary of the multipole expansion of the local electrostatic potential up to dipolar terms. But quadrupolarization of the medium should not be neglected if the molecules which build up the medium possess large quadrupole moment or do not have any dipole moment. If we include the quadrupolar terms in Maxwell equations we obtain the quadrupolar analogue of Poisson's equation: $\nabla^2 \phi - L^2_Q\nabla^4 \phi = - \rho / \varepsilon$. This equation is of the fourth order and it requires not only the two classical boundary conditions but also two additional ones: continuous electric field and the relation of the jump of the normal quadrupolarizability at the surface to the intrinsic normal surface dipole moment. The account of the quadrupole moment of the molecules leads to significant differences compared to the classical electrostatic theory.


title: On holomorphic Artin L-functions

abstract: Let $K/\mathbb Q$ be a finite Galois extension, $s_0\in \mathbb C\setminus \{1\}$, ${\it Hol}(s_0)$ the semigroup of Artin L-functions holomorphic at $s_0$. If the Galois group is almost monomial then Artin's L-functions are holomorphic at $s_0$ if and only if $ {\it Hol}(s_0)$ is factorial. This holds also if $s_0$ is a zero of an irreducible L-function of dimension $\leq 2$, without any condition on the Galois group.


title: Scale-invariant radio jets and varying black hole spin

abstract: Compact radio cores associated with relativistic jets are often observed in both active galactic nuclei and X-ray binaries. Their radiative properties follow some general scaling laws which primarily depend on their masses and accretion rates. However, it has been suggested that the black hole spin can also strongly influence the power and radio flux of these. Here, we attempt to estimate the dependency of the radio luminosity of steady jets launched by accretion disks on black hole mass, accretion rate and spin using numerical simulations. We make use of 3D GRMHD simulations of accretion disks around low-luminosity black holes in which the jet radio emission is produced by the jet sheath. We find that the radio flux increases roughly by a factor of 6 as the back hole spin increases from a~0 to a=0.98. This is comparable to the increase in accretion power with spin, meaning that the ratio between radio jet and accretion power is hardly changing. Although our jet spine power scales as expected for the Blandford-Znajek process, the dependency of jet radio luminosity on the black hole spin is somewhat weaker. Also weakly rotating black holes can produce visible radio jets. The overall scaling of the radio emission with black hole mass and accretion rate is consistent with the scale-invariant analytical models used to explain the fundamental plane of black hole activity. Spin does not introduce a significant scatter in this model. The jet-sheath model can describe well resolved accreting systems, such as SgrA* and M87, as well as the general scaling behavior of low-luminosity black holes. Hence the model should be applicable to a wide range of radio jets in sub-Eddington black holes. The black hole spin has an effect on the production of visible radio jet, but it may not be the main driver to produce visible radio jets. An extension of our findings to powerful quasars remains speculative.


title: Detection of Very Low-Frequency Quasi-Periodic Oscillations in the 2015 Outburst of V404 Cygni

abstract: In June 2015, the black hole X-ray binary (BHXRB) V404 Cygni went into outburst for the first time since 1989. Here, we present a comprehensive search for quasi-periodic oscillations (QPOs) of V404 Cygni during its recent outburst, utilizing data from six instruments on board five different X-ray missions: Swift/XRT, Fermi/GBM, Chandra/ACIS, INTEGRAL's IBIS/ISGRI and JEM-X, and NuSTAR. We report the detection of a QPO at 18 mHz simultaneously with both Fermi/GBM and Swift/XRT, another example of a rare but slowly growing new class of mHz-QPOs in BHXRBs linked to sources with a high orbital inclination. Additionally, we find a duo of QPOs in a Chandra/ACIS observation at 73 mHz and 1.03 Hz, as well as a QPO at 136 mHz in a single Swift/XRT observation that can be interpreted as standard Type-C QPOs. Aside from the detected QPOs, there is significant structure in the broadband power, with a strong feature observable in the Chandra observations between 0.1 and 1 Hz. We discuss our results in the context of current models for QPO formation.


title: A Four-Body Convex Central Configuration with Perpendicular Diagonals Is Necessarily a Kite

abstract: We prove that any four-body convex central configuration with perpendicular diagonals must be a kite configuration. The result extends to general power-law potential functions, including the planar four-vortex problem.


title: Pacific Neutrinos: Towards a High Precision Measurement of CP Violation ?

abstract: The application of deep sea low energy neutrino detection techniques to long baseline neutrino physics is investigated, with a focus on a possible configuration based on a FNAL neutrino beam impinging a detector hosted by the NEPTUNE/OOI submarine observatories offshore of British Columbia.


title: Coherence depletion in the Grover quantum search algorithm

abstract: We investigate the role of quantum coherence depletion (QCD) in Grover search algorithm (GA) by using several typical measures of quantum coherence and quantum correlations. By using the relative entropy of coherence measure ($\mathcal{C}_r$), we show that the success probability depends on the QCD. The same phenomenon is also found by using the $l_1$ norm of coherence measure ($\mathcal{C}_{l_1}$). In the limit case, the cost performance is defined to characterize the behavior about QCD in enhancing the success probability of GA, which is only related to the number of searcher items and the scale of database, no matter using $\mathcal{C}_r$ or $\mathcal{C}_{l_1}$. In generalized Grover search algorithm (GGA), the QCD for a class of states increases with the required optimal measurement time. In comparison, the quantification of other quantum correlations in GA, such as pairwise entanglement, multipartite entanglement, pairwise discord and genuine multipartite discord, cannot be directly related to the success probability or the optimal measurement time. Additionally, we do not detect pairwise nonlocality or genuine tripartite nonlocality in GA since Clauser-Horne-Shimony-Holt inequality and Svetlichny's inequality are not violated.


title: Critical properties of the eight-vertex model in a field

abstract: The general eight-vertex model on a square lattice is studied numerically by using the Corner Transfer Matrix Renormalization Group method. The method is tested on the symmetric (zero-field) version of the model, the obtained dependence of critical exponents on model's parameters is in agreement with Baxter's exact solution and weak universality is verified with a high accuracy. It was suggested longtime ago that the symmetric eight-vertex model is a special exceptional case and in the presence of external fields the eight-vertex model falls into the Ising universality class. We confirm numerically this conjecture in a subspace of vertex weights, except for two specific combinations of vertical and horizontal fields for which the system still exhibits weak universality.


title: Extended Particles and the Exterior Calculus

abstract: These notes were delivered as a series of NIMROD lectures at the Rutherford Appleton Laboratory by the author in February 1976 (RL-76-022). The purpose of these lectures was primarily two-fold: to discuss the classical theory of free point particles, free strings and free membranes from a unified viewpoint; and to present in the process of doing this the rudiments of an intrinsic geometrical calculus that the author has found of immense value in investigating these systems. It is shown how the equations of motion for such classically extended relativistic systems arise in a very simple manner from a principle of stationary action and furthermore how the boundary conditions for finite systems may be derived in a gauge invariant way. Momenta are naturally introduced and the primary constraints that exist in a Hamiltonian description follow simply. Calculations may proceed in an index-free manner until components are required. It is at this stage that one can, if one desires, impose gauge conditions and remove non-independent degrees of freedom. Such methods can be applied in any spacetime of any dimension and metric, and examples are given throughout.


title: A rational use of BCA code MARLOWE for ballistic effects of ion beam irradiation in the ion mixing formalism: comparison to Molecular Dynamics

abstract: Understanding ballistic effects caused by ion beam irradiation, and linking them with induced structure can be a key point for controlling and predicting the microstructure of irradiated materials. For this purpose, we have investigated ballistic effects from an ion mixing formalism point of view. The displacement cascades in copper and AgCu alloy were obtained using binary collision approximation (BCA) and molecular dynamics (MD) simulations. We employed BCA-based code MARLOWE for its ability to simulate high energy displacement cascades. A first set of simulations was performed using both methods on pure copper for energies ranging from 0.5 keV to 20 keV. The results of BCA and MD simulations are compared, evidencing rationally parametrized MARLOWE to be predictive. A second set of simulations was then carried out using BCA only. Following experimental studies, AgCu alloy was subjected to 1 MeV krypton ions. MARLOW simulations are found to be in good agreement with experimental results.


title: Spin transfer driven resonant expulsion of a magnetic vortex core for efficient rf detector

abstract: Spin transfer magnetization dynamics have led to considerable advances in Spintronics, including opportunities for new nanoscale radiofrequency devices. Among the new functionalities is the radiofrequency(rf) detection using the spin diode rectification effect in spin torque nano-oscillators (STNOs). In this study, we focus on a new phenomenon, the resonant expulsion of a magnetic vortex in STNOs. This effect is observed when the excitation vortex radius, due to spin torques associated to rf currents, becomes larger than the actual radius of the STNO. This vortex expulsion is leading to a sharp variation of the voltage at the resonant frequency. Here we show that the detected frequency can be tuned by different parameters; furthermore, a simultaneous detection of different rf signals can be achieved by real time measurements with several STNOs having different diameters. This result constitutes a first proof-of-principle towards the development of a new kind of nanoscale rf threshold detector.


title: Adiabatic physics of an exchange-coupled spin-dimer system: magnetocaloric effect, zero-point fluctuations, and possible two-dimensional universal behavior

abstract: We present the magnetic and thermal properties of the bosonic-superfluid phase in a spin-dimer network using both quasistatic and rapidly-changing pulsed magnetic fields. The entropy derived from a heat-capacity study reveals that the pulsed-field measurements are strongly adiabatic in nature and are responsible for the onset of a significant magnetocaloric effect (MCE). In contrast to previous predictions we show that the MCE is not just confined to the critical regions, but occurs for all fields greater than zero at sufficiently low temperatures. We explain the MCE using a model of the thermal occupation of exchange-coupled dimer spin-states and highlight that failure to take this effect into account inevitably leads to incorrect interpretations of experimental results. In addition, the heat capacity in our material is suggestive of an extraordinary contribution from zero-point fluctuations and appears to indicate universal behavior with different critical exponents at the two field-induced critical points. The data are consistent with a two-dimensional nature of spin excitations in the system.


title: Functional limit theorems for the number of busy servers in a $G/G/\infty$ queue

abstract: We discuss weak convergence of the number of busy servers in a $G/G/\infty$ queue in the $J_1$-topology on the Skorokhod space. We prove two functional limit theorems, with random and nonrandom centering, respectively, thereby solving two open problems stated in Mikosch and Resnick (2006}. A new integral representation for the limit Gaussian process is given.


title: Regularization parameter selection in indirect regression by residual based bootstrap

abstract: Residual-based analysis is generally considered a cornerstone of statistical methodology. For a special case of indirect regression, we investigate the residual-based empirical distribution function and provide a uniform expansion of this estimator, which is also shown to be asymptotically most precise. This investigation naturally leads to a completely data-driven technique for selecting a regularization parameter used in our indirect regression function estimator. The resulting methodology is based on a smooth bootstrap of the model residuals. A simulation study demonstrates the effectiveness of our approach.


title: A random version of principal component analysis in data clustering

abstract: Principal component analysis (PCA) is a widespread technique for data analysis that relies on the covariance-correlation matrix of the analyzed data. However to properly work with high-dimensional data, PCA poses severe mathematical constraints on the minimum number of different replicates or samples that must be included in the analysis. Here we show that a modified algorithm works not only on well dimensioned datasets, but also on degenerated ones.


title: Direct semi-parametric estimation of the state price density implied in option prices

abstract: We present a model for direct semi-parametric estimation of the State Price Density (SPD) implied in quoted option prices. We treat the observed prices as expected values of possible pay-offs at maturity, weighted by the unknown probability density function. We model the logarithm of the latter as a smooth function while matching the expected values of the potential pay-offs with the observed prices. This leads to a special case of the penalized composite link model. Our estimates do not rely on any parametric assumption on the underlying asset price dynamics and are consistent with no-arbitrage conditions. The model shows excellent performance in simulations and in application to real data.


title: Performance Assessment of Variational Integrators for Thermomechanical Problems

abstract: Structure-preserving integrators are in the focus of ongoing research because of their distinguished features of robustness and long time stability. In particular, their formulation for coupled problems that include dissipative mechanisms is still an active topic. Conservative formulations, such as the thermo-elastic case without heat conduction, fit well into a variational framework and have been solved with variational integrators, whereas the inclusions of viscosity and heat transfer are still under investigation. To encompass viscous forces and heat transfer, an extension of Hamilton's principle is required. In this contribution we derive variational integrators for thermo-viscoelastic systems with classical heat transfer. Their results are compared for two discrete model problems vs. Energy-Entropy-Momentum methods. Such comparisons allow to draw conclusions about their relative performance, weaknesses and strengths.


title: Heat flux of a granular gas with homogeneous temperature

abstract: A steady state of a granular gas with homogeneous granular temperature, no mass flow, and nonzero heat flux is studied. The state is created by applying an external position--dependent force or by enclosing the grains inside a curved two--dimensional silo. At a macroscopic level, the state is identified with one solution to the inelastic Navier--Stokes equations, due to the coupling between the heat flux induced by the density gradient and the external force. On the contrary, at the mesoscopic level, by exactly solving a BGK model or the inelastic Boltzmann equation in an approximate way, a one--parametric family of solutions is found. Molecular dynamics simulations of the system in the quasi--elastic limit are in agreement with the theoretical results.


title: Real-time optimal control via Deep Neural Networks: study on landing problems

abstract: Recent research on deep learning, a set of machine learning techniques able to learn deep architectures, has shown how robotic perception and action greatly benefits from these techniques. In terms of spacecraft navigation and control system, this suggests that deep architectures may be considered now to drive all or part of the on-board decision making system. In this paper this claim is investigated in more detail training deep artificial neural networks to represent the optimal control action during a pinpoint landing, assuming perfect state information. It is found to be possible to train deep networks for this purpose and that the resulting landings, driven by the trained networks, are close to simulated optimal ones. These results allow for the design of an on-board real time optimal control system able to cope with large sets of possible initial states while still producing an optimal response.


title: Experimental Determination of Dynamical Lee-Yang Zeros

abstract: Statistical physics provides the concepts and methods to explain the phase behavior of interacting many-body systems. Investigations of Lee-Yang zeros --- complex singularities of the free energy in systems of finite size --- have led to a unified understanding of equilibrium phase transitions. The ideas of Lee and Yang, however, are not restricted to equilibrium phenomena. Recently, Lee-Yang zeros have been used to characterize non-equilibrium processes such as dynamical phase transitions in quantum systems after a quench or dynamic order-disorder transitions in glasses. Here, we experimentally realize a scheme for determining Lee-Yang zeros in such non-equilibrium settings. We extract the dynamical Lee-Yang zeros of a stochastic process involving Andreev tunneling between a normal-state island and two superconducting leads from measurements of the dynamical activity along a trajectory. From the short-time behavior of the Lee-Yang zeros, we predict the large-deviation statistics of the activity which is typically difficult to measure. Our method paves the way for further experiments on the statistical mechanics of many-body systems out of equilibrium.


title: Efficient fiber-coupled single-photon source based on quantum dots in a photonic-crystal waveguide

abstract: Many photonic quantum information processing applications would benefit from a high brightness, fiber-coupled source of triggered single photons. Here, we present a fiber-coupled photonic-crystal waveguide single-photon source relying on evanescent coupling of the light field from a tapered out-coupler to an optical fiber. A two-step approach is taken where the performance of the tapered out-coupler is recorded first on an independent device containing an on-chip reflector. Reflection measurements establish that the chip-to-fiber coupling efficiency exceeds 80 %. The detailed characterization of a high-efficiency photonic-crystal waveguide extended with a tapered out-coupling section is then performed. The corresponding overall single-photon source efficiency is 10.9 % $\pm$ 2.3 %, which quantifies the success probability to prepare an exciton in the quantum dot, couple it out as a photon in the waveguide, and subsequently transfer it to the fiber. The applied out-coupling method is robust, stable over time, and broadband over several tens of nanometers, which makes it a highly promising pathway to increase the efficiency and reliability of planar chip-based single-photon sources.


title: Explicit equivalence of quadratic forms over $\mathbb{F}_q(t)$

abstract: We propose a randomized polynomial time algorithm for computing nontrivial zeros of quadratic forms in 4 or more variables over $\mathbb{F}_q(t)$, where $\mathbb{F}_q$ is a finite field of odd characteristic. The algorithm is based on a suitable splitting of the form into two forms and finding a common value they both represent. We make use of an effective formula for the number of fixed degree irreducible polynomials in a given residue class. We apply our algorithms for computing a Witt decomposition of a quadratic form, for computing an explicit isometry between quadratic forms and finding zero divisors in quaternion algebras over quadratic extensions of $\mathbb{F}_q(t)$.


title: Topological aspects of the multi-language phases of the Naming Game on community-based networks

abstract: The Naming Game is an agent-based model where individuals communicate to name an initially unnamed object. On a large class of networks continual pairwise interactions lead the system to an ultimate consensus state, in which agents converge on a globally shared name. Soon after the introduction of the model, it was observed in literature that on community-based networks the path to consensus passes through metastable multi-language states. Subsequently, it was proposed to use this feature as a mean to discover communities in a given network. In this paper we show that metastable states correspond to genuine multi-language phases, emerging in the thermodynamic limit when the fraction of links connecting communities drops below critical thresholds. In particular, we study the transition to multi-language states in the stochastic block model and on networks with community overlap. We also examine the scaling of critical thresholds under variations of topological properties of the network, such as the number and relative size of communities and the structure of intra-/inter-community links. Our results provide a theoretical justification for the proposed use of the model as a community-detection algorithm.


title: A fast solution method for time dependent multidimensional Schr\"odinger equations

abstract: In this paper we propose fast solution methods for the Cauchy problem for the multidimensional Schr\"odinger equation. Our approach is based on the approximation of the data by the basis functions introduced in the theory of approximate approximations. We obtain high-order approximations also in higher dimensions up to a small saturation error, which is negligible in computations, and we prove error estimates in mixed Lebesgue spaces for the inhomogeneous equation. The proposed method is very efficient in high dimensions if the densities allow separated representations. We illustrate the efficiency of the procedure on different examples, up to approximation order 6 and space dimension 200.


title: Seismic diagnosis from gravity modes strongly affected by rotation

abstract: Most of the information we have about the internal rotation of stars comes from modes that are weakly affected by rotation, for example by using rotational splittings. In contrast, we present here a method, based on the asymptotic theory of Prat et al. (2016), which allows us to analyse the signature of rotation where its effect is the most important, that is in low-frequency gravity modes that are strongly affected by rotation. For such modes, we predict two spectral patterns that could be confronted to observed spectra and those computed using fully two-dimensional oscillation codes.


title: The p-radical closure of local noetherian rings

abstract: Given a local noetherian ring $R$ whose formal completion is integral, we introduce and study the $p$-radical closure $R^\text{prc}$. Roughly speaking, this is the largest purely inseparable $R$-subalgebra inside the formal completion $\hat{R}$. It turns out that the finitely generated intermediate rings $R\subset A\subset R^\text{prc}$ have rather peculiar properties. They can be used in a systematic way to provide examples of integral local rings whose normalization is non-finite, that do not admit a resolution of singularities, and whose formal completion is non-reduced.


title: $\kappa$-generalized models of income and wealth distributions: A survey

abstract: The paper provides a survey of results related to the "$\kappa$-generalized distribution", a statistical model for the size distribution of income and wealth. Topics include, among others, discussion of basic analytical properties, interrelations with other statistical distributions as well as aspects that are of special interest in the income distribution field, such as the Gini index and the Lorenz curve. An extension of the basic model that is most able to accommodate the special features of wealth data is also reviewed. The survey of empirical applications given in this paper shows the $\kappa$-generalized models of income and wealth to be in excellent agreement with the observed data in many cases.


title: Simplification of inclusion-exclusion on intersections of unions with application to network systems reliability

abstract: Reliability of safety-critical systems is an important issue in system engineering and in most practical situations the reliability of a non series-parallel network system has to be calculated. Some methods for calculating reliability use the probability principle of inclusion-exclusion. When dealing with complex networks, this leads to very long mathematical expressions which are usually computationally very expensive to calculate. In this paper, we provide a new expression to simplify the probability principle of inclusion-exclusion's formula for intersections of unions, which appear when calculating reliability on non series parallel network systems. This new expression has much less terms, which reduces enormously the computational cost. We also show that the general form of the probability principle of inclusion-exclusion's formula has double exponential complexity whereas the simplified form has only exponential complexity with a linear exponent. Finally, we illustrate how to use this result when calculating the reliability of a door management system in aircraft engineering.


title: $\Lambda$ polarization in peripheral collisions at moderate relativistic energies

abstract: The polarization of $\Lambda$ hyperons from relativistic flow vorticity is studied in peripheral heavy ion reactions at FAIR and NICA energies, just above the threshold of the transition to the Quark-Gluon Plasma. Previous calculations at higher energies with larger initial angular momentum, predicted significant $\Lambda$ polarization based on the classical vorticity term in the polarization, while relativistic modifications decreased the polarization and changed its structure in the momentum space. At the lower energies studied here, we see the same effect namely that the relativistic modifications decrease the polarization arising from the initial shear flow vorticity.


title: Geometric phase of neutrinos: differences between Dirac and Majorana neutrinos

abstract: We analize the non-cyclic geometric phase for neutrinos. We find that the geometric phase and the total phase associated to the mixing phenomenon provide a tool to distinguish between Dirac and Majorana neutrinos. Our results hold for neutrinos propagating in vacuum and through the matter. Future experiments, based on interferometry, could reveal the nature of neutrinos.


title: Weight-dependent commutation relations and combinatorial identities

abstract: We derive combinatorial identities for variables satisfying specific systems of commutation relations, in particular elliptic commutation relations. The identities thus obtained extend corresponding ones for $q$-commuting variables $x$ and $y$ satisfying $yx=qxy$. In particular, we obtain weight-dependent binomial theorems, functional equations for generalized exponential functions, we derive results for an elliptic derivative of elliptic commuting variables, and finally study weight-dependent extensions of the Weyl algebra which we connect to rook theory.


title: New numerical scheme of Atangana-Baleanu fractional integral: an application to groundwater flow within leaky aquifer

abstract: Many physical problems can be described using the integral equations known as Volterra equations. There exists quite a number of analytical methods that can handle these equations for linear cases. For non-linear cases, a numerical scheme is needed to obtain an approximation. Recently a new concept of fractional differentiation was introduced using the Mittag-Leffler function as kernel, and the associate fractional integral was also presented. Up to this point there is no numerical approximation of this new integral in the literature. Therefore, to accommodate researchers working in the field of numerical analysis, we propose in this paper a new numerical scheme for the new fractional integral. To test the accuracy of the new numerical scheme, we first revisit the groundwater model within a leaky aquifer by reverting the time classical derivative with the Atangana-Baleanu fractional derivative in Caputo sense. The new model is solved numerically by using this new scheme.


title: Localization in periodically modulated speckle potentials

abstract: Disorder in a 1D quantum lattice induces Anderson localization of the eigenstates and drastically alters transport properties of the lattice. In the original Anderson model, the addition of a periodic driving increases, in a certain range of the driving's frequency and amplitude, localization length of the appearing Floquet eigenstates. We go beyond the uncorrelated disorder case and address the experimentally relevant situation when spatial correlations are present in the lattice potential. Their presence induces the creation of an effective mobility edge in the energy spectrum of the system. We find that a slow driving leads to resonant hybridization of the Floquet states, by increasing both the participation numbers and effective widths of the states in the strongly localized band and decreasing values of these characteristics for the states in the quasi-extended band. Strong driving homogenizes the bands, so that the Floquet states loose compactness and tend to be spatially smeared. In the basis of the stationary Hamiltonian, these states retain localization in terms of participation number but become de-localized and spectrum-wide in term of their effective widths. Signatures of thermalization are also observed.


title: Isospin-violating dark matter in the light of recent data

abstract: In scenarios where the dark matter interacts differently with protons and neutrons (isospin-violating dark matter), the interpretation of the experimental limits on the dark matter spin-independent cross section may be significantly modified. On the one hand, the direct detection constraints are shifted depending on the target nucleus, possibly changing the hierarchy among different experiments. On the other hand, the relative strength between the bounds from neutrino detectors and those from direct detection experiments is altered, allowing the former to be more competitive. In this paper, the status of isospin-violating dark matter is assessed in the light of recent data, and the prospects for its detection in the near future are analyzed. We find, for example, that there are regions in the parameter space where IceCube currently provides the most stringent limits on the spin-independent cross section, or others where the expected sensitivity of DEAP-3600 is well above the LUX exclusion limit. Our results highlight the complementarity among different targets in direct detection experiments, and between direct detection and neutrino searches in the quest for a dark matter signal.


title: Formal Abel-Jacobi maps

abstract: We realize the infinitesimal Abel-Jacobi map as a morphism of formal deformation theories, realized as a morphism in the homotopy category of differential graded Lie algebras. The whole construction is carried out in a general setting, of which the classical Abel-Jacobi map is a special example.


title: A Framework for Certified Self-Stabilization

abstract: We propose a general framework to build certified proofs of distributed self-stabilizing algorithms with the proof assistant Coq. We first define in Coq the locally shared memory model with composite atomicity, the most commonly used model in the self-stabilizing area. We then validate our framework by certifying a non trivial part of an existing silent self-stabilizing algorithm which builds a $k$-clustering of the network. We also certify a quantitative property related to the output of this algorithm. Precisely, we show that the computed $k$-clustering contains at most $\lfloor \frac{n-1}{k+1} \rfloor + 1$ clusterheads, where $n$ is the number of nodes in the network. To obtain these results, we also developed a library which contains general tools related to potential functions and cardinality of sets.


title: Polarized User and Topic Tracking in Twitter

abstract: Digital traces of conversations in micro-blogging platforms and OSNs provide information about user opinion with a high degree of resolution. These information sources can be exploited to under- stand and monitor collective behaviors. In this work, we focus on polarization classes, i.e., those topics that require the user to side exclusively with one position. The proposed method provides an iterative classification of users and keywords: first, polarized users are identified, then polarized keywords are discovered by monitoring the activities of previously classified users. This method thus allows tracking users and topics over time. We report several experiments conducted on two Twitter datasets during political election time-frames. We measure the user classification accuracy on a golden set of users, and analyze the relevance of the extracted keywords for the ongoing political discussion.


title: A class of quasi-linear Allen-Cahn type equations with dynamic boundary conditions

abstract: In this paper, we consider a class of coupled systems of PDEs, denoted by (ACE)$_{\varepsilon}$ for $ \varepsilon \geq 0 $. For each $ \varepsilon \geq 0 $, the system (ACE)$_{\varepsilon}$ consists of an Allen-Cahn type equation in a bounded spacial domain $ \Omega $, and another Allen-Cahn type equation on the smooth boundary $ \Gamma := \partial \Omega $, and besides, these coupled equations are transmitted via the dynamic boundary conditions. In particular, the equation in $ \Omega $ is derived from the non-smooth energy proposed by Visintin in his monography "Models of phase transitions": hence, the diffusion in $ \Omega $ is provided by a quasilinear form with singularity. The objective of this paper is to build a mathematical method to obtain meaningful $ L^2 $-based solutions to our systems, and to see some robustness of (ACE)$_\varepsilon$ with respect to $ \varepsilon \geq 0 $. On this basis, we will prove two Main Theorems 1 and 2, which will be concerned with the well-posedness of (ACE)$_\varepsilon$ for each $ \varepsilon \geq 0 $, and the continuous dependence of solutions to (ACE)$_\varepsilon$ for the variations of $ \varepsilon \geq 0 $, respectively.


title: A simple model to describe intrinsic stellar noise for exoplanet detection around red giants

abstract: In spite of the huge advances in exoplanet research provided by the NASA Kepler Mission, there remain only a small number of transit detections around evolved stars. Here we present a reformulation of the noise properties of red-giant stars, where the intrinsic stellar granulation, and the stellar oscillations described by asteroseismology play a key role. The new noise model is a significant improvement on the current Kepler results for evolved stars. Our noise model may be used to help understand planet detection thresholds for the ongoing K2 and upcoming TESS missions, and serve as a predictor of stellar noise for these missions. As an application of our noise model, we explore the minimum detectable planet radii for red giant stars, and find that Neptune sized planets should be detectable around low luminosity red giant branch stars.


title: Variational principles and symmetries on fibered multisymplectic manifolds

abstract: The standard techniques of variational calculus are geometrically stated in the ambient of fiber bundles endowed with a (pre)multisymplectic structure. Then, for the corresponding variational equations, conserved quantities (or, what is equivalent, conservation laws), symmetries, Cartan (Noether) symmetries, gauge symmetries and different versions of Noether's theorem are studied in this ambient. In this way, this constitutes a general geometric framework for all these topics that includes, as special cases, first and higher order field theories and (non-autonomous) mechanics.


title: The earliest drawings of datable auroras and a two-tail comet from the Syriac Chronicle of Z\=uqn\=in

abstract: People have probably been watching the sky since the beginning of human history. Observers in pre-telescopic ages recorded anomalous events and these astronomical records in the historical documents provide uniquely valuable information for modern scientists. Records with drawings are particularly useful, as the verbal expressions recorded by pre-telescopic observers, who did not know the physical nature of the phenomena, are often ambiguous. However, drawings for specific datable events in the historical documents are much fewer than the verbal records. Therefore, in this paper, we show the possible earliest drawings of datable auroras and a two-tail comet in a manuscript of the Chronicle of Z\=uqn\=in, a Syriac chronicle up to 775/776 CE to interpret their nature. They provide not only the historical facts in the realm around Amida but also information about low-latitude aurora observations due to extreme space weather events and the existence of sun-grazing comets.


title: Type oriented parallel programming for Exascale

abstract: Whilst there have been great advances in HPC hardware and software in recent years, the languages and models that we use to program these machines have remained much more static. This is not from a lack of effort, but instead by virtue of the fact that the foundation that many programming languages are built on is not sufficient for the level of expressivity required for parallel work. The result is an implicit trade-off between programmability and performance which is made worse due to the fact that, whilst many scientific users are experts within their own fields, they are not HPC experts. Type oriented programming looks to address this by encoding the complexity of a language via the type system. Most of the language functionality is contained within a loosely coupled type library that can be flexibly used to control many aspects such as parallelism. Due to the high level nature of this approach there is much information available during compilation which can be used for optimisation and, in the absence of type information, the compiler can apply sensible default options thus supporting both the expert programmer and novice alike. We demonstrate that, at no performance or scalability penalty when running on up to 8196 cores of a Cray XE6 system, codes written in this type oriented manner provide improved programmability. The programmer is able to write simple, implicit parallel, HPC code at a high level and then explicitly tune by adding additional type information if required.


title: Remembering Guido Altarelli

abstract: In memory of Guido Altarelli I present my personal recollections of the early times and his major role played in the development of QCD.


title: Buffer-Aided Relaying For The Two-Hop Full-Duplex Relay Channel With Self-Interference

abstract: In this paper, we investigate the two-hop full-duplex (FD) relay channel with self-interference and fading, which is comprised of a source, an FD relay, and a destination, where a direct source-destination link does not exist and the FD relay is impaired by self-interference. For this channel, we propose three buffer-aided relaying schemes with adaptive reception-transmission at the FD relay for the cases when the source and the relay both perform adaptive-power allocation, fixed-power allocation, and fixed-rate transmission, respectively. The proposed buffer-aided relaying schemes significantly improve the achievable rate and the throughput of the considered relay channel by enabling the FD relay to adaptively select to either receive, transmit, or simultaneously receive and transmit in a given time slot based on the qualities of the receiving, transmitting, and self-interference channels. Our numerical results show that significant performance gains are achieved using the proposed buffer-aided relaying schemes compared to conventional FD relaying, where the FD relay is forced to always simultaneously receive and transmit, and to buffer-aided half-duplex relaying, where the half-duplex relay cannot simultaneously receive and transmit.


title: CogALex-V Shared Task: LexNET - Integrated Path-based and Distributional Method for the Identification of Semantic Relations

abstract: We present a submission to the CogALex 2016 shared task on the corpus-based identification of semantic relations, using LexNET (Shwartz and Dagan, 2016), an integrated path-based and distributional method for semantic relation classification. The reported results in the shared task bring this submission to the third place on subtask 1 (word relatedness), and the first place on subtask 2 (semantic relation classification), demonstrating the utility of integrating the complementary path-based and distributional information sources in recognizing concrete semantic relations. Combined with a common similarity measure, LexNET performs fairly good on the word relatedness task (subtask 1). The relatively low performance of LexNET and all other systems on subtask 2, however, confirms the difficulty of the semantic relation classification task, and stresses the need to develop additional methods for this task.


title: Implementation of traveling odd Schr$\"o$dinger cat states in circuit-QED

abstract: We propose a realistic scheme of generating a traveling odd Schr$\"o$dinger cat state and a generalized entangled coherent state in circuit quantum electrodynamics (circuit-QED). A squeezed vacuum state is used as initial resource of nonclassical states, which can be created through a Josephson traveling-wave parametric amplifier, and travels through a transmission line. Because a single-photon subtraction from the squeezed vacuum gives with very high fidelity an odd Schr$\"o$dinger cat state, we consider a specific circuit-QED setup consisting of the Josephson amplifier creating the traveling resource in a line, a beam-splitter coupling two transmission lines, and a single photon detector located at the end of the other line. When a single microwave photon is detected by measuring the excited state of a superconducting qubit in the detector, a heralded cat state is generated with high fidelity in the opposite line. For example, we show that the high fidelity of the outcome with the ideal cat state can be achieved with appropriate squeezing parameters theoretically. As its extended setup, we suggest that generalized entangled coherent states can be also built probabilistically and useful for microwave quantum information processing for error-correctable qudits in circuit-QED.


title: Learning Bound for Parameter Transfer Learning

abstract: We consider a transfer-learning problem by using the parameter transfer approach, where a suitable parameter of feature mapping is learned through one task and applied to another objective task. Then, we introduce the notion of the local stability and parameter transfer learnability of parametric feature mapping,and thereby derive a learning bound for parameter transfer algorithms. As an application of parameter transfer learning, we discuss the performance of sparse coding in self-taught learning. Although self-taught learning algorithms with plentiful unlabeled data often show excellent empirical performance, their theoretical analysis has not been studied. In this paper, we also provide the first theoretical learning bound for self-taught learning.


title: Proximal primal-dual best approximation algorithm with memory

abstract: We propose a new modified primal-dual proximal best approximation method for solving convex not necessarily differentiable optimization problems. The novelty of the method relies on introducing memory by taking into account iterates computed in previous steps in the formulas defining current iterate. To this end we consider projections onto intersections of halfspaces generated on the basis of the current as well as the previous iterates. To calculate these projections we are using recently obtained closed-form expressions for projectors onto polyhedral sets. The resulting algorithm with memory inherits strong convergence properties of the original best approximation proximal primal-dual algorithm. Additionally, we compare our algorithm with the original (non-inertial) one with the help of the so called attraction property defined below. Extensive numerical experimental results on image reconstruction problems illustrate the advantages of including memory into the original algorithm.


title: The HADES RV Programme with HARPS-N@TNG II. Data treatment and simulations

abstract: The distribution of exoplanets around low-mass stars is still not well understood. Such stars, however, present an excellent opportunity of reaching down to the rocky and habitable planet domains. The number of current detections used for statistical purposes is still quite modest and different surveys, using both photometry and precise radial velocities, are searching for planets around M dwarfs. Our HARPS-N red dwarf exoplanet survey is aimed at the detection of new planets around a sample of 78 selected stars, together with the subsequent characterization of their activity properties. Here we investigate the survey performance and strategy. From 2700 observed spectra, we compare the radial velocity determinations of the HARPS-N DRS pipeline and the HARPS-TERRA code, we calculate the mean activity jitter level, we evaluate the planet detection expectations, and we address the general question of how to define the strategy of spectroscopic surveys in order to be most efficient in the detection of planets. We find that the HARPS-TERRA radial velocities show less scatter and we calculate a mean activity jitter of 2.3 m/s for our sample. For a general radial velocity survey with limited observing time, the number of observations per star is key for the detection efficiency. In the case of an early M-type target sample, we conclude that approximately 50 observations per star with exposure times of 900 s and precisions of about 1 m/s maximizes the number of planet detections.


title: Topological rigidity fails for quotients of the Davis complex

abstract: A Coxeter group acts properly and cocompactly by isometries on the Davis complex for the group; we call the quotient of the Davis complex under this action the Davis orbicomplex for the group. We prove the set of finite covers of the Davis orbicomplexes for the set of one-ended Coxeter groups is not topologically rigid. We exhibit a quotient of a Davis complex by a one-ended right-angled Coxeter group which has two finite covers that are homotopy equivalent but not homeomorphic. We discuss consequences for the abstract commensurability classification of Coxeter groups.


title: Admissibility in Positive Logics

abstract: The paper studies admissibility of multiple-conclusion rules in the positive logics. Using modification of a method used by M.~Wajsberg in the proof of the separation theorem, it is shown that the problem of admissibility in positive logics is equivalent to the problem of admissibility in intermediate logics defined by positive additional axioms.


title: Resolving Three Single-Photon Emitters By Only Detecting Third-Order Correlation Functions

abstract: Rayleigh criterion has been continuously broken through with modern detecting instruments, many methods on resolving two emitters below the classical diffraction limit has been developed, one of these methods is detecting high order correlation functions of the emitters. We propose a method to resolve the single-photon emitters with only detecting the highest order correlation functions. We find that the variables that characterize the emitters are contained in the highest order correlation functions. Especially when there are more than two emitters, the spatial distribution of the emitters can be deduced only from the highest order correlation functions, which gives a new vision and a new method on beating classical diffraction limit using high order correlation functions.


title: On the flat geometry of the cuspidal edge

abstract: We study the geometry of the cuspidal edge $M$ in $\mathbb R^3$ derived from its contact with planes and lines (referred to as flat geometry). The contact of $M$ with planes is measured by the singularities of the height functions on $M$. We classify submersions on a model of $M$ by diffeomorphisms and recover the contact of $M$ with planes from that classification. The contact of $M$ with lines is measured by the singularities of orthogonal projections of $M$. We list the generic singularities of the projections and obtain the generic deformations of the apparent contour (profile) when the direction of projection varies locally in $S^2$. We also relate the singularities of the height functions and of the projections to some geometric invariants of the cuspidal edge.


title: Identification of Fully Physical Consistent Inertial Parameters using Optimization on Manifolds

abstract: This paper presents a new condition, the fully physical consistency for a set of inertial parameters to determine if they can be generated by a physical rigid body. The proposed condition ensure both the positive definiteness and the triangular inequality of 3D inertia matrices as opposed to existing techniques in which the triangular inequality constraint is ignored. This paper presents also a new parametrization that naturally ensures that the inertial parameters are fully physical consistency. The proposed parametrization is exploited to reformulate the inertial identification problem as a manifold optimization problem, that ensures that the identified parameters can always be generated by a physical body. The proposed optimization problem has been validated with a set of experiments on the iCub humanoid robot.


title: Negative refractive index induced by percolation in disordered metamaterials

abstract: An effective medium model is developed for disordered metamaterials containing a spatially random distribution of dielectric spheres. Similar to effective medium models for ordered metamaterials, this model predicts resonances in the effective permeability and permittivity arising from electric- and magnetic-dipole Mie resonances in the spheres. In addition, the model predicts a redshift of the electric resonance with increasing particle loading. Interestingly, when the particle loading exceeds the percolation threshold of 33\%, the model predicts that the electric resonance overlaps with the magnetic resonance, resulting in a negative refractive index.


title: Accelerating BLAS and LAPACK via Efficient Floating Point Architecture Design

abstract: Basic Linear Algebra Subprograms (BLAS) and Linear Algebra Package (LAPACK) form basic building blocks for several High Performance Computing (HPC) applications and hence dictate performance of the HPC applications. Performance in such tuned packages is attained through tuning of several algorithmic and architectural parameters such as number of parallel operations in the Directed Acyclic Graph of the BLAS/LAPACK routines, sizes of the memories in the memory hierarchy of the underlying platform, bandwidth of the memory, and structure of the compute resources in the underlying platform. In this paper, we closely investigate the impact of the Floating Point Unit (FPU) micro-architecture for performance tuning of BLAS and LAPACK. We present theoretical analysis for pipeline depth of different floating point operations like multiplier, adder, square root, and divider followed by characterization of BLAS and LAPACK to determine several parameters required in the theoretical framework for deciding optimum pipeline depth of the floating operations. A simple design of a Processing Element (PE) is presented and shown that the PE outperforms the most recent custom realizations of BLAS and LAPACK by 1.1X to 1.5X in Gflops/W, and 1.9X to 2.1X in Gflops/mm^2.


title: Remarks on Local Lie algebras of pairs of functions

abstract: We study local Lie algebras of pairs of functions which generate infinitesimal symmetries of almost-cosymplectic-contact structures of odd dimensional manifolds.


title: Exact solutions of an anisotropic universe in a modified teleparallel gravity model via the Noether and B.N.S. approaches

abstract: In this paper, we present the Noether symmetries of locally rotationally symmetric Bianchi type I (LRS BI), an anisotropic model, in the context of the teleparallel gravity. We study a certain modified teleparallel theory based on the action that, in particular, contains a coupling between the scalar field and field strength (magnetism part). We derive the symmetry generators and show that, by means of cyclic variables approach, we can not obtain a suitable solution for field equations. Hence by the use of B.N.S. approach, we solve the equations which carry Noether currents as well. By data analysis of the obtained results, we show compatible results with observational data at the last half the age of universe which is accelerating.


title: Percolation in real multiplex networks

abstract: We present an exact mathematical framework able to describe site-percolation transitions in real multiplex networks. Specifically, we consider the average percolation diagram valid over an infinite number of random configurations where nodes are present in the system with given probability. The approach relies on the locally treelike ansatz, so that it is expected to accurately reproduce the true percolation diagram of sparse multiplex networks with negligible number of short loops. The performance of our theory is tested in social, biological, and transportation multiplex graphs. When compared against previously introduced methods, we observe improvements in the prediction of the percolation diagrams in all networks analyzed. Results from our method confirm previous claims about the robustness of real multiplex networks, in the sense that the average connectedness of the system does not exhibit any significant abrupt change as its individual components are randomly destroyed.


title: Tracking in dense environments and its inefficiency measurement using pixel $dE/dx$

abstract: We present a measurement of the charged particle reconstruction inefficiency inside of jet cores, using data collected by the ATLAS experiment in 2015 of $pp$ collisions produced at the LHC, at a center-of-mass energy of 13 TeV. The determination of this inefficiency is important for jet energy scale and mass calibration, as well as multiple other performance studies and analyses. A data driven method is used, where the fraction of lost particle tracks is determined from energy deposition $dE/dx$ in the pixel detector. The fraction of lost tracks is found to be less than 5%, which is an improvement since the previous study, and agrees well within systematic uncertainties with a Monte Carlo simulation.


title: On the Stability of a General Magnetic Field Topology in Stellar Radiative Zones

abstract: This paper provides a brief overview of the formation of stellar fossil magnetic fields and what potential instabilities may occur given certain configurations of the magnetic field. In particular, a purely magnetic instability can occur for poloidal, toroidal, and mixed poloidal-toroidal axisymmetric magnetic field configurations as originally studied in Tayler (1973), Markey & Tayler (1973), and Tayler (1980). However, most of the magnetic field configurations observed at the surface of massive stars are non- axisymmetric. Thus, extending earlier studies of the axisymmetric Tayler instability in spherical geometry (Goossens, 1980), we introduce a formulation for the global change in the potential energy contained in a convectively-stable region given an arbitrary Lagrangian perturbation, which permits the inclusion of both axisymmetric and non-axisymmetric magnetic fields. With this tool in hand, a path is shown by which more general stability criterion can be established.


title: Evolution of CMB spectral distortion anisotropies and tests of primordial non-Gaussianity

abstract: Anisotropies in distortions to the frequency spectrum of the cosmic microwave background (CMB) can be created through spatially varying heating processes in the early Universe. For instance, the dissipation of small-scale acoustic modes does create distortion anisotropies, in particular for non-Gaussian primordial perturbations. In this work, we derive approximations that allow describing the associated distortion field. We provide a systematic formulation of the problem using Fourier-space window functions, clarifying and generalizing previous approximations. Our expressions highlight the fact that the amplitudes of the spectral-distortion fluctuations induced by non-Gaussianity depend also on the homogeneous value of those distortions. Absolute measurements are thus required to obtain model-independent distortion constraints on primordial non-Gaussianity. We also include a simple description for the evolution of distortions through photon diffusion, showing that these corrections can usually be neglected. Our formulation provides a systematic framework for computing higher order correlation functions of distortions with CMB temperature anisotropies and can be extended to describe correlations with polarization anisotropies.


title: Dynamical Properties of Internal Shocks Revisited

abstract: Internal shocks between propagating plasma shells, originally ejected at different times with different velocities are believed to play a major role in dissipating the kinetic energy, thereby explaining the observed lightcurve and spectra in a large range of transient objects. Even if initially the colliding plasmas are cold, following the first collision the plasma shells are substantially heated, implying that in a scenario of multiple collisions, most collisions take place between plasmas of non-zero temperatures. Here, we calculate the dynamical properties of plasmas resulting from collision between arbitrarily hot plasma shells, moving at arbitrary speeds. We provide simple analytical expressions valid for both the ultra-relativistic and Newtonian velocities, for both hot and cold plasmas. We derive the minimum criteria required for the formation of the two-shock wave system, and show that in the relativistic limit, the minimum Lorentz factor is proportional to the square root of the ratio of the initial plasmas enthalpies. We provide basic scaling laws of synchrotron emission from both the forward and reverse shock waves, and show how these can be used to deduce the properties of the colliding shells. Finally, we discuss the implications of these results in the study of several astronomical transients, such as x-ray binaries, radio load quasars and gamma-ray bursts.


title: The Probabilistic Model Checker Storm (Extended Abstract)

abstract: We present a new probabilistic model checker Storm. Using state-of-the-art libraries, we aim for both high performance and versatility. This extended abstract gives a brief overview of the features of Storm.


title: Coulomb spin liquid in anion-disordered pyrochlore Tb$_2$Hf$_2$O$_7$

abstract: The charge ordered structure of ions and vacancies characterizing rare-earth pyrochlore oxides serves as a model for the study of geometrically frustrated magnetism. The organization of magnetic ions into networks of corner-sharing tetrahedra gives rise to highly correlated magnetic phases with strong fluctuations, including spin liquids and spin ices. It is an open question how these ground states governed by local rules are affected by disorder. In the pyrochlore Tb$_2$Hf$_2$O$_7$, we demonstrate that the vicinity of the disordering transition towards a defective fluorite structure translates into a tunable density of anion Frenkel disorder while cations remain ordered. Quenched random crystal fields and disordered exchange interactions can therefore be introduced into otherwise perfect pyrochlore lattices of magnetic ions. We show that disorder can play a crucial role in preventing long-range magnetic order at low temperatures, and instead induces a strongly-fluctuating Coulomb spin liquid with defect-induced frozen magnetic degrees of freedom.


title: The procedure of excluding of the nonlinear trend fir the models described by stochastic differential and difference equations

abstract: We consider the diffusion process and its approximation by Markov chain with nonlinear increasing trends. The usual parametrix method is not appliable because these models have unbounded trends. We describe a procedure that allows to exclude nonlinear growing trend and move to stochastic differential equation with reduced drift and diffusion coefficients. A similar procedure is considered for a Markov chain


title: Uncertainty propagation with functionally correlated quantities

abstract: Many uncertainty propagation software exist, written in different programming languages, but not all of them are able to handle functional correlation between quantities. In this paper we review one strategy to deal with uncertainty propagation of quantities that are functionally correlated, and introduce a new software offering this feature: the Julia package Measurements.jl. It supports real and complex numbers with uncertainty, arbitrary-precision calculations, mathematical and linear algebra operations with matrices and arrays.


title: Reins to the Cloud: Compromising Cloud Systems via the Data Plane

abstract: Virtual switches have become popular among cloud operating systems to interconnect virtual machines in a more flexible manner. However, this paper demonstrates that virtual switches introduce new attack surfaces in cloud setups, whose effects can be disastrous. Our analysis shows that these vulnerabilities are caused by: (1) inappropriate security assumptions (privileged virtual switch execution in kernel and user space), (2) the logical centralization of such networks (e.g., OpenStack or SDN), (3) the presence of bi-directional communication channels between data plane systems and the centralized controller, and (4) non-standard protocol parsers. Our work highlights the need to accommodate the data plane(s) in our threat models. In particular, it forces us to revisit today's assumption that the data plane can only be compromised by a sophisticated attacker: we show that compromising the data plane of modern computer networks can actually be performed by a very simple attacker with limited resources only and at low cost (i.e., at the cost of renting a virtual machine in the Cloud). As a case study, we fuzzed only 2\% of the code-base of a production quality virtual switch's packet processor (namely OvS), identifying serious vulnerabilities leading to unauthenticated remote code execution. In particular, we present the "rein worm" which allows us to fully compromise test-setups in less than 100 seconds. We also evaluate the performance overhead of existing mitigations such as ASLR, PIEs, and unconditional stack canaries on OvS. We find that while applying these countermeasures in kernel-space incurs a significant overhead, in user-space the performance overhead is negligible.


title: Predicting seasonal influenza transmission using Regression Models with Temporal Dependence

abstract: In this manuscript, we use meteorological information in Galicia (Spain) to propose a novel approach to predict the incidence of influenza. Our approach extends the GLS methods in the multivariate framework to functional regression models with dependent errors. A simulation study shows that the GLS estimators render better estimations of the parameters associated with the regression model and obtain extremely good results from the predictive point of view. Thus they improve the classical linear approach. It proposes an iterative version of the GLS estimator (called iGLS) that can help to model complicated dependence structures, uses the distance correlation measure $\mathcal{R}$ to select relevant information to predict influenza rate and applies the GLS procedure to the prediction of the influenza rate using readily available functional variables. These kinds of models are extremely useful to health managers in allocating resources in advance for an epidemic outbreak.


title: Spinning particles coupled to gravity

abstract: Recent experimental work has determined that free falling $^{87}$Rb atoms on Earth, with vertically aligned spins, follow geodesics, thus apparently ruling out spin--gravitation interactions. It is showed that while some spinning matter models coupled to gravitation referenced to in that work seem to be ruled out by the experiment, those same experimental results confirm theoretical results derived from a Lagrangian description of spinning particles coupled to gravity constructed over forty years ago. A proposal to carry out (similar but) different experiments which will help to test the validity of the Universality of Free Fall as opposed to the correctness of the aforementioned Lagrangian theory, is presented.


title: Sensitivity of the mixing current technique to detect nano-mechanical motion

abstract: Detection of nano-mechanical displacement by transport techniques has reached high level of sensitivity and versatility. In order to detect the amplitude of oscillation of nano-mechanical oscillator a widely used technique consists to couple this motion capacitively to a single-electron transistor and to detect the high-frequency modulation of the current through the non-linear mixing with an electric signal at a slighltly detuned frequency. The method known as current-mixing technique is employed in particular for the detection of suspended carbon nanotubes. In this paper we study theoretically the limiting conditions on the sensitivity of this method. The sensitivity is increased by increasing the response function to the signal, but also by reducing the noise. For these reasons we study systematically the response function, the effect of current-and displacement-fluctuations, and finally the case where the tunnelling rate of the electrons are of the same order or larger of the resonating frequency. We find thus upper bounds to the sensitivity of the detection technique.


title: Cross section measurement of alpha particle induced nuclear reactions on natural cadmium up to 52 MeV

abstract: Cross sections of alpha particle induced nuclear reactions have been measured on thin natural cadmium targets foils in the energy range from 11 to 51.2 MeV. This work was a part of our systematic study on excitation functions of light ion induced nuclear reactions on different target materials. Regarding the cross sections, the alpha induced reactions are not deeply enough investigated. Some of the produced isotopes are of medical interest, others have application in research and industry. The radioisotope $^{117m}$Sn is a very important theranostic (therapeutic + diagnostic) radioisotope, so special care was taken to the results for that isotope. The well-established stacked foil technique followed by gamma-spectrometry with HPGe gamma spectrometers were used. The target and monitor foils in the stack were commercial high purity metal foils. From the irradiated targets $^{117m}$Sn, $^{113}$Sn, $^{110}$Sn, $^{117m,g}$In, $^{116m}$In, $^{115m}$In, $^{114m}$In, $^{113m}$In, $^{111}$In, $^{110m,g}$In, $^{109m}$In, $^{108m}$In, $^{115g}$Cd and $^{111m}$Cd were identified and their excitation functions were derived. The results were compared with the data of the previous measurements from the literature and with the results of the theoretical nuclear reaction model code calculations TALYS 1.8 (TENDL-2015) and EMPIRE 3.2 (Malta). From the cross section curves thick target yields were calculated and compared with the available literature data.


title: Improving PageRank for Local Community Detection

abstract: Community detection is a classical problem in the field of graph mining. While most algorithms work on the entire graph, it is often interesting in practice to recover only the community containing some given set of seed nodes. In this paper, we propose a novel approach to this problem, using some low-dimensional embedding of the graph based on random walks starting from the seed nodes. From this embedding, we propose some simple yet efficient versions of the PageRank algorithm as well as a novel algorithm, called WalkSCAN, that is able to detect multiple communities, possibly overlapping. We provide insights into the performance of these algorithms through the theoretical analysis of a toy network and show that WalkSCAN outperforms existing algorithms on real networks.


title: Path dependent equations driven by H\"older processes

abstract: This paper investigates existence results for path-dependent differential equations driven by a H{\"o}lder function where the integrals are understood in the Young sense. The two main results are proved via an application of Schauder theorem and the vector field is allowed to be unbounded. The H{\"o}lder function is typically the trajectory of a stochastic process.


title: Bifurcation induced by the aspect ratio in a turbulent Von-K\'arm\'an swirling flow

abstract: We evaluate the effect of two experimental parameters on the slow dynamics of a Von-K\'arm\'an swirling flow driven by two propellers in a closed cylinder. The first parameter is the inertia mo- mentum of the propellers, and the second parameter is the aspect ratio, i.e. the distance between the propellers $H$ divided by the diameter $D$. We use a cell with a fixed diameter $D$ but where the distance between the propellers can be turned continuously and where the inertia from the pro- pellers can also be changed using different gears. No change on the dynamics is observed when the momentum of inertia is modified. Some dramatic changes of the shear layer position are observed modifying the aspect ratio $\Gamma=H/D$. A bifurcation of the shear layer position appears. Whereas for low $\Gamma$ the shear layer position has a smooth evolution when turning the asymmetry between the rotation frequency of the propellers, for high $\Gamma$ the transition becomes abrupt and a symmetry breaking appears. Secondly we observe that the spontaneous reversals already observed in this ex- periment for $\Gamma=1$ [de la Torre & Burguete PRL 99, 054101 (2007)] exist only in a strait window of aspect ratio. We show using an experimental study of the mean flow structure and a numerical approach based on a Langevin equation with coloured noise that the shear layer position seems to be decided by the mean flow structure whereas the reversals are linked to the spatial distribution of the turbulent fluctuations in the cell.


title: Constraints on the Primordial Black Hole Abundance from the First Advanced LIGO Observation Run Using the Stochastic Gravitational-Wave Background

abstract: Advanced LIGO's discovery of gravitational-wave events is stimulating extensive studies on the origin of binary black holes. Assuming that the gravitational-wave events can be explained by binary primordial black hole mergers, we utilize the upper limits on the stochastic gravitational-wave background given by Advanced LIGO as a new observational window to independently constrain the abundance of primordial black holes in dark matter. We show that Advanced LIGO's first observation run gives the best constraint on the primordial black hole abundance in the mass range $1 M_\odot \lesssim M_\text{PBH}\lesssim 100 M_\odot$, pushing the previous microlensing and dwarf galaxy dynamics constraints tighter by 1 order of magnitude. Moreover, we discuss the possibility to detect the stochastic gravitational-wave background from primordial black holes, in particular from subsolar mass primordial black holes, by Advanced LIGO in the near future.


title: Wilf's conjecture for numerical semigroups

abstract: Let $S\subseteq \mathbb{N}$ be a numerical semigroup with multiplicity $m$, embedding dimension $\nu$ and conductor $c=f+1=qm-\rho$ for some $q,\rho\in\mathbb{N}$ with $\rho<m$. Let Ap$(S,m) = \{w\_0<w_1 < \ldots < w_{m-1}\}$ be the Ap\'ery set of $S$. The aim of this paper is to prove Wilf's Conjecture in some special cases. First, we prove that if $w_{m-1}\geq w_1+w_\alpha$ and $(2+\frac{\alpha-3}{q})\nu\geq m$ for some $1<\alpha<m-1$, then $S$ satisfies Wilf's Conjecture. Then, we prove the conjecture in the following cases: $(2+\frac{1}{q})\nu\geq m$, $m-\nu\leq 5$ and $m=9$. Finally, the conjecture is proved if $w_{m-1} \geq w_{\alpha-1} + w_\alpha$ and $(\frac{\alpha+3}{3})\nu\geq m$ for some $1<\alpha<m-1$.


title: Combining Usability and Privacy Protection in Free-Access Public Cloud Storage Servers: Review of the Main Threats and Challenges

abstract: The 21st century belongs to the world of computing, specially as a result of the so-called cloud computing. This technology enables ubiquitous information management and thus people can access all their data from any place and at any time. In this landscape, the emergence of cloud storage has had an important role in the last five years. Nowadays, several free-access public cloud storage services make it possible for users to have a free backup of their assets and to manage and share them, representing a low-cost opportunity for Small and Medium Companies (SME). However, the adoption of cloud storage involves data outsourcing, so a user does not have the guarantee about the way her data will be processed and protected. Therefore, it seems necessary to endow public cloud storage with a set of means to protect users' confidentiality and privacy, to assess data integrity and to guarantee a proper backup of information assets. Along this paper we discuss the main challenges to achieve such a goal, underlining the set of functionalities already implemented in the most popular public cloud storage services.


title: Two-dimensional topological superconducting phases emerged from d-wave superconductors in proximity to antiferromagnets

abstract: Motivated by the recent observations of nodeless superconductivity in the monolayer CuO$_{2}$ grown on the Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8+\delta }$ substrates, we study the two-dimensional superconducting (SC) phases described by the two-dimensional $t$-$J$ model in proximity to an antiferromagnetic (AF) insulator. We found that (i) the nodal d-wave SC state can be driven via a continuous transition into a nodeless d-wave pairing state by the proximity induced AF field. (ii) The energetically favorable pairing states in the strong field regime have extended s-wave symmetry and can be nodal or nodeless. (iii) Between the pure d-wave and s-wave paired phases, there emerge two topologically distinct SC phases with ($s+$i$d$) symmetry, i.e., the weak and strong pairing phases, and the weak pairing phase is found to be a $Z_{2}$ topological superconductor protected by valley symmetry, exhibiting robust gapless non-chiral edge modes. These findings strongly suggest that the high-$T_{c}$ superconductors in proximity to antiferromagnets can realize fully gapped symmetry protected topological SC.


title: Shadow lines in the arithmetic of elliptic curves

abstract: Let E/Q be an elliptic curve and p a rational prime of good ordinary reduction. For every imaginary quadratic field K/Q satisfying the Heegner hypothesis for E we have a corresponding line in E(K)\otimes Q_p, known as a shadow line. When E/Q has analytic rank 2 and E/K has analytic rank 3, shadow lines are expected to lie in E(Q)\otimes Q_p. If, in addition, p splits in K/Q, then shadow lines can be determined using the anticyclotomic p-adic height pairing. We develop an algorithm to compute anticyclotomic p-adic heights which we then use to provide an algorithm to compute shadow lines. We conclude by illustrating these algorithms in a collection of examples.


title: Smooth models for the Coulomb potential

abstract: Smooth model potentials with parameters selected to reproduce the spectrum of one-electron atoms are used to approximate the singular Coulomb potential. Even when the potentials do not mimic the Coulomb singularity, much of the spectrum is reproduced within the chemical accuracy. For the Hydrogen atom, the smooth approximations to the Coulomb potential are more accurate for higher angular momentum states. The transferability of the model potentials from an attractive interaction (Hydrogen atom) to a repulsive one (Harmonium and the uniform electron gas) is discussed.


title: Weighted Bergman Projections on the Hartogs Triangle: Exponential Decay

abstract: We study weighted Bergman projections on the Hartogs triangle in $\mathbb{C}^2$. We show that projections corresponding to exponentially vanishing weights have degenerate $L^p$ mapping properties.


title: On exponential functionals of processes with independent increments

abstract: In this paper we study the exponential functionals of the processes $X$ with independent increments , namely $$I_t= \int _0^t\exp(-X_s)ds, _,\,\, t\geq 0,$$ and also $$I_{\infty}= \int _0^{\infty}\exp(-X_s)ds.$$ When $X$ is a semi-martingale with absolutely continuous characteristics, we derive recurrent integral equations for Mellin transform ${\bf E}( I_t^{\alpha})$, $\alpha\in\mathbb{R}$, of the integral functional $I_t$. Then we apply these recurrent formulas to calculate the moments. We present also the corresponding results for the exponential functionals of Levy processes, which hold under less restrictive conditions then in the paper of Bertoin, Yor (2005). In particular, we obtain an explicit formula for the moments of $I_t$ and $I_{\infty}$, and we precise the exact number of finite moments of $I_{\infty}$.


title: GPflow: A Gaussian process library using TensorFlow

abstract: GPflow is a Gaussian process library that uses TensorFlow for its core computations and Python for its front end. The distinguishing features of GPflow are that it uses variational inference as the primary approximation method, provides concise code through the use of automatic differentiation, has been engineered with a particular emphasis on software testing and is able to exploit GPU hardware.


title: High quality electron beam generation in a proton-driven hollow plasma wakefield accelerator

abstract: Simulations of proton-driven plasma wakefield accelerators have demonstrated substantially higher accelerating gradients compared to conventional accelerators and the viability of accelerating electrons to the energy frontier in a single plasma stage. However, due to the strong intrinsic transverse fields varying both radially and in time, the witness beam quality is still far from suitable for practical application in future colliders. Here we demonstrate efficient acceleration of electrons in proton-driven wakefields in a hollow plasma channel. In this regime, the witness bunch is positioned in the region with a strong accelerating field, free from plasma electrons and ions. We show that the electron beam carrying the charge of about 10% of 1 TeV proton driver charge can be accelerated to 0.6 TeV with preserved normalized emittance in a single channel of 700 m. This high quality and high charge beam may pave the way for the development of future plasma-based energy frontier colliders.


title: Stratification of patient trajectories using covariate latent variable models

abstract: Standard models assign disease progression to discrete categories or stages based on well-characterized clinical markers. However, such a system is potentially at odds with our understanding of the underlying biology, which in highly complex systems may support a (near-)continuous evolution of disease from inception to terminal state. To learn such a continuous disease score one could infer a latent variable from dynamic "omics" data such as RNA-seq that correlates with an outcome of interest such as survival time. However, such analyses may be confounded by additional data such as clinical covariates measured in electronic health records (EHRs). As a solution to this we introduce covariate latent variable models, a novel type of latent variable model that learns a low-dimensional data representation in the presence of two (asymmetric) views of the same data source. We apply our model to TCGA colorectal cancer RNA-seq data and demonstrate how incorporating microsatellite-instability (MSI) status as an external covariate allows us to identify genes that stratify patients on an immune-response trajectory. Finally, we propose an extension termed Covariate Gaussian Process Latent Variable Models for learning nonparametric, nonlinear representations. An R package implementing variational inference for covariate latent variable models is available at http://github.com/kieranrcampbell/clvm.


title: Ekeland's inverse function theorem in graded Fr\'echet spaces revisited for multifunctions

abstract: In this paper, we present some implicit function theorems for set-valued mappings between Fr\'echet spaces. The proof relies on Lebesgue's Dominated Convergence Theorem and on Ekeland's variational principle. An application to the existence of solutions of differential equations in Fr\'echet spaces with non-smooth data is given.


title: A nuclear quantum effect with pure anharmonicity and the anomalous thermal expansion of silicon

abstract: Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well-understood. Adapting harmonic phonons to the specific volume at temperature, the quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test atomistic mechanisms, we performed inelastic neutron scattering experiments from 100-1500K on a single-crystal of silicon to measure the changes in phonon frequencies. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity and nuclear quantum effects, reproduced the measured shifts of individual phonons with temperature, whereas quasiharmonic shifts were mostly of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a fortuitous cancellation of contributions from individual phonons.


title: Compressive K-means

abstract: The Lloyd-Max algorithm is a classical approach to perform K-means clustering. Unfortunately, its cost becomes prohibitive as the training dataset grows large. We propose a compressive version of K-means (CKM), that estimates cluster centers from a sketch, i.e. from a drastically compressed representation of the training dataset. We demonstrate empirically that CKM performs similarly to Lloyd-Max, for a sketch size proportional to the number of cen-troids times the ambient dimension, and independent of the size of the original dataset. Given the sketch, the computational complexity of CKM is also independent of the size of the dataset. Unlike Lloyd-Max which requires several replicates, we further demonstrate that CKM is almost insensitive to initialization. For a large dataset of 10^7 data points, we show that CKM can run two orders of magnitude faster than five replicates of Lloyd-Max, with similar clustering performance on artificial data. Finally, CKM achieves lower classification errors on handwritten digits classification.


title: Finding Largest Common Substructures of Molecules in Quadratic Time

abstract: Finding the common structural features of two molecules is a fundamental task in cheminformatics. Most drugs are small molecules, which can naturally be interpreted as graphs. Hence, the task is formalized as maximum common subgraph problem. Albeit the vast majority of molecules yields outerplanar graphs this problem remains NP-hard. We consider a variation of the problem of high practical relevance, where the rings of molecules must not be broken, i.e., the block and bridge structure of the input graphs must be retained by the common subgraph. We present an algorithm for finding a maximum common connected induced subgraph of two given outerplanar graphs subject to this constraint. Our approach runs in time $\mathcal{O}(\Delta n^2)$ in outerplanar graphs on $n$ vertices with maximum degree $\Delta$. This leads to a quadratic time complexity in molecular graphs, which have bounded degree. The experimental comparison on synthetic and real-world datasets shows that our approach is highly efficient in practice and outperforms comparable state-of-the-art algorithms.


title: Polarization-free Quantization of Linear Field Theories

abstract: It is well-known that there exist infinitely-many inequivalent representations of the canonical (anti)-commutation relations of Quantum Field Theory (QFT). A way out, suggested by Algebraic QFT, is to instead define the quantum theory as encompassing all possible (abstract) states. In the present paper, we describe a quantization scheme for general linear (aka. free) field theories that can be seen as intermediate between traditional Fock quantization and full Algebraic QFT, in the sense that: * it provides a constructive, explicit description of the resulting space of quantum states; * it does not require the choice of a polarization, aka. the splitting of classical solutions into positive vs. negative-frequency modes: in fact, any Fock representation corresponding to a "reasonable" choice of polarization is naturally embedded; * it supports the implementation of a "large enough" class of linear symplectomorphisms of the classical, infinite-dimensional phase space. The proposed quantization (like Algebraic QFT) is notably meant for use on curved spacetimes, where the lack of a preferred choice of polarization makes the introduction of a Fock representation problematic. Accordingly, we illustrate it in the simple case of a free Klein-Gordon field on a spatially-compact, cosmological spacetime.


title: Field-Induced Quantum Phases in Frustrated Spin-Dimer Model: A Sign-Problem-Free Quantum Monte Carlo Study

abstract: The magnetization process of a frustrated spin-1/2 spin-dimer model on a square lattice is investigated by means of a sign-problem-free quantum Monte Carlo algorithm developed recently. Rich field-induced quantum phases are discovered. We find two spin superfluids satisfying different symmetries upon layer permutation. There exist as well two solid phases with distinct checkerboard patterns. Besides, two kinds of spin supersolids are observed over a finite regime of magnetic fields. The latter two phases can be stabilized when small but nonzero spin anisotropy is present. Various field-induced transitions among these phases are explored. Our findings may provide guidance to the search of these interesting quantum phases in real frustrated spin-dimer compounds.