title: Practical aspects of X-ray imaging polarimetry of supernova remnants and other extended sources

abstract: The new generation of X-ray polarisation detectors, gas pixel detectors, which will be employed by the future space missions IXPE and eXTP, allows for spatially resolved X-ray polarisation studies. This will be of particular interest for X-ray synchrotron emission from extended sources like young supernova remnants and pulsar wind nebulae. Here we report on employing a polarisation statistic that can be used to makes maps in the Stokes I, Q, and U parameters, a method that we expand by correcting for the energy-dependent instrumental modulation factor, using optimal weighting of the signal. In order to explore the types of Stokes maps that can be obtained, we present a Monte Carlo simulation program called xpolim, with which different polarisation weighting schemes are explored. We illustrate its use with simulations of polarisation maps of young supernova remnants, after having described the general science case for polarisation studies of supernova remnants, and its connection to magnetic-field turbulence. We use xpolim simulations to show that in general deep, ~2 Ms observations are needed to recover polarisation signals, in particular for Cas A, for which in the polarisation fraction may be as low as 5%.


title: Extended Homotopy Quantum Field Theories and their Orbifoldization

abstract: We present a precise definition of extended homotopy quantum field theories and develop an orbifold construction for these theories when the target space is the classifying space of a finite group $G$, i.e. for $G$-equivariant topological field theories. More precisely, we use a bicategorical version of the parallel section functor to associate to an extended equivariant topological field theory an ordinary extended topological field theory. Thereby, we give a unification, geometric underpinning and vast generalization of algebraic concepts of orbifoldization. In the special case of 3-2-1-dimensional equivariant topological field theories, we investigate the equivariant modular structure on the categories that such theories yield upon evaluation on the circle. By means of our orbifold construction this equivariant modular structure will be related to the modular structure on the orbifold category. We also generalize our orbifold construction to a pushforward operation along an arbitrary morphism of finite groups and hence provide a valuable tool for the construction of extended homotopy quantum field theories.


title: Fast and Sample Near-Optimal Algorithms for Learning Multidimensional Histograms

abstract: We study the problem of robustly learning multi-dimensional histograms. A $d$-dimensional function $h: D \rightarrow \mathbb{R}$ is called a $k$-histogram if there exists a partition of the domain $D \subseteq \mathbb{R}^d$ into $k$ axis-aligned rectangles such that $h$ is constant within each such rectangle. Let $f: D \rightarrow \mathbb{R}$ be a $d$-dimensional probability density function and suppose that $f$ is $\mathrm{OPT}$-close, in $L_1$-distance, to an unknown $k$-histogram (with unknown partition). Our goal is to output a hypothesis that is $O(\mathrm{OPT}) + \epsilon$ close to $f$, in $L_1$-distance. We give an algorithm for this learning problem that uses $n = \tilde{O}_d(k/\epsilon^2)$ samples and runs in time $\tilde{O}_d(n)$. For any fixed dimension, our algorithm has optimal sample complexity, up to logarithmic factors, and runs in near-linear time. Prior to our work, the time complexity of the $d=1$ case was well-understood, but significant gaps in our understanding remained even for $d=2$.


title: Silicon Vibrating Wires at Low Temperatures

abstract: Nowadays microfabrication techniques originating from micro-electronics enable to create mechanical objects of micron-size. The field of Micro-Electro-Mechanical devices (MEMS) is continuously expanding, with an amazingly broad range of applications at room temperature. Vibrating objects (torsional oscillators, vibrating wires) widely used at low temperatures to study quantum fluids, can be replaced advantageously by Silicon MEMS. In this letter we report on the study of Silicon vibrating wire devices. A goal-post structure covered with a metal layer is driven at resonance by the Laplace force acting on a current in a magnetic field, while the induced voltage arising from the cut magnetic flux allows to detect the motion. The characteristics of the resonance have been studied from 10 mK to 30 K, in vacuum and in $^4$He gas. In this article, we focus on the results obtained above 1.5 K, in vacuum and gas, and introduce some features observed at lower temperatures. The resonant properties can be quantitatively understood by means of simple models, from the linear regime to a highly non-linear response at strong drives. We demonstrate that the non-linearity is mostly due to the geometry of the vibrators. We also show that in our device the friction mechanisms originate in the metallic layers, and can be fully characterized. The interaction with $^4$He gas is fit to theory without adjustable parameters.


title: Closed-form solution to cooperative visual-inertial structure from motion

abstract: This paper considers the problem of visual-inertial sensor fusion in the cooperative case and it provides new theoretical contributions, which regard its observability and its resolvability in closed form. The case of two agents is investigated. Each agent is equipped with inertial sensors (accelerometer and gyroscope) and with a monocular camera. By using the monocular camera, each agent can observe the other agent. No additional camera observations (e.g., of external point features in the environment) are considered. All the inertial sensors are assumed to be affected by a bias. First, the entire observable state is analytically derived. This state includes the absolute scale, the relative velocity between the two agents, the three Euler angles that express the rotation between the two agent frames and all the accelerometer and gyroscope biases. Second, the paper provides the extension of the closed-form solution given in [19] (which holds for a single agent) to the aforementioned cooperative case. The impact of the presence of the bias on the performance of this closed-form solution is investigated. As in the case of a single agent, this performance is significantly sensitive to the presence of a bias on the gyroscope, while, the presence of a bias on the accelerometer is negligible. Finally, a simple and effective method to obtain the gyroscope bias is proposed. Extensive simulations clearly show that the proposed method is successful. It is amazing that, it is possible to automatically retrieve the absolute scale and simultaneously calibrate the gyroscopes not only without any prior knowledge (as in [13]), but also without external point features in the environment.


title: 6D Pose Estimation using an Improved Method based on Point Pair Features

abstract: The Point Pair Feature (Drost et al. 2010) has been one of the most successful 6D pose estimation method among model-based approaches as an efficient, integrated and compromise alternative to the traditional local and global pipelines. During the last years, several variations of the algorithm have been proposed. Among these extensions, the solution introduced by Hinterstoisser et al. (2016) is a major contribution. This work presents a variation of this PPF method applied to the SIXD Challenge datasets presented at the 3rd International Workshop on Recovering 6D Object Pose held at the ICCV 2017. We report an average recall of 0.77 for all datasets and overall recall of 0.82, 0.67, 0.85, 0.37, 0.97 and 0.96 for hinterstoisser, tless, tudlight, rutgers, tejani and doumanoglou datasets, respectively.


title: Exclusive vector meson photoproduction with proton dissociation in photon-hadron interactions at the LHC

abstract: At forward rapidities and high energies we expect to probe the non-linear regime of Quantum Chromodynamics (QCD). One of the most promising observables to constrain the QCD dynamics in this regime is exclusive vector meson photoproduction (EVMP). We study the EVMP in association with a leading baryon (product of the proton dissociation) in photon-hadron interactions that take place in $pp$ and $p$Pb collisions at large impact parameters. We present the rapidity distributions for $\rho$ and $J/\psi$ photoproduction in association with a leading baryon (neutron and delta states) at LHC run 2 energies. Our results show that the $V + \Delta$ cross section is almost 30 % of the $V + n$ one. Our results also show that a future experimental analysis of these processes is, in principle, feasible and can be useful to study the leading particle production.


title: Hardy spaces for the Dunkl harmonic oscillator

abstract: Let $\Delta$ and $L=\Delta -\|\mathbf x\|^2$ be the Dunkl Laplacian and the Dunkl harmonic oscillator respectively. We define the Hardy space $\mathcal H^1$ associated with the Dunkl harmonic oscillator by means of the nontangential maximal function with respect to the semigroup $e^{tL}$. We prove that the space $\mathcal H^1$ admits characterizations by relevant Riesz transforms and atomic decompositions. The atoms which occur in the atomic decompositions are of local type.


title: Simple advecting structures and the edge of chaos in subcritical tokamak plasmas

abstract: In tokamak plasmas, sheared flows perpendicular to the driving temperature gradients can strongly stabilize linear modes. While the system is linearly stable, regimes with persistent nonlinear turbulence may develop, i.e. the system is subcritical. A perturbation with small but finite amplitude may be sufficient to push the plasma into a regime where nonlinear effects are dominant and thus allow sustained turbulence. The minimum threshold for nonlinear instability to be triggered provides a criterion for assessing whether a tokamak is likely to stay in the quiescent (laminar) regime. At the critical amplitude, instead of transitioning to the turbulent regime or decaying to a laminar state, the trajectory will map out the edge of chaos. Surprisingly, a quasi-traveling-wave solution is found as an attractor on this edge manifold. This simple advecting solution is qualitatively similar to, but simpler than, the avalanche-like bursts seen in earlier turbulent simulations and provides an insight into how turbulence is sustained in subcritical plasma systems. For large flow shearing rate, the system is only convectively unstable, and given a localised initial perturbation, will eventually return to a laminar state at a fixed spatial location.


title: Non-Uniqueness of Stationary Solutions in Extremum Seeking Control

abstract: Extremum seeking control (ESC) is a classical adaptive control method for steady-state optimization, purely based on output feedback. It is well known that the extremum seeking control loop, under certain mild conditions on the controller, has a stable stationary periodic solution in the vicinity of an extremum point of the steady-state input-output map of the plant. However, this is a local result only and this paper investigates whether this solution is necessarily unique given that the underlying optimization problem is convex. We first derive a necessary condition that any stationary solution of the ESC loop must satisfy. For plants in which the extremum point is due to a purely static nonlinearity, such as in Hammerstein or Wiener plants, the condition involves the steady-state gradient. However, for more general plants the necessary condition involves the phase lag of the locally linearized plant, indicating the possible existence of solutions without any relationship to optimality. Combining the derived necessary condition with the existence of a local solution close to the optimum, we employ bifurcation theory to trace out branches of stationary solutions. We derive conditions on when these branches bifurcate, resulting in multiple stationary solutions. The results show that cyclic fold bifurcations may exist, resulting in the existence of multiple stationary period one solutions, of which only one is related to the optimality conditions. We illustrate the results through an example in which the conversion in a chemical reactor is optimized using ESC. We show that at least five stationary solutions may exist simultaneously for realistic control parameters, and that several of these solutions are stable. One consequence of the non-uniqueness is that one in general needs to start close to the optimum to ensure convergence to the near-optimal solution.


