title: Exploring circular polarization in the CMB due to conventional sources of cosmic birefringence

abstract: The circular polarization of the cosmic microwave background (CMB) is usually taken to be zero since it is not generated by Thomson scattering. Here we explore the actual level of circular polarization in the CMB generated by conventional cosmological sources of birefringence. We consider two classes of mechanisms for birefringence. One is alignment of the matter to produce an anisotropic susceptibility tensor: the hydrogen spins can be aligned either by density perturbations or CMB anisotropies themselves. The other is anisotropy of the radiation field coupled to the non-linear response of the medium to electromagnetic fields: this can occur either via photon-photon scattering (non-linear response of the vacuum); atomic hyperpolarizability (non-linear response of neutral atoms); or plasma delay (non-linear response of free electrons). The strongest effect comes from photon-photon scattering from recombination at a level of $\sim 10^{-14}$ K. Our results are consistent with a negligible circular polarization of the CMB in comparison with the linear polarization or the sensitivity of current and near-term experiments.


title: Graph Algorithms for Topology Identification using Power Grid Probing

abstract: To perform any meaningful optimization task, power distribution operators need to know the topology and line impedances of their electric networks. Nevertheless, distribution grids currently lack a comprehensive metering infrastructure. Although smart inverters are widely used for control purposes, they have been recently advocated as the means for an active data acquisition paradigm: Reading the voltage deviations induced by intentionally perturbing inverter injections, the system operator can potentially recover the electric grid topology. Adopting inverter probing for feeder processing, a suite of graph-based topology identification algorithms is developed here. If the grid is probed at all leaf nodes but voltage data are metered at all nodes, the entire feeder topology can be successfully recovered. When voltage data are collected only at probing buses, the operator can find a reduced feeder featuring key properties and similarities to the actual feeder. To handle modeling inaccuracies and load non-stationarity, noisy probing data need to be preprocessed. If the suggested guidelines on the magnitude and duration of probing are followed, the recoverability guarantees carry over from the noiseless to the noisy setup with high probability.


title: Thermoelectric current in topological insulator nanowires with impurities

abstract: In this paper we consider charge current generated by maintaining a temperature difference over a nanowire at zero voltage bias. For topological insulator nanowires in a perpendicular magnetic field the current can change sign as the temperature of one end is increased. Here we study how this thermoelectric current sign reversal depends on magnetic field and how impurities affect the size of the thermoelectric current. We consider both scalar and magnetic impurities and show that their influence on the current are quite similar, although the magnetic impurities seem to be more effective in reducing the effect. For moderate impurity concentration the sign reversal persists.


title: Construction of a class of scalar quantum field theory models in any dimensions

abstract: It is observed that certain convex envelopes of Wightman type functionals corresponding to scalar, stochastically positive quantum fields consist of Wightman type functionals only .This leads to the construction of a large classes of not quasi-free scalar quantum field theory models obeying all Wightman axioms in any dimensions


title: Sorting by Swaps with Noisy Comparisons

abstract: We study sorting of permutations by random swaps if each comparison gives the wrong result with some fixed probability $p<1/2$. We use this process as prototype for the behaviour of randomized, comparison-based optimization heuristics in the presence of noisy comparisons. As quality measure, we compute the expected fitness of the stationary distribution. To measure the runtime, we compute the minimal number of steps after which the average fitness approximates the expected fitness of the stationary distribution. We study the process where in each round a random pair of elements at distance at most $r$ are compared. We give theoretical results for the extreme cases $r=1$ and $r=n$, and experimental results for the intermediate cases. We find a trade-off between faster convergence (for large $r$) and better quality of the solution after convergence (for small $r$).


title: Discontinuity propagation in delay differential-algebraic equations

abstract: The propagation of primary discontinuities in initial value problems for linear delay differential-algebraic equations (DDAEs) is discussed. Based on the (quasi-) Weierstra{\ss} form for regular matrix pencil, a complete characterization of the different propagation types is given and algebraic criteria in terms of the matrices are developed. The analysis, which is based on the method of steps, takes into account all possible inhomogeneities and history functions and thus serves as the worst-case scenario. Moreover, it reveals possible hidden delays in the DDAE. The new classification for DDAEs is compared to existing approaches in the literature and the impact of splicing conditions on the classification is studied.


title: On the continuity of entropy of Lorenz maps

abstract: We consider a one parameter family of Lorenz maps indexed by their point of discontinuity $p$ and constructed from a pair of bilipschitz functions. We prove that their topological entropies vary continuously as a function of $p$ and discuss Milnor's monotonicity conjecture in this setting.


title: Predicting polarisation hemisphere switch in C60 caused by the motion of an internal point charge with an electrostatic and quantum chemistry solutions

abstract: A classical electrostatic solution for polarisation charge on the interface of a dielectric cavity interacting with an internal point charge is presented. This solution comes from the Gauss electrostatic potential as a sum of two terms, the cavity and the point charge, expanded with Legendre polynomials. Subsequent application of the Dietrich-Newman boundary conditions defines the problem of emptiness in the inside and isotropic dielectric medium on the outside of the spherical interface to obtain an equation which describes the surface charge density on the interface. These results are compared with quantum chemical calculations using density functional theory for neutral C60 fullerenes. Comparison showed that there was good qualitative agreement between the classical electrostatic theory and quantum calculations. The polarisation effect that occurs as a result of the motion of the trapped particle inside the C60 molecule shows potential for a polarisable nanoswitch which might be used in nanotechnology as electronic component.


title: Effects of Topology Knowledge and Relay Depth on Asynchronous Consensus

abstract: Consider a point-to-point message-passing network. We are interested in the asynchronous crash-tolerant consensus problem in incomplete networks. We study the feasibility and efficiency of approximate consensus under different restrictions on topology knowledge and the relay depth, i.e., the maximum number of hops any message can be relayed. These two constraints are common in large-scale networks, and are used to avoid memory overload and network congestion respectively. Specifically, for different values of integers k, k , we consider that each node knows all its neighbors of at most k-hop distance (k-hop topology knowledge), and the relay depth is k . We consider both directed and undirected graphs. More concretely, we answer the following main question in asynchronous systems: What is a tight condition on the underlying communication graphs for achieving approximate consensus if each node has only a k-hop topology knowledge and relay depth k? To prove that the necessary conditions presented in the paper are also sufficient, we have developed algorithms that achieve consensus in graphs satisfying those conditions: -The first class of algorithms requires k-hop topology knowledge and relay depth k. Unlike prior algorithms, these algorithms do not flood the network, and each node does not need the full topology knowledge. We show how the convergence time and the message complexity of those algorithms is affected by k, providing the respective upper bounds. -The second set of algorithms requires only one-hop neighborhood knowledge, i.e., immediate incoming and outgoing neighbors, but needs to flood the network (i.e., relay depth is n, where n is the number of nodes). One result that may be of independent interest is a topology discovery mechanism to learn and "estimate" the topology in asynchronous directed networks with crash faults.


title: Similar but Different: Exploiting Users' Congruity for Recommendation Systems

abstract: The pervasive use of social media provides massive data about individuals' online social activities and their social relations. The building block of most existing recommendation systems is the similarity between users with social relations, i.e., friends. While friendship ensures some homophily, the similarity of a user with her friends can vary as the number of friends increases. Research from sociology suggests that friends are more similar than strangers, but friends can have different interests. Exogenous information such as comments and ratings may help discern different degrees of agreement (i.e., congruity) among similar users. In this paper, we investigate if users' congruity can be incorporated into recommendation systems to improve it's performance. Experimental results demonstrate the effectiveness of embedding congruity related information into recommendation systems.


title: Unveiling the physical conditions of the youngest disks: A warm embedded disk in L1527

abstract: [Abridged] Protoplanetary disks have been studied extensively, both physically and chemically, to understand the environment in which planets form. However, the first steps of planet formation are likely to occur already when the protostar and disk are still embedded in their natal envelope. The initial conditions for planet formation may thus be provided by these young embedded disks, of which the physical and chemical structure is poorly characterized. We aim to constrain the midplane temperature structure, one of the critical unknowns, of the embedded disk around L1527. In particular, we set out to determine whether there is an extended cold outer region where CO is frozen out, as is the case for Class II disks. We use archival ALMA data to directly observe the midplane of the near edge-on L1527 disk. Optically thick $^{13}$CO ($J=2-1$) and C$^{18}$O ($J=2-1$) emission is observed throughout the disk and inner envelope, while N$_2$D$^+ (J=3-2$), which can only be abundant when CO is frozen out, is not detected. Both CO isotopologues have brightness temperatures $\gtrsim$ 25 K along the midplane. Disk and envelope emission can be disentangled kinematically, because the largest velocities are reached in the disk. A power law radial temperature profile constructed using the highest midplane temperature at these velocities suggest that the temperature is above 20 K out to at least 75 AU, and possibly throughout the entire 125 AU disk. Radiative transfer models show that a model without CO freeze-out in the disk matches the C$^{18}$O observations better than a model with the CO snowline at $\sim$70 AU. In addition, there is no evidence for a large (order of magnitude) depletion of CO. The disk around L1527 is likely to be warm enough to have CO present in the gas phase throughout the disk, suggesting that young embedded disks can indeed be warmer than the more evolved Class II disks.


title: Explicit inverse of tridiagonal matrix with applications in autoregressive modeling

abstract: We present the explicit inverse of a class of symmetric tridiagonal matrices which is almost Toeplitz, except that the first and last diagonal elements are different from the rest. This class of tridiagonal matrices are of special interest in complex statistical models which uses the first order autoregression to induce dependence in the covariance structure, for instance, in econometrics or spatial modeling. They also arise in interpolation problems using the cubic spline. We show that the inverse can be expressed as a linear combination of Chebyshev polynomials of the second kind and present results on the properties of the inverse, such as bounds on the row sums, the trace of the inverse and its square, and their limits as the order of the matrix increases.


title: T-duality and high-derivative gravity theories: the BTZ black hole/string paradigm

