title: Optical injection of a spin current into a zigzag nanoribbon of monolayer MoS2 with antiferromagnetic Kekule distortion

abstract: The Kekule pattern of the (anti)ferromagnetic exchange field on monolayer MoS2 can be induced by proximity to the (111) surface of BiFeO3 on both sides. The three-band tight binding model of the MoS2 layer with Kekule patterned exchange field is applied to describe the heterostructures. The tight binding model is justified by the first principle calculations. The magnetization orientations of the substrates control the pattern of the exchange field, which then switches the band structures of the lowest zigzag edge states between being metallic and insulating. The lowest four zigzag edge bands provide conducting channels with a spin-polarized current. Optical excitation of carriers in these bands generates sizable spin and charge currents, which are theoretically modeled by the perturbation solution of the semiconductor Bloch equation. The injected spin currents have multiple resonant peaks at a few frequencies, which can be switched off by rotating the magnetization orientations of the substrates.


title: Stability of topologically protected edge states in nonlinear quantum walks: Additional bifurcations unique to Floquet systems

abstract: Recently, effects of nonlinearity on topologically nontrivial systems have attracted attention and the stability of topologically protected edge states has been studied for a quantum walk with nonlinear effects, which is akin to time-periodically driven systems (Floquet systems). In the previous work, it has been found that the edge states can be stable attractors or unstable repellers depending on their intrinsic topological property, while the stability is not affected by the strength of nonlinearity. In the present work, we find additional bifurcations at which edge states change from stable attractors to unstable repellers with increasing the strength of nonlinearity in nonlinear quantum walks, for the first time. The new bifurcations are unique to Floquet systems, since we take dynamical properties of Floquet systems into consideration by directly applying the time-evolution operator of the quantum walks to the linear stability analysis. Our results shed new light on nonlinear effects on topological edge states in Floquet systems.

categories: quant-ph cond-mat.other math-ph math.MP nlin.PS


title: Elastic ripening and inhibition of liquid-liquid phase separation

abstract: Phase separation has recently emerged as an important organizational principle in the dense and heterogeneous environment within living cells. Here, we use a synthetic system to show that compressive stresses in a polymer network suppress phase separation of the solvent that swells it. These stresses create a barrier to droplet nucleation that leads to robust, stabilized mixtures well beyond the liquid-liquid phase separation boundary. Network stresses not only alter the stability of mixtures, but they also have a dramatic effect on the ripening of droplets. Gradients in network stresses can drive a new form of ripening, where solute is transported down stiffness gradients. This elastic ripening can be much faster than conventional surface tension driven Ostwald ripening.


title: Common-path generation of stable cylindrical perfect vector vortex beams with arbitrary order

abstract: A Highly flexible and efficient method of generating stable radially and Azimuthally polarized perfect optical vortex beams and all higher-order cylindrical vector vortex beams are proposed and demonstrated. The method is the most convenient implementation of the superposition of two orthogonally circularly polarized optical vortex beams of arbitrary integer topological charges. By simply controlling the relative amplitude and phase between vertical and horizontal polarization component of an input perfect vortex beam on a phase-sensitive Spatial Light Modulator (SLM) and using a common-path interferometry, all arbitrary order Poincar\'e beams are prepared. Calibration curve of relative phase vs grey-scale of the phase function on the SLM is drawn to facilitate the determination of required phase-offset using pre-calculated grey-scale colour-map. Generated beams states of polarization are uniformly distributed throughout the beam cross-section. The interference pattern of the two-beam in a projected linear polarization state also gives the cylindrically symmetric petal beams with control over the number and orientation of petals.


title: Accelerating Reinforcement Learning through GPU Atari Emulation

abstract: We introduce CuLE (CUDA Learning Environment), a CUDA port of the Atari Learning Environment (ALE) which is used for the development of deep reinforcement algorithms. CuLE overcomes many limitations of existing CPU-based emulators and scales naturally to multiple GPUs. It leverages GPU parallelization to run thousands of games simultaneously and it renders frames directly on the GPU, to avoid the bottleneck arising from the limited CPU-GPU communication bandwidth. CuLE generates up to 155M frames per hour on a single GPU, a finding previously achieved only through a cluster of CPUs. Beyond highlighting the differences between CPU and GPU emulators in the context of reinforcement learning, we show how to leverage the high throughput of CuLE by effective batching of the training data, and show accelerated convergence for A2C+V-trace. CuLE is available at https://github.com/NVLabs/cule .


title: Shaped On-Off Keying Using Polar Codes

abstract: The probabilistic shaping scheme from Honda and Yamamoto (2013) for polar codes is used to enable power-efficient signaling for on-off keying (OOK). As OOK has a non-symmetric optimal input distribution, shaping approaches that are based on the concatenation of a distribution matcher followed by systematic encoding do not result in optimal signaling. Instead, these approaches represent a time sharing scheme where only a fraction of the codeword symbols is shaped. The proposed scheme uses a polar code for joint distribution matching and forward error correction which enables asymptotically optimal signaling. Numerical simulations show a gain of 1.8 dB compared to uniform transmission at a spectral efficiency of 0.25 bits/channel use for a blocklength of 65,536 bits.


title: Exploring sentence informativeness

abstract: This study is a preliminary exploration of the concept of informativeness -how much information a sentence gives about a word it contains- and its potential benefits to building quality word representations from scarce data. We propose several sentence-level classifiers to predict informativeness, and we perform a manual annotation on a set of sentences. We conclude that these two measures correspond to different notions of informativeness. However, our experiments show that using the classifiers' predictions to train word embeddings has an impact on embedding quality.


title: Provenance Analysis for Logic and Games

abstract: A model checking computation checks whether a given logical sentence is true in a given finite structure. Provenance analysis abstracts from such a computation mathematical information on how the result depends on the atomic data that describe the structure. In database theory, provenance analysis by interpretations in commutative semirings has been rather succesful for positive query languages (such a unions of conjunctive queries, positive relational algebra, or datalog). However, it did not really offer an adequate treatment of negation or missing information. Here we propose a new approach for the provenance analysis of logics with negation, such as first-order logic and fixed-point logics. It is closely related to a provenance analysis of the associated model-checking games, and based on new semirings of dual-indeterminate polynomials or dual-indeterminate formal power series. These are obtained by taking quotients of traditional provenance semirings by congruences that are generated by products of positive and negative provenance tokens. Beyond the use for model-checking problems in logics, provenance analysis of games is of independent interest. Provenance values in games provide detailed information about the number and properties of the strategies of the players, far beyond the question whether or not a player has a winning strategy from a given position.


title: Effect of nucleation on icy pebble growth in protoplanetary discs

abstract: Solid particles in protoplanetary discs can grow by direct vapour deposition outside of ice lines. The presence of microscopic silicate particles may nevertheless hinder growth into large pebbles, since the available vapour is deposited predominantly on the small grains that dominate the total surface area. Experiments on heterogeneous ice nucleation, performed to understand ice clouds in the Martian atmosphere, show that the formation of a new ice layer on a silicate surface requires a substantially higher water vapour pressure than the deposition of water vapour on an existing ice surface. In this paper, we investigate how the difference in partial vapour pressure needed for deposition of vapour on water ice versus heterogeneous ice nucleation on silicate grains influences particle growth close to the water ice line. We developed and tested a dynamical 1D deposition and sublimation model, where we include radial drift, sedimentation, and diffusion in a turbulent protoplanetary disc. We find that vapour is deposited predominantly on already ice-covered particles, since the vapour pressure exterior of the ice line is too low for heterogeneous nucleation on bare silicate grains. Icy particles can thus grow to centimetre-sized pebbles in a narrow region around the ice line, whereas silicate particles stay dust-sized and diffuse out over the disc. The inhibition of heterogeneous ice nucleation results in a preferential region for growth into planetesimals close to the ice line where we find large icy pebbles. The suppression of heterogeneous ice nucleation on silicate grains may also be the mechanism behind some of the observed dark rings around ice lines in protoplanetary discs, as the presence of large ice pebbles outside ice lines leads to a decrease in the opacity there.


title: Time-independent free energies from metadynamics via Mean Force Integration

abstract: Inspired by thermodynamic integration, we propose a method for the calculation of time-independent free energy profiles from history-dependent biased simulations via Mean Force Integration (MFI). MFI circumvents the need for computing the ensemble average of the bias acting on the system c(t) and can be applied to different variants of metadynamics. Moreover, MFI naturally extends to aggregate information obtained from independent metadynamics simulations, allowing to converge free energy surfaces without the need to sample recrossing events in a single continuous trajectory. We validate MFI against one- and two-dimensional analytical potentials and by computing the conformational free energy landscape of ibuprofen in the bulk of its most common crystal phase.


title: The Static Local Field Correction of the Warm Dense Electron Gas: An ab Initio Path Integral Monte Carlo Study and Machine Learning Representation

abstract: The study of matter at extreme densities and temperatures as they occur in astrophysical objects and state-of-the art experiments with high-intensity lasers is of high current interest for many applications. While no overarching theory for this regime exists, accurate data for the density response of correlated electrons to an external perturbation are of paramount importance. In this context, the key quantity is given by the local field correction (LFC), which provides a wave-vector resolved description of exchange-correlation effects. In this work, we present extensive new path integral Monte Carlo (PIMC) results for the static LFC of the uniform electron gas, which are subsequently used to train a fully connected deep neural network. This allows us to present a continuous representation of the LFC with respect to wave-vector, density, and temperature covering the entire warm dense matter regime. Both the PIMC data and neural-net results are available online. Moreover, we expect the presented combination of ab initio calculations with machine-learning methods to be a promising strategy for many applications.


