title: A guided network propagation approach to identify disease genes that combines prior and new information

abstract: A major challenge in biomedical data science is to identify the causal genes underlying complex genetic diseases. Despite the massive influx of genome sequencing data, identifying disease-relevant genes remains difficult as individuals with the same disease may share very few, if any, genetic variants. Protein-protein interaction networks provide a means to tackle this heterogeneity, as genes causing the same disease tend to be proximal within networks. Previously, network propagation approaches have spread signal across the network from either known disease genes or genes that are newly putatively implicated in the disease (e.g., found to be mutated in exome studies or linked via genome-wide association studies). Here we introduce a general framework that considers both sources of data within a network context. Specifically, we use prior knowledge of disease-associated genes to guide random walks initiated from genes that are newly identified as perhaps disease-relevant. In large-scale testing across 24 cancer types, we demonstrate that our approach for integrating both prior and new information not only better identifies cancer driver genes than using either source of information alone but also readily outperforms other state-of-the-art network-based approaches. To demonstrate the versatility of our approach, we also apply it to genome-wide association data to identify genes functionally relevant for several complex diseases. Overall, our work suggests that guided network propagation approaches that utilize both prior and new data are a powerful means to identify disease genes.


title: Robust Traffic Control Using a First Order Macroscopic Traffic Flow Model

abstract: Traffic control is at the core of research in transportation engineering because it is one of the best practices for reducing traffic congestion. It has been shown in recent years that the traffic control problem involving Lighthill-Whitham-Richards (LWR) model can be formulated as a Linear Programming (LP) problem given that the corresponding initial conditions and the model parameters in the fundamental diagram are fixed. However, the initial conditions can be uncertain when studying actual control problems. This paper presents a stochastic programming formulation of the boundary control problem involving chance constraints, to capture the uncertainty in the initial conditions. Different objective functions are explored using this framework, and case studies for both a single highway link and a small network are conducted. In addition, the optimal results are validated with Monte Carlo simulation.


title: Graph Inference Learning for Semi-supervised Classification

abstract: In this work, we address semi-supervised classification of graph data, where the categories of those unlabeled nodes are inferred from labeled nodes as well as graph structures. Recent works often solve this problem via advanced graph convolution in a conventionally supervised manner, but the performance could degrade significantly when labeled data is scarce. To this end, we propose a Graph Inference Learning (GIL) framework to boost the performance of semi-supervised node classification by learning the inference of node labels on graph topology. To bridge the connection between two nodes, we formally define a structure relation by encapsulating node attributes, between-node paths, and local topological structures together, which can make the inference conveniently deduced from one node to another node. For learning the inference process, we further introduce meta-optimization on structure relations from training nodes to validation nodes, such that the learnt graph inference capability can be better self-adapted to testing nodes. Comprehensive evaluations on four benchmark datasets (including Cora, Citeseer, Pubmed, and NELL) demonstrate the superiority of our proposed GIL when compared against state-of-the-art methods on the semi-supervised node classification task.


title: Sequent calculi for a unity of logic

abstract: We present a novel unity of logic, viz., a single sequent calculus that embodies classical, intuitionistic and linear logics. Concretely, we define classical linear logic negative (CLL$^-$), a new logic that is classical and linear yet discards the polarities and the strict De Morgan laws in classical linear logic (CLL). Then, we define unlinearisation and classicalisation on sequent calculi such that unlinearisation maps CLL$^-$ (resp. intuitionistic linear logic (ILL)) to classical logic (CL) (resp. intuitionistic logic (IL)), and classicalisation maps IL (resp. ILL) to CL (resp. CLL$^-$) modulo conservative extensions. By these two maps, only a sequent calculus for a conservative extension of ILL suffices for ILL, IL, CLL$^-$ and CL. This result achieves a simple, highly systematic unity of logic by discarding the polarities and the strict De Morgan laws, which (arguably) only CLL has, and consisting of the uniform classicalisation and unlinearisation, which commute. Previous methods do not satisfy these points. Our unity also clarifies the dichotomies between intuitionisity and classicality, and between linearity and non-linearity of logic, which are symmetric.


title: FRaZ: A Generic High-Fidelity Fixed-Ratio Lossy Compression Framework for Scientific Floating-point Data

abstract: With ever-increasing volumes of scientific floating-point data being produced by high-performance computing applications, significantly reducing scientific floating-point data size is critical, and error-controlled lossy compressors have been developed for years. None of the existing scientific floating-point lossy data compressors, however, support effective fixed-ratio lossy compression. Yet fixed-ratio lossy compression for scientific floating-point data not only compresses to the requested ratio but also respects a user-specified error bound with higher fidelity. In this paper, we present FRaZ: a generic fixed-ratio lossy compression framework respecting user-specified error constraints. The contribution is twofold. (1) We develop an efficient iterative approach to accurately determine the appropriate error settings for different lossy compressors based on target compression ratios. (2) We perform a thorough performance and accuracy evaluation for our proposed fixed-ratio compression framework with multiple state-of-the-art error-controlled lossy compressors, using several real-world scientific floating-point datasets from different domains. Experiments show that FRaZ effectively identifies the optimum error setting in the entire error setting space of any given lossy compressor. While fixed-ratio lossy compression is slower than fixed-error compression, it provides an important new lossy compression technique for users of very large scientific floating-point datasets.


title: Higher order curvature flows of plane curves with generalised Neumann boundary conditions

abstract: We consider the parabolic polyharmonic diffusion and $L^2$-gradient flows of the $m$-th arclength derivative of curvature for regular closed curves evolving with generalised Neumann boundary conditions. In the polyharmonic case, we prove that if the curvature of the initial curve is small in $L^2$, then the evolving curve converges exponentially in the $C^\infty$ topology to a straight horizontal line segment. The same behaviour is shown for the $L^2$-gradient flow provided the energy of the initial curve is sufficiently small. In each case the smallness conditions depend only on $m$.


title: The radial acceleration relation in galaxy clusters

abstract: Recently, the discovery of the radial acceleration relation (RAR) in galaxies has been regarded as an indirect support of alternative theories of gravity such as Modified Newtonian Dynamics (MOND) and modified gravity. This relation indicates a tight correlation between dynamical mass and baryonic mass in galaxies with different sizes and morphology. However, if the RAR relation is scale-independent and could be explained by alternative theories of gravity, this relation should be universal and true for galaxy clusters as well. In this article, by using the x-ray data of a sample of galaxy clusters, we investigate if there exists any tight correlation between dynamical mass and baryonic mass in galaxy clusters, assuming hot gas mass distribution almost representing baryonic distribution and that the galaxy clusters are virialized. We show that the resulting RAR of 52 non-cool-core galaxy clusters scatters in a large parameter space, possibly due to our simplifying assumptions and unclear matter content in galaxy clusters. This might indicate that the RAR is unlikely to be universal and scale-independent.


title: Demonstration of second harmonic generationin gallium phosphide nano-waveguides

abstract: We designed, fabricated and tested gallium phosphide (GaP) nano-waveguides for second harmonic generation (SHG). We demonstrate SHG in the visible range around 655 nm using low power continuous-wave pump in the optical communication O-band. Our structures utilize modal phase matching, such that lower order eigenmodes of the pump are phase matched to higher order eigenmodes of the second harmonic. We observe phase matched SHG for different combinations of interacting modes by varying the widths of the waveguides and tuning the wavelength of the pump. The presented results contribute to the development of integrated photonic platforms with efficient nonlinear wave-mixing processes for classical and quantum applications.


title: On the weak Lefschetz property for almost complete intersections generated by uniform powers of general linear forms

abstract: In 2012, Migliore, the first author, and Nagel conjectured that, for all $n\geq 4$, the artinian ideal $I=(L_0^d,\ldots,L_{2n+1}^d) \subset R=k[x_0,\ldots,x_{2n}]$ generated by the $d$-th powers of $2n+2$ general linear forms fails to have the weak Lefschetz property if and only if $d>1$. This paper is entirely devoted to prove partially this conjecture. More precisely, we prove that $R/I$ fails to have the weak Lefschetz property, provided $4\leq n\leq 8,\ d\geq 4$ or $d=2r,\ 1\leq r\leq 8,\ 4\leq n\leq 2r(r+2)-1$.


title: Learning to Augment Expressions for Few-shot Fine-grained Facial Expression Recognition

abstract: Affective computing and cognitive theory are widely used in modern human-computer interaction scenarios. Human faces, as the most prominent and easily accessible features, have attracted great attention from researchers. Since humans have rich emotions and developed musculature, there exist a lot of fine-grained expressions in real-world applications. However, it is extremely time-consuming to collect and annotate a large number of facial images, of which may even require psychologists to correctly categorize them. To the best of our knowledge, the existing expression datasets are only limited to several basic facial expressions, which are not sufficient to support our ambitions in developing successful human-computer interaction systems. To this end, a novel Fine-grained Facial Expression Database - F2ED is contributed in this paper, and it includes more than 200k images with 54 facial expressions from 119 persons. Considering the phenomenon of uneven data distribution and lack of samples is common in real-world scenarios, we further evaluate several tasks of few-shot expression learning by virtue of our F2ED, which are to recognize the facial expressions given only few training instances. These tasks mimic human performance to learn robust and general representation from few examples. To address such few-shot tasks, we propose a unified task-driven framework - Compositional Generative Adversarial Network (Comp-GAN) learning to synthesize facial images and thus augmenting the instances of few-shot expression classes. Extensive experiments are conducted on F2ED and existing facial expression datasets, i.e., JAFFE and FER2013, to validate the efficacy of our F2ED in pre-training facial expression recognition network and the effectiveness of our proposed approach Comp-GAN to improve the performance of few-shot recognition tasks.


title: A Derivative-Free Method for Solving Elliptic Partial Differential Equations with Deep Neural Networks

