title: Methods for removal of unwanted signals from gravity time-series: comparison using linear techniques complemented with analysis of system dynamics

abstract: The presence of undesirable dominating signals in geophysical experimental data is a challenge in many subfields. One remarkable example is surface gravimetry, where frequencies from Earth tides correspond to time-series fluctuations up to a thousand times larger than the phenomena of major interest, such as hydrological gravity effects or co-seismic gravity changes. This work discusses general methods for removal of unwanted dominating signals by applying them to 8 long-period gravity time-series of the International Geodynamics and Earth Tides Service, equivalent to the acquisition from 8 instruments in 5 locations representative of the network. We compare three different conceptual approaches for tide removal: frequency filtering, physical modelling and data-based modelling. Each approach reveals a different limitation to be considered depending on the intended application. Vestiges of tides remain in the residues for the modelling procedures, whereas the signal was distorted in different ways by the filtering and data-based procedures. The linear techniques employed were power spectral density, spectrogram, cross-correlation and classical harmonics decomposition, while the system dynamics was analysed by state-space reconstruction and estimation of the largest Lyapunov exponent. Although the tides could not be completely eliminated, they were sufficiently reduced to allow observation of geophysical events of interest above the $10 \text{ nm s}^{-2}$ level, exemplified by a hydrology-related event of $60 \text{ nm s}^{-2}$. The implementations adopted for each conceptual approach are general, so that their principles could be applied to other kinds of data affected by undesired signals composed mainly by periodic or quasi-periodic components.


title: The full automorphism group of $\overline{T}$

abstract: Let $\overline G$ be the wonderful compactification of a simple affine algebraic group $G$ of adjoint type defined over $\mathbb C.$ Let ${\overline T}\subset \overline G$ be the closure of a maximal torus $T\subset G.$ We prove that the group of all automorphisms of the variety $\overline T$ is the semi-direct product $N_G(T)\rtimes D,$ where $N_G(T)$ is the normalizer of $T$ in $G$ and $D$ is the group of all automorphisms of the Dynkin diagram, if $G\not= {\rm PSL}(2,\mathbb{C})$. Note that if $G = {\rm PSL}(2,\mathbb{C})$, then $\overline{T} = {\mathbb C}{\mathbb P}^1$ and so in this case $\text{Aut}(\overline T)= {\rm PSL}(2,\mathbb{C})$.


title: Low-energy theory for strained graphene: an approach up to second-order in the strain tensor

abstract: An analytical study of low-energy electronic excited states in an uniformly strained graphene is carried out up to second-order in the strain tensor. We report an new effective Dirac Hamiltonian with an anisotropic Fermi velocity tensor, which reveals the graphene trigonal symmetry being absent in low-energy theories to first-order in the strain tensor. In particular, we demonstrate the dependence of the Dirac-cone elliptical deformation on the stretching direction respect to graphene lattice orientation. We further analytically calculate the optical conductivity tensor of strained graphene and its transmittance for a linearly polarized light with normal incidence. Finally, the obtained analytical expression of the Dirac point shift allows a better determination and understanding of pseudomagnetic fields induced by nonuniform strains.