title: A Practical Fixed-Parameter Algorithm for Constructing Tree-Child Networks from Multiple Binary Trees

abstract: We present the first fixed-parameter algorithm for constructing a tree-child phylogenetic network that displays an arbitrary number of binary input trees and has the minimum number of reticulations among all such networks. The algorithm uses the recently introduced framework of cherry picking sequences and runs in $O((8k)^k \mathrm{poly}(n, t))$ time, where $n$ is the number of leaves of every tree, $t$ is the number of trees, and $k$ is the reticulation number of the constructed network. Moreover, we provide an efficient parallel implementation of the algorithm and show that it can deal with up to $100$ input trees on a standard desktop computer, thereby providing a major improvement over previous phylogenetic network construction methods.


title: Representational Capacity of Deep Neural Networks -- A Computing Study

abstract: There is some theoretical evidence that deep neural networks with multiple hidden layers have a potential for more efficient representation of multidimensional mappings than shallow networks with a single hidden layer. The question is whether it is possible to exploit this theoretical advantage for finding such representations with help of numerical training methods. Tests using prototypical problems with a known mean square minimum did not confirm this hypothesis. Minima found with the help of deep networks have always been worse than those found using shallow networks. This does not directly contradict the theoretical findings---it is possible that the superior representational capacity of deep networks is genuine while finding the mean square minimum of such deep networks is a substantially harder problem than with shallow ones.


title: A conformal Yamabe problem with potential on the euclidean space

abstract: We consider, in the Euclidean setting, a conformal Yamabe-type equation related to a potential generalization of the classical constant scalar curvature problem and which naturally arises in the study of Ricci solitons structures. We prove existence and nonexistence results, focusing on the radial case, under some general hypothesis on the potential.


title: A polynomial bound for the number of maximal systems of imprimitivity of a finite transitive permutation group

abstract: We show that, there exists a constant $a$ such that, for every subgroup $H$ of a finite group $G$, the number of maximal subgroups of $G$ containing $H$ is bounded above by $a|G:H|^{3/2}$. In particular, a transitive permutation group of degree $n$ has at most $an^{3/2}$ maximal systems of imprimitivity. When $G$ is soluble, generalizing a classic result of Tim Wall, we prove a much stroger bound, that is, the number of maximal subgroups of $G$ containing $H$ is at most $|G:H|-1$.


title: Interactive Learning of Environment Dynamics for Sequential Tasks

abstract: In order for robots and other artificial agents to efficiently learn to perform useful tasks defined by an end user, they must understand not only the goals of those tasks, but also the structure and dynamics of that user's environment. While existing work has looked at how the goals of a task can be inferred from a human teacher, the agent is often left to learn about the environment on its own. To address this limitation, we develop an algorithm, Behavior Aware Modeling (BAM), which incorporates a teacher's knowledge into a model of the transition dynamics of an agent's environment. We evaluate BAM both in simulation and with real human teachers, learning from a combination of task demonstrations and evaluative feedback, and show that it can outperform approaches which do not explicitly consider this source of dynamics knowledge.


title: An approximate version of Jackson's conjecture

abstract: In 1981 Jackson showed that the diregular bipartite tournament (a complete bipartite graph whose edges are oriented so that every vertex has the same in- and outdegree) contains a Hamilton cycle, and conjectured that in fact the edge set of it can be partitioned into Hamilton cycles. We prove an approximate version of this conjecture: For every $c>1/2$ and $\varepsilon>0$ there exists $n_0$ such that every $cn$-regular bipartite digraph on $2n\geq n_0$ vertices contains $(1-\varepsilon)cn$ edge-disjoint Hamilton cycles.


title: Microscopic reweighting for non-equilibrium steady states dynamics

abstract: Computer simulations generate trajectories at a single, well-defined thermodynamic state point. Statistical reweighting offers the means to reweight static and dynamical properties to different equilibrium state points by means of analytic relations. We extend these ideas to non-equilibrium steady states by relying on a maximum path entropy formalism subject to physical constraints. Stochastic thermodynamics analytically relates the forward and backward probabilities of any pathway through the external non-conservative force, enabling reweighting both in and out of equilibrium. We avoid the combinatorial explosion of microtrajectories by systematically constructing pathways through Markovian transitions. We further identify a quantity that is invariant to dynamical reweighting, analogous to the density of states in equilibrium reweighting.


title: A spatially extended model to assess the role of landscape structure on the pollination service of Apis mellifera

abstract: Apis mellifera plays a crucial role as pollinator of the majority of crops linked to food production and thus its presence is currently fundamental to our health and survival. The composition and configuration of the landscape in which Apis mellifera lives will likely determine the well-being of the hives and the pollination service that their members can provide to the crops. Here we present a spatially explicit model that predicts the spatial distribution of visits by Apis mellifera to crops, by simulating daily trips of honey bees, the demographical dynamic of each hive and their honey production. This model goes beyond existing approaches by including 1) a flower resource affected by the feedback interaction between nectar extraction, pollination, blossoming and repeated visits, 2) a pollinators dynamic that allows competition through short term resource depletion, 3) a probabilistic approach of the foraging behavior, modeling the fact that the pollinators have only partial knowledge of the resource on their surroundings, and 4) the specific and systematic foraging behavior and strategies of Apis mellifera at the moment of choosing foraging sites, as opposed to those adopted by solitary and wild pollinators. With a balance between simplicity and realism we show the importance of keeping a minimal fraction of natural habitat in an agricultural landscape. We also evaluate the effects of the landscape's structure on pollination, and demonstrate that there exists an optimal size of natural habitat patches that maximizes the pollination service for a fixed fraction of natural habitat.


title: Impact of main ion pressure anisotropy on stellarator impurity transport

abstract: Main ions influence impurity dynamics through a variety of mechanisms; in particular, via impurity-ion collisions. To lowest order in an expansion in the main ion mass over the impurity mass, the impurity-ion collision operator only depends on the component of the main ion distribution that is odd in the parallel velocity. These lowest order terms give the parallel friction of the impurities with the main ions, which is typically assumed to be the main cause of collisional impurity transport. Next-order terms in the mass ratio expansion of the impurity-ion collision operator, proportional to the component of the main ion distribution that is even in the parallel velocity, are usually neglected. However, in stellarators, the even component of the main ion distribution can be very large. In this article, such next-order terms in the mass ratio expansion of the impurity-ion collision operator are retained, and analytical expressions for the neoclassical radial flux of trace impurities are calculated in the Pfirsch-Schl\"uter, plateau and $1/\nu$ regimes. The new terms provide a drive for impurity transport that is physically very different from parallel friction: they are associated to anisotropy in the pressure of the main ions, which translates into impurity pressure anisotropy. It is argued that main ion pressure anisotropy must be taken into account for a correct description of impurity transport in certain realistic stellarator plasmas. Examples are given by numerically evaluating the analytical expressions for the impurity flux.


title: Inferring Accurate Bus Trajectories from Noisy Estimated Arrival Time Records

abstract: Urban commuting data has long been a vital source of understanding population mobility behaviour and has been widely adopted for various applications such as transport infrastructure planning and urban anomaly detection. While individual-specific transaction records (such as smart card (tap-in, tap-out) data or taxi trip records) hold a wealth of information, these are often private data available only to the service provider (e.g., taxicab operator). In this work, we explore the utility in harnessing publicly available, albeit noisy, transportation datasets, such as noisy "Estimated Time of Arrival" (ETA) records (commonly available to commuters through transit Apps or electronic signages). We first propose a framework to extract accurate individual bus trajectories from such ETA records, and present results from both a primary city (Singapore) and a secondary city (London) to validate the techniques. Finally, we quantify the upper bound on the spatiotemporal resolution, of the reconstructed trajectory outputs, achieved by our proposed technique.


title: Trojans in the Solar Neighborhood

abstract: About 20% of stars in the solar vicinity are in the Hercules stream, a bundle of stars that move together with a velocity distinct from the Sun. Its origin is still uncertain. Here, we explore the possibility that Hercules is made of trojans, stars captured at L4, one the Lagrangian points of the stellar bar. Using GALAKOS--a high-resolution N-body simulation of the Galactic disk--we follow the motions of stars in the co-rotating frame of the bar and confirm previous studies on Hercules being formed by stars in co-rotation resonance with the bar. Unlike previous work, we demonstrate that the retrograde nature of trojan orbits causes the asymmetry in the radial velocity distribution, typical of Hercules in the solar vicinity. We show that trojans remain at capture for only a finite amount of time, before escaping L4 without being captured again. We anticipate that in the kinematic plane the Hercules stream will de-populate along the bar major axis and be visible at azimuthal angles behind the solar vicinity with a peak towards L4. This test can exclude the OLR origin of the Hercules stream and be validated by Gaia DR3 and DR4.


title: Invariant measure for stochastic Schr\"odinger equations

abstract: Quantum trajectories are Markov processes that describe the time-evolution of a quantum system undergoing continuous indirect measurement. Mathematically, they are defined as solutions of the so-called "Stochastic Schr\"odinger Equations", which are nonlinear stochastic differential equations driven by Poisson and Wiener processes. This paper is devoted to the study of the invariant measures of quantum trajectories. Particularly, we prove that the invariant measure is unique under an ergodicity condition on the mean time evolution, and a "purification" condition on the generator of the evolution. We further show that quantum trajectories converge in law exponentially fast towards this invariant measure. We illustrate our results with examples where we can derive explicit expressions for the invariant measure.


title: Discovery of a giant and luminous Lya+CIV+HeII nebula at z=3.326 with extreme emission line ratios

