title: A Global Hartman-Grobman Theorem

abstract: We showed that for any bounded neighborhood of a hyperbolic equilibrium point $x_0$, there is a transformation which is locally homeomorphism, such that the system is changed into a linear system in this neighborhood. If the eigenvalues of $Df(x_0)$ are all located in the left-half complex plane, then there is a homeomorphism on the whole region of attraction such that the nonlinear system on the region of attraction is changed into a linear system under such a coordinate change.


title: Traffic Modelling and Prediction via Symbolic Regression on Road Sensor Data

abstract: The continuous expansion of the urban traffic sensing infrastructure has led to a surge in the volume of widely available road related data. Consequently, increasing effort is being dedicated to the creation of intelligent transportation systems, where decisions on issues ranging from city-wide road maintenance planning to improving the commuting experience are informed by computational models of urban traffic instead of being left entirely to humans. The automation of traffic management has received substantial attention from the research community, however, most approaches target highways, produce predictions valid for a limited time window or require expensive retraining of available models in order to accurately forecast traffic at a new location. In this article, we propose a novel and accurate traffic flow prediction method based on symbolic regression enhanced with a lag operator. Our approach produces robust models suitable for the intricacies of urban roads, much more difficult to predict than highways. Additionally, there is no need to retrain the model for a period of up to 9 weeks. Furthermore, the proposed method generates models that are transferable to other segments of the road network, similar to, yet geographically distinct from the ones they were initially trained on. We demonstrate the achievement of these claims by conducting extensive experiments on data collected from the Darmstadt urban infrastructure.


title: Spin-wave diode and circulator based on unidirectional coupling

abstract: In magnonics, an emerging branch of wave physics characterized by low-energy consumption, it is highly desirable to realize circuit elements within the scope of spin-wave computing. Here, based on numerical simulations, we demonstrate the functionality of the spin-wave diode and the circulator to steer and manipulate spin waves over a wide range of frequency in the GHz regime. They take advantage of the unidirectional coupling induced by the interfacial Dzyaloshinskii-Moriya interaction to transfer the spin wave between thin ferromagnetic layers in only one direction of propagation. Using the multilayered structure consisting of Py and Co in direct contact with Pt, we obtain sub-micrometer-size devices of high efficiency. In the diode, the power loss ratio between forward and reverse direction reaches 22 dB, while in the four-port circulator, the efficiency exceeds 13 dB. Thus, our work contributes to the emerging branch of energy-efficient magnonic logic devices, where, thanks to short wavelength of spin waves, it is possible to realize nanoscale devices.


title: On the gauge group of Galois objects

abstract: We study the Ehresmann--Schauenburg bialgebroid of a noncommutative principal bundle as a quantization of the classical gauge groupoid of a principal bundle. When the base algebra is in the centre of the total space algebra, the gauge group of the noncommutative principal bundle is isomorphic to the group of bisections of the bialgebroid. In particular we consider Galois objects (non-trivial noncommutative bundles over a point in a sense) for which the bialgebroid is a Hopf algebra. For these we give a crossed module structure for the bisections and the automorphisms of the bialgebroid. Examples include Galois objects of group Hopf algebras and of Taft algebras.


title: Consistency Analysis of Replication-Based Probabilistic Key-Value Stores

abstract: Partial quorum systems are widely used in distributed key-value stores due to their latency benefits at the expense of providing weaker consistency guarantees. The probabilistically bounded staleness framework (PBS) studied the latency-consistency trade-off of Dynamo-style partial quorum systems through Monte Carlo event-based simulations. In this paper, we study the latency-consistency trade-off for such systems analytically and derive a closed-form expression for the inconsistency probability. Our approach allows fine-tuning of latency and consistency guarantees in key-value stores, which is intractable using Monte Carlo event-based simulations.


title: Understanding the role of predictive time delay and biased propagator in RAVE

abstract: In this work, we revisit our recent iterative machine learning (ML) -- molecular dynamics (MD) technique "Reweighted autoencoded variational Bayes for enhanced sampling (RAVE)" (Ribeiro, Bravo, Wang, Tiwary, J. Chem. Phys. 149 072301 (2018) and Wang, Ribeiro, Tiwary, Nature Commun. 10 3573 (2019)) and analyze as well as formalize some of its approximations. These including: (a) the choice of a predictive time-delay, or how far into the future should the ML try to predict the state of a given system output from MD, and (b) for short time-delays, how much of an error is made in approximating the biased propagator for the dynamics as the unbiased propagator. We demonstrate through a master equation framework as to why the exact choice of time-delay is irrelevant as long as a small non-zero value is adopted. We also derive a correction to reweight the biased propagator, and somewhat to our dissatisfaction but also to our reassurance, find that it barely makes a difference to the intuitive picture we had previously derived and used.


title: Analyzing Differentiable Fuzzy Logic Operators

abstract: The AI community is increasingly putting its attention towards combining symbolic and neural approaches, as it is often argued that the strengths and weaknesses of these approaches are complementary. One recent trend in the literature are weakly supervised learning techniques that employ operators from fuzzy logics. In particular, these use prior background knowledge described in such logics to help the training of a neural network from unlabeled and noisy data. By interpreting logical symbols using neural networks, this background knowledge can be added to regular loss functions, hence making reasoning a part of learning. We study, both formally and empirically, how a large collection of logical operators from the fuzzy logic literature behave in a differentiable learning setting. We find that many of these operators, including some of the most well-known, are highly unsuitable in this setting. A further finding concerns the treatment of implication in these fuzzy logics, and shows a strong imbalance between gradients driven by the antecedent and the consequent of the implication. Furthermore, we introduce a new family of fuzzy implications (called sigmoidal implications) to tackle this phenomenon. Finally, we empirically show that it is possible to use Differentiable Fuzzy Logics for semi-supervised learning, and compare how different operators behave in practice. We find that, to achieve the largest performance improvement over a supervised baseline, we have to resort to non-standard combinations of logical operators which perform well in learning, but no longer satisfy the usual logical laws.


title: Revisiting the reactivity between HCO and CH$_3$ on interstellar grain surfaces

abstract: Formation of interstellar complex organic molecules is currently thought to be dominated by the barrierless coupling between radicals on the interstellar icy grain surfaces. Previous standard DFT results on the reactivity between CH$_3$ and HCO on amorphous water surfaces, showed that formation of CH$_4$ + CO by H transfer from HCO to CH$_3$ assisted by water molecules of the ice was the dominant channel. However, the adopted description of the electronic structure of the biradical (i.e., CH$_3$/HCO) system was inadequate (without the broken-symmetry (BS) approach). In this work, we revisit the original results by means of BS-DFT both in gas phase and with one water molecule simulating the role of the ice. Results indicate that adoption of BS-DFT is mandatory to describe properly biradical systems. In the presence of the single water molecule, the water-assisted H transfer exhibits a high energy barrier. In contrast, CH$_3$CHO formation is found to be barrierless. However, direct H transfer from HCO to CH$_3$ to give CO and CH$_4$ presents a very low energy barrier, hence being a potential competitive channel to the radical coupling and indicating, moreover, that the physical insights ofthe original work remain valid.


title: Time-Varying Gaussian-Cauchy Mixture Models for Financial Risk Management

abstract: There are various metrics for financial risk, such as value at risk (VaR), expected shortfall, expected/unexpected loss, etc. When estimating these metrics, it was very common to assume Gaussian distribution for the asset returns, which may underestimate the real risk of the market, especially during the financial crisis. In this paper, we propose a series of time-varying mixture models for risk analysis and management. These mixture models contain two components: one component with Gaussian distribution, and the other one with a fat-tailed Cauchy distribution. We allow the distribution parameters and component weights to change over time to increase the flexibility of the models. Monte Carlo Expectation-Maximization algorithm is utilized to estimate the parameters. To verify the good performance of our models, we conduct some simulation studies, and implement our models to the real stock market. Based on these studies, our models are appropriate under different economic conditions, and the component weights can capture the correct pattern of the market volatility.


title: Multivariate Probabilistic Time Series Forecasting via Conditioned Normalizing Flows

abstract: Time series forecasting is often fundamental to scientific and engineering problems and enables decision making. With ever increasing data set sizes, a trivial solution to scale up predictions is to assume independence between interacting time series. However, modeling statistical dependencies can improve accuracy and enable analysis of interaction effects. Deep learning methods are well suited for this problem, but multivariate models often assume a simple parametric distribution and do not scale to high dimensions. In this work we model the multivariate temporal dynamics of time series via an autoregressive deep learning model, where the data distribution is represented by a conditioned normalizing flow. This combination retains the power of autoregressive models, such as good performance in extrapolation into the future, with the flexibility of flows as a general purpose high-dimensional distribution model, while remaining computationally tractable. We show that it improves over the state-of-the-art for standard metrics on many real-world data sets with several thousand interacting time-series.


title: A Numerical Model for the Multi-wavelength Lightcurves of PSR J0030+0451

abstract: Recent modeling of Neutron Star Interior Composition Explorer(NICER) observations of the millisecond pulsar PSR J0030+0451 suggests that the magnetic field of the pulsar is non-dipolar. We construct a magnetic field configuration where foot points of the open field lines closely resemble the hotspot configuration from NICER observations. Using this magnetic field as input, we perform force-free simulations of the magnetosphere of PSR J0030+0451, showing the three-dimensional structure of its plasma-filled magnetosphere. Making simple and physically motivated assumptions about the emitting regions, we are able to construct the multi-wavelength lightcurves that qualitatively agree with the corresponding observations. The agreement suggests that multipole magnetic structures are the key to modeling this type of pulsars, and can be used to constrain the magnetic inclination angle and the location of radio emission.