title: Topological and metric spaces are full subcategories of the category of simplicial objects of the category of filters

abstract: We observe that the category of topological space, uniform spaces, and simplicial sets are all, in a natural way, full subcategories of the same larger category, namely the simplicial category of filters; this is, moreover, implicit in the definitions of a topological and uniform space. We use these embeddings to rewrite the notions of completeness, precompactness, compactness, Cauchy sequence, and equicontinuity in the language of category theory, which we hope might be of use in formalisation of mathematics and tame topology. We formulate some arising open questions.


title: SimCommSys: Taking the errors out of error-correcting code simulations

abstract: In this paper we present SimCommSys, a Simulator of Communication Systems that we are releasing under an open source license. The core of the project is a set of C++ libraries defining communication system components and a distributed Monte Carlo simulator. Of principal interest is the error-control coding component, where various kinds of binary and non-binary codes are implemented, including turbo, LDPC, repeat-accumulate, and Reed-Solomon. The project also contains a number of ready-to-build binaries implementing various stages of the communication system (such as the encoder and decoder), a complete simulator, and a system benchmark. Finally, SimCommSys also provides a number of shell and python scripts to encapsulate routine use cases. As long as the required components are already available in SimCommSys, the user may simulate complete communication systems of their own design without any additional programming. The strict separation of development (needed only to implement new components) and use (to simulate specific constructions) encourages reproducibility of experimental work and reduces the likelihood of error. Following an overview of the framework, we provide some examples of how to use the framework, including the implementation of a simple codec, the specification of communication systems and their simulation.


title: Orienting lipid-coated graphitic micro-particles in solution using AC electric fields: A new theoretical dual-ellipsoid Laplace model for electro-orientation

abstract: Graphitic micro-particles coated with thin layers in solution are technologically interesting as they can be manipulated with electric fields. Modeling the electrical manipulation of submerged layered micro-particles analytically or numerically is not straight forward. In particular, the generation of reliable quantitative torque predictions for electro-orientation experiments has been elusive. The traditional Laplace model approximates the coated particle by an ellipsoid with a confocal ellipsoidal layer and solves Laplace's equation to produce convenient analytical predictions. However, due to the non-uniformity of the layer thickness around the ellipsoid, this method can lead to incorrect torque predictions. Here we present a new theoretical dual-ellipsoid Laplace model that corrects the effect of the non-uniform layer thickness by calculating two layered ellipsoids, each accounting for the correct layer thickness along each relevant direction for the torque. Our model describes the electro-orientation of submerged lipid-coated graphitic micro-particles in the presence of an alternating current (AC) electric field and is valid for ellipsoids with moderate aspect ratios and coated with thin shells. It is one of the first models to generate correct quantitative electric torque predictions. We present model results for the torque versus frequency and compare them to our measurements for lipid-coated highly ordered pyrolytic graphite (HOPG) micro-flakes in aqueous NaCl solution at MHz frequencies. The results show how the lipid shell changes the overall electrical properties of the micro-flakes so that the torque is low at low frequencies and increases at higher frequencies into the MHz regime. The torque depends critically on the lipid-shell thickness, the solution conductivity and the shape of the particle, all of which can be used as handles to control the response of the particles.


title: The Cram\'er conjecture holds with a positive probability

abstract: We prove that a positive proportion of the intervals of any fixed scalar multiple of $\log(X)$ in the dyadic interval $[X,2X]$ contain a prime number. We also show that a positive proportion of the congruence classes modulo $q$ contain a prime number smaller than any fixed scalar multiple of $\varphi(q)\log(q).$


title: Wormhole Potentials and Throats from Quasi-Normal Modes

abstract: Exotic compact objects refer to a wide class of black hole alternatives or effective models to describe phenomenologically quantum gravitational effects on the horizon scale. In this work we show how the knowledge of the quasi-normal mode spectrum of non-rotating wormhole models can be used to reconstruct the effective potential that appears in the perturbation equations. From this it is further possible to obtain the parameters that characterize the specific wormhole model, which in this paper was chosen to be the one by Damour and Solodukhin. We also address the question whether one can distinguish such type of wormholes from ultra compact stars, if only the quasi-normal mode spectrum is known. We have proven that this is not possible by using the trapped modes only, but requires additional information. The here presented inverse method is an extension of work that has previously been developed and applied to the oscillation spectra of ultra compact stars and gravastars. However, it is not limited to the study of exotic compact objects, but applicable to symmetric double barrier potentials that appear in one-dimensional wave equations. Therefore we think it can be of interest for other fields too.


title: The Weighted Kendall and High-order Kernels for Permutations

abstract: We propose new positive definite kernels for permutations. First we introduce a weighted version of the Kendall kernel, which allows to weight unequally the contributions of different item pairs in the permutations depending on their ranks. Like the Kendall kernel, we show that the weighted version is invariant to relabeling of items and can be computed efficiently in $O(n \ln(n))$ operations, where $n$ is the number of items in the permutation. Second, we propose a supervised approach to learn the weights by jointly optimizing them with the function estimated by a kernel machine. Third, while the Kendall kernel considers pairwise comparison between items, we extend it by considering higher-order comparisons among tuples of items and show that the supervised approach of learning the weights can be systematically generalized to higher-order permutation kernels.


title: Reductions of points on algebraic groups, II

abstract: Let $A$ be the product of an abelian variety and a torus over a number field $K$, and let $m$ be a positive integer. If $\alpha \in A(K)$ is a point of infinite order, we consider the set of primes $\mathfrak p$ of $K$ such that the reduction $(\alpha \bmod \mathfrak p)$ is well defined and has order coprime to $m$. This set admits a natural density, which we are able to express as a finite sum of products of $\ell$-adic integrals, where $\ell$ varies in the set of prime divisors of $m$. We deduce that the density is a rational number, whose denominator is bounded (up to powers of $m$) in a very strong sense. This extends the results of the paper "Reductions of points on algebraic groups" by Davide Lombardo and the second author, where the case $m$ prime is established.


title: Sinking during earthquakes: Critical acceleration criteria control drained soil liquefaction

abstract: This article focuses on liquefaction of saturated granular soils, triggered by earthquakes. Liquefaction is definedhere as the transition from a rigid state, in which the granular soil layer supports structures placed on its surface, toa fluidlike state, in which structures placed initially on the surface sink to their isostatic depth within the granularlayer.We suggest a simple theoretical model for soil liquefaction and show that buoyancy caused by the presence ofwater inside a granular medium has a dramatic influence on the stability of an intruder resting at the surface of themedium.We confirm this hypothesis by comparison with laboratory experiments and discrete-element numericalsimulations. The external excitation representing ground motion during earthquakes is simulated via horizontalsinusoidal oscillations of controlled frequency and amplitude. In the experiments, we use particles only slightlydenser than water, which as predicted theoretically increases the effect of liquefaction and allows clear depth-of-sinkingmeasurements. In the simulations, a micromechanical model simulates grains using molecular dynamicswith friction between neighbors. The effect of the fluid is captured by taking into account buoyancy effects onthe grains when they are immersed. We show that the motion of an intruder inside a granular medium is mainlydependent on the peak acceleration of the ground motion and establish a phase diagram for the conditions underwhich liquefaction happens, depending on the soil bulk density, friction properties, presence of water, and peak acceleration of the imposed large-scale soil vibrations.We establish that in liquefaction conditions, most cases relaxtoward an equilibrium position following an exponential in time.We also show that the equilibrium position itself,for most liquefaction regimes, corresponds to the isostatic equilibrium of the intruder inside a medium of effectivedensity. The characteristic time to relaxation is shown to be essentially a function of the peak ground velocity.


title: A $q$-analogue for Euler's $\zeta(6)=\pi^6/945$

abstract: We give a $q$-analogue of $\zeta(6)=\pi^6/945$. Our main results are stated in Theorems 2.1 and 2.2 below.


title: Training wide residual networks for deployment using a single bit for each weight

abstract: For fast and energy-efficient deployment of trained deep neural networks on resource-constrained embedded hardware, each learned weight parameter should ideally be represented and stored using a single bit. Error-rates usually increase when this requirement is imposed. Here, we report large improvements in error rates on multiple datasets, for deep convolutional neural networks deployed with 1-bit-per-weight. Using wide residual networks as our main baseline, our approach simplifies existing methods that binarize weights by applying the sign function in training; we apply scaling factors for each layer with constant unlearned values equal to the layer-specific standard deviations used for initialization. For CIFAR-10, CIFAR-100 and ImageNet, and models with 1-bit-per-weight requiring less than 10 MB of parameter memory, we achieve error rates of 3.9%, 18.5% and 26.0% / 8.5% (Top-1 / Top-5) respectively. We also considered MNIST, SVHN and ImageNet32, achieving 1-bit-per-weight test results of 0.27%, 1.9%, and 41.3% / 19.1% respectively. For CIFAR, our error rates halve previously reported values, and are within about 1% of our error-rates for the same network with full-precision weights. For networks that overfit, we also show significant improvements in error rate by not learning batch normalization scale and offset parameters. This applies to both full precision and 1-bit-per-weight networks. Using a warm-restart learning-rate schedule, we found that training for 1-bit-per-weight is just as fast as full-precision networks, with better accuracy than standard schedules, and achieved about 98%-99% of peak performance in just 62 training epochs for CIFAR-10/100. For full training code and trained models in MATLAB, Keras and PyTorch see https://github.com/McDonnell-Lab/1-bit-per-weight/ .


title: Soft gluon evolution and non-global logarithms

