title: Dispersion Analysis of Finite Difference and Discontinuous Galerkin Schemes for Maxwell's Equations in Linear Lorentz Media

abstract: In this paper, we consider Maxwell's equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes studied in our previous research [5,6]. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. The results for the numerical dispersion analysis can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwell's equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures of numerical dispersion errors. Finally, we highlight the limitations of the second order accurate temporal discretizations considered.


title: A Nonparametric Approach to High-dimensional k-sample Comparison Problems

abstract: High-dimensional k-sample comparison is a common applied problem. We construct a class of easy-to-implement nonparametric distribution-free tests based on new tools and unexplored connections with spectral graph theory. The test is shown to possess various desirable properties along with a characteristic exploratory flavor that has practical consequences. The numerical examples show that our method works surprisingly well under a broad range of realistic situations.


title: Benchmarking Nonequilibrium Green's Functions against Configuration Interaction for time-dependent Auger decay processes

abstract: We have recently proposed a Nonequilibrium Green's Function (NEGF) approach to include Auger decay processes in the ultrafast charge dynamics of photoionized molecules. Within the so called Generalized Kadanoff-Baym Ansatz the fundamental unknowns of the NEGF equations are the reduced one-particle density matrix of bound electrons and the occupations of the continuum states. Both unknowns are one-time functions like the density in Time-Dependent Functional Theory (TDDFT). In this work we assess the accuracy of the approach against Configuration Interaction (CI) calculations in one-dimensional model systems. Our results show that NEGF correctly captures qualitative and quantitative features of the relaxation dynamics provided that the energy of the Auger electron is much larger than the Coulomb repulsion between two holes in the valence shells. For the accuracy of the results dynamical electron-electron correlations or, equivalently, memory effects play a pivotal role. The combination of our NEGF approach with the Sham-Schl\"uter equation may provide useful insights for the development of TDDFT exchange-correlation potentials with a history dependence.


title: Fault Tolerant and Fully Dynamic DFS in Undirected Graphs: Simple Yet Efficient

abstract: We present an algorithm for a fault tolerant Depth First Search (DFS) Tree in an undirected graph. This algorithm is drastically simpler than the current state-of-the-art algorithms for this problem, uses optimal space and optimal preprocessing time, and still achieves better time complexity. This algorithm also leads to a better time complexity for maintaining a DFS tree in a fully dynamic environment.


title: The bispectra of galactic CMB foregrounds and their impact on primordial non-Gaussianity estimation

abstract: We use the binned bispectrum estimator to determine the bispectra of the dust, free-free, synchrotron, and AME galactic foregrounds using maps produced by the Commander component separation method from Planck 2015 data. We find that all of these peak in the squeezed configuration, allowing for potential confusion with in particular the local primordial shape. Applying an additional functionality implemented in the binned bispectrum estimator code, we then use these galactic bispectra as templates in an $f_\mathrm{NL}$ analysis of other maps. After testing and validating the method and code with simulations, we show that we detect the dust in the raw 143 GHz map with the expected amplitude (the other galactic foregrounds are too weak at 143 GHz to be detected) and that no galactic residuals are detected in the cleaned CMB map. We also investigate the effect of the mask on the templates and the effect of the choice of binning on a joint dust-primordial $f_\mathrm{NL}$ analysis.


title: BSDE Representation and Randomized Dynamic Programming Principle for Stochastic Control Problems of Infinite-Dimensional Jump-Diffusions

abstract: We consider a general class of stochastic optimal control problems, where the state process lives in a real separable Hilbert space and is driven by a cylindrical Brownian motion and a Poisson random measure; no special structure is imposed on the coefficients, which are also allowed to be path-dependent; in addition, the diffusion coefficient can be degenerate. For such a class of stochastic control problems, we prove, by means of purely probabilistic techniques based on the so-called randomization method, that the value of the control problem admits a probabilistic representation formula (known as non-linear Feynman-Kac formula) in terms of a suitable backward stochastic differential equation. This probabilistic representation considerably extends current results in the literature on the infinite-dimensional case, and it is also relevant in finite dimension. Such a representation allows to show, in the non-path-dependent (or Markovian) case, that the value function satisfies the so-called randomized dynamic programming principle. As a consequence, we are able to prove that the value function is a viscosity solution of the corresponding Hamilton-Jacobi-Bellman equation, which turns out to be a second-order fully non-linear integro-differential equation in Hilbert space.


title: Can everyday AI be ethical. Fairness of Machine Learning Algorithms

abstract: Combining big data and machine learning algorithms, the power of automatic decision tools induces as much hope as fear. Many recently enacted European legislation (GDPR) and French laws attempt to regulate the use of these tools. Leaving aside the well-identified problems of data confidentiality and impediments to competition, we focus on the risks of discrimination, the problems of transparency and the quality of algorithmic decisions. The detailed perspective of the legal texts, faced with the complexity and opacity of the learning algorithms, reveals the need for important technological disruptions for the detection or reduction of the discrimination risk, and for addressing the right to obtain an explanation of the auto- matic decision. Since trust of the developers and above all of the users (citizens, litigants, customers) is essential, algorithms exploiting personal data must be deployed in a strict ethical framework. In conclusion, to answer this need, we list some ways of controls to be developed: institutional control, ethical charter, external audit attached to the issue of a label.


title: Submodular Stochastic Probing with Prices

abstract: We introduce Stochastic Probing with Prices (SPP), a variant of the Stochastic Probing (SP) model in which we must pay a price to probe an element. A SPP problem involves two set systems $(N,\mathcal{I}_{in})$ and $(N,\mathcal{I}_{out})$ where each $e\in N$ is active with probability $p_e$. To discover whether $e$ is active, it must be probed by paying the price $\Delta_e$. If it is probed and active, then it is irrevocably added to the solution. Moreover, at all times, the set of probed elements must lie in $\mathcal{I}_{out}$, and the solution (the set of probed and active elements) must lie in $\mathcal{I}_{in}$. The goal is to maximize a set function $f$ minus the cost of the probes. We give a bi-criteria approximation algorithm to the online version of this problem, in which the elements are shown to the algorithm in a possibly adversarial order. Our results translate to state-of-the-art approximations for the traditional (online) stochastic probing problem.


title: Judiciously 3-partitioning 3-uniform hypergraphs

abstract: Bollob\'as, Reed and Thomason proved every $3$-uniform hypergraph $\mathcal{H}$ with $m$ edges has a vertex-partition $V(\mathcal{H})=V_1 \sqcup V_2 \sqcup V_3$ such that each part meets at least $\frac{1}{3}(1-\frac{1}{e})m$ edges, later improved to $0.6m$ by Halsegrave and improved asymptotically to $0.65m+o(m)$ by Ma and Yu. We improve this asymptotic bound to $\frac{19}{27}m+o(m)$, which is best possible up to the error term, resolving a special case of a conjecture of Bollob\'as and Scott.


title: 2018 Low-Power Image Recognition Challenge