abstract: We show that the temperature and entropy of a BTZ black hole are invariant under T-duality to next to leading order in $M_\star^{-2}$, $M_\star$ being the scale suppressing higher-curvature/derivative terms in the Lagrangian. We work in the framework of a two-parameter family of theories exhibiting T-duality, which includes (but goes beyond) String Theory. Interestingly enough, the AdS/CFT correspondence enforces quantization conditions on these parameters. In the particular case of bosonic/heterotic string theory, our results extend those of a classical paper by Horowitz and Welch. For generic (albeit quantized) values of the parameters, it suggests that T-duality might be an interesting tool to constrain consistent low-energy effective actions while entailing physical equivalences outside String Theory. Moreover, it generates a new family of regular asymptotically flat black string solutions in three-dimensions.


title: On multivariate modifications of Cramer Lundberg risk model with constant intensities

abstract: The paper considers very general multivariate modifications of Cramer-Lundberg risk model. The claims can be of different types and can arrive in groups. The groups arrival processes within a type have constant intensities. The counting groups processes are dependent multivariate compound Poisson processes of type I. We allow empty groups and show that in that case we can find stochastically equivalent Cramer-Lundberg model with non-empty groups. The investigated model generalizes the risk model with common shocks, the Poisson risk process of order k, the Poisson negative binomial, the Polya-Aeppli, the Polya-Aeppli of order k among others. All of them with one or more types of polices. The relations between the numerical characteristics and distributions of the components of the risk processes are proven to be corollaries of the corresponding formulae of the Cramer-Lundberg risk model.


title: Analytical results for the quantum non-Markovianity of spin ensembles undergoing pure dephasing dynamics

abstract: We study analytically the non-Markovianity of a spin ensemble, with arbitrary number of spins and spin quantum number, undergoing a pure dephasing dynamics. The system is considered as a part of a larger spin ensemble of any geometry with pairwise interactions. We derive exact formulas for the reduced dynamics of the system and for its non-Markovianity as assessed by the witness of Lorenzo \emph{et al.} [Phys.~Rev.~A \textbf{88}, 020102(R) (2013)]. The non-Markovianity is further investigated in the thermodynamic limit when the environment's size goes to infinity. In this limit and for finite-size systems, we find that the Markovian's character of the system's dynamics crucially depends on the range of the interactions. We also show that, when the system and its environment are initially in a separable state, the appearance of non-Markovianity is independent of the entanglement generation between the system and its environment.


title: Is the security of quantum cryptography guaranteed by the laws of physics?

abstract: It is often claimed that the security of quantum key distribution (QKD) is guaranteed by the laws of physics. However, this claim is content-free if the underlying theoretical definition of QKD is not actually compatible with the laws of physics. This paper observes that (1) the laws of physics pose serious obstacles to the security of QKD and (2) the same laws are ignored in all QKD "security proofs".


title: A note on some sub-Gaussian random variables

abstract: In [8] the author of this paper continued the research on the complex-valued discrete random variables $X_l(m,N)$ ($0\le l\le N-1$, $1\le M\le N)$ recently introduced and studied in [24]. Here we extend our results by considering $X_l(m,N)$ as sub-Gaussian random variables. Our investigation is motivated by the known fact thatthe so-called Restricted Isometry Property (RIP) introduced in [4] holds with high probability for any matrix generated by a sub-Gaussian random variable. Notice that sensing matrices with the RIP play a crucial role in Theory of compressive sensing. Our main results concern the proofs of the lower and upper bound estimates of the expected values of the random variables $|X_l(m,N)|$, $|U_l(m,N)|$ and $|V_l(m,N)|$, where $U_l(m,N)$ and $U_l(m,N)$ are the real and the imaginary part of $X_l(m,N)$, respectively. These estimates are also given in terms of related sub-Gaussian norm $\Vert \cdot\Vert_{\psi_2}$ considered in [28]. Moreover, we prove a refinement of the mentioned upper bound estimates for the real and the imaginary part of $X_l(m,N)$.


title: Some limit laws for quantum walks with applications to a version of the Parrondo paradox

abstract: A quantum walker moves on the integers with four extra degrees of freedom, performing a coin-shift operation to alter its internal state and position at discrete units of time. The time evolution is described by a unitary process. We focus on finding the limit probability law for the position of the walker and study it by means of Fourier analysis. The quantum walker exhibits both localization and a ballistic behavior. Our two results are given as limit theorems for a 2-period time-dependent walk and they describe the location of the walker after it has repeated the unitary process a large number of times. The theorems give an analytical tool to study some of the Parrondo type behavior in a quantum game which was studied by J. Rajendran and C. Benjamin by means of very nice numerical simulations [1]. With our analytical tools at hand we can easily explore the "phase space" of parameters of one of the games, similar to the winning game in their papers. We include numerical evidence that our two games, similar to theirs, exhibit a Parrondo type paradox.


title: Event-based Moving Object Detection and Tracking

abstract: Event-based vision sensors, such as the Dynamic Vision Sensor (DVS), are ideally suited for real-time motion analysis. The unique properties encompassed in the readings of such sensors provide high temporal resolution, superior sensitivity to light and low latency. These properties provide the grounds to estimate motion extremely reliably in the most sophisticated scenarios but they come at a price - modern event-based vision sensors have extremely low resolution and produce a lot of noise. Moreover, the asynchronous nature of the event stream calls for novel algorithms. This paper presents a new, efficient approach to object tracking with asynchronous cameras. We present a novel event stream representation which enables us to utilize information about the dynamic (temporal) component of the event stream, and not only the spatial component, at every moment of time. This is done by approximating the 3D geometry of the event stream with a parametric model; as a result, the algorithm is capable of producing the motion-compensated event stream (effectively approximating egomotion), and without using any form of external sensors in extremely low-light and noisy conditions without any form of feature tracking or explicit optical flow computation. We demonstrate our framework on the task of independent motion detection and tracking, where we use the temporal model inconsistencies to locate differently moving objects in challenging situations of very fast motion.


title: A note on non-inclusion "interval" root solvers

abstract: In this note we give comments on recent interval extensions of the King-Ostrowski multipoint methods for solving nonlinear equations published in at least five scientific journals offering faulty "self-validated" methods. We show that they do not possess essential feature of genuine interval methods -- inclusion property.


title: Well-posedness of Hamilton-Jacobi equations in population dynamics and applications to large deviations

abstract: We prove Freidlin-Wentzell type large deviation principles for various rescaled models in populations dynamics that have immigration and possibly harvesting: birth-death processes, Galton-Watson trees, epidemic SI models, and prey-predator models. The proofs are carried out using a general analytic approach based on the well-posedness of a class of associated Hamilton-Jacobi equations. The notable feature for these Hamilton-Jacobi equations is that the Hamiltonian can be discontinuous at the boundary. We prove a well-posedness result for a large class of Hamilton-Jacobi equations corresponding to one-dimensional models, and give partial results for the multi-dimensional setting.


title: Kepler Data Validation I -- Architecture, Diagnostic Tests, and Data Products for Vetting Transiting Planet Candidates

abstract: The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for ~200,000 individual stellar targets in at least one observing quarter in the primary Kepler Mission. Light curves for target stars are extracted in the Kepler Science Data Processing Pipeline, and are searched for transiting planet signatures. A Threshold Crossing Event is generated in the transit search for targets where the transit detection threshold is exceeded and transit consistency checks are satisfied. These targets are subjected to further scrutiny in the Data Validation (DV) component of the Pipeline. Transiting planet candidates are characterized in DV, and light curves are searched for additional planets after transit signatures are modeled and removed. A suite of diagnostic tests is performed on all candidates to aid in discrimination between genuine transiting planets and instrumental or astrophysical false positives. Data products are generated per target and planet candidate to document and display transiting planet model fit and diagnostic test results. These products are exported to the Exoplanet Archive at the NASA Exoplanet Science Institute, and are available to the community. We describe the DV architecture and diagnostic tests, and provide a brief overview of the data products. Transiting planet modeling and the search for multiple planets on individual targets are described in a companion paper. The final revision of the Kepler Pipeline code base is available to the general public through GitHub. The Kepler Pipeline has also been modified to support the TESS Mission which will commence in 2018.


title: Rise of the First Super-Massive Stars

abstract: We use high resolution adaptive mesh refinement simulations to model the formation of massive metal-free stars in the early Universe. By applying Lyman-Werner (LW) backgrounds of 100 J$_{21}$ and 1000 J$_{21}$ respectively we construct environments conducive to the formation of massive stars. We find that only in the case of the higher LW backgrounds are super-critical accretion rates realised that are necessary for super-massive star (SMS) formation. Mild fragmentation is observed for both backgrounds. Violent dynamical interactions between the stars that form in the more massive halo formed (1000 J$_{21}$ background) results in the eventual expulsion of the two most massive stars from the halo. In the smaller mass halo (100 J$_{21}$ background) mergers of stars occur before any multibody interactions and a single massive Pop III star is left at the centre of the halo at the end of our simulation. Feedback from the very massive Pop III stars is not effective in generating a large HII region with ionising photons absorbed within a few thousand AU of the star. In all cases a massive black hole seed is the expected final fate of the most massive objects. The seed of the massive Pop III star which remained at the centre of the less massive halo experiences steady accretion rates of almost $10^{-2}$ M$_{\odot}$/yr and if these rates continue could potentially experience super-Eddington accretion rates in the immediate aftermath of collapsing into a black hole.


title: On sufficient conditions for mixed monotonicity

abstract: Mixed monotone systems form an important class of nonlinear systems that have recently received attention in the abstraction-based control design area. Slightly different definitions exist in the literature, and it remains a challenge to verify mixed monotonicity of a system in general. In this paper, we first clarify the relation between different existing definitions of mixed monotone systems, and then give two sufficient conditions for mixed monotone functions defined on Euclidean space. These sufficient conditions are more general than the ones from the existing control literature, and they suggest that mixed monotonicity is a very generic property. Some discussions are provided on the computational usefulness of the proposed sufficient conditions.


