title: Stratospheric Aerosol Particles Size Distribution Based on Multi-Color Polarization Measurements of the Twilight Sky

abstract: Polarization measurements of the twilight background with Wide-Angle Polarization Camera (WAPC) are used to detect the depolarization effect caused by stratospheric aerosol near the altitude of 20 km. Based on a number of observations in central Russia in spring and summer 2016, we found the parameters of lognormal size distribution of aerosol particles. This confirmed the previously published results of the colorimetric method as applied to the same twilights. The mean particle radius (about 0.1 micrometers) and size distribution are also in agreement with the recent data of in situ and space-based remote sensing of stratospheric aerosol. Methods considered here provide two independent techniques of the stratospheric aerosol study based on the twilight sky analysis.


title: Maximizing Wiener Index for Trees with Given Vertex Weight and Degree Sequences

abstract: The Wiener index is maximized over the set of trees with the given vertex weight and degree sequences. This model covers the traditional "unweighed" Wiener index, the terminal Wiener index, and the vertex distance index. It is shown that there exists an optimal caterpillar. If weights of internal vertices increase in their degrees, then an optimal caterpillar exists with weights of internal vertices on its backbone monotonously increasing from some central point to the ends of the backbone, and the same is true for pendent vertices. A tight upper bound of the Wiener index value is proposed and an efficient greedy heuristics is developed that approximates well the optimal index value. Finally, a branch and bound algorithm is built and tested for the exact solution of this NP-complete problem.


title: A family of charged compact objects with anisotropic pressure

abstract: Utilizing an ansatz developed by Maurya and co-workers we present a class of exact solutions of the Einstein-Maxwell field equations describing a spherically symmetric compact object. A detailed physical analysis of these solutions in terms of stability, compactness and regularity indicate that these solutions may be used to model strange star candidates. In particular, we model the strange star candidate Her X-1 and show that our solution conforms to observational data to an excellent degree of accuracy. An interesting and novel phenomenon which arises in this model is the fact that the relative difference between the electromagnetic force and the force due to pressure anisotropy changes sign within the stellar interior. This may be a an additional mechanism required for stability against cracking of the stellar object.


title: Bayesian Approaches to Distribution Regression

abstract: Distribution regression has recently attracted much interest as a generic solution to the problem of supervised learning where labels are available at the group level, rather than at the individual level. Current approaches, however, do not propagate the uncertainty in observations due to sampling variability in the groups. This effectively assumes that small and large groups are estimated equally well, and should have equal weight in the final regression. We account for this uncertainty with a Bayesian distribution regression formalism, improving the robustness and performance of the model when group sizes vary. We frame our models in a neural network style, allowing for simple MAP inference using backpropagation to learn the parameters, as well as MCMC-based inference which can fully propagate uncertainty. We demonstrate our approach on illustrative toy datasets, as well as on a challenging problem of predicting age from images.


title: Asymptotics of Nonlinear LSE Precoders with Applications to Transmit Antenna Selection

abstract: This paper studies the large-system performance of Least Square Error (LSE) precoders which~minimize~the~input-output distortion over an arbitrary support subject to a general penalty function. The asymptotics are determined via the replica method in a general form which encloses the Replica Symmetric (RS) and Replica Symmetry Breaking (RSB) ans\"atze. As a result, the "marginal decoupling property" of LSE precoders for $b$-steps of RSB is derived. The generality of the studied setup enables us to address special cases in which the number of active transmit antennas are constrained. Our numerical investigations depict that the computationally efficient forms of LSE precoders based on "$\ell_1$-norm" minimization perform close to the cases with "zero-norm" penalty function which have a considerable improvements compared to the random antenna selection. For the case with BPSK signals and restricted number of active antennas, the results show that RS fails to predict the performance while the RSB ansatz is consistent with theoretical bounds.


title: Linear dichroism in angle-resolved core-level photoemission spectra reflecting 4f ground-state symmetry of strongly correlated cubic Pr compounds

abstract: We report experimentally observed linear dichroism in angle-resolved core-level photoemission spectra of PrIr2Zn20 and PrB6 in cubic symmetry. The different anisotropic 4f charge distributions between the compounds due to the crystalline-electric-field splitting are responsible for the difference in the linear dichroism, which has been verified by spectral simulations with the full multiplet theory for a single-site Pr3+ ion in cubic symmetry. The observed linear dichroism and polarization-dependent spectra in two different photoelectron directions for PrIr2Zn20 are reproduced by theoretical analysis for the Gamma_3 ground state, whereas those of the Pr 3d and 4d core levels indicate the Gamma_5 ground state for PrB6.


title: Displayed Categories

abstract: We introduce and develop the notion of *displayed categories*. A displayed category over a category C is equivalent to "a category D and functor F : D --> C", but instead of having a single collection of "objects of D" with a map to the objects of C, the objects are given as a family indexed by objects of C, and similarly for the morphisms. This encapsulates a common way of building categories in practice, by starting with an existing category and adding extra data/properties to the objects and morphisms. The interest of this seemingly trivial reformulation is that various properties of functors are more naturally defined as properties of the corresponding displayed categories. Grothendieck fibrations, for example, when defined as certain functors, use equality on objects in their definition. When defined instead as certain displayed categories, no reference to equality on objects is required. Moreover, almost all examples of fibrations in nature are, in fact, categories whose standard construction can be seen as going via displayed categories. We therefore propose displayed categories as a basis for the development of fibrations in the type-theoretic setting, and similarly for various other notions whose classical definitions involve equality on objects. Besides giving a conceptual clarification of such issues, displayed categories also provide a powerful tool in computer formalisation, unifying and abstracting common constructions and proof techniques of category theory, and enabling modular reasoning about categories of multi-component structures. As such, most of the material of this article has been formalised in Coq over the UniMath library, with the aim of providing a practical library for use in further developments.


title: A Derived Decomposition for Equivariant $D$-modules

abstract: We show that the $G$-equivariant coherent derived category of $D$-modules on $\mathfrak{g}$ admits an orthogonal decomposition in to blocks indexed by cuspidal data (in the sense of Lusztig). Each block admits a monadic description in terms a certain differential graded algebra related to the homology of Steinberg varieties, which resembles a "triple affine" Hecke algebra. Our results generalize the work of Rider and Rider--Russell on constructible complexes on the nilpotent cone, and the earlier work of the author on the abelian category of equivariant $D$-modules on $\mathfrak{g}$. However, the algebra controlling the entire derived category of $D$-modules appears to be substantially more complicated than either of these special cases, as evidenced by the non-splitting of the Mackey filtration on the monad controlling each block.


title: Iterative wave-front reconstruction in the Fourier domain

abstract: The use of Fourier methods in wave-front reconstruction can significantly reduce the computation time for large telescopes with a high number of degrees of freedom. However, Fourier algorithms for discrete data require a rectangular data set which conform to specific boundary requirements, whereas wave-front sensor data is typically defined over a circular domain (the telescope pupil). Here we present an iterative Gerchberg routine modified for the purposes of discrete wave-front reconstruction which adapts the measurement data (wave-front sensor slopes) for Fourier analysis, fulfilling the requirements of the Fast Fourier Transform (FFT) and providing accurate reconstruction. The routine is used in the adaptation step only and can be coupled to any other Wiener-like or least-squares method. We compare simulations using this method with previous Fourier methods and show an increase in performance in terms of Strehl ratio and a reduction in noise propagation for a 40x40 SPHERE-like adaptive optics system. For closed loop operation with minimal iterations the Gerchberg method provides an improvement in Strehl, from 95.4% to 96.9% in K-band. This corresponds to ~40 nm improvement in rms, and avoids the high spatial frequency errors present in other methods, providing an increase in contrast towards the edge of the correctable band.


title: Maximum principle for a stochastic delayed system involving terminal state constraints

abstract: We investigate a stochastic optimal control problem where the controlled system is depicted as a stochastic differential delayed equation; however, at the terminal time, the state is constrained in a convex set. We firstly introduce an equivalent backward delayed system depicted as a time-delayed backward stochastic differential equation. Then a stochastic maximum principle is obtained by virtue of Ekeland's variational principle. Finally, applications to a state constrained stochastic delayed linear-quadratic control model and a production-consumption choice problem are studied to illustrate the main obtained result.


title: Challenges in Monocular Visual Odometry: Photometric Calibration, Motion Bias and Rolling Shutter Effect

abstract: Monocular visual odometry (VO) and simultaneous localization and mapping (SLAM) have seen tremendous improvements in accuracy, robustness and efficiency, and have gained increasing popularity over recent years. Nevertheless, not so many discussions have been carried out to reveal the influences of three very influential yet easily overlooked aspects: photometric calibration, motion bias and rolling shutter effect. In this work, we evaluate these three aspects quantitatively on the state of the art of direct, feature-based and semi-direct methods, providing the community with useful practical knowledge both for better applying existing methods and developing new algorithms of VO and SLAM. Conclusions (some of which are counter-intuitive) are drawn with both technical and empirical analyses to all of our experiments. Possible improvements on existing methods are directed or proposed, such as a sub-pixel accuracy refinement of ORB-SLAM which boosts its performance.


title: A Feature Embedding Strategy for High-level CNN representations from Multiple ConvNets

