title: Quasi-phases and pseudo-transitions in one-dimensional models with nearest neighbor interactions

abstract: There are some particular one-dimensional models, such as the Ising-Heisenberg spin models with a variety of chain structures, which exhibit unexpected behaviors quite similar to the first and second order phase transition, which could be confused naively with an authentic phase transition. Through the analysis of the first derivative of free energy, such as entropy, magnetization, and internal energy, a "sudden" jump that closely resembles a first-order phase transition at finite temperature occurs. However, by analyzing the second derivative of free energy, such as specific heat and magnetic susceptibility at finite temperature, it behaves quite similarly to a second-order phase transition exhibiting an astonishingly sharp and fine peak. The correlation length also confirms the evidence of this pseudo-transition temperature, where a sharp peak occurs at the pseudo-critical temperature. We also present the necessary conditions for the emergence of these quasi-phases and pseudo-transitions.


title: Fastron: An Online Learning-Based Model and Active Learning Strategy for Proxy Collision Detection

abstract: We introduce the Fastron, a configuration space (C-space) model to be used as a proxy to kinematic-based collision detection. The Fastron allows iterative updates to account for a changing environment through a combination of a novel formulation of the kernel perceptron learning algorithm and an active learning strategy. Our simulations on a 7 degree-of-freedom arm indicate that proxy collision checks may be performed at least 2 times faster than an efficient polyhedral collision checker and at least 8 times faster than an efficient high-precision collision checker. The Fastron model provides conservative collision status predictions by padding C-space obstacles, and proxy collision checking time does not scale poorly as the number of workspace obstacles increases. All results were achieved without GPU acceleration or parallel computing.


title: Computing optimal experimental designs with respect to a compound Bayes risk criterion

abstract: We consider the problem of computing optimal experimental design on a finite design space with respect to a compound Bayes risk criterion, which includes the linear criterion for prediction in a random coefficient regression model. We show that the problem can be restated as constrained A-optimality in an artificial model. This permits using recently developed computational tools, for instance the algorithms based on the second-order cone programming for optimal approximate design, and mixed-integer second-order cone programming for optimal exact designs. We demonstrate the use of the proposed method for the problem of computing optimal designs of a random coefficient regression model with respect to an integrated mean squared error criterion.


title: Lattice Surgery with a Twist: Simplifying Clifford Gates of Surface Codes

abstract: We present a planar surface-code-based scheme for fault-tolerant quantum computation which eliminates the time overhead of single-qubit Clifford gates, and implements long-range multi-target CNOT gates with a time overhead that scales only logarithmically with the control-target separation. This is done by replacing hardware operations for single-qubit Clifford gates with a classical tracking protocol. Inter-qubit communication is added via a modified lattice surgery protocol that employs twist defects of the surface code. The long-range multi-target CNOT gates facilitate magic state distillation, which renders our scheme fault-tolerant and universal.


title: Calculating the Expected Value of Sample Information using Efficient Nested Monte Carlo: A Tutorial

abstract: Objective: The Expected Value of Sample Information (EVSI) quantifies the economic benefit of reducing uncertainty in a health economic model by collecting additional information. This has the potential to improve the allocation of research budgets. Despite this, practical EVSI evaluations are limited, partly due to the computational cost of estimating this value using the "gold-standard" nested simulation methods. Recently, however, Heath et al developed an estimation procedure that reduces the number of simulations required for this "gold-standard" calculation. Up to this point, this new method has been presented in purely technical terms. Study Design: This study presents the practical application of this new method to aid its implementation. We use a worked example to illustrate the key steps of the EVSI estimation procedure before discussing its optimal implementation using a practical health economic model. Methods: The worked example is based on a three parameter linear health economic model. The more realistic model evaluates the cost-effectiveness of a new chemotherapy treatment which aims to reduce the number of side effects experienced by patients. We use a Markov Model structure to evaluate the health economic profile of experiencing side effects. Results: This EVSI estimation method offers accurate estimation within a feasible computation time, seconds compared to days, even for more complex model structures. The EVSI estimation is more accurate if a greater number of nested samples are used, even for a fixed computational cost. Conclusions: This new method reduces the computational cost of estimating the EVSI by nested simulation.


title: The formation of a Spitzer bubble RCW79 triggered by cloud-cloud collision

abstract: Understanding the mechanism of O star formation is one of the most important issues in current astrophysics. It is also an issue of keen interest how O stars affect their surroundings and trigger secondary star formation. An H\,\emissiontype{II} region RCW79 is one of the typical Spitzer bubbles alongside of RCW120. New observations of CO $J=$ 1--0 emission with Mopra and NANTEN2 revealed that molecular clouds are associated with RCW79 in four velocity components over a velocity range of 20 km s$^{-1}$. We hypothesize that two of the clouds collided with each other and the collision triggered the formation of 12 O stars inside of the bubble and the formation of 54 low mass young stellar objects along the bubble wall. The collision is supported by observational signatures of bridges connecting different velocity components in the colliding clouds. The whole collision process happened in a timescale of $\sim$1 Myr. RCW79 has a larger size by a factor of 30 in the projected area than RCW120 with a single O star, and the large size favored formation of the 12 O stars due to the larger accumulated gas in the collisional shock compression.


title: From Electrodiffusion Theory to the Electrohydrodynamics of Leaky Dielectrics through the Weak Electrolyte Limit

abstract: The Taylor-Melcher (TM) model is the standard model for describing the dynamics of poorly conducting leaky dielectric fluids under an electric field. The TM model treats the fluids as Ohmic conductors, without modeling the underlying ion dynamics. On the other hand, electrodiffusion models, which have been successful in describing electrokinetic phenomena, incorporate ionic concentration dynamics. Mathematical reconciliation of the electrodiffusion picture and the TM model has been a major issue for electrohydrodynamic theory. Here, we derive the TM model from an electrodiffusion model in which we explicitly model the electrochemistry of ion dissociation. We introduce salt dissociation reaction terms in the bulk electrodiffusion equations and take the limit in which the salt dissociation is weak; the assumption of weak dissociation corresponds to the fact that the TM model describes poor conductors. Together with the assumption that the Debye length is small, we derive the TM model with or without the surface charge convection term depending upon the scaling of relevant dimensionless parameters. An important quantity that emerges is the Galvani potential (GP), the jump in voltage across the liquid-liquid interface between the two leaky dielectric media; the GP arises as a natural consequence of the interfacial boundary conditions for the ionic concentrations, and is absent under certain parametric conditions. When the GP is absent, we recover the TM model. Our analysis also reveals the structure of the Debye layer at the liquid-liquid interface, which suggests how interfacial singularities may arise under strong imposed electric fields. In the presence of a non-zero GP, our model predicts that the liquid droplet will drift under an imposed electric field, the velocity of which is computed explicitly to leading order.


title: Competing effects of Hund's splitting and symmetry-breaking perturbations on electronic order in Pb$_{1-x}$Sn$_{x}$Te

abstract: We study the effect of a uniform external magnetization on p-wave superconductivity on the (001) surface of the crystalline topological insulator(TCI) Pb$_{1-x}$Sn$_{x}$Te. It was shown by us in an earlier work that a chiral p-wave finite momentum pairing (FFLO) state can be stabilized in this system in the presence of weak repulsive interparticle interactions. In particular, the superconducting instability is very sensitive to the Hund's interaction in the multiorbital TCI, and no instabilities are found to be possible for the "wrong" sign of the Hund's splitting. Here we show that for a finite Hund's splitting of interactions, a significant value of the external magnetization is needed to degrade the surface superconductivity, while in the absence of the Hund's interaction, an arbitrarily small external magnetization can destroy the superconductivity. This implies that multiorbital effects in this system play an important role in stabilizing electronic order on the surface.


title: A computational hierarchy in human cortex

abstract: Hierarchies feature prominently in anatomical accounts of cortical organisation. An open question is which computational (algorithmic) processes are implemented by these hierarchies. One renowned hypothesis is that cortical hierarchies implement a model of the world's causal structure and serve to infer environmental states from sensory inputs. This view, which casts perception as hierarchical Bayesian inference, has become a highly influential concept in both basic and clinical neuroscience. So far, however, a direct correspondence between the predicted order of hierarchical Bayesian computations and the sequence of evoked neuronal activity has not been demonstrated. Here, we present evidence for this correspondence from neuroimaging and electrophysiological data in healthy volunteers. Trial-wise sequences of hierarchical computations were inferred from participants' behaviour during a social learning task that required multi-level inference about intentions. We found that the temporal sequence of neuronal activity matched the order of computations as predicted by the theory. These findings provide strong evidence for the operation of hierarchical Bayesian inference in human cortex. Furthermore, our approach offers a novel strategy for the combined computational-physiological phenotyping of patients with disorders of perception, such as schizophrenia or autism.


title: Measurements of Physical Parameters of White Dwarfs: A Test of the Mass-Radius Relation

abstract: We present a detailed spectroscopic and photometric analysis of 219 DA and DB white dwarfs for which trigonometric parallax measurements are available. Our aim is to compare the physical parameters derived from the spectroscopic and photometric techniques, and then to test the theoretical mass-radius relation for white dwarfs using these results. The agreement between spectroscopic and photometric parameters is found to be excellent, especially for effective temperatures, showing that our model atmospheres and fitting procedures provide an accurate, internally consistent analysis. Values of surface gravity and solid angle, obtained respectively from spectroscopy and photometry, are combined with parallax measurements in various ways to study the validity of the mass-radius relation from an empirical point of view. After a thorough examination of our results, we find that 73% and 92% of the white dwarfs are consistent within 1 and 2$\sigma$ confidence levels, respectively, with the predictions of the mass-radius relation, thus providing strong support to the theory of stellar degeneracy. Our analysis also allows us to identify 15 stars that are better interpreted in terms of unresolved double degenerate binaries. Atmospheric parameters for both components in these binary systems are obtained using a novel approach. We further identify a few white dwarfs that are possibly composed of an iron core rather than a carbon/oxygen core, since they are consistent with Fe-core evolutionary models.