title: The $r$-derangement numbers

abstract: The classical derangement numbers count fixed point-free permutations. In this paper we study the enumeration problem of generalized derangements, when some of the elements are restricted to be in distinct cycles in the cycle decomposition. We find exact formula, combinatorial relations for these numbers as well as analytic and asymptotic description. Moreover, we study deeper number theoretical properties, like modularity, $p$-adic valuations, and diophantine problems.


title: Phase transitions and magnetization of the mixed-spin Ising-Heisenberg double sawtooth frustrated ladder

abstract: The mixed spin-(1,1/2) Ising-Heisenberg double sawtooth ladder containing mixture of both spin-1 and spin-1/2 nodal atoms, and the spin-1/2 interstitial dimers is approximately solved by the transfer-matrix method. Here, we study in detail the ground-state phase diagrams, also influences of the bilinear exchange coupling on the rungs and cyclic four-spin exchange interaction in square plaquette of each block on the magnetization and magnetic susceptibility of the suggested ladder at low temperature. Such a double sawtooth ladder may be found in a Shastry-Sutherland Lattice-type. In spite of odd and even blocks spin ordering are different from each other, but due to the commutation relation between all different block Hamiltonians, phase diagrams, magnetization behavior and thermodynamic properties of the model are the same for odd and even blocks. We show that at low temperature, both exchange couplings can change the quality and quantity of the magnetization plateaus versus the magnetic field changes. Specially, we find a new magnetization plateau M/Ms = 5/6 for this model. Besides, we examine the magnetic susceptibility and specific heat of the model in detail. It is proven that behaviors of the magnetization and the magnetic susceptibility coincide at low temperature. The specific heat displays diverse temperature dependencies, which include a Schottky-type peak at a special temperature interval. We observe that with increase of the bilinear exchange coupling on the rungs, second peak temperature dependence grows.


title: Twisted argyle quivers and Higgs bundles

abstract: Ordinarily, quiver varieties are constructed as moduli spaces of quiver representations in the category of vector spaces. It is also natural to consider quiver representations in a richer category, namely that of vector bundles on some complex variety equipped with a fixed sheaf that twists the morphisms. Representations of A-type quivers in this twisted category --- known in the literature as "holomorphic chains" --- have practical use in questions concerning the topology of the moduli space of Higgs bundles. In that problem, the variety is a Riemann surface of genus at least 2, and the twist is its canonical line bundle. We extend the treatment of twisted A-type quiver representations to any genus using the Hitchin stability condition induced by Higgs bundles and computing their deformation theory. We then focus in particular on so-called "argyle quivers", where the rank labelling alternates between 1 and integers $r_i\geq1$. We give explicit geometric identifications of moduli spaces of twisted representations of argyle quivers on $\mathbb{P}^1$ using invariant theory for a non-reductive action via Euclidean reduction on polynomials. This leads to a stratification of the moduli space by change of bundle type, which we identify with "collision manifolds" of invariant zeroes of polynomials. We also relate the present work to Bradlow-Daskalopoulos stability and Thaddeus' pullback maps to stable tuples. We apply our results to computing $\mathbb{Q}$-Betti numbers of low-rank twisted Higgs bundle moduli spaces on $\mathbb{P}^1$, where the Higgs fields take values in an arbitrary ample line bundle. Our results agree with conjectural Poincar\'e series arising from the ADHM recursion formula.


title: Minimising the expectation value of the procurement cost in electricity markets based on the prediction error of energy consumption

abstract: In this paper, we formulate a method for minimising the expectation value of the procurement cost of electricity in two popular spot markets: {\it day-ahead} and {\it intra-day}, under the assumption that expectation value of unit prices and the distributions of prediction errors for the electricity demand traded in two markets are known. The expectation value of the total electricity cost is minimised over two parameters that change the amounts of electricity. Two parameters depend only on the expected unit prices of electricity and the distributions of prediction errors for the electricity demand traded in two markets. That is, even if we do not know the predictions for the electricity demand, we can determine the values of two parameters that minimise the expectation value of the procurement cost of electricity in two popular spot markets. We demonstrate numerically that the estimate of two parameters often results in a small variance of the total electricity cost, and illustrate the usefulness of the proposed procurement method through the analysis of actual data.


title: A note on $p$-adic locally analytic functions with application to behavior of the $p$-adic valuations of Stirling numbers

abstract: The aim of this paper is to prove conjectures concerning $p$-adic valuations of Stirling numbers of the second kind $S(n,k)$, $n,k\in\mathbb{N}_+$, stated by Amdeberhan, Manna and Moll and Berrizbeitia et al., where $p$ is a prime number. The proof is based on elementary facts from $p$-adic analysis.


title: Anomaly Detection in Road Networks Using Sliding-Window Tensor Factorization

abstract: Anomaly detection in road networks is vital for traffic management and emergency response. However, existing approaches do not directly address multiple anomaly types. We propose a tensor-based spatio-temporal model for detecting multiple types of anomalies in road networks. First, we represent network traffic data as a 3rd-order tensor. Next, we acquire spatial and multi-scale temporal patterns of traffic variations via a novel, computationally efficient tensor factorization algorithm: sliding window tensor factorization. Then, from the factorization results, we can identify different anomaly types by measuring deviations from different spatial and temporal patterns. Finally, we discover path-level anomalies by formulating anomalous path inference as a linear program that solves for the best matched paths of anomalous links. We evaluate the proposed methods via both synthetic experiments and case studies based on a real-world vehicle trajectory dataset, demonstrating advantages of our approach over baselines.


title: Computational methods in cardiovascular mechanics

abstract: The introduction of computational models in cardiovascular sciences has been progressively bringing new and unique tools for the investigation of the physiopathology. Together with the dramatic improvement of imaging and measuring devices on one side, and of computational architectures on the other one, mathematical and numerical models have provided a new, clearly noninvasive, approach for understanding not only basic mechanisms but also patient-specific conditions, and for supporting the design and the development of new therapeutic options. The terminology in silico is, nowadays, commonly accepted for indicating this new source of knowledge added to traditional in vitro and in vivo investigations. The advantages of in silico methodologies are basically the low cost in terms of infrastructures and facilities, the reduced invasiveness and, in general, the intrinsic predictive capabilities based on the use of mathematical models. The disadvantages are generally identified in the distance between the real cases and their virtual counterpart required by the conceptual modeling that can be detrimental for the reliability of numerical simulations.


title: Low-cost orbital-based linear-scaling \emph{ab initio} molecular dynamics for weakly-interacting systems

abstract: Within the framework of linear-scaling Kohn-Sham density functional theory, a robust method for maintaining compact localized orbitals close to the ground state is coupled with nuclear dynamics. This allows to obviate the commonly employed optimization of the one-electron density matrix and thus create an efficient orbital-only molecular dynamics method for weakly-interacting systems. An application to liquid water demonstrates that the low computational overhead of the method makes it well-suited for routine simulations whereas its linear-scaling complexity allows to extend first-principle dynamical studies of molecular systems to previously inaccessible length scales.


title: Massive Multiple Input Massive Multiple Output for 5G Wireless Backhauling

abstract: In this paper, we propose a new technique for the future fifth generation cellular network wireless backhauling. We show that hundreds of bits per second per Hertz (bits per second per Hz) of spectral efficiency can be attained at a high carrier frequency (such as 26 GHz) between large antenna arrays deployed along structures (such as lamp posts) that are close and roughly parallel to each other. Hundreds of data streams are spatially multiplexed through a short range and line of sight massive multiple input massive multiple output propagation channel thanks to a new low complexity spatial multiplexing scheme, called block discrete Fourier transform based spatial multiplexing with maximum ratio transmission. Its performance in real and existing environments is assessed using accurate ray-tracing tools and antenna models. In the best simulated scenario, 1.6 kbits per second per Hz of spectral efficiency is attained, corresponding to 80% of Singular Value Decomposition performance, with a transmitter and a receiver that are 200 and 10,000 times less complex, respectively.


title: Detection of a population of carbon-enhanced metal-poor stars in the Sculptor dwarf spheroidal galaxy

abstract: The study of the chemical abundances of metal-poor stars in dwarf galaxies provides a venue to constrain paradigms of chemical enrichment and galaxy formation. Here we present metallicity and carbon abundance measurements of 100 stars in Sculptor from medium-resolution (R ~ 2000) spectra taken with the Magellan/Michigan Fiber System mounted on the Magellan-Clay 6.5m telescope at Las Campanas Observatory. We identify 24 extremely metal-poor star candidates ([Fe/H] < -3.0) and 21 carbon-enhanced metal-poor (CEMP) star candidates. Eight carbon-enhanced stars are classified with at least 2$\sigma$ confidence and five are confirmed as such with follow-up R~6000 observations using the Magellan Echellette Spectrograph on the Magellan-Baade 6.5m telescope. We measure a CEMP fraction of 36% for stars below [Fe/H] = -3.0, indicating that the prevalence of carbon-enhanced stars in Sculptor is similar to that of the halo (~43%) after excluding likely CEMP-s and CEMP-r/s stars from our sample. However, we do not detect that any CEMP stars are strongly enhanced in carbon (e.g., [C/Fe] > 1.0). The existence of a large number of CEMP stars both in the halo and in Sculptor suggests that some halo CEMP stars may have originated from accreted early analogs of dwarf galaxies.


title: Modular Linear Optical Circuits

abstract: We propose and demonstrate a modular architecture for reconfigurable on-chip linear-optical circuits. Each module contains 10 independent phase-controlled Mach-Zehnder interferometers; several such modules can be connected to each other to build large reconfigurable interferometers. With this architecture, large interferometers are easier to build and characterize than with traditional, bespoke, monolithic designs. We demonstrate our approach by fabricating three modules in the form of UV-written silica-on-silicon chips. We characterize these chips, connect them to each other, and implement a wide range of linear optical transformations. We envisage that this architecture will enable many future experiments in quantum optics.