abstract: The Low-Power Image Recognition Challenge (LPIRC, https://rebootingcomputing.ieee.org/lpirc) is an annual competition started in 2015. The competition identifies the best technologies that can classify and detect objects in images efficiently (short execution time and low energy consumption) and accurately (high precision). Over the four years, the winners' scores have improved more than 24 times. As computer vision is widely used in many battery-powered systems (such as drones and mobile phones), the need for low-power computer vision will become increasingly important. This paper summarizes LPIRC 2018 by describing the three different tracks and the winners' solutions.


title: A deep learning pipeline for product recognition on store shelves

abstract: Recognition of grocery products in store shelves poses peculiar challenges. Firstly, the task mandates the recognition of an extremely high number of different items, in the order of several thousands for medium-small shops, with many of them featuring small inter and intra class variability. Then, available product databases usually include just one or a few studio-quality images per product (referred to herein as reference images), whilst at test time recognition is performed on pictures displaying a portion of a shelf containing several products and taken in the store by cheap cameras (referred to as query images). Moreover, as the items on sale in a store as well as their appearance change frequently over time, a practical recognition system should handle seamlessly new products/packages. Inspired by recent advances in object detection and image retrieval, we propose to leverage on state of the art object detectors based on deep learning to obtain an initial productagnostic item detection. Then, we pursue product recognition through a similarity search between global descriptors computed on reference and cropped query images. To maximize performance, we learn an ad-hoc global descriptor by a CNN trained on reference images based on an image embedding loss. Our system is computationally expensive at training time but can perform recognition rapidly and accurately at test time.


title: Wind LiDAR measurement campaign at CESAR Observatory in Cabauw: preliminary results

abstract: A two-year verification campaign of the ZephIR~300 vertical profiling wind LiDAR has recently been started by the Royal Netherlands Meteorological Institute (KNMI) at the CESAR Observatory (Cabauw Experimental Site for Atmospheric Research) in Cabauw, The Netherlands. In this paper we present preliminary results of the first six months of this campaign (February 14 - July 31, 2018), focusing on the (height-dependent) data availability of the ZephIR~300 under various meteorological conditions (precipitation intensity, cloud base height and visibility) and the data quality via a comparison with in situ wind measurements at several levels in the 213-m tall meteorological mast.


title: Fractional Generalizations of Gradient Mechanics

abstract: This short chapter provides a fractional generalization of gradient mechanics, an approach (originally advanced by the author in the mid 80s) that has gained world-wide attention in the last decades due to its capability of modeling pattern forming instabilities and size effects in materials, as well as eliminating undesired elastic singularities. It is based on the incorporation of higher-order gradients (in the form of Laplacians) in the classical constitutive equations multiplied by appropriate internal lengths accounting for the geometry/topology of underlying micro/nano structures. This review will focus on the fractional generalization of the gradient elasticity equations (GradEla) an extension of classical elasticity to incorporate the Laplacian of Hookean stress by replacing the standard Laplacian by its fractional counterpart. On introducing the resulting fractional constitutive equation into the classical static equilibrium equation for the stress, a fractional differential equation is obtained whose fundamental solutions are derived by using the Greens function procedure. As an example, Kelvins problem is analyzed within the aforementioned setting. Then, an extension to consider constitutive equations for a restrictive class of nonlinear elastic deformations and deformation theory of plasticity is pursued. Finally, the methodology is applied for extending the authors higher-order diffusion theory from the integer to the fractional case.


title: Auction Theory Adaptations for Real Life Applications

abstract: We develop extensions to auction theory results that are useful in real life scenarios. 1. Since valuations are generally positive we first develop approximations using the log-normal distribution. This would be useful for many finance related auction settings since asset prices are usually non-negative. 2. We formulate a positive symmetric discrete distribution, which is likely to be followed by the total number of auction participants, and incorporate this into auction theory results. 3. We develop extensions when the valuations of the bidders are interdependent and incorporate all the results developed into a final combined realistic setting. 4. Our methods can be a practical tool for bidders and auction sellers to maximize their profits. The models developed here could be potentially useful for inventory estimation and for wholesale procurement of financial instruments and also non-financial commodities. All the propositions are new results and they refer to existing results which are stated as Lemmas.


title: Zero-one generation laws for finite simple groups

abstract: Let $G$ be a simple algebraic group over the algebraic closure of $GF(p)$ ($p$ prime), and let $G(q)$ denote a corresponding finite group of Lie type over $GF(q)$, where $q$ is a power of $p$. Let $X$ be an irreducible subvariety of $G^r$ for some $r\ge 2$. We prove a zero-one law for the probability that $G(q)$ is generated by a random $r$-tuple in $X(q) = X\cap G(q)^r$: the limit of this probability as $q$ increases (through values of $q$ for which $X$ is stable under the Frobenius morphism defining $G(q)$) is either 1 or 0. Indeed, to ensure that this limit is 1, one only needs $G(q)$ to be generated by an $r$-tuple in $X(q)$ for two sufficiently large values of $q$. We also prove a version of this result where the underlying characteristic is allowed to vary. In our main application, we apply these results to the case where $r=2$ and the irreducible subvariety $X = C\times D$, a product of two conjugacy classes of elements of finite order in $G$. This leads to new results on random $(2,3)$-generation of finite simple groups $G(q)$ of exceptional Lie type: provided $G(q)$ is not a Suzuki group, we show that the probability that a random involution and a random element of order 3 generate $G(q)$ tends to $1$ as $q \rightarrow \infty$. Combining this with previous results for classical groups, this shows that finite simple groups (apart from Suzuki groups and $PSp_4(q)$) are randomly $(2,3)$-generated. Our tools include algebraic geometry, representation theory of algebraic groups, and character theory of finite groups of Lie type.


title: ELT Spectroscopy of Extragalactic Massive Stars

abstract: Quantitative spectroscopy of O-type stars in galaxies beyond 1 Mpc requires the unprecedented sensitivity of the so-called Extremely Large Telescopes. Visible spectroscopy with these exciting new facilities will give us the first chance to explore the properties of large samples of massive stars in very metal-poor galaxies, as well as the broad range of environments in the galaxies of the Sculptor Group and beyond. In this context, we introduce performance simulations of visible spectroscopy of extragalactic massive stars with the MOSAIC instrument, which is a multi-object spectrograph in development for the European Extremely Large Telescope (ELT). Preliminary results demonstrate that we can obtain sufficient signal-to-noise from one night of observations in galaxies out to distances of 3 to 4 Mpc (V = 22-23 mag). We also highlight the goal of obtaining ultraviolet spectroscopy of the same individual stars beyond 1 Mpc; this will require a future large-aperture space telescope with UV capabilities such as the LUVOIR concept currently under study.


title: An Experimental Evaluation of the Generalized Sinusoidal Frequency Modulated Waveform for Active Sonar Systems

abstract: This paper experimentally evaluates the Generalized Sinusoidal Frequency Modulated (GSFM) waveform, a generalization of the Sinusoidal Frequency Modulated (SFM) waveform. The Instantaneous Frequency (IF) of the GSFM resembles the time/voltage characteristic of a Linear FM (LFM) chirp waveform. Consequently, the GSFM possesses an Ambiguity Function (AF) that resembles a thumbtack shape. Practical sonar system design must consider two factors beyond the AF. The spectral efficiency (SE), defined as the ratio of energy in an operational frequency band to the total waveform energy, is another important metric for waveform design. The Peak-to-Average-Power Ratio (PAPR) quantifies how close the waveform is to constant amplitude. These measures predict a waveform's energy efficiency and ability to be accurately replicated on practical piezoelectric transducers, which have limits on both their bandwidth and maximum transmit power. This paper explores these design considerations for the GSFM waveform and evaluates its performance against a host of other well established waveforms using simulated and experimental acoustic data. The GSFM possesses superior SE, PAPR, and overall energy efficiency when compared to thumbtack waveforms.


title: Functional Dynamics by Intention Recognition in Iterated Games

abstract: Intention recognition is an important characteristic of intelligent agents. In their interactions with others, they try to read others' intentions and make an image of others to choose their actions accordingly. While the way in which players choose their actions depending on such intentions has been investigated in game theory, how dynamic changes in intentions by mutually reading others' intentions are incorporated into game theory has not been explored. We present a novel formulation of game theory in which players read others' intentions and change their own through an iterated game. Here, intention is given as a function of the other's action and the own action to be taken accordingly as the dependent variable, while the mutual recognition of intention is represented as the functional dynamics. It is shown that a player suffers no disadvantage when he/she recognizes the other's intention, whereas the functional dynamics reach equilibria in which both players' intentions are optimized. These cover a classical Nash and Stackelberg equilibria but we extend them in this study: Novel equilibria exist depending on the degree of mutual recognition. Moreover, the degree to which each player recognizes the other can also differ. This formulation is applied to resource competition, duopoly, and prisoner's dilemma games. For example, in the resource competition game with player-dependent capacity on gaining the resource, the superior player's recognition leads to the exploitation of the other, while the inferior player's recognition leads to cooperation through which both players' payoffs increase.


title: Quantum Berezinskii-Kosterltz-Thouless Transition for Topological Insulator

abstract: We consider the interacting helical liquid system at the one-dimensional edge of a two-dimensional topological insulator, coupled to an external magnetic field and s-wave superconductor and map it to an XYZ spin chain system. This model undergoes quantum Berezinskii-Kosterlitz-Thouless (BKT) transition with two limiting conditions. We derive the renormalization group (RG) equations explicitly and also present the flow lines behavior. We also present the behavior of RG flow lines based on the exact solution. We observe that the physics of Majorana fermion zero modes and the gaped Ising-ferromagnetic phase, which appears in a different context. We observe that the evidence of gapless helical Luttinger liquid phase as a common non-topological quantum phase for both quantum BKT transitions. We explain analytically and physically that there is no Majorana-Ising transition. In the presence of chemical potential, the system shows the commensurate to incommensurate transition.


title: A characterization of the Edge of Criticality in Binary Echo State Networks

abstract: Echo State Networks (ESNs) are simplified recurrent neural network models composed of a reservoir and a linear, trainable readout layer. The reservoir is tunable by some hyper-parameters that control the network behaviour. ESNs are known to be effective in solving tasks when configured on a region in (hyper-)parameter space called \emph{Edge of Criticality} (EoC), where the system is maximally sensitive to perturbations hence affecting its behaviour. In this paper, we propose binary ESNs, which are architecturally equivalent to standard ESNs but consider binary activation functions and binary recurrent weights. For these networks, we derive a closed-form expression for the EoC in the autonomous case and perform simulations in order to assess their behavior in the case of noisy neurons and in the presence of a signal. We propose a theoretical explanation for the fact that the variance of the input plays a major role in characterizing the EoC.


title: The Physical Essence of Pulsar Glitch

abstract: Based on the magnetic dipole radiation from the 3P2 neutron superfluid vortices (3P2NSFV) in neutron stars, we propose a model of glitch for young pulsars by oscillation between B phase and A phase of 3P2 Neutron superfluid. The main behavior of glitches of pulsars may be naturally explained by our model.


title: A Study of the High-Luminosity Quasar HS 1946+7658

abstract: We study the variability of the quasar HS 1946+7658 on intra-night time scale based on both our own optical and archival X-ray data. We find the quasar non-variable during about 11 hours of optical monitoring. This is in accordance with the low intra-night variability duty cycle of radio-quiet quasars. Regarding the X-rays, we cannot make a firm conclusion about the quasar variability owing to the controversial results of the light curves statistical analysis. In addition, we calibrated Johnson-Cousins $BVRI$ magnitudes of 7 field stars that are to be used as secondary standards.


title: Laying the Groundwork for the Development of the Data Archive of the New Robotic Telescope

abstract: The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope which will be online in the early-2020s; and with the development of a new data archive facility for a range of telescopes at the National Astronomical Research Institute of Thailand. The Newton Fund enabled us to collaborate in designing a new versatile generic system that serves all purposes. In the end, we conclude that a single system would not meet the needs of all parties and only adopt similar front-ends while the back-ends are bespoke to our respective systems and data-flows.


title: Comment on `An elementary argument for the magnetic field outside a solenoid'

abstract: An alternative elementary argument for the magnetic field outside a solenoid is described.


title: Topological exploration of artificial neuronal network dynamics

abstract: One of the paramount challenges in neuroscience is to understand the dynamics of individual neurons and how they give rise to network dynamics when interconnected. Historically, researchers have resorted to graph theory, statistics, and statistical mechanics to describe the spatiotemporal structure of such network dynamics. Our novel approach employs tools from algebraic topology to characterize the global properties of network structure and dynamics. We propose a method based on persistent homology to automatically classify network dynamics using topological features of spaces built from various spike-train distances. We investigate the efficacy of our method by simulating activity in three small artificial neural networks with different sets of parameters, giving rise to dynamics that can be classified into four regimes. We then compute three measures of spike train similarity and use persistent homology to extract topological features that are fundamentally different from those used in traditional methods. Our results show that a machine learning classifier trained on these features can accurately predict the regime of the network it was trained on and also generalize to other networks that were not presented during training. Moreover, we demonstrate that using features extracted from multiple spike-train distances systematically improves the performance of our method.


title: Hives Determined by Pairs in the Affine Grassmannian over Discrete Valuation Rings

abstract: Let ${\mathcal O}$ be a discrete valuation ring with quotient field ${\cal K}$. The affine Grassmannian ${\cal G}r$ is the set of full-rank ${\mathcal O}$-modules contained in ${\cal K}^n$. Given $\Lambda \in {\cal G}r$, invariant factors $inv(\Lambda)=\lambda \in {\mathbb Z}^n$ stratify ${\cal G}r$. Left-multiplication by $GL_{n}({\cal K})$ stratifies ${\cal G}r \times {\cal G}r$ where $inv(N,\Lambda) = \mu$ if $(N,\Lambda)$ and $(I_{n} ,M)$ are in the same $GL_{n}({\cal K})$ orbit, and $inv(M) = \mu$. We present an elementary map from ${\cal G}r \times {\cal G}r$ to hives (in the sense of Knutson and Tao) of type $(\mu,\nu,\lambda)$ where $inv(N,\Lambda) = \mu$, $inv(N) = \nu$, and $inv(\Lambda) = \lambda$. Earlier work by the authors determined Littlewood-Richardson fillings from matrix pairs over certain rings ${\mathcal O}$, and later Kamnitzer utilized properties of MV polytopes to define a map from ${\cal G}r\times {\cal G}r$ to hives over ${\mathcal O} = {\mathbb C}[[t]]$. Our proof uses only linear algebra methods over any discrete valuation ring, where hive entries are minima of sums of orders of invariant factors over certain submodules. Our map is analogous to a conjectured construction of hives from Hermitian matrix pairs due to Danilov and Koshevoy.


title: Highly confined Love waves modes by defect states in a holey SiO2 /quartz phononic crystal

abstract: Highly confined Love modes are demonstrated in a phononic crystal based on a square array of etched holes in SiO 2 deposited on the ST-cut quartz. An optimal choice of the geometrical parameters contributes to a wide stop-band for shear waves' modes. The introduction of a defect by removing lines of holes leads to the nearly flat modes within the band gap and consequently paves the way to implement advanced designs of electroacoustic filters and high-performance cavity resonators. The calculations are based on the finite element method in considering the elastic and piezoelectric properties of the materials. Interdigital transducers are employed to measure the transmission spectra. The geometrical parameters enabling the appearance of confined cavity modes within the band gap and the efficiency of the electric excitation were investigated.


title: Relaxed micromorphic model of transient wave propagation in anisotropic band-gap metastructures

abstract: In this paper, we show that the transient waveforms arising from several localised pulses in a micro-structured material can be reproduced by a corresponding generalised continuum of the relaxed micromorphic type. Specifically, we compare the dynamic response of a bounded micro-structured material to that of bounded continua with special kinematic properties: (i) the relaxed micromorphic continuum and (ii) an equivalent Cauchy linear elastic continuum. We show that, while the Cauchy theory is able to describe the overall behaviour of the metastructure only at low frequencies, the relaxed micromorphic model goes far beyond by giving a correct description of the pulse propagation in the frequency band-gap and at frequencies intersecting the optical branches. In addition, we observe a computational time reduction associated with the use of the relaxed micromorphic continuum, compared to the sensible computational time needed to perform a transient computation in a micro-structured domain.


title: Long-term optical monitoring of TeV emitting Blazars

abstract: We present 10 years of R-band monitoring data of 31 northern blazars which were either detected at very high energy (VHE) gamma rays or listed as potential VHE gamma-ray emitters. The data comprise 11820 photometric data points in the R-band obtained in 2002-2012. We analyze the light curves by determining their power spectral density (PSD) slopes assuming a power-law dependence with a single slope $\beta$ and a Gaussian probability density function (PDF). We use the multiple fragments variance function (MFVF) combined with a forward-casting approach and likelihood analysis to determine the slopes and perform extensive simulations to estimate the uncertainties of the derived slopes. We also look for periodic variations via Fourier analysis and quantify the false alarm probability through a large number of simulations. Comparing the obtained PSD slopes to values in the literature, we find the slopes in the radio band to be steeper than those in the optical and gamma rays. Our periodicity search yielded one target, Mrk 421, with a significant (p<5%) period. Finding one significant period among 31 targets is consistent with the expected false alarm rate, but the period found in Mrk~421 is very strong and deserves further consideration}.


title: Unitary (g,K) modules of SU(2,1)

abstract: Let G=SU(2,1). In this paper we parametrize irreducible unitary (g,K) modules of G. The parametrization is done in two steps. Firstly, we parametrize irreducible (g,K) modules (Theorem 3). In the second step we find unitary (g,K) modules (Theorem 5). One can compare our results with results of Kraljevic.


title: Stochastic Thermodynamics of oscillators networks

abstract: We apply the stochastic thermodynamics formalism to describe the dynamics of systems of complex Langevin and Fokker-Planck equations. We provide in particular a simple and general recipe to calculate thermodynamical currents, dissipated and propagating heat for networks of nonlinear oscillators. By using the Hodge decomposition of thermodynamical forces and fluxes, we derive a formula for entropy production that generalises the notion of non-potential forces and makes trans- parent the breaking of detailed balance and of time reversal symmetry for states arbitrarily far from equilibrium. Our formalism is then applied to describe the off-equilibrium thermodynamics of a few examples, notably a continuum ferromagnet, a network of classical spin-oscillators and the Frenkel-Kontorova model of nano friction.


title: Schr\"{o}dinger-type equations with sign-changing nonlinearities: a survey

abstract: We are looking for solutions to nonlinear Schr\"odinger-type equations of the form $$ (-\Delta)^{\alpha / 2} u (x) + V(x) u(x) = h (x,u(x)), \quad x \in \mathbb{R}^N, $$ where $V : \mathbb{R}^N \rightarrow \mathbb{R}$ is an external potential (bounded, unbounded or singular) and the nonlinearity $h : \mathbb{R}^N \times \mathbb{R} \rightarrow \mathbb{R}$ may be sign-changing. The aim of this paper is to give an introduction to the problem and the variational approach based on the Nehari manifold technique, and to survey recent results concerning the existence of ground and bound states.


title: Modeling Inner Boundary Values at 18 Solar Radii During Solar Quiet time for Global Three-dimensional Time-Dependent Magnetohydrodynamic Numerical Simulation

abstract: The solar wind speed plays a key role in the transport of CME out of the Sun and ultimately determines the arrival time of CME-driven shocks in the heliosphere. Here, we develop an empirical model of the solar wind parameters at the inner boundary (18 solar radii, Rs) used in our global, 3D MHD model (G3DMHD) or other equivalent ones. The model takes solar magnetic field maps at 2.5 Rs (which is based on the Potential Field Source Surface, PFSS model) and interpolates the solar wind plasma and field out to 18 Rs using the algorithm of Wang and Sheeley [1990a]. A formula V_{18Rs} = V1 + V2 fs^{\alpha} is used to calculate the solar wind speed at 18 Rs, where V1 is in a range of 150-350 km/s, V2 is in the range of 250-500 km/s, and fs is an expansion factor, which was derived from the Wang and Sheeley (WS) algorithm at 2.5 Rs. To estimate the solar wind density and temperature at 18 Rs, we assume an incompressible solar wind and a constant total pressure. The three free parameters are obtained by adjusting simulation results to match in-situ observations (Wind) for more than 54 combination of V1, V2 and {\alpha} during a quiet solar wind interval, CR2082. We found V18Rs = (150 +/- 50) + (500 +/- 100) fs^-0.4 km/s performs reasonably well in predicting solar wind parameters at 1 AU not just for CR 2082 but other quiet solar period. Comparing results from the present study with those from WSA [Arge et al. 2000; 2004] we conclude that i) Results of using V_{18Rs} with the full rotation data (FR) as input to drive G3DMHD model is better than the results of WSA using FR, or daily updated. ii) When using a modified daily updated 4-day-advanced solar wind speed predictions WSA performs slightly better than our G3DMHD. iii) When using V_{18Rs} as input, G3DMHD model performs much better than the WSA formula. We argue the necessity of the extra angular width ({\theta}b) parameter used in WSA.


title: Structure of Gel Phase DPPC Determined by X-ray Diffraction

abstract: High resolution low angle x-ray data are reported for the gel phase of DPPC lipid bilayers, extending the previous q range of 1.0 {\AA}-1 to 1.3 {\AA}-1, and employing a new technique to obtain more accurate intensities and form factors |F(q)| for the highest orders of diffraction. Combined with previous wide angle x-ray and volumetric data, a space filling model is employed to obtain gel phase structure at a mesoscopic level. A new conclusion from this analysis is that the hydrocarbon chains from opposing monolayers are mini-interdigitated; this would help explain the previously well-established result that the opposing monolayers are strongly coupled with respect to their chain tilt directions. Even more detailed structural features are described that have not been obtained from experiment but that could, in principle, be obtained from simulations that would first be validated by agreement with the wide angle and the new low angle |F(q)| x-ray data.


title: Magnetoresistance in Co-hBN-NiFe tunnel junctions enhanced by resonant tunneling through single defects in ultrathin hBN barriers

abstract: Hexagonal boron nitride (hBN) is a prototypical high-quality two-dimensional insulator and an ideal material to study tunneling phenomena, as it can be easily integrated in vertical van der Waals devices. For spintronic devices, its potential has been demonstrated both for efficient spin injection in lateral spin valves and as a barrier in magnetic tunnel junctions (MTJs). Here we reveal the effect of point defects inevitably present in mechanically exfoliated hBN on the tunnel magnetoresistance of Co-hBN-NiFe MTJs. We observe a clear enhancement of both the conductance and magnetoresistance of the junction at well-defined bias voltages, indicating resonant tunneling through magnetic (spin-polarized) defect states. The spin polarization of the defect states is attributed to exchange coupling of a paramagnetic impurity in the few-atomic-layer thick hBN to the ferromagnetic electrodes. This is confirmed by excellent agreement with theoretical modelling. Our findings should be taken into account in analyzing tunneling processes in hBN-based magnetic devices. More generally, our study shows the potential of using atomically thin hBN barriers with defects to engineer the magnetoresistance of MTJs and to achieve spin filtering, opening the door towards exploiting the spin degree of freedom in current studies of point defects as quantum emitters.


title: A Learning-based Power Management for Networked Microgrids Under Incomplete Information

abstract: This paper presents an approximate Reinforcement Learning (RL) methodology for bi-level power management of networked Microgrids (MG) in electric distribution systems. In practice, the cooperative agent can have limited or no knowledge of the MG asset behavior and detailed models behind the Point of Common Coupling (PCC). This makes the distribution systems unobservable and impedes conventional optimization solutions for the constrained MG power management problem. To tackle this challenge, we have proposed a bi-level RL framework in a price-based environment. At the higher level, a cooperative agent performs function approximation to predict the behavior of entities under incomplete information of MG parametric models; while at the lower level, each MG provides power-flow-constrained optimal response to price signals. The function approximation scheme is then used within an adaptive RL framework to optimize the price signal as the system load and solar generation change over time. Numerical experiments have verified that, compared to previous works in the literature, the proposed privacy-preserving learning model has better adaptability and enhanced computational speed.


title: $Q$-Series and Quantum Spin Networks

abstract: The tail of a quantum spin network in the two-sphere is a $q$-series associated to the network. We study the existence of the head and tail functions of quantum spin networks colored by $2n$. We compute the $q$-series for an infinite family of quantum spin networks and give the relation between the tail of these networks and the tail of the colored Jones polynomial. Finally, we show that the family of quantum spin networks under study satisfies a natural product structure, making these networks satisfy a natural product structure.


title: A Coordinated Volt-Var Control Scheme for Distribution Systems with High DER Penetration

abstract: In this paper, a new Volt/Var Control (VVC) scheme is proposed to facilitate the coordination between the conventional VVC devices and the new smart PV inverters to provide an effective voltage control on a system with high PV penetration. The proposed scheme decomposes the problem into two levels. The first level uses Load Tap Changer (LTC) and Voltage Regulators (VRs) to adjust the voltage level on the circuit to keep the voltages along the circuit within the desired range. The second level determines Var support needed from smart inverters to smooth the fast voltage variations while providing effective power factor correction to keep the power losses at minimum. The case study shows that the proposed VVC method is very effective in maintaining acceptable voltages on the system under various operating conditions while meeting the operational constrains. The results also show the computational efficiency of the method.


title: Simulation of elliptic and hypo-elliptic conditional diffusions

abstract: Suppose $X$ is a multidimensional diffusion process. Assume that at time zero the state of $X$ is fully observed, but at time $T>0$ only linear combinations of its components are observed. That is, one only observes the vector $L X_T$ for a given matrix $L$. In this paper we show how samples from the conditioned process can be generated. The main contribution of this paper is to prove that guided proposals, introduced in Schauer et al. (2017), can be used in a unified way for both uniformly and hypo-elliptic diffusions, also when $L$ is not the identity matrix. This is illustrated by excellent performance in two challenging cases: a partially observed twice integrated diffusion with multiple wells and the partially observed FitzHugh-Nagumo model.


title: On the spectral radius of compact operator cocycles

abstract: We extend the notions of joint and generalized spectral radii to cocycles acting on Banach spaces and obtain a version of Berger-Wang's formula when restricted to the space of cocycles taking values in the space of compact operators. Moreover, we observe that the previous quantities depends continuously on the underlying cocycle.


title: Algorithms for Destructive Shift Bribery

abstract: We study the complexity of Destructive Shift Bribery. In this problem, we are given an election with a set of candidates and a set of voters (each ranking the candidates from the best to the worst), a despised candidate $d$, a budget $B$, and prices for shifting $d$ back in the voters' rankings. The goal is to ensure that $d$ is not a winner of the election. We show that this problem is polynomial-time solvable for scoring protocols (encoded in unary), the Bucklin and Simplified Bucklin rules, and the Maximin rule, but is NP-hard for the Copeland rule. This stands in contrast to the results for the constructive setting (known from the literature), for which the problem is polynomial-time solvable for $k$-Approval family of rules, but is NP-hard for the Borda, Copeland, and Maximin rules. We complement the analysis of the Copeland rule showing W-hardness for the parameterization by the budget value, and by the number of affected voters. We prove that the problem is W-hard when parameterized by the number of voters even for unit prices. From the positive perspective we provide an efficient algorithm for solving the problem parameterized by the combined parameter the number of candidates and the maximum bribery price (alternatively the number of different bribery prices).


title: $\Delta \left( 27\right)$ flavor singlet-triplet Higgs model for fermion masses and mixings

abstract: We propose a multiscalar singlet extension of the singlet-triplet Higgs model capable of explaining the SM fermion mass spectrum and mixing parameters. Our model is based on the $\Delta \left( 27\right) $ family symmetry, supplemented by cyclic symmetries, which are spontaneously broken thus yielding the observed hierarchy of the SM charged fermion masses and quark mixing angles. The masses of the light active neutrinos are produced by type-II seesaw mechanism mediated by the neutral component of the $SU(2)_{L}$ scalar triplet. The model symmetries lead to the extended Gatto-Sartori-Tonin relations between the quark masses and mixing angles.


title: Predicting Factuality of Reporting and Bias of News Media Sources

abstract: We present a study on predicting the factuality of reporting and bias of news media. While previous work has focused on studying the veracity of claims or documents, here we are interested in characterizing entire news media. These are under-studied but arguably important research problems, both in their own right and as a prior for fact-checking systems. We experiment with a large list of news websites and with a rich set of features derived from (i) a sample of articles from the target news medium, (ii) its Wikipedia page, (iii) its Twitter account, (iv) the structure of its URL, and (v) information about the Web traffic it attracts. The experimental results show sizable performance gains over the baselines, and confirm the importance of each feature type.


title: Collective Effects and Performance of Algorithmic Electric Vehicle Charging Strategies

abstract: We combine the power flow model with the proportionally fair optimization criterion to study the control of congestion within a distribution electric grid network. The form of the mathematical optimization problem is a convex second order cone that can be solved by modern non-linear interior point methods and constitutes the core of a dynamic simulation of electric vehicles (EV) joining and leaving the charging network. The preferences of EV drivers, represented by simple algorithmic strategies, are conveyed to the optimizing component by real-time adjustments to user-specific weighting parameters that are then directly incorporated into the objective function. The algorithmic strategies utilize a small number of parameters that characterize the user's budgets, expectations on the availability of vehicles and the charging process. We investigate the collective behaviour emerging from individual strategies and evaluate their performance by means of computer simulation.


title: Mass Functions, Luminosity Functions, and Completeness Measurements from Clustering Redshifts

abstract: This paper presents stellar mass functions and i-band luminosity functions for Sloan Digital Sky Survey (SDSS) galaxies at $i < 21$ using clustering redshifts, from which we also compute targeting completeness measurements for the Baryon Oscillation Spectroscopic Survey (BOSS). Clustering redshifts is a method of obtaining the redshift distribution of a sample of galaxies with only photometric information by measuring the angular crosscorrelation with a spectroscopic sample in different redshift bins. We construct a spectroscopic sample containing data from the BOSS + eBOSS surveys, allowing us to recover redshift distributions from photometric data out to $z\simeq 2.5$. We produce k-corrected i-band luminosity functions and stellar mass functions by applying clustering redshifts to SDSS DR8 galaxies in small bins of colour and magnitude. There is little evolution in the mass function between $0.2 < z < 0.8$, implying the most massive galaxies form most of their mass before $z = 0.8$. These mass functions are used to produce stellar mass completeness estimates for the Baryon Oscillation Spectroscopic Survey (BOSS), giving a stellar mass completeness of $80\%$ above $M_{\star} > 10^{11.4}$ between $0.2 < z < 0.7$, with completeness falling significantly at redshifts higher than 0.7, and at lower masses. Large photometric datasets will be available in the near future (DECaLS, DES, Euclid), so this, and similar techniques will become increasingly useful in order to fully utilise this data.


title: Three remarkable properties of the Normal distribution

abstract: In this paper, we present three remarkable properties of the normal distribution: first that if two independent variables's sum is normally distributed, then each random variable follows a normal distribution (which is referred to as the Levy Cramer theorem), second a variation of the Levy Cramer theorem that states that if two independent symmetric random variables with finite variance have their sum and their difference independent, then each random variable follows a standard normal distribution, and third that the normal distribution is characterized by the fact that it is the only distribution for which the sample mean and variance are independent (which is a central property for deriving the Student distribution and referred as the Geary theorem). The novelty of this paper is to provide new, quicker or self contained proofs of theses theorems.


title: Quantum circuits of c -- Z and SWAP gates optimization and entanglement

abstract: We have studied the algebraic structure underlying the quantum circuits composed by $\cZ$ and $\Swap$ gates. Our results are applied to optimize the circuits and to understand the emergence of entanglement when a circuit acts on a fully factorized state.


title: Description of Glass Transition kinetics in 3D XY-model in terms of Gauge Field Theory

abstract: We consider a gauge theory of the glass transition in the frustrated XY model being simplest model containing topologically nontrivial excitations. We describe the transition kinetics and find that the three-dimensional system exhibits the Vogel-Fulcher-Tamman criticality heralding its freezing into a spin glass. We analytically show that the system demonstrates all glass transition properties, like the logarithmic relaxation, and corresponding behavior of linear and non-linear susceptibility. The mode-coupling theory equation in the Zwanziger-Mori representation also is derived in framework of our approach. Our findings provide insights into the topological origin of glass formation, that allows to make progress in understanding glass-transition processes in more intricate systems.


title: SAVOIAS: A Diverse, Multi-Category Visual Complexity Dataset

abstract: Visual complexity identifies the level of intricacy and details in an image or the level of difficulty to describe the image. It is an important concept in a variety of areas such as cognitive psychology, computer vision and visualization, and advertisement. Yet, efforts to create large, downloadable image datasets with diverse content and unbiased groundtruthing are lacking. In this work, we introduce Savoias, a visual complexity dataset that compromises of more than 1,400 images from seven image categories relevant to the above research areas, namely Scenes, Advertisements, Visualization and infographics, Objects, Interior design, Art, and Suprematism. The images in each category portray diverse characteristics including various low-level and high-level features, objects, backgrounds, textures and patterns, text, and graphics. The ground truth for Savoias is obtained by crowdsourcing more than 37,000 pairwise comparisons of images using the forced-choice methodology and with more than 1,600 contributors. The resulting relative scores are then converted to absolute visual complexity scores using the Bradley-Terry method and matrix completion. When applying five state-of-the-art algorithms to analyze the visual complexity of the images in the Savoias dataset, we found that the scores obtained from these baseline tools only correlate well with crowdsourced labels for abstract patterns in the Suprematism category (Pearson correlation r=0.84). For the other categories, in particular, the objects and advertisement categories, low correlation coefficients were revealed (r=0.3 and 0.56, respectively). These findings suggest that (1) state-of-the-art approaches are mostly insufficient and (2) Savoias enables category-specific method development, which is likely to improve the impact of visual complexity analysis on specific application areas, including computer vision.


title: Measurement of the top quark mass in the $t\bar{t}\to$ lepton+jets channel from $\sqrt{s}=8$ TeV ATLAS data and combination with previous results

abstract: The top quark mass is measured using a template method in the $t\bar{t}\to$ lepton+jets channel (lepton is $e$ or $\mu$) using ATLAS data recorded in 2012 at the LHC. The data were taken at a proton-proton centre-of-mass energy of $\sqrt{s}=8$ TeV and correspond to an integrated luminosity of 20.2 fb$^{-1}$. The $t\bar{t}\to$ lepton+jets channel is characterized by the presence of a charged lepton, a neutrino and four jets, two of which originate from bottom quarks ($b$). Exploiting a three-dimensional template technique, the top quark mass is determined together with a global jet energy scale factor and a relative $b$-to-light-jet energy scale factor. The mass of the top quark is measured to be $m_{top}= 172.08 \pm 0.39 (stat) \pm 0.82 (syst)$ GeV. A combination with previous ATLAS $m_{top}$ measurements gives $m_{top}= 172.69 \pm 0.25 (stat) \pm 0.41 (syst)$ GeV.


title: Optomechanics of chiral dielectric metasurfaces

abstract: Electromagnetic fields coupled with mechanical degrees of freedom have recently shown exceptional and innovative applications, ultimately leading to mesoscopic optomechanical devices operating in the quantum regime of motion. Simultaneously, micromechanical elements have provided new ways to enhance and manipulate the optical properties of passive photonic elements. Following this concept, in this article we show how combining a chiral metasurface with a GaAs suspended micromembrane can offer new scenarios for controlling the polarization state of near-infrared light beams. Starting from the uncommon properties of chiral metasurface to statically realize target polarization states and circular and linear dichroism, we report mechanically induced, ~300 kHz polarization modulation, which favorably compares, in terms of speed, with liquid-crystals commercial devices. Moreover, we demonstrate how the mechanical resonance can be non-trivially affected by the input light polarization (and chiral state) via a thermoelastic effect triggered by intracavity photons. This work inaugurates the field of Polarization Optomechanics, which could pave the way to fast polarimetric devices, polarization modulators and dynamically tunable chiral state generators and detectors, as well as giving access to new form of polarization nonlinearities and control.


title: On typical representations for depth-zero components of split classical groups

abstract: Let ${\bf G}$ be a split classical group over a non-Archimedean local field $F$ with the cardinality of the residue field $q_F>5$. Let $M$ be the group of $F$-points of a Levi factor of a proper $F$-parabolic subgroup of ${\bf G}$. Let $[M, \sigma_M]_M$ be an inertial class such that $\sigma_M$ contains a depth-zero Moy--Prasad type of the form $(K_M, \tau_M)$, where $K_M$ is a hyperspecial maximal compact subgroup of $M$. Let $K$ be a hyperspecial maximal compact subgroup of ${\bf G}(F)$ such that $K$ contains $K_M$. In this article, we classify $\mathfrak{s}$-typical representations of $K$. In particular, we show that the $\mathfrak{s}$-typical representations of $K$ are precisely the irreducible subrepresentations of $\ind_J^K\lambda$, where $(J, \lambda)$ is a level-zero $G$-cover of $(K\cap M, \tau_M)$.


title: A rigidity result for the Holm-Staley b-family of equations with application to the asymptotic stability of the Degasperis-Procesi peakon

abstract: We prove that the peakons are asymptotically H 1-stable, under the flow of the Degasperis-Procesi equation, in the class of functions with a momentum density that belongs to M + (R). The key argument is a rigidity result for uniformly in time exponentially decaying global solutions that is shared by the Holm-Staley b-family of equations for b $\ge$ 1. This extends previous results obtained for the Camassa-Holm equation (b = 2).


title: Efficient Fastest-Path Computations in Road Maps

abstract: In the age of real-time online traffic information and GPS-enabled devices, fastest-path computations between two points in a road network modeled as a directed graph, where each directed edge is weighted by a "travel time" value, are becoming a standard feature of many navigation-related applications. To support this, very efficient computation of these paths in very large road networks is critical. Fastest paths may be computed as minimal-cost paths in a weighted directed graph, but traditional minimal-cost path algorithms based on variants of the classic Dijkstra algorithm do not scale well, as in the worst case they may traverse the entire graph. A common improvement, which can dramatically reduce the number of traversed graph vertices, is the A* algorithm, which requires a good heuristic lower bound on the minimal cost. We introduce a simple, but very effective, heuristic function based on a small number of values assigned to each graph vertex. The values are based on graph separators and computed efficiently in a preprocessing stage. We present experimental results demonstrating that our heuristic provides estimates of the minimal cost which are superior to those of other heuristics. Our experiments show that when used in the A* algorithm, this heuristic can reduce the number of vertices traversed by an order of magnitude compared to other heuristics.


title: Acoustic flat lensing using an indefinite medium

abstract: Acoustic flat lensing is achieved here by tuning a phononic array to have indefinite medium behaviour in a narrow frequency spectral region along the acoustic branch. This is confirmed by the occurrence of a flat band along an unusual path in the Brillouin zone and by interpreting the intersection point of isofrequency contours on the corresponding isofrequency surface; coherent directive beams are formed whose reflection from the array surfaces create lensing. Theoretical predictions are corroborated by time-domain experiments, airborne acoustic waves generated by a source with a frequency centered about $10.6$ kHz, placed at three different distances from one side of a finite phononic crystal slab, constructed from polymeric spheres, yield distinctive focal spots on the other side. These experiments evaluate the pressure field using optical feedback interferometry and demonstrate precise control of the three-dimensional wave trajectory through a sonic crystal.


title: A Bayesian model for sparse graphs with flexible degree distribution and overlapping community structure

abstract: We consider a non-projective class of inhomogeneous random graph models with interpretable parameters and a number of interesting asymptotic properties. Using the results of Bollob\'as et al. [2007], we show that i) the class of models is sparse and ii) depending on the choice of the parameters, the model is either scale-free, with power-law exponent greater than 2, or with an asymptotic degree distribution which is power-law with exponential cut-off. We propose an extension of the model that can accommodate an overlapping community structure. Scalable posterior inference can be performed due to the specific choice of the link probability. We present experiments on five different real-world networks with up to 100,000 nodes and edges, showing that the model can provide a good fit to the degree distribution and recovers well the latent community structure.


title: Stochastic Methods for the Neutron Transport Equation I: Linear Semigroup asymptotics

abstract: The Neutron Transport Equation (NTE) describes the flux of neutrons through an inhomogeneous fissile medium. In this paper, we reconnect the NTE to the physical model of the spatial Markov branching process which describes the process of nuclear fission, transport, scattering, and absorption. By reformulating the NTE in its mild form and identifying its solution as an expectation semigroup, we use modern techniques to develop a Perron-Frobenius (PF) type decomposition, showing that growth is dominated by a leading eigenfunction and its associated left and right eigenfunctions. In the spirit of results for spatial branching and fragmentation processes, we use our PF decomposition to show the existence of an intrinsic martingale and associated spine decomposition. Moreover, we show how criticality in the PF decomposition dictates the convergence of the intrinsic martingale. The mathematical difficulties in this context come about through unusual piecewise linear motion of particles coupled with an infinite type-space which is taken as neutron velocity. The fundamental nature of our PF decomposition also plays out in accompanying work [20, 9].


title: RadioPropa - A Modular Raytracer for In-Matter Radio Propagation

abstract: Experiments for radio detection of UHE particles such as e.g. ARA/ARIANNA or NuMoon require detailed understanding of the propagation of radio waves in the surrounding matter. The index of refraction in e.g. polar ice or lunar rock may have a complex spatial structure that makes detailed simulations of the radio propagation necessary to design the respective experiments and analyze their data. Here, we present RadioPropa as a new modular ray tracing code that solves the eikonal equation with a Runge-Kutta method in arbitrary refractivity fields. RadioPropa is based on the cosmic ray propagation code CRPropa, which has been forked to allow efficient incorporation of the required data structures for ray tracing while retaining its modular design. This allows for the setup of versatile simulation geometries as well as the easy inclusion of additional physical effects such as e.g. partial reflection on boundary layers in the simulations. We discuss the principal design of the code as well as its performance in example applications.


title: Constraining neutrino mass with tomographic weak lensing peak counts

abstract: Massive cosmic neutrinos change the structure formation history by suppressing perturbations on small scales. Weak lensing data from galaxy surveys probe the structure evolution and thereby can be used to constrain the total mass of the three active neutrinos. However, much of the information is at small scales where the dynamics are nonlinear. Traditional approaches with second order statistics thus fail to fully extract the information in the lensing field. In this paper, we study constraints on the neutrino mass sum using lensing peak counts, a statistic which captures non-Gaussian information beyond the second order. We use the ray-traced weak lensing mocks from the Cosmological Massive Neutrino Simulations (MassiveNuS), and apply LSST-like noise. We discuss the effects of redshift tomography, multipole cutoff $\ell_{\rm max}$ for the power spectrum, smoothing scale for the peak counts, and constraints from peaks of different heights. We find that combining peak counts with the traditional lensing power spectrum can improve the constraint on neutrino mass sum, $\Omega_m$, and $A_s$ by 39%, 32%, and 60%, respectively, over that from the power spectrum alone.


title: Approximation and Loewner Theory of holomorphic covering mappings

abstract: We give conditions in order to approximate locally uniformly holomorphic covering mappings of the unit ball of $\mathbb{C}^n$ with respect to an arbitrary norm, with entire holomorphic covering mappings. The results rely on a generalization of the Loewner theory for covering mappings which we develop in the paper.


title: A Direct Proof of the Reflection Principle for Brownian Motion

abstract: We present a self-contained proof of the reflection principle for Brownian Motion.


title: Determining r-Robustness of Digraphs Using Mixed Integer Linear Programming

abstract: Convergence guarantees of many resilient consensus algorithms are based on the graph theoretic properties of $r$- and $(r,s)$-robustness. These algorithms guarantee consensus of normally behaving agents in the presence of a bounded number of arbitrarily misbehaving agents if the values of the integers $r$ and $s$ are sufficiently high. However, determining the largest integer $r$ for which an arbitrary digraph is $r$-robust is highly nontrivial. This paper introduces a novel method for calculating this value using mixed integer linear programming. The method only requires knowledge of the graph Laplacian matrix, and can be formulated with affine objective and constraints, except for the integer constraint. Integer programming methods such as branch-and-bound can allow both lower and upper bounds on $r$ to be iteratively tightened. Simulations suggest the proposed method demonstrates greater efficiency than prior algorithms.


title: Weighted dynamic finger in binary search trees

abstract: It is shown that the online binary search tree data structure GreedyASS performs asymptotically as well on a sufficiently long sequence of searches as any static binary search tree where each search begins from the previous search (rather than the root). This bound is known to be equivalent to assigning each item $i$ in the search tree a positive weight $w_i$ and bounding the search cost of an item in the search sequence $s_1,\ldots,s_m$ by $$O\left(1+ \log \frac{\displaystyle \sum_{\min(s_{i-1},s_i) \leq x \leq \max(s_{i-1},s_i)}w_x}{\displaystyle \min(w_{s_i},w_{s_{i-1}})} \right)$$ amortized. This result is the strongest finger-type bound to be proven for binary search trees. By setting the weights to be equal, one observes that our bound implies the dynamic finger bound. Compared to the previous proof of the dynamic finger bound for Splay trees, our result is significantly shorter, stronger, simpler, and has reasonable constants.


title: Towards Uniform Online Spherical Tessellations

abstract: The problem of uniformly placing N points onto a sphere finds applications in many areas. For example, points on the sphere correspond to unit quaternions as well as to the group of rotations SO(3) and the online version of generating uniform rotations (known as "incremental generation") plays a crucial role in a large number of engineering applications ranging from robotics and aeronautics to computer graphics. An online version of this problem was recently studied with respect to the gap ratio as a measure of uniformity. The first online algorithm of Chen et al. was upper-bounded by 5.99 and later improved to 3.69, which is achieved by considering a circumscribed dodecahedron followed by a recursive decomposition of each face. In this paper we provide a more efficient tessellation technique based on the regular icosahedron, which improves the upper-bound for the online version of this problem, decreasing it to approximately 2.84. Moreover, we show that the lower bound for the gap ratio of placing at least three points is the golden ratio, approx. 1.618, and for at least four points is no less than 1.726.


title: Electronic quantum effect on hydrogen bond geometry in water dimer

abstract: During compression of a water dimer calculated with high-precision first-principles methods, the trends of H-bond and O-H bond lengths show quantum effect of the electronic structure. We found that the H-bond length keeps decreasing, while the O-H bond length increases up to the stable point and decreases beyond it when the water dimer is further compressed. The remarkable properties are different from those observed in most previous researches which can be understood and extrapolated through classical simulation. The observations can be explained by the decrease in orbital overlap and change in the exchange repulsion interaction between water monomers. The dominant interaction between water monomers changes from electrostatic interaction to exchange repulsion at the turning point of the O-H bond length when the O...O distance is decreased. These findings highlight the quantum effect on the hydrogen bond in extreme conditions and play an important role in the recognition of the hydrogen bond structure and mechanism.


title: The s process in AGB stars as constrained by a large sample of Barium stars

abstract: Context. Barium (Ba) stars are dwarf and giant stars enriched in elements heavier than iron produced by the slow neutron-capture process (s process). They belong to binary systems where the primary star evolved through the asymptotic giant branch (AGB) phase,during which it produced the s-process elements and transferred them onto the secondary, now observed as a Ba star. Aims. We compare the largest homogeneous set of Ba giant star observations of the s-process elements Y, Zr, La, Ce, and Nd with AGB nucleosynthesis models to reach a better understanding of the s process in AGB stars. Methods. By considering the light-s (ls: Y and Zr) heavy-s (hs: La, Ce, and Nd) and elements individually, we computed for the first time quantitative error bars for the different hs-element/ls-element abundance ratios, and for each of the sample stars. We compared these ratios to low-mass AGB nucleosynthesis models. We excluded La from our analysis because the strong La lines in some of the sample stars cause an overestimation and unreliable abundance determination, as compared to the other observed hs-type elements. Results. All the computed hs-type to ls-type element ratios show a clear trend of increasing with decreasing metallicity with a small spread (less than a factor of 3). This trend is predicted by low-mass AGB models where 13C is the main neutron source. The comparison with rotating AGB models indicates the need for the presence of an angular momentum transport mechanism that should not transport chemical species, but significantly reduce the rotational speed of the core in the advanced stellar evolutionary stages. This is an independent confirmation of asteroseismology observations of the slow down of core rotation in giant stars, and of rotational velocities of white dwarfs lower than predicted by models without an extra angular momentum transport mechanism.


title: Eliminating the effect of acoustic noise on cantilever spring constant calibration

abstract: A common use for atomic force microscopy is to quantify local forces through tip-sample interactions between the probe tip and a sample surface. The accuracy of these measurements depends on the accuracy to which the cantilever spring constant is known. Recent work has demonstrated that the measured spring constant of a cantilever can vary up to a factor of two, even for the exact same cantilever measured by different users on different microscopes. Here we demonstrate that a standard method for calibrating the spring constant (using the oscillations due to thermal energy) is susceptible to ambient noise, which can alter the result significantly. We demonstrate a new step-by-step method to measure the spring constant by actively driving the cantilever to measure the resonance frequency and quality factor, giving results that are unaffected by acoustic noise. Our method can be performed rapidly on any atomic force microscope without any expensive additional hardware.


title: A note on the quantum query complexity of permutation symmetric functions

abstract: It is known since the work of [AA14] that for any permutation symmetric function $f$, the quantum query complexity is at most polynomially smaller than the classical randomized query complexity, more precisely that $R(f) = \widetilde{O}\left(Q^7(f)\right)$. In this paper, we improve this result and show that $R(f) = {O}\left(Q^3(f)\right)$ for a more general class of symmetric functions. Our proof is constructive and relies largely on the quantum hardness of distinguishing a random permutation from a random function with small range from Zhandry [Zha15].


title: FixMiner: Mining Relevant Fix Patterns for Automated Program Repair

abstract: Patching is a common activity in software development. It is generally performed on a source code base to address bugs or add new functionalities. In this context, given the recurrence of bugs across projects, the associated similar patches can be leveraged to extract generic fix actions. While the literature includes various approaches leveraging similarity among patches to guide program repair, these approaches often do not yield fix patterns that are tractable and reusable as actionable input to APR systems. In this paper, we propose a systematic and automated approach to mining relevant and actionable fix patterns based on an iterative clustering strategy applied to atomic changes within patches. The goal of FixMiner is thus to infer separate and reusable fix patterns that can be leveraged in other patch generation systems. Our technique, FixMiner, leverages Rich Edit Script which is a specialized tree structure of the edit scripts that captures the AST-level context of the code changes. FixMiner uses different tree representations of Rich Edit Scripts for each round of clustering to identify similar changes. These are abstract syntax trees, edit actions trees, and code context trees. We have evaluated FixMiner on thousands of software patches collected from open source projects. Preliminary results show that we are able to mine accurate patterns, efficiently exploiting change information in Rich Edit Scripts. We further integrated the mined patterns to an automated program repair prototype, PARFixMiner, with which we are able to correctly fix 26 bugs of the Defects4J benchmark. Beyond this quantitative performance, we show that the mined fix patterns are sufficiently relevant to produce patches with a high probability of correctness: 81% of PARFixMiner's generated plausible patches are correct.


title: Detection of planetary signals by reflected light of the host star using the autocorrelation of spectra

abstract: We consider an alternative to the cross-correlation function (CCF), that uses the autocor- relation function (ACF), to measure in spectra the reflected light of the stars by their planetary companion. The objective of this work is to assess and analyse the efficiency of the ACF in detecting planetary signals by a detection of reflected light. To do so, we first compare the ACF and the CCF using artificial spectra containing a planetary signal. We then use the ACF to analyse spectra of the 51 Peg + 51 Peg b system and compare our results with those obtained by Martins et al. (2015) who previously analysed the same system using the CCF. The functionalities of the ACF and the way it is implemented are similar to that of the CCF. One of the main advantages of the ACF is the fact that, unlike the CCF, no weighted binary mask is required. This makes the ACF easier to use. The results related to simulated spectra showed that the ACF allowed us to decrease the boundary condition threshold for the use of spectra compared to the CCF so that more data could be used to recover the planetary signal. For the 51 Peg b planet, we achieved a detection significance of 5.52 {\sigma}noise with the ACF compared to 3.70 {\sigma}noise with the CCF (Martins et al. 2015). We conclude that the ACF has the potential to become a prominent technique in detecting exoplanets considering its efficiency, ease of use and rapid execution time.


title: Length spectra of flat metrics coming from q-differentials

abstract: When geometric structures on surfaces are determined by the lengths of curves, it is natural to ask: which curves' lengths do we really need to know? It is a result of Duchin--Leininger--Rafi that any flat metric induced by a unit-norm quadratic differential is determined by its marked simple length spectrum. We generalize the notion of simple curves to that of q-simple curves, for any positive integer q, and show that the lengths of q-simple curves suffice to determine a non-positively curved Euclidean cone metric induced by a q-differential.


title: VStore: A Data Store for Analytics on Large Videos

abstract: We present VStore, a data store for supporting fast, resource-efficient analytics over large archival videos. VStore manages video ingestion, storage, retrieval, and consumption. It controls video formats along the video data path. It is challenged by i) the huge combinatorial space of video format knobs; ii) the complex impacts of these knobs and their high profiling cost; iii) optimizing for multiple resource types. It explores an idea called backward derivation of configuration: in the opposite direction along the video data path, VStore passes the video quantity and quality expected by analytics backward to retrieval, to storage, and to ingestion. In this process, VStore derives an optimal set of video formats, optimizing for different resources in a progressive manner. VStore automatically derives large, complex configurations consisting of more than one hundred knobs over tens of video formats. In response to queries, VStore selects video formats catering to the executed operators and the target accuracy. It streams video data from disks through decoder to operators. It runs queries as fast as 362x of video realtime.


title: Phase Separation Dynamics Induced by an Interaction Quench of a Correlated Fermi-Fermi Mixture in a Double Well

abstract: We explore the interspecies interaction quench dynamics of ultracold spin-polarized few-body mass balanced Fermi-Fermi mixtures confined in a double-well with an emphasis on the beyond Hartree-Fock correlation effects. It is shown that the ground state of particle imbalanced mixtures exhibits a symmetry breaking of the single-particle density for strong interactions in the Hartree-Fock limit, which is altered within the many-body approach. Quenching the interspecies repulsion towards the strongly interacting regime the two species phase separate within the Hartree-Fock approximation while remaining miscible in the many-body treatment. Despite their miscible character on the one-body level the two species are found to be strongly correlated and exhibit a phase separation on the two-body level that suggests the anti-ferromagnetic like behavior of the few-body mixture. For particle balanced mixtures we show that an intrawell fragmentation (filamentation) of the density occurs both for the ground state as well as upon quenching from weak to strong interactions, a result that is exclusively caused by the presence of strong correlations. Inspecting the two-body correlations a phase separation of the two species is unveiled being a precursor towards an anti-ferromagnetic state. Finally, we simulate in-situ single-shot measurements and showcase how our findings can be retrieved by averaging over a sample of single-shot images.


title: Accretion Disk MHD Winds and Blazar Classification

abstract: The Fermi Gamma-Ray Space Telescope observations of blazars show a strong correlation between the spectral index of their gamma-ray spectra and their synchrotron peak frequency $\nu_{\rm{pk}}^{\rm{syn}}$; additionally, the rate of Compton Dominance of these sources also seems to be a function of $\nu_{\rm{pk}}^{\rm{syn}}$. In this work, we adopt the assumption that the nonthermal emission of blazars is primarily due to radiation by a population of Fermi-accelerated electrons in a relativistic outflow (jet) along the symmetry axis of the blazar's accretion disk. Furthermore, we assume that the Compton component is related to an external photon field of photons, which are scattered from particles of the magnetohydrodynamic (MHD) wind emanating from the accretion disk. Our results reproduce well the aforementioned basic observational trends of blazar classification by varying just one parameter, namely the mass accretion rate onto the central black hole.


title: Parametric Feedback Cooling of Rigid Body Nanodumbbells in Levitated Optomechanics

abstract: We theoretically investigate the rigid body dynamics of an optically levitated nanodumbbell under parametric feedback cooling and provide a simplified model for describing the motion. Differing from previous studies, the spin of the nanoparticle about its symmetry axis is considered non-negligible. Simulations reveal that standard parametric feedback cooling can extract energy from two of the five rotational degrees of freedom when the nanoparticle is levitated using a linearly polarized laser beam. The dynamics after feedback cooling are characterized by a normal mode describing precession about the laser polarization axis together with spin about the nanoparticle's symmetry axis. Cooling the remaining mode requires an asymmetry in the two librational frequencies associated with motion about the polarization axis as well as information about the two frequencies of rotation about the polarization axis. Introducing an asymmetric potential allows full cooling of the librational coordinates if the frequencies of both are used in the feedback modulation and is an avenue for entering the librational quantum regime. The asymmetry in the potential needs to be large enough for practical cooling times as the cooling rate of the system depends non-linearly on the degree of asymmetry, a condition that is easily achieved experimentally.

categories: quant-ph cond-mat.mes-hall physics.app-ph physics.optics


title: The White Dwarf Luminosity Functions from the Pan-STARRS 1 3{\pi} Steradian Survey

abstract: A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3{\pi} Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3{\pi} sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100.


title: Simplest phonons and pseudo-phonons in field theory

abstract: We study the emergence of Nambu-Goldstone modes due to broken translation symmetry in field theory. Purely spontaneous breaking yields a massless phonon which develops a mass upon introducing a perturbative explicit breaking. The pseudo-phonon mass agrees with Gell Mann-Oakes-Renner relations. We analyze the simplest possible theories featuring gradient Mexican hats and describing space-dependent order parameters. We comment on homogeneous translation breaking and the connections with holographic Q-lattices.


title: Recent Developments in Gluon Fusion Higgs Calculations

abstract: During recent years perturbative fixed order and resummation calculations have decreased uncertainties on predictions for gluon fusion Higgs production cross sections tremendously. Most exciting results have been published just this year. In these proceedings I present an overview of recent and most recent developments of these calculations that allow theory predictions to compete with the experimental precision reached by future collider upgrades.


title: Symmetry-Breaking Interlayer Dzyaloshinskii-Moriya Interactions in Synthetic Antiferromagnets

abstract: The magnetic interfacial Dzyaloshinskii-Moriya interaction (DMI) in multi-layered thin films can lead to exotic chiral spin states, of paramount importance for future spintronic technologies. Interfacial DMI is normally manifested as an intralayer interaction, mediated via a paramagnetic heavy metal in systems lacking inversion symmetry. Here we show how, by designing synthetic antiferromagnets with canted magnetization states, it is also possible to observe interfacial interlayer-DMI at room temperature. The interlayer-DMI breaks the symmetry of the magnetic reversal process via the emergence of noncollinear spin states, which results in chiral exchange-biased hysteresis loops. This work opens up yet unexplored avenues for the development of new chiral spin textures in multi-layered thin film systems.


title: Superselection Sectors of Gravitational Subregions

abstract: Motivated by the problem of defining the entanglement entropy of the graviton, we study the division of the phase space of general relativity across subregions. Our key requirement is demanding that the separation into subregions is imaginary---i.e., that entangling surfaces are not physical. This translates into a certain condition on the symplectic form. We find that gravitational subregions that satisfy this condition are bounded by surfaces of extremal area. We characterise the 'centre variables' of the phase space of the graviton in such subsystems, which can be taken to be the conformal class of the induced metric in the boundary, subject to a constraint involving the traceless part of the extrinsic curvature. We argue that this condition works to discard local deformations of the boundary surface to infinitesimally nearby extremal surfaces, that are otherwise available for generic codimension-2 extremal surfaces of dimension $\geq$ 2.


title: An alternative approach to the Finger of God in large scale structures

abstract: It is generally accepted that linear theory of growth of structure under gravity produces a squashed structure in the two-point correlation function (2PCF) along the line of sight (LoS). On the other hand, the observed radial spread out structure known as Finger of God (FoG) is attributed to non-linear effects, like virial relaxation or random motions in the inner regions of galaxy clusters, or non-linear terms in collapse theory. In this paper we argue that the squashed structure associated with the redshift-space ($s-$) linear theory 2PCF is obtained only when this function is displayed in real-space ($r-$), or when the mapping from $r-$ to $s-$space is trivialized by one of a series of approximations: the "wide separation" sometimes mistaken by the "distant observer", the $\mu(r) \sim \mu(s)$ in the cosine of the position vector with the LoS, or plainly by approximating $s \sim r$. To display the $s-$space 2PCF properly in a grid in $s-$space requires one to solve for the mapping function $\bm s(\bm r)$ and its inverse, which exists in the small gravitational disturbance case. We find this mapping by using plane of the sky projections of the $r-$ and $s-$ 2PCFs in the simplest case of the linear Kaiser spectrum with a conservative power-law $r-$space 2PCF. We show that even in this simple case, a structure quite similar to the FoG is observed in the small scale region, while in the large scale the expected squashed structure is obtained. We show that the observed structure is a function of three parameters: the dimensionless growth rate for visible matter ($\beta$), the power-law exponent ($\gamma$) and the cosmological distortion given by the Alcock-Paczy{\'n}ski parameter ($AP$). We conclude that there appears to be the way of discerning between cosmological and gravitational redshift distortions avoiding known degeneracies.


title: Two-stage stochastic approximation for dynamic rebalancing of shared mobility systems

abstract: Mobility systems featuring shared vehicles are often unable to serve all potential customers, as the distribution of demand does not coincide with the positions of vehicles at any given time. System operators often choose to reposition these shared vehicles (such as bikes, cars, or scooters) actively during the course of the day to improve service rate. They face a complex dynamic optimization problem in which many integer-valued decisions must be made, using real-time state and forecast information, and within the tight computation time constraints inherent to real-time decision-making. We first present a novel nested-flow formulation of the problem, and demonstrate that its linear relaxation is significantly tighter than one from existing literature. We then adapt a two-stage stochastic approximation scheme from the generic SPAR algorithm due to Powell et al., in which rebalancing plans are optimized against a value function representing the expected cost (in terms of fulfilled and unfulfilled customer demand) of the future evolution of the system. The true value function is approximated by a separable function of contributions due to the rebalancing actions carried out at each station and each time step of the planning horizon. The new algorithm requires surprisingly few iterations to yield high-quality solutions, and is suited to real-time use as it can be terminated early if required. We provide insight into this good performance by examining the mathematical properties of our new flow formulation, and perform rigorous tests on standardized benchmark networks to explore the effect of system size. We then use data from Philadelphia's public bike sharing scheme to demonstrate that the approach also yields performance gains for real systems.


title: Random triangular Burnside groups

abstract: We introduce a model for random groups in varieties of $n$-periodic groups as $n$-periodic quotients of triangular random groups. We show that for an explicit $d_{\mathrm{crit}}\in(1/3,1/2)$, for densities $d\in(1/3,d_{\mathrm{crit}})$ and for $n$ large enough, the model produces \emph{infinite} $n$-periodic groups. As an application, we obtain, for every fixed large enough $n$, for every $p\in (1,\infty)$ an infinite $n$-periodic group with fixed points for all isometric actions on $L^p$-spaces. Our main contribution is to show that certain random triangular groups are uniformly acylindrically hyperbolic.


title: A Boundary-Value Problem for the 3-D Fractional Wave Equation with Singularity

abstract: In this paper, a boundary-value problem for the 3-D wave equation with Caputo and Bessel operators is investigated. Sufficient conditions on the initial data are established for the existence of a unique solution to the considered problem.


title: Disambiguating Music Artists at Scale with Audio Metric Learning

abstract: We address the problem of disambiguating large scale catalogs through the definition of an unknown artist clustering task. We explore the use of metric learning techniques to learn artist embeddings directly from audio, and using a dedicated homonym artists dataset, we compare our method with a recent approach that learn similar embeddings using artist classifiers. While both systems have the ability to disambiguate unknown artists relying exclusively on audio, we show that our system is more suitable in the case when enough audio data is available for each artist in the train dataset. We also propose a new negative sampling method for metric learning that takes advantage of side information such as music genre during the learning phase and shows promising results for the artist clustering task.


title: A Neural Transition-based Model for Nested Mention Recognition

abstract: It is common that entity mentions can contain other mentions recursively. This paper introduces a scalable transition-based method to model the nested structure of mentions. We first map a sentence with nested mentions to a designated forest where each mention corresponds to a constituent of the forest. Our shift-reduce based system then learns to construct the forest structure in a bottom-up manner through an action sequence whose maximal length is guaranteed to be three times of the sentence length. Based on Stack-LSTM which is employed to efficiently and effectively represent the states of the system in a continuous space, our system is further incorporated with a character-based component to capture letter-level patterns. Our model achieves the state-of-the-art results on ACE datasets, showing its effectiveness in detecting nested mentions.


title: On tangential transversality

abstract: This is the first of two closely related papers on transversality. Here we introduce the notion of tangential transversality of two closed subsets of a Banach space. It is an intermediate property between transversality and subtransversality. Using it, we obtain a variety of known results and some new ones in a unified way. Our proofs do not use variational principles and we are concentrated mainly on tangential conditions in the primal space.


title: Logarithmic Hyperseries

abstract: We define the field $\mathbb{L}$ of logarithmic hyperseries, construct on $\mathbb{L}$ natural operations of differentiation, integration, and composition, establish the basic properties of these operations, and characterize these operations uniquely by such properties.


abstract: In this paper, we introduce McTorch, a manifold optimization library for deep learning that extends PyTorch. It aims to lower the barrier for users wishing to use manifold constraints in deep learning applications, i.e., when the parameters are constrained to lie on a manifold. Such constraints include the popular orthogonality and rank constraints, and have been recently used in a number of applications in deep learning. McTorch follows PyTorch's architecture and decouples manifold definitions and optimizers, i.e., once a new manifold is added it can be used with any existing optimizer and vice-versa. McTorch is available at https://github.com/mctorch .


title: Spectral and timing analysis of the bursting pulsar GRO J1744-28 with RXTE observations

abstract: We analyzed RXTE/PCA observations of the bursting pulsar GRO J1744-28 during its two outbursts in 1995-1997. We have found a significant transition, i.e., a sharp change, of both spectral and timing properties around a flux of $\rm 10^{-8}\ ergs\ cm^{-2}\ s^{-1}$ (3-30 keV), which corresponds to a luminosity of 2--8$\rm \times 10^{37}\ ergs\ s^{-1}$ for a distance of 4--8 kpc. We define the faint (bright) state when the flux is smaller (larger) than this threshold. In the faint state, the spectral hardness increases significantly with the increasing flux, while remaining almost constant in the bright state. In addition, we find that the pulsed fraction is positively related to both the flux and the energy (<30keV) in both states. Thanks to the very stable pulse profile shape in all energy bands, a hard X-ray lag could be measured. This lag is only significant in the faint state (reaching $\sim$ 20\,ms between 3-4 keV and 16-20 keV), and much smaller in the bright state ($\lesssim$2ms). We speculate that the discovered transition might be caused by changes of the accretion structure on the surface of the neutron star, and the hard X-ray lags are due to a geometry effect.


title: The AMBRE Project: searching for the closest solar siblings

abstract: Finding solar siblings, that is, stars that formed in the same cluster as the Sun, will yield information about the conditions at the Sun's birthplace. We search for solar sibling candidates in AMBRE, the very large spectra database of solar vicinity stars. Since the ages and chemical abundances of solar siblings are very similar to those of the Sun, we carried out a chemistry- and age-based search for solar sibling candidates. We used high-resolution spectra to derive precise stellar parameters and chemical abundances of the stars. We used these spectroscopic parameters together with Gaia DR2 astrometric data to derive stellar isochronal ages. Gaia data were also used to study the kinematics of the sibling candidates. From the about 17000 stars that are characterized within the AMBRE project, we first selected 55 stars whose metallicities are closest to the solar value (-0.1 < [Fe/H] < 0.1 dex). For these stars we derived precise chemical abundances of several iron-peak, alpha- and neutron-capture elements, based on which we selected 12 solar sibling candidates with average abundances and metallicities between -0.03 to 0.03 dex. Our further selection left us with 4 candidates with stellar ages that are compatible with the solar age within observational uncertainties. For the 2 of the hottest candidates, we derived the carbon isotopic ratios, which are compatible with the solar value. HD186302 is the most precisely characterized and probably the most probable candidate of our 4 best candidates. Very precise chemical characterization and age estimation is necessary to identify solar siblings. We propose that in addition to typical chemical tagging, the study of isotopic ratios can give further important information about the relation of sibling candidates with the Sun. Ideally, asteroseismic age determinations of the candidates could solve the problem of imprecise isochronal ages.


title: On strong tangential transversality

abstract: This is the second of two closely related papers on transversality. Here we introduce the notion of strong tangential transversality of two closed subsets of a Banach space which is a natural sufficient condition for tangential transversality. Some properties of uniform tangent sets are obtained. A new sum rule for the Clarke subdifferential is proven.


title: Rubble Pile Asteroids

abstract: The moniker rubble pile is typically applied to all solar system bodies with Diameter between 200m and 10km - where in this size range there is an abundance of evidence that nearly every object is bound primarily by self-gravity with significant void space or bulk porosity between irregularly shaped constituent particles. The understanding of this population is derived from wide-ranging population studies of derived shape and spin, decades of observational studies in numerous wavelengths, evidence left behind from impacts on planets and moons and the in situ study of a few objects via spacecraft flyby or rendezvous. The internal structure, however, which is responsible for the name rubble pile, is never directly observed, but belies a violent history. Many or most of the asteroids on near-Earth orbits, and the ones most accessible for rendezvous and in situ study, are likely byproducts of the continued collisional evolution of the Main Asteroid Belt.