title: Multilayer Brain Networks

abstract: The field of neuroscience is facing an unprecedented expanse in the volume and diversity of available data. Traditionally, network models have provided key insights into the structure and function of the brain. With the advent of big data in neuroscience, both more sophisticated models capable of characterizing the increasing complexity of the data and novel methods of quantitative analysis are needed. Recently multilayer networks, a mathematical extension of traditional networks, have gained increasing popularity in neuroscience due to their ability to capture the full information of multi-model, multi-scale, spatiotemporal data sets. Here, we review multilayer networks and their applications in neuroscience, showing how incorporating the multilayer framework into network neuroscience analysis has uncovered previously hidden features of brain networks. We specifically highlight the use of multilayer networks to model disease, structure-function relationships, network evolution, and link multi-scale data. Finally, we close with a discussion of promising new directions of multilayer network neuroscience research and propose a modified definition of multilayer networks designed to unite and clarify the use of the multilayer formalism in describing real-world systems.


title: Finding bases of uncountable free abelian groups is usually difficult

abstract: We investigate effective properties of uncountable free abelian groups. We show that identifying free abelian groups and constructing bases for such groups is often computationally hard, depending on the cardinality. For example, we show, under the assumption $V=L$, that there is a first-order definable free abelian group with no first-order definable basis.


title: Feature selection in high-dimensional dataset using MapReduce

abstract: This paper describes a distributed MapReduce implementation of the minimum Redundancy Maximum Relevance algorithm, a popular feature selection method in bioinformatics and network inference problems. The proposed approach handles both tall/narrow and wide/short datasets. We further provide an open source implementation based on Hadoop/Spark, and illustrate its scalability on datasets involving millions of observations or features.


title: Plasmonic nanolaser for intracavity spectroscopy and sensorics

abstract: We demonstrate intracavity plasmonic laser spectroscopy using a plasmonic laser created from a periodically-perforated silver film with a liquid gain medium. An active zone of the laser is formed by a highly elongated spot of pumping. This results in an 80x decrease in the threshold pumping level; in a significantly more efficient diffusive mixing of dye molecules, which substantially suppresses the effect of their bleaching; and in the ability to reduce the volume of the gain medium to as little as 400 nL. We use this design for a stable plasmonic laser in multiple measurements and demonstrate that it is highly effective as a spaser spectroscopy sensor for intracavity detection of an absorptive dye at 70 ppb. This work provides an opportunity to develop applications of intracavity plasmonic laser spectroscopy in biological label detection and other fields.


title: On regularity theory for n/p-harmonic maps into manifolds

abstract: In this paper we continue the investigation of the regularity of the so-called weak $\frac{n}{p}$-harmonic maps in the critical case. These are critical points of the following nonlocal energy \[ {\mathcal{L}}_s(u)=\int_{\mathbb{R}^n}| ( {-\Delta})^{\frac{s}{2}} u(x)|^p dx\,, \] where $u\in \dot{H}^{s,p}(\mathbb{R}^n,\mathcal{N})$ and ${\mathcal{N}}\subset\mathbb{R}^N$ is a closed $k$ dimensional smooth manifold and $s=\frac{n}{p}$. We prove H\"older continuity for such critical points for $p \leq 2$. For $p > 2$ we obtain the same under an additional Lorentz-space assumption. The regularity theory is in the two cases based on regularity results for nonlocal Schr\"odinger systems with an antisymmetric potential.


title: Singularities of Whitham flows for hyperelliptic spectral curves

abstract: We consider the Whitham equations for deformations of hyperelliptic spectral curves, which preserve all periods of a meromorphic differential. If the meromorphic differential has a root at a fixed point of the hyperelliptic involution, then the Whitham flow has a singularity. We prove that the stable and unstable manifolds are non-empty and extend the Whitham flow continuously through the singularity.


title: A hybrid six-dimensional muon cooling channel using gas filled rf cavities

abstract: An alternative cooling approach to prevent rf breakdown in magnetic fields is described that simultaneously reduces all six phase-space dimensions of a muon beam. In this process, cooling is accomplished by reducing the beam momentum through ionization energy loss in discrete absorbers and replenishing the momentum loss only in the longitudinal direction through gas-filled rf cavities. The advantage of gas filled cavities is that they can run at high gradients in magnetic fields without breakdown. With this approach, we show that our channel can achieve a decrease of the 6-dimensional phase-space volume by several orders of magnitude. With the aid of numerical simulations, we demonstrate that the transmission of our proposed channel is comparable to that of an equivalent channel with vacuum rf cavities. Finally, we discuss the sensitivity of the channel performance to the choice of gas and operating pressure.


title: The neutral Higgs self-couplings in the (h)MSSM

abstract: We consider the Minimal Supersymmetric extension of the Standard Model in the regime where the supersymmetric breaking scale is extermely large. In this MSSM, not only the Higgs masses will be affected by large radiative corrections, the dominant part of which is provided by the third generation quark/squark sector, but also the various self-couplings among the Higgs states. In this note, assuming that squarks are extremely heavy, we evaluate the next-to-leading order radiative corrections to the two neutral CP-even Higgs self-couplings $\lambda_{H h h}$ and $\lambda_{h h h}$ and to the partial decay width $\Gamma(H \to h h)$ that are most relevant at the LHC. The calculation is performed using an effective field theory approach that resums the large logarithmic squark contributions and allows to keep under control the perturbative expansion. Since the direct loop vertex corrections are generally missing in this effective approach, we have properly renormalised the effective theory to take them into account. Finally, we perform a comparison of the results in this effective MSSM with those obtained in a much simpler way in the so-called hMSSM approach in which the mass value for the lightest Higgs boson $M_{h}=125$ GeV is used as an input. We show that the hMSSM provides a reasonably good approximation of the corrected self-couplings and $H \to hh$ decay rate and, hence, it can be used also in these cases.


title: Partial preservation of frequencies and Floquet exponents of invariant tori in the reversible KAM context 2

abstract: We consider the persistence of smooth families of invariant tori in the reversible context 2 of KAM theory under various weak nondegeneracy conditions via Herman's method. The reversible KAM context 2 refers to the situation where the dimension of the fixed point manifold of the reversing involution is less than half the codimension of the invariant torus in question. The nondegeneracy conditions we employ ensure the preservation of any prescribed subsets of the frequencies of the unperturbed tori and of their Floquet exponents (the eigenvalues of the coefficient matrix of the variational equation along the torus).


title: Nearest Embedded and Embedding Self-Nested Trees

abstract: Self-nested trees present a systematic form of redundancy in their subtrees and thus achieve optimal compression rates by DAG compression. A method for quantifying the degree of self-similarity of plants through self-nested trees has been introduced by Godin and Ferraro in 2010. The procedure consists in computing a self-nested approximation, called the nearest embedding self-nested tree, that both embeds the plant and is the closest to it. In this paper, we propose a new algorithm that computes the nearest embedding self-nested tree with a smaller overall complexity, but also the nearest embedded self-nested tree. We show from simulations that the latter is mostly the closest to the initial data, which suggests that this better approximation should be used as a privileged measure of the degree of self-similarity of plants.


title: Optimal velocity control of a viscous Cahn-Hilliard system with convection and dynamic boundary conditions

abstract: In this paper, we investigate a distributed optimal control problem for a convective viscous Cahn-Hilliard system with dynamic boundary conditions. Such systems govern phase separation processes between two phases taking place in an incompressible fluid in a container and, at the same time, on the container boundary. The cost functional is of standard tracking type, while the control is exerted by the velocity of the fluid in the bulk. In this way, the coupling between the state (given by the associated order parameter and chemical potential) and control variables in the governing system of nonlinear partial differential equations is bilinear, which presents an additional difficulty for the analysis. The nonlinearities in the bulk and surface free energies are of logarithmic type, which entails that the thermodynamic forces driving the phase separation process may become singular. We show existence for the optimal control problem under investigation, prove the Fr\'echet differentiability of the associated control-to-state mapping in suitable Banach spaces and derive the first-order necessary optimality conditions in terms of a variational inequality and the associated adjoint system. Due to the strong nonlinear couplings between state variables and control, the corresponding proofs require a considerable analytical effort.


title: Sufficiently collapsed irreducible Alexandrov 3-spaces are geometric

abstract: We prove that sufficiently collapsed, closed and irreducible three-dimensional Alexandrov spaces are modeled on one of the eight three-dimensional Thurston geometries. This extends a result of Shioya and Yamaguchi, originally formulated for Riemannian manifolds, to the Alexandrov setting.


title: Morse theory methods for a class of quasi-linear elliptic systems of higher order

abstract: We develop the local Morse theory for a class of non-twice continuously differentiable functionals on Hilbert spaces, including a new generalization of the Gromoll-Meyer's splitting theorem and a weaker Marino-Prodi perturbation type result. They are applicable to a wide range of multiple integrals with quasi-linear elliptic Euler equations and systems of higher order.


title: Tunable $\varphi$-Josephson junction with a quantum anomalous Hall insulator

abstract: We theoretically study the Josephson current in a superconductor/quantum anomalous Hall insulator/superconductor junction by using the lattice Green function technique. When an in-plane external Zeeman field is applied to the quantum anomalous Hall insulator, the Josephson current $J$ flows without a phase difference across the junction $\theta$. The phase shift $\varphi$ appealing in the current-phase relationship $J\propto \sin(\theta-\varphi$) is proportional to the amplitude of Zeeman fields and depends on the direction of Zeeman fields. A phenomenological analysis of the Andreev reflection processes explains the physical origin of $\varphi$. A quantum anomalous Hall insulator breaks time-reversal symmetry and mirror reflection symmetry simultaneously. However it preserves magnetic mirror reflection symmetry. Such characteristic symmetry property enable us to have a tunable $\varphi$-junction with a quantum Hall insulator.