abstract: We present the discovery of HLock01-LAB, a luminous and large Lya nebula at z=3.326. Medium-band imaging and long-slit spectroscopic observations with the Gran Telescopio Canarias reveal extended emission in the Lya 1215\AA, CIV1550\AA, and HeII 1640\AA lines over ~100kpc, and a total luminosity L(Lya)=(6.4+/-0.1)x10^44 erg s^-1. HLock01-LAB presents an elongated morphology aligned with two faint radio sources contained within the central ~8kpc of the nebula. The radio structures are consistent to be faint radio jets or lobes of a central galaxy, whose spectrum shows nebular emission characteristic of a type-II active galactic nucleus (AGN). The continuum emission of the AGN at short wavelengths is, however, likely dominated by stellar emission of the host galaxy, for which we derive a stellar mass M* = 2.3x10^11 Msun. The detection of extended emission in CIV and CIII] indicates that the gas within the nebula is not primordial. Feedback may have enriched the halo at at least 50 kpc from the nuclear region. Using rest-frame UV emission-line diagnostics, we find that the gas in the nebula is likely heated by the AGN. Nevertheless, at the center of the nebula we find extreme emission line ratios of Lya/CIV~60 and Lya/HeII~80, one of the highest values measured to date, and well above the standard values of photoionization models (Lya/HeII~30 for case B photoionization). Our data suggest that jet-induced shocks are likely responsible for the increase of the electron temperature and, thus, the observed Lya enhancement in the center of the nebula. This scenario is further supported by the presence of radio structures and perturbed kinematics in this region. The large Lya luminosity in HLock01-LAB is likely due to a combination of AGN photoionization and jet-induced shocks, highlighting the diversity of sources of energy powering Lya nebulae. [abridged]


title: A Graph Neural Network Approach for Scalable Wireless Power Control

abstract: Deep neural networks have recently emerged as a disruptive technology to solve NP-hard wireless resource allocation problems in a real-time manner. However, the adopted neural network structures, e.g., multi-layer perceptron (MLP) and convolutional neural network (CNN), are inherited from deep learning for image processing tasks, and thus are not tailored to problems in wireless networks. In particular, the performance of these methods deteriorates dramatically when the wireless network size becomes large. In this paper, we propose to utilize graph neural networks (GNNs) to develop scalable methods for solving the power control problem in $K$-user interference channels. Specifically, a $K$-user interference channel is first modeled as a complete graph, where the quantitative information of wireless channels is incorporated as the features of the graph. We then propose an interference graph convolutional neural network (IGCNet) to learn the optimal power control in an unsupervised manner. It is shown that one-layer IGCNet is a universal approximator to continuous set functions, which well matches the permutation invariance property of interference channels and it is robust to imperfect channel state information (CSI). Extensive simulations will show that the proposed IGCNet outperforms existing methods and achieves significant speedup over the classic algorithm for power control, namely, WMMSE. The code is available on https://github.com/yshenaw/Globecom2019.


title: An Optimal Control Approach to Early Stopping Variational Methods for Image Restoration

abstract: We investigate a well-known phenomenon of variational approaches in image processing, where typically the best image quality is achieved when the gradient flow process is stopped before converging to a stationary point. This paradox originates from a tradeoff between optimization and modelling errors of the underlying variational model and holds true even if deep learning methods are used to learn highly expressive regularizers from data. In this paper, we take advantage of this paradox and introduce an optimal stopping time into the gradient flow process, which in turn is learned from data by means of an optimal control approach. As a result, we obtain highly efficient numerical schemes that achieve competitive results for image denoising and image deblurring. A nonlinear spectral analysis of the gradient of the learned regularizer gives enlightening insights about the different regularization properties.


title: Empowering A* Search Algorithms with Neural Networks for Personalized Route Recommendation

abstract: Personalized Route Recommendation (PRR) aims to generate user-specific route suggestions in response to users' route queries. Early studies cast the PRR task as a pathfinding problem on graphs, and adopt adapted search algorithms by integrating heuristic strategies. Although these methods are effective to some extent, they require setting the cost functions with heuristics. In addition, it is difficult to utilize useful context information in the search procedure. To address these issues, we propose using neural networks to automatically learn the cost functions of a classic heuristic algorithm, namely A* algorithm, for the PRR task. Our model consists of two components. First, we employ attention-based Recurrent Neural Networks (RNN) to model the cost from the source to the candidate location by incorporating useful context information. Instead of learning a single cost value, the RNN component is able to learn a time-varying vectorized representation for the moving state of a user. Second, we propose to use a value network for estimating the cost from a candidate location to the destination. For capturing structural characteristics, the value network is built on top of improved graph attention networks by incorporating the moving state of a user and other context information. The two components are integrated in a principled way for deriving a more accurate cost of a candidate location. Extensive experiment results on three real-world datasets have shown the effectiveness and robustness of the proposed model.


title: Automatic Generation of Acceptance Test Cases from Use Case Specifications: an NLP-based Approach

abstract: Acceptance testing is a validation activity performed to ensure the conformance of software systems with respect to their functional requirements. In safety critical systems, it plays a crucial role since it is enforced by software standards. Test engineers need to identify all the representative test execution scenarios from requirements, determine the runtime conditions that trigger these scenarios, and finally provide the input data that satisfy these conditions. Given that requirements specifications are typically large and often provided in natural language, the generation of acceptance test cases tends to be expensive and error-prone. In this paper, we present UMTG, an approach that supports the generation of executable, system-level, acceptance test cases from requirements specifications in natural language, with the goal of reducing the manual effort required to generate test cases and ensuring requirements coverage. More specifically, UMTG automates the generation of acceptance test cases based on use case specifications and a domain model for the system under test, which are commonly produced in many development environments. Unlike existing approaches, it does not impose strong restrictions on the expressiveness of use case specifications. We rely on recent advances in natural language processing to automatically identify test scenarios and to generate formal constraints that capture conditions triggering the execution of the scenarios, thus enabling the generation of test data. In two industrial case studies, UMTG automatically and correctly translated 95% of the use case specification steps into formal constraints required for test data generation; furthermore, it generated test cases that exercise not only all the test scenarios manually implemented by experts, but also some critical scenarios not previously considered.


title: Are Parametric Markov Chains Monotonic?

abstract: This paper presents a simple algorithm to check whether reachability probabilities in parametric Markov chains are monotonic in (some of) the parameters. The idea is to construct - only using the graph structure of the Markov chain and local transition probabilities - a pre-order on the states. Our algorithm cheaply checks a sufficient condition for monotonicity. Experiments show that monotonicity in several benchmarks is automatically detected, and monotonicity can speed up parameter synthesis up to orders of magnitude faster than a symbolic baseline.


title: A Hermite-like basis for faster matrix-free evaluation of interior penalty discontinuous Galerkin operators

abstract: This work proposes a basis for improved throughput of matrix-free evaluation of discontinuous Galerkin symmetric interior penalty discretizations on hexahedral elements. The basis relies on ideas of Hermite polynomials. It is used in a fully discontinuous setting not for higher order continuity but to minimize the effective stencil width, namely to limit the neighbor access of an element to one data point for the function value and one for the derivative. The basis is extended to higher orders with nodal contributions derived from roots of Jacobi polynomials and extended to multiple dimensions with tensor products, which enable the use of sum factorization. The beneficial effect of the reduced data access on modern processors is shown. Furthermore, the viability of the basis in the context of multigrid solvers is analyzed. While a plain point-Jacobi approach is less efficient than with the best nodal polynomials, a basis change via sum-factorization techniques enables the combination of the fast matrix-vector products with effective multigrid constituents. The basis change is essentially for free on modern hardware because these computations can be hidden behind the cost of the data access.


title: A note on the local Lipschitz triviality of values of complex polynomial functions

abstract: We address the question of the bi-Lipschitz local triviality of a complex polynomial function over a complex value. Our main result state that a non constant complex polynomial admits a locally bi-Lipschitz trivial value if and only if it is a polynomial in a single complex variable.


title: Performance evaluation of THz wireless systems under the joint impact of misalignment fading and phase noise

abstract: In this paper, we investigate the joint impact of misalignment fading and local oscillator (LO) phase noise (PHN) in multi-carrier terahertz (THz) wireless systems. In more detail, after establishing a suitable system model that takes into account the particularities of the THz channel, as well as the transceivers characteristics, we present simulation results that quantify the joint impact of misalignment fading and PHN in terms of average signal-to-interference-plus-noise-ratio (SINR) and outage probability (OP).


title: Criteria for Rigor in Visualization Design Study

abstract: We develop a new perspective on research conducted through visualization design study that emphasizes design as a method of inquiry and the broad range of knowledge-contributions achieved through it as multiple, subjective, and socially constructed. From this interpretivist position we explore the nature of visualization design study and develop six criteria for rigor. We propose that rigor is established and judged according to the extent to which visualization design study research and its reporting are INFORMED, REFLEXIVE, ABUNDANT, PLAUSIBLE, RESONANT, and TRANSPARENT. This perspective and the criteria were constructed through a four-year engagement with the discourse around rigor and the nature of knowledge in social science, information systems, and design. We suggest methods from cognate disciplines that can support visualization researchers in meeting these criteria during the planning, execution, and reporting of design study. Through a series of deliberately provocative questions, we explore implications of this new perspective for design study research in visualization, concluding that as a discipline, visualization is not yet well positioned to embrace, nurture, and fully benefit from a rigorous, interpretivist approach to design study. The perspective and criteria we present are intended to stimulate dialogue and debate around the nature of visualization design study and the broader underpinnings of the discipline.


title: Combinatorial parametrised spectra