abstract: We consider soft-gluon evolution at the amplitude level. Our evolution includes Coulomb exchanges and applies to generic hard-scattering processes involving any number of coloured partons. We emphasise the special role played by a Lorentz-invariant evolution variable, which coincides with the transverse momentum of the latest emission in a suitably defined dipole zero-momentum frame. We also relate the evolution algorithm, which was used originally in the derivation of super-leading logarithms, to renormalization group evolution equations that have been encountered recently. Handling large colour matrices presents the most significant challenge to numerical implementations and we present a means to expand systematically about the leading colour approximation.


title: ZpL: a p-adic precision package

abstract: We present a new package ZpL for the mathematical software system SM. It implements a sharp tracking of precision on p-adic numbers, following the theory of ultrametric precision introduced in [4]. The underlying algorithms are mostly based on automatic dierentiation techniques. We introduce them, study their complexity and discuss our design choices. We illustrate the bene-ts of our package (in comparison with previous implementations) with a large sample of examples coming from linear algebra, com-mutative algebra and dierential equations.


title: On Emergent Gravity, Black Hole Entropy and Galactic Rotation Curves

abstract: In this work we derive a generalized Newtonian gravitational force and show that it can account for the anomalous galactic rotation curves. We derive the entropy-area relationship applying the Feynman-Hibbs procedure to the supersymmetric Wheeler-DeWitt equation of the Schwarzschild black hole. We obtain the modifications to the Newtonian gravitational force from the entropic formulation of gravity.


title: Weighted Double Deep Multiagent Reinforcement Learning in Stochastic Cooperative Environments

abstract: Recently, multiagent deep reinforcement learning (DRL) has received increasingly wide attention. Existing multiagent DRL algorithms are inefficient when facing with the non-stationarity due to agents update their policies simultaneously in stochastic cooperative environments. This paper extends the recently proposed weighted double estimator to the multiagent domain and propose a multiagent DRL framework, named weighted double deep Q-network (WDDQN). By utilizing the weighted double estimator and the deep neural network, WDDQN can not only reduce the bias effectively but also be extended to scenarios with raw visual inputs. To achieve efficient cooperation in the multiagent domain, we introduce the lenient reward network and the scheduled replay strategy. Experiments show that the WDDQN outperforms the existing DRL and multiaent DRL algorithms, i.e., double DQN and lenient Q-learning, in terms of the average reward and the convergence rate in stochastic cooperative environments.


title: Can Neural Networks Understand Logical Entailment?

abstract: We introduce a new dataset of logical entailments for the purpose of measuring models' ability to capture and exploit the structure of logical expressions against an entailment prediction task. We use this task to compare a series of architectures which are ubiquitous in the sequence-processing literature, in addition to a new model class---PossibleWorldNets---which computes entailment as a "convolution over possible worlds". Results show that convolutional networks present the wrong inductive bias for this class of problems relative to LSTM RNNs, tree-structured neural networks outperform LSTM RNNs due to their enhanced ability to exploit the syntax of logic, and PossibleWorldNets outperform all benchmarks.


title: Continuous-time Markov games with asymmetric information

abstract: We study a two-player zero-sum stochastic differential game with asymmetric information where the payoff depends on a controlled continuous-time Markov chain X with finite state space which is only observed by player 1. This model was already studied in Cardaliaguet et al. (2015) through an approximating sequence of discrete-time games. Our first contribution is the proof of the existence of the value in the continuous-time model based on duality techniques. This value is shown to be the unique solution of the same Hamilton-Jacobi equation with convexity constraints which characterized the limit value obtained in Cardaliaguet et al. (2015). Our second main contribution is to provide a simpler equivalent formulation for this Hamilton-Jacobi equation using directional derivatives and exposed points, which we think is interesting for its own sake as the associated comparison principle has a very simple proof which avoids all the technical machinery of viscosity solutions.


title: On the images of Galois representations attached to low weight Siegel modular forms

abstract: Let $\pi$ be a cuspidal automorphic representation of $\mathrm{GSp}_4(\mathbf{A_Q})$, whose archimedean component is a holomorphic discrete series or limit of discrete series representation. If $\pi$ is not CAP or endoscopic, then we show that its associated $\ell$-adic Galois representations are irreducible and crystalline for $100\%$ of primes $\ell$. If, moreover, $\pi$ is neither an automorphic induction nor a symmetric cube lift, then we show that, for $100\%$ of primes $\ell$, the image of its mod $\ell$ Galois representation contains $\mathrm{Sp}_4(\mathbf{F}_{\ell})$.


title: Thermal algebraic-decay charge liquid driven by competing short-range Coulomb repulsion

abstract: We explore the possibility of a Berezinskii-Kosterlitz-Thouless-like critical phase for the charge degrees of freedom in the intermediate-temperature regime between the charge-ordered and disordered phases in two-dimensional systems with competing short-range Coulomb repulsion. As the simplest example, we investigate the extended Hubbard model with on-site and nearest-neighbor Coulomb interactions on a triangular lattice at half filling in the atomic limit by using a classical Monte Carlo method, and find a critical phase, characterized by algebraic decay of the charge correlation function, belonging to the universality class of the two-dimensional XY model with a $\mathbb{Z}_6$ anisotropy. Based on the results, we discuss possible conditions for the critical phase in materials.


title: Computation of optimal transport and related hedging problems via penalization and neural networks

abstract: This paper presents a widely applicable approach to solving (multi-marginal, martingale) optimal transport and related problems via neural networks. The core idea is to penalize the optimization problem in its dual formulation and reduce it to a finite dimensional one which corresponds to optimizing a neural network with smooth objective function. We present numerical examples from optimal transport, martingale optimal transport, portfolio optimization under uncertainty and generative adversarial networks that showcase the generality and effectiveness of the approach.


title: Optimal Stochastic Delivery Planning in Full-Truckload and Less-Than-Truckload Delivery

abstract: With an increasing demand from emerging logistics businesses, Vehicle Routing Problem with Private fleet and common Carrier (VRPPC) has been introduced to manage package delivery services from a supplier to customers. However, almost all of existing studies focus on the deterministic problem that assumes all parameters are known perfectly at the time when the planning and routing decisions are made. In reality, some parameters are random and unknown. Therefore, in this paper, we consider VRPPC with hard time windows and random demand, called Optimal Delivery Planning (ODP). The proposed ODP aims to minimize the total package delivery cost while meeting the customer time window constraints. We use stochastic integer programming to formulate the optimization problem incorporating the customer demand uncertainty. Moreover, we evaluate the performance of the ODP using test data from benchmark dataset and from actual Singapore road map.


title: Chaos, decoherence and emergent extradimensions in D-brane dynamics with fluctuations

abstract: We study, by using tools of the dynamical system theory, a fermionic string streched from a non-commutative D2-brane (stack of D0-branes in the BFSS model) to a probe D0-brane as a quantum system driven by a chaotic system, the (classical and quantum) chaos being induced in the D2-brane dynamics by quantum fluctuations. We show that this dynamics with fluctuations induces a decoherence phenomenon on the reduced density matrix of the fermionic string which is characteristic of the chaotic behaviour since it presents an horizon of coherence. Moreover we show that this dynamics is associated with an invariant torus which involves extradimensions emerging from the fluctuations for the viewpoint of the fermionic string, extending a three dimensional space by six compact dimensions. The situation studied can be considered as a model of qubit (supported by the fermionic string) in interaction with a quantum black hole (modelized by the non-commutative D2-brane).


title: Nanocircuits in loop structures: continuous waves preclude gauge invariant wavelengths

abstract: Tunnel junctions for quantum computing require discrete spectra from continuous waves on a doubly connected coordinate or loop. For an electron on a metal ring discrete spectra follow from discontinuous Bloch waves. Can both propositions be true? We find using a gauge function originating in the Lagrangian that continuity on a ring or loop violates gauge invariance of the de Broglie wavelength. This same gauge function shows that Lagrangians for the electron on a ring and the charge on a junction are mutual transforms. Thus persistent current on a metal ring and the Coulomb blockade on a tunnel junction seem to be the same dynamical theory based on discontinuous Bloch waves on the compact perimeter of a circle


title: Large-$N$ $\mathbb{CP}^{N-1}$ sigma model on a finite interval: general Dirichlet boundary conditions

abstract: This is the third of the series of articles on the large-$N$ two-dimensional $\mathbb{CP}^{N-1}$ sigma model, defined on a finite space interval $L$ with Dirichlet boundary conditions. Here the cases of the general Dirichlet boundary conditions are studied, where the relative $\mathbb{CP}^{N-1}$ orientations at the two boundaries are generic, and numerical solutions are presented. Distinctive features of the $\mathbb{CP}^{N-1}$ sigma model, as compared e.g., to an $O(N)$ model, which were not entirely evident in the basic properties studied in the first two articles in the large $N$ limit, manifest themselves here. It is found that the total energy is minimized when the fields are aligned in the same direction at the two boundaries.


title: Geometrical Equivalence and Action Type Geometrical Equivalence of Group Representations

abstract: In this paper we present the example which proves that we can not conclude the geometrical equivalence of group representations from the corresponding action-type geometrical equivalence and group geometrical equivalence.


title: Unsupervised Grammar Induction with Depth-bounded PCFG

abstract: There has been recent interest in applying cognitively or empirically motivated bounds on recursion depth to limit the search space of grammar induction models (Ponvert et al., 2011; Noji and Johnson, 2016; Shain et al., 2016). This work extends this depth-bounding approach to probabilistic context-free grammar induction (DB-PCFG), which has a smaller parameter space than hierarchical sequence models, and therefore more fully exploits the space reductions of depth-bounding. Results for this model on grammar acquisition from transcribed child-directed speech and newswire text exceed or are competitive with those of other models when evaluated on parse accuracy. Moreover, gram- mars acquired from this model demonstrate a consistent use of category labels, something which has not been demonstrated by other acquisition models.


title: Acceleration of particles to high energy via gravitational repulsion in the Schwarzschild field