title: When Labels Fall Short: Property Graph Simulation via Blending of Network Structure and Vertex Attributes

abstract: Property graphs can be used to represent heterogeneous networks with labeled (attributed) vertices and edges. Given a property graph, simulating another graph with same or greater size with the same statistical properties with respect to the labels and connectivity is critical for privacy preservation and benchmarking purposes. In this work we tackle the problem of capturing the statistical dependence of the edge connectivity on the vertex labels and using the same distribution to regenerate property graphs of the same or expanded size in a scalable manner. However, accurate simulation becomes a challenge when the attributes do not completely explain the network structure. We propose the Property Graph Model (PGM) approach that uses a label augmentation strategy to mitigate the problem and preserve the vertex label and the edge connectivity distributions as well as their correlation, while also replicating the degree distribution. Our proposed algorithm is scalable with a linear complexity in the number of edges in the target graph. We illustrate the efficacy of the PGM approach in regenerating and expanding the datasets by leveraging two distinct illustrations. Our open-source implementation is available on GitHub.


title: Secure Full-Duplex Device-to-Device Communication

abstract: This paper considers full-duplex (FD) device-to-device (D2D) communications in a downlink MISO cellular system in the presence of multiple eavesdroppers. The D2D pair communicate sharing the same frequency band allocated to the cellular users (CUs). Since the D2D users share the same frequency as the CUs, both the base station (BS) and D2D transmissions interfere each other. In addition, due to limited processing capability, D2D users are susceptible to external attacks. Our aim is to design optimal beamforming and power control mechanism to guarantee secure communication while delivering the required quality-of-service (QoS) for the D2D link. In order to improve security, artificial noise (AN) is transmitted by the BS. We design robust beamforming for secure message as well as the AN in the worst-case sense for minimizing total transmit power with imperfect channel state information (CSI) of all links available at the BS. The problem is strictly non-convex with infinitely many constraints. By discovering the hidden convexity of the problem, we derive a rank-one optimal solution for the power minimization problem.


title: Deep Active Inference

abstract: This work combines the free energy principle from cognitive neuroscience and the ensuing active inference dynamics with recent advances in variational inference in deep generative models, and evolution strategies to introduce the "deep active inference" agent. This agent minimises a variational free energy bound on the average surprise of its sensations, which is motivated by a homeostatic argument. It does so by optimising the parameters of a generative latent variable model of its sensory inputs, together with a variational density approximating the posterior distribution over the latent variables, given its observations, and by acting on its environment to actively sample input that is likely under this generative model. The internal dynamics of the agent are implemented using deep and recurrent neural networks, as used in machine learning, making the deep active inference agent a scalable and very flexible class of active inference agent. Using the mountain car problem, we show how goal directed behaviour can be implemented by defining appropriate priors on the latent states in the agent's model. Furthermore, we show that the deep active inference agent can learn a generative model of the environment, which can be sampled from to understand the agent's beliefs about the environment and its interaction therewith.


title: Exploring RNA structure and dynamics through enhanced sampling simulations

abstract: RNA function is intimately related to its structural dynamics. Molecular dynamics simulations are useful for exploring biomolecular flexibility but are severely limited by the accessible timescale. Enhanced sampling methods allow this timescale to be effectively extended in order to probe biologically-relevant conformational changes and chemical reactions. Here, we review the role of enhanced sampling techniques in the study of RNA systems. We discuss the challenges and promises associated with the application of these methods to force-field validation, exploration of conformational landscapes and ion/ligand-RNA interactions, as well as catalytic pathways. Important technical aspects of these methods, such as the choice of the biased collective variables and the analysis of multi-replica simulations, are examined in detail. Finally, a perspective on the role of these methods in the characterization of RNA dynamics is provided.

categories: q-bio.BM physics.bio-ph physics.chem-ph physics.comp-ph


title: TIPS: Mining Top-K Locations to Minimize User-Inconvenience for Trajectory-Aware Services

abstract: Facility location problems aim to identify the best locations to set up new services. Majority of the existing works typically assume that the users are static. However, there exists a wide array of services such as fuel stations, ATMs, food joints, etc., that are widely accessed by mobile users besides the static ones. Such trajectory-aware services should, therefore, factor in the trajectories of its users rather than simply their static locations. In this work, we introduce the problem of optimal placement of facility locations for such trajectory-aware services that minimize the user inconvenience. The inconvenience of a user is the extra distance traveled by her from her regular path to avail a service. We call this the TIPS problem (Trajectory-aware Inconvenience-minimizing Placement of Services) and consider two variants of it. The goal of the first variant, MAXTIPS, is to minimize the maximum inconvenience faced by any user, while that of the second, AVGTIPS, is to minimize the average inconvenience over all the users. We show that both these problems are NP-hard, and propose multiple efficient heuristics to solve them. Empirical evaluation on real urban-scale road networks validate the efficiency and effectiveness of the proposed heuristics.


title: Tunable Optoelectronic Properties of Triply-Bonded Carbon Molecules with Linear and Graphyne Substructures