abstract: We introduce a deep neural network based method for solving a class of elliptic partial differential equations. We approximate the solution of the PDE with a deep neural network which is trained under the guidance of a probabilistic representation of the PDE in the spirit of the Feynman-Kac formula. The solution is given by an expectation of a martingale process driven by a Brownian motion. As Brownian walkers explore the domain, the deep neural network is iteratively trained using a form of reinforcement learning. Our method is a 'Derivative-Free Loss Method' since it does not require the explicit calculation of the derivatives of the neural network with respect to the input neurons in order to compute the training loss. The advantages of our method are showcased in a series of test problems: a corner singularity problem, an interface problem, and an application to a chemotaxis population model.


title: Wettability and surface energy of parylene F

abstract: Parylenes are barrier materials employed as protective layers. However, many parylenes are unsuitable for applications under harsh conditions. A new material, parylene F, demonstrates considerable potential for a wide range of applications due to its high temperature and UV resistance. For the first time, the wettability and surface energy of parylene F were investigated to determine the feasibility of parylene F as an alternative to the commonly employed parylene C. The results show that parylene F has a hydrophobic surface with a water contact angle of 109.63 degrees. We found that 3.5 ul probe liquid is an optimal value for the contact angle measurement of parylene F. Moreover, we found that the Owens-Wendt-Kaelble and the Lifshitz-van der Waals/acid-base approaches are unsuitable for determining the surface energy of parylene F, whereas an approach based on the limitless liquid-solid interface wetting system is compatible. Furthermore, the results show that parylene F has a surface energy of 39.05 mJ/m2. Considering the improved resistance, relatively low cost, and the desirable properties, parylene F can replace parylene C for applications under harsh conditions.


title: Enhanced pair production in frequency modulated Sauter potential wells

abstract: Electron-positron pair production in frequency modulated Sauter potential wells is investigated in the framework of the computational quantum field theory. In combined potential wells with a static Sauter potential well and a frequency modulated oscillating one, the modulated amplitude has a large effect on the number of created pairs. The optimal modulation amplitude of frequency at different center frequencies is obtained, which increases the number of electrons at about two times. However, for a single oscillating potential well with frequency modulation, chirp effect is sensitive to the center frequency, and the number of electrons can be enhanced even to four orders of magnitude at a regime of low center frequency. It implies that for a slowly oscillating Sauter potential well, the chirp effect through the frequency modulation is better than adding a static potential well to improve the pair production.


title: Bouncing compact objects I: Quantum extension of the Oppenheimer-Snyder collapse

abstract: This article proposes a generalization of the Oppenheimer-Snyder model which describes a bouncing compact object. The corrections responsible for the bounce are parameterized in a general way so as to remain agnostic about the specific mechanism of singularity resolution at play. It thus develops an effective theory based on a thin shell approach, inferring generic properties of such a UV complete gravitational collapse. The main result comes in the form of a strong constraint applicable to general UV models : if the dynamics of the collapsing star exhibits a bounce, it always occurs below, or at most at the energy threshold of horizon formation, so that only an instantaneous trapping horizon may be formed while a trapped region never forms. This conclusion relies solely on i) the assumption of continuity of the induced metric across the time-like surface of the star and ii) the assumption of a classical Schwarzschild geometry describing the (vacuum) exterior of the star. In particular, it is completely independent of the choice of corrections inside the star which leads to singularity-resolution. The present model provides thus a general framework to discuss bouncing compact objects, for which the interior geometry is modeled either by a classical or a quantum bounce. In the later case, our no-go result regarding the formation of trapped region suggests that additional structure, such as the formation of an inner horizon, is needed to build consistent models of matter collapse describing black-to-white hole bounces. Indeed, such additional structure is needed to keep quantum gravity effects confined to the high curvature regime, in the deep interior region, providing thus a new challenge for current constructions of quantum black-to-white hole bounce models.


title: Nonlinearity in the Lorentz Oscillator Model

abstract: The magnetic force is retained in the Lorentz Oscillator Model and a perturbation solution is derived beyond the dipole approximation. Perturbation series for the electric dipole, magnetic dipole, and electric quadrupole moments in addition to their macroscopic, spatially-averaged counterparts are also obtained. These expressions are shown to describe many optical phenomena: the normal, linear Zeeman, Faraday, and magneto-optic Kerr effects, the linear Stark, Pockels, and Kerr effects (when the restoring force is anharmonic), the momentum and angular momentum of a photon in matter (contributing to the Abraham-Minkowski debate), the photon drag and inverse Faraday effects, and sum-frequency and second-harmonic generation in centrosymmetric media. The qualitative and quantitative accuracy is assessed by comparison to both experiment and more fundamental theory (e.g. quantum mechanics). New insights into the relationship between the Stark and Pockels/Kerr effects also emerge.


title: Implication Zroupoids and Birkhoff Systems

abstract: An algebra $A = \langle A, \to, 0 \rangle$, where $\to$ is binary and $0$ is a constant, is called an implication zroupoid (I-zroupoid, for short) if A satisfies the identities: $(x \to y) \to z \approx ((z' \to x) \to (y \to z)')'$, where $x' := x \to 0$, and $0'' \approx 0$. These algebras generalize De Morgan algebras and $\lor$-semilattices with zero. Let I denote the variety of implication zroupoids. For details on the motivation leading to these algebras, we refer the reader to [San12] (or the relevant papers mentioned at the end of this paper). The investigations into the structure of the lattice of subvarieties of I, begun in [San12], have continued in [CS16a, CS16b, CS17a, CS17b, CS18a, CS18b, CS19] and [GSV19]. The present paper is a sequel to this series of papers and is devoted to making further contributions to the theory of implication zroupoids. The identity (BR): $x \land (x \lor y) \approx x \lor (x \land y)$ is called the Birkhoff's identity. The main purpose of this paper is to prove that if A is an algebra in the variety I, then the derived algebra $A_{mj} := \langle A; \land, \lor \rangle$, where $a \land b := (a \to b')'$ and $a \lor b := (a' \land b')'$, satisfies the Birkhoff's identity. As a consequence, we characterize the implication zroupoids A whose derived algebras $A_{mj}$ are Birkhoff systems. It also follows from the main result that there are bisemigroups that are not bisemilattices but satisfy the Birkhoff's identity, which suggests a more general notion, than Birkhoff systems, of "Birkhoff bisemigroups" as bisemigroups satisfying the Birkhoff's identity. The paper concludes with an open problem on Birkhoff bisemigroups.


title: Interpreting Galaxy Deblender GAN from the Discriminator's Perspective

abstract: Generative adversarial networks (GANs) are well known for their unsupervised learning capabilities. A recent success in the field of astronomy is deblending two overlapping galaxy images via a branched GAN model. However, it remains a significant challenge to comprehend how the network works, which is particularly difficult for non-expert users. This research focuses on behaviors of one of the network's major components, the Discriminator, which plays a vital role but is often overlooked, Specifically, we enhance the Layer-wise Relevance Propagation (LRP) scheme to generate a heatmap-based visualization. We call this technique Polarized-LRP and it consists of two parts i.e. positive contribution heatmaps for ground truth images and negative contribution heatmaps for generated images. Using the Galaxy Zoo dataset we demonstrate that our method clearly reveals attention areas of the Discriminator when differentiating generated galaxy images from ground truth images. To connect the Discriminator's impact on the Generator, we visualize the gradual changes of the Generator across the training process. An interesting result we have achieved there is the detection of a problematic data augmentation procedure that would else have remained hidden. We find that our proposed method serves as a useful visual analytical tool for a deeper understanding of GAN models.


title: Interpolation and Sampling in Analytic Tent Spaces

abstract: Introduced by Coifman, Meyer, and Stein, the tent spaces have seen wide applications in Harmonic Analysis. Their analytic cousins have seen some applications involving the derivatives of Hardy space functions. Moreover, the tent spaces have been a recent focus of research. We introduce the concept of interpolating and sampling sequences for analytic tent spaces analogously to the same concepts for Bergman spaces. We then characterize such sequences in terms of Seip's upper and lower uniform density. We accomplish this by exploiting a kind of M\"obius invariance for the tent spaces.


title: Bouncing compact objects. Part II: Effective theory of a pulsating Planck star

abstract: This article presents an effective quantum extension of the seminal Oppenheimer-Snyder (OS) collapse in which the singularity resolution is modeled using the effective dynamics of the spatially closed loop quantum cosmology. Imposing the minimal junction conditions, namely the Israel-Darmois conditions, we glue this bouncing LQC geometry to the classical vacuum exterior Schwarzschild geometry across a time-like thin-shell. Consistency of the construction leads to several major deviations from the classical OS collapse model. Firstly, no trapped region can form and the bounce occurs always above, or at most at the Schwarzschild radius. Secondly, the bouncing star discussed here admits an IR cut-off, additionally to the UV cut-off and corresponds therefore to a pulsating compact object. Thirdly, the scale at which quantum gravity effects become non-negligible is encoded in the ratio between the UV cut-off of the quantum theory and the IR cut-off, which in turn, encodes the minimal energy density $\rho_{\text{min}}$ of the star prior to collapse. This energy density is no more fixed by the mass and maximal radius as in the classical OS model, but is now a free parameter of the model. In the end, while the present model cannot describe a black-to-white hole bounce as initially suggested by the Planck star model, it provides a concrete realization of a pulsating compact object based on LQC techniques. Consistency of the model shows that its regime of applicability is restricted to Planckian relics while macroscopic stellar objects are excluded. This first minimal construction should serve as a platform for further investigations in order to explore the physics of bouncing compact objects within the framework of loop quantum cosmology.


title: Quantum decoherence by Coulomb interaction