title: The second boundary value problem of the prescribed affine mean curvature equation and related linearized Monge-Amp\`ere equation

abstract: These lecture notes are concerned with the solvability of the second boundary value problem of the prescribed affine mean curvature equation and related regularity theory of the Monge-Amp\`ere and linearized Monge-Amp\`ere equations. The prescribed affine mean curvature equation is a fully nonlinear, fourth order, geometric partial differential equation of the following form $$\sum_{i, j=1}^n U^{ij}\frac{\partial^2}{\partial {x_i}\partial{x_j}}\left[(\det D^2 u)^{-\frac{n+1}{n+2}}\right]=f$$ where $(U^{ij})$ is the cofactor matrix of the Hessian matrix $D^2 u$ of a locally uniformly convex function $u$. Its variant is related to the problem of finding K\"ahler metrics of constant scalar curvature in complex geometry. We first introduce the background of the prescribed affine mean curvature equation which can be viewed as a coupled system of Monge-Amp\`ere and linearized Monge-Amp\`ere equations. Then we state key open problems and present the solution of the second boundary value problem that prescribes the boundary values of the solution $u$ and its Hessian determinant $\det D^2 u$. Its proof uses important tools from the boundary regularity theory of the Monge-Amp\`ere and linearized Monge-Amp\`ere equations that we will present in the lecture notes.


title: Differentiable Learning of Logical Rules for Knowledge Base Reasoning

abstract: We study the problem of learning probabilistic first-order logical rules for knowledge base reasoning. This learning problem is difficult because it requires learning the parameters in a continuous space as well as the structure in a discrete space. We propose a framework, Neural Logic Programming, that combines the parameter and structure learning of first-order logical rules in an end-to-end differentiable model. This approach is inspired by a recently-developed differentiable logic called TensorLog, where inference tasks can be compiled into sequences of differentiable operations. We design a neural controller system that learns to compose these operations. Empirically, our method outperforms prior work on multiple knowledge base benchmark datasets, including Freebase and WikiMovies.


title: The Local Limit of Random Sorting Networks

abstract: A sorting network is a geodesic path from $12 \cdots n$ to $n \cdots 21$ in the Cayley graph of $S_n$ generated by adjacent transpositions. For a uniformly random sorting network, we establish the existence of a local limit of the process of space-time locations of transpositions in a neighbourhood of $an$ for $a\in[0,1]$ as $n\to\infty$. Here time is scaled by a factor of $1/n$ and space is not scaled. The limit is a swap process $U$ on $\mathbb{Z}$. We show that $U$ is stationary and mixing with respect to the spatial shift and has time-stationary increments. Moreover, the only dependence on $a$ is through time scaling by a factor of $\sqrt{a(1-a)}$. To establish the existence of $U$, we find a local limit for staircase-shaped Young tableaux. These Young tableaux are related to sorting networks through a bijection of Edelman and Greene.


title: Galerkin Least-Squares Stabilization in Ice Sheet Modeling - Accuracy, Robustness, and Comparison to other Techniques

abstract: We investigate the accuracy and robustness of one of the most common methods used in glaciology for the discretization of the $\mathfrak{p}$-Stokes equations: equal order finite elements with Galerkin Least-Squares (GLS) stabilization. Furthermore we compare the results to other stabilized methods. We find that the vertical velocity component is more sensitive to the choice of GLS stabilization parameter than horizontal velocity. Additionally, the accuracy of the vertical velocity component is especially important since errors in this component can cause ice surface instabilities and propagate into future ice volume predictions. If the element cell size is set to the minimum edge length and the stabilization parameter is allowed to vary non-linearly with viscosity, the GLS stabilization parameter found in literature is a good choice on simple domains. However, near ice margins the standard parameter choice may result in significant oscillations in the vertical component of the surface velocity. For these cases, other stabilization techniques, such as the interior penalty method, result in better accuracy and are less sensitive to the choice of the stabilization parameter. During this work we also discovered that the manufactured solutions often used to evaluate errors in glaciology are not reliable due to high artificial surface forces at singularities. We perform our numerical experiments in both FEniCS and Elmer/Ice.


title: Divisible sandpile on Sierpinski gasket graphs

abstract: The divisible sandpile model is a growth model on graphs that was introduced by Levine and Peres as a tool to study internal diffusion limited aggregation. In this work we investigate the shape of the divisible sandpile model on the graphical Sierpinski gasket SG. We show that the shape is a ball in the graph metric of SG. Moreover we give an exact representation of the odometer function of the divisible sandpile.


title: The Robot Crawler Model on Complete k-Partite and Erd\H{o}s-R\'enyi Random Graphs

abstract: Web crawlers are used by internet search engines to gather information about the web graph. In this paper we investigate a simple process which models such software by walking around the vertices of a graph. Once initial random vertex weights have been assigned, the robot crawler traverses the graph deterministically following a greedy algorithm, always visiting the neighbour of least weight and then updating this weight to be the highest overall. We consider the maximum, minimum and average number of steps taken by the crawler to visit every vertex of firstly, complete k-partite graphs and secondly, sparse Erd\H{o}s-R\'enyi random graphs. Our work follows on from a paper of Bonato et. al. who introduced the model.


title: Convex Relaxations of Chance Constrained AC Optimal Power Flow

abstract: High penetration of renewable energy sources and the increasing share of stochastic loads require the explicit representation of uncertainty in tools such as the optimal power flow (OPF). Current approaches follow either a linearized approach or an iterative approximation of non-linearities. This paper proposes a semidefinite relaxation of a chance constrained AC-OPF which is able to provide guarantees for global optimality. Using a piecewise affine policy, we can ensure tractability, accurately model large power deviations, and determine suitable corrective control policies for active power, reactive power, and voltage. We state a tractable formulation for two types of uncertainty sets. Using a scenario-based approach and making no prior assumptions about the probability distribution of the forecast errors, we obtain a robust formulation for a rectangular uncertainty set. Alternatively, assuming a Gaussian distribution of the forecast errors, we propose an analytical reformulation of the chance constraints suitable for semidefinite programming. We demonstrate the performance of our approach on the IEEE 24 and 118 bus system using realistic day-ahead forecast data and obtain tight near-global optimality guarantees.


title: Asymptotic enumeration of graphs by degree sequence, and the degree sequence of a random graph

abstract: In this paper we relate a fundamental parameter of a random graph, its degree sequence, to a simple model of nearly independent binomial random variables. This confirms a conjecture made in 1997. As a result, many interesting functions of the joint distribution of graph degrees, such as the distribution of the median degree, become amenable to estimation. Our result is established by proving an asymptotic formula conjectured in 1990 for the number of graphs with given degree sequence. In particular, this gives an asymptotic formula for the number of $d$-regular graphs for all $d$, as $n\to\infty$.


title: Dense blowup for parabolic SPDEs

abstract: The main result of this paper is that there are examples of stochastic partial differential equations [hereforth, SPDEs] of the type $$ \partial_t u=\frac12\Delta u +\sigma(u)\eta \qquad\text{on $(0\,,\infty)\times\mathbb{R}^3$}$$ such that the solution exists and is unique as a random field in the sense of Dalang and Walsh, yet the solution has unbounded oscillations in every open neighborhood of every space-time point. We are not aware of the existence of such a construction in spatial dimensions below $3$. En route, it will be proved that there exist a large family of parabolic SPDEs whose moment Lyapunov exponents grow at least sub exponentially in its order parameter in the sense that there exist $A_1,\beta\in(0\,,1)$ such that \[ \underline{\gamma}(k) := \liminf_{t\to\infty}t^{-1}\inf_{x\in\mathbb{R}^3} \log\mathbb{E}\left(|u(t\,,x)|^k\right) \ge A_1\exp(A_1 k^\beta) \qquad\text{for all $k\ge 2$}. \] This sort of "super intermittency" is combined with a local linearization of the solution, and with techniques from Gaussian analysis in order to establish the unbounded oscillations of the sample functions of the solution to our SPDE.


title: Amplitude reconstruction from complete electroproduction experiments and truncated partial-wave expansions

abstract: We compare the methods of amplitude reconstruction, for a complete experiment and a truncated partial-wave analysis, applied to the electroproduction of pseudoscalar mesons. We give examples which show, in detail, how the amplitude reconstruction (observables measured at a single energy and angle) is related to a truncated partial-wave analysis (observables measured at a single energy and a number of angles). A connection is made to existing data.


title: Admittance Control Parameter Adaptation for Physical Human-Robot Interaction

abstract: In physical human-robot interaction, the coexistence of robots and humans in the same workspace requires the guarantee of a stable interaction, trying to minimize the effort for the operator. To this aim, the admittance control is widely used and the appropriate selection of the its parameters is crucial, since they affect both the stability and the ability of the robot to interact with the user. In this paper, we present a strategy for detecting deviations from the nominal behavior of an admittance-controlled robot and for adapting the parameters of the controller while guaranteeing the passivity. The proposed methodology is validated on a KUKA LWR 4+.


title: Optimized Secure Position Sharing with Non-trusted Servers

abstract: Today, location-based applications and services such as friend finders and geo-social networks are very popular. However, storing private position information on third-party location servers leads to privacy problems. In our previous work, we proposed a position sharing approach for secure management of positions on non-trusted servers, which distributes position shares of limited precision among servers of several providers. In this paper, we propose two novel contributions to improve the original approach. First, we optimize the placement of shares among servers by taking their trustworthiness into account. Second, we optimize the location update protocols to minimize the number of messages between mobile device and location servers.


title: Unusual enhancement of effective magnetic anisotropy with decreasing particle size in maghemite nanoparticles

abstract: Experimental results and a model are presented to explain the observed unusual enhancement of the effective magnetic anisotropy Keff with decreasing particle size D from 15 nm to 2.5 nm in {\gamma}-Fe2O3 nanoparticles (NPs). The samples include oleic acid-coated NPs with D = 2.5, 3.4, 6.3 and 7.0 nm investigated here, with the results on other sizes taken from literature. Keff is determined from the analysis of the frequency dependence of the blocking temperature TB after considering the effects of interparticle interactions on TB. The data of Keff vs. D is fit to the derived core-shell-based relation: Keff = Kb + (6KS/D) +Ksh{[1-(2d/D)]^(-3) -1}, with Kb = 1.9 x10^5 ergs/cm^3 as the bulk-like contribution of the core, KS = 0.035 ergs/cm^2 as the contribution of the surface layer, and Ksh = 1.057 x10^4 ergs/cm^3 as the contributions of spins in the shell of thickness d = 1.1 nm. This equation represents an extension of the often-used Eq.: Keff = Kb + (6KS/D). Significance of this new result is discussed.


title: Revealing Hidden Potentials of the q-Space Signal in Breast Cancer

abstract: Mammography screening for early detection of breast lesions currently suffers from high amounts of false positive findings, which result in unnecessary invasive biopsies. Diffusion-weighted MR images (DWI) can help to reduce many of these false-positive findings prior to biopsy. Current approaches estimate tissue properties by means of quantitative parameters taken from generative, biophysical models fit to the q-space encoded signal under certain assumptions regarding noise and spatial homogeneity. This process is prone to fitting instability and partial information loss due to model simplicity. We reveal unexplored potentials of the signal by integrating all data processing components into a convolutional neural network (CNN) architecture that is designed to propagate clinical target information down to the raw input images. This approach enables simultaneous and target-specific optimization of image normalization, signal exploitation, global representation learning and classification. Using a multicentric data set of 222 patients, we demonstrate that our approach significantly improves clinical decision making with respect to the current state of the art.


title: Exploring Increasing-Chord Paths and Trees

abstract: A straight-line drawing $\Gamma$ of a graph $G=(V,E)$ is a drawing of $G$ in the Euclidean plane, where every vertex in $G$ is mapped to a distinct point, and every edge in $G$ is mapped to a straight line segment between their endpoints. A path $P$ in $\Gamma$ is called increasing-chord if for every four points (not necessarily vertices) $a,b,c,d$ on $P$ in this order, the Euclidean distance between $b,c$ is at most the Euclidean distance between $a,d$. A spanning tree $T$ rooted at some vertex $r$ in $\Gamma$ is called increasing-chord if $T$ contains an increasing-chord path from $r$ to every vertex in $T$. In this paper we prove that given a vertex $r$ in a straight-line drawing $\Gamma$, it is NP-complete to determine whether $\Gamma$ contains an increasing-chord spanning tree rooted at $r$. We conjecture that finding an increasing-chord path between a pair of vertices in $\Gamma$, which is an intriguing open problem posed by Alamdari et al., is also NP-complete, and show a (non-polynomial) reduction from the 3-SAT problem.


title: Optical appearance of a compact binary system in the neighbourhood of supermassive black hole

abstract: Optical appearance of a compact binary star in the field of a supermassive black hole is modeled in a strong field regime. Expressions for the redshift, magnification coefficient and pulsar extinction time are derived. By using the vierbein formalism we have derived equations of motion of compact binary star in the external gravitational field. We have analysed both the evolution of redshift of the optical ray from the usual star or white dwarf, and the times of arrival of pulses of pulsar. The results are illustrated by a calculation for a model binary system for the case of external gravitational field of a Schwarzschild black hole. The obtained results can be used for fitting timing data from the X-ray pulsars that moves in the neighbourhood of the Galactic Center (Sgr A*).


title: Scheduling Post-Disaster Repairs in Electricity Distribution Networks

abstract: Natural disasters, such as hurricanes, earthquakes and large wind or ice storms, typically require the repair of a large number of components in electricity distribution networks. Since power cannot be restored before these repairs have been completed, optimally scheduling the available crews to minimize the cumulative duration of the customer interruptions reduces the harm done to the affected community. Considering the radial network structure of the distribution system, this repair and restoration process can be modeled as a scheduling problem with soft precedence constraints. As a benchmark, we first formulate this problem as a time-indexed ILP with valid inequalities. Three practical methods are then proposed to solve the problem: (i) an LP-based list scheduling algorithm, (ii) a single to multi-crew repair schedule conversion algorithm, and (iii) a dispatch rule based on $\rho$-factors which can be interpreted as Component Importance Measures. We show that the first two algorithms are $2$ and $\left(2 - \frac{1}{m}\right)$ approximations respectively. We also prove that the latter two algorithms are equivalent. Numerical results validate the effectiveness of the proposed methods.


title: Large-amplitude rapid X-ray variability in the narrow-line Seyfert 1 galaxy PG 1404$+$226

abstract: We present the first results from a detailed analysis of a new, long ($\sim100$ ks) XMM-Newton observation of the narrow-line Seyfert 1 galaxy PG 1404$+$226 which showed a large-amplitude, rapid X-ray variability by a factor of $\sim7$ in $\sim10$ ks with an exponential rise and a sharp fall in the count rate. We investigate the origin of the soft X-ray excess emission and rapid X-ray variability in the source through time-resolved spectroscopy and fractional root-mean-squared (rms) spectral modeling. The strong soft X-ray excess below 1 keV observed both in the time-averaged and time-resolved spectra is described by the intrinsic disk Comptonization model as well as the relativistic reflection model where the emission is intensive merely in the inner regions ($r_{\rm in}<1.7 r_{\rm g}$) of an ionized accretion disk. We detected no significant UV variability while the soft X-ray excess flux varies together with the primary power-law emission (as $F_{{\rm primary}}\propto F_{{\rm excess}}^{1.54}$), although with a smaller amplitude, as expected in the reflection scenario. The observed X-ray fractional rms spectrum is approximately constant with a drop at $\sim0.6$ keV and is described by a non-variable emission line component with the observed energy of $\sim0.6$ keV and two variable spectral components: a more variable primary power-law emission and a less variable soft excess emission. Our results suggest the `lamppost geometry' for the primary X-ray emitting hot corona which illuminates the innermost accretion disk due to strong gravity and gives rise to the soft X-ray excess emission.


title: On Hidden Sums Compatible with A Given Block Cipher Diffusion Layer

abstract: Sometimes it is possible to embed an algebraic trapdoor into a block cipher. Building on previous research, in this paper we investigate an especially dangerous algebraic structure, which is called a hidden sum and which is related to some regular subgroups of the affine group. Mixing group theory arguments and cryptographic tools, we pass from characterizing our hidden sums to designing an efficient algorithm to perform the necessary preprocessing for the exploitation of the trapdoor.


title: Runge--Kutta convolution coercivity and its use for time-dependent boundary integral equations

abstract: A coercivity property of temporal convolution operators is an essential tool in the analysis of time-dependent boundary integral equations and their space and time discretisations. It is known that this coercivity property is inherited by convolution quadrature time discretisation based on A-stable multistep methods, which are of order at most two. Here we study the question as to which Runge--Kutta-based convolution quadrature methods inherit the convolution coercivity property. It is shown that this holds without any restriction for the third-order Radau IIA method, and on permitting a shift in the Laplace domain variable, this holds for all algebraically stable Runge--Kutta methods and hence for methods of arbitrary order. As an illustration, the discrete convolution coercivity is used to analyse the stability and convergence properties of the time discretisation of a non-linear boundary integral equation that originates from a non-linear scattering problem for the linear wave equation. Numerical experiments illustrate the error behaviour of the Runge--Kutta convolution quadrature time discretisation.