abstract: In this paper we present a detailed computational study of the electronic structure and optical properties of triply-bonded hydrocarbons with linear, and graphyne substructures, with the aim of identifying their potential in opto-electronic device applications. For the purpose, we employed a correlated electron methodology based upon the Pariser-Parr-Pople model Hamiltonian, coupled with the configuration interaction (CI) approach, and studied structures containing up to 42 carbon atoms. Our calculations, based upon large-scale CI expansions, reveal that the linear structures have intense optical absorption at the HOMO-LUMO gap, while the graphyne ones have those at higher energies. Thus, the opto-electronic properties depend on the topology of the {graphyne substructures, suggesting that they can be tuned by means of structural modifications. Our results are in very good agreement with the available experimental data.


title: Higher Spin Chern-Simons Theory and the Super Boussinesq hierarchy

abstract: In this paper we construct a map between a solution of supersymmetric Chern-Simons higher spin gravity based on the superalgebra $sl(3|2)$ with Lifshitz scaling and the $N=2$ super Boussinesq hierarchy. We show that under this map the time evolution equations of both theories coincide. In addition, we identify the Poisson structure of the Chern-Simons theory induced by gauge transformation with the second Hamiltonian structure of the super Boussinesq hierarchy.


title: Towards Runtime Adaptation of Actor Systems

abstract: In this dissertation we focus on providing effective adaptations that can be localised and applied to specific concurrent actors, thereby only causing a temporary disruption to the parts of the system requiring mitigation, while leaving the rest of the system intact. We make the application of localised adaptations efficient through incremental synchronisation, whereby the specifier can strategically suspend specific parts of the system, whenever this is strictly required for ensuring that adaptations are effectively applied. We also study static analysis techniques to determine whether the specified incremental synchronisation is in some sense adequate for local adaptations to be carried out. We thus identify a number of generic adaptations that can be applied to any actor system, regardless of its design and the code that it executes. We implement the identified adaptations as an extension of an existing Runtime Verification tool for actor-systems, thereby creating a RA framework for monitoring and mitigating actor systems. In parallel to our implementation we also develop a formal model of our RA framework that further serves to guide our implementation. This model also enables us to better understand the subtle errors that erroneously specified adaptation scripts may introduce. We thus develop a static type system for detecting and rejecting erroneous adaptation scripts prior to deployment, thereby providing the specifier with assistance for writing valid scripts. Although the static typesystem analyses scripts with respect to certain assumptions, we do not assume that the monitored system abides by these assumptions. We therefore augment our RA framework with dynamic checks for halting monitoring whenever the system deviates from our assumption. Based on this dynamically checked model of our RA framework, we prove type soundness for our static type system.


title: Local well-posedness of the Hall-MHD system in $H^s(\mathbb {R}^n)$ with $s>\frac n2$

abstract: We establish local well-posedness of the Hall-magneto-hydrodynamics (Hall-MHD) system in the Sobolev space $\left(H^s(\mathbb{R}^n)\right)^2$ with $s>\frac n2$. The previously known local well-posedness space was $\left(H^s(\mathbb{R}^n)\right)^2$ with $s>\frac n2+1$. Thus the result presented here is an improvement.


title: Ping-pong configurations and circular orders on free groups

abstract: We discuss actions of free groups on the circle with "ping-pong" dynamics; these are dynamics determined by a finite amount of combinatorial data, analogous to Schottky domains or Markov partitions. Using this, we show that the free group $F_n$ admits an isolated circular order if and only if n is even, in stark contrast with the case for linear orders. This answers a question from (Mann, Rivas, 2016). Inspired by work of Alvarez, Barrientos, Filimonov, Kleptsyn, Malicet, Menino and Triestino, we also exhibit examples of "exotic" isolated points in the space of all circular orders on $F_2$. Analogous results are obtained for linear orders on the groups $F_n \times \mathbb{Z}$.


title: A Deep Reinforcement Learning Chatbot

abstract: We present MILABOT: a deep reinforcement learning chatbot developed by the Montreal Institute for Learning Algorithms (MILA) for the Amazon Alexa Prize competition. MILABOT is capable of conversing with humans on popular small talk topics through both speech and text. The system consists of an ensemble of natural language generation and retrieval models, including template-based models, bag-of-words models, sequence-to-sequence neural network and latent variable neural network models. By applying reinforcement learning to crowdsourced data and real-world user interactions, the system has been trained to select an appropriate response from the models in its ensemble. The system has been evaluated through A/B testing with real-world users, where it performed significantly better than many competing systems. Due to its machine learning architecture, the system is likely to improve with additional data.


title: Interaction anisotropy and the KPZ to KPZQ transition in particle deposition at the edges of drying drops

abstract: The deposition process at the edge of evaporating colloidal drops varies with the shape of suspended particles. Experiments with prolate ellipsoidal particles suggest that the spatiotemporal properties of the deposit depend strongly on particle aspect ratio. As the aspect ratio increases, the particles form less densely-packed deposits and the statistical behavior of the deposit interface crosses over from the Kardar-Parisi-Zhang (KPZ) universality class to another universality class which was suggested to be consistent with the KPZ plus quenched disorder. Here, we numerically study the effect of particle interaction anisotropy on deposit growth. In essence, we model the ellipsoids, at the interface, as disk-like particles with two types of interaction patches that correspond to specific features at the poles and equator of the ellipsoid. The numerical results corroborate experimental observations and further suggest that the deposition transition can stem from interparticle interaction anisotropy. Possible extensions of our model to other systems are also discussed.


title: Fourier Method for Approximating Eigenvalues of Indefinite Stekloff Operator

abstract: We introduce an efficient method for computing the Stekloff eigenvalues associated with the Helmholtz equation. In general, this eigenvalue problem requires solving the Helmholtz equation with Dirichlet and/or Neumann boundary condition repeatedly. We propose solving the related constant coefficient Helmholtz equation with Fast Fourier Transform (FFT) based on carefully designed extensions and restrictions of the equation. The proposed Fourier method, combined with proper eigensolver, results in an efficient and clear approach for computing the Stekloff eigenvalues.


title: From large deviations to semidistances of transport and mixing: coherence analysis for finite Lagrangian data

abstract: One way to analyze complicated non-autonomous flows is through trying to understand their transport behavior. In a quantitative, set-oriented approach to transport and mixing, finite time coherent sets play an important role. These are time-parametrized families of sets with unlikely transport to and from their surroundings under small or vanishing random perturbations of the dynamics. Here we propose, as a measure of transport and mixing for purely advective (i.e., deterministic) flows, (semi)distances that arise under vanishing perturbations in the sense of large deviations. Analogously, for given finite Lagrangian trajectory data we derive a discrete-time and space semidistance that comes from the "best" approximation of the randomly perturbed process conditioned on this limited information of the deterministic flow. It can be computed as shortest path in a graph with time-dependent weights. Furthermore, we argue that coherent sets are regions of maximal farness in terms of transport and mixing, hence they occur as extremal regions on a spanning structure of the state space under this semidistance---in fact, under any distance measure arising from the physical notion of transport. Based on this notion we develop a tool to analyze the state space (or the finite trajectory data at hand) and identify coherent regions. We validate our approach on idealized prototypical examples and well-studied standard cases.


title: Boltzmann-type models with uncertain binary interactions

abstract: In this paper we study binary interaction schemes with uncertain parameters for a general class of Boltzmann-type equations with applications in classical gas and aggregation dynamics. We consider deterministic (i.e., a priori averaged) and stochastic kinetic models, corresponding to different ways of understanding the role of uncertainty in the system dynamics, and compare some thermodynamic quantities of interest, such as the mean and the energy, which characterise the asymptotic trends. Furthermore, via suitable scaling techniques we derive the corresponding deterministic and stochastic Fokker-Planck equations in order to gain more detailed insights into the respective asymptotic distributions. We also provide numerical evidences of the trends estimated theoretically by resorting to recently introduced structure preserving uncertainty quantification methods.


title: Conformal Window 2.0: The Large $N_f$ Safe Story

abstract: We extend the phase diagram of SU(N) gauge-fermion theories as function of number of flavours and colours to the region in which asymptotic freedom is lost. We argue, using large $N_f$ results, for the existence of an ultraviolet interacting fixed point at sufficiently large number of flavours opening up to a second ultraviolet conformal window in the number of flavours vs colours phase diagram. We first review the state-of-the-art for the large $N_f$ beta function and then estimate the lower boundary of the ultraviolet window. The theories belonging to this new region are examples of safe non-abelian quantum electro dynamics, termed here {\it safe QCD}. Therefore, according to Wilson, they are fundamental. An important critical quantity is the fermion mass anomalous dimension at the ultraviolet fixed point that we determine at leading order in $1/N_f$. We discover that its value is comfortably below the bootstrap bound. We also investigate the abelian case and find that at the potential ultraviolet fixed point the related fermion mass anomalous dimension has a singular behaviour suggesting that a more careful investigation of its ultimate fate is needed.


title: Quantum Computation of Scattering Amplitudes in Scalar Quantum Electrodynamics

abstract: We present a quantum algorithm for the calculation of scattering amplitudes of massive charged scalar particles in scalar quantum electrodynamics. Our algorithm is based on continuous-variable quantum computing architecture resulting in expo- nential speedup over classical methods. We derive a simple form of the Hamiltonian including interactions, and a straightforward implementation of the constraint due to gauge invariance.


title: Constraining the dark energy models with H(z) data: an approach independent of $H_{0}$

abstract: We study the performance of the latest $H(z)$ data in constraining the cosmological parameters of different cosmological models, including that of Chevalier-Polarski-Linder $w_{0}w_{1}$ parametrization. First, we introduce a statistical procedure in which the chi-square estimator is not affected by the value of the Hubble constant. As a result, we find that the $H(z)$ data do not rule out the possibility of either non-flat models or dynamical dark energy cosmological models. However, we verify that the time varying equation of state parameter $w(z)$ is not constrained by the current expansion data. Combining the $H(z)$ and the Type Ia supernova data we find that the $H(z)$/SNIa overall statistical analysis provides a substantial improvement of the cosmological constraints with respect to those of the $H(z)$ analysis. Moreover, the $w_{0}-w_{1}$ parameter space provided by the $H(z)$/SNIa joint analysis is in a very good agreement with that of Planck 2015, which confirms that the present analysis with the $H(z)$ and SNIa probes correctly reveals the expansion of the Universe as found by the team of Planck. Finally, we generate sets of Monte Carlo realizations in order to quantify the ability of the $H(z)$ data to provide strong constraints on the dark energy model parameters. The Monte Carlo approach shows significant improvement of the constraints, when increasing the sample to 100 $H(z)$ measurements. Such a goal can be achieved in the future, especially in the light of the next generation of surveys.


title: Learning from lions: inferring the utility of agents from their trajectories

abstract: We build a model using Gaussian processes to infer a spatio-temporal vector field from observed agent trajectories. Significant landmarks or influence points in agent surroundings are jointly derived through vector calculus operations that indicate presence of sources and sinks. We evaluate these influence points by using the Kullback-Leibler divergence between the posterior and prior Laplacian of the inferred spatio-temporal vector field. Through locating significant features that influence trajectories, our model aims to give greater insight into underlying causal utility functions that determine agent decision-making. A key feature of our model is that it infers a joint Gaussian process over the observed trajectories, the time-varying vector field of utility and canonical vector calculus operators. We apply our model to both synthetic data and lion GPS data collected at the Bubye Valley Conservancy in southern Zimbabwe.


title: Partial Differential Chow Forms and a Type of Partial Differential Chow varieties

abstract: We first present an intersection theory of partial differential varieties with quasi-generic differential hypersurfaces. Then based on the generic intersection theory, we define the partial differential Chow form for an irreducible partial differential variety $V$ of Kolchin polynomial $\omega_V(t)=(d+1){t+m\choose m}-{t+m-s\choose m}$. And we establish for the partial differential Chow form most of the basic properties of the ordinary differential Chow form. Furthermore, we prove the existence of a type of partial differential Chow varieties.


title: Stopping Times of Random Walks on a Hypercube

abstract: A random walk on a $N$-dimensional hypercube is a discrete time stochastic process whose state space is the set $\{-1,+1\}^{N}$, which has uniform probability of reaching any neighbour state, and probability zero of reaching a non-neighbour state, in one step. This random walk is often studied as a process associated with the Ehrenfest Urn Model. This paper aims to present results about the time that such random walk takes to self-intersect and to return to a set of states. We also present results about the time that the random walk on a hypercube takes to visit a given set and a random set of states. Asymptotic distributions and bounds are presented for these times. The coupling of random walks is widely used as a tool to prove the results.


title: Deterministic physical systems under uncertain initial conditions: the case of maximum entropy applied to projectile motion

abstract: The kinematics and dynamics of deterministic physical systems have been a foundation of our understanding of the world since Galileo and Newton. For real systems, however, uncertainty is largely present via external forces such as friction or lack of precise knowledge about the initial conditions of the system. In this work we focus in the latter case and describe the use of inference methodologies in solving the statistical properties of classical systems subject to uncertain initial conditions. In particular we describe the application of the formalism of Maximum Entropy (MaxEnt) inference to the problem of projectile motion given information about the average horizontal range over many realizations. By using MaxEnt we can invert the problem and use the provided information on the average range to reduce the original uncertainty in the initial conditions, while also achieving additional insights based on the shape of the posterior probabilities for the initial conditions probabilities and the projectile path distribution itself. The wide applicability of this procedure, as well as its ease of use, reveals a useful tool by which to revisit a large number of physics problems, from classrooms to frontier research.


title: Transversal magnetoresistance and Shubnikov-de Haas oscillations in Weyl semimetals

abstract: We explore theoretically the magnetoresistance of Weyl semimetals in transversal magnetic fields away from charge neutrality. The analysis within the self-consistent Born approximation is done for the two different models of disorder: (i) short-range impurties and (ii) charged (Coulomb) impurities. For these models of disorder, we calculate the conductivity away from charge neutrality point as well as the Hall conductivity, and analyze the transversal magnetoresistance (TMR) and Shubnikov-de Haas oscillations for both types of disorder. We further consider a model with Weyl nodes shifted in energy with respect to each other (as found in various materials) with the chemical potential corresponding to the total charge neutrality. In the experimentally most relevant case of Coulomb impurities, we find in this model a large TMR in a broad range of quantizing magnetic fields. More specifically, in the ultra-quantum limit, where only the zeroth Landau level is effective, the TMR is linear in magnetic field. In the regime of moderate (but still quantizing) magnetic fields, where the higher Landau levels are relevant, the rapidly growing TMR is supplemented by strong Shubnikov-de Haas oscillations, consistent with experimental observations.


title: Random Coin Tossing with unknown bias

abstract: Consider a coin tossing experiment which consists of tossing one of two coins at a time, according to a renewal process. The first coin is fair and the second has probability $1/2 + \theta$, $\theta \in [-1/2,1/2]$, $\theta$ unknown but fixed, of head. The biased coin is tossed at the renewal times of the process, and the fair one at all the other times. The main question about this experiment is whether or not it is possible to determine $\theta$ almost surely as the number of tosses increases, given only the probabilities of the renewal process and the observed sequence of heads and tails. We will construct a confidence interval for $\theta$ and determine conditions on the process for its almost sure convergence. It will be shown that recurrence is in fact a necessary condition for the almost sure convergence of the interval, although the convergence still holds if the process is null recurrent but the expected number of renewals up to and including time $N$ is $O(N^{1/2+\alpha}), 0 \leq \alpha < 1/2$. It solves an open problem presented by Harris and Keane (1997). We also generalize this experiment for random variables on $L^{2}$ which are sampled according to a renewal process from either one of two distributions.


title: Topologically Charged Nodal Surface

abstract: We report the existence of topologically charged nodal surface, a band degeneracy on a two-dimensional surface in momentum space that is topologically charged. We develop a Hamiltonian for such charged nodal surface, and show that such a Hamiltonian can be implemented in a tight-binding model as well as in an acoustic meta-material. We also identify a topological phase transition, through which the charges of the nodal surface changes by absorbing or emitting an integer number of Weyl points. Our result indicates that in the band theory, topologically charged objects are not restrict to zero dimension as in a Weyl point, and thus pointing to previously unexplored opportunities for the design of topological materials.

categories: physics.optics cond-mat.mes-hall cond-mat.other physics.class-ph


title: Strongly Enhanced Sensitivity in Planar Microwave Sensors Based on Metamaterial coupling

abstract: Limited sensitivity and sensing range are arguably the greatest challenges in microwave sensor design. Recent attempts to improve these properties have relied on metamaterial- (MTM-) inspired open-loop resonators (OLRs) coupled to transmission lines (TLs). Although the strongly resonant properties of the OLR sensitively reflect small changes in the environment through a shift in its resonance frequency, the resulting sensitivities remain ultimately limited by the level of coupling between the OLR and the TL. This work introduces a novel solution to this problem that employs negative-refractiveindex TL (NRI-TL) MTMs to substantially improve this coupling so as to fully exploit its resonant properties. A MTM-infused planar microwave sensor is designed for operation at 2.5 GHz, and is shown to exhibit a significant improvement in sensitivity and linearity. A rigorous signal-flow analysis (SFA) of the sensor is proposed and shown to provide a fully analytical description of all salient features of both the conventional and MTM-infused sensors. Full-wave simulations confirm the analytical predictions, and all data demonstrate excellent agreement with measurements of a fabricated prototype. The proposed device is shown to be especially useful in the characterization of commonly-available high-permittivity liquids as well as in sensitively distinguishing concentrations of ethanol/methanol in water.


title: The exceptional TeV flaring activity of the blazar 1ES 1959+650 in 2015 and 2016 as observed with VERITAS

abstract: The high-synchrotron-peaked blazar 1ES 1959+650 was among the first extragalactic sources detected in the very high energy gamma ray band (VHE, E > 100 GeV). In October 2015, the source entered an extended period of activity that continued through July 2016, during which several strong VHE flares were observed. This flaring activity in the TeV band was accompanied by a strong increase in the optical, X-ray, and GeV gamma-ray flux of the source, surpassing its brightest recorded flux states. The VERITAS telescope array performed observations of 1ES 1959+650 between October 2015 and June 2016, and detected the source multiple times at a flux higher than the Crab nebula flux in the TeV band, representing the brightest flares of this object since 2002. We here present results from the analysis of 32 hours of VERITAS observations obtained during this period and as well as a contemporaneous multi-wavelength observations in the optical, X-ray, and GeV gamma-ray bands.


title: Six Higgs Doublets Model for Dark Matter

abstract: We consider an extension of the Higgs sector in the standard model (SM) with six Higgs doublets. The gauge couplings are unified without supersymmetry in this model. The lightest of the extra Higgs particles, being stabilized by a discrete symmetry imposed from the outset, presents a plausible candidate for dark matter. For a specific acceptable benchmark point, we show that the model is viable regarding the constraints of relic density, direct detection and invisible Higgs decay. We comment on the mean free path of the dark matter candidate.


title: Cosmic acceleration in dust only Universe via energy-momentum powered gravity

abstract: We propose a modified theory of gravitation constructed by the addition of the term $f(T_{\mu\nu}T^{\mu\nu})$ to the Einstein-Hilbert action, and elaborate a particular case $f(T_{\mu\nu}T^{\mu\nu})=\alpha(T_{\mu\nu}T^{\mu\nu})^{\eta}$, where $\alpha$ and $\eta$ are real constants, dubbed as energy-momentum powered gravity (EMPG). We search for viable cosmologies arising from EMPG especially in the context of the late-time accelerated expansion of the Universe. We investigate the ranges of the EMPG parameters $(\alpha,\eta)$ on theoretical as well as observational grounds leading to the late-time acceleration of the Universe with pressureless matter only, while keeping the successes of standard general relativity at early times. We find that $\eta=0$ corresponds to the $\Lambda$CDM model, whereas $\eta\neq 0$ leads to a $w$CDM-type model. However, the underlying physics of the EMPG model is entirely different in the sense that the energy in the EMPG Universe is sourced by pressureless matter only. Moreover, the energy of the pressureless matter is not conserved, namely, in general it does not dilute as $\rho\propto a^{-3}$ with the expansion of the Universe. Finally, we constrain the parameters of an EMPG-based cosmology with a recent compilation of 28 Hubble parameter measurements, and find that this model describes an evolution of the Universe similar to that in the $\Lambda$CDM model. We briefly discuss that EMPG can be unified with Starobinsky gravity to describe the complete history of the Universe including the inflationary era.


title: Brane world models with bulk perfect fluid and broken 4D Poincare invariance

abstract: We consider 5D brane world models with broken global 4D Poincar\a'{e} invariance (4D part of the spacetime metric is not conformal to the Minkowski spacetime). The bulk is filled with the negative cosmological constant and may contain a perfect fluid. In the case of empty bulk (the perfect fluid is absent), it is shown that one brane solution always has either a physical or a coordinate singularity in the bulk. We cut off these singularities in the case of compact two brane model and obtain regular exact solutions for both 4D Poincar\a'{e} broken and restored invariance. When the perfect fluid is present in the bulk, we get the master equation for the metric coefficients in the case of arbitrary bulk perfect fluid equation of state (EoS) parameters. In two particular cases of EoS, we obtain the analytic solutions for thin and thick branes. First one generalizes the well known Randall-Sundrum model with one brane to the case of the bulk anisotropic perfect fluid. In the second solution, the 4D Poincar\a'{e} invariance is restored. Here, the spacetime goes asymptotically to the anti-de Sitter one far from the thick brane.


title: Brain structural connectivity atrophy in Alzheimer's disease

abstract: Analysis and quantification of brain structural changes, using Magnetic resonance imaging (MRI), are increasingly used to define novel biomarkers of brain pathologies, such as Alzheimer's disease (AD). Network-based models of the brain have shown that both local and global topological properties can reveal patterns of disease propagation. On the other hand, intra-subject descriptions cannot exploit the whole information context, accessible through inter-subject comparisons. To address this, we developed a novel approach, which models brain structural connectivity atrophy with a multiplex network and summarizes it within a classification score. On an independent dataset multiplex networks were able to correctly segregate, from normal controls (NC), AD patients and subjects with mild cognitive impairment that will convert to AD (cMCI) with an accuracy of, respectively, $0.86 \pm 0.01$ and $0.84 \pm 0.01$. The model also shows that illness effects are maximally detected by parceling the brain in equal volumes of $3000$ $mm^3$ ("patches"), without any $a$ $priori$ segmentation based on anatomical features. A direct comparison to standard voxel-based morphometry on the same dataset showed that the multiplex network approach had higher sensitivity. This method is general and can have twofold potential applications: providing a reliable tool for clinical trials and a disease signature of neurodegenerative pathologies.


title: A Nonparametric Statistical Approach to Content Analysis of Items

abstract: In order to use psychometric instruments to assess a multidimensional construct, we may decompose it in dimensions and, in order to assess each dimension, develop a set of items, so one may assess the construct as a whole, by assessing its dimensions. In this scenario, the content analysis of items aims to verify if the developed items are assessing the dimension they are supposed to. In the content analysis process, it is customary to request the judgement of specialists in the studied construct about the dimension that the developed items assess, what makes it a subjective process as it relies upon the personal opinion of the specialists. This paper aims to develop a nonparametric statistical approach to the content analysis of items in order to present a practical method to assess the consistency of the content analysis process, by the development of a statistical test that seeks to determine if all the specialists have the same capability to judge the items. A simulation study is conducted to assess the consistency of the test and it is applied to a real validation process.


title: PWC-Net: CNNs for Optical Flow Using Pyramid, Warping, and Cost Volume

abstract: We present a compact but effective CNN model for optical flow, called PWC-Net. PWC-Net has been designed according to simple and well-established principles: pyramidal processing, warping, and the use of a cost volume. Cast in a learnable feature pyramid, PWC-Net uses the cur- rent optical flow estimate to warp the CNN features of the second image. It then uses the warped features and features of the first image to construct a cost volume, which is processed by a CNN to estimate the optical flow. PWC-Net is 17 times smaller in size and easier to train than the recent FlowNet2 model. Moreover, it outperforms all published optical flow methods on the MPI Sintel final pass and KITTI 2015 benchmarks, running at about 35 fps on Sintel resolution (1024x436) images. Our models are available on https://github.com/NVlabs/PWC-Net.


title: Collisional and thermal dissociation of $J/\psi$ and $\Upsilon$ states at the LHC

abstract: We present new results for the suppression of high transverse momentum charmonium [$J/\psi, \psi(2S)$] and bottomonium [$\Upsilon(1S),\Upsilon(2S),\Upsilon(3S)$] states in Pb+Pb collisions at the Large Hadron Collider. Our theoretical formalism combines the collisional dissociation of quarkonia, as they propagate in the quark-gluon plasma, with the thermal wavefunction effects due to the screening of the $Q\bar{Q}$ attractive potential in the medium. We find that a good description of the relative suppression of the ground and higher excited quarkonium states, transverse momentum and centrality distributions is achieved, when comparison to measurements at a center-of-mass energy of 2.76 TeV is performed. Theoretical predictions for the highest Pb+Pb center-of-mass energy of 5.02 TeV at the LHC, where new experimental results are being finalized, are also presented.


title: Adaptive PCA for Time-Varying Data

abstract: In this paper, we present an online adaptive PCA algorithm that is able to compute the full dimensional eigenspace per new time-step of sequential data. The algorithm is based on a one-step update rule that considers all second order correlations between previous samples and the new time-step. Our algorithm has O(n) complexity per new time-step in its deterministic mode and O(1) complexity per new time-step in its stochastic mode. We test our algorithm on a number of time-varying datasets of different physical phenomena. Explained variance curves indicate that our technique provides an excellent approximation to the original eigenspace computed using standard PCA in batch mode. In addition, our experiments show that the stochastic mode, despite its much lower computational complexity, converges to the same eigenspace computed using the deterministic mode.


title: Nonreciprocal transmission of sound in viscous fluid with asymmetric scatterers

abstract: Two common concepts of nonreciprocity in sound propagation are based on nonlinear effects [1, 2] and on local circulation of fluid [3, 4]. They originate from two known methods of breaking a time reversal symmetry, that is necessary for observation of nonreciprocal effects. Both concepts require additional devices to be installed with their own power sources. Recently it was demonstrated that acoustical losses may serve as a source of T-symmetry violation, thus leading to nonreciprocity in reflection of sound from gradient-index metasurface [5]. Here, we explore viscosity of fluid as a natural factor of T-symmetry breaking. We report experimental observation of the nonreciprocal transmission of ultrasound through a water-submerged phononic crystal consisting of asymmetric rods. Asymmetry, or broken P-symmetry, is the second necessary factor for nonreciprocity. Experimental results are in agreement with numerical simulations based on the Navier-Stokes equation. This passive nonreciprocal linear device is cheap, robust and does not require an energy source.


title: A practical, unitary simulator for non-Markovian complex processes

abstract: Stochastic processes are as ubiquitous throughout the quantitative sciences as they are notorious for being difficult to simulate and predict. In this letter we propose a unitary quantum simulator for discrete-time stochastic processes which requires less internal memory than any classical analogue throughout the simulation. The simulator's internal memory requirements equal those of the best previous quantum models. However, in contrast to previous models it only requires a (small) finite-dimensional Hilbert space. Moreover, since the simulator operates unitarily throughout, it avoids any unnecessary information loss. We provide a stepwise construction for simulators for a large class of stochastic processes hence directly opening the possibility for experimental implementations with current platforms for quantum computation. The results are illustrated for an example process.


title: NMR in an electric field: A bulk probe of the hidden spin and orbital polarizations

abstract: Recent theoretical work has established the presence of hidden spin and orbital textures in non-magnetic materials with inversion symmetry. Here, we propose that these textures can be detected by nuclear magnetic resonance (NMR) measurements carried out in the presence of an electric field. In crystals with hidden polarizations, a uniform electric field produces a staggered magnetic field that points to opposite directions at atomic sites related by spatial inversion. As a result, the NMR resonance peak corresponding to inversion partner nuclei is split into two peaks. The magnitude of the splitting is proportional to the electric field and depends on the orientation of the electric field with respect to the crystallographic axes and the external magnetic field. As a case study, we present a theory of electric-field-induced splitting of NMR peaks for $^{77}$Se, $^{125}$Te and $^{209}$Bi in Bi$_2$Se$_3$ and Bi$_2$Te$_3$. In conducting samples with current densities of $\simeq 10^6\, {\rm A/cm}^2$, the splitting for Bi can reach $100\, {\rm kHz}$, which is comparable to or larger than the intrinsic width of the NMR lines. In order to observe the effect experimentally, the peak splitting must also exceed the linewidth produced by the Oersted field. In Bi$_2$Se$_3$, this requires narrow wires of radius $\lesssim 1\, \mu{\rm m}$. We also discuss other potentially more promising candidate materials, such as SrRuO$_3$ and BaIr$_2$Ge$_2$, whose crystal symmetry enables strategies to suppress the linewidth produced by the Oersted field.


title: Rank-Select Indices Without Tears

abstract: A rank-select index for a sequence $B=(b_1,\ldots,b_n)$ of $n$ bits is a data structure that, if provided with an operation to access $\Theta(\log n)$ arbitrary consecutive bits of $B$ in constant time (thus $B$ is stored outside of the data structure), can compute $\mathit{rank}_B(j)=\sum_{i=1}^j b_i$ for given $j\in\{0,\ldots,n\}$ and $\mathit{select}_B(k)=\min\{j:\mathit{rank}_B(j)\ge k\}$ for given $k\in\{1,\ldots,\sum_{i=1}^n b_i\}$. We describe a new rank-select index that, like previous rank-select indices, occupies $O(n\log\log n/\log n)$ bits and executes $\mathit{rank}$ and $\mathit{select}$ queries in constant time. Its derivation is intended to be particularly easy to follow and largely free of tedious low-level detail, its operations are given by straight-line code, and we show that it can be constructed in $O(n/\log n)$ time.


title: Partial accretion in the propeller stage of low mass X-ray binary Aql X--1

abstract: Aql X--1 is one of the most prolific low mass X-ray binary transients (LMXBTs) showing outbursts almost annually. We present the results of our spectral analyses of RXTE/PCA observations of the 2000 and the 2011 outbursts. We investigate the spectral changes related to the changing disk-magnetosphere interaction modes of Aql X--1. The X-ray light curves of the outbursts of LMXBTs typically show phases of fast rise and exponential decay. The decay phase shows a "knee" where the flux goes from the slow decay to the rapid decay stage. We assume that the rapid decay corresponds to a weak propeller stage at which a fraction of the inflowing matter in the disk accretes onto the star. We introduce a novel method for inferring, from the light curve, the fraction of the inflowing matter in the disk that accretes onto the NS depending on the fastness parameter. We determine the fastness parameter range within which the transition from the accretion to the partial propeller stage is realized. This fastness parameter range is a measure of the scale-height of the disk in units of the inner disk radius. We applied the method to a sample of outbursts of Aql X--1 with different maximum flux and duration times. We show that different outbursts with different maximum luminosity and duration follow a similar path in the parameter space of accreted/inflowing mass flux fraction versus fastness parameter.


title: Optimizing collective fieldtaxis of swarming agents through reinforcement learning

abstract: Swarming of animal groups enthralls scientists in fields ranging from biology to physics to engineering. Complex swarming patterns often arise from simple interactions between individuals to the benefit of the collective whole. The existence and success of swarming, however, nontrivially depend on microscopic parameters governing the interactions. Here we show that a machine-learning technique can be employed to tune these underlying parameters and optimize the resulting performance. As a concrete example, we take an active matter model inspired by schools of golden shiners, which collectively conduct phototaxis. The problem of optimizing the phototaxis capability is then mapped to that of maximizing benefits in a continuum-armed bandit game. The latter problem accepts a simple reinforcement-learning algorithm, which can tune the continuous parameters of the model. This result suggests the utility of machine-learning methodology in swarm-robotics applications.


title: Independent Experimentation on the Activation of Deuterium-Loaded Materials by X-Ray Exposure

abstract: In an attempt to replicate work performed elsewhere [1], we have searched for x-ray induced nuclear activation of deuterated materials. Our first results, reported in September 2015, showed no evidence of nuclear activation, contrary to the results reported in [1]. The primary shortcoming of our first attempt, however, was that the x-ray tube we used was capable of a maximum accelerating potential of only 160kV, less than the 200kV used on the comparison test sample from [1]. Thus, our results could not exclude processes initiated by x-rays near 200 keV in energy. A second potential shortcoming was that our x-ray tube did not have a microfocus beam, like the one used in [1]. Recently, we have irradiated additional samples with an x-ray tube matching the specifications of the tube used in [1]. By matching the test conditions with essentially identical equipment, we have overcome both of the shortcomings of our previous measurements. As before, we find no evidence of activation. The bulk of this report, pages 2 through 35, contain the results and analysis from the first round of testing. The extension to 200kV is found in the final appendix, namely Appendix A5.


title: Read Mapping Near Non-Volatile Memory

abstract: DNA sequencing is the physical/biochemical process of identifying the location of the four bases (Adenine, Guanine, Cytosine, Thymine) in a DNA strand. As semiconductor technology revolutionized computing, modern DNA sequencing technology (termed Next Generation Sequencing, NGS)revolutionized genomic research. As a result, modern NGS platforms can sequence hundreds of millions of short DNA fragments in parallel. The sequenced DNA fragments, representing the output of NGS platforms, are termed reads. Besides genomic variations, NGS imperfections induce noise in reads. Mapping each read to (the most similar portion of) a reference genome of the same species, i.e., read mapping, is a common critical first step in a diverse set of emerging bioinformatics applications. Mapping represents a search-heavy memory-intensive similarity matching problem, therefore, can greatly benefit from near-memory processing. Intuition suggests using fast associative search enabled by Ternary Content Addressable Memory (TCAM) by construction. However, the excessive energy consumption and lack of support for similarity matching (under NGS and genomic variation induced noise) renders direct application of TCAM infeasible, irrespective of volatility, where only non-volatile TCAM can accommodate the large memory footprint in an area-efficient way. This paper introduces GeNVoM, a scalable, energy-efficient and high-throughput solution. Instead of optimizing an algorithm developed for general-purpose computers or GPUs, GeNVoM rethinks the algorithm and non-volatile TCAM-based accelerator design together from the ground up. Thereby GeNVoM can improve the throughput by up to 113.5 times (3.6); the energy consumption, by up to 210.9 times (1.36), when compared to a GPU (accelerator) baseline, which represents one of the highest-throughput implementations known.


title: An Approach to Differential Invariants of $G$-Structures

abstract: In a natural way, the local diffeomorphisms of a manifold onto itself act on the reference frame bundles of any order and on the bundles associated with them. Due to the transitivity, the invariants by diffeomorphisms of an associated bundle correspond to the real functions on the orbit space of the action of the jet group on the typical fiber, which enables to construct the said associated bundle. Scalar differential invariants are functions on bundles of jets of sections of associated bundles, which are invariant by the standard action of local diffeomorphisms. We show that these jet bundles are, in turn, associated bundles with reference frame bundles of higher order. By the previous result, this will allow to recognize scalar differential invariants as real functions over a quotient space, regardless of the base manifold. We apply the above to $G$-structures, and describe their scalar differential invariants as functions over certain quotient spaces. Each quotient space contains a dense differentiable manifold, whose dimension is "the number of functionally independent scalar differential invariants". We will obtain a lower bound of that number and calculate it for some cases, in particular for the scalar conformal differential invariants.


title: Four-point probe measurements using current probes with voltage feedback to measure electric potentials

abstract: We present a four-point probe resistance measurement technique which uses four equivalent current measuring units, resulting in minimal hardware requirements and corresponding sources of noise. Local sample potentials are measured by a software feedback loop which adjusts the corresponding tip voltage such that no current flows to the sample. The resulting tip voltage is then equivalent to the sample potential at the tip position. We implement this measurement method into a multi-tip scanning tunneling microscope setup such that potentials can also be measured in tunneling contact, allowing in principle truly non-invasive four-probe measurements. The resulting measurement capabilities are demonstrated for BiSbTe$_3$ and Si$(111)-(7\times7)$ samples.


title: Observational constraints on dark matter-dark energy scattering cross section

abstract: In this letter, we report precise and robust observational constraints on dark matter-dark energy scattering cross section, using the latest data from cosmic microwave background (CMB) Planck temperature and polarization, baryon acoustic oscillations (BAO) measurements and weak gravitational lensing data from Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). The scattering scenario consists of a pure momentum exchange between the dark components, and we find $\sigma_d < 10^{-29} \, {\rm cm^2}$ at 95\% CL from the joint analysis (CMB + BAO + CFHTLenS), for typical dark matter particle mass of the order 1-10 ${\rm GeV}/c^2$. We notice that the scattering among the dark components may influence the growth of large scale structure in the Universe, leaving the background cosmology unaltered.