abstract: Following the rapidly growing digital image usage, automatic image categorization has become preeminent research area. It has broaden and adopted many algorithms from time to time, whereby multi-feature (generally, hand-engineered features) based image characterization comes handy to improve accuracy. Recently, in machine learning, pre-trained deep convolutional neural networks (DCNNs or ConvNets) have been that the features extracted through such DCNN can improve classification accuracy. Thence, in this paper, we further investigate a feature embedding strategy to exploit cues from multiple DCNNs. We derive a generalized feature space by embedding three different DCNN bottleneck features with weights respect to their Softmax cross-entropy loss. Test outcomes on six different object classification data-sets and an action classification data-set show that regardless of variation in image statistics and tasks the proposed multi-DCNN bottleneck feature fusion is well suited to image classification tasks and an effective complement of DCNN. The comparisons to existing fusion-based image classification approaches prove that the proposed method surmounts the state-of-the-art methods and produces competitive results with fully trained DCNNs as well.


title: Towards Verifying Nonlinear Integer Arithmetic

abstract: We eliminate a key roadblock to efficient verification of nonlinear integer arithmetic using CDCL SAT solvers, by showing how to construct short resolution proofs for many properties of the most widely used multiplier circuits. Such short proofs were conjectured not to exist. More precisely, we give n^{O(1)} size regular resolution proofs for arbitrary degree 2 identities on array, diagonal, and Booth multipliers and quasipolynomial- n^{O(\log n)} size proofs for these identities on Wallace tree multipliers.


title: On a product of certain primes

abstract: We study the properties of the product, which runs over the primes, $$\mathfrak{p}_n = \prod_{s_p(n) \, \geq \, p} p \quad (n \geq 1),$$ where $s_p(n)$ denotes the sum of the base-$p$ digits of $n$. One important property is the fact that $\mathfrak{p}_n$ equals the denominator of the Bernoulli polynomial $B_n(x) - B_n$, where we provide a short $p$-adic proof. Moreover, we consider the decomposition $\mathfrak{p}_n = \mathfrak{p}_n^- \cdot \mathfrak{p}_n^+$, where $\mathfrak{p}_n^+$ contains only those primes $p > \sqrt{n}$. Let $\omega( \cdot )$ denote the number of prime divisors. We show that $\omega( \mathfrak{p}_n^+ ) < \sqrt{n}$, while we raise the explicit conjecture that $$\omega( \mathfrak{p}_n^+ ) \, \sim \, \kappa \, \frac{\sqrt{n}}{\log n} \quad \text{as $n \to \infty$}$$ with a certain constant $\kappa > 1$, supported by several computations.


title: A Deep Reinforced Model for Abstractive Summarization

abstract: Attentional, RNN-based encoder-decoder models for abstractive summarization have achieved good performance on short input and output sequences. For longer documents and summaries however these models often include repetitive and incoherent phrases. We introduce a neural network model with a novel intra-attention that attends over the input and continuously generated output separately, and a new training method that combines standard supervised word prediction and reinforcement learning (RL). Models trained only with supervised learning often exhibit "exposure bias" - they assume ground truth is provided at each step during training. However, when standard word prediction is combined with the global sequence prediction training of RL the resulting summaries become more readable. We evaluate this model on the CNN/Daily Mail and New York Times datasets. Our model obtains a 41.16 ROUGE-1 score on the CNN/Daily Mail dataset, an improvement over previous state-of-the-art models. Human evaluation also shows that our model produces higher quality summaries.


title: Irreversible proliferation of magnetic moments at cleaved surfaces of the topological Kondo insulator SmB6

abstract: The compound SmB$_6$ is the best established realization of a topological Kondo insulator, in which a topological insulator state is obtained through Kondo coherence. Recent studies have found evidence that the surface of SmB$_6$ hosts ferromagnetic domains, creating an intrinsic platform for unidirectional ballistic transport at the domain boundaries. Here, surface-sensitive X-ray absorption (XAS) and bulk-sensitive resonant inelastic X-ray scattering (RIXS) spectra are measured at the Sm N$_{4,5}$-edge, and used to evaluate electronic symmetries, excitations and temperature dependence near the surface of cleaved samples. The XAS data show that the density of large-moment atomic multiplet states on a cleaved surface grows irreversibly over time, to a degree that likely exceeds a related change that has recently been observed in the surface 4f orbital occupation.


title: Characterization of the Lp-range of the Poisson transform on the Octonionic Hyperbolic Plane

abstract: We give an image characterization of the Poisson transform of L-p functions on the boundary of the exceptional symmetric space by using the Hardy-type norm.


title: Modeling observers as physical systems representing the world from within: Quantum theory as a physical and self-referential theory of inference

abstract: In 1929 Szilard pointed out that the physics of the observer may play a role in the analysis of experiments. The same year, Bohr pointed out that complementarity appears to arise naturally in psychology where both the objects of perception and the perceiving subject belong to 'our mental content'. Here we argue that the formalism of quantum theory can be derived from two related intuitive principles: (i) inference is a physical process performed by physical systems, observers, which are part of the experimental setup---this implies non-commutativity and imaginary-time quantum mechanics; (ii) experiments must be described from a first-person perspective---this leads to self-reference, complementarity, and real-time quantum dynamics. This approach sheds new light on the foundations of quantum theory and suggests fundamental equations in physics are typically of second order due to the physical nature of the observer. It also suggests some experimental conjectures: (i) the quantum of action could be understood as the result of the additional energy required to transition from unconscious to conscious perception; (ii) humans can observe a single photon of visible light; (iii) self-aware systems and the neural correlates of the self should be composed of two complementary sub-systems, much like the DNA molecule is composed of two strands---this may help explain the double-hemisphere architecture of the brain. Moreover, this approach may help bridge the gap between science and human experience. We discuss the potential implications of these ideas for the modern research programs on consciousness and contemplative science. As side results: (i) we show that message-passing algorithms and stochastic processes can be written in a quantum-like manner; (ii) we provide evidence that non-stoquasticity, a quantum computational resource, may be related to non-equilibrium phenomena.


title: A Sobering Assessment of Small-Molecule Force Field Methods for Low Energy Conformer Predictions

abstract: We have carried out a large scale computational investigation to assess the utility of common small-molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely-used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn-Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening.


title: Quantum Emulation of Extreme Non-equilibrium Phenomena with Trapped Atoms

abstract: Ultracold atomic physics experiments offer a nearly ideal context for the investigation of quantum systems far from equilibrium. We describe three related emerging directions of research into extreme non-equilibrium phenomena in atom traps: quantum emulation of ultrafast atom-light interactions, coherent phasonic spectroscopy in tunable quasicrystals, and realization of Floquet matter in strongly-driven lattice systems. We show that all three should enable quantum emulation in parameter regimes inaccessible in solid-state experiments, facilitating a complementary approach to open problems in non-equilibrium condensed matter.


title: Categorical structures for type theory in univalent foundations

abstract: In this paper, we analyze and compare three of the many algebraic structures that have been used for modeling dependent type theories: categories with families, split type-categories, and representable maps of presheaves. We study these in univalent type theory, where the comparisons between them can be given more elementarily than in set-theoretic foundations. Specifically, we construct maps between the various types of structures, and show that assuming the Univalence axiom, some of the comparisons are equivalences. We then analyze how these structures transfer along (weak and strong) equivalences of categories, and, in particular, show how they descend from a category (not assumed univalent/saturated) to its Rezk completion. To this end, we introduce relative universes, generalizing the preceding notions, and study the transfer of such relative universes along suitable structure. We work throughout in (intensional) dependent type theory; some results, but not all, assume the univalence axiom. All the material of this paper has been formalized in Coq, over the UniMath library.


title: Infrared Divergences in QED, Revisited

abstract: Recently it has been shown that the vacuum state in QED is infinitely degenerate. Moreover a transition among the degenerate vacua is induced in any nontrivial scattering process and determined from the associated soft factor. Conventional computations of scattering amplitudes in QED do not account for this vacuum degeneracy and therefore always give zero. This vanishing of all conventional QED amplitudes is usually attributed to infrared divergences. Here we show that if these vacuum transitions are properly accounted for, the resulting amplitudes are nonzero and infrared finite. Our construction of finite amplitudes is mathematically equivalent to, and amounts to a physical reinterpretation of, the 1970 construction of Faddeev and Kulish.


title: FDR-Corrected Sparse Canonical Correlation Analysis with Applications to Imaging Genomics

abstract: Reducing the number of false discoveries is presently one of the most pressing issues in the life sciences. It is of especially great importance for many applications in neuroimaging and genomics, where datasets are typically high-dimensional, which means that the number of explanatory variables exceeds the sample size. The false discovery rate (FDR) is a criterion that can be employed to address that issue. Thus it has gained great popularity as a tool for testing multiple hypotheses. Canonical correlation analysis (CCA) is a statistical technique that is used to make sense of the cross-correlation of two sets of measurements collected on the same set of samples (e.g., brain imaging and genomic data for the same mental illness patients), and sparse CCA extends the classical method to high-dimensional settings. Here we propose a way of applying the FDR concept to sparse CCA, and a method to control the FDR. The proposed FDR correction directly influences the sparsity of the solution, adapting it to the unknown true sparsity level. Theoretical derivation as well as simulation studies show that our procedure indeed keeps the FDR of the canonical vectors below a user-specified target level. We apply the proposed method to an imaging genomics dataset from the Philadelphia Neurodevelopmental Cohort. Our results link the brain connectivity profiles derived from brain activity during an emotion identification task, as measured by functional magnetic resonance imaging (fMRI), to the corresponding subjects' genomic data.


title: Extended Red Emission in IC59 and IC63