abstract: The performance of modern quantum devices in communication, metrology or microscopy relies on the quantum-classical interaction which is generally described by the theory of decoherence. Despite the high relevance for long coherence times in quantum electronics, decoherence mechanisms mediated by the Coulomb force are not well understood yet and several competing theoretical models exist. Here, we present an experimental study of the Coulomb-induced decoherence of free electrons in a superposition state in a biprism electron interferometer close to a semiconducting and metallic surface. The decoherence was determined through a contrast loss at different beam path separations, surface distances and conductibilities. To clarify the current literature discussion, four theoretical models were compared to our data. We could rule out three of them and got good agreement with a theory based on macroscopic quantum electrodynamics. The results will enable the determination and minimization of specific decoherence channels in the design of novel quantum instruments.


title: Positively curved K\"ahler metrics on tube domains and their applications to optimal transport

abstract: In this article, we study a class of K\"ahler manifolds defined on tube domains in $\mathbb{C}^n$, and in particular those which have $O(n) \times \mathbb{R}^n$ symmetry. For these, we prove a uniqueness result showing that any such manifold which is complete and has non-negative orthogonal bisectional curvature ($n \geq 3$) or non-negative bisectional curvature ($n \geq 2$) is biholomorphically isometric to $\mathbb{C}^n$. We also consider another curvature tensor called the \emph{orthogonal anti-bisectional curvature}. We find necessary and sufficient conditions for a complete $O(n)$-symmetric tube domain to have non-negative orthogonal anti-bisectional curvature and provide several examples of complete metrics which satisfy this condition. Finally, we discuss some applications of these spaces within optimal transport. In particular, we study "synthetic" curvature bounds for non-smooth geometries and how they can be applied to the rough geometry induced by the Monge cost $c(x,y)=\|x-y\|$.


title: A Reliable Gravity Compensation Control Strategy for dVRK Robotic Arms With Nonlinear Disturbance Forces

abstract: External disturbance forces caused by nonlinear springy electrical cables in the Master Tool Manipulator (MTM) of the da Vinci Research Kit (dVRK) limits the usage of the existing gravity compensation methods. Significant motion drifts at the MTM tip are often observed when the MTM is located far from its identification trajectory, preventing the usage of these methods for the entire workspace reliably. In this paper, we propose a general and systematic framework to address the problems of the gravity compensation for the MTM of the dVRK. Particularly, high order polynomial models were used to capture the highly nonlinear disturbance forces and integrated with the Multi-step Least Square Estimation (MLSE) framework. This method allows us to identify the parameters of both the gravitational and disturbance forces for each link sequentially, preventing residual error passing among the links of the MTM with uneven mass distribution. A corresponding gravity compensation controller was developed to compensate the gravitational and disturbance forces. The method was validated with extensive experiments in the majority of the manipulator's workspace, showing significant performance enhancements over existing methods. Finally, a deliverable software package in MATLAB and C++ was integrated with dVRK and published in the dVRK community for open-source research and development.


title: Gagliardo-Nirenberg, Trudinger-Moser and Morrey inequalities on Dirichlet spaces

abstract: With a view towards Riemannian or sub-Riemannian manifolds, RCD metric spaces and specially fractals, this paper makes a step further in the development of a theory of heat semigroup based $(1,p)$ Sobolev spaces in the general framework of Dirichlet spaces. Under suitable assumptions that are verified in a variety of settings, the tools developed by D. Bakry, T. Coulhon, M. Ledoux and L. Saloff-Coste in the paper "Sobolev inequalities in disguise" allow us to obtain the whole family of Gagliardo-Nirenberg and Trudinger-Moser inequalities with optimal exponents. The latter depend not only on the Hausdorff and walk dimensions of the space but also on other invariants. In addition, we prove Morrey type inequalities and apply them to study the infimum of the exponents that ensure continuity of Sobolev functions. The results are illustrated for fractals using the Vicsek set, whereas several conjectures are made for nested fractals and the Sierpinski carpet.


title: On the sets of lengths of Puiseux monoids generated by multiple geometric sequences

abstract: In this paper, we study some of the factorization aspects of rational multicyclic monoids, that is, additive submonoids of the nonnegative rational numbers generated by multiple geometric sequences. In particular, we provide a complete description of the rational multicyclic monoids $M$ that are hereditarily atomic (i.e., every submonoid of $M$ is atomic). Additionally, we show that the sets of lengths of certain rational multicyclic monoids are finite unions of multidimensional arithmetic progressions, while their unions satisfy the Structure Theorem for Unions of Sets of Lengths. Finally, we realize arithmetic progressions as the sets of distances of some additive submonoids of the nonnegative rational numbers.


title: A New Fairness Model based on User's Objective for Multi-user Multi-processor Online Scheduling

abstract: Resources of a multi-user system in multi-processor online scheduling are shared by competing users in which fairness is a major performance criterion for resource allocation. Fairness ensures equality in resource sharing among the users. According to our knowledge, fairness based on the user's objective has neither been comprehensively studied nor a formal fairness model has been well defined in the literature. This motivates us to explore and define a new model to ensure algorithmic fairness with quantitative performance measures based on optimization of the user's objective. In this paper, we propose a new model for fairness in Multi-user Multi-processor Online Scheduling Problem(MUMPOSP). We introduce and formally define quantitative fairness measures based on user's objective by optimizing makespan for individual user in our proposed fairness model. We also define the unfairness of deprived users and absolute fairness of an algorithm. We obtain lower bound results for the absolute fairness for m identical machines with equal length jobs. We show that our proposed fairness model can serve as a framework for measuring algorithmic fairness by considering various optimality criteria such as flow time and sum of completion times.


title: Selective C-C Coupling by Spatially Confined Dimeric Metal Centers

abstract: Direct conversion of carbon dioxide (CO2) to high-energy fuels and high-value chemicals is a fascinating sustainable strategy. For most of the current electrocatalysts for CO2 reduction, however, multi-carbon products are inhibited by large overpotentials and low selectivity. For practical applications, there remains a big gap of knowledge in proper manipulation of the C-C coupling process. Herein, we exploit dispersed 3d transition metal dimers as spatially confined dual reaction centers for selective reduction of CO2 to liquid fuels. Various nitrogenated holey carbon monolayers are shown to be promising templates to stabilize these metal dimers and dictate their electronic structures, allowing precise control of the catalytic activity and product selectivity. By comprehensive first-principles calculations, we screen the suitable transition metal dimers that universally have high activity for ethanol (C2H5OH). Furthermore, remarkable selectivity for C2H5OH against other C1 and C2 products is found for Fe2 dimer anchored on C2N monolayer. The correlation between the activity and d band center of the supported metal dimer as well as the role of electronic coupling between the metal dimer and the carbon substrates are thoroughly elucidated.


title: Prediction for structure stability and ultrahigh hydrogen evolution performance of monolayer 2H-CrS2

abstract: By a combination of the first-principles calculations and climbing image nudged elastic band method (ciNEB) we investigate structure stabilities and hydrogen evolution reaction (HER) performance of monolayer 2H-CrS2. The results suggest the free energy for the Volmer reaction in the monolayer 2H-CrS2 with S vacancy is 0.07 eV, comparable with Pt-based catalyst, and HER on the surface of the monolayer is prone to the Volmer-Heyrovsky mechanism with no energy barrier. We propose that high HER performance stems from the reduction of the energy level of d-band center. Additionally, the S vacancy leads to defect states in the middle of electronic bandgap and the reduction of potential barrier between the S atom layer and the vacuum, which is conducive to improve HER performance.


title: On Integer Sequences Associated To Two Distinct Sums

abstract: In this note, we show the existence of integer sequences of lengths at least 3 (except 7) such that for every integer in position $i\equiv 1\pmod{4}$ (respectively position $j\equiv 3\pmod{4}$), counting from left to right, the sum of the integer and the adjacent integer(s) has a constant sum $x$ (respectively $y$) with $x\ne y$.


title: Proceedings 16th Workshop on Quantitative Aspects of Programming Languages and Systems

abstract: This EPTCS volume contains the proceedings of the 16th Workshop on Quantitative Aspects of Programming Languages and Systems (QAPL 2019) held in Prague, Czech Republic, on Sunday 7 April 2019. QAPL 2019 was a satellite event of the European Joint Conferences on Theory and Practice of Software (ETAPS 2019). QAPL focuses on quantitative aspects of computations, which may refer to the use of physical quantities (time, bandwidth, etc.) as well as mathematical quantities (e.g., probabilities) for the characterisation of the behaviour and for determining the properties of systems. Such quantities play a central role in defining both the model of systems (architecture, language design, semantics) and the methodologies and tools for the analysis and verification of system properties. The aim of the QAPL workshop series is to discuss the explicit use of time and probability and general quantities either directly in the model or as a tool for the analysis or synthesis of systems. The 16th edition of QAPL also focuses on discussing the developments, challenges and results in this area covered by our workshop in its nearly 20-year history.


title: New Normal Forms For Degree Three Polynomials and Rational Functions

abstract: When studying families in the moduli space of dynamical systems, choosing an appropriate representative function for a conjugacy class can be a delicate task. The most delicate questions surround rationality of the conjugacy class compared to rationality of the defining polynomials of the representation. We give a normal form for degree three polynomials which has the property that the set of fixed points is equal to the set of fixed point multipliers. This normal form is given in terms of moduli space invariants and, hence, has nice rationality properties. We further classify all degree three rational maps which can be conjugated to have a similar relationship between the fixed points and the fixed point multipliers.


title: On the stability of the Lions-Peetre method of real interpolation with functional parameter