title: Sub-pm/$\mathrm{\mathbf{\sqrt{\rm Hz}}}$ non-reciprocal noise in the LISA backlink fiber

abstract: The future space-based gravitational wave detector Laser Interferometer Space Antenna (LISA) requires bidirectional exchange of light between its two optical benches on board of each of its three satellites. The current baseline foresees a polarization-maintaining single-mode fiber for this backlink connection. Phase changes which are common in both directions do not enter the science measurement, but differential ("non-reciprocal") phase fluctuations directly do and must thus be guaranteed to be small enough. We have built a setup consisting of a Zerodur$^{\rm TM}$ baseplate with fused silica components attached to it using hydroxide-catalysis bonding and demonstrated the reciprocity of a polarization-maintaining single-mode fiber at the 1 pm/$\sqrt{\textrm{Hz}}$ level as is required for LISA. We used balanced detection to reduce the influence of parasitic optical beams on the reciprocity measurement and a fiber length stabilization to avoid nonlinear effects in our phase measurement system (phase meter). For LISA, a different phase meter is planned to be used that does not show this nonlinearity. We corrected the influence of beam angle changes and temperature changes on the reciprocity measurement in post-processing.


title: The large scale impact of offshore wind farm structures on pelagic primary productivity in the southern North Sea

abstract: The increasing demand for renewable energy is projected to result in a 40-fold increase in offshore wind electricity in the European Union by 2030. Despite a great number of local impact studies for selected marine populations, the regional ecosystem impacts of offshore wind farm structures are not yet well assessed nor understood. Our study investigates whether the accumulation of epifauna, dominated by the filter feeder Mytilus edulis (blue mussel), on turbine structures affects pelagic primary productivity and ecosystem functioning in the southern North Sea. We estimate the anthropogenically increased potential distribution based on the current projections of turbine locations and reported patterns of M. edulis settlement. This distribution is integrated through the Modular Coupling System for Shelves and Coasts to state-of-the-art hydrodynamic and ecosystem models. Our simulations reveal non-negligible potential changes in regional annual primary productivity of up to 8% within the offshore wind farm area, and induced maximal increases of the same magnitude in daily productivity also far from the wind farms. Our setup and modular coupling are effective tools for system scale studies of other environmental changes arising from large-scale offshore wind-farming such as ocean physics and distributions of pelagic top predators.