abstract: Gravitational repulsion is an inherent aspect of the Schwarzschild solution of the Einstein-Hilbert field equations of general relativity. We show that this circumstance means that it is possible to gravitationally accelerate particles to the highest cosmic ray energies.


title: SmartUnit: Empirical Evaluations for Automated Unit Testing of Embedded Software in Industry

abstract: In this paper, we aim at the automated unit coverage-based testing for embedded software. To achieve the goal, by analyzing the industrial requirements and our previous work on automated unit testing tool CAUT, we rebuild a new tool, SmartUnit, to solve the engineering requirements that take place in our partner companies. SmartUnit is a dynamic symbolic execution implementation, which supports statement, branch, boundary value and MC/DC coverage. SmartUnit has been used to test more than one million lines of code in real projects. For confidentiality motives, we select three in-house real projects for the empirical evaluations. We also carry out our evaluations on two open source database projects, SQLite and PostgreSQL, to test the scalability of our tool since the scale of the embedded software project is mostly not large, 5K-50K lines of code on average. From our experimental results, in general, more than 90% of functions in commercial embedded software achieve 100% statement, branch, MC/DC coverage, more than 80% of functions in SQLite achieve 100% MC/DC coverage, and more than 60% of functions in PostgreSQL achieve 100% MC/DC coverage. Moreover, SmartUnit is able to find the runtime exceptions at the unit testing level. We also have reported exceptions like array index out of bounds and divided-by-zero in SQLite. Furthermore, we analyze the reasons of low coverage in automated unit testing in our setting and give a survey on the situation of manual unit testing with respect to automated unit testing in industry.


title: On VOAs Associated to Jordan Algebras of Type $C$

abstract: In this article, we construct certain universal VOAs whose Greiss algebras are type C Jordan algebras. We also prove the corresponding simplicity result.


title: DeePCG: constructing coarse-grained models via deep neural networks

abstract: We introduce a general framework for constructing coarse-grained potential models without ad hoc approximations such as limiting the potential to two- and/or three-body contributions. The scheme, called Deep Coarse-Grained Potential (abbreviated DeePCG), exploits a carefully crafted neural network to construct a many-body coarse-grained potential. The network is trained with full atomistic data in a way that preserves the natural symmetries of the system. The resulting model is very accurate and can be used to sample the configurations of the coarse-grained variables in a much faster way than with the original atomistic model. As an application we consider liquid water and use the oxygen coordinates as the coarse-grained variables, starting from a full atomistic simulation of this system at the ab-initio molecular dynamics level. We found that the two-body, three-body and higher order oxygen correlation functions produced by the coarse-grained and full atomistic models agree very well with each other, illustrating the effectiveness of the DeePCG model on a rather challenging task.


title: Morrey spaces related to certain nonnegative potentials and fractional integrals on the Heisenberg groups

abstract: Let $\mathcal L=-\Delta_{\mathbb H^n}+V$ be a Schr\"odinger operator on the Heisenberg group $\mathbb H^n$, where $\Delta_{\mathbb H^n}$ is the sub-Laplacian on $\mathbb H^n$ and the nonnegative potential $V$ belongs to the reverse H\"older class $RH_s$ with $s\geq Q/2$. Here $Q=2n+2$ is the homogeneous dimension of $\mathbb H^n$. For given $\alpha\in(0,Q)$, the fractional integrals associated to the Schr\"odinger operator $\mathcal L$ is defined by $\mathcal I_{\alpha}={\mathcal L}^{-{\alpha}/2}$. In this article, we first introduce the Morrey space $L^{p,\kappa}_{\rho,\infty}(\mathbb H^n)$ and weak Morrey space $WL^{p,\kappa}_{\rho,\infty}(\mathbb H^n)$ related to the nonnegative potential $V$. Then we establish the boundedness of fractional integrals ${\mathcal L}^{-{\alpha}/2}$ on these new spaces. Furthermore, in order to deal with certain extreme cases, we also introduce the spaces $\mathrm{BMO}_{\rho,\infty}(\mathbb H^n)$ and $\mathcal{C}^{\beta}_{\rho,\infty}(\mathbb H^n)$ with exponent $\beta\in(0,1]$.


title: On the Chow groups of some hyperkaehler fourfolds with a non-symplectic involution II

abstract: This article is about hyperk\"ahler fourfolds $X$ admitting a non-symplectic involution $\iota$. The Bloch-Beilinson conjectures predict the way $\iota$ should act on certain pieces of the Chow groups of $X$. The main result is a verification of this prediction for Fano varieties of lines on certain cubic fourfolds. This has some interesting consequences for the Chow ring of the quotient $X/\iota$.


title: Hadron Spectroscopy and Dynamics from Light-Front Holography and Superconformal Algebra

abstract: QCD is not supersymmetrical in the traditional sense -- the QCD Lagrangian is based on quark and gluonic fields, not squarks nor gluinos. However, its hadronic eigensolutions conform to a representation of superconformal algebra, reflecting the underlying conformal symmetry of chiral QCD and its Pauli matrix representation. The eigensolutions of superconformal algebra provide a unified Regge spectroscopy of meson, baryon, and tetraquarks in the same 4-plet representation with a universal Regge slope. The pion $q \bar q$ eigenstate has zero mass for $m_q=0.$ The superconformal relations also can be extended to heavy-light quark mesons and baryons. The combined approach of light-front holography and superconformal algebra also provides insight into the origin of the QCD mass scale and color confinement. A key observation is the remarkable dAFF principle which shows how a mass scale can appear in the Hamiltonian and the equations of motion while retaining the conformal symmetry of the action. When one applies the dAFF procedure to chiral QCD, a mass scale $\kappa$ appears which determines universal Regge slopes, hadron masses in the absence of the Higgs coupling, and the mass parameter underlying the form of the nonperturbative QCD running coupling: $\alpha_s(Q^2) \propto \exp{-{Q^2/4 \kappa^2}}$, in agreement with the effective charge determined from measurements of the Bjorken sum rule. The mass scale $\kappa$ underlying hadron masses can be connected to the parameter $\Lambda_{\overline {MS}}$ in the QCD running coupling by matching its predicted nonperturbative form to the perturbative QCD regime. One also obtains predictions for spacelike and timelike hadronic form factors, structure functions, distribution amplitudes, and transverse momentum distributions.


title: Dark Energy Accretion onto Van der Waal's Black Hole

abstract: We consider the most general static spherically symmetric black hole metric. The accretion of the fluid flow around the Van der Waal's black hole is investigated and we calculate the fluid's four-velocity, the critical point and the speed of sound during the accretion process. We also analyze the nature of the universe's density and the mass of the black hole during accretion of the fluid flow. The density of the fluid flow is also taken into account. We observe that the mass is related to redshift. We compare the accreting power of the Van der Waal's black hole with Schwarzschild black hole for different accreting fluid.


title: Semantic Vector Spaces for Broadening Consideration of Consequences

abstract: Reasoning systems with too simple a model of the world and human intent are unable to consider potential negative side effects of their actions and modify their plans to avoid them (e.g., avoiding potential errors). However, hand-encoding the enormous and subtle body of facts that constitutes common sense into a knowledge base has proved too difficult despite decades of work. Distributed semantic vector spaces learned from large text corpora, on the other hand, can learn representations that capture shades of meaning of common-sense concepts and perform analogical and associational reasoning in ways that knowledge bases are too rigid to perform, by encoding concepts and the relations between them as geometric structures. These have, however, the disadvantage of being unreliable, poorly understood, and biased in their view of the world by the source material. This chapter will discuss how these approaches may be combined in a way that combines the best properties of each for understanding the world and human intentions in a richer way.


title: Phase Transition In $SU(N)\times U(1)$ Gauge Theory With Many Fundamental Bosons

abstract: Here we study the Renormalization group flow of $SU(N)\times U(1)$ gauge theory with $M$-fundamental bosons in $4-\epsilon$ dimension by calculating the beta functions. We found a new stable fixed point in the zero mass plane for $M>M_\text{crit}$ by expanding upto $O(\epsilon)$. This indicates a second order phase transition. We also calculated the critical exponents in both $\epsilon$ expansion and also in the large-$M$ expansion.


title: Complexity of finding near-stationary points of convex functions stochastically

abstract: In a recent paper, we showed that the stochastic subgradient method applied to a weakly convex problem, drives the gradient of the Moreau envelope to zero at the rate $O(k^{-1/4})$. In this supplementary note, we present a stochastic subgradient method for minimizing a convex function, with the improved rate $\widetilde O(k^{-1/2})$.


title: Simultaneous Solving of Batched Linear Programs on a GPU

abstract: Linear Programs (LPs) appear in a large number of applications and offloading them to a GPU is viable to gain performance. Existing work on offloading and solving an LP on a GPU suggests that there is performance gain generally on large sized LPs (typically 500 constraints, 500 variables and above). In order to gain performance from a GPU, for applications involving small to medium sized LPs, we propose batched solving of a large number of LPs in parallel. In this paper, we present the design and implementation of a batched LP solver in CUDA, keeping memory coalescent access, low CPU-GPU memory transfer latency and load balancing as the goals. The performance of the batched LP solver is compared against sequential solving in the CPU using the open source solver GLPK (GNU Linear Programming Kit) and the CPLEX solver from IBM. The evaluation on selected LP benchmarks from the Netlib repository displays a maximum speed-up of 95x and 5x with respect to CPLEX and GLPK solver respectively, for a batch of 1e5 LPs. We demonstrate the application of our batched LP solver to enhance performance in the domain of state-space exploration of mathematical models of control systems design.


title: Moore: Interval Arithmetic in C++20

abstract: This article presents the Moore library for interval arithmetic in C++20. It gives examples of how the library can be used, and explains the basic principles underlying its design.


title: Profinite commensurability of S-arithmetic groups

abstract: Given an S-arithmetic group, we ask how much information on the ambient algebraic group, number field of definition, and set of places S is encoded in the commensurability class of the profinite completion. As a first step, we show that the profinite commensurability class of a higher rank S-arithmetic group determines the number field up to arithmetical equivalence and the places in S above unramified primes. We include applications to profiniteness questions of group invariants.