abstract: Let $\vec{X}=(X_0, X_1)$ be a compatible couple of Banach spaces, $ 1\le p \le \infty$ and let $ \varphi$ be positive quasi-concave function. Denote by $\overline{X}_{\varphi,p}=(X_0,X_1)_{\varphi,p}$ the real interpolation spaces defined by S. Janson (1981). We give necessary and sufficient conditions on $ \varphi_{0}$, $\varphi_{1}$ and $\varphi$ for the validity of \begin{equation*} \left(\overline{X}_{\varphi_{0},1},\overline{X}_{\varphi_{1},1} \right) _{\varphi,p}= \left(\overline{X}_{\varphi_{0},\infty},\overline{X}_{\varphi_{1},\infty}\right)_{\varphi,p} \end{equation*} for all $ 1\le p\le \infty$, and all Banach couples $\overline{X}. $


title: Comparing School Choice and College Admission Mechanisms By Their Immunity to Strategic Admissions

abstract: Recently dozens of school districts and college admissions systems around the world have reformed their admission rules. As a main motivation for these reforms the policymakers cited strategic flaws of the rules: students had strong incentives to game the system, which caused dramatic consequences for non-strategic students. However, almost none of the new rules were strategy-proof. We explain this puzzle. We show that after the reforms the rules became more immune to strategic admissions: each student received a smaller set of schools that he can get in using a strategy, weakening incentives to manipulate. Simultaneously, the admission to each school became strategy-proof to a larger set of students, making the schools more available for non-strategic students. We also show that the existing explanation of the puzzle due to Pathak and S\"onmez (2013) is incomplete.


title: Spectra of the Rarita-Schwinger operator on some symmetric spaces

abstract: We give a method to calculate spectra of the square of the Rarita-Schwinger operator on compact symmetric spaces. According to Weitzenb\"{o}ck formulas, the operator can be written by the Laplace operator, which is the Casimir operator on compact symmetric spaces. Then we can obtain the spectra by using the Freudenthal's formula and branching rules. As examples, we calculate the spectra on the sphere, the complex projective space, and the quaternionic projective space.


title: Optimal Crossover Designs for Generalized Linear Models

abstract: We identify locally $D$-optimal crossover designs for generalized linear models. We use generalized estimating equations to estimate the model parameters along with their variances. To capture the dependency among the observations coming from the same subject, we propose six different correlation structures. We identify the optimal allocations of units for different sequences of treatments. For two-treatment crossover designs, we show via simulations that the optimal allocations are reasonably robust to different choices of the correlation structures. We discuss a real example of multiple treatment crossover experiments using Latin square designs. Using a simulation study, we show that a two-stage design with our locally $D$-optimal design at the second stage is more efficient than the uniform design, especially when the responses from the same subject are correlated.


title: Lower bounds for the maximum number of runners that cause loneliness, and its application to Isolation

abstract: We consider (n+1) runners with given constant unique integer speeds running along the circumference of a circle whose circumferential length is one, and all runners starting from the same point. We define and give lower bounds to a first problem PMAX of finding, for every runner r, the maximum number of runners that can be simultaneously separated from runner r by a distance of atleast d. For d=1/(2^(floor(lg(n)))), a lower bound for PMAX is ( n - ((n-1)/floor(lg(n))) ), which makes the fraction of simultaneously separated runners tend to 1 as n tends to infinity. Next, we define and give upper bounds to a second problem ISOLATE of finding, for every runner r, the minimum number of steps needed to isolate r, assuming that the runners that can be simultaneously separated from r by atleast d, are removed at each step. For d=1/(2^(floor(lg(n)))), an upper bound for ISOLATE is ( lg(n - 1)/lg(floor(lg(n))) ).


title: Precision predictions for $B \rightarrow \rho \tau \nu_\tau$ and $B \rightarrow \omega \tau \nu_\tau$ in the SM and beyond

abstract: We present new precision predictions for semitauonic decays involving $\rho$ and $\omega$ final state mesons. These decay channels offer an interesting orthogonal probe to study the existing B anomalies in semitauonic transitions and are accessible with the Belle II experiment. The predictions are based on combining existing light-cone sum-rule calculations for the form factors with measured experimental spectra from the BaBar and Belle collaborations. This allows us to reliably extrapolate the light-lepton form factor predictions to large values of the four-momentum transfer squared, $q^2$, and in turn to derive precise predictions for $R_{\rho}$ and $R_{\omega}$, the ratio of the total decay rates of $B \rightarrow \rho \tau \nu_\tau$ and $B \rightarrow \omega \tau \nu_\tau$ for $\tau$ final states with respect to light leptons in the SM. In addition, we investigate the impact of all four-fermi operators on the semitauonic $q^2$ spectra and these ratios.


title: FPCR-Net: Feature Pyramidal Correlation and Residual Reconstruction for Optical Flow Estimation

abstract: Optical flow estimation is an important yet challenging problem in the field of video analytics. The features of different semantics levels/layers of a convolutional neural network can provide information of different granularity. To exploit such flexible and comprehensive information, we propose a semi-supervised Feature Pyramidal Correlation and Residual Reconstruction Network (FPCR-Net) for optical flow estimation from frame pairs. It consists of two main modules: pyramid correlation mapping and residual reconstruction. The pyramid correlation mapping module takes advantage of the multi-scale correlations of global/local patches by aggregating features of different scales to form a multi-level cost volume. The residual reconstruction module aims to reconstruct the sub-band high-frequency residuals of finer optical flow in each stage. Based on the pyramid correlation mapping, we further propose a correlation-warping-normalization (CWN) module to efficiently exploit the correlation dependency. Experiment results show that the proposed scheme achieves the state-of-the-art performance, with improvement by 0.80, 1.15 and 0.10 in terms of average end-point error (AEE) against competing baseline methods - FlowNet2, LiteFlowNet and PWC-Net on the Final pass of Sintel dataset, respectively.


title: Persistence and Stability of the $\mathbb{A}_n$ Quiver

abstract: We introduce two new distances for zigzag persistence modules. The first uses Auslander-Reiten quiver theory, and the second is an extension of the classical interleaving distance. Both are defined over completely general orientations of the $\mathbb{A}_n$ quiver. We compare the first distance to the block distance introduced by M. Botnan and M. Lesnick and obtain the full set of sharp Lipschitz bounds between the two (as bottleneck distances) over pure zigzag orientations. The final portion of the paper presents sharp Lipschitz bounds necessary for the extended interleaving distance to dominate the distance that is created from the Auslander-Reiten quiver. These bounds are obtained for general orientations of the $\mathbb{A}_n$ quiver.


title: Floquet eigenspectra of a nonlinear two-mode system under periodic driving: the emergence of "ring" structures

abstract: We study Floquet eigenspectra of a nonlinear two-mode system under a periodic driving of the off-diagonal coupling. By solving the Gross-Pitaevskii equation numerically, we obtain triangular and loop structures near the crossings of different Floquet branches. At lower driving frequencies, we find "ring" and "double-ring" structures which are distinct from the well-known loop structure. The mechanism of the emergence of these structures is discussed and the parameter windows of their existence are obtained analytically. In addition, we study the evolution of the system under the driving with an adiabatic sweep and find there are some dynamically unstable states in the Floquet eigenspectra which break the quantum adiabaticity.


title: A new method for measurement and quantification of tracer diffusion in nanoconfined liquids

abstract: We report development of a novel instrument to measure tracer diffusion in water under nano-scale confinement. A direct optical access to the confinement region, where water is confined between a tapered fiber and a flat substrate, is made possible by coating the probe with metal and opening a small aperture( 0.1 $\mu$m to 1 $\mu$m) at its end. A well-controlled cut using an ion beam ensures desired lateral confinement area as well as adequate illumination of the confinement gap. The probe is mounted on a tuning-fork based force sensor to control the separation between the probe and the substrate with nm precision. Fluctuations in fluorescence intensity due to diffusion of a dye molecule in water confined between probe and the sample are recorded using a confocal arrangement with a single photon precision. A Monte-Carlo method is developed to determine the diffusion coefficient from the measured autocorrelation of intensity fluctuations which accommodates the specific geometry of confinement and the illumination profile. The instrument allows measurement of diffusion laws under confinement. We found that the diffusion of a tracer molecule is slowed down by more than ten times for the probe-substrate separations of 5 nm and below.


title: Spatiotemporal Camera-LiDAR Calibration: A Targetless and Structureless Approach

abstract: The demand for multimodal sensing systems for robotics is growing due to the increase in robustness, reliability and accuracy offered by these systems. These systems also need to be spatially and temporally co-registered to be effective. In this paper, we propose a targetless and structureless spatiotemporal camera-LiDAR calibration method. Our method combines a closed-form solution with a modified structureless bundle adjustment where the coarse-to-fine approach does not {require} an initial guess on the spatiotemporal parameters. Also, as 3D features (structure) are calculated from triangulation only, there is no need to have a calibration target or to match 2D features with the 3D point cloud which provides flexibility in the calibration process and sensor configuration. We demonstrate the accuracy and robustness of the proposed method through both simulation and real data experiments using multiple sensor payload configurations mounted to hand-held, aerial and legged robot systems. Also, qualitative results are given in the form of a colorized point cloud visualization.


title: A proximal MM method for the zero-norm regularized PLQ composite optimization problem

abstract: This paper is concerned with a class of zero-norm regularized piecewise linear-quadratic (PLQ) composite minimization problems, which covers the zero-norm regularized $\ell_1$-loss minimization problem as a special case. For this class of nonconvex nonsmooth problems, we show that its equivalent MPEC reformulation is partially calm on the set of global optima and make use of this property to derive a family of equivalent DC surrogates. Then, we propose a proximal majorization-minimization (MM) method, a convex relaxation approach not in the DC algorithm framework, for solving one of the DC surrogates which is a semiconvex PLQ minimization problem involving three nonsmooth terms. For this method, we establish its global convergence and linear rate of convergence, and under suitable conditions show that the limit of the generated sequence is not only a local optimum but also a good critical point in a statistical sense. Numerical experiments are conducted with synthetic and real data for the proximal MM method with the subproblems solved by a dual semismooth Newton method to confirm our theoretical findings, and numerical comparisons with a convergent indefinite-proximal ADMM for the partially smoothed DC surrogate verify its superiority in the quality of solutions and computing time.