title: Understanding disk galaxy rotation velocities without dark matter contribution--a physical process for MOND?

abstract: An impact model of gravity designed to emulate Newton's law of gravitation is applied to the radial acceleration of disk galaxies. Based on this model (Wilhelm et al. 2013), the rotation velocity curves can be understood without the need to postulate any dark matter contribution. The increased acceleration in the plane of the disk is a consequence of multiple interactions of gravitons (called "quadrupoles" in the original paper) and the subsequent propagation in this plane and not in three-dimensional space. The concept provides a physical process that relates the fit parameter of the acceleration scale defined by McGaugh et al. (2016) to the mean free path length of gravitons in the disks of galaxies. It may also explain the modification of the gravitational interaction at low acceleration levels in MOND (Milgrom 1983, 1994, 2015, 2016). Three examples are discussed in some detail: The spiral galaxies NGC 7814, NGC 6503 and M 33.


title: Non-Abelian Supertubes

abstract: A supertube is a supersymmetric configuration in string theory which occurs when a pair of branes spontaneously polarizes and generates a new dipole charge extended along a closed curve. The dipole charge of a codimension-2 supertube is characterized by the U-duality monodromy as one goes around the supertube. For multiple codimension-2 supertubes, their monodromies do not commute in general. In this paper, we construct a supersymmetric solution of five-dimensional supergravity that describes two supertubes with such non-Abelian monodromies, in a certain perturbative expansion. In supergravity, the monodromies are realized as the multi-valuedness of the scalar fields, while in higher dimensions they correspond to non-geometric duality twists of the internal space. The supertubes in our solution carry NS5 and 5$_2^2$ dipole charges and exhibit the same monodromy structure as the SU(2) Seiberg-Witten geometry. The perturbative solution has $\mathrm{AdS}_2\times S^2$ asymptotics and vanishing four-dimensional angular momentum. We argue that this solution represents a microstate of four-dimensional black holes with a finite horizon and that it provides a clue for the gravity realization of a pure-Higgs branch state in the dual quiver quantum mechanics.