title: An Observational Study of the Effect of Nike Vaporfly Shoes on Marathon Performance

abstract: We collected marathon performance data from a systematic sample of elite and sub-elite athletes over the period 2015 to 2019, then searched the internet for publicly-available photographs of these performances, identifying whether the Nike Vaporfly shoes were worn or not in each performance. Controlling for athlete ability and race difficulty, we estimated the effect on marathon times of wearing the Vaporfly shoes. Assuming that the effect of Vaporfly shoes is additive, we estimate that the Vaporfly shoes improve men's times between 2.0 and 3.9 minutes, while they improve women's times between 0.8 and 3.5 minutes. Assuming that the effect of Vaporfly shoes is multiplicative, we estimate that they improve men's times between 1.4 and 2.8 percent and women's performances between 0.6 and 2.2 percent. The improvements are in comparison to the shoe the athlete was wearing before switching to Vaporfly shoes, and represents an expected improvement rather than a guaranteed improvement.


title: Three-Dimensional Single-Shot Ptychography

abstract: Here we introduce three-dimensional single-shot ptychography (3DSSP). 3DSSP leverages an additional constraint unique to the single-shot geometry to deconvolve multiple 2D planes of a 3D object. Numeric simulations and analytic calculations demonstrate that 3DSSP reconstructs multiple planes in an extended 3D object with a minimum separation consistent with the depth of field for a conventional microscope. We experimentally demonstrate 3DSSP by reconstructing orthogonal hair strands axially separated by 5 mm. Three-dimensional single-shot ptychography provides a pathway towards volumetric imaging of dynamically evolving systems on ultrafast timescales.


title: Limits of Transmit and Receive Array Gain in Massive MIMO

abstract: In this paper, we consider the transmit and receive antenna array gain of massive MIMO systems. In particular, we look at their dependence on the number of antennas in the array, and the antenna spacing for uniform linear and uniform circular arrays. It is known that the transmit array gain saturates at a certain antenna spacing, but the receive array gain had not been considered. With our physically consistent analysis based on the Multiport Communication Theory, we show that the receive array gain does not saturate, but that there is a peak at a certain antenna spacing when there is no decoupling network at the receiver. As implementing a decoupling network for massive MIMO would be almost impossible, this is a reasonable assumption. Furthermore, we analyze how the array gain changes depending on the antenna spacing and the size of the antenna array and derive design recommendations.


title: Vector perturbations of Kerr-AdS$_5$ and the Painlev\'e VI transcendent

abstract: We analyze the Ansatz of separability for Maxwell equations in generically spinning, five-dimensional Kerr-AdS black holes. We find that the parameter \mu introduced in a previous work by O. Lunin can be interpreted as apparent singularities of the resulting radial and angular equations. Using isomonodromy deformations, we describe a non-linear symmetry of the system, under which \mu is tied to the Painlev\'e VI transcendent. By translating the boundary conditions imposed on the solutions of the equations for quasinormal modes in terms of monodromy data, we find a procedure to fix \mu and study the behavior of the quasinormal modes in the limit of fast spinning small black holes.


title: Diffusion and Consensus in a Weakly Coupled Network of Networks

abstract: We study diffusion and consensus dynamics in a Network of Networks model. In this model, there is a collection of sub-networks, connected to one another using a small number of links. We consider a setting where the links between networks have small weights, or are used less frequently than links within each sub-network. Using spectral perturbation theory, we analyze the diffusion rate and convergence rate of the investigated systems. Our analysis shows that the first order approximation of the diffusion and convergence rates is independent of the topologies of the individual graphs; the rates depend only on the number of nodes in each graph and the topology of the connecting edges. The second order analysis shows a relationship between the diffusion and convergence rates and the information centrality of the connecting nodes within each sub-network. We further highlight these theoretical results through numerical examples.


title: Derived traces of Soergel categories

abstract: We study two kinds of categorical traces of (monoidal) dg categories, with particular interest in categories of Soergel bimodules. First, we explicitly compute the usual Hochschild homology, or derived vertical trace, of the category of Soergel bimodules in arbitrary types. Secondly, we introduce the notion of derived horizontal trace of a monoidal dg category and compute the derived horizontal trace of Soergel bimodules in type A. As an application we obtain a derived annular Khovanov-Rozansky link invariant with an action of full twist insertion, and thus a categorification of the HOMFLY-PT skein module of the solid torus.


title: Shape dependence of renormalized holographic entanglement entropy

abstract: We study the holographic entanglement entropy of deformed entangling regions in three-dimensional CFTs dual to Einstein-AdS gravity, using a renormalization scheme based on the addition of extrinsic counterterms. In this prescription, when even-dimensional manifolds are considered, the universal contribution to the entanglement entropy is identified as the renormalized volume of the Ryu-Takayanagi hypersurface, which is written as the sum of a topological and a curvature term. It is shown that the change in the renormalized entanglement entropy due to the deformation of the entangling surface is encoded purely in the curvature contribution. In turn, as the topological part is given by the Euler characteristic of the Ryu-Takayanagi surface, it remains shape independent. Exploiting the covariant character of the extrinsic counterterms, we apply the renormalization scheme for the case of deformed entangling regions in AdS$_4$/CFT$_3$, recovering the results found in the literature. Finally, we provide a derivation of the relation between renormalized entanglement entropy and Willmore energy. The presence of a lower bound of the latter makes manifest the relation between the AdS curvature of the Ryu-Takayanagi surface and the strong subadditivity property.


title: Formation of morphogenetic patterns in cellular automata

abstract: One of the most important problems in contemporary science, and especially in biology, is to reveal mechanisms of pattern formation. On the level of biological tissues, patterns form due to interactions between cells. These interactions can be long-range if mediated by diffusive molecules or short-range when associated with cell-to-cell contact sites. Mathematical studies of long-range interactions in-volve models based on differential equations while short-range interactions are modelled using discrete type models. In this paper, we use cellular automata (CA) technique to study formation of patterns due to short-range interactions. Namely, we use von Neumann cellular automata represented by a finite set of lattices whose states evolve according to transition rules. Lattices can be considered as representing biological cells (which, in the simplest case, can only be in one of the two different states) while the transition rules define changes in their states due to the cell-to-cell contact interactions. In this model, we identify rules result-ing in the formation of stationary periodic patterns. In our analysis, we distin-guish rules which do not destroy preset patterns and those which cause pattern formation from random initial conditions. Also, we check whether the forming patterns are resistant to noise and analyze the time frame for their formation. Transition rules which allow formation of stationary periodic patterns are then discussed in terms of pattern formation in biology.


title: Diffusion of Chiral Janus Particles in Convection Rolls

abstract: The diffusion of an artificial active particle in a two-dimensional periodic pattern of stationary convection cells is investigated by means of extensive numerical simulations. In the limit of large P\'eclet numbers, i.e., for self-propulsion speeds below a certain depinning threshold and weak roto-translational fluctuations, the particle undergoes asymptotic normal diffusion with diffusion constant proportional to the square root of its diffusion constant at zero flow. Chirality effects in the propulsion mechanism, modeled here by a tunable applied torque, favors particle's jumping between adjacent convection rolls. Roll jumping is signaled by an excess diffusion peak, which appears to separate two distinct active diffusion regimes for low and high chirality. A qualitative interpretation of our simulation results is proposed as a first step toward a fully analytical study of this phenomenon.


title: Inference of solar rotation from perturbations of acoustic mode eigenfunctions

abstract: Today's picture of the internal solar rotation rate profile results essentially from helioseismic analyses of frequency splittings of resonant acoustic waves. Here we present another, complementary estimation of the internal solar rotation rate using the perturbation of the shape of the acoustic waves. For this purpose, we extend a global helioseismic approach developed previously for the investigation of the meridional flow \cite{schad11, schad12, schad13} to work on the components of the differential rotation. We discuss the effect of rotation on mode eigenfunctions and thereon based observables. Based on a numerical study using a simulated rotation rate profile we tailor an inversion approach and also consider the case of the presence of an additional meridional flow. This inversion approach is then applied to data from the MDI (Michelson Doppler Imager aboard the Solar Heliospheric Observatory (SoHO)) instrument and the HMI (Helioseismic and Magnetic Imager aboard the Solar Dynamics Observatory (SDO)) instrument. In the end, rotation rate profiles estimated from eigenfunction perturbation and frequency splittings are compared. The rotation rate profiles from the two different approaches are qualitatively in good agreement, especially for the MDI data. Significant differences are obtained at high latitudes $> 50^{\circ}$ and near the subsurface. The result from HMI data shows larger discrepancies between the different methods. We find that the two global helioseismic approaches provide complementary methods for measuring solar rotation. Comparing the results from different methods may help to reveal systematic influences that affect analyses based on eigenfunction perturbations, like meridional flow measurements.


title: Scalable Neural Methods for Reasoning With a Symbolic Knowledge Base

abstract: We describe a novel way of representing a symbolic knowledge base (KB) called a sparse-matrix reified KB. This representation enables neural modules that are fully differentiable, faithful to the original semantics of the KB, expressive enough to model multi-hop inferences, and scalable enough to use with realistically large KBs. The sparse-matrix reified KB can be distributed across multiple GPUs, can scale to tens of millions of entities and facts, and is orders of magnitude faster than naive sparse-matrix implementations. The reified KB enables very simple end-to-end architectures to obtain competitive performance on several benchmarks representing two families of tasks: KB completion, and learning semantic parsers from denotations.