abstract: We analysed new wide-field, wide- and narrow-band optical images of IC 59 and IC 63, two nebulae which are externally illuminated by the early B-star {\gamma} Cas, with the objective of mapping the extended red emission (ERE), a dust-related photoluminescence process that is still poorly understood, in these two clouds. The spatial distribution of the ERE relative to the direction of the incident radiation and relative to other emission processes, whose carriers and excitation requirements are known, provides important constraints on the excitation of the ERE. In both nebulae, we find the ERE intensity to peak spatially well before the more extended distribution of mid-infrared emission in the unidentified infrared bands, supporting earlier findings that point toward far-ultraviolet (11 eV < E$_\mathrm{{photon}}$ < 13.6 eV) photons as the source of ERE excitation. The band-integrated absolute intensities of the ERE in IC 59 and IC 63 measured relative to the number density of photons available for ERE excitation are lower by about two orders of magnitude compared to ERE intensities observed in the high-latitude diffuse interstellar medium (ISM). This suggests that the lifetime of the ERE carriers is significantly reduced in the more intense radiation field prevailing in IC 59 and IC 63, pointing toward potential carriers that are only marginally stable against photo-processing under interstellar conditions. A model involving isolated molecules or molecular ions, capable of inverse internal conversion and recurrent fluorescence, appears to provide the most likely explanation for our observational results.


title: Compressible-incompressible two-phase flows with phase transition: model problem

abstract: We study the compressible and incompressible two-phase flows separated by a sharp interface with a phase transition and a surface tension. In particular, we consider the problem in $\mathbb{R}^N$, and the Navier-Stokes-Korteweg equations is used in the upper domain and the Navier-Stokes equations is used in the lower domain. We prove the existence of $\mathcal{R}$-bounded solution operator families for a resolvent problem arising from its model problem. According to Shibata \cite{GS2014}, the regularity of $\rho_+$ is $W^1_q$ in space, but to solve the kinetic equation: $\mathbf{u}_\Gamma\cdot\mathbf{n}_t = [[\rho\mathbf{u}]]\cdot\mathbf{n}_t /[[\rho]]$ on $\Gamma_t$ we need $W^{2-1/q}_q$ regularity of $\rho_+$ on $\Gamma_t$, which means the regularity loss. Since the regularity of $\rho_+$ dominated by the Navier-Stokes-Korteweg equations is $W^3_q$ in space, we eliminate the problem by using the Navier-Stokes-Korteweg equations instead of the compressible Navier-Stokes equations.


title: Positive Definiteness of Paired Symmetric Tensors and Elasticity Tensors

abstract: In this paper, we consider higher order paired symmetric tensors and strongly paired symmetric tensors. Elasticity tensors and higher order elasticity tensors in solid mechanics are strongly paired symmetric tensors. A (strongly) paired symmetric tensor is said to be positive definite if the homogeneous polynomial defined by it is positive definite. Positive definiteness of elasticity and higher order elasticity tensors is strong ellipticity in solid mechanics, which plays an important role in nonlinear elasticity theory. We mainly investigate positive definiteness of fourth order three dimensional and sixth order three dimensional (strongly) paired symmetric tensors. We first show that the concerned (strongly) paired symmetric tensor is positive definite if and only if its smallest $M$-eigenvalue is positive. Second, we propose several necessary and sufficient conditions under which the concerned (strongly) paired symmetric tensor is positive definite. Third, we study the conditions under which the homogeneous polynomial defined by a fourth order three dimensional or sixth order three dimensional (strongly) paired symmetric tensor can be written as a sum of squares of polynomials, and further, propose several necessary and/or sufficient conditions to judge whether the concerned (strongly) paired symmetric tensors are positive definite or not. Fourth, by using semidefinite relaxation we propose a sequential semidefinite programming method to compute the smallest $M$-eigenvalue of a fourth order three dimensional (strongly) paired symmetric tensor, by which we can check positive definiteness of the concerned tensor. The preliminary numerical results demonstrate that our method is effective.


title: Summing coincidence in rare event gamma-ray measurements under an ultra-low background environment

abstract: A Monte Carlo method based on the GEANT4 toolkit has been developed to correct the full-energy peak (FEP) efficiencies of a high purity germanium (HPGe) detector equipped with a low background shielding system, and moreover evaluated using summing peaks in a numerical way. It is found that the FEP efficiencies of $^{60}$Co, $^{133}$Ba and $^{152}$Eu can be improved up to 18\% by taking the calculated true summing \mbox{coincidence} factors (TSCFs) correction into account. Counts of summing coincidence $\gamma$ peaks in the spectrum of $^{152}$Eu can be well reproduced using the corrected efficiency curve within an accuracy of 3\%.


title: Management system for the SND experiments

abstract: A new management system for the SND detector experiments (at VEPP-2000 collider in Novosibirsk) is developed. We describe here the interaction between a user and the SND databases. These databases contain experiment configuration, conditions and metadata. The new system is designed in client-server architecture. It has several logical layers corresponding to the users roles. A new template engine is created. A web application is implemented using Node.js framework. At the time the application provides: showing and editing configuration; showing experiment metadata and experiment conditions data index; showing SND log (prototype).


title: On the generalization of Erd\H{o}s-Vincze's theorem about the approximation of closed convex plane curves by polyellipses

abstract: A polyellipse is a curve in the Euclidean plane all of whose points have the same sum of distances from finitely many given points (focuses). The classical version of Erd\H{o}s-Vincze's theorem states that regular triangles can not be presented as the Hausdorff limit of polyellipses even if the number of the focuses can be arbitrary large. In other words the topological closure of the set of polyellipses with respect to the Hausdorff distance does not contain any regular triangle and we have a negative answer to the problem posed by E. V\'azsonyi (Weissfeld) about the approximation of closed convex plane curves by polyellipses. It is the additive version of the approximation of simple closed plane curves by polynomial lemniscates all of whose points have the same product of distances from finitely many given points (focuses). Here we are going to generalize the classical version of Erd\H{o}s-Vincze's theorem for regular polygons in the plane. We will conclude that the error of the approximation tends to zero as the number of the vertices of the regular polygon tends to the infinity. The decreasing tendency of the approximation error gives the idea to construct curves in the topological closure of the set of polyellipses. If we use integration to compute the average distance of a point from a given (focal) set in the plane then the curves all of whose points have the same average distance from the focal set can be given as the Hausdorff limit of polyellipses corresponding to partial sums.


title: Anomalous slowing down of individual human activity due to successive decision-making processes

abstract: Motivated by a host of empirical evidences revealing the bursty character of human dynamics, we develop a model of human activity based on successive switching between an hesitation state and a decision-realization state, with residency times in the hesitation state distributed according to a heavy-tailed Pareto distribution. This model is particularly reminiscent of an individual strolling through a randomly distributed human crowd. Using a stochastic model based on the concept of anomalous and non-Markovian L\'evy walk, we show exactly that successive decision-making processes drastically slow down the progression of an individual faced with randomly distributed obstacles. Specifically, we prove exactly that the average displacement exhibits a sublinear scaling with time that finds its origins in: (i) the intrinsically non-Markovian character of human activity, and (ii) the power law distribution of hesitation times.


title: An Extension of Heron's Formula

abstract: This paper introduces an extension of Heron's formula to approximate area of cyclic n-gons where the error never exceeds $\frac{\pi}{e}-1$


title: CMS Hardware Track Trigger: New Opportunities for Long-Lived Particle Searches at the HL-LHC

abstract: The planned upgrade of the CMS detector for the High Luminosity LHC allows to find tracks in the silicon tracker for every single LHC collision and use them in the first level (hardware) trigger decision. So far, studies by CMS collaboration concentrated on the maintaining the overall trigger performance in the punishing pile up environment. We argue that the potential capabilities of the track trigger are much wider, and may offer groundbreaking opportunities for new physics searches. As an example, and to facilitate community discussion, we use a simple toy simulation to study rare Higgs decays into new particles with lifetime of order of a few mm.


title: [Fe XIV] and [Fe XI] reveal the forward shock in SNR 1E0102.2-7219

abstract: Aims. We study the forward shock in the oxygen-rich young supernova remnant (SNR) 1E0102.2-7219 (1E0102 in short) via optical coronal emission from [Fe XIV] and [Fe XI]: emission lines which offer an alternative method to X-rays to do so. Methods. We have used the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph at the Very Large Telescope (VLT) on Cerro Paranal to obtain deep observations of SNR 1E0102 in the Small Magellanic Cloud. Our observations cover the entire extent of the remnant with a seeing limited spatial resolution of 0.7" = 0.2 pc at the distance of 1E 0102. Results. Our MUSE observations unambiguously reveal the presence of [Fe XIV] and [Fe XI] emission in 1E0102. The emission largely arises from a thin, partial ring of filaments surrounding the fast moving O-rich ejecta in the system. The brightest [Fe XIV] and [Fe XI] emission is found along the eastern and north-western sides of 1E0102, where shocks are driven into denser ISM material, while fainter emission along the northern edge reveals the location of the forward shock in lower density gas, possibly the relic stellar wind cavity. Modeling of the eastern shocks and the photoionization precursor surrounding 1E0102, we derive a pre-shock density $n_H$ = (7.4 +-1.5) cm$^{-3}$, and a shock velocity 330 km/s < $v_s$ < 350 km/s.


title: Survival of the obscuring torus in the most powerful active galactic nuclei