title: Teaching students about informatics and astronomy using real data for detection of asteroids

abstract: In this paper we approach the astronomy teaching process for the students in computer sciences through the controlled investigation method on real astronomical data, using data reduction and quality control of the astrometry of near-Earth asteroids. The method used the data collected on the Isaac Newton Telescope (INT) located at the ORM observatory on the island of La Palma in the Spanish Canary Islands and was successfully tested on a group of students in the second-year of study.


title: Personalized Gaussian Processes for Forecasting of Alzheimer's Disease Assessment Scale-Cognition Sub-Scale (ADAS-Cog13)

abstract: In this paper, we introduce the use of a personalized Gaussian Process model (pGP) to predict per-patient changes in ADAS-Cog13 -- a significant predictor of Alzheimer's Disease (AD) in the cognitive domain -- using data from each patient's previous visits, and testing on future (held-out) data. We start by learning a population-level model using multi-modal data from previously seen patients using a base Gaussian Process (GP) regression. The personalized GP (pGP) is formed by adapting the base GP sequentially over time to a new (target) patient using domain adaptive GPs. We extend this personalized approach to predict the values of ADAS-Cog13 over the future 6, 12, 18, and 24 months. We compare this approach to a GP model trained only on past data of the target patients (tGP), as well as to a new approach that combines pGP with tGP. We find that the new approach, combining pGP with tGP, leads to large improvements in accurately forecasting future ADAS-Cog13 scores.


title: Deep Unsupervised Learning of Visual Similarities

abstract: Exemplar learning of visual similarities in an unsupervised manner is a problem of paramount importance to Computer Vision. In this context, however, the recent breakthrough in deep learning could not yet unfold its full potential. With only a single positive sample, a great imbalance between one positive and many negatives, and unreliable relationships between most samples, training of Convolutional Neural networks is impaired. In this paper we use weak estimates of local similarities and propose a single optimization problem to extract batches of samples with mutually consistent relations. Conflicting relations are distributed over different batches and similar samples are grouped into compact groups. Learning visual similarities is then framed as a sequence of categorization tasks. The CNN then consolidates transitivity relations within and between groups and learns a single representation for all samples without the need for labels. The proposed unsupervised approach has shown competitive performance on detailed posture analysis and object classification.


title: The Parameterized Hardness of the k-Center Problem in Transportation Networks

abstract: In this paper we study the hardness of the $k$-Center problem on inputs that model transportation networks. For the problem, a graph $G=(V,E)$ with edge lengths and an integer $k$ are given and a center set $C\subseteq V$ needs to be chosen such that $|C|\leq k$. The aim is to minimize the maximum distance of any vertex in the graph to the closest center. This problem arises in many applications of logistics, and thus it is natural to consider inputs that model transportation networks. Such inputs are often assumed to be planar graphs, low doubling metrics, or bounded highway dimension graphs. For each of these models, parameterized approximation algorithms have been shown to exist. We complement these results by proving that the $k$-Center problem is W[1]-hard on planar graphs of constant doubling dimension, where the parameter is the combination of the number of centers $k$, the highway dimension $h$, and the pathwidth $p$. Moreover, under the Exponential Time Hypothesis there is no $f(k,p,h)\cdot n^{o(p+\sqrt{k+h})}$ time algorithm for any computable function $f$. Thus it is unlikely that the optimum solution to $k$-Center can be found efficiently, even when assuming that the input graph abides to all of the above models for transportation networks at once! Additionally we give a simple parameterized $(1+\varepsilon)$-approximation algorithm for inputs of doubling dimension $d$ with runtime $(k^k/\varepsilon^{O(kd)})\cdot n^{O(1)}$. This generalizes a previous result, which considered inputs in $D$-dimensional $L_q$ metrics.


title: Improbability of Collisions in $n$-Body Systems

abstract: For a wide class of two-body interactions, including standard examples like gravitational or Coulomb fields, we show that collision orbits in $n$-body systems are of Liouville measure zero for all energies. We use techniques from symplectic geometry to relate the volume of the union of collision orbits to the area of Poincar\'e surfaces surrounding the collision set.


title: Additional contributions from: Nobel Symposium 162 - Microfluidics

abstract: Series of short contributions that are part of Nobel Symposium 162 - Microfluidics arXiv:1712.08369.


title: Improbability of Wandering Orbits Passing Through a Sequence of Poincar\'e Surfaces of Decreasing Size

abstract: Given a volume preserving dynamical system with non-compact phase space, one is sometimes interested in special subsets of its wandering set. One example from celestial mechanics is the set of initial values leading to collision. Another one is the set of initial values of semi-orbits, whose asymptotic velocity does not exist as a limit. We introduce techniques that can be helpful in showing that these sets are of measure zero: by defining a sequence of hypersurfaces, that are eventually hit by each of those semi-orbits and whose total surface area decreases to zero.


title: Adversarial Training for Probabilistic Spiking Neural Networks

abstract: Classifiers trained using conventional empirical risk minimization or maximum likelihood methods are known to suffer dramatic performance degradations when tested over examples adversarially selected based on knowledge of the classifier's decision rule. Due to the prominence of Artificial Neural Networks (ANNs) as classifiers, their sensitivity to adversarial examples, as well as robust training schemes, have been recently the subject of intense investigation. In this paper, for the first time, the sensitivity of spiking neural networks (SNNs), or third-generation neural networks, to adversarial examples is studied. The study considers rate and time encoding, as well as rate and first-to-spike decoding. Furthermore, a robust training mechanism is proposed that is demonstrated to enhance the performance of SNNs under white-box attacks.


title: Indic Handwritten Script Identification using Offline-Online Multimodal Deep Network

abstract: In this paper, we propose a novel approach of word-level Indic script identification using only character-level data in training stage. The advantages of using character level data for training have been outlined in section I. Our method uses a multimodal deep network which takes both offline and online modality of the data as input in order to explore the information from both the modalities jointly for script identification task. We take handwritten data in either modality as input and the opposite modality is generated through intermodality conversion. Thereafter, we feed this offline-online modality pair to our network. Hence, along with the advantage of utilizing information from both the modalities, it can work as a single framework for both offline and online script identification simultaneously which alleviates the need for designing two separate script identification modules for individual modality. One more major contribution is that we propose a novel conditional multimodal fusion scheme to combine the information from offline and online modality which takes into account the real origin of the data being fed to our network and thus it combines adaptively. An exhaustive experiment has been done on a data set consisting of English and six Indic scripts. Our proposed framework clearly outperforms different frameworks based on traditional classifiers along with handcrafted features and deep learning based methods with a clear margin. Extensive experiments show that using only character level training data can achieve state-of-art performance similar to that obtained with traditional training using word level data in our framework.


title: Toward High Temperature Quasi-two-dimensional Superconductivity

abstract: The demonstration of a quasi-two-dimensional electron gas (2DEG) and superconducting properties in LaAlO$_3$/SrTiO$_3$ heterostructures has stimulated intense research activity in recent ten years. The 2DEG has unique properties that are promising for applications in all-oxide electronic devices. The superconductivity in such heterostructures has been observed below 300 mK. For superconductivity applications it is desirable to have more wide temperature of the existence range and the ability to control superconductivity properties by external stimulus. Based on first-principles calculations and theoretical consideration we show that all-oxide heterostructures incorporating ferroelectric constituent, such as BaTiO$_3$/La$_2$CuO$_4$, allow creating 2DEG. We predict a possibility of a high temperature guasi-two-dimensional superconductivity state. This state could be switchable between superconducting and conducting states by ferroelectric polarization reversal. We also discuss that such structures must be more simple for preparation. The proposed concept of ferroelectrically controlled interface superconductivity offers the possibility to design novel electronic devices.


title: Relative hyperbolicity of free-by-cyclic extensions

abstract: Given a finite rank free group $\mathbb{F}$ of $\mathsf{rank}(\mathbb{F})\geq 3$, we show that the mapping torus of $\phi$ is (strongly) relatively hyperbolic if $\phi$ is exponentially growing. We combine our result with the work of Button-Kropholler to answer a question asked by Minasyan-Osin regarding the acylindrical hyperbolicity of such free-by-cyclic extensions. As an application we construct new examples of free-by-free hyperbolic extensions where the elements of the quotient group are not necessarily fully irreducible. We also give a new proof of the Bridson-Groves quadratic isoperimetric inequality theorem.


title: Energy consumption for ion transport in a segmented Paul trap

abstract: There is recent interest in determining energy costs of shortcuts to adiabaticity (STA), but different definitions of "cost" have been used. We demonstrate the importance of taking into account the Control System (CS) for a fair assessment of energy flows and consumptions. We model the energy consumption and power to transport an ion by a STA protocol in a multisegmented Paul trap. The ion is driven by an externally controlled, moving harmonic oscillator. Even if no net ion- energy is gained at destination, setting the time-dependent control parameters is a macroscopic operation that costs energy and results in energy dissipation for the short time scales implied by the intrinsically fast STA processes. The potential minimum is displaced by modulating the voltages on control (dc) electrodes. A secondary effect of the modulation, usually ignored as it does not affect the ion dynamics, is the time- dependent energy shift of the potential minimum. The non trivial part of the energy consumption is due to the electromotive forces to set the electrode voltages through the low-pass filters required to preserve the electronic noise from decohering the ion's motion. The results for the macroscopic CS (the Paul trap) are compared to the microscopic power and energy of the ion alone. Similarities are found -and may be used quantitatively to minimize costs- only when the CS-dependent energy shift of the harmonic oscillator is included in the ion energy.


title: Bayesian Lasso : Concentration and MCMC Diagnosis

abstract: Using posterior distribution of Bayesian LASSO we construct a semi-norm on the parameter space. We show that the partition function depends on the ratio of the l 1 and l 2 norms and present three regimes. We derive the con- centration of Bayesian LASSO, and present MCMC convergence diagnosis. Keywords: LASSO, Bayes, MCMC, log-concave, geometry, incomplete Gamma function