title: Constraints on light Dark Matter fermions from relic density consideration and Tsallis statistics

abstract: The cold dark matter fermions with mass MeV scale, pair produced inside the supernova SN1987A core, can freely stream away from the supernovae and hence contributes to its energy loss rate. Similar type of DM fermions(having similar kind of coupling to the standard model photon), produced from some other sources earlier, could have contributed to the relic density of the Universe. Working in a theory with an effective dark matter-photon coupling (inversely proportional to the scale $\Lambda$) in the formalism of Tsallis statistics, we find the dark matter contribution to the relic density and obtain a upper bound on $\Lambda$ using the experimental bound on the relic density for cold non-baryonic matter i.e. $\Omega h^2 = 0.1186 \pm 0.0020 $. The upper bound obtained from the relic density is shown with the lower bound obtained from the Raffelt's criterion on the emissibity rate of the supernovae SN1987A energy loss $\dot{\varepsilon}(e^+ e^- \to \chi \overline{\chi}) \le 10^{19}~\rm{erg~g^{-1}s^{-1}}$ and the optical depth criteria on the free streaming of the dark matter fermion (produced inside the supernovae core). As the deformation parameter $q$ changes from $1.0$ (undeformed scenario) to $1.1$(deformed scenario), the relic density bound on $\Lambda$ is found to vary from $ \sim 4.9 \times 10^7 $ TeV to $1.6 \times 10^8$ TeV for a fermion dark matter($\chi$) of mass $m_\chi = 30~\rm{MeV}$, which is almost $10$ times more than the lower bound obtained from the SN1987A energy loss rate and the optical depth criteria. \noindent {{\bf Keywords}: Dark matter, Relic density, Supernova cooling, Tsallis statistics, free-streaming, } }


title: Probing secret interactions of eV-scale sterile neutrinos with the diffuse supernova neutrino background

abstract: Sterile neutrinos with mass in the eV-scale and large mixings of order $\theta_0\simeq 0.1$ could explain some anomalies found in short-baseline neutrino oscillation data. Here, we revisit a neutrino portal scenario in which eV-scale sterile neutrinos have self-interactions via a new gauge vector boson $\phi$. Their production in the early Universe via mixing with active neutrinos can be suppressed by the induced effective potential in the sterile sector. We study how different cosmological observations can constrain this model, in terms of the mass of the new gauge boson, $M_\phi$, and its coupling to sterile neutrinos, $g_s$. Then, we explore how to probe part of the allowed parameter space of this particular model with future observations of the diffuse supernova neutrino background by the Hyper-Kamiokande and DUNE detectors. For $M_\phi \sim 5-10$~keV and $g_s \sim 10^{-4}-10^{-2}$, as allowed by cosmological constraints, we find that interactions of diffuse supernova neutrinos with relic sterile neutrinos on their way to the Earth would result in significant dips in the neutrino spectrum which would produce unique features in the event spectra observed in these detectors.


title: Symmetries of Chimera States

abstract: Symmetry broken states arise naturally in oscillatory networks. In this Letter, we investigate chaotic attractors in an ensemble of four mean-coupled Stuart-Landau oscillators with two oscillators being synchronized. We report that these states with partially broken symmetry, so-called chimera states, have different set-wise symmetries in the incoherent oscillators, and in particular some are and some are not invariant under a permutation symmetry on average. This allows for a classification of different chimera states in small networks. We conclude our report with a discussion of related states in spatially extended systems, which seem to inherit the symmetry properties of their counterparts in small networks.


title: Possible Bright Starspots on TRAPPIST-1

abstract: The M8V star TRAPPIST-1 hosts seven roughly Earth-sized planets and is a promising target for exoplanet characterization. Kepler/K2 Campaign 12 observations of TRAPPIST-1 in the optical show an apparent rotational modulation with a 3.3 day period, though that rotational signal is not readily detected in the Spitzer light curve at 4.5 $\mu$m. If the rotational modulation is due to starspots, persistent dark spots can be excluded from the lack of photometric variability in the Spitzer light curve. We construct a photometric model for rotational modulation due to photospheric bright spots on TRAPPIST-1 which is consistent with both the Kepler and Spitzer light curves. The maximum-likelihood model with three spots has typical spot sizes of $R_\mathrm{spot}/R_\star \approx 0.004$ at temperature $T_\mathrm{spot} \gtrsim 5300 \pm 200$ K. We also find that large flares are observed more often when the brightest spot is facing the observer, suggesting a correlation between the position of the bright spots and flare events. In addition, these flares may occur preferentially when the spots are increasing in brightness, which suggests that the 3.3 d periodicity may not be a rotational signal, but rather a characteristic timescale of active regions.


title: Optimal $H_2$ Decentralized Control of Cone Causal Spatially Invariant Systems

abstract: This paper presents an explicit solution to decentralized control of a class of spatially invariant systems. The problem of optimal $H_2$ decentralized control for cone causal systems is formulated. Using Parseval's identity, the optimal $H_2$ decentralized control problem is transformed into an infinite number of model matching problems with a specific structure that can be solved efficiently. In addition, the closed-form expression (explicit formula) of the decentralized controller is derived for the first time. In particular, it is shown that the optimal decentralized controller is given by a specific positive feedback scheme. A constructive procedure to obtain the state-space representation of the decentralized controller is provided. A numerical example is given and compared with previous works which demonstrate the effectiveness of the proposed method.


title: Coverings of rational ruled normal surfaces

abstract: In this work we use arithmetic, geometric, and combinatorial techniques to compute the cohomology of Weil divisors of a special class of normal surfaces, the so-called rational ruled toric surfaces. These computations are used to study the topology of cyclic coverings of such surfaces ramified along Q-normal crossing divisors.


title: The topological biquandle of a link

abstract: To every oriented link $L$, we associate a topologically defined biquandle $\widehat{\mathcal{B}}_{L}$, which we call the topological biquandle of $L$. The construction of $\widehat{\mathcal{B}}_{L}$ is similar to the topological description of the fundamental quandle given by Matveev. We find a presentation of the topological biquandle and explain how it is related to the fundamental biquandle of the link.


title: Analysis of spectral clustering algorithms for community detection: the general bipartite setting

abstract: We consider spectral clustering algorithms for community detection under a general bipartite stochastic block model (SBM). A modern spectral clustering algorithm consists of three steps: (1) regularization of an appropriate adjacency or Laplacian matrix (2) a form of spectral truncation and (3) a k-means type algorithm in the reduced spectral domain. We focus on the adjacency-based spectral clustering and for the first step, propose a new data-driven regularization that can restore the concentration of the adjacency matrix even for the sparse networks. This result is based on recent work on regularization of random binary matrices, but avoids using unknown population level parameters, and instead estimates the necessary quantities from the data. We also propose and study a novel variation of the spectral truncation step and show how this variation changes the nature of the misclassification rate in a general SBM. We then show how the consistency results can be extended to models beyond SBMs, such as inhomogeneous random graph models with approximate clusters, including a graphon clustering problem, as well as general sub-Gaussian biclustering. A theme of the paper is providing a better understanding of the analysis of spectral methods for community detection and establishing consistency results, under fairly general clustering models and for a wide regime of degree growths, including sparse cases where the average expected degree grows arbitrarily slowly.


title: Taking Turing by Surprise? Designing Digital Computers for morally-loaded contexts

abstract: There is much to learn from what Turing hastily dismissed as Lady Lovelace s objection. Digital computers can indeed surprise us. Just like a piece of art, algorithms can be designed in such a way as to lead us to question our understanding of the world, or our place within it. Some humans do lose the capacity to be surprised in that way. It might be fear, or it might be the comfort of ideological certainties. As lazy normative animals, we do need to be able to rely on authorities to simplify our reasoning: that is ok. Yet the growing sophistication of systems designed to free us from the constraints of normative engagement may take us past a point of no-return. What if, through lack of normative exercise, our moral muscles became so atrophied as to leave us unable to question our social practices? This paper makes two distinct normative claims: 1. Decision-support systems should be designed with a view to regularly jolting us out of our moral torpor. 2. Without the depth of habit to somatically anchor model certainty, a computer s experience of something new is very different from that which in humans gives rise to non-trivial surprises. This asymmetry has key repercussions when it comes to the shape of ethical agency in artificial moral agents. The worry is not just that they would be likely to leap morally ahead of us, unencumbered by habits. The main reason to doubt that the moral trajectories of humans v. autonomous systems might remain compatible stems from the asymmetry in the mechanisms underlying moral change. Whereas in humans surprises will continue to play an important role in waking us to the need for moral change, cognitive processes will rule when it comes to machines. This asymmetry will translate into increasingly different moral outlooks, to the point of likely unintelligibility. The latter prospect is enough to doubt the desirability of autonomous moral agents.


title: Supersymmetrising the GSY Soliton

abstract: We supersymmetrise the Hopfion studied in a previous work. This soliton represents a closed semilocal vortex string in $U(1)$ gauge theory. It carries nonzero Hopf number due to the additional winding of a phase modulus as one moves along the closed string. We study this solution in $\mathcal{N}= 2$ supersymmetric QED with two flavours. As a preliminary exercise we compactify one space dimension and consider a straight vortex with periodic boundary conditions. It turns out to be 1/2-BPS saturated. An additional winding along the string can be introduced and it does not spoil the BPS nature of the object. Next, we consider a ring-like vortex in a non-compact space and show that the circumference of the ring $L$ can be stabilised once the previously mentioned winding along the string is introduced. Of course the ring-like vortex is not BPS but its energy becomes close to the BPS bound if $L$ is large, which can be guaranteed in the case that we have a large value of the angular momentum $J$. Thus we arrive at the concept of asymptotically BPS-saturated solitons. BPS saturation is achieved in the limit $J\rightarrow \infty$.


title: Ergodicity in Stationary Graph Processes: A Weak Law of Large Numbers