abstract: Dedicated searches generally find a decreasing fraction of obscured Active Galactic Nuclei (AGN) with increasing AGN luminosity. This has often been interpreted as evidence for a decrease of the covering factor of the AGN torus with increasing luminosity, the so-called receding torus models. Using a complete flux-limited X-ray selected sample of 199 AGN, from the Bright Ultra-hard XMM-Newton Survey, we determine the intrinsic fraction of optical type-2 AGN at 0.05$\leq$z$\leq$1 as a function of rest-frame 2-10 keV X-ray luminosity from 10$^{42}$ to 10$^{45}$ erg/s. We use the distributions of covering factors of AGN tori derived from CLUMPY torus models. Since these distributions combined over the total AGN population need to match the intrinsic type-2 AGN fraction, we reveal a population of X-ray undetected objects with high-covering factor tori, which are increasingly numerous at higher AGN luminosities. When these "missing" objects are included, we find that Compton-thick AGN account at most for 37$_{-10}^{+9}$% of the total population. The intrinsic type-2 AGN fraction is 58$\pm$4% and has a weak, non-significant (less than 2$\sigma$) luminosity dependence. This contradicts the results generally reported by AGN surveys, and the expectations from receding torus models. Our findings imply that the majority of luminous rapidly-accreting supermassive black holes at z<1 reside in highly-obscured nuclear environments but most of them are so deeply embedded that they have so far escaped detection in X-rays in <10 keV wide-area surveys.


title: The discovery of tidal tails around the globular cluster NGC 7492 with Pan-STARRS1

abstract: We report the discovery of tidal tails around the Galactic globular cluster NGC 7492, based on the Data Release 1 of the Pan-STARRS 1 survey. The tails were detected with a version of the matched filter technique applied to the $(g-r,r)$ and $(g-i,i)$ color-magnitude diagrams. Tidal tails emerging from the cluster extend at least $\sim$3.5 degrees in the North-East to South-East direction, equivalent to $\sim1.5$ kpc in projected length.


title: Hydrodynamic charge and heat transport on inhomogeneous curved spaces

abstract: We develop the theory of hydrodynamic charge and heat transport in strongly interacting quasi-relativistic systems on manifolds with inhomogeneous spatial curvature. In solid-state physics, this is analogous to strain disorder in the underlying lattice. In the hydrodynamic limit, we find that the thermal and electrical conductivities are dominated by viscous effects, and that the thermal conductivity is most sensitive to this disorder. We compare the effects of inhomogeneity in the spatial metric to inhomogeneity in the chemical potential, and discuss the extent to which our hydrodynamic theory is relevant for experimentally realizable condensed matter systems, including suspended graphene at the Dirac point.


title: A Herschel/PACS Far Infrared Line Emission Survey of Local Luminous Infrared Galaxies

abstract: We present an analysis of [OI]63, [OIII]88, [NII]122 and [CII]158 far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ~240 local luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). We find pronounced declines -deficits- of line-to-FIR-continuum emission for [NII]122, [OI]63 and [CII]158 as a function of FIR color and infrared luminosity surface density, $\Sigma_{\rm IR}$. The median electron density of the ionized gas in LIRGs, based on the [NII]122/[NII]205 ratio, is $n_{\rm e}$ = 41 cm$^{-3}$. We find that the dispersion in the [CII]158 deficit of LIRGs is attributed to a varying fractional contribution of photo-dissociation-regions (PDRs) to the observed [CII]158 emission, f([CII]PDR) = [CII]PDR/[CII], which increases from ~60% to ~95% in the warmest LIRGs. The [OI]63/[CII]158PDR ratio is tightly correlated with the PDR gas kinetic temperature in sources where [OI]63 is not optically-thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, $n_{\rm H}$, and intensity of the interstellar radiation field, in units of G$_0$, and find G$_0$/$n_{\rm H}$ ratios ~0.1-50 cm$^3$, with ULIRGs populating the upper end of the distribution. There is a relation between G$_0$/$n_{\rm H}$ and $\Sigma_{\rm IR}$, showing a critical break at $\Sigma_{\rm IR}^{\star}$ ~ 5 x 10$^{10}$ Lsun/kpc$^2$. Below $\Sigma_{\rm IR}^{\star}$, G$_0$/$n_{\rm H}$ remains constant, ~0.32 cm$^3$, and variations in $\Sigma_{\rm IR}$ are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above $\Sigma_{\rm IR}^{\star}$, G$_0$/$n_{\rm H}$ increases rapidly with $\Sigma_{\rm IR}$, signaling a departure from the typical PDR conditions found in normal star-forming galaxies towards more intense/harder radiation fields and compact geometries typical of starbursting sources.


title: The Immitigable Nature of Assembly Bias: The Impact of Halo Definition on Assembly Bias

abstract: Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.


title: The Weak Gravity Conjecture and Scalar Fields

abstract: We propose a generalisation of the Weak Gravity Conjecture in the presence of scalar fields. The proposal is guided by properties of extremal black holes in ${\cal N}=2$ supergravity, but can be understood more generally in terms of forbidding towers of stable gravitationally bound states. It amounts to the statement that there must exist a particle on which the gauge force acts more strongly than gravity and the scalar forces combined. We also propose that the scalar force itself should act on this particle stronger than gravity. This implies that generically the mass of this particle decreases exponentially as a function of the scalar field expectation value for super-Planckian variations, which is behaviour predicted by the Refined Swampland Conjecture. In the context of ${\cal N}=2$ supergravity the Weak Gravity Conjecture bound can be tied to bounds on scalar field distances in field space. Guided by this, we present a general proof that for any linear combination of moduli in any Calabi-Yau compactification of string theory the proper field distance grows at best logarithmically with the moduli values for super-Planckian distances.


title: Next Generation Virgo Cluster Survey. XXI. The weak lensing masses of the CFHTLS and NGVS RedGOLD galaxy clusters and calibration of the optical richness