title: Review: Far-Infrared Instrumentation and Technology Development for the Next Decade

abstract: Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.


title: Entanglement Entropy and Decoupling in the Universe

abstract: In the expanding universe, two interacting fields are no longer in thermal contact when the interaction rate becomes smaller than the Hubble expansion rate. After decoupling, two subsystems are usually treated separately in accordance with equilibrium thermodynamics and the thermodynamic entropy gives a fiducial quantity conserved in each subsystem. In this paper, we discuss a correction to this paradigm from quantum entanglement of two coupled fields. The thermodynamic entropy is generalized to the entanglement entropy. We formulate a perturbation theory to derive the entanglement entropy and present Feynman rules in diagrammatic calculations. For specific models to illustrate our formulation, an interacting scalar-scalar system, quantum electrodynamics, and the Yukawa theory are considered. We calculate the entanglement entropy in these models and find a quantum correction to the thermodynamic entropy. The correction is revealed to be important in circumstances of instantaneous decoupling.


title: Chemical enrichment of the planet forming region as probed by accretion

abstract: The chemical conditions in the planet forming regions of protoplanetary discs remain difficult to observe directly. Gas accreting from the disc on to the star provides a way to measure the elemental abundances because even refractory species are in an atomic gaseous form. Here we compare the abundance ratios derived from UV lines probing T Tauri accretion streams to simple models of disc evolution. Although the interpretation of line ratios in terms of abundances is highly uncertain, discs with large cavities in mm images tend to have lower Si emission. Since this can naturally be explained by the suppressed accretion of dust, this suggests that abundance variations are at least partially responsible for the variations seen in the line ratios. Our models of disc evolution due to grain growth, radial drift and the flux of volatile species carried as ices on grain surfaces, give rise to a partial sorting of the atomic species based on the volatility of their dominant molecular carriers. This arises because volatiles are left behind at their snow lines while the grains continue to drift. We show that this reproduces the main features seen in the accretion line ratio data, such as C/N ratios that are a few times solar and the correlation between the Si to volatile ratio with mm-flux. We highlight the fact that developing a more robust linkage between line ratios and abundance ratios and acquiring data for larger samples has the potential to cast considerable light on the chemical history of protoplanetary discs.


title: A Note on the WGC, Effective Field Theory and Clockwork within String Theory