abstract: For stationary signals in time the weak law of large numbers (WLLN) states that ensemble and realization averages are within e of each other with a probability of order O(1/Ne^2) when considering N signal components. The graph WLLN introduced in this paper shows that the same is essentially true for signals supported on graphs. However, the notions of stationarity, ensemble mean, and realization mean are different. Recent papers have defined graph stationary signals as those that satisfy a form of invariance with respect to graph diffusion. The ensemble mean of a graph stationary signal is not a constant but a node-varying signal whose structure depends on the spectral properties of the graph. The realization average of a graph signal is defined here as an average of successive weighted averages of local signal values with signal values of neighboring nodes. The graph WLLN shows that this two node-varying signals are within e of each other with probability of order O(1/Ne^2) in at least some nodes. In stationary time signals, the realization average is not only a consistent estimator of the ensemble mean but also optimal in terms of mean squared error (MSE). This is not true of graph signals. Optimal MSE graph filter designs are also presented. An example problem concerning the estimation of the mean of a Gaussian random field is presented.


title: Multi-Sensor Conflict Measurement and Information Fusion

abstract: In sensing applications where multiple sensors observe the same scene, fusing sensor outputs can provide improved results. However, if some of the sensors are providing lower quality outputs, the fused results can be degraded. In this work, a multi-sensor conflict measure is proposed which estimates multi-sensor conflict by representing each sensor output as interval-valued information and examines the sensor output overlaps on all possible n-tuple sensor combinations. The conflict is based on the sizes of the intervals and how many sensors output values lie in these intervals. In this work, conflict is defined in terms of how little the output from multiple sensors overlap. That is, high degrees of overlap mean low sensor conflict, while low degrees of overlap mean high conflict. This work is a preliminary step towards a robust conflict and sensor fusion framework. In addition, a sensor fusion algorithm is proposed based on a weighted sum of sensor outputs, where the weights for each sensor diminish as the conflict measure increases. The proposed methods can be utilized to (1) assess a measure of multi-sensor conflict, and (2) improve sensor output fusion by lessening weighting for sensors with high conflict. Using this measure, a simulated example is given to explain the mechanics of calculating the conflict measure, and stereo camera 3D outputs are analyzed and fused. In the stereo camera case, the sensor output is corrupted by additive impulse noise, DC offset, and Gaussian noise. Impulse noise is common in sensors due to intermittent interference, a DC offset a sensor bias or registration error, and Gaussian noise represents a sensor output with low SNR. The results show that sensor output fusion based on the conflict measure shows improved accuracy over a simple averaging fusion strategy.


title: Text Data Mining from the Author's Perspective: Whose Text, Whose Mining, and to Whose Benefit?

abstract: Given the many technical, social, and policy shifts in access to scholarly content since the early days of text data mining, it is time to expand the conversation about text data mining from concerns of the researcher wishing to mine data to include concerns of researcher-authors about how their data are mined, by whom, for what purposes, and to whose benefits.


title: Luby--Veli\v{c}kovi\'c--Wigderson revisited: Improved correlation bounds and pseudorandom generators for depth-two circuits

abstract: We study correlation bounds and pseudorandom generators for depth-two circuits that consist of a $\mathsf{SYM}$-gate (computing an arbitrary symmetric function) or $\mathsf{THR}$-gate (computing an arbitrary linear threshold function) that is fed by $S$ $\mathsf{AND}$ gates. Such circuits were considered in early influential work on unconditional derandomization of Luby, Veli\v{c}kovi\'c, and Wigderson [LVW93], who gave the first non-trivial PRG with seed length $2^{O(\sqrt{\log(S/\varepsilon)})}$ that $\varepsilon$-fools these circuits. In this work we obtain the first strict improvement of [LVW93]'s seed length: we construct a PRG that $\varepsilon$-fools size-$S$ $\{\mathsf{SYM},\mathsf{THR}\} \circ\mathsf{AND}$ circuits over $\{0,1\}^n$ with seed length \[ 2^{O(\sqrt{\log S })} + \mathrm{polylog}(1/\varepsilon), \] an exponential (and near-optimal) improvement of the $\varepsilon$-dependence of [LVW93]. The above PRG is actually a special case of a more general PRG which we establish for constant-depth circuits containing multiple $\mathsf{SYM}$ or $\mathsf{THR}$ gates, including as a special case $\{\mathsf{SYM},\mathsf{THR}\} \circ \mathsf{AC^0}$ circuits. These more general results strengthen previous results of Viola [Vio06] and essentially strengthen more recent results of Lovett and Srinivasan [LS11]. Our improved PRGs follow from improved correlation bounds, which are transformed into PRGs via the Nisan--Wigderson "hardness versus randomness" paradigm [NW94]. The key to our improved correlation bounds is the use of a recent powerful \emph{multi-switching} lemma due to H{\aa}stad [H{\aa}s14].


title: Reconstruction Algorithm Design for Mitigating the Orientation Dependent Conspicuity of Fiber-Like signals in Digital Breast Tomosynthesis

abstract: There are a number of clinically relevant tasks in digital breast tomosynthesis (DBT) involving the detection and visual assessment of fiber-like structures such as Cooper's ligaments, blood vessels, and spiculated lesions. Such structures can exhibit orientation dependent variations in conspicuity. This study demonstrates the presence of in-plane orientation-dependent signal conspicuity for fiber-like signals in DBT and shows how reconstruction algorithm design can mitigate this phenomenon. We uncover a tradeoff between minimizing orientation-dependence and preserving depth resolution that is dictated by the regularization strength employed in reconstruction.


title: Reactive Proximity Data Structures for Graphs

abstract: We consider data structures for graphs where we maintain a subset of the nodes called sites, and allow proximity queries, such as asking for the closest site to a query node, and update operations that enable or disable nodes as sites. We refer to a data structure that can efficiently react to such updates as reactive. We present novel reactive proximity data structures for graphs of polynomial expansion, i.e., the class of graphs with small separators, such as planar graphs and road networks. Our data structures can be used in several logistical problems and geographic information systems dealing with real-time data, such as emergency dispatching. We experimentally compare our data structure to Dijkstra's algorithm in a system emulating random queries in a real road network.


title: Measuring Conflict in a Multi-Source Environment as a Normal Measure

abstract: In a multi-source environment, each source has its own credibility. If there is no external knowledge about credibility then we can use the information provided by the sources to assess their credibility. In this paper, we propose a way to measure conflict in a multi-source environment as a normal measure. We examine our algorithm using three simulated examples of increasing conflict and one experimental example. The results demonstrate that the proposed measure can represent conflict in a meaningful way similar to what a human might expect and from it we can identify conflict within our sources.


title: Topological horseshoes for surface homeomorphisms

abstract: In this work we develop a new criterion for the existence of topological horseshoes for surface homeomorphisms in the isotopy class of the identity. Based on our previous work on forcing theory, this new criterion is purely topological and can be expressed in terms of equivariant Brouwer foliations and transverse trajectories. We then apply this new tool in the study of the dynamics of homeomorphisms of surfaces with zero genus and null topological entropy and we obtain several applications. For homeomorphisms of the open annulus $\mathbb{A}$ with zero topological entropy, we show that rotation numbers exists for all points with nonempty omega limits, and that if $\mathbb{A}$ is a generalized region of instability then it admits a single rotation vector. We also offer a new proof of a recent result of Passegi, Potrie and Sambarino, showing that zero entropy dissipative homeomorphisms of the annulus having as an atractor a circloid have a single rotation number. Our work also studies homeomorphisms of the sphere without horseshoes. For these maps we present a structure theorem in terms of fixed point free invariant sub-annuli, as well as a very restricted description of all possible dynamical behavior in the transitive subsets. This description ensures, for instance, that transitive sets can contain at most $2$ distinct periodic orbits and that, in many cases, the restriction of the homeomorphism to the transitive set must be an extension of an odometer. In particular, we show that any nontrivial and stable transitive subset of a dissipative diffeomorphism of the plane is always infinitely renormalizable in the sense of Bonatti-Gambaudo-Lion-Tresser.


title: Empirical Wavelet-based Estimation for Non-linear Additive Regression Models

abstract: Additive regression models are actively researched in the statistical field because of their usefulness in the analysis of responses determined by non-linear relationships with multivariate predictors. In this kind of statistical models, the response depends linearly on unknown functions of predictor variables and typically, the goal of the analysis is to make inference about these functions. In this paper, we consider the problem of Additive Regression with random designs from a novel viewpoint: we propose an estimator based on an orthogonal projection onto a multiresolution space using empirical wavelet coefficients that are fully data driven. In this setting, we derive a mean-square consistent estimator based on periodic wavelets on the interval $[0,1]$. For construction of the estimator, we assume that the joint distribution of predictors is non-zero and bounded on its support; We also assume that the functions belong to a Sobolev space and integrate to zero over the $[0,1]$ interval, which guarantees model identifiability and convergence of the proposed method. Moreover, we provide the $\mathbb{L}_{2}$ risk analysis of the estimator and derive its convergence rate. Theoretically, we show that this approach achieves good convergence rates when the dimensionality of the problem is relatively low and the set of unknown functions is sufficiently smooth. In this approach, the results are obtained without the assumption of an equispaced design, a condition that is typically assumed in most wavelet-based procedures. Finally, we show practical results obtained from simulated data, demonstrating the potential applicability of our method in the problem of additive regression models with random designs.


title: Weighted Bayesian Bootstrap for Scalable Bayes

abstract: We develop a weighted Bayesian Bootstrap (WBB) for machine learning and statistics. WBB provides uncertainty quantification by sampling from a high dimensional posterior distribution. WBB is computationally fast and scalable using only off-theshelf optimization software such as TensorFlow. We provide regularity conditions which apply to a wide range of machine learning and statistical models. We illustrate our methodology in regularized regression, trend filtering and deep learning. Finally, we conclude with directions for future research.


title: How far is the Borel map from being surjective in quasianalytic ultradifferentiable classes?