abstract: We measured stacked weak lensing cluster masses for a sample of 1325 galaxy clusters detected by the RedGOLD algorithm in the Canada-France-Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection of our best richness subsample ($20<\lambda<50$), this is the most comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure optical cluster catalogue in this redshift range, with a total of 346 clusters in $\sim164~deg^2$. We test three different mass models, and our best model includes a basic halo model, with a Navarro Frenk and White profile, and correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass-richness relation. With this model, we obtain a mass-richness relation of $\log{M_{\rm 200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical uncertainties). This result is consistent with other published lensing mass-richness relations. When compared to X-ray masses and mass proxies, we find that on average weak lensing masses are $\sim 10\%$ higher than those derived in the X-ray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm 200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and simulations. We also give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, $L_X$ and $T_X$, and compare them with previous results.


title: M-theoretic Lichnerowicz formula and supersymmetry

abstract: A suitable generalisation of the Lichnerowicz formula can relate the squares of supersymmetric operators to the effective action, the Bianchi identities for fluxes, and some equations of motion. Recently, such formulae have also been shown to underlie the (generalised) geometry of supersymmetric theories. In this paper, we derive an M-theoretic Lichnerowicz formula that describes eleven-dimensional supergravity together with its higher-derivative couplings. The first corrections to the action appear at eight-derivative level, and the construction yields two different supersymmetric invariants, each with a free coefficient. We discuss the restriction of our construction to seven-dimensional internal spaces, and implications for compactifications on manifolds of $G_2$ holonomy. Inclusion of fluxes and computation of contributions with higher than eight derivatives are also discussed.


title: The RedGOLD Cluster Detection Algorithm and its Cluster Candidate Catalogue for the CFHT-LS W1

abstract: RedGOLD searches for red-sequence galaxy overdensities while minimizing contamination from dusty star-forming galaxies. It imposes an NFW profile and calculates cluster detection significance and richness. We optimize these latter two parameters using both simulations and X-ray detected cluster catalogs, and obtain a catalog $\sim 80\%$ pure up to $z \sim 1$, and $\sim 100\%$ ($\sim 70\%$) complete at $z\le 0.6$ ( $z\lesssim1$) for galaxy clusters with $M \gtrsim 10^{14}\ {\rm M_{\odot}}$ at the CFHT-LS Wide depth. In the CFHT-LS W1, we detect 11 cluster candidates per $\rm deg^2$ out to $z\sim1.1$. When we optimize both completeness and purity, RedGOLD obtains a cluster catalog with higher completeness and purity than other public catalogs, obtained using CFHT-LS W1 observations, for $M \gtrsim 10^{14}\ {\rm M_{\odot}}$. We use X-ray detected cluster samples to extend the study of the X-ray temperature-optical richness relation to a lower mass threshold, and find a mass scatter at fixed richness of $\sigma_{lnM|\lambda}=0.39\pm0.07$ and $\sigma_{lnM|\lambda}=0.30\pm0.13$ for the Gozaliasl et al. (2014) and Mehrtens et al. (2012) samples. When considering similar mass ranges as previous work, we recover a smaller scatter in mass at fixed richness. We recover $93\%$ of the redMaPPer detections, and find that its richness estimates is on average $\sim 40-50\%$ larger than ours at $z>0.3$. RedGOLD recovers X-ray cluster spectroscopic redshifts at better than $5\%$ up to $z\sim1$, and the centers within a few tens of arcseconds.


title: The Next Generation Virgo Cluster Survey. XX. RedGOLD Background Galaxy Cluster Detections

abstract: We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.


title: Detection of the Yarkovsky effect for C-type asteroids in the Veritas family

abstract: The age of a young asteroid family can be determined by tracking the orbits of family members backward in time and showing that they converge at some time in the past. Here we consider the Veritas family. We find that the membership of the Veritas family increased enormously since the last detailed analysis of the family. Using backward integration, we confirm the convergence of nodal longitudes $\Omega$, and, for the first time, also obtain a simultaneous convergence of pericenter longitudes $\varpi$. The Veritas family is found to be $8.23^{+0.37}_{-0.31}$~Myr old. To obtain a tight convergence of $\Omega$ and $\varpi$, as expected from low ejection speeds of fragments, the Yarkovsky effect needs to be included in the modeling of the past orbital histories of Veritas family members. Using this method, we compute the Yarkovsky semi-major axis drift rates, ${\rm d}a/{\rm d}t$, for 274 member asteroids. The distribution of ${\rm d}a/{\rm d}t$ values is consistent with a population of C-type objects with low densities and low thermal conductivities. The accuracy of individual ${\rm d}a/{\rm d}t$ measurements is limited by the effect of close encounters of member asteroids to (1) Ceres and other massive asteroids, which cannot be evaluated with confidence.


title: Holography and $AdS_2$ gravity with a dynamical aether

abstract: We study two-dimensional Einstein-aether (or equivalently Ho\v{r}ava-Lifshitz) gravity, which has an $AdS_2$ solution. We examine various properties of this solution in the context of holography. We first show that the asymptotic symmetry group is the full set of time reparametrizations, the one-dimensional conformal group. At the same time there are configurations with finite energy and temperature, which indicate a violation of the Ward identity associated with one-dimensional conformal invariance. These solutions are characterized by a universal causal horizon and we show that the associated entropy of the universal horizon scales with the logarithm of the temperature. We discuss the puzzles associated with this result and argue that the violation of the Ward identity is associated with a type of explicit breaking of time reparametrizations in the hypothetical $0+1$ dimensional dual system.


title: Quantum and private capacities of low-noise channels

abstract: We determine both the quantum and the private capacities of low-noise quantum channels to leading orders in the channel's distance to the perfect channel. It has been an open problem for more than 20 years to determine the capacities of some of these low-noise channels such as the depolarizing channel. We also show that both capacities are equal to the single-letter coherent information of the channel, again to leading orders. We thus find that, in the low noise regime, super-additivity and degenerate codes have negligible benefit for the quantum capacity, and shielding does not improve the private capacity beyond the quantum capacity, in stark contrast to the situation when noisier channels are considered.


title: An Optimal Dimensionality Multi-shell Sampling Scheme with Accurate and Efficient Transforms for Diffusion MRI

abstract: This paper proposes a multi-shell sampling scheme and corresponding transforms for the accurate reconstruction of the diffusion signal in diffusion MRI by expansion in the spherical polar Fourier (SPF) basis. The sampling scheme uses an optimal number of samples, equal to the degrees of freedom of the band-limited diffusion signal in the SPF domain, and allows for computationally efficient reconstruction. We use synthetic data sets to demonstrate that the proposed scheme allows for greater reconstruction accuracy of the diffusion signal than the multi-shell sampling schemes obtained using the generalised electrostatic energy minimisation (gEEM) method used in the Human Connectome Project. We also demonstrate that the proposed sampling scheme allows for increased angular discrimination and improved rotational invariance of reconstruction accuracy than the gEEM schemes.


title: Variational calculation of the ground state of closed-shell nuclei up to $A=40$

abstract: Variational calculations of ground-state properties of $^4$He, $^{16}$O, and $^{40}$Ca are carried out employing realistic phenomenological two- and three-nucleon potentials. The trial wave function includes two- and three-body correlations acting on a product of single-particle determinants. Expectation values are evaluated with a cluster expansion for the spin-isospin dependent correlations considering up to five-body cluster terms. The optimal wave function is obtained by minimizing the energy expectation value over a set of up to 20 parameters by means of a nonlinear optimization library. We present results for the binding energy, charge radius, one- and two-body densities, single-nucleon momentum distribution, charge form factor, and Coulomb sum rule. We find that the employed three-nucleon interaction becomes repulsive for $A\geq16$. In $^{16}$O the inclusion of such a force provides a better description of the properties of the nucleus. In $^{40}$Ca instead, the repulsive behavior of the three-body interaction fails to reproduce experimental data for the charge radius and the charge form factor. We find that the high-momentum region of the momentum distributions, determined by the short-range terms of nuclear correlations, exhibit a universal behavior independent of the particular nucleus. The comparison of the Coulomb sum rules for $^4$He, $^{16}$O, and $^{40}$Ca reported in this work will help elucidate in-medium modifications of the nucleon form factors.


title: A note on the high power Diophantine equations

abstract: In this paper, we solve the simultaneous Diophantine equations(SDE) x_1^u+...+x_n^u=k(y_1^u+...+y_{n/k}); u=1,3, where n >3, and k< n, is a divisor of n , and obtain nontrivial parametric solution for them. Furthermore we present a method for producing another solution for the above Diophantine equation (DE) for the case u = 3, when a solution is given. We work out some examples and find nontrivial parametric solutions for each case in nonzero integers. Also we prove that the other DE p_1x_1^{a_1}+....+p_nx_n^{a_n}=q_1y_1^{b_1}+...+q_my_m^{b_m} , has parametric solution and infinitely many solutions in nonzero integers with the condition that: there is a i such that p_i=1, and (a_i,a_1....a_{i-1}a_{i+1}...b_1...b_m)=1, or there is a j such that q_j=1, and (b_j,a_1...a_nb_1...b_{j-1}b_{j+1}..b_m)=1. Finally we study the DE x^a_y^b=z^c.


title: Dispersion relations for $\eta'\to\eta\pi\pi$

abstract: We present a dispersive analysis of the decay amplitude for $\eta'\to\eta\pi\pi$ that is based on the fundamental principles of analyticity and unitarity. In this framework, final-state interactions are fully taken into account. Our dispersive representation relies only on input for the $\pi\pi$ and $\pi\eta$ scattering phase shifts. Isospin symmetry allows us to describe both the charged and neutral decay channel in terms of the same function. The dispersion relation contains subtraction constants that cannot be fixed by unitarity. We determine these parameters by a fit to Dalitz-plot data from the VES and BES-III experiments. We study the prediction of a low-energy theorem and compare the dispersive fit to variants of chiral perturbation theory.


title: Preconditioning for Accurate Solutions of Linear Systems and Eigenvalue Problems

abstract: This paper develops the preconditioning technique as a method to address the accuracy issue caused by ill-conditioning. Given a preconditioner $M$ for an ill-conditioned linear system $Ax=b$, we show that, if the inverse of the preconditioner $M^{-1}$ can be applied to vectors accurately, then the linear system can be solved accurately. A stability concept called inverse-equivalent accuracy is introduced to describe higher accuracy that is achieved and an error analysis will be presented. As an application, we use the preconditioning approach to accurately compute a few smallest eigenvalues of certain ill-conditioned matrices. Numerical examples are presented to illustrate the error analysis and the performance of the methods.


title: All coherent field states entangle equally

abstract: We analyze the interactions of particle detectors with coherent states of a free scalar field. We find that the eigenvalues of the post-interaction density matrices of i) a single detector, ii) two detectors, and iii) the partial transpose of the latter, are all independent of which coherent state the field was in. A consequence of these results is that a detector pair can harvest the same amount of entanglement from any coherent field state as from the vacuum.


title: An invertibility criterion in a C*-algebra acting on the Hardy space with applications to composition operators

abstract: In this paper we prove an invertibility criterion for certain operators which is given as a linear algebraic combination of Toeplitz operators and Fourier multipliers acting on the Hardy space of the unit disc. Very similar to the case of Toeplitz operators we prove that such operators are invertible if and only if they are Fredholm and their Fredholm index is zero. As an application we prove that for "quasi-parabolic" composition operators the spectra and the essential spectra are equal.


title: Wave-particle duality employing quantum coherence in superposition with non-orthogonal pointers

abstract: We propose the notion of quantum coherence for superpositions over states which are not necessarily mutually orthogonal. This anticipatedly leads to a resource theory of non-orthogonal coherence. We characterize free states and free operations in this theory, and connect the latter with free operations in the resource theory of quantum coherence for orthogonal bases. We show that the concept of non-orthogonal coherence naturally furnishes us with a wave-particle duality in quantum double-slit experiments where the channels beyond the slits are leaky between them. Furthermore, we demonstrate existence of a unique maximally coherent qubit state corresponding to any given purity. In addition, and in contradistinction with the case of orthogonal bases, there appears a non-trivial minimally coherent qubit state for a given purity. We also study the behavior of quantum coherence for some typical configurations of non-orthogonal bases which have no analogs for orthogonal bases. We further investigate the problem of determining the energy cost of creating non-orthogonal coherence, and find that it scales linearly with the non-orthogonal coherence created.


title: A Spectral reciprocity formula and non-vanishing for L-functions on GL(4)xGL(2)

abstract: We develop a reciprocity formula for a spectral sum over central values of L-functions on GL(4)xGL(2). As an application we show that for any self-dual cusp form Pi for SL(4,Z), there exists a Maass form pi for SL(2,Z) such that L(1/2, Pi x pi) is nonvanishing. An important ingredient is a "balanced" Voronoi summation formula involving Kloosterman sums on both sides, which can also be thought of as the functional equation of a certain double Dirichlet series involving Kloosterman sums and GL(4) Hecke eigenvalues.


title: Error estimate for a homogenization problem involving the Laplace-Beltrami operator

abstract: In this paper we prove a concentration result and an error estimate for a model of heat conduction in composite materials having a microscopic structure arranged in a perodic array and thermally active membranes separating the heat conductive phases.


title: The geometric classification of Leibniz algebras

abstract: We describe all rigid algebras and all irreducible components in the variety of four dimensional Leibniz algebras $\mathfrak{Leib}_4$ over $\mathbb{C}.$ In particular, we prove that the Grunewald--O'Halloran conjecture is not valid and the Vergne conjecture is valid for $\mathfrak{Leib}_4.$


title: Averaging method for systems with separatrix crossing

abstract: The averaging method provides a powerful tool for studying evolution in near-integrable systems. Existence of separatrices in the phase space of the underlying integrable system is an obstacle for application of standard results that justify using of averaging. We establish estimates that allow to use averaging method when the underlying integrable system is a system with one rotating phase, and the evolution leads to separatrix crossings.