abstract: It has been recently argued that Higgsing of theories with $U(1)^n$ gauge interactions consistent with the Weak Gravity Conjecture (WGC) may lead to effective field theories parametrically violating WGC constraints. The minimal examples typically involve Higgs scalars with a large charge with respect to a $U(1)$ (e.g. charges $(Z,1)$ in $U(1)^2$ with $Z\gg 1$). This type of Higgs multiplets play also a key role in clockwork $U(1)$ theories. We study these issues in the context of heterotic string theory and find that, while indeed there is no new physics at the standard magnetic WGC scale $\Lambda\sim g_{IR} M_P$, the string scale is just slightly above, at a scale $\sim \sqrt{k_{IR}}\Lambda$. Here $k_{IR}$ is the level of the IR $U(1)$ worldsheet current. We show that, unlike the standard magnetic cutoff, this bound is insensitive to subsequent Higgsing. One may argue that this constraint gives rise to no bound at the effective field theory level since $k_0$ is model dependent and in general unknown. However there is an additional constraint to be taken into account, which is that the Higgsing scalars with large charge $Z$ should be part of the string massless spectrum, which becomes an upper bound $k_{IR}\leq k_0^2$, where $k_0$ is the level of the UV currents. Thus, for fixed $k_0$, $Z$ cannot be made parametrically large. The upper bound on the charges $Z$ leads to limitations on the size and structure of hierarchies in an iterated $U(1)$ clockwork mechanism.


title: The Outer Halos of Very Massive Galaxies: BCGs and their DSC in the Magneticum Simulations

abstract: Recent hydrodynamic cosmological simulations cover volumes up to Gpc^3 and resolve halos across a wide range of masses and environments, from massive galaxy clusters down to normal galaxies, while following a large variety of physical processes (star formation, chemical enrichment, AGN feedback) to allow a self-consistent comparison to observations at multiple wavelengths. Using the Magneticum simulations, we investigate the buildup of the diffuse stellar component (DSC) around massive galaxies within group and cluster environments. The DSC in our simulations reproduces the spatial distribution of the observed intracluster light (ICL) as well as its kinematic properties remarkably well. For galaxy clusters and groups we find that, although the DSC in almost all cases shows a clear separation from the brightest cluster galaxy (BCG) with regard to its dynamic state, the radial stellar density distribution in many halos is often characterized by a single Sersic profile, representing both the BCG component and the DSC, very much in agreement with current observational results. Interestingly, even in those halos that clearly show two components in both the dynamics and the spatial distribution of the stellar component, no correlation between them is evident.


title: Rotational strain in Weyl semimetals: A continuum approach

abstract: The coupling of lattice deformations to the low energy electronic excitations of Dirac matter involve novel types of electron--phonon couplings as the celebrated elastic gauge fields first analyzed in graphene. In the continuum low energy approach, lattice deformations coupling to the electronic degrees of freedom are characterized by the (symmetric) strain tensor defined in elasticity theory. We review these couplings in Weyl semimetals and examine the coupling of electronic excitations to the antisymmetric part of the deformation gradient tensor associated to rotational strain. The new couplings, absent in the two dimensional materials, have important physical implications: they give rise to new elastic gauge fields, contribute to the deformation potential, tilt the cones and generate new pseudo--Zeeman couplings.


title: On the Penrose inequality in anti-deSitter space

abstract: For asymptotically flat spacetimes the Penrose inequality gives an initial data test for the weak cosmic censorship hypothesis. We give a formulation of this inequality for asymptotically anti-deSitter (AAdS) spacetimes, and show that the inequality holds for time asymmetric data in spherical symmetry. Our analysis is motivated by the constant-negative-spatial-curvature form of the AdS black hole metric.


title: Resolving API Mentions in Informal Documents

abstract: Developer forums contain opinions and information related to the usage of APIs. API names in forum posts are often not explicitly linked to their official resources. Automatic linking of an API mention to its official resources can be challenging for various reasons, such as, name overloading. We present a technique, ANACE, to automatically resolve API mentions in the textual contents of forum posts. Given a database of APIs, we first detect all words in a forum post that are potential references to an API. We then use a combination of heuristics and machine learning to eliminate false positives and to link true positives to the actual APIs and their resources.


title: Effective potential at three loops

abstract: I present the effective potential at three-loop order for a general renormalizable theory, using the \MSbar renormalization scheme and Landau gauge fixing. As applications and illustrative points of reference, the results are specialized to the supersymmetric Wess-Zumino model and to the Standard Model. In each case, renormalization group scale invariance provides a consistency check. In the Wess-Zumino model, the required vanishing of the minimum vacuum energy yields an additional check. For the Standard Model, I carry out the resummation of Goldstone boson contributions, which provides yet more opportunities for non-trivial checks, and obtain the minimization condition for the Higgs vacuum expectation value at full three-loop order. An infrared divergence due to doubled photon propagators appears in the three-loop Standard Model effective potential, but it does not affect the minimization condition or physical observables and can be eliminated by resummation.


title: Stripe order from the perspective of the Hubbard model

abstract: A microscopic understanding of the strongly correlated physics of the cuprates must account for the translational and rotational symmetry breaking that is present across all cuprate families, commonly in the form of stripes. Here we investigate emergence of stripes in the Hubbard model, a minimal model believed to be relevant to the cuprate superconductors, using determinant quantum Monte Carlo (DQMC) simulations at finite temperatures and density matrix renormalization group (DMRG) ground state calculations. By varying temperature, doping, and model parameters, we characterize the extent of stripes throughout the phase diagram of the Hubbard model. Our results show that including the often neglected next-nearest-neighbor hopping leads to the absence of spin incommensurability upon electron-doping and nearly half-filled stripes upon hole-doping. The similarities of these findings to experimental results on both electron and hole-doped cuprate families support a unified description across a large portion of the cuprate phase diagram.


title: Synthesis of mixed hypermetallic oxide BaOCa$^+$ from laser-cooled reagents in an atom-ion hybrid trap

abstract: Hypermetallic alkaline earth (M) oxides of formula MOM have been studied under plasma conditions that preclude insight into their formation mechanism. We present here the application of emerging techniques in ultracold physics to the synthesis of a mixed hypermetallic oxide, BaOCa$^+$. These methods, augmented by high-level electronic structure calculations, permit detailed investigation of the bonding and structure, as well as the mechanism of its formation via the barrierless reaction of Ca $(^3P_J)$ with BaOCH$_{3}^+$. Further investigations of the reaction kinetics as a function of collision energy over the range 0.005 K to 30 K and of individual Ca fine-structure levels compare favorably with calculations based on long-range capture theory.


title: Mean ergodicity vs weak almost periodicity

abstract: We provide explicit examples of positive and power-bounded operators on $c_0$ and $\ell^\infty$ which are mean ergodic but not weakly almost periodic. As a consequence we prove that a countably order complete Banach lattice on which every positive and power-bounded mean ergodic operator is weakly almost periodic is necessarily a KB-space. This answers several open questions from the literature. Finally, we prove that if $T$ is a positive mean ergodic operator with zero fixed space on an arbitrary Banach lattice, then so is every power of $T$.


title: 2-FAST: Fast and accurate computation of projected two-point functions

abstract: We present the 2-point function from Fast and Accurate Spherical Bessel Transformation (2-FAST) algorithm for a fast and accurate computation of integrals involving one or two spherical Bessel functions. These types of integrals occur when projecting the galaxy power spectrum $P(k)$ onto the configuration space, $\xi_\ell^\nu(r)$, or spherical harmonic space, $C_\ell(\chi,\chi')$. First, we employ the FFTlog transformation of the power spectrum to divide the calculation into $P(k)$-dependent coefficients and $P(k)$-independent integrations of basis functions multiplied by spherical Bessel functions. We find analytical expressions for the latter integrals in terms of special functions, for which recursion provides a fast and accurate evaluation. The algorithm, therefore, circumvents direct integration of highly oscillating spherical Bessel functions.


title: High Coherence Plane Breaking Packaging for Superconducting Qubits

abstract: We demonstrate a pogo pin package for a superconducting quantum processor specifically designed with a nontrivial layout topology (e.g., a center qubit that cannot be accessed from the sides of the chip). Two experiments on two nominally identical superconducting quantum processors in pogo packages, which use commercially available parts and require modest machining tolerances, are performed at low temperature (10 mK) in a dilution refrigerator and both found to behave comparably to processors in standard planar packages with wirebonds where control and readout signals come in from the edges. Single- and two-qubit gate errors are also characterized via randomized benchmarking. More detailed crosstalk measurements indicate levels of crosstalk less than -40 dB at the qubit frequencies, opening the possibility of integration with extensible qubit architectures.


title: Power Network Regulation Benchmark for Switched-Mode Optimal Control

abstract: Power network regulation is presented as a benchmark problem for assessing and developing switched-mode optimal control approaches like mode scheduling, sliding window scheduling and modal design. Power network evolution modeled by the swing equations and coupled with controllable switching components is a nonlinear, high-dimensional problem. The proposed benchmark problem is the 54 generator IEEE 118 Bus Test Case composed of 106 states. Open questions include scalability in state and number of modes of operation, as well as real-time implementation, reliability, hysteresis, and timing constraints. Can the entire North American power network be regulated? Can every transmission line have independent switching control authority?


title: EMD-regression for modelling multi-scale relationships, and application to weather-related cardiovascular mortality

abstract: In a number of environmental studies, relationships between natural processes are often assessed through regression analyses, using time series data. Such data are often multi-scale and non-stationary, leading to a poor accuracy of the resulting regression models and therefore to results with moderate reliability. To deal with this issue, the present paper introduces the EMD-regression methodology consisting in applying the empirical mode decomposition (EMD) algorithm on data series and then using the resulting components in regression models. The proposed methodology presents a number of advantages. First, it accounts of the issues of non-stationarity associated to the data series. Second, this approach acts as a scan for the relationship between a response variable and the predictors at different time scales, providing new insights about this relationship. To illustrate the proposed methodology it is applied to study the relationship between weather and cardiovascular mortality in Montreal, Canada. The results shed new knowledge concerning the studied relationship. For instance, they show that the humidity can cause excess mortality at the monthly time scale, which is a scale not visible in classical models. A comparison is also conducted with state of the art methods which are the generalized additive models and distributed lag models, both widely used in weather-related health studies. The comparison shows that EMD-regression achieves better prediction performances and provides more details than classical models concerning the relationship.