title: Cathodo- and radioluminescence of Tm$^{3+}$:YAG and Nd$^{3+}$:YAG in an extended wavelength range

abstract: We have studied the cathodo- and radioluminescence of Nd:YAG and of Tm:YAG single crystals in an extended wavelength range up to $\approx 5\,\mu$m in view of developing a new kind of detector for low-energy, low-rate energy deposition events. Whereas the light yield in the visible range is as large as $\approx 10^{4}\,$photons/MeV, in good agreement with literature results, in the infrared range we have found a light yield $\approx 5\times 10^{4}\,$photons/MeV, thereby proving that ionizing radiation is particularly efficient in populating the low lying levels of rare earth doped crystals.


title: On the affine random walk on the torus

abstract: Let $\mu$ be a borelian probability measure on $\mathbf{G}:=\mathrm{SL}_d(\mathbb{Z}) \ltimes \mathbb{T}^d$. Define, for $x\in \mathbb{T}^d$, a random walk starting at $x$ denoting for $n\in \mathbb{N}$, \[ \left\{\begin{array}{rcl} X_0 &=&x\\ X_{n+1} &=& a_{n+1} X_n + b_{n+1} \end{array}\right. \] where $((a_n,b_n))\in \mathbf{G}^\mathbb{N}$ is an iid sequence of law $\mu$. Then, we denote by $\mathbb{P}_x$ the measure on $(\mathbb{T}^d)^\mathbb{N}$ that is the image of $\mu^{\otimes \mathbb{N}}$ by the map $\left((g_n) \mapsto (x,g_1 x, g_2 g_1 x, \dots , g_n \dots g_1 x, \dots)\right)$ and for any $\varphi \in \mathrm{L}^1((\mathbb{T}^d)^\mathbb{N}, \mathbb{P}_x)$, we set $\mathbb{E}_x \varphi((X_n)) = \int \varphi((X_n)) \mathrm{d}\mathbb{P}_x((X_n))$. Bourgain, Furmann, Lindenstrauss and Mozes studied this random walk when $\mu$ is concentrated on $\mathrm{SL}_d(\mathbb{Z}) \ltimes\{0\}$ and this allowed us to study, for any h\"older-continuous function $f$ on the torus, the sequence $(f(X_n))$ when $x$ is not too well approximable by rational points. In this article, we are interested in the case where $\mu$ is not concentrated on $\mathrm{SL}_d(\mathbb{Z}) \ltimes \mathbb{Q}^d/\mathbb{Z}^d$ and we prove that, under assumptions on the group spanned by the support of $\mu$, the Lebesgue's measure $\nu$ on the torus is the only stationary probability measure and that for any h\"older-continuous function $f$ on the torus, $\mathbb{E}_x f(X_n)$ converges exponentially fast to $\int f\mathrm{d}\nu$. Then, we use this to prove the law of large numbers, a non-concentration inequality, the functional central limit theorem and it's almost-sure version for the sequence $(f(X_n))$. In the appendix, we state a non-concentration inequality for products of random matrices without any irreducibility assumption.


title: Political Homophily in Independence Movements: Analysing and Classifying Social Media Users by National Identity

abstract: Social media and data mining are increasingly being used to analyse political and societal issues. Here we undertake the classification of social media users as supporting or opposing ongoing independence movements in their territories. Independence movements occur in territories whose citizens have conflicting national identities; users with opposing national identities will then support or oppose the sense of being part of an independent nation that differs from the officially recognised country. We describe a methodology that relies on users' self-reported location to build large-scale datasets for three territories -- Catalonia, the Basque Country and Scotland. An analysis of these datasets shows that homophily plays an important role in determining who people connect with, as users predominantly choose to follow and interact with others from the same national identity. We show that a classifier relying on users' follow networks can achieve accurate, language-independent classification performances ranging from 85% to 97% for the three territories.


title: Equivariance Through Parameter-Sharing

abstract: We propose to study equivariance in deep neural networks through parameter symmetries. In particular, given a group $\mathcal{G}$ that acts discretely on the input and output of a standard neural network layer $\phi_{W}: \Re^{M} \to \Re^{N}$, we show that $\phi_{W}$ is equivariant with respect to $\mathcal{G}$-action iff $\mathcal{G}$ explains the symmetries of the network parameters $W$. Inspired by this observation, we then propose two parameter-sharing schemes to induce the desirable symmetry on $W$. Our procedures for tying the parameters achieve $\mathcal{G}$-equivariance and, under some conditions on the action of $\mathcal{G}$, they guarantee sensitivity to all other permutation groups outside $\mathcal{G}$.


title: Parametric Analysis of Cherenkov Light LDF from EAS for High Energy Gamma Rays and Nuclei: Ways of Practical Application

abstract: In this paper we propose a 'knee-like' approximation of the lateral distribution of the Cherenkov light from extensive air showers in the energy range 30-3000 TeV and study a possibility of its practical application in high energy ground-based gamma-ray astronomy experiments (in particular, in TAIGA-HiSCORE). The approximation has a very good accuracy for individual showers and can be easily simplified for practical application in the HiSCORE wide angle timing array in the condition of a limited number of triggered stations.


title: Economic inequality and mobility for stochastic models with multiplicative noise

abstract: In this article, we discuss a dynamical stochastic model that represents the time evolution of income distribution of a population, where the dynamics develop from an interplay of multiple economic exchanges in the presence of multiplicative noise. The model remit stretches beyond the conventional framework of a Langevin-type kinetic equation in that our model dynamics is self-consistently constrained by dynamical conservation laws emerging from population and wealth conservation. This model is numerically solved and analyzed to interpret the inequality of income as a function of relevant dynamical parameters like the {\it mobility} $M$ and the {\it total income} $\mu$. In our model, inequality is quantified by the {\it Gini index} $G$. In particular, correlations between any two of the mobility index $M$ and/or the total income $\mu$ with the Gini index $G$ are investigated and compared with the analogous correlations resulting from an equivalent additive noise model. Our findings highlight the importance of a multiplicative noise based economic modeling structure in the analysis of inequality. The model also depicts the nature of correlation between mobility and total income of a population from the perspective of inequality measure.


title: Combining the $k$-CNF and XOR Phase-Transitions

abstract: The runtime performance of modern SAT solvers on random $k$-CNF formulas is deeply connected with the 'phase-transition' phenomenon seen empirically in the satisfiability of random $k$-CNF formulas. Recent universal hashing-based approaches to sampling and counting crucially depend on the runtime performance of SAT solvers on formulas expressed as the conjunction of both $k$-CNF and XOR constraints (known as $k$-CNF-XOR formulas), but the behavior of random $k$-CNF-XOR formulas is unexplored in prior work. In this paper, we present the first study of the satisfiability of random $k$-CNF-XOR formulas. We show empirical evidence of a surprising phase-transition that follows a linear trade-off between $k$-CNF and XOR constraints. Furthermore, we prove that a phase-transition for $k$-CNF-XOR formulas exists for $k = 2$ and (when the number of $k$-CNF constraints is small) for $k > 2$.


title: Radio emission and mass loss rate limits of four young solar-type stars

abstract: Observations of free-free continuum radio emission of four young main-sequence solar-type stars (EK Dra, Pi1 UMa, Chi1 Ori; and Kappa1 Cet) are studied to detect stellar winds or at least to place upper limits on their thermal radio emission, which is dominated by the ionized wind. These stars are excellent proxies for representing the young Sun. Upper limits on mass loss rates are calculated using their observational radio emission. Our aim is to re-examine the faint young Sun paradox by assuming that the young Sun was more massive in its past, and hence to find a possible solution for this famous problem. The observations of our sample are performed with the Karl G. Jansky VLA with excellent sensitivity, using the C-band and the Ku-band. ALMA observations are performed at 100 GHz. For the estimation of the mass loss limits, spherically symmetric winds and stationary, anisotropic, ionized winds are assumed. We compare our results to 1) mass loss rate estimates of theoretical rotational evolution models, and 2) to results of the indirect technique of determining mass loss rates: Lyman-alpha absorption. We are able to derive the most stringent direct upper limits on mass loss so far from radio observations. Two objects, EK Dra and Chi1 Ori, are detected at 6 and 14 GHz down to an excellent noise level. These stars are very active and additional radio emission identified as non-thermal emission was detected, but limits for the mass loss rates of these objects are still derived. The stars Pi1 UMa and Kappa1 Cet were not detected in either C-band or in Ku-band. For these objects we give upper limits to their radio free-free emission and calculate upper limits to their mass loss rates. Finally, we reproduce the evolution of the Sun and derive an estimate for the solar mass of the Sun at a younger age.


title: Magnetization reversal by superconducting current in $\varphi_0$ Josephson junctions

abstract: We study magnetization reversal in a $\varphi_0$ Josephson junction with direct coupling between magnetic moment and Josephson current. Our simulations of magnetic moment dynamics show that by applying an electric current pulse, we can realize the full magnetization reversal. We propose different protocols of full magnetization reversal based on the variation of the Josephson junction and pulse parameters, particularly, electric current pulse amplitude, damping of magnetization and spin-orbit interaction. We discuss experiments which can probe the magnetization reversal in $\varphi_0$-junctions.


title: Backhaul-aware Robust 3D Drone Placement in 5G+ Wireless Networks

abstract: Using drones as flying base stations is a promising approach to enhance the network coverage and area capacity by moving supply towards demand when required. However deployment of such base stations can face some restrictions that need to be considered. One of the limitations in drone base stations (drone-BSs) deployment is the availability of reliable wireless backhaul link. This paper investigates how different types of wireless backhaul offering various data rates would affect the number of served users. Two approaches, namely, network-centric and user-centric, are introduced and the optimal 3D backhaul-aware placement of a drone-BS is found for each approach. To this end, the total number of served users and sum-rates are maximized in the network-centric and user-centric frameworks, respectively. Moreover, as it is preferred to decrease drone-BS movements to save more on battery and increase flight time and to reduce the channel variations, the robustness of the network is examined as how sensitive it is with respect to the users displacements.


title: Learning Hierarchical Features from Generative Models