title: Generalized Kloosterman Sums from M2-branes

abstract: Kloosterman sums play a special role in analytic number theory, for expressing the integer Fourier coefficients of modular forms as an infinite sum of Bessel functions, also known as Rademacher formula. The generalization to vector-valued modular forms is known as generalized Kloosterman sums. In the paper arxiv:1404.0033, a remarkable connection between these arithmetic sums and quantum black hole entropy was found. Nevertheless, the computation was particular for one-eighth BPS black holes in $\mathcal{N}=8$ string theory, which have a simple counting formula. Here, we review this construction and extend it to the case of one-quarter BPS black holes in $\mathcal{N}=4$ string theory, which are counted by (mock) Jacobi forms of arbitrary index. The main result is an holographic derivation of the Kloosterman sums which includes the intricate sum over phases and depends exactly on the spectral flow sectors and the spectrum of polar states. On the microscopic side we derive an analytic formula for the Kloosterman sums valid for any index, whereas from the macroscopic side we reproduce the same formula from the M-theory path integral on $\mathbb{Z}_c$ orbifolds of $AdS_2\times S^1$. A key aspect of the derivation is the identification of the spectral flow sectors with the contribution of M2 branes wrapping cycles on the compactification manifold. After a careful treatment of the measure, the sum over orbifolds results in the sum over Bessels, in perfect agreement with the Rademacher expansion at any order in perturbation theory.


title: Chemical abundances of 1111 FGK stars from the HARPS GTO planet search program II: Cu, Zn, Sr, Y, Zr, Ba, Ce, Nd and Eu

abstract: To understand the formation and evolution of the different stellar populations within our Galaxy it is essential to combine detailed kinematical and chemical information for large samples of stars. We derive chemical abundances of Cu, Zn, Sr, Y, Zr, Ba, Ce, Nd and Eu for a large sample of more than 1000 FGK dwarf stars with high-resolution ($R \sim$\,115000) and high-quality spectra from the HARPS-GTO program. The abundances are derived by a standard Local Thermodinamyc Equilibrium (LTE) analysis using measured Equivalent Widths (EWs) injected to the code MOOG and a grid of Kurucz ATLAS9 atmospheres. We find that thick disk stars are chemically disjunct for Zn and Eu and also show on average higher Zr but lower Ba and Y when compared to the thin disk stars. We also discovered that the previously identified high-$\alpha$ metal-rich population is also enhanced in Cu, Zn, Nd and Eu with respect to the thin disk but presents Ba and Y abundances lower on average, following the trend of thick disk stars towards higher metallities and further supporting the different chemical composition of this population. The ratio of heavy-s to light-s elements of thin disk stars presents the expected behaviour (increasing towards lower metallicities) and can be explained by a major contribution of low-mass AGB stars for s-process production at disk metallicities. However, the opposite trend found for thick disk stars suggests that intermediate-mass AGB stars played an important role in the enrichment of the gas from where these stars formed. Previous works in the literature also point to a possible primary production of light-s elements at low metallicities to explain this trend. Finally, we also find an enhancement of light-s elements in the thin disk at super solar metallicities which could be caused by the contribution of metal-rich AGB stars. (short version)


title: Imagination improves Multimodal Translation

abstract: We decompose multimodal translation into two sub-tasks: learning to translate and learning visually grounded representations. In a multitask learning framework, translations are learned in an attention-based encoder-decoder, and grounded representations are learned through image representation prediction. Our approach improves translation performance compared to the state of the art on the Multi30K dataset. Furthermore, it is equally effective if we train the image prediction task on the external MS COCO dataset, and we find improvements if we train the translation model on the external News Commentary parallel text.


title: A rational analysis of curiosity

abstract: We present a rational analysis of curiosity, proposing that people's curiosity is driven by seeking stimuli that maximize their ability to make appropriate responses in the future. This perspective offers a way to unify previous theories of curiosity into a single framework. Experimental results confirm our model's predictions, showing how the relationship between curiosity and confidence can change significantly depending on the nature of the environment. Please refer to https://psyarxiv.com/wg5m6/ for a more updated version of this manuscript with a more detailed modeling section with extensive experiments.


title: Reconfiguring the Imaging Pipeline for Computer Vision

abstract: Advancements in deep learning have ignited an explosion of research on efficient hardware for embedded computer vision. Hardware vision acceleration, however, does not address the cost of capturing and processing the image data that feeds these algorithms. We examine the role of the image signal processing (ISP) pipeline in computer vision to identify opportunities to reduce computation and save energy. The key insight is that imaging pipelines should be designed to be configurable: to switch between a traditional photography mode and a low-power vision mode that produces lower-quality image data suitable only for computer vision. We use eight computer vision algorithms and a reversible pipeline simulation tool to study the imaging system's impact on vision performance. For both CNN-based and classical vision algorithms, we observe that only two ISP stages, demosaicing and gamma compression, are critical for task performance. We propose a new image sensor design that can compensate for skipping these stages. The sensor design features an adjustable resolution and tunable analog-to-digital converters (ADCs). Our proposed imaging system's vision mode disables the ISP entirely and configures the sensor to produce subsampled, lower-precision image data. This vision mode can save ~75% of the average energy of a baseline photography mode while having only a small impact on vision task accuracy.


title: On the records

abstract: World record setting has long attracted public interest and scientific investigation. Extremal records summarize the limits of the space explored by a process, and the historical progression of a record sheds light on the underlying dynamics of the process. Existing analyses of prediction, statistical properties, and ultimate limits of record progressions have focused on particular domains. However, a broad perspective on how record progressions vary across different spheres of activity needs further development. Here we employ cross-cutting metrics to compare records across a variety of domains, including sports, games, biological evolution, and technological development. We find that these domains exhibit characteristic statistical signatures in terms of rates of improvement, "burstiness" of record-breaking time series, and the acceleration of the record breaking process. Specifically, sports and games exhibit the slowest rate of improvement and a wide range of rates of "burstiness." Technology improves at a much faster rate and, unlike other domains, tends to show acceleration in records. Many biological and technological processes are characterized by constant rates of improvement, showing less burstiness than sports and games. It is important to understand how these statistical properties of record progression emerge from the underlying dynamics. Towards this end, we conduct a detailed analysis of a particular record-setting event: elite marathon running. In this domain, we find that studying record-setting data alone can obscure many of the structural properties of the underlying process. The marathon study also illustrates how some of the standard statistical assumptions underlying record progression models may be inappropriate or commonly violated in real-world datasets.


title: HAT-P-26b: A Neptune-Mass Exoplanet with a Well Constrained Heavy Element Abundance

abstract: A correlation between giant-planet mass and atmospheric heavy elemental abundance was first noted in the past century from observations of planets in our own Solar System, and has served as a cornerstone of planet formation theory. Using data from the Hubble and Spitzer Space Telescopes from 0.5 to 5 microns, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent H2O absorption bands with a maximum base-to-peak amplitude of 525ppm in the transmission spectrum. Using the water abundance as a proxy for metallicity, we measured HAT-P-26b's atmospheric heavy element content [4.8 (-4.0 +21.5) times solar]. This likely indicates that HAT-P-26b's atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals.


title: Lieb-Robinson bounds on $n$-partite connected correlation functions

abstract: Lieb and Robinson provided bounds on how fast bipartite connected correlations can arise in systems with only short-range interactions. We generalize Lieb-Robinson bounds on bipartite connected correlators to multipartite connected correlators. The bounds imply that an $n$-partite connected correlator can reach unit value in constant time. Remarkably, the bounds also allow for an $n$-partite connected correlator to reach a value that is exponentially large with system size in constant time, a feature which stands in contrast to bipartite connected correlations. We provide explicit examples of such systems.


title: Comparing probabilistic predictive models applied to football

abstract: We propose two Bayesian multinomial-Dirichlet models to predict the final outcome of football (soccer) matches and compare them to three well-known models regarding their predictive power. All the models predicted the full-time results of 1710 matches of the first division of the Brazilian football championship and the comparison used three proper scoring rules, the proportion of errors and a calibration assessment. We also provide a goodness of fit measure. Our results show that multinomial-Dirichlet models are not only competitive with standard approaches, but they are also well calibrated and present reasonable goodness of fit.


title: Fitting phase--type scale mixtures to heavy--tailed data and distributions

abstract: We consider the fitting of heavy tailed data and distribution with a special attention to distributions with a non--standard shape in the "body" of the distribution. To this end we consider a dense class of heavy tailed distributions introduced recently, employing an EM algorithm for the the maximum likelihood estimates of its parameters. We present methods for fitting to observed data, histograms, censored data, as well as to theoretical distributions. Numerical examples are provided with simulated data and a benchmark reinsurance dataset. We empirically demonstrate that our model can provide excellent fits to heavy--tailed data/distributions with minimal assumptions


title: Object-Level Context Modeling For Scene Classification with Context-CNN

abstract: Convolutional Neural Networks (CNNs) have been used extensively for computer vision tasks and produce rich feature representation for objects or parts of an image. But reasoning about scenes requires integration between the low-level feature representations and the high-level semantic information. We propose a deep network architecture which models the semantic context of scenes by capturing object-level information. We use Long Short Term Memory(LSTM) units in conjunction with object proposals to incorporate object-object relationship and object-scene relationship in an end-to-end trainable manner. We evaluate our model on the LSUN dataset and achieve results comparable to the state-of-art. We further show visualization of the learned features and analyze the model with experiments to verify our model's ability to model context.


title: An Extreme Rotating Black Hole in New Massive Gravity Theory