title: Higher order Seiberg-Witten functionals and their associated gradient flows

abstract: We define functionals generalising the Seiberg-Witten functional on closed $spin^c$ manifolds, involving higher order derivatives of the curvature form and spinor field. We then consider their associated gradient flows and, using a gauge fixing technique, are able to prove short time existence for the flows. We then prove energy estimates along the flow, and establish local $L^2$-derivative estimates. These are then used to show long time existence of the flow in sub-critical dimensions. In the critical dimension, we are able to show that long time existence is obstructed by an $L^{k+2}$ curvature concentration phenomenon.


title: Strong measure zero and meager-additive sets through the prism of fractal measures

abstract: We develop a theory of \emph{sharp measure zero} sets that parallels Borel's \emph{strong measure zero}, and prove a theorem analogous to Galvin-Myscielski-Solovay Theorem, namely that a set of reals has sharp measure zero if and only if it is meager-additive. Some consequences: A subset of $2^\omega$ is meager-additive if and only if it is $\mathcal E$-additive; if $f:2^\omega\to2^\omega$ is continuous and $X$ is meager-additive, then so is $f(X)$.


title: Complexity, Centralization, and Fragility in Economic Networks

abstract: Trade networks, across which countries distribute their products, are crucial components of the globalized world economy. Their structure is strongly heterogeneous across products, given the different features of the countries which buy and sell goods. By using a diversified pool of indicators from network science and product complexity theory, we quantitatively confirm the intuition that, overall, products with higher complexity -- i.e., with larger technological content and number of components -- are traded through a more centralized network -- i.e., with a small number of countries concentrating most of the export flow. Since centralized networks are known to be more vulnerable, we argue that the current composition of production and trading is associated to high fragility at the level of the most complex -- thus strategic -- products.


title: Volume average regularization for the Wheeler-DeWitt equation

abstract: In this article, I present a volume average regularization for the second functional derivative operator that appears in the metric-basis Wheeler-DeWitt equation. Naively, the second functional derivative operator in the Wheeler-DeWitt equation is infinite, since it contains terms with a factor of a delta function or derivatives of the delta function. More precisely, the second functional derivative contains terms that are only well defined as a distribution---these terms only yield meaningful results when they appear within an integral. The second functional derivative may, therefore, be regularized by performing an integral average of the distributional terms over some finite volume; I argue that such a regularization is appropriate if one regards quantum general relativity (from which the Wheeler-DeWitt equation may be derived) to be the low-energy effective field theory of a full theory of quantum gravity. I also show that a volume average regularization can be viewed as a natural generalization of the same-variable second partial derivative for an ordinary multivariable function. Using the regularized second functional derivative operator, I construct an approximate solution to the Wheeler-DeWitt equation in the low-curvature, long-distance limit.


title: On Streaming Algorithms for the Steiner Cycle and Path Cover Problem on Interval Graphs and Falling Platforms in Video Games

abstract: We introduce a simplified model for platform game levels with falling platforms based on interval graphs and show that solvability of such levels corresponds to finding Steiner cycles or Steiner paths in the corresponding graphs. Linear time algorithms are obtained for both of these problems. We also study these algorithms as streaming algorithms and analyze the necessary memory with respect to the maximum number of intervals contained in another interval. This corresponds to understanding which parts of a level have to be visible at each point to allow the player to make optimal deterministic decisions.


title: "Isocrater" impacts: Conditions and mantle dynamical responses for different impactor types

abstract: Impactors of different types and sizes can produce a final crater of the same diameter on a planet under certain conditions. We derive the condition for such "isocrater impacts" from scaling laws, as well as relations that describe how the different impactors affect the interior of the target planet; these relations are also valid for impacts that are too small to affect the mantle. The analysis reveals that in a given isocrater impact, asteroidal impactors produce anomalies in the interior of smaller spatial extent than cometary or similar impactors. The differences in the interior could be useful for characterizing the projectile that formed a given crater on the basis of geophysical observations and potentially offer a possibility to help constrain the demographics of the ancient impactor population. A series of numerical models of basin-forming impacts on Mercury, Venus, the Moon, and Mars illustrates the dynamical effects of the different impactor types on different planets. It shows that the signature of large impacts may be preserved to the present in Mars, the Moon, and Mercury, where convection is less vigorous and much of the anomaly merges with the growing lid. On the other hand, their signature will long have been destroyed in Venus, whose vigorous convection and recurring lithospheric instabilities obliterate larger coherent anomalies.


title: Nonparametric Estimation of a distribution function from doubly truncated data under dependence

abstract: The NPMLE of a distribution function from doubly truncated data was introduced in the seminal paper of Efron and Petrosian. The consistency of the Efron-Petrosian estimator depends however on the assumption of independent truncation. In this work we introduce an extension of the Efron-Petrosian NPMLE when the variable of interest and the truncation variables may be dependent. The proposed estimator is constructed on the basis of a copula function which represents the dependence structure between the variable of interest and the truncation variables. Two different iterative algorithms to compute the estimator in practice are introduced, and their performance is explored through an intensive Monte Carlo simulation study. We illustrate the use of the estimators on two real data examples.


title: Quantum Nonlocality and Indistinguishability

abstract: It is shown how a "meddlesome" photon indistinguishable from another photon of an entangled pair can affect the result of an Einstein- Podolsky-Rosen (EPR) experiment. This makes it clear the importance of the notion of field over that of particle.


title: X-ray study of a sample of FR0 radio galaxies: unveiling the nature of the central engine

abstract: FR0s are compact radio sources that represent the bulk of the Radio-Loud (RL) AGN population, but they are still poorly understood. Pilot studies on these sources have been already performed at radio and optical wavelengths: here we present the first X-ray study of a sample of 19 FR0 radio galaxies selected from the SDSS/NVSS/FIRST sample of Best & Heckman (2012), with redshift $\leq$ 0.15, radio size $\leq$ 10 kpc and optically classified as low-excitation galaxies (LEG). The X-ray spectra are modeled with a power-law component absorbed by Galactic column density with, in some cases, a contribution from thermal extended gas. The X-ray photons are likely produced by the jet as attested by the observed correlation between X-ray (2-10 keV) and radio (5 GHz) luminosities, similar to FRIs. The estimated Eddington-scaled luminosities indicate a low accretion rate. Overall, we find that the X-ray properties of FR0s are indistinguishable from those of FRIs, thus adding another similarity between AGN associated with compact and extended radio sources. A comparison between FR0s and low luminosity BL Lacs, rules out important beaming effects in the X-ray emission of the compact radio galaxies. FR0s have different X-ray properties with respect to young radio sources (e.g. GPS/CSS sources), generally characterized by higher X-ray luminosities and more complex spectra. In conclusion, the paucity of extended radio emission in FR0s is probably related to the intrinsic properties of their jets that prevent the formation of extended structures, and/or to intermittent activity of their engines.


title: Magnetic field generation in plasma waves driven by co-propagating intense twisted lasers

abstract: We present a new magnetic field generation mechanism in underdense plasmas driven by the beating of two, co-propagating, Laguerre-Gaussian (LG) orbital angular momentum (OAM) laser pulses with different frequencies and also different twist indices. The resulting twisted ponderomotive force drives up an electron plasma wave with a helical rotating structure. To second order, there is a nonlinear rotating current leading to the onset of an intense, static axial magnetic field, which persists over a long time in the plasma (ps scale) after the laser pulses have passed by. The results are confirmed in three-dimensional particle-in-cell simulations and also theoretical analysis. For the case of 300 fs duration, 3.8x10^17 W/cm^2 peak laser intensity we observe magnetic field of up to 0.4 MG. This new method of magnetic field creation may find applications in charged beam collimation and controlled fusion.


title: Feedback Control of Scalar Conservation Laws with Application to Density Control in Freeways by Means of Variable Speed Limits

abstract: The paper provides results for the stabilization of a spatially uniform equilibrium profile for a scalar conservation law that arises in the study of traffic dynamics under variable speed limit control. Two different control problems are studied: the problem with free speed limits at the inlet and the problem with no speed limits at the inlet. Explicit formulas are provided for respective feedback laws that guarantee stabilization of the desired equilibrium profile. For the first problem, global asymptotic stabilization is achieved; while for the second problem, regional exponential stabilization is achieved. Moreover, the solutions for the corresponding closed-loop systems are guaranteed to be classical solutions, i.e., there are no shocks. The obtained results are illustrated by means of a numerical example.


title: Few- and many-nucleon systems with semilocal coordinate-space regularized chiral nucleon-nucleon forces

abstract: We employ a variety of ab initio methods including Faddeev-Yakubovsky equations, No-Core Configuration Interaction Approach, Coupled-Cluster Theory and In-Medium Similarity Renormalization Group to perform a comprehensive analysis of the nucleon-deuteron elastic and breakup reactions and selected properties of light and medium-mass nuclei up to 48Ca using the recently constructed semilocal coordinate-space regularized chiral nucleon-nucleon potentials. We compare the results with those based on selected phenomenological and chiral EFT two-nucleon potentials, discuss the convergence pattern of the chiral expansion and estimate the achievable theoretical accuracy at various chiral orders using the novel approach to quantify truncation errors of the chiral expansion without relying on cutoff variation. We also address the robustness of this method and explore alternative ways to estimate the theoretical uncertainty from the truncation of the chiral expansion.


title: The statistical properties of solar wind temperature parameters near 1 AU