abstract: The Borel map takes germs at 0 of smooth functions to the sequence of iterated partial derivatives at 0. In the literature, it is well known that the restriction of this mapping to the germs of quasianalytic ultradifferentiable classes which are strictly containing the real analytic functions can never be onto the corresponding sequence space. In this paper, we are interested in studying how large the image of the Borel map is and we investigate the size and the structure of this image by using different approaches (Baire residuality, prevalence and lineability). We give an answer to this question in the very general setting of quasianalytic ultradifferentiable classes defined by weight matrices, which contains as particular cases the classes defined by a single weight sequence or by a weight function.


title: Optimizing power oscillations in an ellipsometric system

abstract: Ellipsometry is a powerful and well-established optical technique used in the characterisation of materials. It works by combining the components of elliptically polarized light in order to draw information about the optical system. We propose an ellipsometric experimental set up to study polarization interference in the total internal reflection regime for Gaussian laser beams. The relative phase between orthogonal states can be measured as a power oscillation of the optical beam transmitted through a dielectric block and whose orthogonal components are then mixed by a polarizer. We show under which conditions the plane wave analysis is valid and when the power oscillation can be optimized to reproduce a full pattern of oscillation and to simulate quarter and half wave plates.


title: Bias in OLAP Queries: Detection, Explanation, and Removal

abstract: On line analytical processing (OLAP) is an essential element of decision-support systems. OLAP tools provide insights and understanding needed for improved decision making. However, the answers to OLAP queries can be biased and lead to perplexing and incorrect insights. In this paper, we propose HypDB, a system to detect, explain, and to resolve bias in decision-support queries. We give a simple definition of a \emph{biased query}, which performs a set of independence tests on the data to detect bias. We propose a novel technique that gives explanations for bias, thus assisting an analyst in understanding what goes on. Additionally, we develop an automated method for rewriting a biased query into an unbiased query, which shows what the analyst intended to examine. In a thorough evaluation on several real datasets we show both the quality and the performance of our techniques, including the completely automatic discovery of the revolutionary insights from a famous 1973 discrimination case.


title: On The Nature of ultra-faint Dwarf Galaxy Candidates II: The case of Cetus II

abstract: We obtained deep Gemini GMOS-S $g,r$ photometry of the ultra-faint dwarf galaxy candidate Cetus II with the aim of providing stronger constraints on its size, luminosity and stellar population. Cetus II is an important object in the size-luminosity plane as it occupies the transition zone between dwarf galaxies and star clusters. All known objects smaller than Cetus II ($r_h \sim 20$ pc) are reported to be star clusters, while most larger objects are likely dwarf galaxies. We found a prominent excess of main-sequence stars in the colour-magnitude diagram of Cetus II, best described by a single stellar population with an age of 11.2 Gyr, metallicity of [Fe/H] = $-1.28$ dex, an [$\alpha$/Fe] = 0.0 dex at a heliocentric distance of 26.3$\pm$1.2 kpc. As well as being spatially located within the Sagittarius dwarf tidal stream, these properties are well matched to the Sagittarius galaxy's Population B stars. Interestingly, like our recent findings on the ultra-faint dwarf galaxy candidate Tucana V, the stellar field in the direction of Cetus II shows no evidence of a concentrated overdensity despite tracing the main sequence for over six magnitudes. These results strongly support the picture that Cetus II is not an ultra-faint stellar system in the Milky Way halo, but made up of stars from the Sagittarius tidal stream.


title: NASA's Asteroid Grand Challenge: Strategy, Results and Lessons Learned

abstract: Beginning in 2012, NASA utilized a strategic process to identify broad societal questions, or grand challenges, that are well suited to the aerospace sector and align with national priorities. This effort generated NASA's first grand challenge, the Asteroid Grand Challenge, a large scale effort using multidisciplinary collaborations and innovative engagement mechanisms focused on finding and addressing asteroid threats to human populations. In April 2010, President Barack Obama announced a mission to send humans to an asteroid by 2025. This resulted in the agency's Asteroid Redirect Mission to leverage and maximize existing robotic and human efforts to capture and reroute an asteroid, with the goal of eventual human exploration. The AGC, initiated in 2013, complemented ARM by expanding public participation, partnerships, and other approaches to find, understand, and overcome these potentially harmful asteroids. This paper describes a selection of AGC activities implemented from 2013 to 2017 and their results, excluding those conducted by NASA's Near Earth Object Observations Program and other organizations. The strategic development of the initiative is outlined as well as initial successes, strengths, and weaknesses resulting from the first four years of AGC activities and approaches. Finally, we describe lesson learned and areas for continued work and study. The AGC lessons learned and strategies could inform the work of other agencies and organizations seeking to conduct a global scientific investigation with matrixed organizational support, multiple strategic partners, and numerous internal and external open innovation approaches and audiences.


title: Learning to recognize Abnormalities in Chest X-Rays with Location-Aware Dense Networks

abstract: Chest X-ray is the most common medical imaging exam used to assess multiple pathologies. Automated algorithms and tools have the potential to support the reading workflow, improve efficiency, and reduce reading errors. With the availability of large scale data sets, several methods have been proposed to classify pathologies on chest X-ray images. However, most methods report performance based on random image based splitting, ignoring the high probability of the same patient appearing in both training and test set. In addition, most methods fail to explicitly incorporate the spatial information of abnormalities or utilize the high resolution images. We propose a novel approach based on location aware Dense Networks (DNetLoc), whereby we incorporate both high-resolution image data and spatial information for abnormality classification. We evaluate our method on the largest data set reported in the community, containing a total of 86,876 patients and 297,541 chest X-ray images. We achieve (i) the best average AUC score for published training and test splits on the single benchmarking data set (ChestX-Ray14), and (ii) improved AUC scores when the pathology location information is explicitly used. To foster future research we demonstrate the limitations of the current benchmarking setup and provide new reference patient-wise splits for the used data sets. This could support consistent and meaningful benchmarking of future methods on the largest publicly available data sets.


title: Compact Convolutional Neural Networks for Classification of Asynchronous Steady-state Visual Evoked Potentials

abstract: Steady-State Visual Evoked Potentials (SSVEPs) are neural oscillations from the parietal and occipital regions of the brain that are evoked from flickering visual stimuli. SSVEPs are robust signals measurable in the electroencephalogram (EEG) and are commonly used in brain-computer interfaces (BCIs). However, methods for high-accuracy decoding of SSVEPs usually require hand-crafted approaches that leverage domain-specific knowledge of the stimulus signals, such as specific temporal frequencies in the visual stimuli and their relative spatial arrangement. When this knowledge is unavailable, such as when SSVEP signals are acquired asynchronously, such approaches tend to fail. In this paper, we show how a compact convolutional neural network (Compact-CNN), which only requires raw EEG signals for automatic feature extraction, can be used to decode signals from a 12-class SSVEP dataset without the need for any domain-specific knowledge or calibration data. We report across subject mean accuracy of approximately 80% (chance being 8.3%) and show this is substantially better than current state-of-the-art hand-crafted approaches using canonical correlation analysis (CCA) and Combined-CCA. Furthermore, we analyze our Compact-CNN to examine the underlying feature representation, discovering that the deep learner extracts additional phase and amplitude related features associated with the structure of the dataset. We discuss how our Compact-CNN shows promise for BCI applications that allow users to freely gaze/attend to any stimulus at any time (e.g., asynchronous BCI) as well as provides a method for analyzing SSVEP signals in a way that might augment our understanding about the basic processing in the visual cortex.


title: Convolutional Neural Networks and Language Embeddings for End-to-End Dialect Recognition

abstract: Dialect identification (DID) is a special case of general language identification (LID), but a more challenging problem due to the linguistic similarity between dialects. In this paper, we propose an end-to-end DID system and a Siamese neural network to extract language embeddings. We use both acoustic and linguistic features for the DID task on the Arabic dialectal speech dataset: Multi-Genre Broadcast 3 (MGB-3). The end-to-end DID system was trained using three kinds of acoustic features: Mel-Frequency Cepstral Coefficients (MFCCs), log Mel-scale Filter Bank energies (FBANK) and spectrogram energies. We also investigated a dataset augmentation approach to achieve robust performance with limited data resources. Our linguistic feature research focused on learning similarities and dissimilarities between dialects using the Siamese network, so that we can reduce feature dimensionality as well as improve DID performance. The best system using a single feature set achieves 73% accuracy, while a fusion system using multiple features yields 78% on the MGB-3 dialect test set consisting of 5 dialects. The experimental results indicate that FBANK features achieve slightly better results than MFCCs. Dataset augmentation via speed perturbation appears to add significant robustness to the system. Although the Siamese network with language embeddings did not achieve as good a result as the end-to-end DID system, the two approaches had good synergy when combined together in a fused system.


title: On $L^1$-estimates for probability solutions to Fokker-Planck-Kolmogorov equations

abstract: We prove two new results connected with elliptic Fokker-Planck-Kolmogorov equations with drifts integrable with respect to solutions. The first result answers negatively a long-standing question and shows that a density of a probability measure satisfying the Fokker-Planck-Kolmogorov equation with a drift integrable with respect to this density can fail to belong to the Sobolev class~$W^{1,1}(\mathbb{R}^d)$. There is also a version of this result for densities with respect to Gaussian measures. The second new result gives some positive information about properties of such solutions: the solution density is proved to belong to certain fractional Sobolev classes.


title: Colour fields of the quark-antiquark excited flux tube

abstract: We present colour field density profiles for some of the first SU(3) gluonic excitations of the flux tube in the presence of a static quark-antiquark pair. The results are obtained from a large set of gluonic operators.


title: Sharp inequalities for linear combinations of orthogonal martingales

abstract: For any two real-valued continuous-path martingales $X=\{X_t\}_{t\geq 0}$ and $Y=\{Y_t\}_{t\geq 0}$, with $X$ and $Y$ being orthogonal and $Y$ being differentially subordinate to $X$, we obtain sharp $L^p$ inequalities for martingales of the form $aX+bY$ with $a, b$ real numbers. The best $L^p$ constant is equal to the norm of the operator $aI+bH$ from $L^p$ to $L^p$, where $H$ is the Hilbert transform on the circle or real line. The values of these norms were found by Hollenbeck, Kalton and Verbitsky \cite{HKV}.