abstract: New Massive Gravity is an alternative theory to General Relativity that is used to describe the gravitational field in a (2+1)-dimensional spacetime. Black hole solutions have been found in this theory, in particular an asymptotically anti-de Sitter rotating black hole. We analyse some features of this solution as its event horizon, black hole area and distance to the horizon, specially in the rotating extreme case, showing that they have shared features with extreme black holes in 4-dimensional General relativity. This limit case is interesting in the search of geometric inequalities as the ones found for the Kerr black hole in (3+1)-General Relativity.


title: Group and round quadratic forms

abstract: We offer some elementary characterisations of group and round quadratic forms. These characterisations are applied to establish new (and recover existing) characterisations of Pfister forms. We establish "going-up" results for group and anisotropic round forms with respect to iterated Laurent series fields, which contrast with the established results with respect to rational function field extensions. For forms of two-power dimension, we determine when there exists a field extension over which the form becomes an anisotropic group form that is not round.


title: The MUSCLES Treasury Survey IV: Scaling Relations for Ultraviolet, Ca II K, and Energetic Particle Fluxes from M Dwarfs

abstract: Characterizing the UV spectral energy distribution (SED) of an exoplanet host star is critically important for assessing its planet's potential habitability, particularly for M dwarfs as they are prime targets for current and near-term exoplanet characterization efforts and atmospheric models predict that their UV radiation can produce photochemistry on habitable zone planets different than on Earth. To derive ground-based proxies for UV emission for use when Hubble Space Telescope observations are unavailable, we have assembled a sample of fifteen early-to-mid M dwarfs observed by Hubble, and compared their non-simultaneous UV and optical spectra. We find that the equivalent width of the chromospheric Ca II K line at 3933 Angstroms, when corrected for spectral type, can be used to estimate the stellar surface flux in ultraviolet emission lines, including H I Lyman alpha. In addition, we address another potential driver of habitability: energetic particle fluxes associated with flares. We present a new technique for estimating soft X-ray and >10 MeV proton flux during far-UV emission line flares (Si IV and He II) by assuming solar-like energy partitions. We analyze several flares from the M4 dwarf GJ 876 observed with Hubble and Chandra as part of the MUSCLES Treasury Survey and find that habitable zone planets orbiting GJ 876 are impacted by large Carrington-like flares with peak soft X-ray fluxes >1e-3 W m-2 and possible proton fluxes ~100-1000 pfu, approximately four orders of magnitude more frequently than modern-day Earth.


title: The minimal ideal in multiplier algebras

abstract: Let $\mathcal A$ be a simple, $\sigma$-unital, non-unital, non-elementary C*-algebra and let $I_{min}$ be the intersection of all the ideals of $\mathcal M(\mathcal A)$ that properly contain $\mathcal A$. $I_{min}$ coincides with the ideal defined by Lin (Simple C*-algebras with continuous scales and simple corona algebras. 112, (1991) Proc. Amer.Math. Soc) in terms of approximate units of $\mathcal A$ and $I_{min}/\mathcal A$ is purely infinite and simple. If $\mathcal A$ is separable, or if $\mathcal A$ has the (SP) property and its dimension semigroup $D(\mathcal A)$ of Murray-von Neumann equivalence classes of projections of $\mathcal A$ is order separable, or if $\mathcal A$ has strict comparison of positive elements by traces, then $\mathcal A\ne I_{min}$. If the tracial simplex $ \mathcal T(\mathcal A)$ is nonempty, let $ I_{con}$ be the closure of the linear span of the elements $A\in\mathcal M(\mathcal A)_+$ such that the evaluation map $\hat A(\tau)=\tau(A)$ is continuous. If $\mathcal A$ has strict comparison of positive element by traces then $I_{min}= I_{con}$. Furthermore, $I_{min}$ too has strict comparison of positive elements in the sense that if $A, B\in (I_{min})_+$, $B\not \in\mathcal A$ and $d_\tau(A)< d_\tau(B)$ for all $\tau\in \mathcal T(\mathcal A)$ for which $d_\tau(B)< \infty$, then $A\preceq B$. However if $\mathcal A$ does not have strict comparison of positive elements by traces then $I_{min}\ne I_{con}$ can occur: a counterexample is provided by Villadsen's AH algebras without slow dimension growth. If the dimension growth is flat, $ I_{con}$ is the largest proper ideal of $\mathcal M(\mathcal A)$.


title: Quantum graphs with the Bethe-Sommerfeld property

abstract: In contrast to the usual quantum systems which have at most a finite number of open spectral gaps if they are periodic in more than one direction, periodic quantum graphs may have gaps arbitrarily high in the spectrum. This property of graph Hamiltonians, being generic in a sense, inspires the question about the existence of graphs with a finite and nonzero number of spectral gaps. We show that the answer depends on the vertex couplings together with commensurability of the graph edges. A finite and nonzero number of gaps is excluded for graphs with scale invariant couplings; on the other hand, we demonstrate that graphs featuring a finite nonzero number of gaps do exist, illustrating the claim on the example of a rectangular lattice with a suitably tuned $\delta$-coupling at the vertices.


title: Macroscopic Quantum Tunneling Escape of Bose-Einstein Condensates

abstract: Recent experiments on macroscopic quantum tunneling reveal a non-exponential decay of the number of atoms trapped in a quasibound state behind a potential barrier. Through both experiment and theory, we demonstrate this non-exponential decay results from interactions between atoms. Quantum tunneling of tens of thousands of 87 Rb atoms in a Bose-Einstein condensate is modeled by a modified Jeffreys-Wentzel-Kramers-Brillouin model, taking into account the effective time-dependent barrier induced by the mean-field. Three-dimensional Gross-Pitaevskii simulations corroborate a mean-field result when compared with experiments. However, with one-dimensional modeling using time-evolving block decimation, we present an effective renormalized mean-field theory that suggests many-body dynamics for which a bare mean-field theory may not apply.


title: Multiparticle production in the large lambda n limit: Realising Higgsplosion in a scalar QFT

abstract: In a scalar theory which we use as a simplified model for the Higgs sector, we adopt the semiclassical formalism of Son for computations of $n$-particle production cross-sections in the high-multiplicity $n\to \infty$ weak-coupling $\lambda \to 0$ regime with the value of $\lambda n$ held fixed and large. The approach relies on the use of singular classical solutions to a certain boundary value problem. In the past this formalism has been successfully used and verified in computations of perturbative multi-particle processes at tree-level, and also at the next-to-leading order level in the small $\lambda n$ expansion near the multi-particle mass threshold. We apply this singular solutions formalism in the regime of ultra-high multiplicities where $\lambda n \gg 1$, and compute the leading positive $\sim n\,\sqrt{\lambda n}$ contribution to the exponent of the multi-particle rate in this large $\lambda n$ limit. The computation is carried out near the multi-particle mass threshold where the multiplicity $n$ approaches its maximal value allowed by kinematics. This calculation relies on the idea of Gorsky and Voloshin to use a thin wall approximation for the singular solutions that resemble critical bubbles. This approximation is justified in precisely the high-multiplicity $\sqrt{\lambda n} \to \infty$ regime of interest. Based on our results we show that the scalar theory with a spontaneous symmetry breaking used here as a simplified model for the Higgs sector, is very likely to realise the high-energy Higgsplosion phenomenon.


title: Assessment of Bayesian Expected Power via Bayesian Bootstrap

abstract: The Bayesian expected power (BEP) has become increasingly popular in sample size determination and assessment of the probability of success (POS) for a future trial. The BEP takes into consideration the uncertainty around the parameters assumed by a power analysis and is thus more robust compared to the traditional power that assumes a single set of parameters. Current methods for assessing BEP are often based in a parametric framework by imposing a model on the pilot data to derive and sample from the posterior distributions of the parameters. Implementation of the model-based approaches can be analytically challenging and computationally costly especially for multivariate data sets; it also runs the risk of generating misleading BEP if the model is mis-specified. We propose an approach based on the Bayesian bootstrap technique (BBS) to simulate future trials in the presence of individual-level pilot data, based on which the empirical BEP can be calculated. The BBS approach is model-free with no assumptions about the distribution of the prior data and circumvents the analytical and computational complexity associated with obtaining the posterior distribution of the parameters. Information from multiple pilot studies is also straightforward to combine. We also propose the double bootstrap (BS2), a frequentist counterpart to the BBS, that shares similar properties and achieves the same goal as the BBS for BEP assessment. Simulation studies and case studies are presented to demonstrate the implementation of the BBS and BS2 techniques and to compare the BEP results with model-based approaches.


title: Recognising Axionic Dark Matter by Compton and de-Broglie Scale Modulation of Pulsar Timing

abstract: Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background.


title: Lattice Gas with Molecular Dynamics Collision Operator

abstract: We introduce a lattice gas implementation that is based on coarse-graining a Molecular Dynamics (MD) simulation. Such a lattice gas is similar to standard lattice gases, but its collision operator is informed by an underlying MD simulation. This can be considered an optimal lattice gas implementation because it allows for the representation of any system that can be simulated with MD. We show here that equilibrium behavior of the popular lattice Boltzmann algorithm is consistent with this optimal lattice gas. This comparison allows us to make a more accurate identification of the expressions for temperature and pressure in lattice Boltzmann simulations which turn out to be related not only to the physical temperature and pressure but also to the lattice discretization. We show that for any spatial discretization we need to choose a particular temporal discretization to recover the lattice Boltzmann equilibrium.


title: Analysis of the finite element method for the Laplace--Beltrami equation on surfaces with regions of high curvature using graded meshes