abstract: We obtain combinatorial model categories of parametrised spectra, together with systems of base change Quillen adjunctions associated to maps of parameter spaces. We work with simplicial objects and use Hovey's sequential and symmetric stabilisation machines. By means of a Grothendieck construction for model categories, we produce combinatorial model categories controlling the totality of parametrised stable homotopy theory. The global model category of parametrised symmetric spectra is equipped with a symmetric monoidal model structure (the external smash product) inducing pairings in twisted cohomology groups. As an application of our results we prove a tangent prolongation of Simpson's theorem, characterising tangent $\infty$-categories of presentable $\infty$-categories as accessible localisations of $\infty$-categories of presheaves of parametrised spectra. Applying these results to the homotopy theory of smooth $\infty$-stacks produces well-behaved (symmetric monoidal) model categories of smooth parametrised spectra. These models provide a concrete foundation for studying twisted differential cohomology, incorporating previous work of Bunke and Nikolaus.


title: Accurate Sampling with Noisy Forces from Approximate Computing

abstract: In scientific computing, the acceleration of atomistic computer simulations by means of custom hardware is finding ever growing application. A major limitation, however, is that the high efficiency in terms of performance and low power consumption entails the massive usage of low-precision computing units. Here, based on the approximate computing paradigm, we present an algorithmic method to rigorously compensate for numerical inaccuracies due to low-accuracy arithmetic operations, yet still obtaining exact expectation values using a properly modified Langevin-type equation.

categories: physics.comp-ph cond-mat.stat-mech cs.NA math.NA physics.chem-ph


title: Dynamic multi-focus laser writing with acousto-optofluidics

abstract: Laser writing of materials is normally performed by the sequential scanning of a single focused beam across a sample. This process is time-consuming and it can severely limit the throughput of laser systems in key applications such as surgery, microelectronics, or manufacturing. Here we report a parallelization strategy based on ultrasound waves in a liquid to diffract light into multiple beamlets. Adjusting amplitude, frequency, or phase of ultrasound allows tunable multi-focus distributions with sub-microsecond control. When combined with sample translation, the dynamic splitting of light leads to high-throughput laser processing, as demonstrated by locally modifying the morphological and wettability properties of metals, polymers, and ceramics. The results illustrate how acousto-optofluidic systems are universal tools for fast multi-focus generation, with potential impact in fields such as imaging or optical trapping.


title: L\'evy-Ito Models in Finance

abstract: We present an overview of the broad class of financial models in which the prices of assets are L\'evy-Ito processes driven by an $n$-dimensional Brownian motion and an independent Poisson random measure. The Poisson random measure is associated with an $n$-dimensional L\'evy process. Each model consists of a pricing kernel, a money market account, and one or more risky assets. We show how the excess rate of return above the interest rate can be calculated for risky assets in such models, thus showing the relationship between risk and return when asset prices have jumps. The framework is applied to a variety of asset classes, allowing one to construct new models as well as interesting generalizations of familiar models.


title: Network-Assisted D2D Relay Selection Under the Presence of Dynamic Obstacles

abstract: Millimeter wave (\texttt{mmWave}) channels in device to device (\texttt{D2D}) communication are susceptible to blockages in spite of using directional beams from multi-input multi-output (\texttt{MIMO}) antennas to compensate for high propagation loss. This motivates one to look for the presence of obstacles while forming \texttt{D2D} links among user equipments (\texttt{UEs}) which are in motion. In \texttt{D2D} communication, moving \texttt{UEs} also act as relays to forward data from one \texttt{UE} to another which introduces the problem of relay selection. The problem becomes more challenging when the obstacles are also in motion (dynamic obstacles) along with the moving \texttt{UEs}. First we have developed a probabilistic model for relay selection which considers both moving \texttt{UEs} and dynamic obstacles. Then we have analyzed the probability of dynamic obstacles blocking a link in 3D Euclidean space by exploiting the information from \texttt{MIMO} radar connected to the base station. Finally, using this information, we have developed unique strategies based on simple geometry to find the best relay which maximizes the expected data rate. Through simulations we have shown that our proposed strategy gives a significant improvement in packet loss due to mobility of nodes and dynamic obstacles in a \texttt{mmWave} channel over traditional approaches which do not consider dynamic obstacle's presence.


title: A Comparative Evaluation of Visual and Natural Language Question Answering Over Linked Data

abstract: With the growing number and size of Linked Data datasets, it is crucial to make the data accessible and useful for users without knowledge of formal query languages. Two approaches towards this goal are knowledge graph visualization and natural language interfaces. Here, we investigate specifically question answering (QA) over Linked Data by comparing a diagrammatic visual approach with existing natural language-based systems. Given a QA benchmark (QALD7), we evaluate a visual method which is based on iteratively creating diagrams until the answer is found, against four QA systems that have natural language queries as input. Besides other benefits, the visual approach provides higher performance, but also requires more manual input. The results indicate that the methods can be used complementary, and that such a combination has a large positive impact on QA performance, and also facilitates additional features such as data exploration.


title: Particle-In-Cell simulations of the parallel proton firehose instability influenced by the electron temperature anisotropy in solar wind conditions

abstract: In situ observations of the solar wind show a limited level of particle temperature anisotropy with respect to the interplanetary magnetic field direction. Kinetic electromagnetic instabilities are efficient to prevent the excessive growth of the anisotropy of particle velocity distribution functions. Among them, the firehose instabilities are often considered to prevent the increase of the parallel temperature and hence to shape the velocity distribution functions of electrons and protons in the solar wind. We present a non-linear modeling of the parallel firehose instability, retaining a kinetic description for both the electrons and protons. One-dimensional (1D) fully kinetic Particle-In-Cell simulations using the Energy Conserving semi-implicit method (ECsim) are performed to clarify the role of the electron temperature anisotropy in the development of the parallel proton firehose instability. We found that in the presence of an electron temperature anisotropy, such that the temperature parallel to the background magnetic field is higher than the temperature in the perpendicular direction, the onset of the parallel proton firehose instability occurs earlier and its growth rate is faster. The enhanced wave fluctuations contribute to the particle scattering reducing the temperature anisotropy to a stable, nearly isotropic state. The simulation results compare well with linear theory. A test case of 1D simulations at oblique angles with respect to the magnetic field is also considered, as a first step to study the cumulative effect of protons and electrons on the full spectrum of instabilities.


title: Unified phase diagram of reversible-irreversible, jamming and yielding transitions in cyclically sheared soft sphere packings

abstract: Self-organization, and transitions from reversible to irreversible behaviour, of interacting particle assemblies driven by externally imposed stresses or deformation is of interest in comprehending diverse phenomena in soft matter. They have been investigated in a wide range of systems, such as colloidal suspensions, glasses, and granular matter. In different density and driving regimes, such behaviour is related to yielding of amorphous solids, jamming, and memory formation, \emph{etc.} How these phenomena are related to each other has not, however, been much studied. In order to obtain a unified view of the different regimes of behaviour, and transitions between them, we investigate computationally the response of soft sphere assemblies to athermal cyclic shear deformation over a wide range of densities and amplitudes of shear deformation. Cyclic shear deformation induces transitions from reversible to irreversible behaviour in both unjammed and jammed soft sphere packings. Well above isotropic jamming density ($\bf{\phi_J}$), this transition corresponds to yielding. In the vicinity of the jamming point, up to a higher density limit we designate ${\bf \phi_J^{cyc}}$, an unjammed phase emerges between a localised, \emph{absorbing} phase, and a diffusive, {\emph irreversible} phase. The emergence of the unjammed phase signals the shifting of the jamming point to higher densities as a result of annealing, and opens a window where shear jamming becomes possible for frictionless packings. Below $\bf{\phi_J}$, two distinct localised states, termed point and loop reversibile, are observed. We characterise in detail the different regimes and transitions between them, and obtain a unified density-shear amplitude phase diagram.


title: Time-Optimal Control of Collisional $\sqrt{\text{SWAP}}$ Gates in Ultracold Atomic Systems

abstract: We use quantum optimal control to identify fast collision-based two-qubit $\sqrt{\text{SWAP}}$ gates in ultracold atoms. We show that a significant speed up can be achieved by optimizing the full gate instead of separately optimizing the merge-wait-separate sequence of the trapping potentials. Our optimal strategy does not rely on the atoms populating the lowest eigenstates of the merged potential, and it crucially includes accumulation of quantum phases before the potentials are fully merged. Our analyses transcend the particular trapping geometry, but to compare with previous works, we present systematic results for an optical lattice and find greatly improved gate durations and fidelities.


title: Tidal disruption events as a site of an evolving relativistic spectral~line

abstract: In nuclei of galaxies strong tidal forces can destroy stars passing within a critical distance from the central super-massive black hole (SMBH). Observational signatures of tidal disruption events (TDEs) depend on the environment around the SBMH horizon and the level of its accretion activity. Evidence for optical and UV spectral features has been reported in TDE flares; nevertheless, to test the effects of General Relativity in the immediate vicinity of SMBH, the relativistically broadened and skewed X-ray line would tell us significantly more useful information. This will require to proceed beyond inactive nuclei. To this end we consider a system where the material from a disrupted star forms a gaseous ring that circularises near the tidal radius around SMBH, it gradually spreads in radius by viscous processes, and resides embedded within a hot corona. In our test calculation the remnant trail is assumed to be fully circularized and embedded in a hot environment and illuminated by X-rays from a surrounding corona or a jet base of (mild) AGN activity. We show the expected effects on the observed profile and the centroid energy. In future the evolving spectral features can enhance the diagnostic capability and provide a novel way to reveal the parameters of TDEs in such sources; namely, the distance of the remnant gas from the SMBH, the radial extent of the gaseous trail, and the spin of the SMBH could be measured.


title: Language Modelling for Sound Event Detection with Teacher Forcing and Scheduled Sampling