title: Experimental quantum polarimetry using heralded single photons

abstract: We perform experimental quantum polarimetry using a heralded single photon to analyze the optical activity of linearly polarized light traversing a chiral medium. Three kinds of estimators are considered to estimate the concentrations of sucrose solutions from measuring the rotation angle of the linear polarization of the output photons. Through repetition of independent and identical measurements performed for each individual scheme and different concentration sucrose solutions, we compare the estimation uncertainty among the three schemes. The results are also compared to classical benchmarks for which a coherent state of light is taken into account. The quantum enhancement in the estimation uncertainty is evaluated and the impact of experimental and technical imperfections is discussed. In this work, we lay out a route for future applications relying on quantum polarimetry.


title: DNNs as Layers of Cooperating Classifiers

abstract: A robust theoretical framework that can describe and predict the generalization ability of deep neural networks (DNNs) in general circumstances remains elusive. Classical attempts have produced complexity metrics that rely heavily on global measures of compactness and capacity with little investigation into the effects of sub-component collaboration. We demonstrate intriguing regularities in the activation patterns of the hidden nodes within fully-connected feedforward networks. By tracing the origin of these patterns, we show how such networks can be viewed as the combination of two information processing systems: one continuous and one discrete. We describe how these two systems arise naturally from the gradient-based optimization process, and demonstrate the classification ability of the two systems, individually and in collaboration. This perspective on DNN classification offers a novel way to think about generalization, in which different subsets of the training data are used to train distinct classifiers; those classifiers are then combined to perform the classification task, and their consistency is crucial for accurate classification.


title: Hyper-positive definite functions II: A complete study of branching-Toeplitz operators

abstract: We introduce and give a more or less complete study of a family of branching-Toeplitz operators on the Hilbert space $\ell^2(T_q)$ indexed by a rooted homogeneous tree $T_q$ of degree $q\ge 2$. The finite dimensional analogues of such operators form a very natural family of structured sparse matrices called branching-Toeplitz matrices and will also be investigated. The branching-Toeplitz operators/matrices in this paper should be viewed as natural generalizations of the standard Toeplitz operators/matrices. We will apply our results to construct a family of determinantal point processes on homogeneous trees which are branching-type strong stationary stochastic processes.


title: MAGMO: Polarimetry of 1720-MHz OH Masers towards Southern Star Forming Regions