title: Laser Cooling of Optically Trapped Molecules

abstract: Calcium monofluoride (CaF) molecules are loaded into an optical dipole trap (ODT) and subsequently laser cooled within the trap. Starting with magneto-optical trapping, we sub-Doppler cool CaF and then load $150(30)$ CaF molecules into an ODT. Enhanced loading by a factor of five is obtained when sub-Doppler cooling light and trapping light are on simultaneously. For trapped molecules, we directly observe efficient sub-Doppler cooling to a temperature of $60(5)$ $\mu\text{K}$. The trapped molecular density of $8(2)\times10^7$ cm$^{-3}$ is an order of magnitude greater than in the initial sub-Doppler cooled sample. The trap lifetime of 750(40) ms is dominated by background gas collisions.


title: COPA: Constrained PARAFAC2 for Sparse & Large Datasets

abstract: PARAFAC2 has demonstrated success in modeling irregular tensors, where the tensor dimensions vary across one of the modes. An example scenario is modeling treatments across a set of patients with the varying number of medical encounters over time. Despite recent improvements on unconstrained PARAFAC2, its model factors are usually dense and sensitive to noise which limits their interpretability. As a result, the following open challenges remain: a) various modeling constraints, such as temporal smoothness, sparsity and non-negativity, are needed to be imposed for interpretable temporal modeling and b) a scalable approach is required to support those constraints efficiently for large datasets. To tackle these challenges, we propose a {\it CO}nstrained {\it PA}RAFAC2 (COPA) method, which carefully incorporates optimization constraints such as temporal smoothness, sparsity, and non-negativity in the resulting factors. To efficiently support all those constraints, COPA adopts a hybrid optimization framework using alternating optimization and alternating direction method of multiplier (AO-ADMM). As evaluated on large electronic health record (EHR) datasets with hundreds of thousands of patients, COPA achieves significant speedups (up to 36 times faster) over prior PARAFAC2 approaches that only attempt to handle a subset of the constraints that COPA enables. Overall, our method outperforms all the baselines attempting to handle a subset of the constraints in terms of speed, while achieving the same level of accuracy. Through a case study on temporal phenotyping of medically complex children, we demonstrate how the constraints imposed by COPA reveal concise phenotypes and meaningful temporal profiles of patients. The clinical interpretation of both the phenotypes and the temporal profiles was confirmed by a medical expert.


title: Identification of GSM and LTE Signals Using Their Second-order Cyclostationarity

abstract: Automatic signal identification (ASI) has various millitary and commercial applications, such as spectrum surveillance and cognitive radio. In this paper, a novel ASI algorithm is proposed for the identification of GSM and LTE signals, which is based on the pilot-induced second-order cyclostationarity. The proposed algorithm provides a very good performance at low signal-to-noise ratios and short observation times, with no need for channel estimation, and timing and frequency synchronization. Simulations and off-the-air signals acquired with the ThinkRF WSA4000 receiver are used to confirm the findings.


title: Boosting computational power through spatial multiplexing in quantum reservoir computing

abstract: Quantum reservoir computing provides a framework for exploiting the natural dynamics of quantum systems as a computational resource. It can implement real-time signal processing and solve temporal machine learning problems in general, which requires memory and nonlinear mapping of the recent input stream using the quantum dynamics in computational supremacy region, where the classical simulation of the system is intractable. A nuclear magnetic resonance spin-ensemble system is one of the realistic candidates for such physical implementations, which is currently available in laboratories. In this paper, considering these realistic experimental constraints for implementing the framework, we introduce a scheme, which we call a spatial multiplexing technique, to effectively boost the computational power of the platform. This technique exploits disjoint dynamics, which originate from multiple different quantum systems driven by common input streams in parallel. Accordingly, unlike designing a single large quantum system to increase the number of qubits for computational nodes, it is possible to prepare a huge number of qubits from multiple but small quantum systems, which are operationally easy to handle in laboratory experiments. We numerically demonstrate the effectiveness of the technique using several benchmark tasks and quantitatively investigate its specifications, range of validity, and limitations in detail.


title: Feasibility of rational shape design of single-polymer micelle using spontaneous surface curvature

abstract: Polymeric micelles are used in a variety of applications, with the micelle's shape often playing an important role. Consequently, a scheme to design micelles of arbitrary shape is desirable. In this paper, we consider micelles formed from a single, linear, multiblock copolymer, and we study how easily the micelle's shape can be controlled by altering the copolymer block lengths. Using a rational design scheme, we identify a few aspects of the multiblock composition that are expected to have a well-behaved, predictable effect on micelle shape. Starting from a reference micelle composition, itself already exhibiting a nonstandard shape having a moderately sized dimple, we alter these aspects of the multiblock composition and observe the regularity of the micelle shape response. The response of the shape is found to be somewhat smooth, but significantly nonlinear and sometimes nonmonotonic, suggesting that sophisticated techniques may be required to aid in micelle design.


title: Plasmonic Superconductivity in Layered Materials

abstract: Plasmonic excitations behave fundamentally different in layered materials in comparison to bulk systems. They form gapless modes, which in turn couple at low energies to the electrons. Thereby they can strongly influence superconducting instabilities. Here, we show how these excitations can be controlled from the outside via changes in the dielectric environment or in the doping level, which allows for external tuning of the superconducting transition temperature. By solving the gap equation for an effective system, we find that the plasmonic influence can both strongly enhance or reduce the transition temperature, depending on the details of the plasmon-phonon interplay. We formulate simple experimental guidelines to find plasmon- induced elevated transition temperatures in layered materials.


title: Bayesian optimization for computationally extensive probability distributions

abstract: An efficient method for finding a better maximizer of computationally extensive probability distributions is proposed on the basis of a Bayesian optimization technique. A key idea of the proposed method is to use extreme values of acquisition functions by Gaussian processes for the next training phase, which should be located near a local maximum or a global maximum of the probability distribution. Our Bayesian optimization technique is applied to the posterior distribution in the effective physical model estimation, which is a computationally extensive probability distribution. Even when the number of sampling points on the posterior distributions is fixed to be small, the Bayesian optimization provides a better maximizer of the posterior distributions in comparison to those by the random search method, the steepest descent method, or the Monte Carlo method. Furthermore, the Bayesian optimization improves the results efficiently by combining the steepest descent method and thus it is a powerful tool to search for a better maximizer of computationally extensive probability distributions.


title: Spanning Trees With Edge Conflicts and Wireless Connectivity

abstract: We introduce the problem of finding a spanning tree along with a partition of the tree edges into fewest number of feasible sets, where constraints on the edges define feasibility. The motivation comes from wireless networking, where we seek to model the irregularities seen in actual wireless environments. Not all node pairs may be able to communicate, even if geographically close --- thus, the available pairs are modeled with a link graph $\mathcal{L}=(V,E)$. Also, signal attenuation need not follow a nice geometric formulas --- hence, interference is modeled by a conflict (hyper)graph $\mathcal{C}=(E,F)$ on the links. The objective is to maximize the efficiency of the communication, or equivalently minimizing the length of a schedule of the tree edges in the form of a coloring. We find that in spite of all this generality, the problem can be approximated linearly in terms of a versatile parameter, the inductive independence of the interference graph. Specifically, we give a simple algorithm that attains a $O(\rho \log n)$-approximation, where $n$ is the number of nodes and $\rho$ is the inductive independence, and show that near-linear dependence on $\rho$ is also necessary. We also treat an extension to Steiner trees, modeling multicasting, and obtain a comparable result. Our results suggest that several canonical assumptions of geometry, regularity and "niceness" in wireless settings can sometimes be relaxed without a significant hit in algorithm performance.


title: It was the training data pruning too!

abstract: We study the current best model (KDG) for question answering on tabular data evaluated over the WikiTableQuestions dataset. Previous ablation studies performed against this model attributed the model's performance to certain aspects of its architecture. In this paper, we find that the model's performance also crucially depends on a certain pruning of the data used to train the model. Disabling the pruning step drops the accuracy of the model from 43.3% to 36.3%. The large impact on the performance of the KDG model suggests that the pruning may be a useful pre-processing step in training other semantic parsers as well.


title: Mechanisms of pyroelectricity in three- and two-dimensional materials

abstract: Pyroelectricity is a very promising phenomenon in three- and two-dimensional (2D) materials, but first-principles calculations have not so far been used to elucidate the underlying mechanisms. Here we report density-functional theory (DFT) calculations based on the Born-Szigeti theory of pyroelectricity, by combining fundamental thermodynamics and the modern theory of polarization. We find satisfactory agreement with experimental data in the case of bulk benchmark materials, showing that the so-called electron-phonon renormalization, whose contribution has been viewed as negligible, is important. We predict out-of-plane pyroelectricity in the recently synthesized Janus MoSSe monolayer and in-plane pyroelectricity in the group-IV monochalcogenide GeS monolayer. It is notable that the so-called secondary pyroelectricity is found to be dominant in GeS monolayer. The present work opens a theoretical route to study the pyroelectric effect using DFT and provides a valuable tool in the search for new candidates for pyroelectric applications.


title: Nonribbon 2-links all of whose components are trivial knots and some of whose band-sums are nonribbon knots

abstract: There is a nonribbon 2-link all of whose components are trivial 2-knots and one of whose band-sums is a nonribbon 2-knot.


title: Deep Bayesian Supervised Learning given Hypercuboidally-shaped, Discontinuous Data, using Compound Tensor-Variate & Scalar-Variate Gaussian Processes