abstract: A sound event detection (SED) method typically takes as an input a sequence of audio frames and predicts the activities of sound events in each frame. In real-life recordings, the sound events exhibit some temporal structure: for instance, a "car horn" will likely be followed by a "car passing by". While this temporal structure is widely exploited in sequence prediction tasks (e.g., in machine translation), where language models (LM) are exploited, it is not satisfactorily modeled in SED. In this work we propose a method which allows a recurrent neural network (RNN) to learn an LM for the SED task. The method conditions the input of the RNN with the activities of classes at the previous time step. We evaluate our method using F1 score and error rate (ER) over three different and publicly available datasets; the TUT-SED Synthetic 2016 and the TUT Sound Events 2016 and 2017 datasets. The obtained results show an increase of 9% and 2% at the F1 (higher is better) and a decrease of 7% and 2% at ER (lower is better) for the TUT Sound Events 2016 and 2017 datasets, respectively, when using our method. On the contrary, with our method there is a decrease of 4% at F1 score and an increase of 7% at ER for the TUT-SED Synthetic 2016 dataset.


title: A Short Proof of Bernoulli Disjointness via the Local Lemma

abstract: Recently, Glasner, Tsankov, Weiss, and Zucker showed that if $\Gamma$ is an infinite discrete group, then every minimal $\Gamma$-flow is disjoint from the Bernoulli shift $2^\Gamma$. Their proof is somewhat involved; in particular, it invokes separate arguments for different classes of groups. In this note, we give a short and self-contained proof of their result using purely combinatorial methods applicable to all groups at once. Our proof relies on the Lov\'{a}sz Local Lemma, an important tool in probabilistic combinatorics that has recently found several applications in the study of dynamical systems.


title: An Inverse Mass Expansion for the Mutual Information in Free Scalar QFT at Finite Temperature

abstract: We study the entanglement entropy and the mutual information in coupled harmonic systems at finite temperature. Interestingly, we find that the mutual information does not vanish at infinite temperature, but it rather reaches a specific finite value, which can be attributed to classical correlations solely. We further obtain high and low temperature expansions for both quantities. Then, we extend the analysis performed in the seminal paper by Srednicki (Phys. Rev. Lett. 71, 666 (1993)) for free real scalar field theories in Minkowski space-time in 3+1 dimensions at a thermal state. We find that the mutual information obeys an area law, similar to that obeyed by the entanglement entropy at vanishing temperature. The coefficient of this area law does not vanish at infinite temperature. Then, we calculate this coefficient perturbatively in an $1/\mu$ expansion, where $\mu$ is the mass of the scalar field. Finally, we study the high and low temperature behaviour of the area law term.


title: A fibre Smart Displacement Based (FSDB) beam element for the nonlinear analysis of R/C members

abstract: Beam finite elements for non linear plastic analysis of beam-like structures are formulated according to Displacement Based (DB) or Force Based (FB) approaches. DB formulations rely on modelling the displacement field by means of displacement shape functions. Despite the greater simplicity of DB over FB approaches, the latter provide more accurate responses requiring a coarser mesh. To fill the existent gap between the two approaches the improvement of the DB formulation without the introduction of mesh refinement is needed. To this aim the authors recently provided a contribution to the improvement of the DB approach by proposing new enriched adaptive displacement shape functions leading to the Smart Displacement Based (SDB) beam element. In this paper the SDB element is extended to include the axial force-bending moment interaction, crucial for the analysis of r/c cross sections. The proposed extension requires the formulation of discontinuous axial displacement shape functions dependent on the diffusion of plastic deformations. The stiffness matrix of the extended smart element is provided explicitly and shown to be dependent on the displacement shape functions updating. The axial force-bending moment interaction is approached by means of a fibre discretisation of the r/c cross section. The extended element, addressed as Fibre Smart Displacement Based (FSDB) beam element, is shown to be accurate and furthermore accompanied by the proposal of an optional procedure to verify strong equilibrium of the axial force along the beam element, which is usually not accomplished by DB beam elements. Given a fixed mesh discretisation the performance of the FSDB beam element is compared with the DB approach to show the better accuracy of the proposed element.


title: Using Muon Rings for the Optical Throughput Calibration of the Cherenkov Telescope Array

abstract: Muon ring images observed with Imaging Atmospheric Cherenkov Telescopes (IACTs) provide a powerful means to calibrate the optical throughput of IACTs and monitor their optical point spread function. We investigate whether muons ring images can be used as the primary optical throughput calibration method for the telescopes of the future Cherenkov Telescope Array (CTA) and find several additional systematic effects in comparison to previous works. To ensure that the method achieves the accuracy required by CTA, these systematic effects need to be taken into account and minor modifications to the hardware and analysis are necessary. We derive analytic estimates for the expected muon data rates to be used for optical throughput calibration, monitoring of the optical point spread function, with achievable statistical and systematic uncertainties, and explore the potential of muon ring images as a secondary method of camera pixel flat-fielding.


title: Matrix cofactorization for joint spatial-spectral unmixing of hyperspectral images

abstract: Hyperspectral unmixing aims at identifying a set of elementary spectra and the corresponding mixture coefficients for each pixel of an image. As the elementary spectra correspond to the reflectance spectra of real materials, they are often very correlated yielding an ill-conditioned problem. To enrich the model and to reduce ambiguity due to the high correlation, it is common to introduce spatial information to complement the spectral information. The most common way to introduce spatial information is to rely on a spatial regularization of the abundance maps. In this paper, instead of considering a simple but limited regularization process, spatial information is directly incorporated through the newly proposed context of spatial unmixing. Contextual features are extracted for each pixel and this additional set of observations is decomposed according to a linear model. Finally the spatial and spectral observations are unmixed jointly through a cofactorization model. In particular, this model introduces a coupling term used to identify clusters of shared spatial and spectral signatures. An evaluation of the proposed method is conducted on synthetic and real data and shows that results are accurate and also very meaningful since they describe both spatially and spectrally the various areas of the scene.


title: Fluctuations of the product of random matrices and generalized Lyapunov exponent

abstract: I present a general framework allowing to carry out explicit calculation of the moment generating function of random matrix products $\Pi_n=M_nM_{n-1}\cdots M_1$, where $M_i$'s are i.i.d.. Following Tutubalin [Theor. Probab. Appl. {\bf 10}, 15 (1965)], the calculation of the generating function is reduced to finding the largest eigenvalue of a certain transfer operator associated with a family of representations of the group. The formalism is illustrated by considering products of random matrices from the group $\mathrm{SL}(2,\mathbb{R})$ where explicit calculations are possible. For concreteness, I study in detail transfer matrix products for the one-dimensional Schr\"odinger equation where the random potential is a L\'evy noise (derivative of a L\'evy process). In this case, I obtain a general formula for the variance of $\ln||\Pi_n||$ and for the variance of $\ln|\psi(x)|$, where $\psi(x)$ is the wavefunction, in terms of a single integral involving the Fourier transform of the invariant density of the matrix product. Finally I discuss the continuum limit of random matrix products (matrices close to the identity ). In particular, I investigate a simple case where the spectral problem providing the generalized Lyapunov exponent can be solved exactly.


title: From probabilistic mechanics to quantum theory

abstract: We show that quantum theory (QT) is a substructure of classical probabilistic physics. The central quantity of the classical theory is Hamilton's function, which determines canonical equations, a corresponding flow, and a Liouville equation for a probability density. We extend this theory in two respects: (1) The same structure is defined for arbitrary observables. Thus we have all of the above entities generated not only by Hamilton's function but by every observable. (2) We introduce for each observable a phase space function representing the classical action. This is a redundant quantity in a classical context but indispensable for the transition to QT. The basic equations of the resulting theory take a "quantum-like" form, which allows for a simple derivation of QT by means of a projection to configuration space reported previously [Quantum Stud.:Math. Found. (2018) 5:219-227]. We obtain the most important relations of QT, namely the form of operators, Schr\"odinger's equation, eigenvalue equations, commutation relations, expectation values, and Born's rule. Implications for the interpretation of QT are discussed, as well as an alternative projection method allowing for a derivation of spin.


title: Predicting Visual Memory Schemas with Variational Autoencoders

abstract: Visual memory schema (VMS) maps show which regions of an image cause that image to be remembered or falsely remembered. Previous work has succeeded in generating low resolution VMS maps using convolutional neural networks. We instead approach this problem as an image-to-image translation task making use of a variational autoencoder. This approach allows us to generate higher resolution dual channel images that represent visual memory schemas, allowing us to evaluate predicted true memorability and false memorability separately. We also evaluate the relationship between VMS maps, predicted VMS maps, ground truth memorability scores, and predicted memorability scores.


title: Cosmological information content in redshift-space power spectrum of SDSS-like galaxies in the quasi-nonlinear regime up to $k=0.3\,h\,$Mpc$^{-1}$

abstract: Clustering properties and peculiar velocities of halos in large-scale structure carry a wealth of cosmological information over a wide range of scales from linear to nonlinear scales. We use halo catalogs in a suite of high-resolution $N$-body simulations to construct mock catalogs of galaxies that resemble the SDSS-like luminous early-type galaxies at three redshift bins in the range $0.15\le z\le 0.7$. To do this we include 10 nuisance parameters to model variations in halo-galaxy connections for each redshift bin; the halo occupation distribution, and the spatial and velocity distributions of galaxies in the host halos. We evaluate the Fisher information matrix for the redshift-space power spectrum of SDSS-like galaxies using different sets of the mock catalogs that are generated from changes in each of model parameters; cosmological parameters ($\sigma_8$ and $\Omega_{\rm m}$), the halo-galaxy connection parameters, and the cosmological distances ($D_{\rm A}$ and $H$ parameters at each redshift bin) for modeling an apparent geometrical distortion of the redshift-space power spectrum (the Alcock-Paczynski effect). We show that combining the monopole and quadrupole power spectra of galaxies allows for accurate estimations of the cosmological parameters and the cosmological distances, even after marginalization over the halo-galaxy parameters, by lifting the parameter degeneracies that are otherwise inevitable if either of the two spectra alone is used. When including the galaxy power spectrum information up to $k=0.3\,h\,$Mpc$^{-1}$, we find about factor of 6 gain in the cosmological information content of ($\sigma_8,\Omega_{\rm m},D_{\rm A}$'s and $H$'s) compared to $k=0.2\,h\,$Mpc$^{-1}$.