abstract: We derive error estimates for the piecewise linear finite element approximation of the Laplace--Beltrami operator on a bounded, orientable, $C^3$, surface without boundary on general shape regular meshes. As an application, we consider a problem where the domain is split into two regions: one which has relatively high curvature and one that has low curvature. Using a graded mesh we prove error estimates that do not depend on the curvature on the high curvature region. Numerical experiments are provided.


title: On the local well-posedness of Lovelock and Horndeski theories

abstract: We investigate local well-posedness of the initial value problem for Lovelock and Horndeski theories of gravity. A necessary condition for local well-posedness is strong hyperbolicity of the equations of motion. Even weak hyperbolicity can fail for strong fields so we restrict to weak fields. The Einstein equation is known to be strongly hyperbolic in harmonic gauge so we study Lovelock theories in harmonic gauge. We show that the equation of motion is always weakly hyperbolic for weak fields but, in a generic weak-field background, it is not strongly hyperbolic. For Horndeski theories, we prove that, for weak fields, the equation of motion is always weakly hyperbolic in any generalized harmonic gauge. For some Horndeski theories there exists a generalized harmonic gauge for which the equation of motion is strongly hyperbolic in a weak-field background. This includes "k-essence" like theories. However, for more general Horndeski theories, there is no generalized harmonic gauge for which the equation of motion is strongly hyperbolic in a generic weak-field background. Our results show that the standard method used to establish local well-posedness of the Einstein equation does not extend to Lovelock or general Horndeski theories. This raises the possibility that these theories may not admit a well-posed initial value problem even for weak fields.


title: Trajectory Planning of Automated Vehicles in Tube-like Road Segments

abstract: This paper presents a method based on linear programming for trajectory planning of automated vehicles, combining obstacle avoidance, time scheduling for the reaching of waypoints and time-optimal traversal of tube-like road segments. System modeling is conducted entirely spatial-based. Kinematic vehicle dynamics as well as time are expressed in a road-aligned coordinate frame with path along the road centerline serving as the dependent variable. We elaborate on control rate constraints in the spatial domain. A vehicle dimension constraint heuristic is proposed to constrain vehicle dimensions inside road boundaries. It is outlined how friction constraints are accounted for. The discussion is extended to dynamic vehicle models. The benefits of the proposed method are illustrated by a comparison to a time-based method.


title: On covering systems of integers

abstract: A covering system of the integers is a finite collection of modular residue classes $\{a_m \bmod{m}\}_{m \in S}$ whose union is all integers. Given a finite set $S$ of moduli, it is often difficult to tell whether there is a choice of residues modulo elements of $S$ covering the integers. Hough has shown that if the smallest modulus in $S$ is at least $10^{16}$, then there is none. However, the question of whether there is a covering of the integers with all odd moduli remains open. We consider multiplicative restrictions on the set of moduli to generalize Hough's negative solution to the minimum modulus problem. In particular, we find that every covering system of the integers has a modulus divisible by a prime number less than or equal to $19$. Hough and Nielsen have shown that every covering system has a modulus divisible by either $2$ or $3$.


title: Theory of Dark Matter

abstract: We discuss the hypothesis that the constituents of dark matter in the galactic halo are Primordial Intermediate-Mass Black Holes (PIMBHs). The status of axions and WIMPs is discussed, as are the methods for detecting PIMBHs with emphasis on microlensing. The role of the angular momentum J of the PIMBHs in their escaping previous detection is considered.


title: Optimal fidelity multi-level Monte Carlo for quantification of uncertainty in simulations of cloud cavitation collapse

abstract: We quantify uncertainties in the location and magnitude of extreme pressure spots revealed from large scale multi-phase flow simulations of cloud cavitation collapse. We examine clouds containing 500 cavities and quantify uncertainties related to their initial spatial arrangement. The resulting 2000-dimensional space is sampled using a non-intrusive and computationally efficient Multi-Level Monte Carlo (MLMC) methodology. We introduce novel optimal control variate coefficients to enhance the variance reduction in MLMC. The proposed optimal fidelity MLMC leads to more than two orders of magnitude speedup when compared to standard Monte Carlo methods. We identify large uncertainties in the location and magnitude of the peak pressure pulse and present its statistical correlations and joint probability density functions with the geometrical characteristics of the cloud. Characteristic properties of spatial cloud structure are identified as potential causes of significant uncertainties in exerted collapse pressures.


title: Stable small bubble clusters in two-dimensional foams

abstract: Key features of the mechanical response of amorphous particulate materials, such as foams, emulsions, and granular media, to applied stress are determined by the frequency and size of particle rearrangements that occur as the system transitions from one mechanically stable state to another. This work describes coordinated experimental and computational studies of bubble rafts, which are quasi-two dimensional systems of bubbles confined to the air-water interface. We focus on small mechanically stable clusters of four, five, six, and seven bubbles with two different sizes with diameter ratio $\sigma_L/\sigma_S = 1.4$. Focusing on small bubble clusters, which can be viewed as subsystems of a larger system, allows us to investigate the full ensemble of clusters that form, measure the respective frequencies with which the clusters occur, and determine the form of the bubble-bubble interactions. We emphasize several important results. First, for clusters with N > 5 bubbles, we find using discrete element simulations that short-range attractive interactions between bubbles give rise to a larger ensemble of distinct mechanically stable clusters compared to that generated by long-range attractive interactions. The additional clusters in systems with short-range attractions possess larger gaps between pairs of neighboring bubbles on the periphery of the clusters. The ensemble of bubble clusters observed in experiments is similar to the ensemble of clusters with long-range attractive interactions. We also compare the frequency with which each cluster occurs in simulations and experiments. We find that the cluster frequencies are extremely sensitive to the protocol used to generate them and only weakly correlated to the energy of the clusters.


title: Cellular Automaton-Like Model of Arising Physical-Like Properties

abstract: A simple relation of the order of $n$ abstract objects generates an $n-2$ dimensional basis of three dimensional vectors. A cellular automaton-like model of evolution of this system is postulated. During this evolution, some quantities stabilize with time and form a discrete spectrum of values. The presented model may have some general aspects in common with a cellular automaton representation of a quantum system.


title: Linear and non-linear flow mode in Pb-Pb collisions at $\sqrt{s_{\rm NN}} =$ 2.76 TeV

abstract: The second and the third order anisotropic flow, $V_{2}$ and $V_3$, are mostly determined by the corresponding initial spatial anisotropy coefficients, $\varepsilon_{2}$ and $\varepsilon_{3}$, in the initial density distribution. In addition to their dependence on the same order initial anisotropy coefficient, higher order anisotropic flow, $V_n$ ($n > 3$), can also have a significant contribution from lower order initial anisotropy coefficients, which leads to mode-coupling effects. In this Letter we investigate the linear and non-linear modes in higher order anisotropic flow $V_n$ for $n=4$, $5$, $6$ with the ALICE detector at the Large Hadron Collider. The measurements are done for particles in the pseudorapidity range $|\eta| < 0.8$ and the transverse momentum range $0.2 < p_{\rm T} < 5.0$ GeV/$c$ as a function of collision centrality. The results are compared with theoretical calculations and provide important constraints on the initial conditions, including initial spatial geometry and its fluctuations, as well as the ratio of the shear viscosity to entropy density of the produced system.


title: An overview and comparative analysis of Recurrent Neural Networks for Short Term Load Forecasting

abstract: The key component in forecasting demand and consumption of resources in a supply network is an accurate prediction of real-valued time series. Indeed, both service interruptions and resource waste can be reduced with the implementation of an effective forecasting system. Significant research has thus been devoted to the design and development of methodologies for short term load forecasting over the past decades. A class of mathematical models, called Recurrent Neural Networks, are nowadays gaining renewed interest among researchers and they are replacing many practical implementation of the forecasting systems, previously based on static methods. Despite the undeniable expressive power of these architectures, their recurrent nature complicates their understanding and poses challenges in the training procedures. Recently, new important families of recurrent architectures have emerged and their applicability in the context of load forecasting has not been investigated completely yet. In this paper we perform a comparative study on the problem of Short-Term Load Forecast, by using different classes of state-of-the-art Recurrent Neural Networks. We test the reviewed models first on controlled synthetic tasks and then on different real datasets, covering important practical cases of study. We provide a general overview of the most important architectures and we define guidelines for configuring the recurrent networks to predict real-valued time series.


title: The Network Nullspace Property for Compressed Sensing of Big Data over Networks

abstract: We present a novel condition, which we term the net- work nullspace property, which ensures accurate recovery of graph signals representing massive network-structured datasets from few signal values. The network nullspace property couples the cluster structure of the underlying network-structure with the geometry of the sampling set. Our results can be used to design efficient sampling strategies based on the network topology.


title: ROCKER: A Refinement Operator for Key Discovery

abstract: The Linked Data principles provide a decentral approach for publishing structured data in the RDF format on the Web. In contrast to structured data published in relational databases where a key is often provided explicitly, finding a set of properties that allows identifying a resource uniquely is a non-trivial task. Still, finding keys is of central importance for manifold applications such as resource deduplication, link discovery, logical data compression and data integration. In this paper, we address this research gap by specifying a refinement operator, dubbed ROCKER, which we prove to be finite, proper and non-redundant. We combine the theoretical characteristics of this operator with two monotonicities of keys to obtain a time-efficient approach for detecting keys, i.e., sets of properties that describe resources uniquely. We then utilize a hash index to compute the discriminability score efficiently. Therewith, we ensure that our approach can scale to very large knowledge bases. Results show that ROCKER yields more accurate results, has a comparable runtime, and consumes less memory w.r.t. existing state-of-the-art techniques.