abstract: Deep neural networks have been shown to be very successful at learning feature hierarchies in supervised learning tasks. Generative models, on the other hand, have benefited less from hierarchical models with multiple layers of latent variables. In this paper, we prove that hierarchical latent variable models do not take advantage of the hierarchical structure when trained with existing variational methods, and provide some limitations on the kind of features existing models can learn. Finally we propose an alternative architecture that do not suffer from these limitations. Our model is able to learn highly interpretable and disentangled hierarchical features on several natural image datasets with no task specific regularization or prior knowledge.


title: Forward Event-Chain Monte Carlo: Fast sampling by randomness control in irreversible Markov chains

abstract: Irreversible and rejection-free Monte Carlo methods, recently developed in Physics under the name Event-Chain and known in Statistics as Piecewise Deterministic Monte Carlo (PDMC), have proven to produce clear acceleration over standard Monte Carlo methods, thanks to the reduction of their random-walk behavior. However, while applying such schemes to standard statistical models, one generally needs to introduce an additional randomization for sake of correctness. We propose here a new class of Event-Chain Monte Carlo methods that reduces this extra-randomization to a bare minimum. We compare the efficiency of this new methodology to standard PDMC and Monte Carlo methods. Accelerations up to several magnitudes and reduced dimensional scalings are exhibited.


title: McGan: Mean and Covariance Feature Matching GAN

abstract: We introduce new families of Integral Probability Metrics (IPM) for training Generative Adversarial Networks (GAN). Our IPMs are based on matching statistics of distributions embedded in a finite dimensional feature space. Mean and covariance feature matching IPMs allow for stable training of GANs, which we will call McGan. McGan minimizes a meaningful loss between distributions.


title: Dynamic principle for ensemble control tools

abstract: Dynamical equations describing physical systems at statistical equilibrium are commonly extended by mathematical tools called "thermostats". These tools are designed for sampling ensembles of statistical mechanics. We propose a dynamic principle for derivation of stochastic and deterministic thermostats. It is based on fundamental physical assumptions such that the canonical measure is invariant for the thermostat dynamics. This is a clear advantage over a range of recently proposed and widely discussed in the literature mathematical thermostat schemes. Following justification of the proposed principle we show its generality and usefulness for modeling a wide range of natural systems.


title: Asymmetric Tri-training for Unsupervised Domain Adaptation

abstract: Deep-layered models trained on a large number of labeled samples boost the accuracy of many tasks. It is important to apply such models to different domains because collecting many labeled samples in various domains is expensive. In unsupervised domain adaptation, one needs to train a classifier that works well on a target domain when provided with labeled source samples and unlabeled target samples. Although many methods aim to match the distributions of source and target samples, simply matching the distribution cannot ensure accuracy on the target domain. To learn discriminative representations for the target domain, we assume that artificially labeling target samples can result in a good representation. Tri-training leverages three classifiers equally to give pseudo-labels to unlabeled samples, but the method does not assume labeling samples generated from a different domain.In this paper, we propose an asymmetric tri-training method for unsupervised domain adaptation, where we assign pseudo-labels to unlabeled samples and train neural networks as if they are true labels. In our work, we use three networks asymmetrically. By asymmetric, we mean that two networks are used to label unlabeled target samples and one network is trained by the samples to obtain target-discriminative representations. We evaluate our method on digit recognition and sentiment analysis datasets. Our proposed method achieves state-of-the-art performance on the benchmark digit recognition datasets of domain adaptation.


title: Super Congruences Involving Multiple Harmonic Sums and Bernoulli Numbers

abstract: Let $m$, $r$ and $n$ be positive integers. We denote by ${\bf k}\vdash n$ any tuple of odd positive integers ${\bf k}=(k_1,\dots,k_t)$ such that $k_1+\dots+k_t=n$ and $k_j\ge 3$ for all $j$. In this paper we prove that for every sufficiently large prime $p$ $$ \sum_{\substack{l_1+l_2+\cdots+l_n=mp^r p\nmid l_1 l_2 \cdots l_n }} \frac1{l_1l_2\cdots l_n} \equiv p^{r-1} \sum_{{\bf k}\vdash n} C_{m,{\bf k}} B_{p-{\bf k}} \pmod{p^r} $$ where $B_{p-{\bf k}}=B_{p-k_1}B_{p-k_2}\cdots B_{p-k_t}$ are products of Bernoulli numbers and the coefficients $C_{m,{\bf k}}$ are polynomials of $m$ independent of $p$ and $r$. This generalizes previous results by many different authors and confirms a conjecture by the authors and their collaborators.


title: A Mutually-Dependent Hadamard Kernel for Modelling Latent Variable Couplings

abstract: We introduce a novel kernel that models input-dependent couplings across multiple latent processes. The pairwise joint kernel measures covariance along inputs and across different latent signals in a mutually-dependent fashion. A latent correlation Gaussian process (LCGP) model combines these non-stationary latent components into multiple outputs by an input-dependent mixing matrix. Probit classification and support for multiple observation sets are derived by Variational Bayesian inference. Results on several datasets indicate that the LCGP model can recover the correlations between latent signals while simultaneously achieving state-of-the-art performance. We highlight the latent covariances with an EEG classification dataset where latent brain processes and their couplings simultaneously emerge from the model.


title: A Modular Andre-Oort Statement with Derivatives

abstract: In unpublished notes, Pila discussed some theory surrounding the modular function $j$ and its derivatives. A focal point of these notes was the statement of two conjectures regarding $j$, $j'$ and $j"$: a Zilber-Pink type statement incorporating $j$, $j'$ and $j"$, which was an extension of an apparently weaker conjecture of Andre-Oort type. In this paper, I first cover some background regarding $j$, $j'$ and $j"$, mostly covering the work already done by Pila. Then I use a seemingly novel adaptation of the o-minimal Pila-Zannier strategy to prove a weakened version of Pila's "Modular Andre-Oort with Derivatives" conjecture. Under the assumption of a Schanuel-type conjecture, the central theorem of the paper implies Pila's conjecture in full generality, as well as a more precise statement on the same lines.


title: Higher derivatives in Type II and M-theory on Calabi-Yau threefolds

abstract: The four- and five-dimensional effective actions of Calabi-Yau threefold compactifications are derived with a focus on terms involving up to four space-time derivatives. The starting points for these reductions are the ten- and eleven-dimensional supergravity actions supplemented with the known eight-derivative corrections that have been inferred from Type II string amplitudes. The corrected background solutions are determined and the fluctuations of the Kahler structure of the compact space and the form-field background are discussed. It is concluded that the two-derivative effective actions for these fluctuations only takes the expected supergravity form if certain additional ten- and eleven-dimensional higher-derivative terms for the form-fields are included. The main results on the four-derivative terms include a detailed treatment of higher-derivative gravity coupled to Kahler structure deformations. This is supplemented by a derivation of the vector sector in reductions to five dimensions. While the general result is only given as an expansion in the fluctuations, a complete treatment of the one-Kahler modulus case is presented for both Type II theories and M-theory.


title: Game-Theoretic Semantics for ATL+ with Applications to Model Checking

abstract: We develop game-theoretic semantics (GTS) for the fragment ATL+ of the full Alternating-time Temporal Logic ATL*, essentially extending a recently introduced GTS for ATL. We first show that the new game-theoretic semantics is equivalent to the standard semantics of ATL+ (based on perfect recall strategies). We then provide an analysis, based on the new semantics, of the memory and time resources needed for model checking ATL+. Based on that, we establish that strategies that use only a very limited amount of memory suffice for ATL+. Furthermore, using the GTS we provide a new algorithm for model checking of ATL+ and identify a natural hierarchy of tractable fragments of ATL+ that extend ATL.


title: Stability of Electrodeposition at Solid-Solid Interfaces and Implications for Metal Anodes

abstract: We generalize the conditions for stable electrodeposition at isotropic solid-solid interfaces using a kinetic model which incorporates the effects of stresses and surface tension at the interface. We develop a stability diagram that shows two regimes of stability: previously known pressure-driven mechanism and a new density-driven stability mechanism that is governed by the relative density of metal in the two phases. We show that inorganic solids and solid polymers generally do not lead to stable electrodeposition, and provide design guidelines for achieving stable electrodeposition.


title: Photonic implementation of Majorana-based Berry phases

abstract: Geometric phases, generated by cyclic evolutions of quantum systems, offer an inspiring playground for advancing fundamental physics and technologies, alike. Intriguingly, the exotic statistics of anyons realised in physical systems can be interpreted as a topological version of geometric phases. However, non-Abelian statistics has not yet been demonstrated in the laboratory. Here we employ an all optical quantum system that simulates the statistical evolution of Majorana fermions. As a result we experimentally realise non-Abelian Berry phases with the topological characteristic that they are invariant under continues deformations of their control parameters. We implement a universal set of Majorana inspired gates by performing topological and non-topological evolutions and investigate their resilience against perturbative errors. Our photonic experiment, while it is not scalable, it suggests the intriguing possibility of experimentally simulating Majorana statistics with scalable technologies.


title: Current Induced Damping of Nanosized Quantum Moments in the Presence of Spin-Orbit Interaction

abstract: Motivated by the need to understand current-induced magnetization dynamics at the nanoscale, we have developed a formalism, within the framework of Keldysh Green function approach, to study the current-induced dynamics of a ferromagnetic (FM) nanoisland overlayer on a spin-orbit-coupling (SOC) Rashba plane. In contrast to the commonly employed classical micromagnetic LLG simulations the magnetic moments of the FM are treated {\it quantum mechanically}. We obtain the density matrix of the whole system consisting of conduction electrons entangled with the local magnetic moments and calculate the effective damping rate of the FM. We investigate two opposite limiting regimes of FM dynamics: (1) The precessional regime where the magnetic anisotropy energy (MAE) and precessional frequency are smaller than the exchange interactions, and (2) The local spin-flip regime where the MAE and precessional frequency are comparable to the exchange interactions. In the former case, we show that due to the finite size of the FM domain, the \textquotedblleft Gilbert damping\textquotedblright does not diverge in the ballistic electron transport regime, in sharp contrast to Kambersky's breathing Fermi surface theory for damping in metallic FMs. In the latter case, we show that above a critical bias the excited conduction electrons can switch the local spin moments resulting in demagnetization and reversal of the magnetization. Furthermore, our calculations show that the bias-induced antidamping efficiency in the local spin-flip regime is much higher than that in the rotational excitation regime.


