title: Combining Symbolic Execution and Model Checking to Verify MPI Programs

abstract: Message passing is the standard paradigm of programming in high-performance computing. However, verifying Message Passing Interface (MPI) programs is challenging, due to the complex program features (such as non-determinism and non-blocking operations). In this work, we present MPI symbolic verifier (MPI-SV), the first symbolic execution based tool for automatically verifying MPI programs with non-blocking operations. MPI-SV combines symbolic execution and model checking in a synergistic way to tackle the challenges in MPI program verification. The synergy improves the scalability and enlarges the scope of verifiable properties. We have implemented MPI-SV (footnote: https://mpi-sv.github.io) and evaluated it with 111 real-world MPI verification tasks. The pure symbolic execution-based technique successfully verifies 61 out of the 111 tasks (55\%) within one hour, while in comparison, MPI-SV verifies 100 tasks (90\%). On average, compared with pure symbolic execution, MPI-SV achieves 19x speedups on verifying the satisfaction of the critical property and 5x speedups on finding violations.


title: Improved Part Segmentation Performance by Optimising Realism of Synthetic Images using Cycle Generative Adversarial Networks

abstract: In this paper we report on improved part segmentation performance using convolutional neural networks to reduce the dependency on the large amount of manually annotated empirical images. This was achieved by optimising the visual realism of synthetic agricultural images.In Part I, a cycle consistent generative adversarial network was applied to synthetic and empirical images with the objective to generate more realistic synthetic images by translating them to the empirical domain. We first hypothesise and confirm that plant part image features such as color and texture become more similar to the empirical domain after translation of the synthetic images.Results confirm this with an improved mean color distribution correlation with the empirical data prior of 0.62 and post translation of 0.90. Furthermore, the mean image features of contrast, homogeneity, energy and entropy moved closer to the empirical mean, post translation. In Part II, 7 experiments were performed using convolutional neural networks with different combinations of synthetic, synthetic translated to empirical and empirical images. We hypothesised that the translated images can be used for (i) improved learning of empirical images, and (ii) that learning without any fine-tuning with empirical images is improved by bootstrapping with translated images over bootstrapping with synthetic images. Results confirm our second and third hypotheses. First a maximum intersection-over-union performance was achieved of 0.52 when bootstrapping with translated images and fine-tuning with empirical images; an 8% increase compared to only using synthetic images. Second, training without any empirical fine-tuning resulted in an average IOU of 0.31; a 55% performance increase over previous methods that only used synthetic images.


title: Electron-mediated phonon-phonon coupling drives the vibrational relaxation of CO on Cu(100)

abstract: We bring forth a consistent theory for the electron-mediated vibrational intermode coupling that clarifies the microscopic mechanism behind the vibrational relaxation of adsorbates on metal surfaces. Our analysis points out the inability of state-of-the-art nonadiabatic theories to quantitatively reproduce the experimental linewidth of the CO internal stretch mode on Cu(100) and it emphasizes the crucial role of the electron-mediated phonon-phonon coupling in this regard. The results demonstrate a strong electron-mediated coupling between the internal stretch and low-energy CO modes, but also a significant role of surface motion. Our nonadiabatic theory is also able to explain the temperature dependence of the internal stretch phonon linewidth, thus far considered a sign of the direct anharmonic coupling.


title: Optimal control of light propagation and exciton transfer in arrays of molecular-like noble-metal clusters

abstract: We demonstrate theoretically the possibility of optimal control of light propagation and exciton transfer in arrays constructed of subnanometer sized noble-metal clusters by using phase-shaped laser pulses and analyze the mechanism underlying this process. The theoretical approach for simulation of light propagation in the arrays is based on the numerical solution of the coupled time-dependent Schr\"odinger equation and the classical electric field propagation in an iterative self-consistent manner. The electronic eigenstates of individual clusters and the dipole couplings are obtained from ab initio TDDFT calculations. The total electric field is propagated along the array by coupling an external excitation electric field with the electric fields produced by all clusters. A genetic algorithm is used to determine optimal pulse shapes which drive the excitation in a desired direction. The described theoretical approach is applied to control the light propagation and exciton transfer dynamics into a T-shaped structure built of seven Ag8 clusters. We demonstrate that a selective switching of light localization is possible in $\sim$5 nm sized cluster arrays which might serve as a building block for plasmonic devices with an ultrafast operation regime


title: "Cultural additivity" and how the values and norms of Confucianism, Buddhism, and Taoism co-exist, interact, and influence Vietnamese society: A Bayesian analysis of long-standing folktales, using R and Stan

abstract: Every year, the Vietnamese people reportedly burned about 50,000 tons of joss papers, which took the form of not only bank notes, but iPhones, cars, clothes, even housekeepers, in hope of pleasing the dead. The practice was mistakenly attributed to traditional Buddhist teachings but originated in fact from China, which most Vietnamese were not aware of. In other aspects of life, there were many similar examples of Vietnamese so ready and comfortable with adding new norms, values, and beliefs, even contradictory ones, to their culture. This phenomenon, dubbed "cultural additivity", prompted us to study the co-existence, interaction, and influences among core values and norms of the Three Teachings--Confucianism, Buddhism, and Taoism--as shown through Vietnamese folktales. By applying Bayesian logistic regression, we evaluated the possibility of whether the key message of a story was dominated by a religion (dependent variables), as affected by the appearance of values and anti-values pertaining to the Three Teachings in the story (independent variables).


title: C-LSTM: Enabling Efficient LSTM using Structured Compression Techniques on FPGAs

abstract: Recently, significant accuracy improvement has been achieved for acoustic recognition systems by increasing the model size of Long Short-Term Memory (LSTM) networks. Unfortunately, the ever-increasing size of LSTM model leads to inefficient designs on FPGAs due to the limited on-chip resources. The previous work proposes to use a pruning based compression technique to reduce the model size and thus speedups the inference on FPGAs. However, the random nature of the pruning technique transforms the dense matrices of the model to highly unstructured sparse ones, which leads to unbalanced computation and irregular memory accesses and thus hurts the overall performance and energy efficiency. In contrast, we propose to use a structured compression technique which could not only reduce the LSTM model size but also eliminate the irregularities of computation and memory accesses. This approach employs block-circulant instead of sparse matrices to compress weight matrices and reduces the storage requirement from $\mathcal{O}(k^2)$ to $\mathcal{O}(k)$. Fast Fourier Transform algorithm is utilized to further accelerate the inference by reducing the computational complexity from $\mathcal{O}(k^2)$ to $\mathcal{O}(k\text{log}k)$. The datapath and activation functions are quantized as 16-bit to improve the resource utilization. More importantly, we propose a comprehensive framework called C-LSTM to automatically optimize and implement a wide range of LSTM variants on FPGAs. According to the experimental results, C-LSTM achieves up to 18.8X and 33.5X gains for performance and energy efficiency compared with the state-of-the-art LSTM implementation under the same experimental setup, and the accuracy degradation is very small.


title: Block minimal bases $\ell$-ifications of matrix polynomials

abstract: The standard way of solving a polynomial eigenvalue problem associated with a matrix polynomial starts by embedding the matrix coefficients of the polynomial into a matrix pencil, known as a strong linearization. This process transforms the problem into an equivalent generalized eigenvalue problem. However, there are some situations in which is more convenient to replace linearizations by other low degree matrix polynomials. This has motivated the idea of a strong $\ell$-ification of a matrix polynomial, which is a matrix polynomial of degree $\ell$ having the same finite and infinite elementary divisors, and the same numbers of left and right minimal indices as the original matrix polynomial. We present in this work a novel method for constructing strong $\ell$-ifications of matrix polynomials of size $m\times n$ and grade $d$ when $\ell< d$, and $\ell$ divides $nd$ or $md$. This method is based on a family called "strong block minimal bases matrix polynomials", and relies heavily on properties of dual minimal bases. We show how strong block minimal bases $\ell$-ifications can be constructed from the coefficients of a given matrix polynomial $P(\lambda)$. We also show that these $\ell$-ifications satisfy many desirable properties for numerical applications: they are strong $\ell$-ifications regardless of whether $P(\lambda)$ is regular or singular, the minimal indices of the $\ell$-ifications are related to those of $P(\lambda)$ via constant uniform shifts, and eigenvectors and minimal bases of $P(\lambda)$ can be recovered from those of any of the strong block minimal bases $\ell$-ifications. In the special case where $\ell$ divides $d$, we introduce a subfamily of strong block minimal bases matrix polynomials named "block Kronecker matrix polynomials", which is shown to be a fruitful source of companion $\ell$-ifications.


title: Exchange-correlation functionals of i-DFT for asymmetrically coupled leads

abstract: A recently proposed density functional approach for steady-state transport through nanoscale systems (called i-DFT) is used to investigate junctions which are asymmetrically coupled to the leads and biased with asymmetric voltage drops. In the latter case, the system can simply be transformed to a physically equivalent one with symmetric voltage drop by a total energy shift of the entire system. For the former case, known exchange correlation gate and bias functionals have to be generalized to take into account the asymmetric coupling to the leads. We show how differential conductance spectra of the constant interaction model evolve with increasing asymmetry of both voltage drops and coupling to the leads.


title: Actions of nonamenable groups on $\mathcal{Z}$-stable $C^*$-algebras

abstract: We study strongly outer actions of discrete groups on C*-algebras in relation to (non)amenability. In contrast to related results for amenable groups, where uniqueness of strongly outer actions on the Jiang-Su algebra is expected, we show that uniqueness fails for all nonamenable groups, and that the failure is drastic. Our main result implies that if $G$ contains a copy of the free group, then there exist uncountable many, non-cocycle conjugate strongly outer actions of $G$ on any Jiang-Su stable tracial C*-algebra. Similar conclusions apply for outer actions on McDuff tracial von Neumann algebras. We moreover show that $G$ is amenable if and only if the Bernoulli shift on a finite strongly self-absorbing C*-algebra absorbs the trivial action on the Jiang-Su algebra. Our methods consist in a careful study of weak containments of the Koopman representations of different Bernoulli-type actions.


title: Modified dipole-dipole interaction and dissipation in an atomic ensemble near surfaces

abstract: We study how the radiative properties of a dense ensemble of atoms can be modified when they are placed near or between metallic or dielectric surfaces. If the average separation between the atoms is comparable or smaller than the wavelength of the scattered photons, the coupling to the radiation field induces long-range coherent interactions based on the interatomic exchange of virtual photons. Moreover, the incoherent scattering of photons back to the electromagnetic field is known to be a many-body process, characterized by the appearance of superradiant and subradiant emission modes. By changing the radiation field properties, in this case by considering a layered medium where the atoms are near metallic or dielectric surfaces, these scattering properties can be dramatically modified. We perform a detailed study of these effects, with focus on experimentally relevant parameter regimes. We finish with a specific application in the context of quantum information storage, where the presence of a nearby surface is shown to increase the storage time of an atomic excitation that is transported across a one-dimensional chain.


title: Yangian Symmetry for the Action of Planar N=4 Super Yang-Mills and N=6 Super Chern-Simons Theories

abstract: In this article we establish the notion of classical Yangian symmetry for planar N=4 supersymmetric Yang-Mills theory and for related planar gauge theories. After revisiting Yangian invariance for the equations of motion, we describe how the bi-local generators act on the action of the model such that the latter becomes exactly invariant. In particular, we elaborate on the relevance of the planar limit and how to act non-linearly with bi-local generators on the cyclic action.


title: Experimental Demonstration of Non-Destructive Discrimination of Arbitrary Set of Orthogonal Quantum States Using 5-qubit IBM Quantum Computer on Cloud

abstract: A protocol for non-destructive descrimination of arbitrary set of orthogonal quantum states was proposed by V. S. Manu et al., using an algorithm based on quantum phase estimation. IBM Corporation has released a superconductivity based 5-qubit (5-qubit transmon bowtie chip 3 and IBM 5-qubit real processor) quantum computer named Quantum Experience and placed it on cloud. In this paper we take advantage of the online availability of those real quantum processors(ibmqx2 amd ibmqx4) and carry out the above protocol that has experimentally demonstrated earlier using NMR quantum processor. Here, we set up experiments for arbitrary one-qubit and two-qubit orthogonal quantum states. The experiment confirmed that the arbitrary orthogonal quantum states can be discriminate in a nondestructive manner with a high fidelity. We compare the outcomes of those experiments which are done by ibmqx2 and ibmqx4 processors. Here, we also show the state tomography for the single qubit experiments.


title: EVA$^2$: Exploiting Temporal Redundancy in Live Computer Vision

abstract: Hardware support for deep convolutional neural networks (CNNs) is critical to advanced computer vision in mobile and embedded devices. Current designs, however, accelerate generic CNNs; they do not exploit the unique characteristics of real-time vision. We propose to use the temporal redundancy in natural video to avoid unnecessary computation on most frames. A new algorithm, activation motion compensation, detects changes in the visual input and incrementally updates a previously-computed output. The technique takes inspiration from video compression and applies well-known motion estimation techniques to adapt to visual changes. We use an adaptive key frame rate to control the trade-off between efficiency and vision quality as the input changes. We implement the technique in hardware as an extension to existing state-of-the-art CNN accelerator designs. The new unit reduces the average energy per frame by 54.2%, 61.7%, and 87.6% for three CNNs with less than 1% loss in vision accuracy.


title: Error Analysis of an Incremental POD Algorithm for PDE Simulation Data

abstract: In our earlier work [Fareed et al., Comput. Math. Appl. 75 (2018), no. 6, 1942-1960], we proposed an incremental SVD algorithm with respect to a weighted inner product to compute the proper orthogonal decomposition (POD) of a set of simulation data for a partial differential equation (PDE) without storing the data. In this work, we perform an error analysis of the incremental SVD algorithm. We also modify the algorithm to incrementally update both the SVD and an error bound when a new column of data is added. We show the algorithm produces the exact SVD of an approximate data matrix, and the operator norm error between the approximate and exact data matrices is bounded above by the computed error bound. This error bound also allows us to bound the error in the incrementally computed singular values and singular vectors. We illustrate our analysis with numerical results for three simulation data sets from a 1D FitzHugh-Nagumo PDE system with various choices of the algorithm truncation tolerances.


title: Potential Impact of Global Navigation Satellite Services on Total Power HI Intensity Mapping Surveys

abstract: Future total-power single-dish HI intensity mapping (HI IM) surveys have the potential to provide unprecedented insight into late time ($z < 1$) cosmology that are competitive with Stage IV dark energy surveys. However, redshifts between $0 < z < 0.2$ lie within the transmission bands of global navigation satellite services (GNSS), and even at higher redshifts out-of-band leakage from GNSS satellites may be problematic. We estimate the impact of GNSS satellites on future single-dish HI IM surveys using realistic estimates of both the total power and spectral structure of GNSS signals convolved with a model SKA beam. Using a simulated SKA HI IM survey covering 30000 sq. deg. of sky and 200 dishes, we compare the integrated GNSS emission on the sky with the expected HI signal. It is found that for frequencies $> 950$ MHz the emission from GNSS satellites will exceed the expected HI signal for all angular scales to which the SKA is sensitive when operating in single-dish mode.


title: Benchmarks for cyber-physical systems: A modular model library for building automation systems (Extended version)

abstract: Building Automation Systems (BAS) are exemplars of Cyber-Physical Systems (CPS), incorporating digital control architectures over underlying continuous physical processes. We provide a modular model library for BAS drawn from expertise developed on a real BAS setup. The library allows to build models comprising of either physical quantities or digital control modules.% which are composable. The structure, operation, and dynamics of the model can be complex, incorporating (i) stochasticity, (ii) non-linearities, (iii) numerous continuous variables or discrete states, (iv) various input and output signals, and (v) a large number of possible discrete configurations. The modular composition of BAS components can generate useful CPS benchmarks. We display this use by means of three realistic case studies, where corresponding models are built and engaged with different analysis goals. The benchmarks, the model library and data collected from the BAS setup at the University of Oxford, are kept on-line at https://github.com/natchi92/BASBenchmarks.


title: Temporal Gaussian Mixture Layer for Videos

abstract: We introduce a new convolutional layer named the Temporal Gaussian Mixture (TGM) layer and present how it can be used to efficiently capture longer-term temporal information in continuous activity videos. The TGM layer is a temporal convolutional layer governed by a much smaller set of parameters (e.g., location/variance of Gaussians) that are fully differentiable. We present our fully convolutional video models with multiple TGM layers for activity detection. The extensive experiments on multiple datasets, including Charades and MultiTHUMOS, confirm the effectiveness of TGM layers, significantly outperforming the state-of-the-arts.


title: A local characterization of crystals for the quantum queer superalgebra

abstract: We define operators on semistandard shifted tableaux and use Stembridge's local characterization for regular graphs to prove they define a crystal structure. This gives a new proof that Schur $P$-polynomials are Schur positive. We define queer crystal operators (also called odd Kashiwara operators) to construct a connected queer crystal on semistandard shifted tableaux of a given shape. Using the tensor rule for queer crystals, this provides a new proof that products of Schur $P$-polynomials are Schur $P$-positive. Finally, to facilitate applications of queer crystals in the context of Schur $P$-positivity, we give local axioms for queer regular graphs, generalizing Stembridge's axioms, that partially characterize queer crystals.


title: Holographic meson decays via worldsheet instantons

abstract: We study meson decays using instanton methods in two string models. The first model is the old string model in flat space which combines strings and massive particles and the second is the holographic, Sakai-Sugimoto model. Using the the old string model, we reproduce the QCD formula for the probability of splitting of the QCD flux tube derived by Casher-Neuberger-Nussinov (CNN). In the holographic model we construct a string worldsheet instanton which interpolates between a single and double string configuration, which determines the decay from one to two dual mesonic particles. The resulting probability for meson decay incorporates both the effects of finite meson size as well as back-reaction of the produced quarks on the QCD flux tube. In the limit of very large strings the probability for a split reduces to the CNN formula. A byproduct of our analysis is the analysis of the moduli space of a generic double concentric Wilson loop with circles which are separated in the holographic direction of the confining background.


title: Overcoming correlation fluctuations in two-photon interference experiments with differently bright and independently blinking remote quantum emitters

abstract: As a fundamental building block for quantum computation and communication protocols, the correct verification of the two-photon interference (TPI) contrast between two independent quantum light sources is of utmost importance. Here, we experimentally demonstrate how frequently present blinking dynamics and changes in emitter brightness critically affect the Hong-Ou-Mandel-type (HOM) correlation histograms of remote TPI experiments measured via the commonly utilized setup configuration. We further exploit this qualitative and quantitative explanation of the observed correlation dynamics to establish an alternative interferometer configuration, which is overcoming the discussed temporal fluctuations, giving rise to an error-free determination of the remote TPI visibility. We prove full knowledge of the obtained correlation by reproducing the measured correlation statistics via Monte-Carlo simulations. As exemplary system, we make use of two pairs of remote semiconductor quantum dots, however, the same conclusions apply for TPI experiments with flying qubits from any kind of remote solid state quantum emitters.


title: Synchronisation of Partial Multi-Matchings via Non-negative Factorisations

abstract: In this work we study permutation synchronisation for the challenging case of partial permutations, which plays an important role for the problem of matching multiple objects (e.g. images or shapes). The term synchronisation refers to the property that the set of pairwise matchings is cycle-consistent, i.e. in the full matching case all compositions of pairwise matchings over cycles must be equal to the identity. Motivated by clustering and matrix factorisation perspectives of cycle-consistency, we derive an algorithm to tackle the permutation synchronisation problem based on non-negative factorisations. In order to deal with the inherent non-convexity of the permutation synchronisation problem, we use an initialisation procedure based on a novel rotation scheme applied to the solution of the spectral relaxation. Moreover, this rotation scheme facilitates a convenient Euclidean projection to obtain a binary solution after solving our relaxed problem. In contrast to state-of-the-art methods, our approach is guaranteed to produce cycle-consistent results. We experimentally demonstrate the efficacy of our method and show that it achieves better results compared to existing methods.


title: A particle-based variational approach to Bayesian Non-negative Matrix Factorization

abstract: Bayesian Non-negative Matrix Factorization (NMF) is a promising approach for understanding uncertainty and structure in matrix data. However, a large volume of applied work optimizes traditional non-Bayesian NMF objectives that fail to provide a principled understanding of the non-identifiability inherent in NMF-- an issue ideally addressed by a Bayesian approach. Despite their suitability, current Bayesian NMF approaches have failed to gain popularity in an applied setting; they sacrifice flexibility in modeling for tractable computation, tend to get stuck in local modes, and require many thousands of samples for meaningful uncertainty estimates. We address these issues through a particle-based variational approach to Bayesian NMF that only requires the joint likelihood to be differentiable for tractability, uses a novel initialization technique to identify multiple modes in the posterior, and allows domain experts to inspect a `small' set of factorizations that faithfully represent the posterior. We introduce and employ a class of likelihood and prior distributions for NMF that formulate a Bayesian model using popular non-Bayesian NMF objectives. On several real datasets, we obtain better particle approximations to the Bayesian NMF posterior in less time than baselines and demonstrate the significant role that multimodality plays in NMF-related tasks.


title: Computing performability measures in Markov chains by means of matrix functions

abstract: We discuss the efficient computation of performance, reliability, and availability measures for Markov chains; these metrics, and the ones obtained by combining them, are often called performability measures. We show that this computational problem can be recasted as the evaluation of a bilinear forms induced by appropriate matrix functions, and thus solved by leveraging the fast methods available for this task. We provide a comprehensive analysis of the theory required to translate the problem from the language of Markov chains to the one of matrix functions. The advantages of this new formulation are discussed, and it is shown that this setting allows to easily study the sensitivities of the measures with respect to the model parameters. Numerical experiments confirm the effectiveness of our approach; the tests we have run show that we can outperform the solvers available in state of the art commercial packages on a representative set of large scale examples.


title: Minimal space with non-minimal square

abstract: We completely solve the problem whether the product of two compact metric spaces admitting minimal maps also admits a minimal map. Recently Boro\'nski, Clark and Oprocha gave a negative answer in the particular case when homeomorphisms rather than continuous maps are considered. In the present paper we show that there is a metric continuum $X$ admitting a minimal map, in fact a minimal homeomorphism, such that $X\times X$ does not admit any minimal map.


title: Fast approximation and exact computation of negative curvature parameters of graphs

abstract: In this paper, we study Gromov hyperbolicity and related parameters, that represent how close (locally) a metric space is to a tree from a metric point of view. The study of Gromov hyperbolicity for geodesic metric spaces can be reduced to the study of graph hyperbolicity. The main contribution of this paper is a new characterization of the hyperbolicity of graphs. This characterization has algorithmic implications in the field of large-scale network analysis. A sharp estimate of graph hyperbolicity is useful, e.g., in embedding an undirected graph into hyperbolic space with minimum distortion [Verbeek and Suri, SoCG'14]. The hyperbolicity of a graph can be computed in polynomial-time, however it is unlikely that it can be done in subcubic time. This makes this parameter difficult to compute or to approximate on large graphs. Using our new characterization of graph hyperbolicity, we provide a simple factor 8 approximation algorithm for computing the hyperbolicity of an $n$-vertex graph $G=(V,E)$ in optimal time $O(n^2)$ (assuming that the input is the distance matrix of the graph). This algorithm leads to constant factor approximations of other graph-parameters related to hyperbolicity (thinness, slimness, and insize). We also present the first efficient algorithms for exact computation of these parameters. All of our algorithms can be used to approximate the hyperbolicity of a geodesic metric space. We also show that a similar characterization of hyperbolicity holds for all geodesic metric spaces endowed with a geodesic spanning tree. Along the way, we prove that any complete geodesic metric space $(X,d)$ has such a geodesic spanning tree. We hope that this fundamental result can be useful in other contexts.


title: Lie algebras and $v_n$-periodic spaces

abstract: We consider a homotopy theory obtained from that of pointed spaces by inverting the maps inducing isomorphisms in $v_n$-periodic homotopy groups. The case n = 0 corresponds to rational homotopy theory. In analogy with Quillen's results in the rational case, we prove that this $v_n$-periodic homotopy theory is equivalent to the homotopy theory of Lie algebras in T(n)-local spectra. We also compare it to the homotopy theory of commutative coalgebras in T(n)-local spectra, where it turns out there is only an equivalence up to a certain convergence issue of the Goodwillie tower of the identity.


title: Unified Phase Diagram of Antiferromagnetic SU(N) Spin Ladders

abstract: Motivated by near-term experiments with ultracold alkaline-earth atoms confined to optical lattices, we establish numerically and analytically the phase diagram of two-leg SU($N$) spin ladders. Two-leg ladders provide a rich and highly non-trivial extension of the single chain case on the way towards the relatively little explored two dimensional situation. Focusing on the experimentally relevant limit of one fermion per site, antiferromagnetic exchange interactions, and $2\leq N \leq 6$, we show that the phase diagrams as a function of the interchain (rung) to intrachain (leg) coupling ratio $J_\perp/J_\Vert$ strongly differ for even vs. odd $N$. For even $N=4$ and 6, we demonstrate that the phase diagram consists of a single valence bond crystal (VBC) with a spatial period of $N/2$ rungs. For odd $N=3$ and 5, we find surprisingly rich phase diagrams exhibiting three distinct phases. For weak rung coupling, we obtain a VBC with a spatial period of $N$ rungs, whereas for strong coupling we obtain a critical phase related to the case of a single chain. In addition, we encounter intermediate phases for odd $N$, albeit of a different nature for $N=3$ as compared to $N=5$. For $N=3$, we find a novel gapless intermediate phase with $J_\perp$-dependent incommensurate spatial fluctuations in a sizeable region of the phase diagram. For $N=5$, there are strong indications for a narrow potentially gapped intermediate phase, whose nature is not entirely clear.


title: A Covariant Approach to Entropic Dynamics

abstract: Entropic Dynamics (ED) is a framework for constructing dynamical theories of inference using the tools of inductive reasoning. A central feature of the ED framework is the special focus placed on time. In previous work a global entropic time was used to derive a quantum theory of relativistic scalar fields. This theory, however, suffered from a lack of explicit or manifest Lorentz symmetry. In this paper we explore an alternative formulation in which the relativistic aspects of the theory are manifest. The approach we pursue here is inspired by the works of Dirac, Kuchar, and Teitelboim in their development of covariant Hamiltonian methods. The key ingredient here is the adoption of a local notion of time, which we call entropic time. This construction allows the expression of arbitrary notion of simultaneity, in accord with relativity. In order to ensure, however, that this local time dynamics is compatible with the background spacetime we must impose a set of Poisson bracket constraints; these constraints themselves result from requiring the dynamcics to be path independent, in the sense of Teitelboim and Kuchar.


title: Nesting Probabilistic Programs

abstract: We formalize the notion of nesting probabilistic programming queries and investigate the resulting statistical implications. We demonstrate that while query nesting allows the definition of models which could not otherwise be expressed, such as those involving agents reasoning about other agents, existing systems take approaches which lead to inconsistent estimates. We show how to correct this by delineating possible ways one might want to nest queries and asserting the respective conditions required for convergence. We further introduce a new online nested Monte Carlo estimator that makes it substantially easier to ensure these conditions are met, thereby providing a simple framework for designing statistically correct inference engines. We prove the correctness of this online estimator and show that, when using the recommended setup, its asymptotic variance is always better than that of the equivalent fixed estimator, while its bias is always within a factor of two.


title: Learning deep structured active contours end-to-end

abstract: The world is covered with millions of buildings, and precisely knowing each instance's position and extents is vital to a multitude of applications. Recently, automated building footprint segmentation models have shown superior detection accuracy thanks to the usage of Convolutional Neural Networks (CNN). However, even the latest evolutions struggle to precisely delineating borders, which often leads to geometric distortions and inadvertent fusion of adjacent building instances. We propose to overcome this issue by exploiting the distinct geometric properties of buildings. To this end, we present Deep Structured Active Contours (DSAC), a novel framework that integrates priors and constraints into the segmentation process, such as continuous boundaries, smooth edges, and sharp corners. To do so, DSAC employs Active Contour Models (ACM), a family of constraint- and prior-based polygonal models. We learn ACM parameterizations per instance using a CNN, and show how to incorporate all components in a structured output model, making DSAC trainable end-to-end. We evaluate DSAC on three challenging building instance segmentation datasets, where it compares favorably against state-of-the-art. Code will be made available.


title: Interior metric and ray-tracing map in the firework black-to-white hole transition

abstract: The possibility that a black hole could tunnel into a white hole has recently received attention. Here we present a metric that improves the 'firework' metric: it describes the entire process and solves the Einstein equations everywhere except on a small transition surface that corresponds to the quantum tunnelling. We compute the corresponding ray-tracing map from past infinity to future infinity explicitly.


title: Measurement of the Dynamical Structure Factor of a 1D Interacting Fermi Gas

abstract: We present measurements of the dynamical structure factor $S(q,\omega)$ of an interacting one-dimensional (1D) Fermi gas for small excitation energies. We use the two lowest hyperfine levels of the $^6$Li atom to form a pseudo-spin-1/2 system whose s-wave interactions are tunable via a Feshbach resonance. The atoms are confined to 1D by a two-dimensional optical lattice. Bragg spectroscopy is used to measure a response of the gas to density ("charge") mode excitations at a momentum $q$ and frequency $\omega$. The spectrum is obtained by varying $\omega$, while the angle between two laser beams determines $q$, which is fixed to be less than the Fermi momentum $k_\textrm{F}$. The measurements agree well with Tomonaga-Luttinger theory.


title: On the MSW neutrino mixing effects in atomic weak interactions and double beta decays

abstract: Matter effects on the mixing of the neutrinos mass eigenstates, also know as the Mikheyev-Smirnov-Wolfenstein effect, seem to be well established in describing the propagation of the neutrino from the source to detecting devices. These effects were mostly considered in bulk matter, but not inside the atoms. Here we consider the effect of the high electron densities existing in the atomic nuclei. We investigate if these effects can affect the known neutrino phenomenology. It was reported that the mixing of the neutrino in high density matter, such as inside a supernova, can affect the Majoron decay probabilities. We investigate if the neutrino mixing effects in the high electron density inside the atomic nuclei can change the neutrinoless double beta decay half-life formula. In both cases we found that the standard results stand. The results look simple, but the road to them is complex and it opens the possibility that the neutrino mixing in atomic nuclei may affect other observables, such as the neutrinoless double beta Majoron decays.


title: Snap ML: A Hierarchical Framework for Machine Learning

abstract: We describe a new software framework for fast training of generalized linear models. The framework, named Snap Machine Learning (Snap ML), combines recent advances in machine learning systems and algorithms in a nested manner to reflect the hierarchical architecture of modern computing systems. We prove theoretically that such a hierarchical system can accelerate training in distributed environments where intra-node communication is cheaper than inter-node communication. Additionally, we provide a review of the implementation of Snap ML in terms of GPU acceleration, pipelining, communication patterns and software architecture, highlighting aspects that were critical for achieving high performance. We evaluate the performance of Snap ML in both single-node and multi-node environments, quantifying the benefit of the hierarchical scheme and the data streaming functionality, and comparing with other widely-used machine learning software frameworks. Finally, we present a logistic regression benchmark on the Criteo Terabyte Click Logs dataset and show that Snap ML achieves the same test loss an order of magnitude faster than any of the previously reported results, including those obtained using TensorFlow and scikit-learn.


title: Forbidden arithmetic progressions in permutations of subsets of the integers

abstract: Permutations of the positive integers avoiding arithmetic progressions of length $5$ were constructed in (Davis et al, 1977), implying the existence of permutations of the integers avoiding arithmetic progressions of length $7$. We construct a permutation of the integers avoiding arithmetic progressions of length $6$. We also prove a lower bound of $\frac{1}{2}$ on the lower density of subsets of positive integers that can be permuted to avoid arithmetic progressions of length $4$, sharpening the lower bound of $\frac{1}{3}$ from (LeSaulnier and Vijay, 2011). In addition, we generalize several results about forbidden arithmetic progressions to construct permutations avoiding generalized arithmetic progressions.


title: Black holes in quartic-order beyond-generalized Proca theories

abstract: The generalized Proca theories with second-order equations of motion can be healthily extended to a more general framework in which the number of propagating degrees of freedom remains unchanged. In the presence of a quartic-order nonminimal coupling to gravity arising in beyond-generalized Proca theories, the speed of gravitational waves $c_t$ on the Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological background can be equal to that of light $c$ under a certain condition. By using this condition alone, we show that the speed of gravitational waves in the vicinity of static and spherically symmetric black holes is also equivalent to $c$ for the propagation of odd-parity perturbations along both radial and angular directions. As a by-product, the black holes arising in our beyond-generalized Proca theories are plagued by neither ghost nor Laplacian instabilities against odd-parity perturbations. We show the existence of both exact and numerical black hole solutions endowed with vector hairs induced by the quartic-order coupling.


title: Applying the Delta method in metric analytics: A practical guide with novel ideas

abstract: During the last decade, the information technology industry has adopted a data-driven culture, relying on online metrics to measure and monitor business performance. Under the setting of big data, the majority of such metrics approximately follow normal distributions, opening up potential opportunities to model them directly without extra model assumptions and solve big data problems via closed-form formulas using distributed algorithms at a fraction of the cost of simulation-based procedures like bootstrap. However, certain attributes of the metrics, such as their corresponding data generating processes and aggregation levels, pose numerous challenges for constructing trustworthy estimation and inference procedures. Motivated by four real-life examples in metric development and analytics for large-scale A/B testing, we provide a practical guide to applying the Delta method, one of the most important tools from the classic statistics literature, to address the aforementioned challenges. We emphasize the central role of the Delta method in metric analytics by highlighting both its classic and novel applications.


title: Electro-rheological fluids under random influences: martingale and strong solutions

abstract: We study generalised Navier--Stokes equations governing the motion of an electro-rheological fluid subject to stochastic perturbation. Stochastic effects are implemented through (i) random initial data, (ii) a forcing term in the momentum equation represented by a multiplicative white noise and (iii) a random character of the variable exponent $p=p(\omega,t,x)$ (as a result of a random electric field). We show the existence of a weak martingale solution provided the variable exponent satisfies $p\geq p^->\frac{3n}{n+2}$ ($p^->1$ in two dimensions). Under additional assumptions we obtain also pathwise solutions.


title: Elliptic flow in a strongly-interacting normal Bose gas

abstract: We study the anisotropic, elliptic expansion of a thermal atomic Bose gas released from an anisotropic trapping potential, for a wide range of interaction strengths across a Feshbach resonance. We show that in our system this hydrodynamic phenomenon is for all interaction strengths fully described by a microscopic kinetic model with no free parameters. The success of this description crucially relies on taking into account the reduced thermalising power of elastic collisions in a strongly interacting gas, for which we derive an analytical theory. We also perform time-resolved measurements that directly reveal the dynamics of the energy transfer between the different expansion axes.

categories: cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph physics.plasm-ph quant-ph


title: A time dependent Stokes interface problem: well-posedness and space-time finite element discretization

abstract: In this paper a time dependent Stokes problem that is motivated by a standard sharp interface model for the fluid dynamics of two-phase flows is studied. This Stokes interface problem has discontinuous density and viscosity coefficients and a pressure solution that is discontinuous across an evolving interface. This strongly simplified two-phase Stokes equation is considered to be a good model problem for the development and analysis of finite element discretization methods for two-phase flow problems. In view of the unfitted finite element methods that are often used for two-phase flow simulations, we are particularly interested in a well-posed variational formulation of this Stokes interface problem in a Euclidean setting. Such well-posed weak formulations, which are not known in the literature, are the main results of this paper. Different variants are considered, namely one with suitable spaces of divergence free functions, a discrete-in-time version of it, and variants in which the divergence free constraint in the solution space is treated by a pressure Lagrange multiplier. The discrete-in-time variational formulation involving the pressure variable for the divergence free constraint is a natural starting point for a space-time finite element discretization. Such a method is introduced and results of numerical experiments with this method are presented.


title: Faces as Lighting Probes via Unsupervised Deep Highlight Extraction

abstract: We present a method for estimating detailed scene illumination using human faces in a single image. In contrast to previous works that estimate lighting in terms of low-order basis functions or distant point lights, our technique estimates illumination at a higher precision in the form of a non-parametric environment map. Based on the observation that faces can exhibit strong highlight reflections from a broad range of lighting directions, we propose a deep neural network for extracting highlights from faces, and then trace these reflections back to the scene to acquire the environment map. Since real training data for highlight extraction is very limited, we introduce an unsupervised scheme for finetuning the network on real images, based on the consistent diffuse chromaticity of a given face seen in multiple real images. In tracing the estimated highlights to the environment, we reduce the blurring effect of skin reflectance on reflected light through a deconvolution determined by prior knowledge on face material properties. Comparisons to previous techniques for highlight extraction and illumination estimation show the state-of-the-art performance of this approach on a variety of indoor and outdoor scenes.


title: Distributed Transactions: Dissecting the Nightmare

abstract: Many distributed storage systems are transactional and a lot of work has been devoted to optimizing their performance, especially the performance of read-only transactions that are considered the most frequent in practice. Yet, the results obtained so far are rather disappointing, and some of the design decisions seem contrived. This paper contributes to explaining this state of affairs by proving intrinsic limitations of transactional storage systems, even those that need not ensure strong consistency but only causality. We first consider general storage systems where some transactions are read-only and some also involve write operations. We show that even read-only transactions cannot be "fast": their operations cannot be executed within one round-trip message exchange between a client seeking an object and the server storing it. We then consider systems (as sometimes implemented today) where all transactions are read-only, i.e., updates are performed as individual operations outside transactions. In this case, read-only transactions can indeed be "fast", but we prove that they need to be "visible". They induce inherent updates on the servers, which in turn impact their overall performance.


title: 3D modelling of macroscopic force-free effects in superconducting thin films and rectangular prisms

abstract: When the magnetic field has a parallel component to the current density ${\bf J}$ there appear force-free effects due to flux cutting and crossing. This results in an anisotropic ${\bf E}({\bf J})$ relation, being ${\bf E}$ the electric field. Understanding force-free effects is interesting not only for the design of superconducting power and magnet applications but also for material characterization. This work develops and applies a fast and accurate computer modeling method based on a variational approach that can handle force-free anisotropic ${\bf E}({\bf J})$ relations and perform fully three dimensional (3D) calculations. We present a systematic study of force-free effects in rectangular thin films and prisms with several finite thicknesses under applied magnetic fields with arbitrary angle $\theta$ with the surface. The results are compared with the same situation with isotropic ${\bf E}({\bf J})$ relation. The thin film situation shows gradual critical current density penetration and a general increase of the magnitude of the magnetization with the angle $\theta$ but a minimum at the remnant state of the magnetization loop. The prism model presents current paths with 3D bending for all angles $\theta$. The average current density over the thickness agrees very well with the thin film model except for the highest angles. The prism hysteresis loops reveal a peak after the remnant state, which is due to the parallel component of the self-magnetic-field and is implicitly neglected for thin films. The presented numerical method shows the capability to take force-free situations into account for general 3D situations with a high number of degrees of freedom. The results reveal new features of force-free effects in thin films and prisms.


title: Forward Modeling of Spectroscopic Galaxy Surveys: Application to SDSS

abstract: Galaxy spectra are essential to probe the spatial distribution of galaxies in our Universe. To better interpret current and future spectroscopic galaxy redshift surveys, it is important to be able to simulate these data sets. We describe Uspec, a forward modeling tool to generate galaxy spectra taking into account some intrinsic galaxy properties as well as instrumental responses of a given telescope. The model for the intrinsic properties of the galaxy population, i.e., the luminosity functions, and size and spectral coefficients distribu- tions, was developed in an earlier work for broad-band imaging surveys [1], and we now aim to test the model further using spectroscopic data. We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We construct selection cuts that match those used to build this LRG sample, which we then apply to data and simulations in the same way. The resulting real and simulated average spectra show a good statistical agreement overall, with residual differences likely coming from a bluer galaxy population of the simulated sam- ple. We also do not explore the impact of non-solar element ratios in our simulations. For a quantitative comparison, we perform Principal Component Analysis (PCA) of the sets of spectra. By comparing the PCs constructed from simulations and data, we find good agree- ment for all components. The distributions of the eigencoefficients also show an appreciable overlap. We are therefore able to properly simulate the LRG sample taking into account the SDSS/BOSS instrumental responses. The differences between the two samples can be ascribed to the intrinsic properties of the simulated galaxy population, which can be reduced by further improvements of our modelling method in the future. We discuss how these results can be useful for the forward modeling of upcoming large spectroscopic surveys.


title: A Multi-Scheme Ensemble Using Coopetitive Soft-Gating With Application to Power Forecasting for Renewable Energy Generation

abstract: In this article, we propose a novel ensemble technique with a multi-scheme weighting based on a technique called coopetitive soft gating. This technique combines both, ensemble member competition and cooperation, in order to maximize the overall forecasting accuracy of the ensemble. The proposed algorithm combines the ideas of multiple ensemble paradigms (power forecasting model ensemble, weather forecasting model ensemble, and lagged ensemble) in a hierarchical structure. The technique is designed to be used in a flexible manner on single and multiple weather forecasting models, and for a variety of lead times. We compare the technique to other power forecasting models and ensemble techniques with a flexible number of weather forecasting models, which can have the same, or varying forecasting horizons. It is shown that the model is able to outperform those models on a number of publicly available data sets. The article closes with a discussion of properties of the proposed model which are relevant in its application.


title: The Role of Cosmic Ray Transport in Shaping the Simulated Circumgalactic Medium

abstract: The majority of galactic baryons reside outside of the galactic disk in the diffuse gas known as the circumgalactic medium (CGM). While state-of-the art simulations excel at reproducing galactic disk properties, many struggle to drive strong galactic winds or to match the observed ionization structure of the CGM using only thermal supernova feedback. To remedy this, recent studies have invoked non-thermal cosmic ray (CR) stellar feedback prescriptions. However, numerical schemes of CR transport are still poorly constrained. We explore how the choice of CR transport affects the multiphase structure of the simulated CGM. We implement anisotropic CR physics in the astrophysical simulation code, {\sc Enzo} and simulate a suite of isolated disk galaxies with varying prescriptions for CR transport: isotropic diffusion, anisotropic diffusion, and streaming. We find that all three transport mechanisms result in strong, metal-rich outflows but differ in the temperature and ionization structure of their CGM. Isotropic diffusion results in a spatially uniform, warm CGM that underpredicts the column densities of low-ions. Anisotropic diffusion develops a reservoir of cool gas that extends further from the galactic center, but disperses rapidly with distance. CR streaming projects cool gas out to radii of 200 kpc, supporting a truly multiphase medium. In addition, we find that streaming is less sensitive to changes in constant parameter values like the CR injection fraction, transport velocity, and resolution than diffusion. We conclude that CR streaming is a more robust implementation of CR transport and motivate the need for detailed parameter studies of CR transport.


title: Improving performance of SEOBNRv3 by $\sim$300x

abstract: When a gravitational wave is detected by Advanced LIGO/Virgo, sophisticated parameter estimation (PE) pipelines spring into action. These pipelines leverage approximants to generate large numbers of theoretical gravitational waveform predictions to characterize the detected signal. One of the most accurate and physically comprehensive classes of approximants in wide use is the "Spinning Effective One Body--Numerical Relativity" (SEOBNR) family. Waveform generation with these approximants can be computationally expensive, which has limited their usefulness in multiple data analysis contexts. In prior work we improved the performance of the aligned-spin approximant SEOBNR version 2 (v2) by nearly 300x. In this work we focus on optimizing the full eight-dimensional, precessing approximant SEOBNR version 3 (v3). While several v2 optimizations were implemented during its development, v3 is far too slow for use in state-of-the-art source characterization efforts for long-inspiral detections. Completion of a PE run after such a detection could take centuries to complete using v3. Here we develop and implement a host of optimizations for v3, calling the optimized approximant v3_Opt. Our optimized approximant is about 340x faster than v3, and generates waveforms that are numerically indistinguishable.


title: Electroweak Gauge Boson Parton Distribution Functions

abstract: Transverse and longitudinal electroweak gauge boson parton distribution functions (PDFs) are computed in terms of deep-inelastic scattering structure functions, following the recently developed method to determine the photon PDF. The calculation provides initial conditions at the electroweak scale for PDF evolution to higher energies. Numerical results for the $W^\pm$ and $Z$ transverse, longitudinal and polarized PDFs, as well as the $\gamma Z$ transverse and polarized PDFs are presented.


title: Likelihood Non-Gaussianity in Large-Scale Structure Analyses

abstract: Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_\ell$) analysis and the Sinha et al. (2017) group multiplicity function ($\zeta$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_\ell$ and $\zeta$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $f\sigma_8$ constraint by $-0.44\sigma$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $\zeta$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $\alpha$, the posteriors are shifted by $+0.43\sigma$ and $-0.51\sigma$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.


title: Two-dimensional Dirac nodal loop magnons in collinear antiferromagnets

abstract: We study the nontrivial linear magnon band crossings in the collinear antiferromagnets on the two-dimensional (2D) CaVO lattice, also realized in some iron-based superconductors such as AFe$_{1.6+x}$Se$_2$ (A = K, Rb, Cs). It is shown that the combination of space-inversion and time-reversal symmetry ($\mathcal{PT}$-symmetry) leads to doubly-degenerate eight magnon branches, which cross each other linearly along a one-dimensional loop in the 2D Brillouin zone. We show that the Dirac nodal loops (DNLs) are not present in the collinear ferromagnet on this lattice. Thus, the current 2D antiferromagnetic DNLs are symmetry-protected and they provide a novel platform to search for their analogs in 2D electronic antiferromagnetic systems.


title: Energetics of turbulence generated by chiral MHD dynamos

abstract: An asymmetry in the number density of left- and right-handed fermions is known to give rise to a new term in the induction equation that can result in a dynamo instability. At high temperatures, when a chiral asymmetry can survive for long enough, this chiral dynamo instability can amplify magnetic fields efficiently, which in turn drive turbulence via the Lorentz force. While it has been demonstrated in numerical simulations that this chiral magnetically driven turbulence exists and strongly affects the dynamics of the magnetic field, the details of this process remain unclear. The goal of this paper is to analyse the energetics of chiral magnetically driven turbulence and its effect on the generation and dynamics of the magnetic field using direct numerical simulations. We study these effects for different initial conditions, including a variation of the initial chiral chemical potential and the magnetic Prandtl number, Pm. In particular, we determine the ratio of kinetic to magnetic energy, Upsilon, in chiral magnetically driven turbulence. Within the parameter space explored in this study, Upsilon reaches a value of approximately 0.064-0.074 -independently of the initial chiral asymmetry and for Pm=1. Our simulations suggest, that Upsilon decreases as a power law when increasing Pm by decreasing the viscosity. While the exact scaling depends on the details of the fitting criteria and the Reynolds number regime, an approximate result of Upsilon(Pm)=0.1 Pm^{-0.4} is reported. Using the findings from our numerical simulations, we analyse the energetics of chiral magnetically driven turbulence in the early Universe.

categories: physics.flu-dyn astro-ph.CO astro-ph.HE physics.plasm-ph


title: Loop Induced Single Top Partner Production and Decay at the LHC

abstract: Most searches for top partners, $T$, are concerned with top partner pair production. However, as these bounds become increasingly stringent, the LHC energy will saturate and single top partner production will become more important. In this paper we study the novel signature of the top partner produced in association with the SM top, $pp\rightarrow T\overline{t}+t\overline{T}$, in a model where the Standard Model (SM) is extended by a vector-like $SU(2)_L$ singlet fermion top partner and a real, SM gauge singlet scalar, $S$. In this model, $pp\rightarrow T\overline{t}+t\overline{T}$ production is possible through loops mediated by the scalar singlet. We find that, with reasonable coupling strengths, the production rate of this channel can dominate top partner pair production at top partner masses of $m_T\gtrsim 1.5$ TeV. In addition, this model allows for the exotic decay modes $T\rightarrow tg$, $T\rightarrow t \gamma$, and $T\rightarrow t S$. In much of the parameter space the loop induced decay $T\rightarrow tg$ dominates and the top partner is quite long lived. New search strategies are necessary to cover these decay modes. We project the the sensitivity of the high luminosity LHC to $pp\rightarrow T\overline{t}+t\overline{T}$ via a realistic collider study. We find with 3 ab$^{-1}$, the LHC is sensitive to this process for masses $m_T\lesssim2$ TeV. In addition, we provide appendices detailing the renormalization of this model.


title: Age-resolved chemistry of red giants in the solar neighbourhood

abstract: In the age of high-resolution spectroscopic stellar surveys of the Milky Way, the number of stars with detailed abundances of multiple elements is rapidly increasing. These elemental abundances are directly influenced by the evolutionary history of the Galaxy, but this can be difficult to interpret without an absolute timeline of the abundance enrichment. We present age-abundance trends for [M/H], [{\alpha}/M], and 17 individual elements using a sample of 721 solar neighbourhood Hipparcos red giant stars observed by APOGEE. These age trends are determined through a Bayesian hierarchical modelling method presented by Feuillet et al. (2016). We confirm that the [{\alpha}/M]- age relation in the solar neighbourhood is steep and relatively narrow (0.20 dex age dispersion), as are the [O/M]- and [Mg/M]-age relations. The age trend of [C/N] is steep and smooth, consistent with stellar evolution. The [M/H]-age relation has a mean age dispersion of 0.28 dex and a complex overall structure. The oldest stars in our sample are those with the lowest and highest metallicities, while the youngest stars are those with solar metallicity. These results provide strong constraints on theoretical models of Galactic chemical evolution (GCE). We compare them to the predictions of one-zone GCE mod- els and multi-zone mixtures, both analytic and numerical. These comparisons support the hypothesis that the solar neighbourhood is composed of stars born at a range of Galactocentric radii, and that the most metal-rich stars likely migrated from a region with earlier and more rapid star formation such as the inner Galaxy.


title: Potentials for Moduli Spaces of A_m-local Systems on Surfaces

abstract: We study properties of potentials on quivers $Q_{\mathcal{T},m}$ arising from cluster coordinates on moduli spaces of $PGL_{m+1}$-local systems on a topological surface with punctures. To every quiver with potential one can associate a $3d$ Calabi-Yau $A_\infty$-category in such a way that a natural notion of equivalence for quivers with potentials (called "right-equivalence") translates to $A_\infty$-equivalence of associated categories. For any quiver one can define a notion of a "primitive" potential. Our first result is the description of the space of equivalence classes of primitive potentials on quivers $Q_{\mathcal{T}, m}$. Then we provide a full description of the space of equivalence classes of all \emph{generic} potentials for the case $m = 2$ (corresponds to $PGL_3$-local systems). In particular, we show that it is finite-dimensional. This claim extends results of Gei\ss, Labardini-Fragoso and Schr\"oer who have proved analogous statement in $m=1$ case. In many cases $3d$ Calabi-Yau $A_\infty$-categories constructed from quivers with potentials are expected to be realized geometrically as Fukaya categories of certain Calabi-Yau $3$-folds. Bridgeland and Smith gave an explicit construction of Fukaya categories for quivers $Q_{\mathcal{T},m=1}$. We propose a candidate for Calabi-Yau $3$-folds that would play analogous role in higher rank cases, $m > 1$. We study their (co)homology and describe a construction of collections of $3$-dimensional spheres that should play a role of generating collections of Lagrangian spheres in corresponding Fukaya categories.


title: Serverless Data Analytics with Flint

abstract: Serverless architectures organized around loosely-coupled function invocations represent an emerging design for many applications. Recent work mostly focuses on user-facing products and event-driven processing pipelines. In this paper, we explore a completely different part of the application space and examine the feasibility of analytical processing on big data using a serverless architecture. We present Flint, a prototype Spark execution engine that takes advantage of AWS Lambda to provide a pure pay-as-you-go cost model. With Flint, a developer uses PySpark exactly as before, but without needing an actual Spark cluster. We describe the design, implementation, and performance of Flint, along with the challenges associated with serverless analytics.


title: A Low-rank Tensor Regularization Strategy for Hyperspectral Unmixing

abstract: Tensor-based methods have recently emerged as a more natural and effective formulation to address many problems in hyperspectral imaging. In hyperspectral unmixing (HU), low-rank constraints on the abundance maps have been shown to act as a regularization which adequately accounts for the multidimensional structure of the underlying signal. However, imposing a strict low-rank constraint for the abundance maps does not seem to be adequate, as important information that may be required to represent fine scale abundance behavior may be discarded. This paper introduces a new low-rank tensor regularization that adequately captures the low-rank structure underlying the abundance maps without hindering the flexibility of the solution. Simulation results with synthetic and real data show that the the extra flexibility introduced by the proposed regularization significantly improves the unmixing results.


title: The Spitzer-IRAC/MIPS Extragalactic Survey (SIMES): II enhanced nuclear accretion rate in galaxy groups at z$\sim$0.2

abstract: For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups.


title: Maximal subalgebras of the exceptional Lie algebras in low characteristic

abstract: In this thesis we consider the maximal subalgebras of the exceptional Lie algebras in algebraically closed fields of positive characteristic. This begins with a quick recap of the article by Herpel and Stewart which considered the Cartan type maximal subalgebras in the exceptional Lie algebras for good characteristic, and then the article by Premet considering non-semisimple maximal subalgebras in good characteristic. For $p=5$ we give an example of what appears to be a new maximal subalgebra in the exceptional Lie algebra of type $E_8$. We show that this maximal subalgebra is isomorphic to the $p$-closure of the non-restricted Witt algebra $W(1;2)$. After this, we focus completely on characteristics $p=2$ and $p=3$ giving examples of new non-semisimple maximal subalgebras in the exceptional Lie algebras. We consider the Weisfeiler filtration associated to these maximal subalgebras and leave many open questions. There are one or two examples of simple maximal subalgebras in $F_4$ for $p=3$ and $E_8$ for $p=2$.


abstract: Self-trapping and acceleration of ions in laser-driven relativistically transparent plasma are investigated with the help of particle-in-cell simulations. A theoretical model based on ion wave breaking is established in describing ion evolution and ion trapping. The threshold for ion trapping is identified. Near the threshold ion trapping is self-regulating and stops when the number of trapped ions is large enough. The model is applied to ion trapping in three-dimensional geometry. Longitudinal distributions of ions and the electric field near the wave breaking point are derived analytically in terms of power-law scalings. The areal density of trapped charge is obtained as a function of the strength of ion wave breaking, which scales with target density for fixed laser intensity. The results of the model are confirmed by the simulations.


title: SMSS J130522.47-293113.0: a high-latitude stellar X-ray source with pc-scale outflow relics?

abstract: We report the discovery of an unusual stellar system, SMSS J130522.47-293113.0. The optical spectrum is dominated by a blue continuum together with emission lines of hydrogen, neutral and ionized helium, and the N III, C III blend at around 4640-4650 Angstrom. The emission line profiles vary in strength and position on timescales as short as 1 day, while optical photometry reveals fluctuations of as much as ~0.2 mag in g on timescales as short as 10-15 min. The system is a weak X-ray source (f_{0.3-10} = (1.2 +/- 0.1) x 10^{-13} ergs cm^{-2} s^{-1} in the 0.3-10 keV band) but is not detected at radio wavelengths (3-sigma upper limit of 50 microJy at 5.5 GHz). The most intriguing property of the system, however, is the existence of two "blobs", a few arcsec in size, that are symmetrically located 3.8 arcmin (2.2 pc for our preferred system distance of ~2 kpc) each side of the central object. The blobs are detected in optical and near-IR broadband images but do not show any excess emission in Halpha images. We discuss the interpretation of the system, suggesting that the central object is most likely a nova-like CV, and that the blobs are relics of a pc-scale accretion-powered collimated outflow.


title: Transport information geometry I: Riemannian calculus on probability simplex

abstract: We formulate the Riemannian calculus of the probability set embedded with $L^2$-Wasserstein metric. This is an initial work of transport information geometry. Our investigation starts with the probability simplex (probability manifold) supported on vertices of a finite graph. The main idea is to embed the probability manifold as a submanifold of the positive measure space with a nonlinear metric tensor. Here the nonlinearity comes from the linear weighted Laplacian operator. By this viewpoint, we establish torsion-free Christoffel symbols, Levi-Civita connections, curvature tensors and volume forms in the probability manifold by Euclidean coordinates. As a consequence, the Jacobi equation, Laplace-Beltrami and Hessian operators on the probability manifold are derived. These geometric computations are also provided in the infinite-dimensional density space (density manifold) supported on a finite-dimensional manifold. In particular, an identity is given connecting the Baker-{\'E}mery $\Gamma_2$ operator (carr{\'e} du champ it{\'e}r{\'e}) by connecting Fisher-Rao information metric and optimal transport metric. Several examples are demonstrated.


title: Halving the bounds for the Markov, Chebyshev, and Chernoff Inequalities using smoothing

abstract: The Markov, Chebyshev, and Chernoff inequalities are some of the most widely used methods for bounding the tail probabilities of random variables. In all three cases, the bounds are tight in the sense that there exists easy examples where the inequalities become equality. Here we will show that through a simple smoothing using auxiliary randomness, that each of the three bounds can be cut in half. In many common cases, the halving can be achieved without the need for the auxiliary randomness.


title: A Candidate Tidal Disruption Event in a Quasar at z=2.359 from Abundance Ratio Variability

abstract: A small fraction of quasars show an unusually high nitrogen-to-carbon ratio (N/C) in their spectra. These "nitrogen-rich" (N-rich) quasars are a long-standing puzzle because their interstellar medium implies stellar populations with abnormally high metallicities. It has recently been proposed that N-rich quasars may result from tidal disruption events (TDEs) of stars by supermassive black holes. The rapid enhancement of nitrogen and the depletion of carbon due to the carbon--nitrogen--oxygen cycle in supersolar mass stars could naturally produce high N/C. However, the TDE hypothesis predicts that the N/C should change with time, which has never hitherto been observed. Here we report the discovery of the first N-rich quasar with rapid N/C variability that could be caused by a TDE. Two spectra separated by 1.7 years (rest-frame) show that the N III]\lambda 1750/C III]\lambda 1909 intensity ratio decayed by ~86%+/-14% (1 \sigma). Optical (rest-frame UV) light-curve and X-ray observations are qualitatively consistent with the TDE hypothesis; though, the time baseline falls short of a definitive proof. Putting the single-object discovery into context, statistical analyses of the ~80 known N-rich quasars with high-quality archival spectra show evidence (at a 5\sigma\ significance level) of a decrease in N/C on timescales of >1 year (rest-frame) and a constant level of ionization (indicated by the C III]\lambda 1909/C IV\lambda 1549 intensity ratio). If confirmed, our results demonstrate the method of identifying TDE candidates in quasars via abundance ratio variability, opening a new window of TDE observations at high redshift (z>2) with upcoming large-scale time-domain spectroscopic surveys.


title: Geometric Adaptive Control for a Quadrotor UAV with Wind Disturbance Rejection

abstract: This paper presents a geometric adaptive control scheme for a quadrotor unmanned aerial vehicle, where the effects of unknown, unstructured disturbances are mitigated by a multilayer neural network that is adjusted online. The stability of the proposed controller is analyzed with Lyapunov stability theory on the special Euclidean group, and it is shown that the tracking errors are uniformly ultimately bounded with an ultimate bound that can be abridged arbitrarily. A mathematical model of wind disturbance on the quadrotor dynamics is presented, and it is shown that the proposed adaptive controller is capable of rejecting the effects of wind disturbances successfully. These are illustrated by numerical examples.


title: Conversion rules for Weyl points and nodal lines in topological media

abstract: According to a widely-held paradigm, a pair of Weyl points with opposite chirality mutually annihilate when brought together. In contrast, we show that such a process is strictly forbidden for Weyl points related by a mirror symmetry, provided that an effective two-band description exists in terms of orbitals with opposite mirror eigenvalue. Instead, such a pair of Weyl points convert into a nodal loop inside a symmetric plane upon the collision. Similar constraints are identified for systems with multiple mirrors, facilitating previously unreported nodal-line and nodal-chain semimetals that exhibit both Fermi-arc and drumhead surface states. We further find that Weyl points in systems symmetric under a $\pi$-rotation composed with time-reversal are characterized by an additional integer charge that we call helicity. A pair of Weyl points with opposite chirality can annihilate only if their helicities also cancel out. We base our predictions on topological crystalline invariants derived from relative homotopy theory, and we test our predictions on simple tight-binding models. The outlined homotopy description can be directly generalized to systems with multiple bands and other choices of symmetry.


title: Inference for case-control studies with incident and prevalent cases

abstract: We propose and study a fully efficient method to estimate associations of an exposure with disease incidence when both, incident cases and prevalent cases, i.e. individuals who were diagnosed with the disease at some prior time point and are alive at the time of sampling, are included in a case-control study. We extend the exponential tilting model for the relationship between exposure and case status to accommodate two case groups, and correct for the survival bias in the prevalent cases through a tilting term that depends on the parametric distribution of the backward time, i.e. the time from disease diagnosis to study enrollment. We construct an empirical likelihood that also incorporates the observed backward times for prevalent cases, obtain efficient estimates of odds ratio parameters that relate exposure to disease incidence and propose a likelihood ratio test for model parameters that has a standard chi-squared distribution. We quantify the changes in efficiency of association parameters when incident cases are supplemented with, or replaced by, prevalent cases in simulations. We illustrate our methods by estimating associations of single nucleotide polymorphisms (SNPs) with breast cancer incidence in a sample of controls and incident and prevalent cases from the U.S. Radiologic Technologists Health Study.


title: Graphs, Ultrafilters and Colourability

abstract: Let $\beta$ be the functor from Set to CHaus which maps each discrete set X to its Stone-Cech compactification, the set $\beta$ X of ultrafilters on X. Every graph G with vertex set V naturally gives rise to a graph $\beta G$ on the set $\beta V$ of ultrafilters on V . In what follows, we interrelate the properties of G and $\beta G$. Perhaps the most striking result is that G can be finitely coloured iff $\beta G$ has no loops.


title: The supernova remnant population in the very-high-energy sky: prospects for CTA

abstract: The detection of very-high-energy gamma rays from supernova remnant shells testifies of the acceleration of particles at strong shocks. Many aspects of the particle acceleration remain however unclear. The study of individual objects is very helpful, but the study of the entire population of SNRs detected in this range and its characteristics can also bring valuable science. Using Monte-Carlo simulations, the population of shells bright in the TeV and multi-TeV range can be simulated. The results of these simulations aim at being compared with observations of in struments operating in these ranges, such as the Cherenkov Telescope Array (CTA). Our results suggest that CTA should be able to effectively constrain the slope of particles accelerated at SNRs and the electron-to-proton ratio.


title: Doppelg\"anger dark energy: modified gravity with non-universal couplings after GW170817

abstract: Gravitational Wave (GW) astronomy severely narrowed down the theoretical space for scalar-tensor theories. We propose a new class of attractor models for Horndeski action in which GWs propagate at the speed of light in the nearby universe but not in the past. To do so we derive new solutions to the interacting dark sector in which the ratio of dark energy and dark matter remains constant, which we refer to as doppelg\"anger dark energy (DDE). We then remove the interaction between dark matter and dark energy by a suitable change of variables. The accelerated expansion that (we) baryons observe is due to a conformal coupling to the dark energy scalar field. We show how in this context it is possible to find a non trivial subset of solutions in which GWs propagate at the speed of light only at low red-shifts. The model is an attractor, thus reaching the limit $c_{T}\to1$ relatively fast. However, the effect of baryons turns out to be non-negligible and severely constrains the form of the Lagrangian. In passing, we found that in the simplest DDE models the no-ghost conditions for perturbations require a non-universal coupling to gravity. In the end, we comment on possible ways to solve the lack of matter domination stage for DDE models.


title: A Family of Minimal and Renormalizable Rectangle Exchange Maps

abstract: A domain exchange map (DEM) is a dynamical system defined on a smooth Jordan domain which is a piecewise translation. We explain how to use cut-and-project sets to construct minimal DEMs. Specializing to the case in which the domain is a square and the cut-and-project set is associated to a Galois lattice, we construct an infinite family of DEMs in which each map is associated to a PV number. We develop a renormalization scheme for these DEMs. Certain DEMs in the family can be composed to create multistage, renormalizable DEMs.


title: First Calorimetric Measurement of Electron Capture in ${}^{193}$Pt with a Transition Edge Sensor

abstract: The neutrino mass can be extracted from a high statistics, high resolution calorimetric spectrum of electron capture in ${}^{163}$Ho. In order to better understand the shape of the calorimetric electron capture spectrum, a second isotope was measured with a close to ideal absorber-source configuration. ${}^{193}$Pt was created by irradiating a ${}^{192}$Pt-enriched platinum foil in a nuclear reactor. This Pt-in-Pt absorber was designed to have a nearly ideal absorber-source configuration. The measured ${}^{193}$Pt calorimetric electron-capture spectrum provides an independent check on the corresponding theoretical calculations, which have thus far been compared only for ${}^{163}$Ho. The first experimental and theoretically-calculated spectra from this ${}^{193}$Pt-in-Pt absorber are presented and overlaid for preliminary comparison of theory with experiment.


title: Ludvig Lorenz (1867) on Light and Electricity

abstract: Independent of Maxwell, in 1867 the Danish physicist L. V. Lorenz proposed a theory in which he identified light with electrical oscillations propagating in a very poor conductor. Lorenz's electrodynamic theory of light, which formally was equivalent to Maxwell's theory but physically quite different from it, was published in well-known journals in German and English but soon fell into oblivion. In 1867 Lorenz also published a paper on his new theory in a semi-popular Danish journal which has generally been overlooked. This other paper is here translated into English and provided with the necessary annotations.


title: Stochastic basins of attraction and generalized committor functions

abstract: We generalize the concept of basin of attraction of a stable state in order to facilitate the analysis of dynamical systems with noise and to assess stability properties of metastable states and long transients. To this end we examine the notions of mean sojourn times and absorption probabilities for Markov chains and study their relation to the basins of attraction. Our approach is applicable to a large variety of problems since in most cases the transfer operator associated to a dynamical system can be approximated by a Markov chain.


title: Adversarial Logit Pairing

abstract: In this paper, we develop improved techniques for defending against adversarial examples at scale. First, we implement the state of the art version of adversarial training at unprecedented scale on ImageNet and investigate whether it remains effective in this setting - an important open scientific question (Athalye et al., 2018). Next, we introduce enhanced defenses using a technique we call logit pairing, a method that encourages logits for pairs of examples to be similar. When applied to clean examples and their adversarial counterparts, logit pairing improves accuracy on adversarial examples over vanilla adversarial training; we also find that logit pairing on clean examples only is competitive with adversarial training in terms of accuracy on two datasets. Finally, we show that adversarial logit pairing achieves the state of the art defense on ImageNet against PGD white box attacks, with an accuracy improvement from 1.5% to 27.9%. Adversarial logit pairing also successfully damages the current state of the art defense against black box attacks on ImageNet (Tramer et al., 2018), dropping its accuracy from 66.6% to 47.1%. With this new accuracy drop, adversarial logit pairing ties with Tramer et al.(2018) for the state of the art on black box attacks on ImageNet.


title: Ab Initio Electron-Phonon Interactions Using Atomic Orbital Wavefunctions

abstract: The interaction between electrons and lattice vibrations determines key physical properties of materials, including their electrical and heat transport, excited electron dynamics, phase transitions, and superconductivity. We present a new ab initio method that employs atomic orbital (AO) wavefunctions to compute the electron-phonon (e-ph) interactions in materials and interpolate the e-ph coupling matrix elements to fine Brillouin zone grids. We detail the numerical implementation of such AO-based e-ph calculations, and benchmark them against direct density functional theory calculations and Wannier function (WF) interpolation. The key advantages of AOs over WFs for e-ph calculations are outlined. Since AOs are fixed basis functions associated with the atoms, they circumvent the need to generate a material-specific localized basis set with a trial-and-error approach, as is needed in WFs. Therefore, AOs are ideal to compute e-ph interactions in chemically and structurally complex materials for which WFs are challenging to generate, and are also promising for high-throughput materials discovery. While our results focus on AOs, the formalism we present generalizes e-ph calculations to arbitrary localized basis sets, with WFs recovered as a special case.


title: Mobile phone records to feed activity-based travel demand models: MATSim for studying a cordon toll policy in Barcelona

abstract: Activity-based models appeared as an answer to the limitations of the traditional trip-based and tour-based four-stage models. The fundamental assumption of activity-based models is that travel demand is originated from people performing their daily activities. This is why they include a consistent representation of time, of the persons and households, time-dependent routing, and microsimulation of travel demand and traffic. In spite of their potential to simulate traffic demand management policies, their practical application is still limited. One of the main reasons is that these models require a huge amount of very detailed input data hard to get with surveys. However, the pervasive use of mobile devices has brought a valuable new source of data. The work presented here has a twofold objective: first, to demonstrate the capability of mobile phone records to feed activity-based transport models, and, second, to assert the advantages of using activity-based models to estimate the effects of traffic demand management policies. Activity diaries for the metropolitan area of Barcelona are reconstructed from mobile phone records. This information is then employed as input for building a transport MATSim model of the city. The model calibration and validation process proves the quality of the activity diaries obtained. The possible impacts of a cordon toll policy applied to two different areas of the city and at different times of the day is then studied. Our results show the way in which the modal share is modified in each of the considered scenario. The possibility of evaluating the effects of the policy at both aggregated and traveller level, together with the ability of the model to capture policy impacts beyond the cordon toll area confirm the advantages of activity-based models for the evaluation of traffic demand management policies.


title: A Generalised Method for Empirical Game Theoretic Analysis

abstract: This paper provides theoretical bounds for empirical game theoretical analysis of complex multi-agent interactions. We provide insights in the empirical meta game showing that a Nash equilibrium of the meta-game is an approximate Nash equilibrium of the true underlying game. We investigate and show how many data samples are required to obtain a close enough approximation of the underlying game. Additionally, we extend the meta-game analysis methodology to asymmetric games. The state-of-the-art has only considered empirical games in which agents have access to the same strategy sets and the payoff structure is symmetric, implying that agents are interchangeable. Finally, we carry out an empirical illustration of the generalised method in several domains, illustrating the theory and evolutionary dynamics of several versions of the AlphaGo algorithm (symmetric), the dynamics of the Colonel Blotto game played by human players on Facebook (symmetric), and an example of a meta-game in Leduc Poker (asymmetric), generated by the PSRO multi-agent learning algorithm.


title: Spread of Information with Confirmation Bias in Cyber-Social Networks

abstract: This paper provides a model to investigate information spreading over cyber-social network of agents communicating with each other. The cyber-social network considered here comprises individuals and news agencies. Each individual holds a belief represented by a scalar. Individuals receive information from news agencies that are closer to their belief, confirmation bias is explicitly incorporated into the model. The proposed dynamics of cyber-social networks is adopted from DeGroot-Friedkin model, where the individual's opinion update mechanism is a convex combination of his innate opinion, his neighbors' opinions at the previous time step (obtained from the social network), and the opinions passed along by news agencies from cyber layer which he follows. The characteristics of the interdependent social and cyber networks are radically different here: the social network relies on trust and hence static while the news agencies are highly dynamic since they are weighted as a function of the distance between an individual state and the state of news agency to account for confirmation bias. The conditions for convergence of the aforementioned dynamics to a unique equilibrium are characterized. The estimation and exact computation of the steady-state values under non-linear and linear state-dependent weight functions are provided. Finally, the impact of polarization in the opinions of news agencies on the public opinion evolution is numerically analyzed in the context of the well-known Krackhardt's advice network.


title: Characterizing Ion Flows Across a Dipolarization Jet

abstract: The structure of dipolarization jets with finite width in the dawn-dusk direction relevant to magnetic reconnection in the Earth's magnetotail is explored with particle-in-cell simulations. We carry out Riemann simulations of the evolution of the jet in the dawn-dusk, north-south plane to investigate the dependence of the jet structure on the jet width in the dawn-dusk direction. We find that the magnetic field and Earth-directed ion flow structure depend on the dawn-dusk width. A reversal in the usual Hall magnetic field near the center of the current sheet on the dusk side of larger jets is observed. For small widths, the maximum velocity of the Earthward flow is significantly reduced below the theoretical limit of the upstream Alfv\'en speed. However, the ion flow speed approaches this limit once the width exceeds the ion Larmor radius based on the normal magnetic field, $B_z$.


title: Photon detector system timing performance in the DUNE 35-ton prototype liquid argon time projection chamber

abstract: The 35-ton prototype for the Deep Underground Neutrino Experiment far detector was a single-phase liquid argon time projection chamber with an integrated photon detector system, all situated inside a membrane cryostat. The detector took cosmic-ray data for six weeks during the period of February 1, 2016 to March 12, 2016. The performance of the photon detection system was checked with these data. An installed photon detector was demonstrated to measure the arrival times of cosmic-ray muons with a resolution better than 32 ns, limited by the timing of the trigger system. A measurement of the timing resolution using closely-spaced calibration pulses yielded a resolution of 15 ns for pulses at a level of 6 photo-electrons. Scintillation light from cosmic-ray muons was observed to be attenuated with increasing distance with a characteristic length of $155 \pm 28$ cm.


title: Distributed Optimization for Second-Order Multi-Agent Systems with Dynamic Event-Triggered Communication

abstract: In this paper, we propose a fully distributed algorithm for second-order continuous-time multi-agent systems to solve the distributed optimization problem. The global objective function is a sum of private cost functions associated with the individual agents and the interaction between agents is described by a weighted undirected graph. We show the exponential convergence of the proposed algorithm if the underlying graph is connected, each private cost function is locally gradient-Lipschitz-continuous, and the global objective function is restricted strongly convex with respect to the global minimizer. Moreover, to reduce the overall need of communication, we then propose a dynamic event-triggered communication mechanism that is free of Zeno behavior. It is shown that the exponential convergence is achieved if the private cost functions are also globally gradient-Lipschitz-continuous. Numerical simulations are provided to illustrate the effectiveness of the theoretical results.


title: An Edge Computing Empowered Radio Access Network With UAV-Mounted FSO Fronthaul and Backhaul: Key Challenges and Approaches

abstract: One promising approach to address the supply-demand mismatch between the terrestrial infrastructure and the temporary and/or unexpected traffic demands is to leverage the unmanned aerial vehicle (UAV) technologies. Motivated by the recent advancement of UAV technologies and retromodulator based free space optical communication, we propose a novel edge-computing empowered radio access network architecture where the fronthaul and backhaul links are mounted on the UAVs for rapid event response and flexible deployment. The implementation of hardware and networking technologies for the proposed architecture are investigated. Due to the limited payload and endurance as well as the high mobility of UAVs, research challenges related to the communication resource management and recent research progress are reported.


title: Harmonic spinors on the Davis hyperbolic 4-manifold

abstract: In this paper we use the G-spin theorem to show that the Davis hyperbolic 4-manifold admits harmonic spinors. This is the first example of a closed hyperbolic 4-manifold that admits harmonic spinors. We also explicitly describe the Spinor bundle of a spin hyperbolic 2- or 4-manifold and show how to calculated the subtle sign terms in the G-spin theorem for an isometry, with isolated fixed points, of a closed spin hyperbolic 2- or 4-manifold.


title: A prediction criterion for working correlation structure selection in GEE

abstract: Generalized estimating equations (GEE) is one of the most commonly used methods for marginal regression analysis of longitudinal data, especially with discrete outcomes. The GEE method models the association among the responses of a subject through a working correlation matrix and correct specification of the working correlation structure ensures efficient estimation of the regression parameters. This study proposes a predicted residual sum of squares (PRESS) statistic as a working correlation selection criterion in GEE. An extensive simulation study is designed to assess the performance of the proposed GEE PRESS criterion and to compare its performance with well-known existing criteria in the literature. The results show that the GEE PRESS criterion has better performance than the weighted error sum of squares SC criterion in all cases and is comparable with that of other existing criteria when the true working correlation structure is AR(1) or exchangeable. Lastly, the working correlation selection criteria are illustrated with the Coronary Artery Risk Development in Young Adults study.