abstract: From targeted observations of ground-state OH masers towards 702 Multibeam (MMB) survey 6.7-GHz methanol masers, between Galactic longitudes 186$^{\circ}$ through the Galactic centre to 20$^{\circ}$, made as part of the `MAGMO' project, we present the physical and polarisation properties of the 1720-MHz OH maser transition, including the identification of Zeeman pairs. We present 10 new and 23 previously catalogued 1720-MHz OH maser sources detected towards star formation regions. In addition, we also detected 16 1720-MHz OH masers associated with supernova remnants and two sites of diffuse OH emission. Towards the 33 star formation masers, we identify 44 Zeeman pairs, implying magnetic field strengths ranging from $-$11.4 to $+$13.2 mG, and a median magnetic field strength of $|B_{LOS}|$ $\sim$ 6 mG. With limited statistics, we present the in-situ magnetic field orientation of the masers and the Galactic magnetic field distribution revealed by the 1720-MHz transition. We also examine the association statistics of 1720-MHz OH SFR masers with other ground-state OH masers, excited-state OH masers, class I and class II methanol masers and water masers, and compare maser positions with mid-infrared images of the parent star forming regions. Of the 33 1720-MHz star formation masers, ten are offset from their central exciting sources, and appear to be associated with outflow activity.


title: Model Predictive Control of Discrete-Continuous Energy Systems via Generalized Disjunctive Programming

abstract: Generalized Disjunctive Programming (GDP) provides an alternative framework to model optimization problems with both discrete and continuous variables. The key idea behind GDP involves the use of logical disjunctions to represent discrete decisions in the continuous space, and logical propositions to denote algebraic constraints in the discrete space. Compared to traditional mixed-integer programming (MIP), the inherent logic structure in GDP yields tighter relaxations that are exploited by global branch and bound algorithms to improve solution quality. In this paper, we present a general GDP model for optimal control of hybrid systems that exhibit both discrete and continuous dynamics. Specifically, we use GDP to formulate a model predictive control (MPC) model for piecewise-affine systems with implicit switching logic. As an example, the GDP-based MPC approach is used as a supervisory control to improve energy efficiency in residential buildings with binary on/off, relay-based thermostats. A simulation study is used to demonstrate the validity of the proposed approach, and the improved solution quality compared to existing MIP-based control approaches.


title: SRPerf: a Performance Evaluation Framework for IPv6 Segment Routing

abstract: Segment Routing is a form of loose source routing. It provides the ability to include a list of instructions (called segments), in the packet headers. The Segment Routing architecture has been first implemented with the MPLS dataplane and then, quite recently, with the IPv6 dataplane (SRv6). IPv6 Segment Routing (SRv6) is a promising solution to support advanced services such as Traffic Engineering, Service Function Chaining, Virtual Private Networks, and Load Balancing. The SRv6 data-plane is supported in many different software forwarding engines including the Linux kernel and VPP software router, as well as in hardware devices. In this paper, we present SRPerf, a performance evaluation framework for software and hardware implementations of SRv6. SRPerf is able to perform different benchmarking tests such as throughput and latency. For throughput tests, we use the Partial Drop Rate (PDR) to characterize a system under test. The architecture of SRPerf can be easily extended to support new benchmarking methodologies as well as different SRv6 implementations. We have used SRPerf to evaluate the performance of the SRv6 implementation in the Linux kernel and in VPP. SRPerf is a valuable tool in the context of software forwarding engines where new features can be added at fast pace, as it helps experimenters to validate their work. In particular, we describe how we have leveraged SRPerf to validate the implementation of some SRv6 behaviors that were missing or wrongly implemented in the Linux kernel mainline.


title: Theory of Higgs Modes in $d$-Wave Superconductors

abstract: By applying a microscopic gauge-invariant kinetic theory in $d$-wave superconductors, we analytically derive the energy spectra of the breathing Higgs mode and, in particular, rotating Higgs mode that is unique for the $d$-wave order parameter. Analytical investigation on their dynamic properties is also revealed. We show that the breathing Higgs mode is optically visible in the second-order regime. Whereas the rotating Higgs mode is optically inactive, we show that this mode can be detected in the pseudogap phase by magnetic resonance experiment. It is interesting to find that with a longitudinal temperature gradient, the charge-neutral rotating Higgs mode generates a thermal Hall current by magnetic field in the pseudogap phase, providing a unique scheme for its detection.


title: Dynamic Phase Transition in the Contact Process with Spatial Disorder: Griffths Phase and Complex Persistence Exponents

abstract: We present a model which displays Griffiths phase i.e. algebraic decay of density with continuously varying exponent in the absorbing phase. In active phase, the memory of initial conditions is lost with continuously varying complex exponent in this model. This is 1-D model where fraction r of sites obey rules leading to directed percolation (DP) class and the rest evolve according to rules leading to compact directed percolation (CDP) class. For infection probability $p < p_c$, the fraction of active sites $\rho(t) = 0$ asymptotically. For $p > p_c$, $\rho(infty) > 0$. At $p = p_c$, $\rho(t)$, the survival probability from single seed and the average number of active sites starting from single seed decay logarithmically. The local persistence $P_l(\infty) > 0$ for $p < p_c$ and $P_l(\infty) = 0$ for $p > p_c$. For $p > p_s$, local persistence $P_l(t)$ decays as a power law with continuously varying exponents. The persistence exponent is clearly complex as $p\rightarrow 1$. The complex exponent implies logarithmic periodic oscillations in persistence. The wavelength and the amplitude of the logarithmic periodic oscillations increases with p. We note that underlying lattice or disorder does not have self-similar structure.


title: Frequency- and Time-Limited Balanced Truncation for Large-Scale Second-Order Systems

abstract: Considering the use of dynamical systems in practical applications, often only limited regions in the time or frequency domain are of interest. Therefor, it usually pays off to compute local approximations of the used dynamical systems in the frequency and time domain. In this paper, we consider a structure-preserving extension of the frequency- and time-limited balanced truncation methods to second-order dynamical systems. We give a full overview about the first-order limited balanced truncation methods and extend those methods to second-order systems by using the different second-order balanced truncation formulas from the literature. Also, we present numerical methods for solving the arising large-scale sparse matrix equations and give numerical modifications to deal with the problematic case of second-order systems. The results are then illustrated on three numerical examples.


title: Structural properties of additive binary hard-sphere mixtures

abstract: An approach to obtain the structural properties of additive binary hard-sphere mixtures is presented. Such an approach, which is a nontrivial generalization of the one recently used for monocomponent hard-sphere fluids [S. Pieprzyk, A. C. Bra\'nka, and D. M. Heyes, Phys. Rev. E 95, 062104 (2017)], combines accurate molecular-dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation. A comparison of the direct correlation functions obtained with the present scheme with those derived from theoretical results stemming from the Percus-Yevick (PY) closure and the so-called rational-function approximation (RFA) is performed. The density dependence of the leading poles of the Fourier transforms of the total correlation functions and the decay of the pair correlation functions of the mixtures are also addressed and compared to the predictions of the two theoretical approximations. A very good overall agreement between the results of the present scheme and those of the RFA is found, thus suggesting that the latter (which is an improvement over the PY approximation) can safely be used to predict reasonably well the long-range behavior, including the structural crossover, of the correlation functions of additive binary hard-sphere mixtures.


title: Exponential contraction in Wasserstein distance on static and evolving manifolds

abstract: In this article, exponential contraction in Wasserstein distance for heat semigroups of diffusion processes on Riemannian manifolds is established under curvature conditions where Ricci curvature is not necessarily required to be non-negative. Compared to the results of Wang (2016), we focus on explicit estimates for the exponential contraction rate. Moreover, we show that our results extend to manifolds evolving under a geometric flow. As application, for the time-inhomogeneous semigroups, we obtain a gradient estimate with an exponential contraction rate under weak curvature conditions, as well as uniqueness of the corresponding evolution system of measures.


title: On Random Matrices Arising in Deep Neural Networks. Gaussian Case

abstract: The paper deals with distribution of singular values of product of random matrices arising in the analysis of deep neural networks. The matrices resemble the product analogs of the sample covariance matrices, however, an important difference is that the population covariance matrices, which are assumed to be non-random in the standard setting of statistics and random matrix theory, are now random, moreover, are certain functions of random data matrices. The problem has been considered in recent work [21] by using the techniques of free probability theory. Since, however, free probability theory deals with population matrices which are independent of the data matrices, its applicability in this case requires an additional justification. We present this justification by using a version of the standard techniques of random matrix theory under the assumption that the entries of data matrices are independent Gaussian random variables. In the subsequent paper [18] we extend our results to the case where the entries of data matrices are just independent identically distributed random variables with several finite moments. This, in particular, extends the property of the so-called macroscopic universality on the considered random matrices.


title: Ab-initio Theory of Fourier-transformed Quasiparticle Interference Maps and Application to the Topological Insulator Bi$_2$Te$_3$

abstract: The quasiparticle interference (QPI) technique is a powerful tool that allows to uncover the structure and properties of electronic structure of a material combined with scattering properties of defects at surfaces. Recently this technique has been pivotal in proving the unique properties of the surface state of topological insulators which manifests itself in the absence of backscattering. In this work we derive a Green function based formalism for the ab initio computation of Fourier-transformed QPI images. We show the efficiency of our new implementation at the examples of QPI that forms around magnetic and non-magnetic defects at the Bi$_2$Te$_3$ surface. This method allows a deepened understanding of the scattering properties of topologically protected electrons off defects and can be a useful tool in the study of quantum materials in the future.


title: Visual Simplified Characters' Emotion Emulator Implementing OCC Model

abstract: In this paper, we present a visual emulator of the emotions seen in characters in stories. This system is based on a simplified view of the cognitive structure of emotions proposed by Ortony, Clore and Collins (OCC Model). The goal of this paper is to provide a visual platform that allows us to observe changes in the characters' different emotions, and the intricate interrelationships between: 1) each character's emotions, 2) their affective relationships and actions, 3) The events that take place in the development of a plot, and 4) the objects of desire that make up the emotional map of any story. This tool was tested on stories with a contrasting variety of emotional and affective environments: Othello, Twilight, and Harry Potter, behaving sensibly and in keeping with the atmosphere in which the characters were immersed.


title: Symmetry-enforced three-dimensional Dirac phononic crystals

abstract: Dirac semimetals, the materials featured with discrete linearly crossing points (called Dirac points) between four bands, are critical states of topologically distinct phases. Such gapless topological states have been accomplished by a band-inversion mechanism, in which the Dirac points can be annihilated pairwise by perturbations without changing the symmetry of the system. Here, we report an experimental observation of Dirac points that are enforced completely by the crystal symmetry, using a nonsymmorphic three-dimensional phononic crystal. Intriguingly, our Dirac phononic crystal hosts four spiral topological surface states, in which the surface states of opposite helicities intersect gaplessly along certain momentum lines, as confirmed by our further surface measurements. The novel Dirac system may release new opportunities for studying the elusive (pseudo)relativistic physics, and also offer a unique prototype platform for acoustic applications.


title: Dynamical blockade in a single mode bosonic system

abstract: We introduce a dynamical blockade phenomenon occurring in a nonlinear bosonic mode induced by a combination of continuous and pulsed excitations. We find that the underlying mechanism for the blockade is general, enhancing antibunching in the strongly nonlinear regime and inducing it in the weakly nonlinear regime, without fine-tuning the system parameters. Moreover, this mechanism shows advantages over existing blockade mechanisms and is suitable for implementation in a wide variety of systems due to its simplicity and universality.

categories: quant-ph cond-mat.mtrl-sci physics.atom-ph physics.optics


title: Status of direct and indirect dark matter searches

abstract: I review the current status of dark matter searches using direct, indirect and accelerator techniques. A detailed review of individual experiments is beyond the scope of these proceedings. I focus instead on the challenges (sometimes limitations) faced by each of the approaches, which is what make them complementary, and the reason we must ensure that they are developed concurrently.


title: Communication-Efficient Distributed Estimator for Generalized Linear Models with a Diverging Number of Covariates

abstract: Distributed statistical inference has recently attracted immense attention. The asymptotic efficiency of the maximum likelihood estimator (MLE), the one-step MLE, and the aggregated estimating equation estimator are established for generalized linear models under the "large $n$, diverging $p_n$" framework, where the dimension of the covariates $p_n$ grows to infinity at a polynomial rate $o(n^\alpha)$ for some $0<\alpha<1$. Then a novel method is proposed to obtain an asymptotically efficient estimator for large-scale distributed data by two rounds of communication. In this novel method, the assumption on the number of servers is more relaxed and thus practical for real-world applications. Simulations and a case study demonstrate the satisfactory finite-sample performance of the proposed estimators.


title: On the conservation of energy in two-dimensional incompressible flows

abstract: We prove the conservation of energy for weak and statistical solutions of the two-dimensional Euler equations, generated as strong (in an appropriate topology) limits of the underlying Navier-Stokes equations and a Monte Carlo-Spectral Viscosity numerical approximation, respectively. We characterize this conservation of energy in terms of a uniform decay of the so-called structure function, allowing us to extend existing results on energy conservation. Moreover, we present numerical experiments with a wide variety of initial data to validate our theory and to observe energy conservation in a large class of two-dimensional incompressible flows.


title: Dissociated dislocation-mediated carbon transport and diffusion in austenitic iron

abstract: Dislocation-solute interaction plays fundamental roles in mechanical properties of alloys. Here, we disclose the essential features of dislocation-carbon interaction in austenitic Fe at the atomistic scale. We show that passage of a Shockley partial dislocation in face-centered cubic iron is able to move carbon atoms on the slip plane forward by one Burgers vector, revealing a novel dissociated dislocation-mediated transport mechanism. This mechanism is induced by shear, which is distinct from the normal thermally activated diffusion process. Furthermore, we show that there exists a fast diffusion channel with significantly reduced diffusion energy barrier in the partial dislocation core, which is highly localized and directional. These inherent geometrical features are crucial for understanding the dependence of the diffusivity of dislocation pipe diffusion on the character of dislocations; most importantly, they can result in unbalanced pinning effect on the leading and trailing partials in a mixed dislocation, consequently facilitating stacking fault formation and deformation twinning. This explains the controversial effects of carbon on deformation twinning observed in various alloys. Our findings pave the road to tune mechanical properties of materials by manipulating dislocation-interstitial interaction.


title: Daugavet- and Delta-points in absolute sums of Banach spaces

abstract: A Daugavet-point (resp.~$\Delta$-point) of a Banach space is a norm one element $x$ for which every point in the unit ball (resp.~element $x$ itself) is in the closed convex hull of unit ball elements that are almost at distance 2 from $x$. A Banach space has the well-known Daugavet property (resp.~diametral local diameter 2 property) if and only if every norm one element is a Daugavet-point (resp.~$\Delta$-point). This paper complements the article "Delta- and Daugavet-points in Banach spaces" by T. A. Abrahamsen, R. Haller, V. Lima, and K. Pirk, where the study of the existence of Daugavet- and $\Delta$-points in absolute sums of Banach spaces was started.


title: Inducing Magnetic Phase Transitions in Monolayer CrI$_3$ via Lattice Deformations

abstract: Atomically thin films of layered chromium triiodide (CrI$_3$) have recently been regarded as suitable candidates to a wide spectrum of technologically relevant applications, mainly owing to the opportunity they offer to achieve a reversible transition between coexisting in-plane ferro- and out-of-plane antiferro-magnetic orders. However, no routes for inducing such a transition have been designed down to the single-layer limit. Here, we address the magnetic response of monolayer CrI$_3$ to in-plane lattice deformations through a combination of isotropic Heisenberg spin Hamiltonians and first-principles calculations. Depending on the magnitude and orientation of the lattice strain exerted, we unveil a series of direction-dependent parallel-to-antiparallel spins crossovers, which yield the emergence of ferromagnetic, N\'eel antiferromagnetic, zigzag and stripy antiferromagnetic ground states. Additionally, we identify a critical point in the magnetic phase diagram whereby the exchange couplings vanish and the magnetism is quenched. Our work establishes guidelines for extensively tailoring the spin interactions in monolayer CrI$_3$ via strain engineering, and further expands the magnetically ordered phases which can be hosted in a two-dimensional crystal.


title: Two-pion scattering amplitude from Bethe-Salpeter wave function at the interaction boundary

abstract: We observe that the ratio of the on-shell scattering amplitude to the Bethe-Salpeter (BS) wave function outside the interaction range is almost independent of time in our quenched calculation of the $I=2$ two-pion scattering with almost zero momentum. In order to discuss the time independence, we present a relation between the two-pion scattering amplitude and the surface term of the BS wave function at the boundary. Using the relation under some assumptions, we show that the ratio is independent of time if the two-pion four-point function in early time is dominated by scattering states with almost zero momentum in addition to the ground state of the two-pion scattering.


title: Fallacies on pairing symmetry and intrinsic electronic Raman spectrum in high-Tc cuprate superconductors

abstract: Certain significant fallacies are involved in discussions of the high-Tc mechanism unsolved for over 30 years in cuprate superconductors. These fallacies are explored with the aim of unravelling this mechanism. Moreover, using polarised electronic Raman scattering in inhomogeneous underdoped cuprate superconductors, the intrinsic nonlinear Raman spectrum is obtained by subtracting the pseudogap characteristic of a nonlinear from the linear Raman spectrum measured in the B2g mode of the node area below the critical temperature. The intrinsic nonlinear behaviour implies the existence of the nodal superconducting gap denying dx2-y2-wave pairing symmetry. An origin of the nodal superconducting gap is discussed.


title: Resonance Line Polarization in Spherically Symmetric Moving Media: a Parametric Study

abstract: In the present paper we consider the problem of resonance line polarization formed in the spherically symmetric expanding atmospheres. For the solution of the concerned polarized transfer equation we use the comoving frame formulation, and apply the Accelerated Lambda Iteration (ALI) method. We restrict ourselves to the non-relativistic regime of velocities wherein mainly Doppler shift effects are significant. For our studies, we consider the scattering on a two-level atom, including the effects of partial frequency redistribution (PFR). We present the dependence of linearly polarized profiles on different atmospheric and atomic parameters.


title: FedVision: An Online Visual Object Detection Platform Powered by Federated Learning

abstract: Visual object detection is a computer vision-based artificial intelligence (AI) technique which has many practical applications (e.g., fire hazard monitoring). However, due to privacy concerns and the high cost of transmitting video data, it is highly challenging to build object detection models on centrally stored large training datasets following the current approach. Federated learning (FL) is a promising approach to resolve this challenge. Nevertheless, there currently lacks an easy to use tool to enable computer vision application developers who are not experts in federated learning to conveniently leverage this technology and apply it in their systems. In this paper, we report FedVision - a machine learning engineering platform to support the development of federated learning powered computer vision applications. The platform has been deployed through a collaboration between WeBank and Extreme Vision to help customers develop computer vision-based safety monitoring solutions in smart city applications. Over four months of usage, it has achieved significant efficiency improvement and cost reduction while removing the need to transmit sensitive data for three major corporate customers. To the best of our knowledge, this is the first real application of FL in computer vision-based tasks.


title: Low-Cost Anti-Copying 2D Barcode by Exploiting Channel Noise Characteristics

abstract: In this paper, for overcoming the drawbacks of the prior approaches, such as low generality, high cost, and high overhead, we propose a Low-Cost Anti-Copying (LCAC) 2D barcode by exploiting the difference between the noise characteristics of legal and illegal channels. An embedding strategy is proposed, and for a variant of it, we also make the corresponding analysis. For accurately evaluating the performance of our approach, a theoretical model of the noise in an illegal channel is established by using a generalized Gaussian distribution. By comparing with the experimental results based on various printers, scanners, and a mobile phone, it can be found that the sample histogram and curve fitting of the theoretical model match well, so it can be concluded that the theoretical model works well. For evaluating the security of the proposed LCAC code, besides the direct-copying (DC) attack, the improved version, which is the synthesized-copying (SC) attack, is also considered in this paper. Based on the theoretical model, we build a prediction function to optimize the parameters of our approach. The parameters optimization incorporates the covertness requirement, the robustness requirement and a tradeoff between the production cost and the cost of illegally-copying attacks together. The experimental results show that the proposed LCAC code with two printers and two scanners can detect the DC attack effectively and resist the SC attack up to the access of 14 legal copies.


title: A Note on Computable Embeddings for Ordinals and Their Reverses

abstract: We continue the study of computable embeddings for pairs of structures, i.e. for classes containing precisely two non-isomorphic structures. Surprisingly, even for some pairs of simple linear orders, computable embeddings induce a non-trivial degree structure. Our main result shows that although $\{\omega \cdot 2, \omega^\star \cdot 2\}$ is computably embeddable in $\{\omega^2, {(\omega^2)}^\star\}$, the class $\{\omega \cdot k,\omega^\star \cdot k\}$ is \emph{not} computably embeddable in $\{\omega^2, {(\omega^2)}^\star\}$ for any natural number $k \geq 3$.


title: Discovering topological surface states of Dirac points

abstract: Dirac materials, unlike the Weyl materials, have not been found in experiments to support intrinsic topological surface states, as the surface arcs in existing systems are unstable against symmetry-preserving perturbations. Utilizing the proposed glide and time-reversal symmetries, we theoretically design and experimentally verify an acoustic crystal of two frequency-isolated three-dimensional Dirac points with $Z_2$ monopole charges and four gapless helicoid surface states.


title: Multi-step Joint-Modality Attention Network for Scene-Aware Dialogue System

abstract: Understanding dynamic scenes and dialogue contexts in order to converse with users has been challenging for multimodal dialogue systems. The 8-th Dialog System Technology Challenge (DSTC8) proposed an Audio Visual Scene-Aware Dialog (AVSD) task, which contains multiple modalities including audio, vision, and language, to evaluate how dialogue systems understand different modalities and response to users. In this paper, we proposed a multi-step joint-modality attention network (JMAN) based on recurrent neural network (RNN) to reason on videos. Our model performs a multi-step attention mechanism and jointly considers both visual and textual representations in each reasoning process to better integrate information from the two different modalities. Compared to the baseline released by AVSD organizers, our model achieves a relative 12.1% and 22.4% improvement over the baseline on ROUGE-L score and CIDEr score.


title: The area minimizing problem in conformal cones

abstract: In this paper we study the area minimizing problem in some kinds of conformal cones. This concept is a generalization of the cones in Eulcidean spaces and the cylinders in product manifolds. We define a non-closed-minimal (NCM) condition for bounded domains. Under this assumption and other necessary conditions we establish the existence of bounded minimal graphs in mean convex conformal cones. Moreover those minimal graphs are the solutions to corresponding area minizing problems. We can solve the area minimizing problem in non-mean convex translating conformal cones if these cones are contained in a larger mean convex conformal cones with the NCM assumption. We give examples to illustrate that this assumption can not be removed for our main results.


title: Causal models for dynamical systems

abstract: A probabilistic model describes a system in its observational state. In many situations, however, we are interested in the system's response under interventions. The class of structural causal models provides a language that allows us to model the behaviour under interventions. It can been taken as a starting point to answer a plethora of causal questions, including the identification of causal effects or causal structure learning. In this chapter, we provide a natural and straight-forward extension of this concept to dynamical systems, focusing on continuous time models. In particular, we introduce two types of causal kinetic models that differ in how the randomness enters into the model: it may either be considered as observational noise or as systematic driving noise. In both cases, we define interventions and therefore provide a possible starting point for causal inference. In this sense, the book chapter provides more questions than answers. The focus of the proposed causal kinetic models lies on the dynamics themselves rather than corresponding stationary distributions, for example. We believe that this is beneficial when the aim is to model the full time evolution of the system and data are measured at different time points. Under this focus, it is natural to consider interventions in the differential equations themselves.


title: Registration made easy -- standalone orthopedic navigation with HoloLens

abstract: In surgical navigation, finding correspondence between preoperative plan and intraoperative anatomy, the so-called registration task, is imperative. One promising approach is to intraoperatively digitize anatomy and register it with the preoperative plan. State-of-the-art commercial navigation systems implement such approaches for pedicle screw placement in spinal fusion surgery. Although these systems improve surgical accuracy, they are not gold standard in clinical practice. Besides economical reasons, this may be due to their difficult integration into clinical workflows and unintuitive navigation feedback. Augmented Reality has the potential to overcome these limitations. Consequently, we propose a surgical navigation approach comprising intraoperative surface digitization for registration and intuitive holographic navigation for pedicle screw placement that runs entirely on the Microsoft HoloLens. Preliminary results from phantom experiments suggest that the method may meet clinical accuracy requirements.


title: Unique continuation property and Poincar\'e inequality for higher order fractional Laplacians with applications in inverse problems

abstract: We prove a unique continuation property for the fractional Laplacian $(-\Delta)^s$ when $s \in (-n/2,\infty)\setminus \mathbb{Z}$. In addition, we study Poincar\'e-type inequalities for the operator $(-\Delta)^s$ when $s\geq 0$. We apply the results to show that one can uniquely recover, up to a gauge, electric and magnetic potentials from the Dirichlet-to-Neumann map associated to the higher order fractional magnetic Schr\"odinger equation. We also study the higher order fractional Schr\"odinger equation with singular electric potential. In both cases, we obtain a Runge approximation property for the equation. Furthermore, we prove a uniqueness result for a partial data problem of the $d$-plane Radon transform in low regularity. Our work extends some recent results in inverse problems for more general operators.


title: Incomplete selected inversion for linear-scaling electronic structure calculations

abstract: Pole Expansion and Selected Inversion (PEXSI) is an efficient scheme for evaluating selected entries of functions of large sparse matrices as required e.g. in electronic structure algorithms. We show that the triangular factorizations computed by the PEXSI scheme exhibit a localization property similar to that of matrix functions, and we present a modified PEXSI algorithm which exploits this observation to achieve linear scaling. To the best of our knowledge, the resulting incomplete PEXSI (iPEXSI) algorithm is the first linear-scaling algorithm which scales provably better than cubically even in the absence of localization, and we hope that this will help to further lower the critical system size where linear-scaling algorithms begin to outperform the diagonalization algorithm.


title: Neural coincidence detection strategies during perception of multi-pitch musical tones

abstract: Multi-pitch perception is investigated in a listening test using 30 recordings of musical sounds with two tones played simultaneously, except for two gong sounds with inharmonic overtone spectrum, judging roughness and separateness as the ability to tell the two tones in each recording apart. 13 sounds were from a Western guitar playing all 13 intervals in one octave, the other sounds were mainly from non-Western instruments comparing familiar with unfamiliar instrument sounds for Western listeners. Additionally the sounds were processed in a cochlear model transferring the mechanical basilar membrane motion into neural spikes followed by a post-processing simulating different degrees of coincidence detection. Separateness perception showed a clear distinction between familiar and unfamiliar sounds, while roughness perception did not. By correlating perception with simulation different perception strategies were found. Familiar sounds correlated strongly positive with high degrees of coincidence detection, where only 3-5 periodicities were left, while unfamiliar sounds correlated with low coincidence levels. This corresponds to an attention to pitch and timbre respectively. Additionally, separateness perception shows an opposite correlation between perception and neural correlates between familiar and unfamiliar sounds. This correlates with the perceptional finding of the distinction between familiar and unfamiliar sounds with separateness.


title: A revisit to periodic continuants

abstract: We give a simple proof of some explicit formulas of periodic continuants by Chebyshev polynomials of the second kind given by Rozsa.


title: Enhanced spectral resolution for broadband coherent anti-Stokes Raman spectroscopy

abstract: The spectral resolution of broadband Fourier-transform coherent anti-Stokes Raman spectroscopy is limited by the maximum optical path length difference that can be scanned within a short time in an interferometer. However, alternatives to the Fourier-transform exist which can bypass this limitation with certain assumptions. We apply one such approach to broadband coherent Raman spectroscopy using interferometers with short delay line (low Fourier spectral resolution) and large delay line (high Fourier spectral resolution). With this method, we demonstrate broadband coherent Raman spectroscopy of closely spaced vibrational bands is possible using a short delay line interferometer, with superior spectral resolution to the longer delay line instrument. We discuss how this approach may be particularly useful for more complex Raman spectra, such as those measured from biological samples.


title: Manifolds with two projective bundle structures

abstract: In this paper we classify varieties of Picard number two having two projective bundle structures of any relative dimension, under the assumption that these structures are mutually uniform. As an application we prove the Campana--Peternell conjecture for varieties of Picard number one admitting $\mathbb C^*$-actions of a certain kind.


title: GraphLIME: Local Interpretable Model Explanations for Graph Neural Networks

abstract: Graph structured data has wide applicability in various domains such as physics, chemistry, biology, computer vision, and social networks, to name a few. Recently, graph neural networks (GNN) were shown to be successful in effectively representing graph structured data because of their good performance and generalization ability. GNN is a deep learning based method that learns a node representation by combining specific nodes and the structural/topological information of a graph. However, like other deep models, explaining the effectiveness of GNN models is a challenging task because of the complex nonlinear transformations made over the iterations. In this paper, we propose GraphLIME, a local interpretable model explanation for graphs using the Hilbert-Schmidt Independence Criterion (HSIC) Lasso, which is a nonlinear feature selection method. GraphLIME is a generic GNN-model explanation framework that learns a nonlinear interpretable model locally in the subgraph of the node being explained. More specifically, to explain a node, we generate a nonlinear interpretable model from its $N$-hop neighborhood and then compute the K most representative features as the explanations of its prediction using HSIC Lasso. Through experiments on two real-world datasets, the explanations of GraphLIME are found to be of extraordinary degree and more descriptive in comparison to the existing explanation methods.


title: Fermi Level Controlled Ultrafast Demagnetization Mechanism in Half-Metallic Heusler Alloy

abstract: The electronic band structure-controlled ultrafast demagnetization mechanism in Co2FexMn1-xSi Heusler alloy is underpinned by systematic variation of composition. We find the spin-flip scattering rate controlled by spin density of states at Fermi level is responsible for non-monotonic variation of ultrafast demagnetization time ({\tau}M) with x with a maximum at x = 0.4. Furthermore, Gilbert damping constant exhibits an inverse relationship with {\tau}M due to the dominance of inter-band scattering mechanism. This establishes a unified mechanism of ultrafast spin dynamics based on Fermi level position.


title: Concentration phenomena in a diffusive aggregation model

abstract: We consider the drift-diffusion equation $$ u_t-\varepsilon \Delta u+\nabla\cdot(u\nabla K\star u)=0 $$ in the whole space with global-in-time bounded solutions. Mass concentration phenomena for radially symmetric solutions of this equation with small diffusivity are studied.


title: Hybrid-functional electronic structure of multilayer graphene

abstract: Multilayer graphene with rhombohedral and Bernal stacking are supposed to be metallic, as predicted by density functional theory calculations using semi-local functionals. However recent angular resolved photoemission and transport data have questioned this point of view. In particular, rhombohedral flakes are suggested to be magnetic insulators. Bernal flakes composed of an even number of layers are insulating, while those composed of an odd number of layers are pseudogapped. Here, by systematically benchmarking with plane waves codes, we develop very accurate all-electron Gaussian basis sets for graphene multilayers. We find that, in agreement with our previous calculations, rhombohedral stacked multilayer graphene are gapped for and magnetic. However, the valence band curvature and the details of the electronic structure depend crucially on the basis set. Only substantially extended basis sets are able to correctly reproduce the effective mass of the valence band top at the K point, while the popular POB-TZVP basis set leads to a severe overestimation. In the case of Bernal stacking, we show that exact exchange gaps the flakes composed by four layers and opens pseudogaps for N = 3, 6, 7, 8. However, the gap or pseudogap size and its behaviour as a function of thickness are not compatible with experimental data. Moreover, hybrid functionals lead to a metallic solution for 5 layers and a magnetic ground state for 5, 6 and 8 layers. Magnetism is very weak with practically no effect on the electronic structure and the magnetic moments are mostly concentrated in the central layers. Our hybrid functional calculations on trilayer Bernal graphene multilayers are in excellent agreement with non-magnetic GW calculations. For thicker multilayers, our calculations are a benchmark for manybody theoretical modeling of the low energy electronic structure.


title: Coulomb-Nuclear Interference:the Latest Modification

abstract: We argue that the account of Coulomb-nuclear interference (CNI) in the differential cross-section of elastic $ pp $ scattering may be easily treated without introduction of intermediate IR regularization ("photon mass"). We also indicate that the parametrization used earlier misses some terms of the second order in $ \alpha $ while it contains a superfluous term of the first order.


title: Left $3$-Engel elements in groups of exponent $60$

abstract: Let $G$ be a group and let $x\in G$ be a left $3$-Engel element of order dividing $60$. Suppose furthermore that $\langle x\rangle^{G}$ has no elements of order $8$, $9$ and $25$. We show that $x$ is then contained in the locally nilpotent radical of $G$. In particular all the left $3$-Engel elements of a group of exponent $60$ are contained in the locally nilpotent radical.


title: Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. VII. Lithium

abstract: Lithium abundances are presented for 91 dwarf and subgiant stars in the Galactic bulge. The analysis is based on line synthesis of the 7Li line at 6707 {\AA} in high-resolution spectra obtained during gravitational microlensing events, when the brightnesses of the targets were highly magnified. Our main finding is that the bulge stars at sub-solar metallicities, and that are older than about eight billion years, does not show any sign of Li production, that is, the Li trend with metallicity is flat (or even slightly declining). This indicates that no lithium was produced during the first few billion years in the history of the bulge. This finding is essentially identical to what is seen for the (old) thick disk stars in the Solar neighbourhood, and adds another piece of evidence for a tight connection between the metal-poor bulge and the Galactic thick disk. For the bulge stars younger than about eight billion years, the sample contains a group of stars at very high metallicities at [Fe/H]~+0.4 that have lithium abundances in the range A(Li)=2.6-2.8. In the Solar neighbourhood the lithium abundances have been found to peak at a A(Li)~3.3 at [Fe/H]~ +0.1 and then decrease by 0.4-0.5 dex when reaching [Fe/H]~+0.4. The few bulge stars that we have at these metallicities, seem to support this declining A(Li) trend. This could indeed support the recent claim that the low A(Li) abundances at the highest metallicities seen in the Solar neighbourhood could be due to stars from the inner disk, or the bulge region, that have migrated to the Solar neighbourhood.


title: Dynamics and Spreading Speed of a Reaction-Diffusion System with Advection Modeling West Nile Virus

abstract: This paper aims to explore the temporal-spatial spreading and asymptotic behaviors of West Nile virus by a reaction-advection-diffusion system with free boundaries, especially considering the impact of advection term on the extinction and persistence of West Nile virus. We define the spatial-temporal risk index $R^{F}_{0}(t)$ with the advection rate and the general basic disease reproduction number $R^D_0$ to get the vanishing-spreading dichotomy regimes of West Nile virus. We show that there exists a threshold value $\mu^{*}$ of the advection rate, and obtain the threshold results of it. When the spreading occurs, we investigate the asymptotic dynamical behaviors of the solution in the long run and first give a sharper estimate that the asymptotic spreading speed of the leftward front is less than the rightward front for $0<\mu<\mu^*$. At last, we give some numerical simulations to identify the significant effects of the advection.


title: Fragmentation favoured in discs around higher mass stars

abstract: We investigate how a protoplanetary disc's susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increase the host star mass. The actual limit for stability is sensitive to the disc temperature, so if the disc is optically thin stellar irradiation can dramatically stabilise discs against gravitational instability. However, even when this is the case we find that discs around $2$M$_{\odot}$ stars are prone to fragmentation, which will act to produce wide-orbit giant planets and brown dwarfs. The consequences of this work are two-fold: that low mass stars could in principle support high disc-to-star mass ratios, and that higher mass stars have discs that are more prone to fragmentation, which is qualitatively consistent with observations that favour high-mass wide-orbit planets around higher mass stars. We also find that the initial masses of these planets depends on the temperature in the disc at large radii, which itself depends on the level of stellar irradiation.


title: Massive discs around low-mass stars

abstract: We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.