title: Query Combinators

abstract: We introduce Rabbit, a combinator-based query language. Rabbit is designed to let data analysts and other accidental programmers query complex structured data. We combine the functional data model and the categorical semantics of computations to develop denotational semantics of database queries. In Rabbit, a query is modeled as a Kleisli arrow for a monadic container determined by the query cardinality. In this model, monadic composition can be used to navigate the database, while other query combinators can aggregate, filter, sort and paginate data; construct compound data; connect self-referential data; and reorganize data with grouping and data cube operations. A context-aware query model, with the input context represented as a comonadic container, can express query parameters and window functions. Rabbit semantics enables pipeline notation, encouraging its users to construct database queries as a series of distinct steps, each individually crafted and tested. We believe that Rabbit can serve as a practical tool for data analytics.


title: Clustering in Discrete Path Planning for Approximating Minimum Length Paths

abstract: In this paper we consider discrete robot path planning problems on metric graphs. We propose a clustering method, Gamma-Clustering for the planning graph that significantly reduces the number of feasible solutions, yet retains a solution within a constant factor of the optimal. By increasing the input parameter Gamma, the constant factor can be decreased, but with less reduction in the search space. We provide a simple polynomial- time algorithm for finding optimal Gamma-Clusters, and show that for a given Gamma, this optimal is unique. We demonstrate the effectiveness of the clustering method on traveling salesman instances, showing that for many instances we obtain significant reductions in computation time with little to no reduction in solution quality.


title: The revival of the Baldwin Effect

abstract: The idea that a genetically fixed behavior evolved from the once differential learning ability of individuals that performed the behavior is known as the Baldwin effect. A highly influential paper [Hinton G.E., Nowlan S.J., 1987. How learning can guide evolution. Complex Syst. 1, 495-502] claimed that this effect can be observed in silico, but here we argue that what was actually shown is that the learning ability is easily selected for. Then we demonstrate the Baldwin effect to happen in the in silico scenario by estimating the probability and waiting times for the learned behavior to become innate. Depending on parameter values, we find that learning can increase the chance of fixation of the learned behavior by several orders of magnitude compared with the non-learning situation.


title: Definition of the relativistic geoid in terms of isochronometric surfaces

abstract: We present a definition of the geoid that is based on the formalism of general relativity without approximations; i.e. it allows for arbitrarily strong gravitational fields. For this reason, it applies not only to the Earth and other planets but also to compact objects such as neutron stars. We define the geoid as a level surface of a time-independent redshift potential. Such a redshift potential exists in any stationary spacetime. Therefore, our geoid is well defined for any rigidly rotating object with constant angular velocity and a fixed rotation axis that is not subject to external forces. Our definition is operational because the level surfaces of a redshift potential can be realized with the help of standard clocks, which may be connected by optical fibers. Therefore, these surfaces are also called isochronometric surfaces. We deliberately base our definition of a relativistic geoid on the use of clocks since we believe that clock geodesy offers the best methods for probing gravitational fields with highest precision in the future. However, we also point out that our definition of the geoid is mathematically equivalent to a definition in terms of an acceleration potential, i.e. that our geoid may also be viewed as a level surface orthogonal to plumb lines. Moreover, we demonstrate that our definition reduces to the known Newtonian and post-Newtonian notions in the appropriate limits. As an illustration, we determine the isochronometric surfaces for rotating observers in axisymmetric static and axisymmetric stationary solutions to Einstein's vacuum field equation, with the Schwarzschild metric, the Erez-Rosen metric, the q-metric and the Kerr metric as particular examples.


title: Entanglement and squeezing in continuous-variable systems

abstract: We introduce a multi-mode squeezing coefficient to characterize entanglement in $N$-partite continuous-variable systems. The coefficient relates to the squeezing of collective observables in the $2N$-dimensional phase space and can be readily extracted from the covariance matrix. Simple extensions further permit to reveal entanglement within specific partitions of a multipartite system. Applications with nonlinear observables allow for the detection of non-Gaussian entanglement.


title: Einstein tori and crooked surfaces

abstract: In hyperbolic space, the angle of intersection and distance classify pairs of totally geodesic hyperplanes. A similar algebraic invariant classifies pairs of hyperplanes in the Einstein universe. In dimension 3, symplectic splittings of a 4-dimensional real symplectic vector space model Einstein hyperplanes and the invariant is a determinant. The classification contributes to a complete disjointness criterion for crooked surfaces in the 3-dimensional Einstein universe.


title: An SDP-Based Algorithm for Linear-Sized Spectral Sparsification

abstract: For any undirected and weighted graph $G=(V,E,w)$ with $n$ vertices and $m$ edges, we call a sparse subgraph $H$ of $G$, with proper reweighting of the edges, a $(1+\varepsilon)$-spectral sparsifier if \[ (1-\varepsilon)x^{\intercal}L_Gx\leq x^{\intercal} L_{H} x\leq (1+\varepsilon) x^{\intercal} L_Gx \] holds for any $x\in\mathbb{R}^n$, where $L_G$ and $L_{H}$ are the respective Laplacian matrices of $G$ and $H$. Noticing that $\Omega(m)$ time is needed for any algorithm to construct a spectral sparsifier and a spectral sparsifier of $G$ requires $\Omega(n)$ edges, a natural question is to investigate, for any constant $\varepsilon$, if a $(1+\varepsilon)$-spectral sparsifier of $G$ with $O(n)$ edges can be constructed in $\tilde{O}(m)$ time, where the $\tilde{O}$ notation suppresses polylogarithmic factors. All previous constructions on spectral sparsification require either super-linear number of edges or $m^{1+\Omega(1)}$ time. In this work we answer this question affirmatively by presenting an algorithm that, for any undirected graph $G$ and $\varepsilon>0$, outputs a $(1+\varepsilon)$-spectral sparsifier of $G$ with $O(n/\varepsilon^2)$ edges in $\tilde{O}(m/\varepsilon^{O(1)})$ time. Our algorithm is based on three novel techniques: (1) a new potential function which is much easier to compute yet has similar guarantees as the potential functions used in previous references; (2) an efficient reduction from a two-sided spectral sparsifier to a one-sided spectral sparsifier; (3) constructing a one-sided spectral sparsifier by a semi-definite program.


title: Transition to turbulence when the Tollmien-Schlichting and bypass routes coexist

abstract: Plane Poiseuille flow, the pressure driven flow between parallel plates, shows a route to turbulence connected with a linear instability to Tollmien-Schlichting (TS) waves, and another one, the bypass transition, that is triggered with finite amplitude perturbation. We use direct numerical simulations to explore the arrangement of the different routes to turbulence among the set of initial conditions. For plates that are a distance $2H$ apart and in a domain of width $2\pi H$ and length $2\pi H$ the subcritical instability to TS waves sets in at $Re_{c}=5815$ that extends down to $Re_{TS}\approx4884$. The bypass route becomes available above $Re_E=459$ with the appearance of three-dimensional finite-amplitude traveling waves. The bypass transition covers a large set of finite amplitude perturbations. Below $Re_c$, TS appear for a tiny set of initial conditions that grows with increasing Reynolds number. Above $Re_c$ the previously stable region becomes unstable via TS waves, but a sharp transition to the bypass route can still be identified. Both routes lead to the same turbulent in the final stage of the transition, but on different time scales. Similar phenomena can be expected in other flows where two or more routes to turbulence compete.


title: Strong couplings and form factors of charmed mesons in holographic QCD

abstract: We extend the two-flavor hard-wall holographic model of Erlich, Katz, Son and Stephanov [Phys.\ Rev.\ Lett.\ {\bf 95}, 261602 (2005)] to four flavors to incorporate strange and charm quarks. The model incorporates chiral and flavor symmetry breaking and provides a reasonable description of masses and weak decay constants of a variety of scalar, pseudoscalar, vector and axial-vector strange and charmed mesons. In particular, we examine flavor symmetry breaking in the strong couplings of the $\rho$ meson to the charmed $D$ and $D^{*}$ mesons. We also compute electromagnetic form factors of the $\pi$, $\rho$, $K$, $K^*$, $D$ and $D^*$ mesons. We compare our results for the $D$ and $D^*$ mesons with lattice QCD data and other nonperturbative approaches.


title: Effect of nuclear matter incompressibility on the 16O+208Pb system

abstract: To analyze the property of nuclear matter in the 16O+208Pb collision system, the internuclear potential of the fusion reaction is calculated by using the Skyrme forces associated with an extensive nuclear matter incompressibility K range in the semiclassical energy density formalism. Comparison of the experimental fusion cross sections and those obtained by using potentials derived from different forces with various K values shows that the incompressibility of nuclear matter changes during the fusion process at different bombarding energies. The results indicate that, as the energy increases, the nuclear matter becomes more incompressible.


title: The effect of deformation and vibration on the alpha decay half-life

abstract: In this work, we expand upon our previous study, the effect of surface vibrations (low-lying vibrational states) on the calculation of penetration probability in alpha decay of spherical isotopes [Phys.Rev.C91,034611(2015)]. To this aim, the Coulomb and proximity potential model, taking in to account the ground state deformations of the involved nuclei along with the surface vibrations in the daughter nucleus, is used to evaluate the alpha decay probability. The results are compared with those obtained by spherical potential barrier, which shows the dramatic effect of employing the ground state deformations in case of deformed nuclei. As well, including of surface vibrations give rise to an increase in the value of tunneling probability in better agreement with experimental data.


title: Embarrassingly Parallel Inference for Gaussian Processes

abstract: Training Gaussian process-based models typically involves an $ O(N^3)$ computational bottleneck due to inverting the covariance matrix. Popular methods for overcoming this matrix inversion problem cannot adequately model all types of latent functions, and are often not parallelizable. However, judicious choice of model structure can ameliorate this problem. A mixture-of-experts model that uses a mixture of $K$ Gaussian processes offers modeling flexibility and opportunities for scalable inference. Our embarrassingly parallel algorithm combines low-dimensional matrix inversions with importance sampling to yield a flexible, scalable mixture-of-experts model that offers comparable performance to Gaussian process regression at a much lower computational cost.


title: The role of the observer in goal-directed behavior