title: Cherenkov Telescope Array Science: Multi-wavelength and multi-messenger perspectives

abstract: The Cherenkov Telescope Array (CTA) will be the major global observatory for VHE gamma-ray astronomy over the next decade and beyond. It will be an explorer of the extreme universe, with a broad scientific potential: from understanding the role of relativistic cosmic particles, to the search for dark matter. Covering photon energies from 20 GeV to 300 TeV, and with an angular resolution unique in the field, of about 1 arc min, CTA will improve on all aspects of the performance with respect to current instruments, surveying the high energy sky hundreds of times faster than previous TeV telescopes, and with a much deeper view. The very large collection area of CTA makes it an important probe of transient phenomena. The first CTA telescope has just been inaugurated in the Canary Islands, Spain, and as more telescopes are added in the coming years, scientific operation will start. It is evident that CTA will have important synergies with many of the new generation astronomical and astroparticle observatories. In this talk we will review the CTA science case from the point of view of its synergies with other instruments and facilities, highlighting the CTA needs in terms of external data, as well as the opportunities and strategies for cooperation to achieve the basic CTA science goals.


title: Random cographs: Brownian graphon limit and asymptotic degree distribution

abstract: We consider uniform random cographs (either labeled or unlabeled) of large size. Our first main result is the convergence towards a Brownian limiting object in the space of graphons. We then show that the degree of a uniform random vertex in a uniform cograph is of order $n$, and converges after normalization to the Lebesgue measure on $[0,1]$. We finally analyze the vertex connectivity (i.e. the minimal number of vertices whose removal disconnects the graph) of random connected cographs, and show that this statistics converges in distribution without renormalization. Unlike for the graphon limit and for the degree of a random vertex, the limiting distribution is different in the labeled and unlabeled settings. Our proofs rely on the classical encoding of cographs via cotrees. We then use mainly combinatorial arguments, including the symbolic method and singularity analysis.


title: Mitigation of readout noise in near-term quantum devices by classical post-processing based on detector tomography

abstract: We propose a simple scheme to reduce readout errors in experiments on quantum systems with finite number of measurement outcomes. Our method relies on performing classical post-processing which is preceded by Quantum Detector Tomography, i.e., the reconstruction of a Positive-Operator Valued Measure (POVM) describing the given quantum measurement device. If the measurement device is affected only by an invertible classical noise, it is possible to correct the outcome statistics of future experiments performed on the same device. To support the practical applicability of this scheme for near-term quantum devices, we characterize measurements implemented in IBM's and Rigetti's quantum processors. We find that for these devices, based on superconducting transmon qubits, classical noise is indeed the dominant source of readout errors. Moreover, we analyze the influence of the presence of coherent errors and finite statistics on the performance of our error-mitigation procedure. Applying our scheme on the IBM's 5-qubit device, we observe a significant improvement of the results of a number of single- and two-qubit tasks including Quantum State Tomography (QST), Quantum Process Tomography (QPT), the implementation of non-projective measurements, and certain quantum algorithms (Grover's search and the Bernstein-Vazirani algorithm). Finally, we present results showing improvement for the implementation of certain probability distributions in the case of five qubits.


title: Disparate energy scaling of trajectory-dependent electronic excitations for slow protons and He ions

abstract: We have simultaneously measured angular distributions and electronic energy loss of helium ions and protons directly transmitted through self-supporting, single-crystalline silicon foils. We have compared the energy loss along channelled and random trajectories for incident ion energies between 50 keV and 200 keV. For all studied cases the energy loss in channelling geometry is found lower than in random geometry. In the case of protons, this difference increases with initial ion energy. This behaviour can be explained by the increasing contribution of excitations of core electrons, which are more likely to happen at small impact parameters accessible only in random geometry. For helium ions we observe a reverse trend - a decrease of the difference between channelled and random energy loss for increasing ion energy. Due to the inefficiency of core-electron excitations even at small impact parameters at such low energies, another mechanism has to be the cause for the observed difference. We provide evidence that the observation originates from reionisation events induced by close collisions of helium ions occurring only along random trajectories.


title: Data Augmentation for Instrument Classification Robust to Audio Effects

abstract: Reusing recorded sounds (sampling) is a key component in Electronic Music Production (EMP), which has been present since its early days and is at the core of genres like hip-hop or jungle. Commercial and non-commercial services allow users to obtain collections of sounds (sample packs) to reuse in their compositions. Automatic classification of one-shot instrumental sounds allows automatically categorising the sounds contained in these collections, allowing easier navigation and better characterisation. Automatic instrument classification has mostly targeted the classification of unprocessed isolated instrumental sounds or detecting predominant instruments in mixed music tracks. For this classification to be useful in audio databases for EMP, it has to be robust to the audio effects applied to unprocessed sounds. In this paper we evaluate how a state of the art model trained with a large dataset of one-shot instrumental sounds performs when classifying instruments processed with audio effects. In order to evaluate the robustness of the model, we use data augmentation with audio effects and evaluate how each effect influences the classification accuracy.


title: Hypersonic Bose-Einstein Condensates in Accelerator Rings

abstract: Some of the most sensitive and precise measurements to date are based on matterwave interferometry with freely falling atomic clouds. Examples include high-precision measurements of inertia, gravity and rotation. In order to achieve these very high sensitivities, the interrogation time has to be very long and consequently the experimental apparatus has to be very tall, in some cases reaching ten or even one hundred meters. Cancelling gravitational acceleration, e.g. in atomtronic circuits6,7 and matterwave guides, will result in compact devices having much extended interrogation times and thus much increased sensitivity both for fundamental and practical measurements. In this letter, we demonstrate extremely smooth and controllable matterwave guides by transporting Bose-Einstein condensates (BEC) over macroscopic distances: We use a novel neutral-atom accelerator ring to bring BECs to very high speeds (16x their velocity of sound) and transport them in a magnetic matterwave guide for 15 cm whilst fully preserving their internal coherence. The high angular momentum of more than 40000h per atom gives access to the higher Landau levels of quantum Hall states. The hypersonic velocities combined with our ability to control the potentials with pico-Kelvin precision open new perspectives in the study of superfluidity and give rise to new regimes of tunnelling and transport. Coherent matterwave guides are expected to enable interaction times of several seconds in highly compact devices. These developments will result in portable guided-atom interferometers for applications such as inertial navigation and gravity mapping.

categories: quant-ph cond-mat.supr-con physics.atom-ph physics.optics


title: A Vine-copula extension for the HAR model

abstract: The heterogeneous autoregressive (HAR) model is revised by modeling the joint distribution of the four partial-volatility terms therein involved. Namely, today's, yesterday's, last week's and last month's volatility components. The joint distribution relies on a (C-) Vine copula construction, allowing to conveniently extract volatility forecasts based on the conditional expectation of today's volatility given its past terms. The proposed empirical application involves more than seven years of high-frequency transaction prices for ten stocks and evaluates the in-sample, out-of-sample and one-step-ahead forecast performance of our model for daily realized-kernel measures. The model proposed in this paper is shown to outperform the HAR counterpart under different models for marginal distributions, copula construction methods, and forecasting settings.


title: Tunable coupling and isolation of single electrons in silicon metal-oxide-semiconductor quantum dots

abstract: Extremely long coherence times, excellent single-qubit gate fidelities and two-qubit logic have been demonstrated with silicon metal-oxide-semiconductor spin qubits, making it one of the leading platforms for quantum information processing. Despite this, a long-standing challenge in this system has been the demonstration of tunable tunnel coupling between single electrons. Here we overcome this hurdle with gate-defined quantum dots and show couplings that can be tuned on and off for quantum operations. We use charge sensing to discriminate between the (2,0) and (1,1) charge states of a double quantum dot and show excellent charge sensitivity. We demonstrate tunable coupling up to 13 GHz, obtained by fitting charge polarization lines, and tunable tunnel rates down to below 1 Hz, deduced from the random telegraph signal. The demonstration of tunable coupling between single electrons in a silicon metal-oxide-semiconductor device provides significant scope for high-fidelity two-qubit logic toward quantum information processing with standard manufacturing.


title: Primordial power spectra for curved inflating universes

abstract: Exact numerical primordial primordial power spectra are computed and plotted for the for the best-fit Planck 2018 curved universe parameters. It is found that the spectra have generic cut-offs and oscillations within the observable window for the level of curvature allowed by current CMB measurements and provide a better fit to current data. Derivations for the Mukhanov-Sasaki equation for curved universes are presented and analysed, and theoretical implications for the quantum and classical initial conditions for inflation are discussed within the curved regime.


title: Effect of black hole--plasma system on light beams