abstract: We undertake Bayesian learning of the high-dimensional functional relationship between a system parameter vector and an observable, that is in general tensor-valued. The ultimate aim is Bayesian inverse prediction of the system parameters, at which test data is recorded. We attempt such learning given hypercuboidally-shaped data that displays strong discontinuities, rendering learning challenging. We model the sought high-dimensional function, with a tensor-variate Gaussian Process (GP), and use three independent ways for learning covariance matrices of the resulting likelihood, which is Tensor-Normal. We demonstrate that the discontinuous data demands that implemented covariance kernels be non-stationary--achieved by modelling each kernel hyperparameter, as a function of the sample function of the invoked tensor-variate GP. Each such function can be shown to be temporally-evolving, and treated as a realisation from a distinct scalar-variate GP, with covariance described adaptively by collating information from a historical set of samples of chosen sample-size. We prove that deep-learning using 2-"layers", suffice, where the outer-layer comprises the tensor-variate GP, compounded with multiple scalar-variate GPs in the "inner-layer", and undertake inference with Metropolis-within-Gibbs. We apply our method to a cuboidally-shaped, discontinuous, real dataset, and subsequently perform forward prediction to generate data from our model, given our results--to perform model-checking.


title: Arithmetic of curves on moduli of local systems

abstract: We investigate the arithmetic of algebraic curves on coarse moduli spaces for special linear rank two local systems on surfaces with fixed boundary traces. We prove a structure theorem for morphisms from the affine line into the moduli space. We show that the set of integral points on any nondegenerate algebraic curve on the moduli space can be effectively determined.


title: Sensitivity Based Thevenin Index for Voltage Stability Assessment Considering N-1 Contingency

abstract: This paper proposes an approach to address the voltage stability assessment (VSA) considering N-1 contingency. The approach leverages the sensitivity based Thevenin index (STI) which involves evaluating the Jacobian matrix at current operating condition. Since the N-1 contingency case is hypothetical, there is no information regarding the operating condition after a foreseen contingency. The proposed approach first estimates the post-contingency operating point as well as possible PV-PQ transitions based on the current operating point. Then the STI for each contingency can be predicted using the estimated operating condition. Numerical results based on IEEE 14-bus system demonstrate the accuracy of the proposed approach in predicting the voltage stability margin under contingency. Moreover, the on-line implementation of the proposed approach is promising since it only involves solving several linear equations.


title: Categorizing Variants of Goodhart's Law

abstract: There are several distinct failure modes for overoptimization of systems on the basis of metrics. This occurs when a metric which can be used to improve a system is used to an extent that further optimization is ineffective or harmful, and is sometimes termed Goodhart's Law. This class of failure is often poorly understood, partly because terminology for discussing them is ambiguous, and partly because discussion using this ambiguous terminology ignores distinctions between different failure modes of this general type. This paper expands on an earlier discussion by Garrabrant, which notes there are "(at least) four different mechanisms" that relate to Goodhart's Law. This paper is intended to explore these mechanisms further, and specify more clearly how they occur. This discussion should be helpful in better understanding these types of failures in economic regulation, in public policy, in machine learning, and in Artificial Intelligence alignment. The importance of Goodhart effects depends on the amount of power directed towards optimizing the proxy, and so the increased optimization power offered by artificial intelligence makes it especially critical for that field.


title: Some properties of ordinary sense slice 1-links: Some answers to the problem (26) of Fox

abstract: We prove that, for any ordinary sense slice 1-link $L$, we can define the Arf invariant and Arf(L)=0. We prove that, for any m-component 1-link L_1, there exists a 3m-component ordinary sense slaice 1-link L_2 of which L_1 is a sublink.


title: Constrained Optimal Consensus in Dynamical Networks

abstract: In this paper, an optimal consensus problem with local inequality constraints is studied for a network of single-integrator agents. The goal is that a group of single-integrator a gents rendezvous at the optimal point of the sum of local convex objective functions. The local objective functions are only available to the corresponding agents that only need to know their relative distances from their neighbors in order to seek the final optimal point. This point is supposed to be confined by some local inequality constraints. To tackle this problem, we integrate the primal dual gradient-based optimization algorithm with a consensus protocol to drive the agents toward the agreed point that satisfies KKT conditions. The asymptotic convergence of the solution of the optimization problem is proven with the help of LaSalle's invariance principle for hybrid systems. A numerical example is presented to show the effectiveness of our protocol.


title: Electrical detection of a skyrmion in a magnetic tunneling junction

abstract: We theoretically investigated a method to detect a single skyrmion in a magnetic tunneling junction (MTJ) geometry. Using the tunneling Hamiltonian approach, we calculated the tunneling magnetoresistance (TMR) ratio of the skyrmion-ferromagnet bilayer system. We show the TMR ratio is determined sorely by the spin profile of the skyrmion and geometrical factor of the device, if only the system is reasonably clean such that the spectral broadening is smaller than the exchange coupling between the local and the itinerant magnetic moment. The TMR ratio in that case can amount to $30\%$ or higher when the diameter of the skyrmion is as large as the size of the device. Since this criterion is easily achievable in real systems, MTJ geometry can be a good candidate of the electrical detection of a single skyrmion i.e., the reading process of the information in the future skyrmionics memory devices.


title: Height functions for motives, II

abstract: This is the Part II of our paper "Height functions for motives". We consider more general period domains and the height functions on more general sets of motives. We also consider the corresponding Hodge theoretic variant of Nevanlinna theory.


title: Causality and Legendrian linking for higher dimensional spacetimes

abstract: Let $(X^{m+1}, g)$ be an $(m+1)$-dimensional globally hyperbolic spacetime with Cauchy surface $M^m$, and let $\widetilde M^m$ be the universal cover of the Cauchy surface. Let $\mathcal N_{X}$ be the contact manifold of all future directed unparameterized light rays in $X$ that we identify with the spherical cotangent bundle $ST^*M.$ Jointly with Stefan Nemirovski we showed when $\widetilde M^m$ is {\bf not\/} a compact manifold, then two points $x, y\in X$ are causally related if and only if the Legendrian spheres $\mathfrak S_x, \mathfrak S_y$ of all light rays through $x$ and $y$ are linked in $\mathcal N_{X}.$ In this short note we use the contact Bott-Samelson theorem of Frauenfelder, Labrousse and Schlenk to show that the same statement is true for all $X$ for which the integral cohomology ring of a closed $\widetilde M$ is {\bf not} the one of the CROSS (compact rank one symmetric space). If $M$ admits a Riemann metric $\overline g$, a point $x$ and a number $\ell>0$ such that all unit speed geodesics starting from $x$ return back to $x$ in time $\ell$, then $(M, \overline g)$ is called a $Y^x_{\ell}$ manifold. Jointly with Stefan Nemirovski we observed that causality in $(M\times \mathbb R, \overline g\oplus -t^2)$ is {\bf not} equivalent to Legendrian linking. Every $Y^x_{\ell}$-Riemann manifold has compact universal cover and its integral cohomology ring is the one of a CROSS. So we conjecture that Legendrian linking is equivalent to causality if and only if one can {\bf not} put a $Y^x_{\ell}$ Riemann metric on a Cauchy surface $M.$


title: A Generalization of the Robust Positive Expectation Theorem for Stock Trading via Feedback Control

abstract: The starting point of this paper is the so-called Robust Positive Expectation (RPE) Theorem, a result which appears in literature in the context of Simultaneous Long-Short stock trading. This theorem states that using a combination of two specially-constructed linear feedback trading controllers, one long and one short, the expected value of the resulting gain-loss function is guaranteed to be robustly positive with respect to a large class of stochastic processes for the stock price. The main result of this paper is a generalization of this theorem. Whereas previous work applies to a single stock, in this paper, we consider a pair of stocks. To this end, we make two assumptions on their expected returns. The first assumption involves price correlation between the two stocks and the second involves a bounded non-zero momentum condition. With known uncertainty bounds on the parameters associated with these assumptions, our new version of the RPE Theorem provides necessary and sufficient conditions on the positive feedback parameter K of the controller under which robust positive expectation is assured. We also demonstrate that our result generalizes the one existing for the single-stock case. Finally, it is noted that our results also can be interpreted in the context of pairs trading.


title: Effects of "stuffing" on the atomic and electronic structure of the pyrochlore Yb$_2$Ti$_2$O$_7$

abstract: There are reasons to believe that the ground state of the magnetic rare earth pyrochlore Yb$_2$Ti$_2$O$_7$ is on the boundary between competing ground states. We have carried out $\mathrm{\textit{ab initio}}$ density functional calculations to determine the most stable chemical formula as a function of the oxygen chemical potential and the likely location of the oxygen atoms in the unit cell of the "stuffed" system. We find that it is energetically favorable in the "stuffed" crystal (with an Yb replacement on a Ti site) to contain oxygen vacancies which dope the Yb 4$\mathrm{\textit{f}}$ orbitals and qualitatively change the electronic properties of the system. In addition, with the inclusion of the contribution of spin-orbit-coupling (SOC) on top of the GGA+U approach, we investigated the electronic structure and the magnetic moments of the most stable "stuffed" system. In our determined "stuffed" structure the valence bands as compared to those of the pure system are pushed down and a change in hybridization between the O 2$\mathrm{\textit{p}}$ orbitals and the metal ion states is found. Our first-principle findings should form a foundation for effective models describing the low-temperature properties of this material whose true ground state remains controversial.


title: How Smart Are `Water Smart Landscapes'?

abstract: Understanding the effectiveness of alternative approaches to water conservation is crucially important for ensuring the security and reliability of water services for urban residents. We analyze data from one of the longest-running "cash for grass" policies - the Southern Nevada Water Authority's Water Smart Landscapes program, where homeowners are paid to replace grass with xeric landscaping. We use a twelve year long panel dataset of monthly water consumption records for 300,000 households in Las Vegas, Nevada. Utilizing a panel difference-in-differences approach, we estimate the average water savings per square meter of turf removed. We find that participation in this program reduced the average treated household's consumption by 18 percent. We find no evidence that water savings degrade as the landscape ages, or that water savings per unit area are influenced by the value of the rebate. Depending on the assumed time horizon of benefits from turf removal, we find that the WSL program cost the water authority about $1.62 per thousand gallons of water saved, which compares favorably to alternative means of water conservation or supply augmentation.