abstract: In goal-directed behavior, a large number of possible initial states end up in the pursued goal. The accompanying information loss implies that goal-oriented behavior is in one-to-one correspondence with an open subsystem whose entropy decreases in time. Yet ultimately, the laws of physics are reversible, so systems capable of yielding goal-directed behavior must transfer the information about initial conditions to other degrees of freedom outside the boundaries of the agent. To operate steadily, they must consume ordered degrees of freedom provided as input, and be dispensed of disordered outputs that act as wastes from the point of view of the aimed objective. Here I argue that a physical system may or may not display goal-directed behavior depending on what exactly is defined as the agent. The borders of the agent must be carefully tailored so as to entail the appropriate information balance sheet. In this game, observers play the role of tailors: They design agents by setting the limits of the system of interest. Brain-guided subjects perform this creative observation task naturally, implying that the observation of goal-oriented behavior is a goal-oriented behavior in itself. Minds evolved to cut out pieces of reality and endow them with intentionality, because ascribing intentionality is an efficient way of modeling the world, and making predictions. One most remarkable agent of whom we have indisputable evidence of its goal-pursuing attitude is the self. Notably, this agent is simultaneously the subject and the object of observation.


title: Relating correlation measures: the importance of the energy gap

abstract: The concept of correlation is central to all approaches that attempt the description of many-body effects in electronic systems. Multipartite correlation is a quantum information theoretical property that is attributed to quantum states independent of the underlying physics. In quantum chemistry, however, the correlation energy (the energy not seized by the Hartree-Fock ansatz) plays a more prominent role. We show that these two different viewpoints on electron correlation are closely related. The key ingredient turns out to be the energy gap within the symmetry-adapted subspace. We then use a few-site Hubbard model and the stretched H$_2$ to illustrate this connection and to show how the corresponding measures of correlation compare.


title: Age Progression/Regression by Conditional Adversarial Autoencoder

abstract: "If I provide you a face image of mine (without telling you the actual age when I took the picture) and a large amount of face images that I crawled (containing labeled faces of different ages but not necessarily paired), can you show me what I would look like when I am 80 or what I was like when I was 5?" The answer is probably a "No." Most existing face aging works attempt to learn the transformation between age groups and thus would require the paired samples as well as the labeled query image. In this paper, we look at the problem from a generative modeling perspective such that no paired samples is required. In addition, given an unlabeled image, the generative model can directly produce the image with desired age attribute. We propose a conditional adversarial autoencoder (CAAE) that learns a face manifold, traversing on which smooth age progression and regression can be realized simultaneously. In CAAE, the face is first mapped to a latent vector through a convolutional encoder, and then the vector is projected to the face manifold conditional on age through a deconvolutional generator. The latent vector preserves personalized face features (i.e., personality) and the age condition controls progression vs. regression. Two adversarial networks are imposed on the encoder and generator, respectively, forcing to generate more photo-realistic faces. Experimental results demonstrate the appealing performance and flexibility of the proposed framework by comparing with the state-of-the-art and ground truth.


title: Deformations, Moduli Stabilisation and Gauge Couplings at One-Loop

abstract: We investigate deformations of $\mathbb{Z}_2$ orbifold singularities on the toroidal orbifold $T^6/(\mathbb{Z}_2\times\mathbb{Z}_6)$ with discrete torsion in the framework of Type IIA orientifold model building with intersecting D6-branes wrapping special Lagrangian cycles. To this aim, we employ the hypersurface formalism developed previously for the orbifold $T^6/(\mathbb{Z}_2\times\mathbb{Z}_2)$ with discrete torsion and adapt it to the $\mathbb{Z}_2\times\mathbb{Z}_6\times\Omega\mathcal{R}$ point group by modding out the remaining $\mathbb{Z}_3$ subsymmetry and the orientifold projection $\Omega\mathcal{R}$. We first study the local behaviour of the $\mathbb{Z}_3\times\Omega\mathcal{R}$ invariant deformation orbits under non-zero deformation and then develop methods to assess the deformation effects on the fractional three-cycle volumes globally. We confirm that D6-branes supporting USp(2N) or SO(2N) gauge groups do not constrain any deformation, while deformation parameters associated to cycles wrapped by D6-branes with U(N) gauge groups are constrained by D-term supersymmetry breaking. These features are exposed in global prototype MSSM, Left-Right symmetric and Pati-Salam models first constructed in arXiv:1509.00048 and arXiv:1409.1236, for which we here count the number of stabilised moduli and study flat directions changing the values of some gauge couplings. Finally, we confront the behaviour of tree-level gauge couplings under non-vanishing deformations along flat directions with the one-loop gauge threshold corrections at the orbifold point and discuss phenomenological implications, in particular on possible LARGE volume scenarios and the corresponding value of the string scale $M_{\text{string}}$, for the same global D6-brane models.


title: xSDK Foundations: Toward an Extreme-scale Scientific Software Development Kit

abstract: Extreme-scale computational science increasingly demands multiscale and multiphysics formulations. Combining software developed by independent groups is imperative: no single team has resources for all predictive science and decision support capabilities. Scientific libraries provide high-quality, reusable software components for constructing applications with improved robustness and portability. However, without coordination, many libraries cannot be easily composed. Namespace collisions, inconsistent arguments, lack of third-party software versioning, and additional difficulties make composition costly. The Extreme-scale Scientific Software Development Kit (xSDK) defines community policies to improve code quality and compatibility across independently developed packages (hypre, PETSc, SuperLU, Trilinos, and Alquimia) and provides a foundation for addressing broader issues in software interoperability, performance portability, and sustainability. The xSDK provides turnkey installation of member software and seamless combination of aggregate capabilities, and it marks first steps toward extreme-scale scientific software ecosystems from which future applications can be composed rapidly with assured quality and scalability.


title: Kac-Moody symmetric spaces

abstract: In the present article we introduce and study a class of topological reflection spaces that we call Kac-Moody symmetric spaces. These generalize Riemannian symmetric spaces of non-compact type. We observe that in a non-spherical Kac-Moody symmetric space there exist pairs of points that do not lie on a common geodesic; however, any two points can be connected by a chain of geodesic segments. We moreover classify maximal flats in Kac-Moody symmetric spaces and study their intersection patterns, leading to a classification of global and local automorphisms. Unlike Riemannian symmetric spaces, non-spherical non-affine irreducible Kac-Moody symmetric spaces also admit an invariant causal structure. For causal and anti-causal geodesic rays with respect to this structure we find a notion of asymptoticity, which allows us to define a future and past boundary of such Kac-Moody symmetric space. We show that these boundaries carry a natural polyhedral structure and are cellularly isomorphic to the halves of the geometric realization of the twin buildings of the underlying split real Kac-Moody group. We also show that every automorphism of the symmetric space is uniquely determined by the induced cellular automorphism of the future and past boundary. The invariant causal structure on a non-spherical non-affine irreducible Kac-Moody symmetric space gives rise to an invariant pre-order on the underlying space, and thus to a subsemigroup of the Kac-Moody group. We conclude that while in some aspects Kac-Moody symmetric spaces closely resemble Riemannian symmetric spaces, in other aspects they behave similarly to ordered affine hovels, their non-Archimedean cousins.


title: Objectivity in non-Markovian spin-boson model

abstract: Objectivity constitutes one of the main features of the macroscopic classical world. An important aspect of the quantum-to-classical transition issue is to explain how such a property arises from the microscopic quantum world. Recently, within the framework of open quantum systems, there has been proposed such a mechanism in terms of the, so-called, Spectrum Broadcast Structures. These are multipartite quantum states of the system of interest and a part of its environment, assumed to be under an observation. This approach requires a departure from the standard open quantum systems methods, as the environment cannot be completely neglected. In the present work we study the emergence of such a state-structures in one of the canonical models of the condensed matter theory: Spin-boson model, describing the dynamics of a two-level system coupled to an environment made up by a large number of harmonic oscillators. We pay much attention to the behavior of the model in the non-Markovian regime, in order to provide a testbed to analyze how the non-Markovian nature of the evolution affects the surfacing of a spectrum broadcast structure.


title: A purity theorem for configuration spaces of smooth compact algebraic varieties

abstract: B. Totaro showed \cite{totaro} that the rational cohomology of configuration spaces of smooth complex projective varieties is isomorphic as an algebra to the $E_2$ term of the Leray spectral sequence corresponding to the open embedding of the configuration space into the Cartesian power. In this note we show that the isomorphism can be chosen to be compatible with the mixed Hodge structures. In particular, we prove that the mixed Hodge structures on the configuration spaces of smooth complex projective varieties are direct sums of pure Hodge structures.


title: HARPO: 1.7 - 74 MeV gamma-ray beam validation of a high angular resolutio n, high linear polarisation dilution, gas time projection chamber telescope and polarimeter

abstract: A presentation at the SciNeGHE conference of the past achievements, of the present activities and of the perspectives for the future of the HARPO project, the development of a time projection chamber as a high-performance gamma-ray telescope and linear polarimeter in the e+e- pair creation regime.


title: keV sterile neutrino Dark Matter

abstract: We give an overview of the current status of keV sterile neutrino Dark Matter. After a short introduction, we start by a general discussion of non-thermal production of Dark Matter, which applies to the three most commonly discussed mechanisms to produce sterile neutrino Dark Matter in the Universe: non-resonant, resonant, and decay production. The main goal in each case is to compute the momentum distribution function $f(p,t)$, which incorporates the full information about the Dark Matter setting under consideration, at least in what concerns its cosmological aspects. While some difficulties lie in the actual computation of this quantity, it is decisive to obtain bounds from cosmic structure formation, which turn out to be the most crucial ones to distinguish different types of production. We will introduce these bounds and we put the resulting limits into a proper context, thereby illustrating that a significant amount of relevant parameter space is available, part of which is testable in particular by Lyman-$\alpha$ data.


title: Boundary-Seeking Generative Adversarial Networks