title: Optimizing Non-Orthogonal Multiple Access in Random Access Networks

abstract: Non-orthogonal multiple access (NOMA) has been considered as a promising solution for improving the spectrum efficiency of next-generation wireless networks. In this paper, the performance of a p-persistent slotted ALOHA system in support of NOMA transmissions is investigated. Specifically, wireless users can choose to use high or low power for data transmissions with certain probabilities. To achieve the maximum network throughput, an analytical framework is developed to analyze the successful transmission probability of NOMA and long term average throughput of users involved in the non-orthogonal transmissions. The feasible region of the maximum number of concurrent users using high and low power to ensure successful NOMA transmissions are quantified. Based on the analysis, an algorithm is proposed to find the optimal transmission probabilities for users to choose high and low power to achieve the maximum system throughput. In addition, the impact of power settings on the network performance is further investigated. Simulations are conducted to validate the analysis.


title: Local continuity of log-concave projection, with applications to estimation under model misspecification

abstract: The log-concave projection is an operator that maps a d-dimensional distribution P to an approximating log-concave density. Prior work by D{\"u}mbgen et al. (2011) establishes that, with suitable metrics on the underlying spaces, this projection is continuous, but not uniformly continuous. In this work we prove a local uniform continuity result for log-concave projection -- in particular, establishing that this map is locally H{\"o}lder-(1/4) continuous. A matching lower bound verifies that this exponent cannot be improved. We also examine the implications of this continuity result for the empirical setting -- given a sample drawn from a distribution P, we bound the squared Hellinger distance between the log-concave projection of the empirical distribution of the sample, and the log-concave projection of P. In particular, this yields interesting statistical results for the misspecified setting, where P is not itself log-concave.


title: Covering of high-dimensional cubes and quantization

abstract: As the main problem, we consider covering of a $d$-dimensional cube by $n$ balls with reasonably large $d$ (10 or more) and reasonably small $n$, like $n=100$ or $n=1000$. We do not require the full coverage but only 90\% or 95\% coverage. We establish that efficient covering schemes have several important properties which are not seen in small dimensions and in asymptotical considerations, for very large $n$. One of these properties can be termed `do not try to cover the vertices' as the vertices of the cube and their close neighbourhoods are very hard to cover and for large $d$ there are far too many of them. We clearly demonstrate that, contrary to a common belief, placing balls at points which form a low-discrepancy sequence in the cube, makes for a very inefficient covering scheme. For a family of random coverings, we are able to provide very accurate approximations to the coverage probability. We then extend our results to the problems of coverage of a cube by smaller cubes and quantization, the latter being also referred to as facility location. Along with theoretical considerations and derivation of approximations, we discuss results of a large-scale numerical investigation.


title: A low cost small reconfigurable autonomous vehicle for research, testing and teaching of guidance, control and positioning systems

abstract: We present a small autonomous vehicle, with the purpose of aiding teachers and researchers to easily deploy different approaches to guidance, control and positioning systems, without the hassle to develop the entire modules required for such endeavor. The design was mainly guided to accomplish a very small cost, along with easy accessibility of the parts used. One of the main problems of using low cost parts is the difficulty of modeling its dynamics, mainly because they lack proper documentation and also because the time varying nature of its characteristics, due to degradation caused by usage. At the same time, any change made to the physical design of the vehicle must be reflected in the mathematical model. We develop a methodology for easy reconfiguration of the vehicle, being able to obtain detailed dynamic mathematical models of the system, as well as a the statistical models needed when implementing stochastic methods such as Kalman Filtering. Baseline example implementations of guidance, control and positioning systems are provided. The hole system is autonomous and self-contained and was tested in a low cost embedded system.


title: A Low-Complexity Framework for Joint User Pairing and Power Control for Cooperative NOMA in 5G and Beyond Cellular Networks

abstract: This paper investigates the performance of cooperative non-orthogonal multiple access (C-NOMA) in a cellular downlink system. The system model consists of a base station (BS) serving multiple users, where users with good channel quality can assist the transmissions between the BS and users with poor channel quality through either half-duplex (HD) or full-duplex (FD) device-to-device (D2D) communications. We formulate and solve a novel optimization problem that jointly determines the optimal D2D user pairing and the optimal power control scheme, where the objective is maximizing the achievable sum rate of the whole system while guaranteeing a certain quality of service (QoS) for all users. The formulated problem is a mixed-integer non-linear program (MINLP) which is generally NPhard. To overcome this issue, we reconstruct the original problem into a bi-level optimization problem that can be decomposed into two sub-problems to be solved independently. The outer problem is a linear assignment problem which can be efficiently handled by the well-known Hungarian method. The inner problem is still a non-convex optimization problem for which finding the optimal solution is challenging. However, we derive the optimal power control policies for both the HD and the FD schemes in closedform expressions, which makes the computational complexity of the inner problems polynomial for every possible pairing configurations. These findings solve ultimately the original MILNP in a timely manner that makes it suitable for real-time and low latency applications. Our simulation results show that the proposed framework outperforms a variety of proposed schemes in the literature and that it can obtain the optimal pairing and power control policies for a network with 100 users in a negligible computational time.


title: The DIRAC code for relativistic molecular calculations

abstract: DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field.


title: Status of the WHIZARD generator for linear colliders

abstract: This summarizes the talk given at the LCWS 2019 conference in Sendai, Japan, on the progress of the WHIZARD event generator in terms of new physics features and technical improvements relevant for the physics programme of future lepton and especially linear colliders. It takes as a reference the version 2.8.2 released in October 2019, and also takes into account the development until version 2.8.3 to be released in February 2020.


title: Light by Light Scattering as a New Probe for Axions

abstract: We study the impact of virtual axions on the polarization of photons inside a cavity during the interaction of high-power laser pulses. A novel detection scheme for measuring the axion-induced ellipticity signal during the Light-by-Light (LBL) scattering process is investigated. We show that a momentum exchange between photons in a probe laser beam and a high-intensity target beam may lead to a resonance at the physical mass of the axion. Consequently, the resonant enhancement of vacuum birefringence gives rise to a large ellipticity signal. This signal enhancement can be applied in order to discriminate between the axion contribution to LBL scattering and the standard model contribution due to electron-positron pairs. The sensitivity of the scheme is studied for experimentally feasible probe light sources and ultrahigh intensity laser backgrounds. It is shown that this technique has the potential to probe the QCD axion in the mass range $10^{-2} \textrm{eV} \lesssim m_{a} \lesssim 1 \textrm{eV}$. In this region the axion induced signal surpasses the standard model background.

categories: hep-ph astro-ph.CO astro-ph.IM hep-ex physics.optics


title: Simulated Annealing with Adaptive Cooling Rates

abstract: As one of the most robust global optimization methods, simulated annealing has received considerable attention, with many variations that attempt to improve the cooling schedule. This paper introduces a variant of simulated annealing that is useful for optimizing atomistic structures, and makes use of the statistical mechanical properties of the system, determined on the fly during optimization, to adaptively control the cooling rate. The adaptive cooling approach is demonstrated to be more computationally efficient than classical simulated annealing, when applied to Lennard-Jones clusters. This increase in efficiency is approximately a factor of two for clusters with 25--40 atoms, and improves with the size of the system.


title: VisMaker: a Question-Oriented Visualization Recommender System for Data Exploration

abstract: The increasingly rapid growth of data production and the consequent need to explore data to obtain answers to the most varied questions have promoted the development of tools to facilitate the manipulation and construction of data visualizations. However, building useful data visualizations is not a trivial task: it may involve a large number of subtle decisions that require experience from their designer. In this paper, we present VisMaker, a visualization recommender tool that uses a set of rules to present visualization recommendations organized and described through questions, in order to facilitate the understanding of the recommendations and assisting the visual exploration process. We carried out two studies comparing our tool with Voyager 2 and analyzed some aspects of the use of tools. We collected feedback from participants to identify the advantages and disadvantages of our recommendation approach. As a result, we gathered comments to help improve the development of tools in this domain.


title: First study of $N_f=2+1+1$ lattice QCD with physical domain-wall quarks

abstract: Using 10-16 units of Nvidia DGX-1, we have generated the first gauge ensemble for $N_f=2+1+1$ lattice QCD with physical $(u/d, s, c)$ domain-wall quarks, on the $ 64^4 $ lattice with lattice spacing $a \sim 0.064 $ fm ($ L > 4 $ fm, and $ M_\pi L > 3 $). The salient feature of this gauge ensemble is that the chiral symmetry is preserved to a high precision and all topological sectors are sampled ergodically. In this paper, we present the first results of the topological susceptibility and the ground-state hadron mass spectra.


title: Tensions in the dark: shedding light on Dark Matter-Dark Energy interactions

abstract: The emergence of an increasingly strong tension between the Hubble rate inferred from early- and late-time observations has reinvigorated interest in nonstandard scenarios, with the aim of reconciling these measurements. One such model involves interactions between Dark Matter and Dark Energy. Here we consider a specific form of the coupling between these two fluids proportional to the Dark Energy energy density, which has been studied extensively in the literature and claimed to substantially alleviate the Hubble tension. We complement the work already discussed in several previous analyses and show that, once all relevant cosmological probes are included simultaneously, the value of the Hubble parameter in this model is $H_0=69.82_{-0.76}^{+0.63}$ km/(s Mpc), which reduces the Hubble tension to $2.5\sigma$. Furthermore, we also perform a statistical model comparison, finding a $\Delta\chi^2$ of $-2.15$ (corresponding to a significance of 1.5$\sigma$) with the inclusion of one additional free parameter, showing no clear preference for this model with respect to $\Lambda$CDM, which is further confirmed with an analysis of the Bayes ratio.


title: Multiplicity-induced-dominancy for delay-differential equations of retarded type

abstract: An important question of ongoing interest for linear time-delay systems is to provide conditions on its parameters guaranteeing exponential stability of solutions. Recent works have explored spectral techniques to show that, for some low-order delay-differential equations of retarded type, spectral values of maximal multiplicity are dominant, and hence determine the asymptotic behavior of the system, a property known as multiplicity-induced-dominancy. This work further explores such a property and shows its validity for general linear delay-differential equations of retarded type of arbitrary order including a single delay in the system's representation. More precisely, an interesting link between characteristic functions with a real root of maximal multiplicity and Kummer's confluent hypergeometric functions is exploited. We also provide examples illustrating our main result.


title: Deploying large fixed file datasets with SquashFS and Singularity

abstract: Shared high-performance computing (HPC) platforms, such as those provided by XSEDE and Compute Canada, enable researchers to carry out large-scale computational experiments at a fraction of the cost of the cloud. Most systems require the use of distributed filesystems (e.g. Lustre) for providing a highly multi-user, large capacity storage environment. These suffer performance penalties as the number of files increases due to network contention and metadata performance. We demonstrate how a combination of two technologies, Singularity and SquashFS, can help developers, integrators, architects, and scientists deploy large datasets (O(10M) files) on these shared systems with minimal performance limitations. The proposed integration enables more efficient access and indexing than normal file-based dataset installations, while providing transparent file access to users and processes. Furthermore, the approach does not require administrative privileges on the target system. While the examples studied here have been taken from the field of neuroimaging, the technologies adopted are not specific to that field. Currently, this solution is limited to read-only datasets. We propose the adoption of this technology for the consumption and dissemination of community datasets across shared computing resources.


title: Thermodynamic Evidence for a Two-Component Superconducting Order Parameter in Sr$_2$RuO$_4$

abstract: Sr$_2$RuO$_4$ has stood as the leading candidate for a spin-triplet superconductor for 26 years. Recent NMR experiments have cast doubt on this candidacy, however, and it is difficult to find a theory of superconductivity that is consistent with all experiments. What is needed are symmetry-based experiments that can rule out broad classes of possible superconducting order parameters. Here we use resonant ultrasound spectroscopy to measure the entire symmetry-resolved elastic tensor of Sr$_2$RuO$_4$ through the superconducting transition. We observe a thermodynamic discontinuity in the shear elastic modulus $c_{66}$, requiring that the superconducting order parameter is two-component. A two-component $p$-wave order parameter, such as $p_x+i p_y$, naturally satisfies this requirement. As this order parameter appears to be precluded by recent NMR experiments, we suggest that two other two-component order parameters, namely $\left\{d_{xz},d_{yz}\right\}$ or $\left\{d_{x^2-y^2},g_{xy(x^2-y^2)}\right\}$, are now the prime candidates for the order parameter of Sr$_2$RuO$_4$.


title: Perturbative BF theory

abstract: We consider a superrenomalizable gauge theory of topological type, in which the structure group is equal to the inhomogeneous group ISU(2). The generating functional of the correlation functions of the gauge fields is derived and its connection with the generating functional of the Chern-Simons theory is discussed. The complete renomalization of this model defined in R^3 is presented. The structure of the ISU(2) conjugacy classes is determined. Gauge invariant observables are defined by means of appropriately normalized traces of ISU(2) holonomies associated with oriented, framed and coloured knots. The perturbative evaluation of the Wilson lines expectation values is investigated and the up-to-third-order contributions to the perturbative expansion of the observables, which correspond to knot invariants, are produced. The general dependence of the knot observables on the framing is worked out.


title: A two-fluid model for black-hole accretion flows: particle acceleration and disc structure

abstract: Hot, tenuous advection-dominated accretion flows around black holes are ideal sites for the Fermi acceleration of relativistic particles at standing shock waves in the accretion disc. Previous work has demonstrated that the shock-acceleration process can be efficient enough to power the observed, strong outflows in radio-loud active galaxies such as M87. However, the dynamical effect (back-reaction) on the flow, exerted by the pressure of the relativistic particles, has not been previously considered, and this effect can have a significant influence on the disc structure. We reexamine the problem by developing a new, two-fluid model for the structure of the accretion disc that includes the dynamical effect of the relativistic particle pressure, combined with the pressure of the background (thermal) gas. The new model is analogous to the two-fluid model of cosmic ray acceleration in supernova-driven shock waves. As part of the model, we also develop a new set of shock jump conditions, which are solved along with the hydrodynamic conservation equations to determine the structure of the accretion disc. The solutions include the formation of a mildly relativistic outflow (jet) at the shock radius, driven by the relativistic particles accelerated in the disc. One of our main conclusions is that in the context of the new two-fluid accretion model, global smooth (shock-free) solutions do not exist, and the disc must always contain a standing shock wave, at least in the inviscid case considered here.


title: Correlation Measurements of Charged Particles and Jets at Mid-Rapidity with Event Activity at Backward-Rapidity in $\sqrt{s_\mathrm{NN}}=200\,\mathrm{GeV}$ $p$+Au Collisions at STAR

abstract: Semi-inclusive charged jet spectra per trigger at STAR are presented binned by event activity (EA) as determined by the Beam Beam Counter (BBC) signal in the Au-going direction. The selected EA determination is motivated by correlations between the number of charged tracks in the Time Projection Chamber (TPC) ($|\eta|<1$) and EA ($\eta_\mathrm{EA}\in[-5,-2]$) which are also presented. The jet spectra per trigger at high EA are suppressed relative to the spectra at low EA. A PYTHIA investigation refutes that the suppression results from a trivial autocorrelation between jet kinematics and the acceptance of the EA and the TPC.


title: Non-equilibrium thermodynamics of quantum processes assisted by transitionless quantum driving: the role of initial state preparation

abstract: Adiabatic evolution is considered to be the ideal situation for most thermodynamic cycles, as it allows for the achievement of maximum efficiency. However, no power output is produced in light of the infinite amount of time required to perform a transformation. Such issue can be overcome through Shortcuts-To-Adiabaticity (STA) protocols, which allows a dynamics to mimic its adiabatic counterpart, but in a finite time. Transitionless quantum driving (TQD) is one form of STA. We develop and study the effects of TQD on a single-qubit system subjected to a changing magnetic field and a general two-qubit Heisenberg model with a time-dependent interaction strength. We establish a quantitative relation between the work produced across the dynamics and the entropy production resulting from the evolution at hand, focusing on the role played by the initial states of the work medium. We identify the states that extremize performance under the TQD protocol. The thermodynamic implications of STA are discussed in the end.


title: Multivariate Monotone Inclusions in Saddle Form

abstract: We propose a novel approach to monotone operator splitting based on the notion of a saddle operator. Under investigation is a highly structured multivariate monotone inclusion problem involving a mix of set-valued, cocoercive, and Lipschitzian monotone operators, as well as various monotonicity-preserving operations among them. This model encompasses most formulations found in the literature. A limitation of existing primal-dual algorithms is that they operate in a product space that is too small to achieve full splitting of our problem in the sense that each operator is used individually. To circumvent this difficulty, we recast the problem as that of finding a zero of a saddle operator that acts on a bigger space. This leads to an algorithm of unprecedented flexibility, which achieves full splitting, exploits the specific attributes of each operator, is asynchronous, and requires to activate only blocks of operators at each iteration, as opposed to activating all of them. The latter feature is of critical importance in large-scale problems. Weak convergence of the main algorithm is established, as well as the strong convergence of a variant. Various applications are discussed, and instantiations of the proposed framework in the context of variational inequalities and minimization problems are presented.


title: On residually finite groups satisfying an Engel type identity

abstract: Let $ n, q $ be positive integers. We show that if $ G $ is a finitely generated residually finite group satisfying the identity $ [x,_ny^q]\equiv 1, $ then there exists a function $ f(n) $ such that $ G $ has a nilpotent subgroup of finite index of class at most $ f(n) $. We also extend this result to locally graded groups.


title: RL agents Implicitly Learning Human Preferences

abstract: In the real world, RL agents should be rewarded for fulfilling human preferences. We show that RL agents implicitly learn the preferences of humans in their environment. Training a classifier to predict if a simulated human's preferences are fulfilled based on the activations of a RL agent's neural network gets .93 AUC. Training a classifier on the raw environment state gets only .8 AUC. Training the classifier off of the RL agent's activations also does much better than training off of activations from an autoencoder. The human preference classifier can be used as the reward function of an RL agent to make RL agent more beneficial for humans.


title: A two-fluid model for black-hole accretion flows: Particle acceleration, outflows, and TeV emission

abstract: The multi-wavelength spectrum observed from M87 extends from radio wavelengths up to TeV gamma-ray energies. The radio through GeV components have been interpreted successfully using SSC models based on misaligned blazar jets, but the origin of the intense TeV emission detected during flares in 2004, 2005, and 2010 remains puzzling. It has been previously suggested that the TeV flares are produced when a relativistic proton jet originating in the core of M87 collides with a molecular cloud (or stellar atmosphere) located less than one parsec from the central black hole. We explore this scenario in detail here using a self-consistent model for the acceleration of relativistic protons in a shocked, two-fluid ADAF accretion disc. The relativistic protons accelerated in the disc escape to power the observed jet outflows. The distribution function for the jet protons is used to compute the TeV emission produced when the jet collides with a cloud or stellar atmosphere. The simulated broadband radiation spectrum includes radio, X-ray, and GeV components generated via synchrotron, as well as TeV emission generated via the production and decay of muons, positrons, and electrons. The self-consistency of the model is verified by computing the relativistic particle pressure using the distribution function, and comparing it with the relativistic particle pressure obtained from the hydrodynamical model. We demonstrate that the model is able to reproduce the multi-wavelength spectrum from M87 observed by VERITAS and HESS during the high-energy flares in 2004, 2005, and 2010.


title: HDG and CG methods for the Indefinite Time-Harmonic Maxwell's Equations under minimal regularity

abstract: We propose to use a hybridizable discontinuous Galerkin (HDG) method combined with the continuous Galerkin (CG) method to approximate Maxwell's equations. We make two contributions in this paper. First, even though there are many papers using HDG methods to approximate Maxwell's equations, to our knowledge they all assume that the coefficients are smooth (or constant). Here, we derive optimal convergence estimates for our HGD-CG approximation when the electromagnetic coefficients are {\em piecewise} smooth. This requires new techniques of analysis. Second, we use CG elements to approximate the Lagrange multiplier used to enforce the divergence condition and we obtain a discrete system in which we can decouple the discrete the Lagrange multiplier. Because we are using a continuous Lagrange multiplier space, the number of degrees of freedom devoted to this are less than for other HDG methods. We present numerical experiments to confirm our theoretical results.


title: Period Relations for Quaternionic Elliptic Functions

abstract: We study elliptic functions in quaternionic analysis, and prove some analogues of classical theorems from the complex case. The main result is a relation between the periods of closed differential 1-forms and 3-forms on H/L where L is a lattice in H. By applying to forms that are regular in the sense of quaternionic analysis, we obtain quaternionic analogues of Riemann's period relations for an abelian surface. We also obtain a quaternionic analogue of Legendre's relation, as an equality between $2\pi^2$ and a linear combination of the four quasi-periods of the quaternionic Weierstrass $\zeta$-function, where the coefficients are rational functions of a basis for L. When L is a left-ideal for a maximal order in a definite quaternion algebra over Q, we obtain three further relations among the quasi-periods.


title: Combining Parametric Land Surface Models with Machine Learning

abstract: A hybrid machine learning and process-based-modeling (PBM) approach is proposed and evaluated at a handful of AmeriFlux sites to simulate the top-layer soil moisture state. The Hybrid-PBM (HPBM) employed here uses the Noah land-surface model integrated with Gaussian Processes. It is designed to correct the model only in climatological situations similar to the training data else it reverts to the PBM. In this way, our approach avoids bad predictions in scenarios where similar training data is not available and incorporates our physical understanding of the system. Here we assume an autoregressive model and obtain out-of-sample results with upwards of a 3-fold reduction in the RMSE using a one-year leave-one-out cross-validation at each of the selected sites. A path is outlined for using hybrid modeling to build global land-surface models with the potential to significantly outperform the current state-of-the-art.


title: Lift and drag forces acting on a particle moving with zero slip velocity near a wall

abstract: The lift and drag forces acting on a small, neutrally-buoyant spherical particle in a single-wall-bounded linear shear flow are examined via numerical computation. The effects of shear rate are isolated from those of slip by setting the particle velocity equal to the local fluid velocity (zero slip), and examining the resulting hydrodynamic forces as a function of separation distance. In contrast to much of the previous numerical literature, low shear Reynolds numbers are considered ($10^{-3} \lesssim Re_{\gamma} \lesssim 10^{-1}$). This shear rate range is relevant when dealing with particulate flows within small channels, for example particle migration in microfluidic devices being used or developed for the biotech industry. We demonstrate a strong dependence of both the lift and drag forces on shear rate. Building on previous theoretical $Re_{\gamma} \ll 1$ studies, a wall-shear based lift correlation is proposed that is applicable when the wall lies both within the inner and outer regions of the disturbed flow. Similarly, we validate an improved drag correlation that includes higher order terms in wall separation distance that more accurately captures the drag force when the particle is close to, but not touching, the wall. Application of the new correlations shows that the examined shear based lift force is as important as the previously examined slip based lift force, highlighting the need to account for shear when predicting the near-wall movement of neutrally-buoyant particles.


title: Are deviations in a gradually varying mean relevant? A testing approach based on sup-norm estimators

abstract: Classical change point analysis aims at (1) detecting abrupt changes in the mean of a possibly non-stationary time series and at (2) identifying regions where the mean exhibits a piecewise constant behavior. In many applications however, it is more reasonable to assume that the mean changes gradually in a smooth way. Those gradual changes may either be non-relevant (i.e., small), or relevant for a specific problem at hand, and the present paper presents statistical methodology to detect the latter. More precisely, we consider the common nonparametric regression model $X_{i} = \mu (i/n) + \varepsilon_{i}$ with possibly non-stationary errors and propose a test for the null hypothesis that the maximum absolute deviation of the regression function $\mu$ from a functional $g (\mu )$ (such as the value $\mu (0)$ or the integral $\int_{0}^{1} \mu (t) dt$) is smaller than a given threshold on a given interval $[x_{0},x_{1}] \subseteq [0,1]$. A test for this type of hypotheses is developed using an appropriate estimator, say $\hat d_{\infty, n}$, for the maximum deviation $ d_{\infty}= \sup_{t \in [x_{0},x_{1}]} |\mu (t) - g( \mu) |$. We derive the limiting distribution of an appropriately standardized version of $\hat d_{\infty,n}$, where the standardization depends on the Lebesgue measure of the set of extremal points of the function $\mu(\cdot)-g(\mu)$. A refined procedure based on an estimate of this set is developed and its consistency is proved. The results are illustrated by means of a simulation study and a data example.


title: Combining Visual and Textual Features for Semantic Segmentation of Historical Newspapers

abstract: The massive amounts of digitized historical documents acquired over the last decades naturally lend themselves to automatic processing and exploration. Research work seeking to automatically process facsimiles and extract information thereby are multiplying with, as a first essential step, document layout analysis. If the identification and categorization of segments of interest in document images have seen significant progress over the last years thanks to deep learning techniques, many challenges remain with, among others, the use of finer-grained segmentation typologies and the consideration of complex, heterogeneous documents such as historical newspapers. Besides, most approaches consider visual features only, ignoring textual signal. In this context, we introduce a multimodal approach for the semantic segmentation of historical newspapers that combines visual and textual features. Based on a series of experiments on diachronic Swiss and Luxembourgish newspapers, we investigate, among others, the predictive power of visual and textual features and their capacity to generalize across time and sources. Results show consistent improvement of multimodal models in comparison to a strong visual baseline, as well as better robustness to high material variance.


title: Purifying Real Images with an Attention-guided Style Transfer Network for Gaze Estimation

abstract: Recently, the progress of learning-by-synthesis has proposed a training model for synthetic images, which can effectively reduce the cost of human and material resources. However, due to the different distribution of synthetic images compared to real images, the desired performance cannot be achieved. Real images consist of multiple forms of light orientation, while synthetic images consist of a uniform light orientation. These features are considered to be characteristic of outdoor and indoor scenes, respectively. To solve this problem, the previous method learned a model to improve the realism of the synthetic image. Different from the previous methods, this paper try to purify real image by extracting discriminative and robust features to convert outdoor real images to indoor synthetic images. In this paper, we first introduce the segmentation masks to construct RGB-mask pairs as inputs, then we design a attention-guided style transfer network to learn style features separately from the attention and bkgd(background) region , learn content features from full and attention region. Moreover, we propose a novel region-level task-guided loss to restrain the features learnt from style and content. Experiments were performed using mixed studies (qualitative and quantitative) methods to demonstrate the possibility of purifying real images in complex directions. We evaluate the proposed method on three public datasets, including LPW, COCO and MPIIGaze. Extensive experimental results show that the proposed method is effective and achieves the state-of-the-art results.


title: On qualitative properties of single-delay linear retarded differential equations: Characteristic roots of maximal multiplicity are necessarily dominant

abstract: This paper presents necessary and sufficient conditions for the existence of a real root of maximal multiplicity in the spectrum of a linear time-invariant single-delay equation of retarded type. We also prove that this root is always strictly dominant, and hence determines the asymptotic behavior of the system. These results are based on improved a priori bounds on the imaginary part of roots on the complex right half-plane.


title: Sharp inequalities on circular and hyperbolic functions using a Bernoulli inequality type

abstract: In this paper, new sharp bounds for circular functions are proved. We provide some improvements of previous results by using infinite products, power series expansions and a generalisation of the so-called Bernoulli inequality. New proofs, refinements as well as new results are offered.


title: Expanding theory testing in general relativity: LIGO and parametrized theories

abstract: The multiple detections of gravitational waves by LIGO (the Laser Interferometer Gravitational-Wave Observatory), operated by Caltech and MIT, have been acclaimed as confirming Einstein's prediction, a century ago, that gravitational waves propagating as ripples in spacetime would be detected. Yunes and Pretorius (2009) investigate whether LIGO's template-based searches encode fundamental assumptions, especially the assumption that the background theory of general relativity is an accurate description of the phenomena detected in the search. They construct the parametrized post-Einsteinian (ppE) framework in response, which broadens those assumptions and allows for wider testing under more flexible assumptions. Their methods are consistent with work on confirmation and testing found in Carnap (1936), Hempel (1969), and Stein (1992, 1994), with the following principles in common: that confirmation is distinct from testing, and that, counterintuitively, revising a theory's formal basis can make it more broadly empirically testable. These views encourage a method according to which theories can be made abstract, to define families of general structures for the purpose of testing. With the development of the ppE framework and related approaches, multi-messenger astronomy is a catalyst for deep reasoning about the limits and potential of the theoretical framework of general relativity.


title: Dual-cyclic polytopes of convex planar polygons with fixed vertex angles

abstract: If we fix the angles at the vertices of a convex planar $n$-gon, the lengths of its edges must satisfy two linear constraints in order for it to close up. If we also require unit perimeter, our vectors of $n$ edge lengths form a convex polytope of dimension $n-3$, each facet of which consists of those $n$-gons in which the length of a particular edge has fallen to zero. Bavard and Ghys require unit area instead, which gives them a hyperbolic polytope. Those two polytopes are combinatorially equivalent, so either is fine for our purposes. Such a fixed-angles polytope is combinatorially richer when the angles are well balanced. We say that fixed external angles are "majority dominant" when every consecutive string of more than half of them sums to more than $\pi$. When $n$ is odd, we show that the fixed-angles polytope for any majority-dominant angles is dual to the cyclic polytope $C_{n-3}(n)$. To extend that result to even $n$, we require that the angles also have "dipole tie-breaking": None of the $n$ strings of length $n/2$ sums to precisely $\pi$, and the $n/2$ that sum to more than $\pi$ overlap as much as possible, all containing a particular angle. Fixing the vertex angles is uncommon, however; people more often fix the edge lengths. That is harder, in part because fixed-lengths $n$-gons may not be convex, but mostly because fixing the lengths constrains the angles nonlinearly -- so the resulting moduli spaces, called "polygon spaces", are curved. Using Schwarz-Christoffel maps, Kapovich and Millson show that the subset of that polygon space in which the $n$-gons are convex and traversed counterclockwise is homeomorphic to the fixed-angles polytope above, for those same fixed values. Each such subset is thus a topological polytope; and it is dual cyclic whenever the fixed lengths are majority dominant and, for even $n$, have dipole tie-breaking.


title: Navier-Stokes equations: an analysis of a possible gap to achieve the energy equality

abstract: The paper is concerned with the IBVP of the Navier-Stokes equations. The goal is the attempt to construct a weak solution enjoying an energy equality. This result is a natural continuation and improvement of the one obtained by the same authors in [2].


title: Theoretical treatment of Coriolis effect using hyper-spherical coordinates, with application to the ro-vibrational spectrum of ozone

abstract: Several alternative methods for description of interaction between rotation and vibration are compared and contrasted using hyper-spherical coordinates for a triatomic molecule. These methods differ by the choice of z-axis and by the assumption of a prolate or oblate rotor shape of the molecule. For each case, block-structure of the rotational-vibrational Hamiltonian matrix is derived and analyzed, and the advantages and disadvantages of each method are made explicit. This theory is then employed to compute ro-vibrational spectra of singly substituted ozone; roughly, 600 vibrational states of 16O18O16O and 16O16O18O isomers combined, with rotational excitations up to J=5 and both inversion parities (21600 coupled ro-vibrational states total). Splittings between the states of different parities, so called K-doublings, are calculated and analyzed. The roles of the asymmetric-top rotor term and the Coriolis coupling term are determined individually, and it is found that they both affect these splittings, but in the opposite directions. Thus, the two effects partially cancel out, and the residual splittings are relatively small. Energies of the ro-vibrational states reported in this work for 16O18O16O and 16O16O18O are in excellent agreement with literature (available for low vibrational excitation). New data obtained here for the highly excited vibrational states enable the first systematic study of the Coriolis effect in symmetric and asymmetric isotopomers of ozone.


title: Analysis of the acoustic waves reflected by a cluster of small holes in the time-domain and the equivalent mass density

abstract: We study the time-domain acoustic scattering problem by a cluster of small holes (i.e. sound-soft obstacles). Based on the retarded boundary integral equation method, we derive the asymptotic expansion of the scattered field as the size of the holes goes to zero. Under certain geometrical constraints on the size and the minimum distance of the holes, we show that the scattered field is approximated by a linear combination of point-sources where the weights are given by the capacitance of each hole and the causal signals (of these point-sources) can be computed by solving a, retarded in time, linear algebraic system. A rigorous justification of the asymptotic expansion and the unique solvability of the linear algebraic system are shown under natural conditions on the cluster of holes. As an application of the asymptotic expansion, we derive, in the limit case when the holes are densely distributed and occupy a bounded domain, the equivalent effective acoustic medium (an equivalent mass density characterized by the capacitance of the holes) that generates, approximately, the same scattered field as the cluster of holes. Conversely, given a locally variable, smooth and positive mass density, satisfying a certain subharmonicity condition, we can design a perforated material with holes, having appropriate capacitances, that generates approximately the same acoustic field as the acoustic medium modelled by the given mass density (and constant speed of propagation). Finally, we numerically verify the asymptotic expansions by comparing the asymptotic approximations with the numerical solutions of the scattered fields via the finite element method.


title: Geodesic rays of the $N$-body problem

abstract: For the Newtonian N-body problem, we study the Jacobi-Maupertuis metric of the nonnegative energy levels. We show that the geodesic rays are expansive, that is to say, all the distances between the bodies must be divergent functions. More precisely, we prove that the evolution of such motions asymptotically decomposes into free particles and subsystems in completely parabolic expansion. The theorem applies in particular to the maximal characteristic curves of any given global viscosity solution of the stationary Hamilton-Jacobi equation $H(x,d_xu) = h$.


title: Epsilon local rigidity and numerical algebraic geometry

abstract: A well-known combinatorial algorithm can decide generic rigidity in the plane by determining if the graph is of Pollaczek-Geiringer-Laman type. Methods from matroid theory have been used to prove other interesting results, again under the assumption of generic configurations. However, configurations arising in applications may not be generic. We present Theorem 5 and its corresponding Algorithm 1 which decide if a configuration is epsilon-locally rigid, a notion we define. A configuration which is epsilon-locally rigid may be locally rigid or flexible, but any continuous deformations remain within a sphere of radius epsilon in configuration space. Deciding epsilon-local rigidity is possible for configurations which are smooth or singular, generic or non-generic. We also present Algorithms 2 and 3 which use numerical algebraic geometry to compute a discrete-time sample of a continuous flex, providing useful visual information for the scientist.


title: U.S. Test System with High Spatial and Temporal Resolution for Renewable Integration Studies

abstract: Planning for power systems with high penetrations of variable renewable energy requires higher spatial and temporal granularity. However, most publicly available test systems are of insufficient fidelity for developing methods and tools for high-resolution planning. This paper presents methods to construct open-access test systems of high spatial resolution to more accurately represent infrastructure and high temporal resolution to represent dynamics of demand and variable resources. To demonstrate, a high-resolution test system representing the United States is created using only publicly available data. This test system is validated by running it in a production cost model, with results compared against historical generation to ensure that they are representative. The resulting open source test system can support power system transition planning and aid in development of tools to answer questions around how best to reach decarbonization goals, using the most effective combinations of transmission expansion, renewable generation, and energy storage.


title: On the role of rotation-vibration coupling in the spectra of ozone isotopomers

abstract: A theoretical framework and computer code are developed for accurate calculations of coupled rotational-vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulae are derived for construction of the Hamiltonian matrix using an efficient combination of the FBR and DVR methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick. Then, accurate calculations of the ro-vibrational spectra are carried out for doubly substituted symmetric (18O16O18O) and asymmetric (18O18O16O) ozone isotopomers for total angular momentum up to J=5. Together with similar data recently reported for the singly substituted symmetric (16O18O16O) and asymmetric (16O16O18O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, by the magnitudes of parity splittings ({\Lambda}-doubling), and by the ratios of ro-vibrational partition functions for asymmetric vs. symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.


title: Generalization and Representational Limits of Graph Neural Networks

abstract: We address two fundamental questions about graph neural networks (GNNs). First, we prove that several important graph properties cannot be computed by GNNs that rely entirely on local information. Such GNNs include the standard message passing models, and more powerful spatial variants that exploit local graph structure (e.g., via relative orientation of messages, or local port ordering) to distinguish neighbors of each node. Our treatment includes a novel graph-theoretic formalism. Second, we provide the first data dependent generalization bounds for message passing GNNs. This analysis explicitly accounts for the local permutation invariance of GNNs. Our bounds are much tighter than existing VC-dimension based guarantees for GNNs, and are comparable to Rademacher bounds for recurrent neural networks.


title: Characterizing low contrast Galactic open clusters with GAIA DR2

abstract: In this study, we characterized 16 objects previously classified as faint or low contrast Galactic open clusters (OCs). We employed parameters associated to the OCs dynamical evolution: core ($r_c$), tidal ($r_t$) and half-mass ($r_{hm}$) radii, age and crossing time ($t_{cr}$). Relations among these parameters were exploited to draw some evolutionary connections. We also included 11 OCs with previous characterizations to provide wider coverage of the parameters space. The investigated sample spans a considerable range in age (log (t/yr) ~7.0 - 9.7) and Galactocentric distance ($R_G$ ~ 6 - 11 kpc). Most of them present solar metallicity. We employed GAIA DR2 astrometry and photometry and selected member stars through a decontamination algorithm which explores the 3D astrometric space ($\mu_{\alpha}$, $\mu_{\delta}$, $\varpi$) to assign membership likelihoods. Previous studies of most of these objects were based mostly on photometric information. All investigated OCs were proved to be real stellar concentrations and relations among their parameters indicate a general disruption scenario in which OCs tend to be more concentrated as they evolve. Internal interactions sucessively drive OCs to develop more dynamically relaxed structures and make them less subject to mass loss due to tidal effects. Tidal radius tends to increase with $R_G$ in accordance with the strength of the Galactic tidal field. Besides, the correlation between the $r_c$ and the dynamical ratio $\tau_{\textrm{dyn}}$ = age/$t_{cr}$ suggests two distinct evolutionary sequences, which may be consequence of different initial formation conditions.


title: Statistical Monitoring of the Covariance Matrix in Multivariate Processes: A Literature Review

abstract: Monitoring several correlated quality characteristics of a process is common in modern manufacturing and service industries. Although a lot of attention has been paid to monitoring the multivariate process mean, not many control charts are available for monitoring the covariance matrix. This paper presents a comprehensive overview of the literature on control charts for monitoring the covariance matrix in a multivariate statistical process monitoring (MSPM) framework. It classifies the research that has previously appeared in the literature. We highlight the challenging areas for research and provide some directions for future research.


title: The Phantom Steering Effect in Q&A Websites

abstract: Badges are commonly used in online platforms as incentives for promoting contributions. It is widely accepted that badges "steer" people's behavior toward increasing their rate of contributions before obtaining the badge. This paper provides a new probabilistic model of user behavior in the presence of badges. By applying the model to data from thousands of users on the Q&A site Stack Overflow, we find that steering is not as widely applicable as was previously understood. Rather, the majority of users remain apathetic toward badges, while still providing a substantial number of contributions to the site. An interesting statistical phenomenon, termed "Phantom Steering," accounts for the interaction data of these users and this may have contributed to some previous conclusions about steering. Our results suggest that a small population, approximately 20%, of users respond to the badge incentives. Moreover, we conduct a qualitative survey of the users on Stack Overflow which provides further evidence that the insights from the model reflect the true behavior of the community. We argue that while badges might contribute toward a suite of effective rewards in an online system, research into other aspects of reward systems such as Stack Overflow reputation points should become a focus of the community.


title: menoci: Lightweight Extensible Web Portal enabling FAIR Data Management for Biomedical Research Projects

abstract: Background: Biomedical research projects deal with data management requirements from multiple sources like funding agencies' guidelines, publisher policies, discipline best practices, and their own users' needs. We describe functional and quality requirements based on many years of experience implementing data management for the CRC 1002 and CRC 1190. A fully equipped data management software should improve documentation of experiments and materials, enable data storage and sharing according to the FAIR Guiding Principles while maximizing usability, information security, as well as software sustainability and reusability. Results: We introduce the modular web portal software menoci for data collection, experiment documentation, data publication, sharing, and preservation in biomedical research projects. Menoci modules are based on the Drupal content management system which enables lightweight deployment and setup, and creates the possibility to combine research data management with a customisable project home page or collaboration platform. Conclusions: Management of research data and digital research artefacts is transforming from individual researcher or groups best practices towards project- or organisation-wide service infrastructures. To enable and support this structural transformation process, a vital ecosystem of open source software tools is needed. Menoci is a contribution to this ecosystem of research data management tools that is specifically designed to support biomedical research projects.


title: The Breakthrough Listen Search for Intelligent Life: A 3.95-8.00 GHz Search for Radio Technosignatures in the Restricted Earth Transit Zone

abstract: We report on a search for artificial narrowband signals of 20 stars within the restricted Earth Transit Zone as a part of the ten-year Breakthrough Listen (BL) search for extraterrestrial intelligence. The restricted Earth Transit Zone is the region of the sky from which an observer would see the Earth transit the Sun with an impact parameter of less than 0.5. This region of the sky is geometrically unique, providing a potential way for an extraterrestrial intelligence to discover the Solar System. The targets were nearby (7-143 pc) and the search covered an electromagnetic frequency range of 3.95-8.00 GHz. We used the Robert C. Byrd Green Bank Telescope to perform these observations with the standard BL data recorder. We searched these data for artificial narrowband ($\sim$Hz) signals with Doppler drift rates of $\pm 20$ Hz s$^{-1}$. We found one set of potential candidate signals on the target HIP 109656 which was then found to be consistent with known properties of anthropogenic radio frequency interference. We find no evidence for radio technosignatures from extraterrestrial intelligence in our observations. The observing campaign achieved a minimum detectable flux which would have allowed detections of emissions that were $10^{-3}$ to $0.88$ times as powerful as the signaling capability of the Arecibo radar transmitter, for the nearest and furthest stars respectively. We conclude that at least $8\%$ of the systems in the restricted Earth Transit Zone within 150 pc do not possess the type of transmitters searched in this survey. To our knowledge, this is the first targeted search for extraterrestrial intelligence of the restricted Earth Transit Zone. All data used in this paper are publicly available via the Breakthrough Listen Public Data Archive (http://seti.berkeley.edu/bldr2).


title: Cleaning Denial Constraint Violations through Relaxation

abstract: Data cleaning is a time-consuming process that depends on the data analysis that users perform. Existing solutions treat data cleaning as a separate offline process that takes place before analysis begins. Applying data cleaning before analysis assumes a priori knowledge of the inconsistencies and the query workload, thereby requiring effort on understanding and cleaning the data that is unnecessary for the analysis. We propose an approach that performs probabilistic repair of denial constraint violations on-demand, driven by the exploratory analysis that users perform. We introduce Daisy, a system that seamlessly integrates data cleaning into the analysis by relaxing query results. Daisy executes analytical query-workloads over dirty data by weaving cleaning operators into the query plan. Our evaluation shows that Daisy adapts to the workload and outperforms traditional offline cleaning on both synthetic and real-world workloads.


title: The positive tropical Grassmannian, the hypersimplex, and the m=2 amplituhedron

abstract: The study of the moment map from the Grassmannian to the hypersimplex, and the relation between torus orbits and matroid polytopes, dates back to the foundational 1987 work of Gelfand-Goresky-MacPherson-Serganova. On the other hand, the amplituhedron is a very new object, defined by Arkani-Hamed-Trnka in connection with scattering amplitudes in $\mathcal{N}=4$ super Yang-Mills theory. In this paper we discover a striking duality between the moment map $\mu:Gr^{\geq0}_{k+1,n}\to\Delta_{k+1,n}$ from the positive Grassmannian $Gr^{\geq0}_{k+1,n}$ to the hypersimplex, and the amplituhedron map $\tilde{Z}:Gr^{\geq0}_{k,n}\to\mathcal{A}_{n,k,2}(Z)$ from $Gr^{\geq0}_{k,n}$ to the $m=2$ amplituhedron. We consider the positroid dissections of both objects, which informally, are subdivisions of $\Delta_{k+1,n}$ (respectively, $\mathcal{A}_{n,k,2}(Z)$) into a disjoint union of images of positroid cells of the positive Grassmannian. At first glance, $\Delta_{k+1,n}$ and $\mathcal{A}_{n,k,2}(Z)$ seem very different - the former is an $(n-1)$-dimensional polytope, while the latter is a $2k$-dimensional non-polytopal subset of $Gr_{k,k+2}$. Nevertheless, we conjecture that positroid dissections of $\Delta_{k+1,n}$ are in bijection with positroid dissections of $\mathcal{A}_{n,k,2}(Z)$ via a map we call T-duality. We prove this conjecture for the (infinite) class of BCFW dissections and give additional experimental evidence. Moreover, we prove that the positive tropical Grassmannian is the secondary fan for the regular positroid subdivisions of the hypersimplex, and propose that it also controls the T-dual positroid subdivisions of the amplituhedron. Along the way, we prove that a matroid polytope is a positroid polytope if and only if all two-dimensional faces are positroid polytopes. Towards the goal of generalizing T-duality for higher $m$, we also define the momentum amplituhedron for any even $m$.


title: Unsupervised Speaker Adaptation using Attention-based Speaker Memory for End-to-End ASR

abstract: We propose an unsupervised speaker adaptation method inspired by the neural Turing machine for end-to-end (E2E) automatic speech recognition (ASR). The proposed model contains a memory block that holds speaker i-vectors extracted from the training data and reads relevant i-vectors from the memory through an attention mechanism. The resulting memory vector (M-vector) is concatenated to the acoustic features or to the hidden layer activations of an E2E neural network model. The E2E ASR system is based on the joint connectionist temporal classification and attention-based encoder-decoder architecture. M-vector and i-vector results are compared for inserting them at different layers of the encoder neural network using the WSJ and TED-LIUM2 ASR benchmarks. We show that M-vectors, which do not require an auxiliary speaker embedding extraction system at test time, achieve similar word error rates (WERs) compared to i-vectors for single speaker utterances and significantly lower WERs for utterances in which there are speaker changes.


title: Lower bounds on Hilbert--Kunz multiplicities and maximal F-signatures

abstract: Hilbert-Kunz multiplicity and F-signature are numerical invariants of commutative rings in positive characteristic that measure severity of singularities: for a regular ring both invariants are equal to one and the converse holds under mild assumptions. A natural question is for what singular rings these invariants are closest to one. For Hilbert--Kunz multiplicity this question was first considered by the last two authors and attracted significant attention. In this paper, we study this question, i.e., an upper bound, for F-signature and revisit lower bounds on Hilbert--Kunz multiplicity.


title: Controlling magnetoresistance by tuning semimetallicity through dimensional confinement and heteroepitaxy

abstract: Controlling the electronic properties via bandstructure engineering is at the heart of modern semiconductor devices. Here, we extend this concept to semimetals where, utilizing LuSb as a model system, we show that quantum confinement lifts carrier compensation and differentially affects the mobility of the electron and hole-like carriers resulting in a strong modification in its large, non-saturating magnetoresistance behavior. Bonding mismatch at the heteroepitaxial interface of a semimetal (LuSb) and a semiconductor (GaSb) leads to the emergence of a novel, two-dimensional, interfacial hole gas and is accompanied by a charge transfer across the interface that provides another avenue to modify the electronic structure and magnetotransport properties in the ultra-thin limit. Our work lays out a general strategy of utilizing confined thin film geometries and heteroepitaxial interfaces to engineer electronic structure in semimetallic systems, which allows control over their magnetoresistance behavior and simultaneously, provides insights into its origin.


title: Spectral Extension and Synchronisation of Microcombs in a Single Microresonator

abstract: We demonstrate that the optical spectrum of a soliton microcomb generated in a microresonator can be coherently extended by using bichromatic pumping: one pump with a frequency in the anomalous dispersion regime of the microresonator is used to generate a bright soliton microcomb while a second pump in the normal dispersion regime both compensates the thermal effect of the microresonator and generates a second frequency comb via cross-phase modulation with the bright soliton, thus extending the spectrum of the microcomb. Experimental results for different microresonators show that the mode spacings of the two combs are synchronised, even when the combs are generated in different mode families. Numerical simulations agree well with experimental results and reveal that a bright optical pulse from the second pump is passively formed in the normal dispersion regime and trapped by the bright soliton. We believe the demonstrated technique provides an alternative way to generate broadband microcombs and multi-coloured frequency combs that can be useful for frequency metrology and spectroscopy. In addition, this method enables the selective enhancement of available optical power in specific parts of the comb spectrum.


title: Learning models of quantum systems from experiments

abstract: An isolated system of interacting quantum particles is described by a Hamiltonian operator. Hamiltonian models underpin the study and analysis of physical and chemical processes throughout science and industry, so it is crucial they are faithful to the system they represent. However, formulating and testing Hamiltonian models of quantum systems from experimental data is difficult because it is impossible to directly observe which interactions the quantum system is subject to. Here, we propose and demonstrate an approach to retrieving a Hamiltonian model from experiments, using unsupervised machine learning. We test our methods experimentally on an electron spin in a nitrogen-vacancy interacting with its spin bath environment, and numerically, finding success rates up to 86%. By building agents capable of learning science, which recover meaningful representations, we can gain further insight on the physics of quantum systems.


title: Transformer on a Diet

abstract: Transformer has been widely used thanks to its ability to capture sequence information in an efficient way. However, recent developments, such as BERT and GPT-2, deliver only heavy architectures with a focus on effectiveness. In this paper, we explore three carefully-designed light Transformer architectures to figure out whether the Transformer with less computations could produce competitive results. Experimental results on language model benchmark datasets hint that such trade-off is promising, and the light Transformer reduces 70% parameters at best, while obtains competitive perplexity compared to standard Transformer. The source code is publicly available.


title: Providing a Link Prediction Model based on Structural and Homophily Similarity in Social Networks

abstract: In recent years, with the growing number of online social networks, these networks have become one of the best markets for advertising and commerce, so studying these networks is very important. Forecasting new edges in online social networks can give us a better understanding of the growth of these networks. There have been many studies of link prediction in the field of engineering and humanities. Scientists attribute the existence of a new relationship between two individuals for two reasons: 1) Proximity to the graph (structure) 2) Similar properties of the two individuals (Homophile law). However, studying the impact of the two approaches working together to create new edges remains an open problem. Similarity metrics can also be divided into two categories; Neighborhood-based and path-based. So far, the above two theoretical approaches (proximity and homophile) have not been found together in the neighborhood-based metrics. In this paper, we first attempt to provide a solution to determine importance of the proximity to the graph and similar features in the connectivity of the graphs. Then obtained weights are assigned to both proximity and homophile. Then the best similarity metric in each approach are obtained. Finally, the selected metric of homophily similarity and structural similarity are combined with the obtained weights. The results of this study were evaluated on two datasets; Zanjan University Graduate School of Social Sciences and Pokec online Social Network. The first data set was collected for this study and then the questionnaires and data collection methods were filled out. Since this dataset is one of the few Iranian datasets that has been compiled with its users' specifications, it can be of great value. In this paper, we have been able to increase the accuracy of Neighborhood-based similarity metric by using two proximity in graph and homophily approaches.


title: Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

abstract: We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.


title: Effects of intrinsic defects and alloying with Fe on the half-metallicity of Co$_2$MnSi

abstract: The electronic structure and half-metallic gap of Co$_{2}$MnSi in the presence of crystallographic defects, partial Fe substitution for Mn, and thermal spin fluctuations are studied using the coherent potential approximation and the disordered local moment method. In the presence of 5\% Co or Mn vacancies the Fermi level shifts down to the minority-spin valence-band maximum. In contrast to NiMnSb, both types of Mn antisite defects in Co$_{2}$MnSi are strongly exchange-coupled to the host magnetization, and thermal spin fluctuations do not strongly affect the half-metallic gap. Partial substitution of Mn by Fe results in considerable changes in the Bloch spectral function near the Fermi level, which strongly deviate from the rigid-band picture. In particular, a light band with the Fe character crosses the Fermi level at about 50\% concentration. At room temperature, Fe substitution of up to 30\% slightly increases the spin polarization at the Fermi level.


title: Classical critical dynamics in quadratically driven Kerr resonators

abstract: Driven-dissipative kerr lattices with two-photon driving are experimentally relevant systems known to exhibit a symmetry-breaking phase transition, which belongs to the universality class of the thermal Ising model for the parameter regime studied here. In this work, we perform finite-size scaling of this system as it is quenched to the transition and the dynamical critical exponent is found to be compatible with $z\approx2.18$ corresponding with metropolis dynamics in classical simulations. Furthermore, we show that the Liouvillian gap scales with the same exponent, similar to scaling of the Hamiltonian gap at quantum phase transitions in closed systems.


title: Improving accuracy of the fifth-order WENO scheme by using the exponential approximation space

abstract: The aim of this study is to develop a novel WENO scheme that improves the performance of the well-known fifth-order WENO methods. The approximation space consists of exponential polynomials with a tension parameter that may be optimized to fit the the specific feature of the data, yielding better results compared to the polynomial approximation space. However, finding an optimal tension parameter is a very important and difficult problem, indeed a topic of active research. In this regard, this study introduces a practical approach to determine an optimal tension parameter by taking into account the relationship between the tension parameter and the accuracy of the exponential polynomial interpolation under the setting of the fifth-order WENO scheme. As a result, the proposed WENO scheme attains an improved order of accuracy (that is, sixth-order) better than other fifth-order WENO methods without loss of accuracy at critical points. A detailed analysis is provided to verify the improved convergence rate. Further, we present modified nonlinear weights based on $L^1$-norm approach along with a new global smoothness indicator. The proposed nonlinear weights reduce numerical dissipation significantly, while attaining better resolution in smooth regions. Some experimental results for various benchmark test problems are presented to demonstrate the ability of the new scheme.


title: Functional Programming in Pattern-Match-Oriented Programming Style

abstract: Throughout the history of functional programming, recursion has emerged as a natural method for describing loops in programs. However, there does often exist a substantial cognitive distance between the recursive definition and the simplest explanation of an algorithm even for the basic list processing functions such as map, concat, or unique; when we explain these functions, we seldom use recursion explicitly as we do in functional programming. For example, map is often explained as follows: the map function takes a function and a list and returns a list of the results of applying the function to all the elements of the list. This paper advocates a new programming paradigm called pattern-match-oriented programming for filling this gap. An essential ingredient of our method is utilizing pattern matching for non-free data types. Pattern matching for non-free data types features non-linear pattern matching with backtracking and extensibility of pattern-matching algorithms. Several non-standard pattern constructs, such as not-patterns, loop patterns, and sequential patterns, are derived from this pattern-matching facility. Based on that result, this paper introduces many programming techniques that replace explicit recursions with an intuitive pattern by confining recursions inside patterns. We classify these techniques as pattern-match-oriented programming design patterns. These programming techniques allow us to redefine not only the most basic functions for list processing such as map, concat, or unique more elegantly than the traditional functional programming style, but also more practical mathematical algorithms and software such as a SAT solver, computer algebra system, and database query language that we had not been able to implement concisely.


title: The Next Decade in AI: Four Steps Towards Robust Artificial Intelligence

abstract: Recent research in artificial intelligence and machine learning has largely emphasized general-purpose learning and ever-larger training sets and more and more compute. In contrast, I propose a hybrid, knowledge-driven, reasoning-based approach, centered around cognitive models, that could provide the substrate for a richer, more robust AI than is currently possible.