abstract: In the paper we discuss propagation of light around Kerr black hole surrounded by non-magnetized cold plasma with infinite conductivity. For that purpose, we use equations for propagation of light rays obtained within Synge's approach in the approximation of geometrical optics. We derive equation of deviation of a ray propagating close to the reference ray, which is a generalization of the well-known Jacobi equation, and use it to calculate the modification of angular distribution of stars observed close to a black hole surrounded by plasma, compared to the uniform star distribution that would be seen without black hole or plasma. We place the observer in the equatorial plane of the Kerr black hole and try various choices of plasma distributions described by mathematically simple formulae. Key features of star distribution on a local sky near the black hole are identified and the influence of plasma on them is discussed.


title: ASYNC: A Cloud Engine with Asynchrony and History for Distributed Machine Learning

abstract: ASYNC is a framework that supports the implementation of asynchrony and history for optimization methods on distributed computing platforms. The popularity of asynchronous optimization methods has increased in distributed machine learning. However, their applicability and practical experimentation on distributed systems are limited because current bulk-processing cloud engines do not provide a robust support for asynchrony and history. With introducing three main modules and bookkeeping system-specific and application parameters, ASYNC provides practitioners with a framework to implement asynchronous machine learning methods. To demonstrate ease-of-implementation in ASYNC, the synchronous and asynchronous variants of two well-known optimization methods, stochastic gradient descent and SAGA, are demonstrated in ASYNC.


title: Regularity results for a class of obstacle problems with $p,q-$growth conditions

abstract: In this paper we prove the local boundedness as well as the local Lipschitz continuity for solutions to a class of obstacle problems of the type $$\min\left\{\int_\Omega {F(x, Dz)}: z\in \mathcal{K}_{\psi}(\Omega)\right\}.$$ Here $\mathcal{K}_{\psi}(\Omega)$ is set of admissible functions $z \in W^{1,p}(\Omega)$ such that $z \ge \psi$ a.e. in $\Omega$, $\psi$ being the obstacle and $\Omega$ being an open bounded set of $\mathbb{R}^n$, $n \ge 2$. The main novelty here is that we are assuming $ F(x, Dz)$ satisfying $(p,q)$-growth conditions {and less restrictive assumptions on the obstacle with respect to the existing regularity results}.


title: Conditional emergence of classical domain and branching of quantum histories

abstract: We outline the Minimalistic Measurement Scheme (MMS) compatible with regular unitary evolution of a closed quantum system. Within this approach, a part of the system becomes informationally isolated (restricted) which leads to a natural emergence of the classical domain. This measurement scenario is a simpler alternative to environment-induced decoherence. In its basic version, MMS involves two ancilla qubits, $A$ and $X$, entangled with each other and with the System $S$. Informational or thermodynamic cost of measurement is represented by $X$-qubit being isolated, i.e. becoming unavailable for future interactions with the rest of the system. Conditional upon this isolation, $A$-qubit, that plays the role of an Apparatus, becomes classical and records the outcome of the measurement. The procedure may be used to perform von Neumann-style projective measurements or generalized ones, that corresponds to Positive-Operator Value Measure (POVM). By repeating the same generalized measurement multiple times with different $A$- and $X$-qubits, one asymptotically approaches the wave function collapse in the basis determined by the premeasurement process. We present a simple result for the total information extracted after $N$ such weak measurements. Building upon MMS, we propose a construction that maps a history of a quantum system onto a set of $A$-qubits. It resembles the Consistent History (CH) formulation of Quantum Mechanics (QM), but is distinct from it, and is built entirely within the conventional QM. In particular, consistency postulate of CH formalism is not automatically satisfied, but rather is an emerging property. Namely, each measurement event corresponds to the branching of mutually exclusive classical realities whose probabilities are additive. In a general case, however, the superposition between different histories is determined by the history density matrix.


title: A note on Machado-Bishop theorem in weighted spaces with applications

abstract: A unified version of Machado-Bishop theorem in weighted spaces is given. A number of applications illustrate its importance.


title: The Cherenkov Telescope Array

abstract: The Cherenkov Telescope Array (CTA) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. It will be capable of detecting gamma rays in the energy range from 20 GeV to more than 300 TeV with unprecedented precision in energy and directional reconstruction. With more than 100 telescopes of three different types it will be located in the northern hemisphere at La Palma, Spain, and in the southern at Paranal, Chile. CTA will be one of the largest astronomical infrastructures in the world with open data access and it will address questions in astronomy, astrophysics and fundamental physics in the next decades. In this presentation we will focus on the status of the CTA construction, the status of the telescope prototypes and highlight some of the physics perspectives.


title: An optimization-based cooperative path-following framework for multiple robotic vehicles

abstract: This paper addresses the design of an optimization-based cooperative path-following control law for multiple robotic vehicles that optimally balances the transient trade-off between coordination and path-following errors. To this end, we formulate a more general multi-agent framework where each agent is associated with (i) a continuous-time dynamical model, which governs the evolution of its state, and (ii) an output equation that is a function of both the state of the agent and a coordination vector. According to a given network topology, each agent can access its state and coordination vector, as well as the coordination vectors of the neighboring agents. In this setup, the goal is to design a distributed control law that steers the output signals to the origin, while simultaneously driving the coordination vectors of the agents of the network to consensus. To solve this, we propose a model predictive control scheme that builds on a pre-existing auxiliary consensus control law to design a performance index that combines the output regulation objective with the consensus objective. Convergence guarantees under which one can solve this coordinated output regulation problem are provided. Numerical simulations display the effectiveness of the proposed scheme applied to a cooperative path following control problem of a network of 3D nonholonomic robotic vehicles.


title: Multi-Granular Text Encoding for Self-Explaining Categorization

abstract: Self-explaining text categorization requires a classifier to make a prediction along with supporting evidence. A popular type of evidence is sub-sequences extracted from the input text which are sufficient for the classifier to make the prediction. In this work, we define multi-granular ngrams as basic units for explanation, and organize all ngrams into a hierarchical structure, so that shorter ngrams can be reused while computing longer ngrams. We leverage a tree-structured LSTM to learn a context-independent representation for each unit via parameter sharing. Experiments on medical disease classification show that our model is more accurate, efficient and compact than BiLSTM and CNN baselines. More importantly, our model can extract intuitive multi-granular evidence to support its predictions.


title: Generating fMRI volumes from T1-weighted volumes using 3D CycleGAN

abstract: Registration between an fMRI volume and a T1-weighted volume is challenging, since fMRI volumes contain geometric distortions. Here we present preliminary results showing that 3D CycleGAN can be used to synthesize fMRI volumes from T1-weighted volumes, and vice versa, which can facilitate registration.


title: Gravitational waves from neutron star mergers and their relation to the nuclear equation of state

abstract: In this article, I introduce ideas and techniques to extract information about the equation of state of matter at very high densities from gravitational waves emitted before, during and after the merger of binary neutron stars. I also review current work and results on the actual use of the first gravitational-wave observation of a neutron-star merger to set constraints on properties of such equation of state. In passing, I also touch on the possibility that what we observe in gravitational waves are not neutron stars, but something more exotic. In order to make this article more accessible, I also review the dynamics and gravitational-wave emission of neutron-star mergers in general, with focus on numerical simulations and on which representations of the equation of state are used for studies on binary systems.


title: End-to-end Learning for GMI Optimized Geometric Constellation Shape

abstract: Autoencoder-based geometric shaping is proposed that includes optimizing bit mappings. Up to 0.2 bits/QAM symbol gain in GMI is achieved for a variety of data rates and in the presence of transceiver impairments. The gains can be harvested with standard binary FEC at no cost w.r.t. conventional BICM.


title: Anomalous transport through algebraically localized states in one-dimension

abstract: Localization in one-dimensional disordered or quasiperiodic non-interacting systems in presence of power-law hopping is very different from localization in short-ranged systems. Power-law hopping leads to algebraic localization as opposed to exponential localization in short-ranged systems. Exponential localization is synonymous with insulating behavior in the thermodynamic limit. Here we show that the same is not true for algebraic localization. We show, on general grounds, that depending on the strength of the algebraic decay, the algebraically localized states can be actually either conducting or insulating in thermodynamic limit. We exemplify this statement with explicit calculations on the Aubry-Andr\'e-Harper model in presence of power-law hopping, with the power-law exponent $\alpha>1$, so that the thermodynamic limit is well-defined. We find a phase of this system where there is a mobility edge separating completely delocalized and algebraically localized states, with the algebraically localized states showing signatures of super-diffusive transport. Thus, in this phase, the mobility edge separates two kinds of conducting states, ballistic and super-diffusive. We trace the occurrence of this behavior to near-resonance conditions of the on-site energies that occur due to the quasi-periodic nature of the potential.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech quant-ph


title: An integral equation based numerical method for the forced heat equation on complex domains

abstract: Integral equation based numerical methods are directly applicable to homogeneous elliptic PDEs, and offer the ability to solve these with high accuracy and speed on complex domains. In this paper, extensions to problems with inhomogeneous source terms and time dependent PDEs, such as the heat equation, have been introduced. One such approach for the heat equation is to first discretise in time, and in each time-step solve a so-called modified Helmholtz equation with a parameter depending on the time step size. The modified Helmholtz equation is then split into two parts: a homogeneous part solved with a boundary integral method and a particular part, where the solution is obtained by evaluating a volume potential over the inhomogeneous source term over a simple domain. In this work, we introduce two components which are critical for the success of this approach: a method to efficiently compute a high-regularity extension of a function outside the domain where it is defined, and a special quadrature method to accurately evaluate singular and nearly singular integrals in the integral formulation of the modified Helmholtz equation for all time step sizes.


title: A survey on Crosstalk and Routing, Modulation Selection, Core and Spectrum Allocation in Elastic Optical Networks