abstract: We present a long-duration ($\sim$10 years) statistical analysis of the temperatures, plasma betas, and temperature ratios for the electron, proton, and alpha-particle populations observed by the \emph{Wind} spacecraft near 1 AU. The mean(median) scalar temperatures are $T{\scriptstyle_{e, tot}}$ $=$ 12.2(11.9) eV, $T{\scriptstyle_{p, tot}}$ $=$ 12.7(8.6) eV, and $T{\scriptstyle_{\alpha, tot}}$ $=$ 23.9(10.8) eV. The mean(median) total plasma betas are $\beta{\scriptstyle_{e, tot}}$ $=$ 2.31(1.09), $\beta{\scriptstyle_{p, tot}}$ $=$ 1.79(1.05), and $\beta{\scriptstyle_{\alpha, tot}}$ $=$ 0.17(0.05). The mean(median) temperature ratios are $\left(T{\scriptstyle_{e}}/T{\scriptstyle_{p}}\right){\scriptstyle_{tot}}$ $=$ 1.64(1.27), $\left(T{\scriptstyle_{e}}/T{\scriptstyle_{\alpha}}\right){\scriptstyle_{tot}}$ $=$ 1.24(0.82), and $\left(T{\scriptstyle_{\alpha}}/T{\scriptstyle_{p}}\right){\scriptstyle_{tot}}$ $=$ 2.50(1.94). We also examined these parameters during time intervals that exclude interplanetary (IP) shocks, times within the magnetic obstacles (MOs) of interplanetary coronal mass ejections (ICMEs), and times that exclude MOs. The only times that show significant alterations to any of the parameters examined are those during MOs. In fact, the only parameter that does not show a significant change during MOs is the electron temperature. Although each parameter shows a broad range of values, the vast majority are near the median. We also compute particle-particle collision rates and compare to effective wave-particle collision rates. We find that, for reasonable assumptions of wave amplitude and occurrence rates, the effect of wave-particle interactions on the plasma is equal to or greater than the effect of Coulomb collisions. Thus, wave-particle interactions should not be neglected when modeling the solar wind.


title: Database Aggregation

abstract: Knowledge can be represented compactly in a multitude ways, from a set of propositional formulas, to a Kripke model, to a database. In this paper we study the aggregation of information coming from multiple sources, each source submitting a database modelled as a first-order relational structure. In the presence of an integrity constraint, we identify classes of aggregators that respect it in the aggregated database, provided all individual databases satisfy it. We also characterise languages for first-order queries on which the answer to queries on the aggregated database coincides with the aggregation of the answers to the query obtained on each individual database. This contribution is meant to be a first step on the application of techniques from rational choice theory to knowledge representation in databases.


title: An ionic impurity in a Bose-Einstein condensate at sub-microkelvin temperatures

abstract: Rydberg atoms immersed in a Bose-Einstein condensate interact with the quantum gas via electron-atom and ion-atom interaction. To suppress the typically dominant electron-neutral interaction, Rydberg states with principal quantum number up to $n = 190$ are excited from a dense and tightly trapped micron-sized condensate. This allows us to explore a regime where the Rydberg orbit exceeds the size of the atomic sample by far. In this case, a detailed lineshape analysis of the Rydberg excitation spectrum provides clear evidence for ion-atom interaction at temperatures well below a microkelvin. Our results may open up ways to enter the quantum regime of ion-atom scattering for the exploration of charged quantum impurities and associated polaron physics.


title: Implicit definability of truth constants in {\L}ukasiewicz logic

abstract: In the framework of propositional {\L}ukasiewicz logic, a suitable notion of implicit definability, tailored to the intended real-valued semantics and referring to the elements of its domain, is introduced. Several variants of implicitly defining each of the rational elements in the standard semantics are explored, and based on that, a faithful interpretation of theories in Rational Pavelka logic in theories in {\L}ukasiewicz logic is obtained. Some of these results were already presented in H\'ajek's "Metamathematics of Fuzzy Logic" as technical statements. A connection to the lack of (deductive) Beth property in {\L}ukasiewicz logic is drawn. Moreover, while irrational elements of the standard semantics are not implicitly definable by finitary means, a parallel development is possible for them in the infinitary {\L}ukasiewicz logic. As an application of definability of the rationals, it is shown how computational complexity results for Rational Pavelka logic can be obtained from analogous results for {\L}ukasiewicz logic. The complexity of the definability notion itself is studied as well. Finally, we review the import of these results for the precision/vagueness discussion for fuzzy logic, and for the general standing of truth constants in {\L}ukasiewicz logic.


title: Asymptotic safety in the dark

abstract: We explore the Renormalization Group flow of massive uncharged fermions -- a candidate for dark matter -- coupled to a scalar field through a Higgs portal. We find that fermionic fluctuations can lower the bound on the scalar mass that arises from vacuum stability. Further, we discuss that despite the perturbative nonrenormalizability of the model, it could be ultraviolet complete at an asymptotically safe fixed point. In our approximation, this simple model exhibits two mechanisms for asymptotic safety: a balance of fermionic and bosonic fluctuations generates a fixed point in the scalar self-interaction; asymptotic safety in the portal coupling is triggered through a balance of canonical scaling and quantum fluctuations. As a consequence of asymptotic safety in the dark sector, the low-energy value of the portal coupling could become a function of the dark fermion mass and the scalar mass, thereby reducing the viable parameter space of the model.


title: Reservoir computing with simple oscillators: Virtual and real networks

abstract: The reservoir computing scheme is a machine learning mechanism which utilizes the naturally occuring computational capabilities of dynamical systems. One important subset of systems that has proven powerful both in experiments and theory are delay-systems. In this work, we investigate the reservoir computing performance of hybrid network-delay systems systematically by evaluating the NARMA10 and the Sante Fe task.. We construct 'multiplexed networks' that can be seen as intermediate steps on the scale from classical networks to the 'virtual networks' of delay systems. We find that the delay approach can be extended to the network case without loss of computational power, enabling the construction of faster reservoir computing substrates.


title: Total Array Gains of Millimeter-Wave Mobile Phone Antennas Under Practical Conditions

abstract: This paper studies a gain of an antenna array embedded on a mobile device operating at a millimeter-wave radio frequency. Assuming that mobile phones at millimeter-wave range operate with a single baseband unit and analog beamforming like phased arrays, we define a total array gain denoting a path gain of the phased antenna array in excess to the omni-directional path gain. The total array gain circumvents the ambiguity of conventional array gain which cannot be uniquely defined as there are multiple choices of a reference single-element antenna in an array. Two types of 8-element patch antenna arrays implemented on a mobile phone chassis, i.e., uniform linear array (ULA) and distributed array (DA) both operating at 60 GHz, are studied. The gain evaluated in a small-cell scenario in an airport shows that DA achieves higher median and outage gain by up to 8 and 6 dB than ULA when different orientations of the mobile phone are considered along with body torso and finger shadowing. There are always postures of the mobile phone where ULA cannot see the line-of-sight due to directionality of the patch antenna and of body and finger shadowing, leading to outage gain of -15 dB in the worst case. The DA has much smaller variation of the gain across different orientations of the phone, even when the human torso shadowing and user's finger effects are considered.


title: On a family of representations of residually finite groups

abstract: For a residually finite group $G$, its normal subgroups $G\supset G_1\supset G_2\cdots$ with $\cap_{n\in\mathbb N}G_n=\{e\}$ and for a growth function $\gamma$ we construct a unitary representation $\pi_\gamma$ of $G$. For the minimal growth, $\pi_\gamma$ is weakly equivalent to the regular representation, and for the maximal growth it is weakly equivalent to the direct sum of the quasiregular representations on the quotients $G/G_n$. In the case of intermediate growth we show two examples of different behaviour of $\pi_\gamma$.


title: Study of phonons in irradiated epitaxial thin films of UO$_2$

abstract: We report experiments to determine the effect of radiation damage on the phonon spectra of the most common nuclear fuel, UO$_2$. We have irradiated thin ($\sim$ 300 nm) epitaxial films of UO$_2$ with 2.1 MeV He$^{2+}$ ions to 0.15 dpa and a lattice swelling of $\Delta$a/a $\sim$ 0.6 %, and then used grazing-incidence inelastic X-ray scattering to measure the phonon spectrum. We succeeded to observe the acoustic modes, both transverse and longitudinal, across the Brillouin zone. The phonon energies, in both the pristine and irradiated samples, are unchanged from those observed in bulk material. On the other hand, the phonon linewidths (inversely proportional to the phonon lifetimes), show a significant broadening when comparing the pristine and irradiated samples. This effect is shown to increase with phonon energy across the Brillouin zone. The decreases in the phonon lifetimes of the acoustic modes are roughly consistent with a 50 % reduction in the thermal conductivity.


title: Links of sandwiched surface singularities and self-similarity

abstract: We characterize sandwiched singularities in terms of their link in two different settings. We first prove that such singularities are precisely the normal surface singularities having self-similar non-archimedean links. We describe this self-similarity both in terms of Berkovich analytic geometry and of the combinatorics of weighted dual graphs. We then show that a complex surface singularity is sandwiched if and only if its complex link can be embedded in a Kato surface in such a way that its complement remains connected.


title: Time-Varying Block Codes for Synchronization Errors: MAP Decoder and Practical Issues

abstract: \begin{abstract} In this paper we consider Time-Varying Block (TVB) codes, which generalize a number of previous synchronization error-correcting codes. We also consider various practical issues related to MAP decoding of these codes. Specifically, we give an expression for the expected distribution of drift between transmitter and receiver due to synchronization errors. We determine an appropriate choice for state space limits based on the drift probability distribution. In turn, we obtain an expression for the decoder complexity under given channel conditions in terms of the state space limits used. For a given state space, we also give a number of optimizations that reduce the algorithm complexity with no further loss of decoder performance. We also show how the MAP decoder can be used in the absence of known frame boundaries, and demonstrate that an appropriate choice of decoder parameters allows the decoder to approach the performance when frame boundaries are known, at the expense of some increase in complexity. Finally, we express some existing constructions as TVB codes, comparing performance with published results, and showing that improved performance is possible by taking advantage of the flexibility of TVB codes.


title: Static memory materials

abstract: We simulate static memory materials on a two-dimensional lattice. The bulk properties of such materials depend on boundary conditions. Considerable information can be stored in various local patterns. We observe local probabilities oscillating with the distance from the boundary. The dependence of the local statistical information on this distance can be described by a linear evolution of classical wave functions, including the superposition principle and classical interference. We speculate that these new phenomena could open new algorithmic possibilities analogous to quantum computing.