abstract: Generative adversarial networks (GANs) are a learning framework that rely on training a discriminator to estimate a measure of difference between a target and generated distributions. GANs, as normally formulated, rely on the generated samples being completely differentiable w.r.t. the generative parameters, and thus do not work for discrete data. We introduce a method for training GANs with discrete data that uses the estimated difference measure from the discriminator to compute importance weights for generated samples, thus providing a policy gradient for training the generator. The importance weights have a strong connection to the decision boundary of the discriminator, and we call our method boundary-seeking GANs (BGANs). We demonstrate the effectiveness of the proposed algorithm with discrete image and character-based natural language generation. In addition, the boundary-seeking objective extends to continuous data, which can be used to improve stability of training, and we demonstrate this on Celeba, Large-scale Scene Understanding (LSUN) bedrooms, and Imagenet without conditioning.


title: Multiphase Porous Electrode Theory

abstract: Porous electrode theory, pioneered by John Newman and collaborators, provides a useful macroscopic description of battery cycling behavior, rooted in microscopic physical models rather than empirical circuit approximations. The theory relies on a separation of length scales to describe transport in the electrode coupled to intercalation within small active material particles. Typically, the active materials are described as solid solution particles with transport and surface reactions driven by concentration fields, and the thermodynamics are incorporated through fitting of the open circuit potential. This approach has fundamental limitations, however, and does not apply to phase-separating materials, for which the voltage is an emergent property of inhomogeneous concentration profiles, even in equilibrium. Here, we present a general theoretical framework for "multiphase porous electrode theory" implemented in an open-source software package called "MPET", based on electrochemical nonequilibrium thermodynamics. Cahn-Hilliard-type phase field models are used to describe the solid active materials with suitably generalized models of interfacial reaction kinetics. Classical concentrated solution theory is implemented for the electrolyte phase, and Newman's porous electrode theory is recovered in the limit of solid-solution active materials with Butler-Volmer kinetics. More general, quantum-mechanical models of Faradaic reactions are also included, such as Marcus-Hush-Chidsey kinetics for electron transfer at metal electrodes, extended for concentrated solutions. The full equations and numerical algorithms are described, and a variety of example calculations are presented to illustrate the novel features of the software compared to existing battery models.


title: Structure of martingale transports in finite dimensions

abstract: We study the structure of martingale transports in finite dimensions. We consider the family $\mathcal{M}(\mu,\nu) $ of martingale measures on $\mathbb{R}^N \times \mathbb{R}^N$ with given marginals $\mu,\nu$, and construct a family of relatively open convex sets $\{C_x:x\in \mathbb{R}^N \}$, which forms a partition of $\mathbb{R}^N$, and such that any martingale transport in $\mathcal{M}(\mu,\nu) $ sends mass from $x$ to within $\overline{C_x}$, $\mu(dx)$--a.e. Our results extend the analogous one-dimensional results of M. Beiglb\"ock and N. Juillet (2016) and M. Beiglb\"ock, M. Nutz, and N. Touzi (2015). We conjecture that the decomposition is canonical and minimal in the sense that it allows to characterise the martingale polar sets, i.e. the sets which have zero mass under all measures in $\mathcal{M}(\mu,\nu)$, and offers the martingale analogue of the characterisation of transport polar sets proved in M. Beiglb\"ock, M. Goldstern, G. Maresch, and W. Schachermayer (2009).


title: Skin Lesion Classification Using Hybrid Deep Neural Networks

abstract: Skin cancer is one of the major types of cancers with an increasing incidence over the past decades. Accurately diagnosing skin lesions to discriminate between benign and malignant skin lesions is crucial to ensure appropriate patient treatment. While there are many computerised methods for skin lesion classification, convolutional neural networks (CNNs) have been shown to be superior over classical methods. In this work, we propose a fully automatic computerised method for skin lesion classification which employs optimised deep features from a number of well-established CNNs and from different abstraction levels. We use three pre-trained deep models, namely AlexNet, VGG16 and ResNet-18, as deep feature generators. The extracted features then are used to train support vector machine classifiers. In the final stage, the classifier outputs are fused to obtain a classification. Evaluated on the 150 validation images from the ISIC 2017 classification challenge, the proposed method is shown to achieve very good classification performance, yielding an area under receiver operating characteristic curve of 83.83% for melanoma classification and of 97.55% for seborrheic keratosis classification.


title: Statistical Anomaly Detection via Composite Hypothesis Testing for Markov Models

abstract: Under Markovian assumptions, we leverage a Central Limit Theorem (CLT) for the empirical measure in the test statistic of the composite hypothesis Hoeffding test so as to establish weak convergence results for the test statistic, and, thereby, derive a new estimator for the threshold needed by the test. We first show the advantages of our estimator over an existing estimator by conducting extensive numerical experiments. We find that our estimator controls better for false alarms while maintaining satisfactory detection probabilities. We then apply the Hoeffding test with our threshold estimator to detecting anomalies in two distinct applications domains: one in communication networks and the other in transportation networks. The former application seeks to enhance cyber security and the latter aims at building smarter transportation systems in cities.


title: HAWC Observations Strongly Favor Pulsar Interpretations of the Cosmic-Ray Positron Excess

abstract: Recent measurements of the Geminga and B0656+14 pulsars by the gamma-ray telescope HAWC (along with earlier measurements by Milagro) indicate that these objects generate significant fluxes of very high-energy electrons. In this paper, we use the very high-energy gamma-ray intensity and spectrum of these pulsars to calculate and constrain their expected contributions to the local cosmic-ray positron spectrum. Among models that are capable of reproducing the observed characteristics of the gamma-ray emission, we find that pulsars invariably produce a flux of high-energy positrons that is similar in spectrum and magnitude to the positron fraction measured by PAMELA and AMS-02. In light of this result, we conclude that it is very likely that pulsars provide the dominant contribution to the long perplexing cosmic-ray positron excess.


title: Least-squares Solutions of Linear Differential Equations

abstract: This study shows how to obtain least-squares solutions to initial and boundary value problems to nonhomogeneous linear differential equations with nonconstant coefficients of any order. However, without loss of generality, the approach has been applied to second order differential equations. The proposed method has two steps. The first step consists of writing a constrained expression, introduced in Ref. \cite{Mortari}, that has embedded the differential equation constraints. These expressions are given in term of a new unknown function, $g (t)$, and they satisfy the constraints, no matter what $g (t)$ is. The second step consists of expressing $g (t)$ as a linear combination of $m$ independent known basis functions, $g (t) = \mathbf{\xi}^T \mathbf{h} (t)$. Specifically, Chebyshev orthogonal polynomials of the first kind are adopted for the basis functions. This choice requires rewriting the differential equation and the constraints in term of a new independent variable, $x\in[-1, +1]$. The procedure leads to a set of linear equations in terms of the unknown coefficients vector, $\mathbf{\xi},$ that is then computed by least-squares. Numerical examples are provided to quantify the solutions accuracy for initial and boundary values problems as well as for a control-type problem, where the state is defined in one point and the costate in another point.


title: Evaluation of the non-elementary integral $\int e^{\lambda x^\alpha} dx, \alpha\ge2$, and other related integrals

abstract: A formula for the non-elementary integral $\int e^{\lambda x^\alpha} dx$ where $\alpha$ is real and greater or equal two, is obtained in terms of the confluent hypergeometric function $_1F_1$. This result is verified by directly evaluating the area under the Gaussian Bell curve, corresponding to $\alpha = 2$, using the asymptotic expression for the confluent hypergeometric function and the Fundamental Theorem of Calculus (FTC). Two different but equivalent expressions, one in terms of the confluent hypergeometric function $_1F_1$ and another one in terms of the hypergeometric function $_1F_2$, are obtained for each of these integrals, $\int \cosh(\lambda x^\alpha)dx$, $\int \sinh(\lambda x^\alpha)dx$, $\int \cos(\lambda x^\alpha)dx$ and $\int \sin(\lambda x^\alpha)dx$, $\lambda\in \mathbb{C}, \alpha\ge2$. And the hypergeometric function $_1F_2$ is expressed in terms of the confluent hypergeometric function $_1F_1$. Some of the applications of the non-elementary integral $\int e^{\lambda x^\alpha}dx,\alpha\ge2$ such as the Gaussian distribution and the Maxwell-Bortsman distribution are given.


title: The spatial statistics of turbulent dissipation rates

abstract: We study the spatial statistics of velocity gradient volatility (i,e., the energy dissipation rate) in turbulent flow. We extend the Kolmogorov-Obukhov theory but also narrow its scope. The models are log normal, with verification from finely resolved large eddy and direct numerical simulations. They are parameterized by a mean and a covariance operator. Addressing applications to large eddy simulations, the mean and the covariance depend on the resolved scale solution. Removing this resolved scale dependence by a locally defined rescaling of the turbulent statistics yields a universal theory for the subgrid statistics, specifically for the mean and variance of turbulent fluctuations, i.e., the log of the energy dissipation rate {\epsilon}. The variance of the velocity gradient statistics is found to be log normal, in accordance with Kolmogorov 1962. Rescaling by the resolved scale mean and variance removes the influence of the resolved scales on the subgrid scales, justifying conceptually the universality we observe. The universality is the basis of new power scaling laws. The new power laws, in turn, allow a simple parameterization of the mean and covariance of the subgrid rescaled dissipation rate, i.e., velocity gradient volatility. We treat the coarse grid resolved space- time location as a random variable. A restriction of the theory to a small range of unresolved but stochastically modelled scales is proposed, with renormalization group ideas offered to overcome this restriction.


title: On some $q$-versions of Ramanujan master theorem

abstract: In this paper, we state some $q$-analogues of the famous Ramanujan's Master Theorem. As applications, some values of Jackson's $q$-integrals involving $q$-special functions are computed.


title: Monte Carlo Action Programming

abstract: This paper proposes Monte Carlo Action Programming, a programming language framework for autonomous systems that act in large probabilistic state spaces with high branching factors. It comprises formal syntax and semantics of a nondeterministic action programming language. The language is interpreted stochastically via Monte Carlo Tree Search. Effectiveness of the approach is shown empirically.


title: The Isostatic Conjecture

abstract: We show that a jammed packing of disks with generic radii, in a generic container, is such that the minimal number of contacts occurs and there is only one dimension of equilibrium stresses. We also point out some connections to packings with different radii and results in the theory of circle packings whose graph forms a triangulation of a given topological surface. We also point out a counterexample, due to F. Nazarov, to a previous conjecture that that triangulated packings with fixed numbers of disks with fixed numbers of disks for each radius claiming that such packings were the most dense.