abstract: Elastic Optical Networks (EON) emerge as a viable solution to supply the current growing demand for bandwidth. With the application of multi-core fibers (MCF) in EON links, it is possible to increase the availability of spectral resources. An EON network with MCF enables Space-Division Multiplexing (SDM), allowing the use of more resources in the fibers and increasing the capacity of attending circuit requests. However, the use of SDM brings some problems of interference between the circuits of a fiber, with greater emphasis on crosstalk interference. In this paper, some important concepts around EON are presented, along with the characterization of SDM technology. The Routing, Modulation, Spectrum and Core Allocation (RMSCA) problem is also characterized, and some solutions currently found in the literature are cited. After, the impact of crosstalk interference is discussed, and which elements are responsible for its occurrence. The paper is concludes with an evaluation of the state of the art, and the discrimination of the main points found from the study of papers related to the SDM-EON scenario.


title: An information-theoretic treatment of quantum dichotomies

abstract: Given two pairs of quantum states, we want to decide if there exists a quantum channel that transforms one pair into the other. The theory of quantum statistical comparison and quantum relative majorization provides necessary and sufficient conditions for such a transformation to exist, but such conditions are typically difficult to check in practice. Here, by building upon work by Matsumoto, we relax the problem by allowing for small errors in one of the transformations. In this way, a simple sufficient condition can be formulated in terms of one-shot relative entropies of the two pairs. In the asymptotic setting where we consider sequences of state pairs, under some mild convergence conditions, this implies that the quantum relative entropy is the only relevant quantity deciding when a pairwise state transformation is possible. More precisely, if the relative entropy of the initial state pair is strictly larger compared to the relative entropy of the target state pair, then a transformation with exponentially vanishing error is possible. On the other hand, if the relative entropy of the target state is strictly larger, then any such transformation will have an error converging exponentially to one. As an immediate consequence, we show that the rate at which pairs of states can be transformed into each other is given by the ratio of their relative entropies. We discuss applications to the resource theories of athermality and coherence.


title: Predicting Human Activities from User-Generated Content

abstract: The activities we do are linked to our interests, personality, political preferences, and decisions we make about the future. In this paper, we explore the task of predicting human activities from user-generated content. We collect a dataset containing instances of social media users writing about a range of everyday activities. We then use a state-of-the-art sentence embedding framework tailored to recognize the semantics of human activities and perform an automatic clustering of these activities. We train a neural network model to make predictions about which clusters contain activities that were performed by a given user based on the text of their previous posts and self-description. Additionally, we explore the degree to which incorporating inferred user traits into our model helps with this prediction task.


title: Resilient Sensor Architecture Design and Tradespace Analysis for Autonomous Vehicle Localization and Mapping

abstract: As autonomous cars are rolled out into new environments, their ability to solve the simultaneous localization and mapping (SLAM) problem becomes critical. In order to tackle this problem, autonomous vehicles rely on sensor suites that provide them with information about their operating environment. When large scale production is taken into consideration, a trade-off between an acceptable sensor suite cost and its resulting performance characteristics arises. Furthermore, guaranteeing the system's performance requires a resilient sensor network design. This work seeks to address such trade-offs by introducing a method that takes into account the performance, cost, and resiliency of distinct sensor selections. As a result, this method is able to offer sensor combination recommendations based on the vehicle's operating environment. It is found that the structure of the environment influences sensor placement, and that the design of a resilient sensor network involves careful consideration of both environmental attributes such as landmark density and location, as well as the available types of complimentary sensors. Demonstration of the proposed approach is shown by evaluating it using sequences from the KITTI Benchmark Suite.


title: Numerical Methods for Simulating Star Formation

abstract: We review the numerical techniques for ideal and non-ideal magneto-hydrodynamics (MHD) used in the context of star formation simulations. We outline the specific challenges offered by modeling star forming environments, which are dominated by supersonic and super-Alfvenic turbulence in a radiative, self-gravitating fluid. These conditions are rather unique in physics and engineering and pose particularly severe restrictions on the robustness and accuracy of numerical codes. One striking aspect is the formation of collapsing fluid elements leading to the formation of singularities that represent point-like objects, namely the proto-stars. Although a few studies have attempted to resolve the formation of the first and second Larson cores, resolution limitations force us to use sink particle techniques, with sub-grid models to compute the accretion rates of mass, momentum and energy, as well as their ejection rate due to radiation and jets from the proto-stars. We discuss the most popular discretisation techniques used in the community, namely smoothed particle hydrodynamics, finite difference and finite volume methods, stressing the importance to maintain a divergence-free magnetic field. We discuss how to estimate the truncation error of a given numerical scheme, and its importance in setting the magnitude of the numerical diffusion. This can have a strong impact on the outcome of these MHD simulations, where both viscosity and resistivity are implemented at the grid scale. We then present various numerical techniques to model non-ideal MHD effects, such as Ohmic and ambipolar diffusion, as well as the Hall effect. These important physical ingredients are posing strong challenges in term of resolution and time stepping. For the latter, several strategies are discussed to overcome the limitations due to prohibitively small time steps (abridged).


title: Localization Principle of the Spectral Expansions of Distributions Connected with Schrodinger Operator

abstract: In this paper the localization properties of the spectral expansions of distributions related to the self adjoint extension of the Schrodinger operator are investigated. Spectral decompositions of the distributions and some classes of distributions are defined. Estimations for Riesz means of the spectral decompositions of the distributions in the norm of the Sobolev classes with negative order are obtained.


title: Post-synaptic potential regularization has potential

abstract: Improving generalization is one of the main challenges for training deep neural networks on classification tasks. In particular, a number of techniques have been proposed, aiming to boost the performance on unseen data: from standard data augmentation techniques to the $\ell_2$ regularization, dropout, batch normalization, entropy-driven SGD and many more.\\ In this work we propose an elegant, simple and principled approach: post-synaptic potential regularization (PSP). We tested this regularization on a number of different state-of-the-art scenarios. Empirical results show that PSP achieves a classification error comparable to more sophisticated learning strategies in the MNIST scenario, while improves the generalization compared to $\ell_2$ regularization in deep architectures trained on CIFAR-10.


title: Positively Hyperbolic Varieties, Tropicalization, and Positroids

abstract: A variety of codimension $c$ in complex affine space is called positively hyperbolic if the imaginary part of any point in it does not lie in any positive linear subspace of dimension $c$. Positively hyperbolic hypersurfaces are defined by stable polynomials. We give a new characterization of positively hyperbolic varieties using sign variations, and show that they are equivalently defined by being hyperbolic with respect to the positive part of the Grassmannian, in the sense of Shamovich and Vinnikov. We prove that positively hyperbolic projective varieties have tropicalizations that are locally subfans of the type $A$ hyperplane arrangement defined by $x_i = x_j$, in which the maximal cones satisfy a non-crossing condition. This gives new proofs of some results of Choe--Oxley--Sokal--Wagner and Br\"and\'en on Newton polytopes and tropicalizations of stable polynomials. We settle the question of which tropical varieties can be obtained as tropicalizations of positively hyperbolic varieties in the case of tropical toric varieties, constant-coefficient tropical curves, and Bergman fans. Along the way, we also give a new characterization of positroids in terms of a non-crossing condition on their Bergman fans.


title: The Dynamical Diquark Model: Fine Structure and Isospin

abstract: We incorporate fine-structure corrections into the dynamical diquark model of multiquark exotic hadrons. These improvements include effects due to finite diquark size, spin-spin couplings within the diquarks, and most significantly, isospin-dependent couplings in the form of pionlike exchanges assumed to occur between the light quarks within the diquarks. Using a simplified two-parameter interaction Hamiltonian, we obtain fits in which the isoscalar $J^{PC} = 1^{++}$ state---identified as the $X(3872)$---appears naturally as the lightest exotic (including all states that are predicted by the model but have not yet been observed), while the $Z_c(3900)$ and $Z_c(4020)$ decay predominantly to $J/\psi$ and $\eta_c$, respectively, in accord with experiment. We explore implications of this model for the excited tetraquark multiplets and the pentaquarks.


title: Topological Weyl-like Phonons and Nodal Line Phonons in Graphene

abstract: By means of first-principles calculations and modeling analysis, we have predicted that the traditional 2D-graphene hosts the topological phononic Weyl-like points (PWs) and phononic nodal line (PNL) in its phonon spectrum. The phonon dispersion of graphene hosts three type-I PWs (both PW1 and PW2 at the BZ corners \emph{K} and \emph{K}', and PW3 locating along the $\Gamma$-\emph{K} line), one type-II PW4 locating along the $\Gamma$-\emph{M} line, and one PNL surrounding the centered $\Gamma$ point in the $q_{x,y}$ plane. The calculations further reveal that Berry curvatures are vanishingly zero throughout the whole BZ, except for the positions of these four pairs of Weyl-like phonons, at which the non-zero singular Berry curvatures appear with the Berry phase of $\pi$ or -$\pi$, confirming its topological non-trivial nature. The topologically protected non-trivial phononic edge states have been also evidenced along both the zigzag-edged and armchair-edged boundaries. These results would pave the ways for further studies of topological phononic properties of graphene, such as phononic destructive interference with a suppression of backscattering and intrinsic phononic quantum Hall-like effects.


title: Some new block designs of dimension three

abstract: The dimension of a block design is the maximum positive integer $d$ such that any $d$ of its points are contained in a proper subdesign. Pairwise balanced designs PBD$(v,K)$ have dimension at least two as long as not all points are on the same line. On the other hand, designs of dimension three appear to be very scarce. We study designs of dimension three with block sizes in $K=\{3,4\}$ or $\{3,5\}$, obtaining several explicit constructions and one nonexistence result in the latter case. As applications, we obtain a result on dimension three triple systems having arbitrary index as well as symmetric latin squares which are covered in a similar sense by proper subsquares.