title: Constraining the optical potential in the search for eta-mesic 4He

abstract: A consistent description of the $ d d \rightarrow$ $^{4}\hspace{-0.03cm}\mbox{He} \eta$ and $ d d \rightarrow$ ($^4$He$\eta$)$_{bound} \rightarrow X $ cross sections was recently proposed with a broad range of real ($V_{0}$) and imaginary (W0), $\eta$-$^4$He optical potential parameters leading to a good agreement with the d d -> 4He eta data. Here we compare the predictions of the model below the $\eta$ production threshold, with the WASA-at-COSY excitation functions for the dd -> 3He N pi reactions to put stronger constraints on (V0, W0). The allowed parameter space (with |V_0| < \sim 60 MeV and |W_0| < \sim 7 MeV estimated at 90% CL ) excludes most optical model predictions of eta-4He nuclei except for some loosely bound narrow states.


title: Learning Weighted Representations for Generalization Across Designs

abstract: Predictive models that generalize well under distributional shift are often desirable and sometimes crucial to building robust and reliable machine learning applications. We focus on distributional shift that arises in causal inference from observational data and in unsupervised domain adaptation. We pose both of these problems as prediction under a shift in design. Popular methods for overcoming distributional shift make unrealistic assumptions such as having a well-specified model or knowing the policy that gave rise to the observed data. Other methods are hindered by their need for a pre-specified metric for comparing observations, or by poor asymptotic properties. We devise a bound on the generalization error under design shift, incorporating both representation learning and sample re-weighting. Based on the bound, we propose an algorithmic framework that does not require any of the above assumptions and which is asymptotically consistent. We empirically study the new framework using two synthetic datasets, and demonstrate its effectiveness compared to previous methods.


title: Evaluating Scoped Meaning Representations

abstract: Semantic parsing offers many opportunities to improve natural language understanding. We present a semantically annotated parallel corpus for English, German, Italian, and Dutch where sentences are aligned with scoped meaning representations in order to capture the semantics of negation, modals, quantification, and presupposition triggers. The semantic formalism is based on Discourse Representation Theory, but concepts are represented by WordNet synsets and thematic roles by VerbNet relations. Translating scoped meaning representations to sets of clauses enables us to compare them for the purpose of semantic parser evaluation and checking translations. This is done by computing precision and recall on matching clauses, in a similar way as is done for Abstract Meaning Representations. We show that our matching tool for evaluating scoped meaning representations is both accurate and efficient. Applying this matching tool to three baseline semantic parsers yields F-scores between 43% and 54%. A pilot study is performed to automatically find changes in meaning by comparing meaning representations of translations. This comparison turns out to be an additional way of (i) finding annotation mistakes and (ii) finding instances where our semantic analysis needs to be improved.


title: Multifrequency superscattering from subwavelength hyperbolic structures

abstract: Superscattering, i.e., a phenomenon of the scattering cross section from a subwavelength object exceeding the single-channel limit, has important prospects in enhanced sensing/spectroscopy, solar cells, and biomedical imaging. Superscattering can be typically constructed only at a single frequency regime, and depends critically on the inescapable material losses. Under such realistic conditions, superscattering has not been predicted nor observed to exist simultaneously at multiple frequency regimes. Here we introduce multifrequency superscattering in a subwavelength hyperbolic structure, which can be made from artificial metamaterials or from naturally-existing materials, such as hexagonal boron nitride (BN), and show the advantage of such hyperbolic materials for reducing structural complexity. The underlying mechanism is revealed to be the multimode resonances at multiple frequency regimes as appear in BN due to the peculiar dispersion of phonon-polaritons. Importantly, the multifrequency superscattering has a high tolerance to material losses and some structural variations, bringing the concept of multifrequency superscattering closer to useful and realistic conditions.


title: 8T SRAM Cell as a Multi-bit Dot Product Engine for Beyond von-Neumann Computing

abstract: Large scale digital computing almost exclusively relies on the von-Neumann architecture which comprises of separate units for storage and computations. The energy expensive transfer of data from the memory units to the computing cores results in the well-known von-Neumann bottleneck. Various approaches aimed towards bypassing the von-Neumann bottleneck are being extensively explored in the literature. Emerging non-volatile memristive technologies have been shown to be very efficient in computing analog dot products in an in-situ fashion. The memristive analog computation of the dot product results in much faster operation as opposed to digital vector in-memory bit-wise Boolean computations. However, challenges with respect to large scale manufacturing coupled with the limited endurance of memristors have hindered rapid commercialization of memristive based computing solutions. In this work, we show that the standard 8 transistor (8T) digital SRAM array can be configured as an analog-like in-memory multi-bit dot product engine. By applying appropriate analog voltages to the read-ports of the 8T SRAM array, and sensing the output current, an approximate analog-digital dot-product engine can be implemented. We present two different configurations for enabling multi-bit dot product computations in the 8T SRAM cell array, without modifying the standard bit-cell structure. Since our proposal preserves the standard 8T-SRAM array structure, it can be used as a storage element with standard read-write instructions, and also as an on-demand analog-like dot product accelerator.


title: GAtor: A First Principles Genetic Algorithm for Molecular Crystal Structure Prediction

abstract: We present the implementation of GAtor, a massively parallel, first principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several sub-populations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a $Z^\prime$=2 structure with P$\bar{1}$ symmetry and a scaffold packing motif, which has not been reported previously.


title: Towards Distributed OPF using ALADIN

abstract: The present paper discusses the application of the recently proposed Augmented Lagrangian Alternating Direction Inexact Newton (ALADIN) method to non-convex AC Optimal Power Flow Problems (OPF) in a distributed fashion. In contrast to the often used Alternating Direction of Multipliers Method (ADMM), ALADIN guarantees locally quadratic convergence for AC OPF. Numerical results for 5 to 300 bus test cases indicate that ALADIN is able to outperform ADMM and to reduce the number of iterations by about one order of magnitude. We compare ALADIN to numerical results for ADMM documented in the literature. The improved convergence speed comes at the cost of increasing the communication effort per iteration. Therefore, we propose a variant of ALADIN that uses inexact Hessians to reduce communication. Additionally, we provide a detailed comparison of these ALADIN variants to ADMM from an algorithmic and communication perspective. Moreover, we prove that ALADIN converges locally at quadratic rate even for the relevant case of suboptimally solved local NLPs.


title: Limiting gaming opportunities on incentive-based demand response programs

abstract: Demand Response (DR) is a program designed to match supply and demand by modifying consumption profile. Some of these programs are based on economic incentives, in which, a user is paid to reduce his energy requirements according to an estimated baseline. Literature review and practice have shown that the counter-factual models of employing baselines are vulnerable for gaming. Classical solutions of mechanism design require that agents communicate their full types which result in greater difficulties for its practical implementation. In this paper, a novel contract is developed to induce individual rationality (voluntary participation) and asymptotic incentive-compatibility (truthfulness) through probability of call, where an agent does not require to report the marginal utility. In this approach, a consumer only announces the baseline and reduction capacity, given a payment scheme that includes cost of electricity, incentive price, and penalty caused by any deviation between self-reported and actual energy consumption. The aggregator decides randomly what users are called to perform the energy reduction. As result, asymptotic truth-telling behavior in incentive-based DR is managed by the aggregator through the probability of call for each agent. Mathematical proofs and numerical studies are provided to demonstrate the properties and advantages of this contract in limiting gaming opportunities and in terms of its implementation.


title: A note on the discrete Cauchy-Kovalevskaya extension

abstract: In this paper we exploit the umbral calculus framework to reformulate the so-called discrete Cauchy-Kovalevskaya extension in the scope of hypercomplex variables. The key idea is to consider not only formal power series representation for the underlying solution, but also integral representations for the Chebyshev polynomials of first and second kind by means of its Cauchy principal values. It turns out that the resulting integral representation associated to our toy problem is a space-time Fourier type inversion formula. Moreover, with the aid of some Laplace transform identities involving the generalized Mittag-Leffler function we are able to establish a link with a Cauchy problem of differential-difference type.


title: Speed of gravitational waves and black hole hair

abstract: The recent detection of GRB 170817A and GW170817 constrains the speed of gravity waves $c_T$ to be that of light, which severely restricts the landscape of modified gravity theories that impact the cosmological evolution of the universe. In this work, we investigate the presence of black hole hair in the remaining viable cosmological theories of modified gravity that respect the constraint $c_T=1$. We focus mainly on scalar-tensor theories of gravity, analyzing static, asymptotically flat black holes in Horndeski, Beyond Horndeski, Einstein-Scalar-Gauss-Bonnet, and Chern-Simons theories. We find that in all of the cases considered here, theories that respect $c_T=1$ do not allow for hair, or have negligible hair. We further comment on vector-tensor theories including Einstein Yang-Mills, Einstein-Aether, and Generalised Proca theories, as well as bimetric theories.


title: Photon Statistics in Collective Strong Coupling: Nano- and Microcavities

abstract: There exists a growing interest in the properties of the light generated by hybrid systems involving a mesoscopic number of emitters as a means of providing macroscopic quantum light sources. In this work, the quantum correlations of the light emitted by a collection of emitters coupled to a generic optical cavity are studied theoretically using an effective Hamiltonian approach. Starting from the single-emitter level, we analyse the persistence of photon antibunching as the ensemble size increases. Not only is the photon blockade effect identifiable, but photon antibunching originated from destructive interference processes (the so-called unconventional antibunching) is also present. We study the dependence of these two types of negative correlations on the spectral detuning between cavity and emitters, as well as its evolution as the time delay between photon detections increases. Throughout this work, the performance of plasmonic nanocavities and dielectric microcavities is compared: despite the distinct energy scales and the differences introduced by their respectively open and closed character, the bunching and antibunching phenomenology presents remarkable similarities in both types of cavities.