title: Elementary Yet Precise Worst-case Analysis of MergeSort, A short version (SV)

abstract: This paper offers two elementary yet precise derivations of an exact formula \[ W(n) = \sum_{i=1} ^{n} \lceil \lg i \rceil = n \lceil \lg n \rceil - 2^{\lceil \lg n \rceil} + 1 \] for the maximum number $ W(n) $ of comparisons of keys performed by $ {\tt MergeSort} $ on an $ n $-element array. The first of the two, due to its structural regularity, is well worth carefully studying in its own right. Close smooth bounds on $ W(n) $ are derived. It seems interesting that $ W(n) $ is linear between the points $ n = 2^{\lfloor \lg n \rfloor} $ and it linearly interpolates its own lower bound $ n \lg n - n + 1 $ between these points.


title: Optimal Allocation of Resources for Suppressing Epidemic Spreading on Networks

abstract: Efficient allocation of limited medical resources is crucial for controlling epidemic spreading on networks. Based on the susceptible-infected-susceptible model, we solve an optimization problem as how best to allocate the limited resources so as to minimize the prevalence, providing that the curing rate of each node is positively correlated to its medical resource. By quenched mean-field theory and heterogeneous mean-field (HMF) theory, we prove that epidemic outbreak will be suppressed to the greatest extent if the curing rate of each node is directly proportional to its degree, under which the effective infection rate $\lambda$ has a maximal threshold $\lambda_c^{opt}=1/\left\langle k \right\rangle$ where $\left\langle k \right\rangle$ is average degree of the underlying network. For weak infection region ($\lambda\gtrsim\lambda_c^{opt}$), we combine a perturbation theory with Lagrange multiplier method (LMM) to derive the analytical expression of optimal allocation of the curing rates and the corresponding minimized prevalence. For general infection region ($\lambda>\lambda_c^{opt}$), the high-dimensional optimization problem is converted into numerically solving low-dimensional nonlinear equations by the HMF theory and LMM. Counterintuitively, in the strong infection region the low-degree nodes should be allocated more medical resources than the high-degree nodes to minimize the prevalence. Finally, we use simulated annealing to validate the theoretical results.


title: Electroweak Phase Transitions in Einstein's Static Universe

abstract: We suggest using Einstein's static universe metric for the metastable state after reheating, instead of the Friedman-Robertson-Walker spacetime. In this case strong static gravitational potential leads to the effective reduction of the Higgs vacuum expectation value, which is found to be compatible with the Standard Model first order electroweak phase transition conditions. Gravity could also increase the CP-violating effects for particles that cross the new phase bubble walls and thus is able to lead to the successful electroweak baryogenesis scenario.


title: Monte Carlo on manifolds: sampling densities and integrating functions

abstract: We describe and analyze some Monte Carlo methods for manifolds in Euclidean space defined by equality and inequality constraints. First, we give an MCMC sampler for probability distributions defined by un-normalized densities on such manifolds. The sampler uses a specific orthogonal projection to the surface that requires only information about the tangent space to the manifold, obtainable from first derivatives of the constraint functions, hence avoiding the need for curvature information or second derivatives. Second, we use the sampler to develop a multi-stage algorithm to compute integrals over such manifolds. We provide single-run error estimates that avoid the need for multiple independent runs. Computational experiments on various test problems show that the algorithms and error estimates work in practice. The method is applied to compute the entropies of different sticky hard sphere systems. These predict the temperature or interaction energy at which loops of hard sticky spheres become preferable to chains.


title: Thermodynamic Limit of Interacting Particle Systems over Time-varying Sparse Random Networks

abstract: We establish a functional weak law of large numbers for observable macroscopic state variables of interacting particle systems (e.g., voter and contact processes) over fast time-varying sparse random networks of interactions. We show that, as the number of agents $N$ grows large, the proportion of agents $\left(\overline{Y}_{k}^{N}(t)\right)$ at a certain state $k$ converges in distribution -- or, more precisely, weakly with respect to the uniform topology on the space of \emph{c\`adl\`ag} sample paths -- to the solution of an ordinary differential equation over any compact interval $\left[0,T\right]$. Although the limiting process is Markov, the prelimit processes, i.e., the normalized macrostate vector processes $\left(\mathbf{\overline{Y}}^{N}(t)\right)=\left(\overline{Y}_{1}^{N}(t),\ldots,\overline{Y}_{K}^{N}(t)\right)$, are non-Markov as they are tied to the \emph{high-dimensional} microscopic state of the system, which precludes the direct application of standard arguments for establishing weak convergence. The techniques developed in the paper for establishing weak convergence might be of independent interest.


title: Simplify proposal for realizing multiqubit tunable phase gate in circuit QED

abstract: We propose a scheme to realize a multiqubit tunable phase gate in a circuit QED setup where two resonators, each coupling with a qudit, are interconnected to a common qudit (d=4). In this proposal, only two levels of each qudit serve as the logical states and the other two levels are used for the gate realization. The proposal is efficient and simple because only a classical microwave pulse is needed, no matter how many qudits are involved, which significantly reduces experimental difficulty. In a nonresonant case, the tunable phase gate can be achieved readily, while under the resonant condition, a {\pi}-phase gate can be realized after a full cycle of Rabi oscillation where the gate speed is rather fast due to the resonant interaction. We have shown that the resulting effective dynamics allows for the creation of a high-fidelity phase gate. The influence of various decoherence processes such as the decay of the resonator mode and the relaxation of the qudits is investigated. Moreover, the proposed scheme can be easily generalized to realize N-qubit phase gate.


title: First Data Release of the Hyper Suprime-Cam Subaru Strategic Program

abstract: The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

categories: astro-ph.IM astro-ph.EP astro-ph.GA astro-ph.HE astro-ph.SR


title: A Knowledge-Based Approach to Word Sense Disambiguation by distributional selection and semantic features

abstract: Word sense disambiguation improves many Natural Language Processing (NLP) applications such as Information Retrieval, Information Extraction, Machine Translation, or Lexical Simplification. Roughly speaking, the aim is to choose for each word in a text its best sense. One of the most popular method estimates local semantic similarity relatedness between two word senses and then extends it to all words from text. The most direct method computes a rough score for every pair of word senses and chooses the lexical chain that has the best score (we can imagine the exponential complexity that returns this comprehensive approach). In this paper, we propose to use a combinatorial optimization metaheuristic for choosing the nearest neighbors obtained by distributional selection around the word to disambiguate. The test and the evaluation of our method concern a corpus written in French by means of the semantic network BabelNet. The obtained accuracy rate is 78 % on all names and verbs chosen for the evaluation.


title: Approches d'analyse distributionnelle pour am\'eliorer la d\'esambigu\"isation s\'emantique

abstract: Word sense disambiguation (WSD) improves many Natural Language Processing (NLP) applications such as Information Retrieval, Machine Translation or Lexical Simplification. WSD is the ability of determining a word sense among different ones within a polysemic lexical unit taking into account the context. The most straightforward approach uses a semantic proximity measure between the word sense candidates of the target word and those of its context. Such a method very easily entails a combinatorial explosion. In this paper, we propose two methods based on distributional analysis which enable to reduce the exponential complexity without losing the coherence. We present a comparison between the selection of distributional neighbors and the linearly nearest neighbors. The figures obtained show that selecting distributional neighbors leads to better results.


title: Cooling of neutron stars in soft X-ray transients

abstract: Thermal states of neutron stars in soft X-ray transients (SXRTs) are thought to be determined by "deep crustal heating" in the accreted matter that drives the quiescent luminosity and cooling via emission of photons and neutrinos from the interior. In this study, we assume a global thermal steady-state of the transient system and calculate the heating curves (quiescent surface luminosity as a function of mean accretion rate) predicted from theoretical models, taking into account variations in the equations of state, superfluidity gaps, thickness of the light element layer and a phenomenological description of the direct Urca threshold. We further provide a statistical analysis on the uncertainties in these parameters, and compare the overall results with observations of several SXRTs, in particular the two sources containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. Interpretation of the observational data indicates that the direct Urca process is required for the most massive stars and also suggests small superfluid gaps.


title: Cosmic strings and other topological defects in nonscaling regimes

abstract: Cosmic strings are topological defects possibly formed in the early Universe, which may be observable due to their gravitational effects on the cosmic microwave background radiation or gravitational wave experiments. To this effect it is important to quantitatively ascertain the network properties, including their density, velocity or the number of strings present, at the various epochs in the observable Universe. Attempts to estimate these numbers often rely on simplistic approximations for the string parameters, such as assuming that the network is scaling. However, in cosmological models containing realistic amounts of radiation, matter and dark energy a string network is never exactly scaling. Here we use the velocity-dependent one-scale model for the evolution of a string network to better quantify how these networks evolve. In particular we obtain new approximate analytic solutions for the behavior of the network during the radiation-to-matter and matter-to-acceleration transitions (assuming, in the latter case, the canonical $\Lambda$ cold dark matter model), and numerically calculate the relevant quantities for a range of possible dark energy models.


title: Composite Reflective/Absorptive IR-Blocking Filters Embedded in Metamaterial Antireflection Coated Silicon

abstract: Infrared (IR) blocking filters are crucial for controlling the radiative loading on cryogenic systems and for optimizing the sensitivity of bolometric detectors in the far-IR. We present a new IR filter approach based on a combination of patterned frequency selective structures on silicon and a thin (50 $\mu \textrm{m}$ thick) absorptive composite based on powdered reststrahlen absorbing materials. For a 300 K blackbody, this combination reflects $\sim$50\% of the incoming light and blocks \textgreater 99.8\% of the total power with negligible thermal gradients and excellent low frequency transmission. This allows for a reduction in the IR thermal loading to negligible levels in a single cold filter. These composite filters are fabricated on silicon substrates which provide excellent thermal transport laterally through the filter and ensure that the entire area of the absorptive filter stays near the bath temperature. A metamaterial antireflection coating cut into these substrates reduces in-band reflections to below 1\%, and the in-band absorption of the powder mix is below 1\% for signal bands below 750 GHz. This type of filter can be directly incorporated into silicon refractive optical elements.