title: On Parallel Solution of Sparse Triangular Linear Systems in CUDA

abstract: The acceleration of sparse matrix computations on modern many-core processors, such as the graphics processing units (GPUs), has been recognized and studied over a decade. Significant performance enhancements have been achieved for many sparse matrix computational kernels such as sparse matrix-vector products and sparse matrix-matrix products. Solving linear systems with sparse triangular structured matrices is another important sparse kernel as demanded by a variety of scientific and engineering applications such as sparse linear solvers. However, the development of efficient parallel algorithms in CUDA for solving sparse triangular linear systems remains a challenging task due to the inherently sequential nature of the computation. In this paper, we will revisit this problem by reviewing the existing level-scheduling methods and proposing algorithms with self-scheduling techniques. Numerical results have indicated that the CUDA implementations of the proposed algorithms can outperform the state-of-the-art solvers in cuSPARSE by a factor of up to $2.6$ for structured model problems and general sparse matrices.


title: Characterization and Enumeration of Complementary Dual Abelian Codes

abstract: Abelian codes and complementary dual codes form important classes of linear codes that have been extensively studied due to their rich algebraic structures and wide applications. In this paper, a family of abelian codes with complementary dual in a group algebra $\mathbb{F}_{p^\nu}[G]$ has been studied under both the Euclidean and Hermitian inner products, where $p$ is a prime, $\nu$ is a positive integer, and $G$ is an arbitrary finite abelian group. Based on the discrete Fourier transform decomposition for semi-simple group algebras and properties of ideas in local group algebras, the characterization of such codes have been given. Subsequently, the number of complementary dual abelian codes in $\mathbb{F}_{p^\nu}[G]$ has been shown to be independent of the Sylow $p$-subgroup of $G$ and it has been completely determined for every finite abelian group $G$. In some cases, a simplified formula for the enumeration has been provided as well. The known results for cyclic complementary dual codes can be viewed as corollaries.


title: Simulating the Ising Model with a Deep Convolutional Generative Adversarial Network

abstract: The deep learning framework is witnessing expansive growth into diverse applications such as biological systems, human cognition, robotics, and the social sciences, thanks to its immense ability to extract essential features from complicated systems. In particular, recent developments of the field have revealed the unique faculty of deep learning to accurately approximate complex physical systems in fluid dynamics, condensed matter physics, etc. The convolutional neural network (CNN) is an efficient approach to represent complex systems with large degrees of freedom. On the other hand, the generative adversarial network (GAN), as an unsupervised learning algorithm, is capable of efficiently imitating the distribution of training data. Here we leverage this unique property of GAN, in conjunction with CNN methodology, to establish an Ising simulator whose generator can produce Ising states given temperature T around the criticality. The generated Ising states well resemble, and essentially replicate, the data from conventional Monte Carlo simulations. Our results demonstrate the universality of GAN as a promising tool in the field of computational and statistical physics.


title: Fermi-momentum dependence of relativistic effective mass below saturation from superscaling of quasielastic electron scattering

abstract: The relativistic effective mass $M^*$, and Fermi momentum, $k_F$, are important ingredients in the determination of the nuclear equation of state, but they have rarely been extracted from experimental data below saturation density where translationally invariant nuclear matter becomes unstable against clusterization into the existing atomic nuclei. Using a novel kind of superscaling analysis of the quasielastic cross section electron scattering data involving a suitable selection criterion and $^{12}$C as a reference nucleus, the global scaling properties of the resulting set of data for 21 nuclei ranging from $^2$H to $^{238}$U are then analyzed. We find that a subset of a third of the about $20000$ data approximately scales to an universal superscaling function with a more constrained uncertainty band than just the reference $^{12}$C case and provides $M^*$ as a function of $k_F$.


title: Complex Word Identification: Challenges in Data Annotation and System Performance

abstract: This paper revisits the problem of complex word identification (CWI) following up the SemEval CWI shared task. We use ensemble classifiers to investigate how well computational methods can discriminate between complex and non-complex words. Furthermore, we analyze the classification performance to understand what makes lexical complexity challenging. Our findings show that most systems performed poorly on the SemEval CWI dataset, and one of the reasons for that is the way in which human annotation was performed.


title: A Revision to the Issue of Frames by Non-minimal Large Field Inflation

abstract: We present an extended study of inflationary models that inflaton field is non-minimally coupled with gravity. We study parameters space of the models up to the second (and in some cases third) order of the slow-roll parameters for usual large field potentials in Jordan and Einstein frames that are connected each other by conformal transformation. We calculate inflationary parameters and the results are compared in both frames and also with observations. By using the recent observational datasets, we present a discussion in order to clarify the physical frame between Jordan and Einstein frames. Also, some suggestions are expressed in order to navigate us for the future works.


title: An NFV and Microservice Based Architecture for On-the-fly Component Provisioning in Content Delivery Networks

abstract: Content Delivery Networks (CDNs) deliver content (e.g. Web pages, videos) to geographically distributed end-users over the Internet. Some contents do sometimes attract the attention of a large group of end-users. This often leads to flash crowds which can cause major issues such as outage in the CDN. Microservice architectural style aims at decomposing monolithic systems into smaller components which can be independently deployed, upgraded and disposed. Network Function Virtualization (NFV) is an emerging technology that aims to reduce costs and bring agility by decoupling network functions from the underlying hardware. This paper leverages the NFV and microservice architectural style to propose an architecture for on-the-fly CDN component provisioning to tackle issues such as flash crowds. In the proposed architecture, CDN components are designed as sets of microservices which interact via RESTFul Web services and are provisioned as Virtual Network Functions (VNFs), which are deployed and orchestrated on-the-fly. We have built a prototype in which a CDN surrogate server, designed as a set of microservices, is deployed on-the-fly. The prototype is deployed on SAVI, a Canadian distributed test bed for future Internet applications. The performance is also evaluated.


title: Statistical symmetry restoration in fully developed turbulence: Renormalization group analysis of two models

abstract: In this paper we consider the model of incompressible fluid described by the stochastic Navier-Stokes equation with finite correlation time of a random force. Inertial-range asymptotic behavior of fully developed turbulence is studied by means of the field theoretic renormalization group within the one-loop approximation. It is corroborated that regardless of the values of model parameters and initial data, the inertial-range behavior of the model is described by limiting case of vanishing correlation time. It indicates that the Galilean symmetry of the model violated by the "colored" random force is restored in the inertial range. This regime corresponds to the only nontrivial fixed point of the renormalization group equation. The stability of this point depends on the relation between the exponents in the energy spectrum ${\cal E} \propto k^{1-y}$ and the dispersion law $\omega \propto k^{2-\eta}$. The second analyzed problem is the passive advection of a scalar field by this velocity ensemble. Correlation functions of the scalar field exhibit anomalous scaling behavior in the inertial-convective range. We demonstrate that in accordance with Kolmogorov's hypothesis of the local symmetry restoration, the main contribution to the operator product expansion is given by the isotropic operator, while anisotropic terms should be considered only as corrections.


title: Scattering of kinks of the sinh-deformed $\varphi^4$ model

abstract: We consider the scattering of kinks of the sinh-deformed $\varphi^4$ model, which is obtained from the well-known $\varphi^4$ model by means of the deformation procedure. Depending on the initial velocity $v_{in}$ of the colliding kinks, different collision scenarios are realized. There is a critical value $v_{cr}$ of the initial velocity, which separates the regime of reflection (at $v_{in}>v_{cr}$) and that of a complicated interaction (at $v_{in}<v_{cr}$) with kinks' capture and escape windows. Besides that, at $v_{in}$ below $v_{cr}$ we observe the formation of a bound state of two oscillons, as well as their escape at some values of $v_{in}$.

categories: hep-th cond-mat.mtrl-sci math-ph math.MP nlin.PS


title: Circumstellar Interaction Models for the Bolometric Light Curve of SN 2017egm

abstract: We explore simple semi-analytic fits to the bolometric light curve of Gaia17biu/SN 2017egm, the most nearby hydrogen-deficient superluminous supernova (SLSN I) yet discovered. SN 2017egm has a quasi-bolometric light curve that is uncharacteristic of other SLSN I by having a nearly linear rise to maximum and decline from peak, with a very sharp transition. Magnetar models have difficulty explaining the sharp peak and may tend to be too bright 20 d after maximum. Light curves powered only by radioactive decay of nickel fail on similar grounds and because they demand greater nickel mass than ejecta mass. Simple models based on circumstellar interaction do have a sharp peak corresponding to the epoch when the forward shock breaks out of the optically-thick circumstellar medium or the reverse shock reaches the inside of the ejecta. We find that models based on circumstellar interaction with a constant-density shell provide an interesting fit to the bolometric light curve from 15 d before to 15 d after peak light of SN 2017egm and that both magnetar and radioactive decay models fail to fit the sharp peak. Future photometric observations should easily discriminate basic CSI models from basic magnetar models. The implications of a CSI model are briefly discussed.


title: Multiple Equivalent Solutions for the Lasso

abstract: Feature selection is an important problem studied in data analytics seeking to identify a minimal-size feature subset that is optimally predictive for an outcome of interest. It is also a powerful tool in Knowledge Discovery as a means for gaining domain insight, e.g., identifying which medical quantities carry unique information for the disease status. It is arguably less recognized however, that the problem may have multiple, equivalent solutions. In that case, it is misleading to domain experts to report only one of them and ignore all other equivalent solutions. In this paper, we extend a well-established single, feature selection algorithm (i.e., reporting a single solution), namely the Lasso algorithm, to the multiple solution problem based on formalized notion of equivalence for both classification and regression tasks. Empirical results are obtained using a fully automated pipeline called Just Add Data Bio or JAD Bio training and selecting multiple, linear as well as nonlinear learners, optimizing hyper-parameter values, and correcting for the bias of multiple inductions (model selection). The results show that multiple solutions do exist in real datasets, as well as the ability of the algorithm to identify a subset of them. A comparison with the Statistical Equivalent Solutions (SES) algorithm shows that Lasso equivalent solutions have better prediction performance at the cost of selecting more features.


title: Merger delay time distribution of extended emission short GRBs

abstract: The most popular progenitor model for short duration Gamma-Ray bursts (sGRBs) is the merger of two compact objects. However, the short GRB population exhibit a certain diversity: some bursts display an extended emission (EE), continuing in soft $\gamma$-rays for a few hundreds of seconds post the initial short pulse. It is currently unclear whether the origin of such bursts is linked to compact object mergers. Within the merger hypothesis, the redshift ($z$) distribution of short GRBs is influenced by the merger delay time, i.e., time elapsed between the merger and the formation of the binary star system, which is dominated by the time-scale for gravitational wave losses during the compact binary phase. We examine redshift distributions of short GRBs with extended emission to see whether their formation channel requires considerable delay post the star formation episode. Our results show that the $z$ distribution of EE bursts is consistent with the merger model. We attempted to compare the delay time distribution of the EE and the non-EE short bursts. However, no statistically significant difference could be seen within the limited sample size.


title: The search for radio emission from exoplanets using LOFAR low-frequency beam-formed observations: Data pipeline and preliminary results for the 55 Cnc system

abstract: Detection of radio emission from exoplanets can provide information on the star-planet system that is difficult to study otherwise, such as the planetary magnetic field, magnetosphere, rotation period, interior structure, atmospheric dynamics and escape, and any star-planet interactions. Such a detection in the radio domain would open up a whole new field in the study of exoplanets. However, currently there are no confirmed detections of an exoplanet at radio frequencies. In this study, we search for non-thermal radio emission from the 55 Cnc system which has 5 known exoplanets. According to theoretical predictions 55 Cnc e, the innermost planet, is among the best targets for this search. We observed for 18 hours with the Low-Frequency Array (LOFAR) Low Band Antenna in the frequency range 26-73 MHz with full-polarization and covered 85% of the orbital phase of 55 Cnc e. During the observations four digital beams within the station beam were recorded simultaneously on 55 Cnc, nearby "empty" sky, a bright radio source, and a pulsar. A pipeline was created to automatically find and mask radio frequency interference, calibrate the time-frequency response of the telescope, and to search for bursty planetary radio signals in our data. Extensive tests and verifications were carried out on the pipeline. Analysis of the first 4 hours of these observations do not contain any exoplanet signal from 55 Cnc but we can confirm that our setup is adequate to detect faint astrophysical signals. We find a 3$\sigma$ upper limit for 55 Cnc of 230 mJy using the pulsar to estimate the sensitivity of the observations and 2.6 Jy using the time-series difference between the target and sky beam. The full data set is still under-going analysis. In the near future we will apply our observational technique and pipeline to the most promising exoplanet candidates for which LOFAR observations have already been obtained.


title: QCD calculations for the LHC: status and prospects

abstract: We briefly review the status of high precision QCD predictions available for LHC processes, focusing on corrections beyond NLO and ways to make the latter available in a convenient format. As a phenomenological example we discuss the two-loop corrections to Higgs boson pair production in gluon fusion.


title: Matrix-Product Constructions for Hermitian Self-Orthogonal Codes

abstract: Self-orthogonal codes have been of interest due to there rich algebraic structures and wide applications. Euclidean self-orthogonal codes have been quite well studied in literature. Here, we have focused on Hermitian self-orthogonal codes. Constructions of such codes have been given based on the well-known matrix-product construction for linear codes. Criterion for the underlying matrix and the input codes required in such constructions have been determined. In many cases, the Hermitian self-orthogonality of the input codes and the assumption that the underlying matrix is unitary can be relaxed. Some special matrices used in the constructions and illustrative examples of good Hermitian self-orthogonal codes have been provided as well.


title: The year-long flux variations in Boyajian's star are asymmetric or aperiodic

abstract: We combine and calibrate publicly available data for Boyajian's star including photometry from ASAS (SN, V, I), Kepler, Gaia, SuperWASP, and citizen scientist observations (AAVSO, HAO and Burke-Gaffney). Precise (mmag) photometry covers the years 2006-2017. We show that the year-long flux variations with an amplitude of ~4% can not be explained with cyclical symmetric or asymmetric models with periods shorter than ten years. If the dips are transits, their period must exceed ten years, or their structure must evolve significantly during each 4-year long cycle.


title: Transformation formulas of a character analogue of $\log\theta_{2}(z)$

abstract: In this paper, transformation formulas for the function \[ A_{1}\left(z,s:\chi\right)=\sum\limits_{n=1}^{\infty}\sum\limits_{m=1}^{\infty}\chi\left(n\right)\chi\left(m\right)\left(-1\right)^{m}n^{s-1}e^{2\pi imnz/k} \] are obtained. Sums that appear in transformation formulas are generalizations of the Hardy--Berndt sums $s_{j}(d,c),$ $j=1,2,5$. As applications of these transformation formulas, reciprocity formulas for these sums are derived and several series relations are presented.


title: A deformation of instanton homology for webs

abstract: A deformation of the authors' instanton homology for webs is constructed by introducing a local system of coefficients. In the case that the web is planar, the rank of the deformed instanton homology is equal to the number of Tait colorings of the web.


title: An FBAR Circulator

abstract: This letter presents the experimental demonstration of a film bulk acoustic resonator (FBAR) circulator at 2.5 GHz. The circulator is based on spatio-temporal modulation of the series resonant frequency of FBARs using varactors and exhibits a large isolation of 76 dB at 2.5 GHz. The FBAR chip (0.25 mm2) consists of three identical FBARs connected in wye configuration. The FBAR0s quality factor (Q) of 1250 and piezoelectric coupling coefficient kt 2 of 3% relaxes the modulation requirements, achieving non-reciprocity with small modulationto- RF frequency ratio bettter than 1:800 (3 MHz:2.5 GHz).


title: Possible quadrupolar nematic phase in the frustrated spin chain LiCuSbO$_4$: an NMR investigation

abstract: The frustrated one-dimensional (1D) quantum magnet LiCuSbO$_4$ is one rare realization of the $J_1-J_2$ spin chain model with an easily accessible saturation field, formerly estimated to 12~T. Exotic multipolar nematic phases were theoretically predicted in such compounds just below the saturation field, but without unambiguous experimental observation so far. In this paper we present extensive experimental research of the compound in the wide temperature (30mK$-$300K) and field (0$-$13.3T) range by muon spin rotation ($\mu$SR), $^7$Li nuclear magnetic resonance (NMR) and magnetic susceptibility (SQUID). $\mu$SR experiments in zero magnetic field demonstrate the absence of long range 3D ordering down to 30mK. Together with former heat capacity data [S.E. Dutton \emph{et al}, Phys. Rev. Lett. 108, 187206 (2012)], magnetic susceptibility measurements suggest short range correlated vector chiral phase in the field range $0-4$T. In the intermediate field values (5$-$12T), the system enters in a 3D ordered spin density wave phase with 0.75$\mu_B$ per copper site at lowest temperatures (125mK), estimated by NMR. At still higher field, the magnetization is found to be saturated above 13T where the spin lattice $T_1^{-1}$ relaxation reveals a spin gap estimated to 3.2(2)K. We narrow down the possibility of observing a multipolar nematic phase to the range 12.5$-$13T.


title: First Principles Free-Energy Theory of Solvation with Atomic Scale Liquid Structure

abstract: Quantum-chemical processes in liquid environments impact broad areas of science, from molecular biology to geology to electrochemistry. While density-functional theory (DFT) has enabled efficient quantum-mechanical calculations which profoundly impact understanding of atomic-scale phenomena, realistic description of the liquid remains a challenge. Here, we present an approach based on joint density-functional theory (JDFT) which addresses this challenge by leveraging the DFT approach not only for the quantum mechanics of the electrons in a solute, but also simultaneously for the statistical mechanics of the molecules in a surrounding equilibrium liquid solvent. Specifically, we develop a new universal description for the interaction of electrons with an arbitrary liquid, providing the missing link to finally transform JDFT into a practical tool for the realistic description of chemical processes in solution. This approach predicts accurate solvation free energies and surrounding atomic-scale liquid structure for molecules and surfaces in multiple solvents without refitting, all at a fraction of the computational cost of methods of comparable detail and accuracy. To demonstrate the potential impact of this method, we determine the structure of the solid/liquid interface, offering compelling agreement with more accurate (but much more computationally intensive) theories and with X-ray reflectivity measurements.


title: Skin Lesion Analysis Toward Melanoma Detection: A Challenge at the 2017 International Symposium on Biomedical Imaging (ISBI), Hosted by the International Skin Imaging Collaboration (ISIC)

abstract: This article describes the design, implementation, and results of the latest installment of the dermoscopic image analysis benchmark challenge. The goal is to support research and development of algorithms for automated diagnosis of melanoma, the most lethal skin cancer. The challenge was divided into 3 tasks: lesion segmentation, feature detection, and disease classification. Participation involved 593 registrations, 81 pre-submissions, 46 finalized submissions (including a 4-page manuscript), and approximately 50 attendees, making this the largest standardized and comparative study in this field to date. While the official challenge duration and ranking of participants has concluded, the dataset snapshots remain available for further research and development.


title: Nonlinear Split Ordered Variational Inequality Problems

abstract: The concept of nonlinear split ordered variational inequality problems on partially ordered vector spaces is a natural extension of linear split vector variational inequality problems on Banach spaces. The results about nonlinear split ordered variational inequality problems are immediately applied to solving nonlinear split vector optimization problems. In this paper, we prove the solvability of some nonlinear split ordered variational inequality problems by using some fixed point theorems on partially ordered spaces, in which the considered mappings may not be required to have any type of continuity and they just satisfy some order-monotonic conditions. As applications of the results about nonlinear split ordered variational inequality problems, we study solvability of some nonlinear split variational inequality problems and linear split variational inequality problems on partially ordered vector spaces and partially ordered Banach spaces.


title: Automatic Detection and Uncertainty Quantification of Landmarks on Elastic Curves

abstract: A population quantity of interest in statistical shape analysis is the location of landmarks, which are points that aid in reconstructing and representing shapes of objects. We provide an automated, model-based approach to inferring landmarks given a sample of shape data. The model is formulated based on a linear reconstruction of the shape, passing through the specified points, and a Bayesian inferential approach is described for estimating unknown landmark locations. The question of how many landmarks to select is addressed in two different ways: (1) by defining a criterion-based approach, and (2) joint estimation of the number of landmarks along with their locations. Efficient methods for posterior sampling are also discussed. We motivate our approach using several simulated examples, as well as data obtained from applications in computer vision and biology; additionally, we explore placements and associated uncertainty in landmarks for various substructures extracted from magnetic resonance image slices.


title: Trapping and Escape in a Turbid Medium

abstract: We investigate the absorption of diffusing molecules in a fluid-filled spherical beaker that contains many small reactive traps. The molecules are absorbed either by hitting a trap or by escaping via the beaker walls. In the physical situation where the number $N$ of traps is large and their radii $a$ are small compared to the beaker radius $R$, the fraction of molecules $E$ that escape to the beaker wall and the complementary fraction $T$ that eventually are absorbed by the traps depend only on the dimensionless parameter combination $\lambda = Na/R$. We compute $E$ and $T$ as a function of $\lambda$ for a spherical beaker and for beakers of other three-dimensional shapes. The asymptotic behavior is found to be universal: $1- E\sim \lambda$ for $\lambda\to 0$ and $E\sim\lambda^{-1/2}$ for $\lambda\to\infty$.


title: Supersymmetric solutions and Borel singularities for N=2 supersymmetric Chern-Simons theories

abstract: In supersymmetric (SUSY) field theory, there exist configurations which formally satisfy SUSY conditions but are not on original path integral contour. We refer to such configurations as complexified supersymmetric solutions (CSS). In this paper we discuss that CSS provide important information on large order behavior of weak coupling perturbative series in SUSY field theories. We conjecture that CSS with a bosonic (fermionic) free parameter give poles (zeroes) of Borel transformation of perturbative series whose locations are uniquely determined by actions of the solutions. We demonstrate this for various SUSY observables in 3d $\mathcal{N}=2$ SUSY Chern-Simons matter theories on sphere. First we construct infinite number of CSS in general 3d $\mathcal{N}=2$ SUSY theory with Lagrangian where adjoint scalar in vector multiplet takes a complex value and matter fields are nontrivial. Then we compare their actions with Borel transformations of perturbative expansions by inverse Chern-Simons levels for the observables and see agreement with our conjecture. It turns out that the CSS explain all the Borel singularities for this case.

categories: hep-th cond-mat.str-el hep-ph math-ph math.MP


title: Fast high-dimensional node generation with variable density

abstract: We present an algorithm for producing discrete distributions with a prescribed nearest-neighbor distance function. Our approach is a combination of quasi-Monte Carlo (Q-MC) methods and weighted Riesz energy minimization: the initial distribution is a stratified Q-MC sequence with some modifications; a suitable energy functional on the configuration space is then minimized to ensure local regularity. The resulting node sets are good candidates for building meshless solvers and interpolants, as well as for other purposes where a point cloud with a controlled separation-covering ratio is required. Applications of a three-dimensional implementation of the algorithm, in particular to atmospheric modeling, are also given.


title: Potential Conditional Mutual Information: Estimators, Properties and Applications

abstract: The conditional mutual information I(X;Y|Z) measures the average information that X and Y contain about each other given Z. This is an important primitive in many learning problems including conditional independence testing, graphical model inference, causal strength estimation and time-series problems. In several applications, it is desirable to have a functional purely of the conditional distribution p_{Y|X,Z} rather than of the joint distribution p_{X,Y,Z}. We define the potential conditional mutual information as the conditional mutual information calculated with a modified joint distribution p_{Y|X,Z} q_{X,Z}, where q_{X,Z} is a potential distribution, fixed airport. We develop K nearest neighbor based estimators for this functional, employing importance sampling, and a coupling trick, and prove the finite k consistency of such an estimator. We demonstrate that the estimator has excellent practical performance and show an application in dynamical system inference.


title: A Case Study Competition Among Methods for Analyzing Large Spatial Data

abstract: The Gaussian process is an indispensable tool for spatial data analysts. The onset of the "big data" era, however, has lead to the traditional Gaussian process being computationally infeasible for modern spatial data. As such, various alternatives to the full Gaussian process that are more amenable to handling big spatial data have been proposed. These modern methods often exploit low rank structures and/or multi-core and multi-threaded computing environments to facilitate computation. This study provides, first, an introductory overview of several methods for analyzing large spatial data. Second, this study describes the results of a predictive competition among the described methods as implemented by different groups with strong expertise in the methodology. Specifically, each research group was provided with two training datasets (one simulated and one observed) along with a set of prediction locations. Each group then wrote their own implementation of their method to produce predictions at the given location and each which was subsequently run on a common computing environment. The methods were then compared in terms of various predictive diagnostics. Supplementary materials regarding implementation details of the methods and code are available for this article online.


title: Tropicalization of Positive Grassmannians

abstract: We introduce combinatorial objects which are parameterized by the positive part of the tropical Grassmannian $Gr(k,n)$. Our method is to relate the Grassmannian to configuration spaces of flags. By work of the first author, and of Goncharov and Shen, configuration spaces of flags naturally tropicalize to give configurations of points in the affine building, which we call higher laminations. We use higher laminations to give two dual objects that are parameterized by the positive tropicalization of $Gr(k,n)$: equivalence classes of higher laminations; or certain restricted subset of higher laminations. This extends results of Speyer and Sturmfels on the tropicalization of $Gr(2,n)$, and of Speyer and Williams on the tropicalization of $Gr(3,6)$ and $Gr(3,7)$. We also analyze the $X$-variety associated to the Grassmannian, and give an interpretation of its positive tropicalization.


title: Coherence of quantum channels

abstract: We investigate the coherence of quantum channels using the Choi-Jamio\l{}kowski isomorphism. The relation between the coherence and the purity of the channel respects a duality relation. It characterizes the allowed values of coherence when the channel has certain purity. This duality has been depicted via the Coherence-Purity (Co-Pu) diagrams. In particular, we study the quantum coherence of the unital and non-unital qubit channels and find out the allowed region of coherence for a fixed purity. We also study coherence of different incoherent channels, namely, incoherent operation (IO), strictly incoherent operation (SIO), physical incoherent operation (PIO) etc. Interestingly, we find that the allowed region for different incoherent operations maintain the relation $PIO\subset SIO \subset IO$. In fact, we find that if PIOs are coherence preserving operations (CPO), its coherence is zero otherwise it has unit coherence and unit purity. Interestingly, different kinds of qubit channels can be distinguished using the Co-Pu diagram. The unital channels generally do not create coherence whereas some nonunital can. All coherence breaking channels are shown to have zero coherence, whereas, this is not usually true for entanglement breaking channels. It turns out that the coherence preserving qubit channels have unit coherence. Although the coherence of the Choi matrix of the incoherent channels might have finite values, its subsystem contains no coherence. This indicates that the incoherent channels can either be unital or nonunital under some conditions.


title: Census of the Local Universe (CLU) Narrow-Band Survey I: Galaxy Catalogs from Preliminary Fields

abstract: We present the Census of the Local Universe (CLU) narrow-band survey to search for emission-line (\ha) galaxies. CLU-\ha~has imaged $\approx$3$\pi$ of the sky (26,470~deg$^2$) with 4 narrow-band filters that probe a distance out to 200~Mpc. We have obtained spectroscopic follow-up for galaxy candidates in 14 preliminary fields (101.6~deg$^2$) to characterize the limits and completeness of the survey. In these preliminary fields, CLU can identify emission lines down to an \ha~flux limit of $10^{-14}$~$\rm{erg~s^{-1}~cm^{-2}}$ at 90\% completeness, and recovers 83\% (67\%) of the \ha~flux from catalogued galaxies in our search volume at the $\Sigma$=2.5 ($\Sigma$=5) color excess levels. The contamination from galaxies with no emission lines is 61\% (12\%) for $\Sigma$=2.5 ($\Sigma$=5). Also, in the regions of overlap between our preliminary fields and previous emission-line surveys, we recover the majority of the galaxies found in previous surveys and identify an additional $\approx$300 galaxies. In total, we find 90 galaxies with no previous distance information, several of which are interesting objects: 7 blue compact dwarfs, 1 green pea, and a Seyfert galaxy; we also identified a known planetary nebula. These objects show that the CLU-\ha~survey can be a discovery machine for objects in our own Galaxy and extreme galaxies out to intermediate redshifts. However, the majority of the CLU-\ha~galaxies identified in this work show properties consistent with normal star-forming galaxies. CLU-\ha~galaxies with new redshifts will be added to existing galaxy catalogs to focus the search for the electromagnetic counterpart to gravitational wave events.


title: The power of sum-of-squares for detecting hidden structures

abstract: We study planted problems---finding hidden structures in random noisy inputs---through the lens of the sum-of-squares semidefinite programming hierarchy (SoS). This family of powerful semidefinite programs has recently yielded many new algorithms for planted problems, often achieving the best known polynomial-time guarantees in terms of accuracy of recovered solutions and robustness to noise. One theme in recent work is the design of spectral algorithms which match the guarantees of SoS algorithms for planted problems. Classical spectral algorithms are often unable to accomplish this: the twist in these new spectral algorithms is the use of spectral structure of matrices whose entries are low-degree polynomials of the input variables. We prove that for a wide class of planted problems, including refuting random constraint satisfaction problems, tensor and sparse PCA, densest-k-subgraph, community detection in stochastic block models, planted clique, and others, eigenvalues of degree-d matrix polynomials are as powerful as SoS semidefinite programs of roughly degree d. For such problems it is therefore always possible to match the guarantees of SoS without solving a large semidefinite program. Using related ideas on SoS algorithms and low-degree matrix polynomials (and inspired by recent work on SoS and the planted clique problem by Barak et al.), we prove new nearly-tight SoS lower bounds for the tensor and sparse principal component analysis problems. Our lower bounds for sparse principal component analysis are the first to suggest that going beyond existing algorithms for this problem may require sub-exponential time.


title: The distribution of symmetry of a naturally reductive nilpotent Lie group

abstract: We show that the distribution of symmetry of a naturally reductive nilpotent Lie group coincides with the invariant distribution induced by the set of fixed vectors of the isotropy. This extends a known result on compact naturally reductive spaces. We also address the study of the quotient by the foliation of symmetry.


title: Tensor product decompositions of II$_1$ factors arising from extensions of amalgamated free product groups

abstract: In this paper we introduce a new family of icc groups $\Gamma$ which satisfy the following product rigidity phenomenon, discovered in [DHI16] (see also [dSP17]): all tensor product decompositions of the II$_1$ factor $L(\Gamma)$ arise only from the canonical direct product decompositions of the underlying group $\Gamma$. Our groups are assembled from certain HNN-extensions and amalgamated free products and include many remarkable groups studied throughout mathematics such as graph product groups, poly-amalgam groups, Burger-Mozes groups, Higman group, various integral two-dimensional Cremona groups, etc. As a consequence, we obtain several new examples of groups that give rise to prime factors.


title: Equation of state of non-relativistic matter from automated perturbation theory and complex Langevin

abstract: We calculate the pressure and density of polarized non-relativistic systems of two-component fermions coupled via a contact interaction at finite temperature. For the unpolarized one-dimensional system with an attractive interaction, we perform a third-order lattice perturbation theory calculation and assess its convergence by comparing with hybrid Monte Carlo. In that regime, we also demonstrate agreement with real Langevin. For the repulsive unpolarized one-dimensional system, where there is a so-called complex phase problem, we present lattice perturbation theory as well as complex Langevin calculations. For our studies, we employ a Hubbard-Stratonovich transformation to decouple the interaction and automate the application of Wick's theorem for perturbative calculations, which generates the diagrammatic expansion at any order. We find excellent agreement between the results from our perturbative calculations and stochastic studies in the weakly interacting regime. In addition, we show predictions for the strong coupling regime as well as for the polarized one-dimensional system. Finally, we show a first estimate for the equation of state in three dimensions where we focus on the polarized unitary Fermi gas.


title: Simultaneous Detection of Signal Regions Using Quadratic Scan Statistics With Applications in Whole Genome Association Studies

abstract: We consider in this paper detection of signal regions associated with disease outcomes in whole genome association studies. Gene- or region-based methods have become increasingly popular in whole genome association analysis as a complementary approach to traditional individual variant analysis. However, these methods test for the association between an outcome and the genetic variants in a pre-specified region, e.g., a gene. In view of massive intergenic regions in whole genome sequencing (WGS) studies, we propose a computationally efficient quadratic scan (Q-SCAN) statistic based method to detect the existence and the locations of signal regions by scanning the genome continuously. The proposed method accounts for the correlation (linkage disequilibrium) among genetic variants, and allows for signal regions to have both causal and neutral variants, and the effects of signal variants to be in different directions. We study the asymptotic properties of the proposed Q-SCAN statistics. We derive an empirical threshold that controls for the family-wise error rate, and show that under regularity conditions the proposed method consistently selects the true signal regions. We perform simulation studies to evaluate the finite sample performance of the proposed method. Our simulation results show that the proposed procedure outperforms the existing methods, especially when signal regions have causal variants whose effects are in different directions, or are contaminated with neutral variants. We illustrate Q-SCAN by analyzing the WGS data from the Atherosclerosis Risk in Communities (ARIC) study.


title: A Grassmann and graded approach to coboundary Lie bialgebras, their classification, and Yang-Baxter equations

abstract: This work pioneers the systematic study and classification (up to Lie algebra automorphisms) of finite-dimensional coboundary Lie bialgebras through Grassmann algebras. Several mathematical structures on Lie algebras, e.g. Killing forms or root decompositions, are extended to the Grassmann algebras of Lie algebras. This simplifies the description of the procedures and tools appearing in the theory of Lie bialgebras and originates novel techniques for its study and classification up to Lie algebra automorphisms. As a particular case, the classification of real three-dimensional coboundary Lie bialgebras is retrieved.


title: Agnostic Stacking of Intergalactic Doublet Absorption: Measuring the NeVIII Population

abstract: We present a blind search for doublet intergalactic metal absorption with a method dubbed `agnostic stacking'. Using a forward-modelling framework we combine this with direct detections in the literature to measure the overall metal population. We apply this novel approach to the search for NeVIII absorption in a set of 26 high-quality COS spectra. We probe to an unprecedented low limit of log N$>$12.3 at 0.47$\leq z \leq$1.34 over a pathlength $\Delta$z = 7.36. This method selects apparent absorption without requiring knowledge of its source. Stacking this mixed population dilutes doublet features in composite spectra in a deterministic manner, allowing us to measure the proportion corresponding to NeVIII absorption. We stack potential NeVIII absorption in two regimes: absorption too weak to be significant in direct line studies (12.3 $<$ log N $<$ 13.7), and strong absorbers (log N $>$ 13.7). We do not detect NeVIII absorption in either regime. Combining our measurements with direct detections, we find that the NeVIII population is reproduced with a power law column density distribution function with slope $\beta = -1.86 \substack{+0.18 \\ -0.26}$ and normalisation log $f_{13.7} = -13.99 \substack{+0.20 \\ -0.23}$, leading to an incidence rate of strong NeVIII absorbers $dn/dz =1.38 \substack{+0.97 \\ -0.82}$. We infer a cosmic mass density for NeVIII gas with 12.3 $<$ log N $<$ 15.0 of $\Omega _{NeVIII} = 2.2 \substack{+1.6 \\ _-1.2} \times 10^{-8}$, a value significantly lower that than predicted by recent simulations. We translate this density into an estimate of the baryon density $\Omega _{b} \approx 1.8 \times 10^{-3}$, constituting 4\% of the total baryonic mass.


title: Equation of State in 2+1 Flavor QCD at High Temperatures

abstract: We calculate the Equation of State at high temperatures in 2+1 flavor QCD using the highly improved staggered quark action. We study the lattice spacing dependence of the pressure at high temperatures using lattices with temporal extent $N_{\tau}=6,~8,~10$ and $12$ and perform continuum extrapolations. We also give a continuum estimate for the Equation of State up to temperatures $T=2$ GeV, which are then compared with results of the weak-coupling calculations. We find a reasonably good agreement with the weak-coupling calculations at the highest temperatures.


title: Connecting traces of galaxy evolution: the missing core mass -- morphological fine structure relation

abstract: Deep exposure imaging of early-type galaxies (ETGs) are revealing the second-order complexity of these objects, which have been long considered uniform, dispersion-supported spheroidals. "Fine structure" features (e.g. ripples, plumes, tidal tails, rings) as well as depleted stellar cores (i.e. central light deficits) characterize a number of massive ETG galaxies, and can be interpreted as the result of galaxy-galaxy interactions. We discuss how the timescale for the evolution of cores and fine structures are comparable, and hence it is expected that they develop in parallel after the major interaction event which shaped the ETG. Using archival data, we compare the "depleted stellar mass" (i.e. the mass missing from the depleted stellar core) against the prominence of the fine structure features, and observe that they correlate inversely. This result confirms our expectation that, while the Super Massive Black Hole (SMBH) binary (constituted by the SMBHs of the merger progenitors) excavates the core via three-body interactions, the gravitational potential of the newborn galaxy relaxes, and the fine structures fade below detection levels. We expect the inverse correlation to hold at least within the first Gyr from the merger which created the SMBH binary; after then, the fine structure evolves independently.


title: Neutrinoless double beta decay matrix elements in light nuclei

abstract: We present the first ab initio calculations of neutrinoless double beta decay matrix elements in $A=6$-$12$ nuclei using Variational Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon potential and Illinois-7 three-nucleon interaction. We study both light Majorana neutrino exchange and potentials arising from a large class of multi-TeV mechanisms of lepton number violation. Our results provide benchmarks to be used in testing many-body methods that can be extended to the heavy nuclei of experimental interest. In light nuclei we have also studied the impact of two-body short range correlations and the use of different forms for the transition operators, such as those corresponding to different orders in chiral effective theory.


title: Real time ridge orientation estimation for fingerprint images

abstract: Fingerprint verification is an important bio-metric technique for personal identification. Most of the automatic verification systems are based on matching of fingerprint minutiae. Extraction of minutiae is an essential process which requires estimation of orientation of the lines in an image. Most of the existing methods involve intense mathematical computations and hence are performed through software means. In this paper a hardware scheme to perform real time orientation estimation is presented which is based on pipelined architecture. Synthesized circuits proved the functionality and accuracy of the suggested method.


title: Dust-trapping vortices and a potentially planet-triggered spiral wake in the pre-transitional disk of V1247 Orionis

abstract: The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.04"-resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply-confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.17"=54 au separation from the star) and the crescent (at 0.38"=120 au) seem smoothly connected through a one-armed spiral arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at $\sim0.3$"$\approx$100 au, where the one-armed spiral arm detected in polarised light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap mm-sized particles, namely the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered light images. At the position of the spiral wake and the crescent we also observe $^{12}$CO (3-2) and H$^{12}$CO$^{+}$ (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.


title: Modelling the line-of-sight contribution in substructure lensing

abstract: We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.


title: Evidence for a Sub-Chandrasekhar Mass Type~Ia Supernova in the Ursa Minor Dwarf Galaxy

abstract: A longstanding problem is identifying the elusive progenitors of Type Ia supernovae (SNe Ia), which can roughly be split into Chandraksekhar and sub-Chandrasekhar mass events. An important difference between these two cases is the nucleosynthetic yield, which is altered by the increased neutron excess in Chandrasekhar progenitors due to their pre-explosion simmering and high central density. From comparison with theoretical nucleosynthesis yields, we show that the chemical composition of the most metal-rich star in the Ursa Minor dwarf galaxy, COS 171, is dominated by nucleosynthesis from a low-metallicity, low-mass, sub-Chandrasekhar mass SN Ia. Key diagnostic abundance ratios include C/Fe, {\alpha}/Fe, Mn/Fe and Ni/Fe ratios, which could not have been produced by Chandrasekhar-mass SNe Ia, Core-Collapse Type II supernovae or Pair-Instsability supernovae. Strong deficiencies of Ni/Fe, Cu/Fe and Zn/Fe also suggest the absence of alpha-rich freeze-out nucleosynthesis, favoring low-mass WD progenitor SNe Ia. Based on comparisons of the measured Mn/Fe and Si/Fe ratios with detonation models, we estimate a WD mass near 0.95Msun. We also compare Mn/Fe and Ni/Fe ratios to the recent theoretical yields predicted by Shen et al., finding consistent results. To explain the COS 171 [Fe/H], at -1.35 dex, requires dilution of the ejecta from a single SNIa event with ~10^4 Msun of material; this is expected for a SN Ia remnant expanding into a warm interstellar medium with n~1 /cm^3. In the future, finding more stars with the unique chemical signatures we highlight here will be important for constraining the rate and environments of sub-Chandrasekhar SNe Ia.


title: Hard X-ray Emission from the M87 AGN Detected with NuSTAR

abstract: M87 hosts a 3-6 billion solar mass black hole with a remarkable relativistic jet that has been regularly monitored in radio to TeV bands. However, hard X-ray emission \gtrsim 10keV, which would be expected to primarily come from the jet or the accretion flow, had never been detected from its unresolved X-ray core. We report NuSTAR detection up to 40 keV from the the central regions of M87. Together with simultaneous Chandra observations, we have constrained the dominant hard X-ray emission to be from its unresolved X-ray core, presumably in its quiescent state. The core spectrum is well fitted by a power law with photon index Gamma=2.11 (+0.15 -0.11). The measured flux density at 40 keV is consistent with a jet origin, although emission from the advection-dominated accretion flow cannot be completely ruled out. The detected hard X-ray emission is significantly lower than that predicted by synchrotron self-Compton models introduced to explain emission above a GeV.


title: Primordial Black Holes as Dark Matter: Constraints From Compact Ultra-Faint Dwarfs

abstract: The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.


title: A non-probabilistic substitute for the Born rule

abstract: The Born rule assigns a probability to any possible outcome of a quantum measurement, but leaves open the question how these probabilities are to be interpreted and, in particular, how they relate to the outcome observed in an actual experiment. We propose to avoid this question by replacing the Born rule with two non-probabilistic postulates: (i) the projector associated to the observed outcome must have a positive overlap with the state of the measured system; (ii) statements about observed outcomes are robust, that is, remain valid under small perturbations of the state. We show that the two postulates suffice to retrieve the interpretations of the Born rule that are commonly used for analysing experimental data.


title: SDSS IV MaNGA - sSFR profiles and the slow quenching of discs in green valley galaxies

abstract: We study radial profiles in H$\alpha$ equivalent width and specific star formation rate (sSFR) derived from spatially-resolved SDSS-IV MaNGA spectroscopy to gain insight on the physical mechanisms that suppress star formation and determine a galaxy's location in the SFR-$\rm M_\star$ diagram. Even within the star-forming `main sequence', the measured sSFR decreases with stellar mass, both in an integrated and spatially-resolved sense. Flat sSFR radial profiles are observed for $\rm log(M_\star/ M_\odot) < 10.5$, while star-forming galaxies of higher mass show a significant decrease in sSFR in the central regions, a likely consequence of both larger bulges and an inside-out growth history. Our primary focus is the green valley, constituted by galaxies lying below the star formation main sequence, but not fully passive. In the green valley we find sSFR profiles that are suppressed with respect to star-forming galaxies of the same mass at all galactocentric distances out to 2 effective radii. The responsible quenching mechanism therefore appears to affect the entire galaxy, not simply an expanding central region. The majority of green valley galaxies of $\rm log(M_\star/ M_\odot) > 10.0$ are classified spectroscopically as central low-ionisation emission-line regions (cLIERs). Despite displaying a higher central stellar mass concentration, the sSFR suppression observed in cLIER galaxies is not simply due to the larger mass of the bulge. Drawing a comparison sample of star forming galaxies with the same $\rm M_\star$ and $\rm \Sigma_{1~kpc}$ (the mass surface density within 1 kpc), we show that a high $\rm \Sigma_{1~kpc}$ is not a sufficient condition for determining central quiescence.


title: Antiferromagnetic Chern insulators in non-centrosymmetric systems

abstract: We investigate a new class of topological antiferromagnetic (AF) Chern insulators driven by electronic interactions in two-dimensional systems without inversion symmetry. Despite the absence of a net magnetization, AF Chern insulators (AFCI) possess a nonzero Chern number $C$ and exhibit the quantum anomalous Hall effect (QAHE). Their existence is guaranteed by the bifurcation of the boundary line of Weyl points between a quantum spin Hall insulator and a topologically trivial phase with the emergence of AF long-range order. As a concrete example, we study the phase structure of the honeycomb lattice Kane-Mele model as a function of the inversion-breaking ionic potential and the Hubbard interaction. We find an easy $z$-axis $C=1$ AFCI phase and a spin-flop transition to a topologically trivial $xy$-plane collinear antiferromagnet. We propose experimental realizations of the AFCI and QAHE in correlated electron materials and cold atom systems.


title: Exploring one particle orbitals in large Many-Body Localized systems

abstract: Strong disorder in interacting quantum systems can give rise to the phenomenon of Many-Body Localization (MBL), which defies thermalization due to the formation of an extensive number of quasi local integrals of motion. The one particle operator content of these integrals of motion is related to the one particle orbitals of the one particle density matrix and shows a strong signature across the MBL transition as recently pointed out by Bera et al. [Phys. Rev. Lett. 115, 046603 (2015); Ann. Phys. 529, 1600356 (2017)]. We study the properties of the one particle orbitals of many-body eigenstates of an MBL system in one dimension. Using shift-and-invert MPS (SIMPS), a matrix product state method to target highly excited many-body eigenstates introduced in [Phys. Rev. Lett. 118, 017201 (2017)], we are able to obtain accurate results for large systems of sizes up to L = 64. We find that the one particle orbitals drawn from eigenstates at different energy densities have high overlap and their occupations are correlated with the energy of the eigenstates. Moreover, the standard deviation of the inverse participation ratio of these orbitals is maximal at the nose of the mobility edge. Also, the one particle orbitals decay exponentially in real space, with a correlation length that increases at low disorder. In addition, we find a 1/f distribution of the coupling constants of a certain range of the number operators of the OPOs, which is related to their exponential decay.


title: Dephasing catastrophe in $4 - \epsilon$ dimensions: A possible instability of the ergodic (many-body-delocalized) phase

abstract: In two dimensions (2D), dephasing by a bath cuts off Anderson localization that would otherwise occur at any energy density for fermions with disorder. For an isolated system with short-range interactions, the system can be its own bath, exhibiting diffusive (non-Markovian) thermal density fluctuations. We recast the dephasing of weak localization due to a diffusive bath as a self-interacting polymer loop. We investigate the critical behavior of the loop in $d=4-\epsilon$ dimensions, and find a nontrivial fixed point corresponding to a temperature $T^* \sim \epsilon >0$ where the dephasing time diverges. Assuming that this fixed point survives to $\epsilon=2$, we associate it to a possible instability of the ergodic phase. Our approach may open a new line of attack against the problem of the ergodic to many-body-localized phase transition in $d > 1$ spatial dimensions.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el


title: Time Evolution of Gamma Rays from Supernova Remnants

abstract: We present a systematic phenomenological study focused on the time evolution of the non-thermal radiation - from radio waves to gamma rays - emitted by typical supernova remnants via hadronic and leptonic mechanisms, for two classes of progenitors: thermonuclear and core-collapse. To this aim, we develop a numerical tool designed to model the evolution of the cosmic-ray spectrum inside a supernova remnant, and compute the associated multi-wavelength emission. We demonstrate the potential of this tool in the context of future population studies based on large collection of high-energy gamma-ray data. We discuss and explore the relevant parameter space involved in the problem, and focus in particular on their impact on the maximum energy of accelerated particles, in order to study the effectiveness and duration of the PeVatron phase. We outline the crucial role of the ambient medium through which the shock propagates during the remnant evolution. In particular, we point out the role of dense clumps in creating a significant hardening in the hadronic gamma-ray spectrum.


title: Virtual homological spectral radii for automorphisms of surfaces

abstract: In this paper, it is shown that any surface automorphism of positive mapping-class entropy possesses a virtual homological eigenvalue which lies outside the unit circle of the complex plane.


title: On nodal solutions of a nonlocal Choquard equation in a bounded domain

abstract: In this paper, we are interested in the least energy nodal solutions to the following nonlocal Choquard equation with a local term \begin{equation*}\left\{\begin{array}{rll} -\Delta u&=\lambda|u|^{p-2}u+\mu \phi(x)|u|^{q-2}u\\ -\Delta \phi&=|u|^q\\ u&=\phi=0 \end{array}\right. \begin{gathered}\begin{array}{rll} &\mbox{in}\ \Omega,\\ &\mbox{in}\ \Omega,\\ &\mbox{on}\ \partial\Omega, \end{array}\end{gathered}\end{equation*} where $ \lambda,\mu>0, p\in [2,6), q\in (1,5)$ and $\Omega\subset \mathbb{R}^3$ is a bounded domain. This problem may be seen as a nonlocal perturbation of the classical Lane-Emden equation $-\Delta u=\lambda|u|^{p-2}u$ in $\Omega.$ The problem has a variational functional with a nonlocal term $\mu\int_{\Omega}\phi|u|^q$. The appearance of the nonlocal term makes the variational functional very different from the local case $\mu=0$, for which the problem has ground state solutions and least energy nodal solutions if $p\in (2,6)$. The problem may also be viewed as a nonlocal Choquard equation with a local pertubation term when $\lambda \not =0$. For $\mu>0$, we show that although ground state solutions always exist, the existence of least energy nodal solution depends on $q$: for $q\in(1,2)$ there does not exist a least energy nodal solution while for $q\in[2,5)$ such a solution exists. Note that $q=2$ is a critical value. In the case of a linear local perturbation, i.e., $p=2,$ if $\lambda<\lambda_1,$ the problem has a positive ground state and a least energy nodal solution. However, if $\lambda\geq \lambda_1,$ the problem has a ground state which changes sign. Hence it is also a least energy nodal solution.


title: The Weizs\"acker-Williams distribution of linearly polarized gluons (and its fluctuations) at small x

abstract: The conventional and linearly polarized Weizs\"acker-Williams gluon distributions at small x are defined from the two-point function of the gluon field in light-cone gauge. They appear in the cross section for dijet production in deep inelastic scattering at high energy. We determine these functions in the small-x limit from solutions of the JIMWLK evolution equations and show that they exhibit approximate geometric scaling. Also, we discuss the functional distributions of these WW gluon distributions over the JIMWLK ensemble at rapidity $Y\sim 1/\alpha_s$. These are determined by a 2d Liouville action for the logarithm of the covariant gauge function $g^2 \mathrm{tr}\,A^+(q)A^+(-q)$. For transverse momenta on the order of the saturation scale we observe large variations across configurations (evolution trajectories) of the linearly polarized distribution up to several times its average, and even to negative values.


title: X-ray mapping of the stellar wind in the binary PSR J2032+4127/MT91 213

abstract: PSR J2032+4127 is a young and rapidly rotating pulsar on a highly eccentric orbit around the high-mass Be star MT91 213. X-ray monitoring of the binary system over a $\sim$4000 day period with Swift has revealed an increase of the X-ray luminosity which we attribute to the synchrotron emission of the shocked pulsar wind. We use Swift X-ray observations to infer a clumpy stellar wind with $r^{-2}$ density profile and constrain the Lorentz factor of the pulsar wind to $10^5<\gamma_w<10^6$. We investigate the effects of an axisymmetric stellar wind with polar gradient on the X-ray emission. Comparison of the X-ray light curve hundreds of days before and after the periastron can be used to explore the polar structure of the wind.


title: Further study on periodic solutions of elliptic equations with a fractional Laplacian

abstract: We obtain some existence theorems for periodic solutions to several linear equations involving fractional Laplacian. We also prove that the lower bound of all periods for semilinear elliptic equations involving fractional Laplacian is not larger than some exact positive constant. Hamiltonian identity, Modica-type inequalities and an estimate of the energy for periodic solutions are also established.


title: ThermSense: Smartphone-based Breathing Sensing Platform using Noncontact Low-Cost Thermal Camera

abstract: The ability of sensing breathing is becoming an increasingly important function for technology that aims at supporting both psychological and physical wellbeing. We demonstrate ThermSense, a new breathing sensing platform based on smartphone technology and low-cost thermal camera, which allows a user to measure his/her breathing pattern in a contact-free manner. With the designed key functions of Thermal Voxel Integration-based breathing estimation and respiration variability spectrogram (RVS, bi-dimensional representation of breathing dynamics), the developed platform provides scalability and flexibility for gathering respiratory physiological measurements ubiquitously. The functionality could be used for a variety of applications from stress monitoring to respiration training.


title: Density split statistics: Cosmological constraints from counts and lensing in cells in DES Y1 and SDSS data

abstract: We derive cosmological constraints from the probability distribution function (PDF) of evolved large-scale matter density fluctuations. We do this by splitting lines of sight by density based on their count of tracer galaxies, and by measuring both gravitational shear around and counts-in-cells in overdense and underdense lines of sight, in Dark Energy Survey (DES) First Year and Sloan Digital Sky Survey (SDSS) data. Our analysis uses a perturbation theory model (see companion paper Friedrich at al.) and is validated using N-body simulation realizations and log-normal mocks. It allows us to constrain cosmology, bias and stochasticity of galaxies w.r.t. matter density and, in addition, the skewness of the matter density field. From a Bayesian model comparison, we find that the data weakly prefer a connection of galaxies and matter that is stochastic beyond Poisson fluctuations on <=20 arcmin angular smoothing scale. The two stochasticity models we fit yield DES constraints on the matter density $\Omega_m=0.26^{+0.04}_{-0.03}$ and $\Omega_m=0.28^{+0.05}_{-0.04}$ that are consistent with each other. These values also agree with the DES analysis of galaxy and shear two-point functions (3x2pt) that only uses second moments of the PDF. Constraints on $\sigma_8$ are model dependent ($\sigma_8=0.97^{+0.07}_{-0.06}$ and $0.80^{+0.06}_{-0.07}$ for the two stochasticity models), but consistent with each other and with the 3x2pt results if stochasticity is at the low end of the posterior range. As an additional test of gravity, counts and lensing in cells allow to compare the skewness $S_3$ of the matter density PDF to its LCDM prediction. We find no evidence of excess skewness in any model or data set, with better than 25 per cent relative precision in the skewness estimate from DES alone.


title: Observational constraints on warm quasi-exponential inflation

abstract: In the present work we study a warm inflationary model defined by a quasi-exponential inflaton potential and an inflaton decay rate proportional to the Hubble rate. The model is characterized by three free parameters. We compute the power spectrum, the scalar spectral index as well as the tensor-to-scalar ratio within the framework of the model, and we compare with the latest Planck data. On the $r-n_s$ plane we show both the theoretical curves and the contour plots allowed by observations, and we constrain the parameters of the model accordingly. The non-linear parameter $f_{NL}$, corresponding to primordial non-Gaussianities, is also discussed and we found that the value predicted by our model is within the bounds imposed by current observational data.


title: A Toolkit For Steady States of Nonlinear Wave Equations: Continuous Time Nesterov and Exponential Time Differencing Schemes

abstract: Several methods exist for finding ground (as well as excited) states of nonlinear waves equations. In this paper we first introduce two modifications of the so-called accelerated imaginary-time evolution method (AITEM). In our first modification, time integration of the underlying gradient flow is done using exponential time differencing instead of using more standard methods. In the second modification, we present a generalization of the gradient flow model, motivated by the work of Nesterov, as well as that of Candes and collaborators. Additionally, we consider combinations of these methods with the so-called spectral renormalization scheme. Finally, we apply these techniques to the so-called Squared Operator Method, enabling convergence to excited states. Various examples are shown to illustrate the effectiveness of these new schemes, comparing them to standard ones established in the literature. In most cases, we find significant reductions in the number of iterations needed to reach convergence.


title: Long-Term Inertial Navigation Aided by Dynamics of Flow Field Features

abstract: A current-aided inertial navigation framework is proposed for small autonomous underwater vehicles in long-duration operations (> 1 hour), where neither frequent surfacing nor consistent bottom-tracking are available. We instantiate this concept through mid-depth, underwater navigation. This strategy mitigates dead-reckoning uncertainty of a traditional inertial navigation system by comparing the estimate of local, ambient flow velocity with preloaded ocean current maps. The proposed navigation system is implemented through a marginalized particle filter where the vehicle's states are sequentially tracked along with sensor bias and local turbulence that is not resolved by general flow prediction. The performance of the proposed approach is first analyzed through Monte Carlo simulations in two artificial background flow fields, resembling real-world ocean circulation patterns, superposed with smaller-scale, turbulent components with Kolmogorov energy spectrum. The current-aided navigation scheme significantly improves the dead-reckoning performance of the vehicle even when unresolved, small-scale flow perturbations are present. For a 6-hour navigation with an automotive-grade inertial navigation system, the current-aided navigation scheme results in positioning estimates with under 3% uncertainty per distance traveled (UDT) in a turbulent, double-gyre flow field, and under 7.3% UDT in a turbulent, meandering jet flow field. Further evaluation with field test data and actual ocean simulation analysis demonstrates consistent performance for a 6-hour mission, positioning result with under 25% UDT for a 24-hour navigation when provided direct heading measurements, and terminal positioning estimate with 16% UDT at the cost of increased uncertainty at an early stage of the navigation.


title: SDSS-IV MaNGA: The Spatial Distribution of Star Formation and its Dependence on Mass, Structure and Environment

abstract: We study the spatially resolved star formation of 1494 galaxies in the SDSSIV-MaNGA Survey. SFRs are calculated using a two-step process, using $H_\alpha$ in star forming regions and $D_n4000$ in regions identified as AGN/LI(N)ER or lineless. The roles of secular and environmental quenching processes are investigated by studying the dependence of the radial profiles of specific star formation rate on stellar mass, galaxy structure and environment. We report on the existence of `Centrally Suppressed' galaxies, which have suppressed SSFR in their cores compared to their disks. The profiles of centrally suppressed and unsuppressed galaxies are distibuted in a bimodal way. Galaxies with high stellar mass and core velocity dispersion are found to be much more likely to be centrally suppressed than low mass galaxies, and we show that this is related to morphology and the presence of AGN/LI(N)ER like emission. Centrally suppressed galaxies also display lower star formation at all radii compared to unsuppressed galaxies. The profiles of central and satellite galaxies are also compared, and we find that satellite galaxies experience lower specific star formation rates at all radii than central galaxies. This uniform suppression could be a signal of the stripping of hot halo gas in the process known as strangulation. We find that satellites are not more likely to be suppressed in their cores than centrals, indicating that the core suppression is an entirely internal process. We find no correlation between the local environment density and the profiles of star formation rate surface density.


title: Learning Independent Features with Adversarial Nets for Non-linear ICA

abstract: Reliable measures of statistical dependence could be useful tools for learning independent features and performing tasks like source separation using Independent Component Analysis (ICA). Unfortunately, many of such measures, like the mutual information, are hard to estimate and optimize directly. We propose to learn independent features with adversarial objectives which optimize such measures implicitly. These objectives compare samples from the joint distribution and the product of the marginals without the need to compute any probability densities. We also propose two methods for obtaining samples from the product of the marginals using either a simple resampling trick or a separate parametric distribution. Our experiments show that this strategy can easily be applied to different types of model architectures and solve both linear and non-linear ICA problems.


title: Device-independent quantum private comparison protocol without a third party

abstract: Since unconditionally secure quantum two-party computations are known to be impossible, most existing quantum private comparison (QPC) protocols adopted a third party. Recently, we proposed a QPC protocol which involves two parties only, and showed that although it is not unconditionally secure, it only leaks an extremely small amount of information to the other party. Here we further propose the device-independent version of the protocol, so that it can be more convenient and dependable in practical applications.


title: Non-perturbative decay of Non-Abelian hair

abstract: We construct a solution of Heterotic supergravity which interpolates between two different AdS$_{3}\times S^{3}\times T^{4}$ geometries corresponding to the near-horizon limits of two 5-dimensional black holes, only one of which has non-Abelian hair. This solution can be used to estimate the amplitude of probability of the non-perturbative decay of the gauge 5-brane responsible for the non-Abelian hair into eight solitonic 5-branes by evaluating its Euclidean action. The Wick rotation of this solution poses several problems which we argue can be overcome by using a non-extremal off-shell (NEOS) deformation of the solution. This NEOS field configuration can be Wick rotated straight away and its Euclidean action can be computed for any value of the deformation parameter. The Euclidean result can then be anti-Wick-rotated and its extremal limit gives the Euclidean action of the original solution, which turns out to be one half of the difference between the entropies of the 5-dimensional black holes.


title: Automated Scalable Bayesian Inference via Hilbert Coresets

abstract: The automation of posterior inference in Bayesian data analysis has enabled experts and nonexperts alike to use more sophisticated models, engage in faster exploratory modeling and analysis, and ensure experimental reproducibility. However, standard automated posterior inference algorithms are not tractable at the scale of massive modern datasets, and modifications to make them so are typically model-specific, require expert tuning, and can break theoretical guarantees on inferential quality. Building on the Bayesian coresets framework, this work instead takes advantage of data redundancy to shrink the dataset itself as a preprocessing step, providing fully-automated, scalable Bayesian inference with theoretical guarantees. We begin with an intuitive reformulation of Bayesian coreset construction as sparse vector sum approximation, and demonstrate that its automation and performance-based shortcomings arise from the use of the supremum norm. To address these shortcomings we develop Hilbert coresets, i.e., Bayesian coresets constructed under a norm induced by an inner-product on the log-likelihood function space. We propose two Hilbert coreset construction algorithms---one based on importance sampling, and one based on the Frank-Wolfe algorithm---along with theoretical guarantees on approximation quality as a function of coreset size. Since the exact computation of the proposed inner-products is model-specific, we automate the construction with a random finite-dimensional projection of the log-likelihood functions. The resulting automated coreset construction algorithm is simple to implement, and experiments on a variety of models with real and synthetic datasets show that it provides high-quality posterior approximations and a significant reduction in the computational cost of inference.


title: Hydrodynamic-to-ballistic crossover in Dirac fluid

abstract: We develop an exactly solvable classical kinetic model of transport in Dirac materials accounting for strong electron-electron (e-e) and electron-hole (e-h) collisions. We use this model to track the evolution of graphene conductivity and properties of its collective excitations across the hydrodynamic-to-ballistic crossover. We find the relaxation rate of electric current by e-e collisions that is possible due to the lack of Galilean invariance, and introduce a universal numerical measure of this non-invariance in arbitrary dimension. We find the two branches of collective excitations in the Dirac fluid: plasmons and electron-hole sound. The sound waves have small viscous damping at the neutrality point both in the hydrodynamic and ballistic regimes, but acquire large damping due to e-h friction even at slight doping. On the contrary, plasmons acquire strong frictional damping at the neutrality point and become well-defined in doped samples.


title: Function spaces obeying a time-varying bandlimit

abstract: Motivated by applications to signal processing and mathematical physics, recent work on the concept of time-varying bandwidth has produced a class of function spaces which generalize the Paley-Wiener spaces of bandlimited functions: any regular simple symmetric linear transformation with deficiency indices $(1,1)$ is naturally represented as multiplication by the independent variable in one of these spaces. We explicitly demonstrate the equivalence of this model for such linear transformations to several other functional models based on the theories of meromorphic model spaces of Hardy space and purely atomic Herglotz measures on the real line, respectively. This theory provides a precise notion of a time-varying or local bandwidth, and we describe how it may be applied to construct signal processing techniques that are adapted to signals obeying a time-varying bandlimit.


title: Compressive three-dimensional super-resolution microscopy with speckle-saturated fluorescence excitation

abstract: Nonlinear structured illumination microscopy (nSIM) is an effective approach for super-resolution wide-field fluorescence microscopy with a theoretically unlimited resolution. In nSIM, carefully designed, highly-contrasted illumination patterns are combined with the saturation of an optical transition to enable sub-diffraction imaging. While the technique proved useful for two-dimensional imaging, extending it to three-dimensions (3D) is challenging due to the fading/fatigue of organic fluorophores under intense cycling conditions. Here, we present a compressed sensing approach that allows for the first time 3D sub-diffraction nSIM of cultured cells by saturating fluorescence excitation. Exploiting the natural orthogonality of transverse speckle illumination planes, 3D probing of the sample is achieved by a single two-dimensional scan. Fluorescence contrast under saturated excitation is ensured by the inherent high density of intensity minima associated with optical vortices in polarized speckle patterns. Compressed speckle microscopy is thus a simple approach that enables 3D super-resolved nSIM imaging with potentially considerably reduced acquisition time and photobleaching.les fast 3D super-resolved imaging with considerably minimized photo-bleaching.


title: A simple discharging method for forbidden subposet problems

abstract: The poset $Y_{k+1, 2}$ consists of $k+2$ distinct elements $x_1$, $x_2$, \dots, $x_{k}$, $y_1$,$y_2$, such that $x_1 \le x_2 \le \dots \le x_{k} \le y_1$,~$y_2$. The poset $Y'_{k+1, 2}$ is the dual of $Y_{k+1, 2}$ Let $\rm{La}^{\sharp}(n,\{Y_{k+1, 2}, Y'_{k+1, 2}\})$ be the size of the largest family $\mathcal{F} \subset 2^{[n]}$ that contains neither $Y_{k+1,2}$ nor $Y'_{k+1,2}$ as an induced subposet. Methuku and Tompkins proved that $\rm{La}^{\sharp}(n, \{Y_{3,2}, Y'_{3,2}\}) = \Sigma(n,2)$ for $n \ge 3$ and they conjectured the generalization that if $k \ge 2$ is an integer and $n \ge k+1$, then $\rm{La}^{\sharp}(n, \{Y_{k+1,2}, Y'_{k+1,2}\}) = \Sigma(n,k)$. In this paper, we introduce a simple discharging approach and prove this conjecture.


title: A note on primes in certain residue classes

abstract: Given positive integers $a_1,\ldots,a_k$, we prove that the set of primes $p$ such that $p \not\equiv 1 \bmod{a_i}$ for $i=1,\ldots,k$ admits asymptotic density relative to the set of all primes which is at least $\prod_{i=1}^k \left(1-\frac{1}{\varphi(a_i)}\right)$, where $\varphi$ is the Euler's totient function. This result is similar to the one of Heilbronn and Rohrbach, which says that the set of positive integer $n$ such that $n \not\equiv 0 \bmod a_i$ for $i=1,\ldots,k$ admits asymptotic density which is at least $\prod_{i=1}^k \left(1-\frac{1}{a_i}\right)$.


title: On weighted approximation with Jacobi weights

abstract: We obtain matching direct and inverse theorems for the degree of weighted $L_p$-approximation by polynomials with the Jacobi weights $(1-x)^\alpha (1+x)^\beta$. Combined, the estimates yield a constructive characterization of various smoothness classes of functions via the degree of their approximation by algebraic polynomials.


title: Functional Decision Theory: A New Theory of Instrumental Rationality

abstract: This paper describes and motivates a new decision theory known as functional decision theory (FDT), as distinct from causal decision theory and evidential decision theory. Functional decision theorists hold that the normative principle for action is to treat one's decision as the output of a fixed mathematical function that answers the question, "Which output of this very function would yield the best outcome?" Adhering to this principle delivers a number of benefits, including the ability to maximize wealth in an array of traditional decision-theoretic and game-theoretic problems where CDT and EDT perform poorly. Using one simple and coherent decision rule, functional decision theorists (for example) achieve more utility than CDT on Newcomb's problem, more utility than EDT on the smoking lesion problem, and more utility than both in Parfit's hitchhiker problem. In this paper, we define FDT, explore its prescriptions in a number of different decision problems, compare it to CDT and EDT, and give philosophical justifications for FDT as a normative theory of decision-making.


title: A New Technique for Reachability of States in Concatenation Automata

abstract: We present a new technique for demonstrating the reachability of states in deterministic finite automata representing the concatenation of two languages. Such demonstrations are a necessary step in establishing the state complexity of the concatenation of two languages, and thus in establishing the state complexity of concatenation as an operation. Typically, ad-hoc induction arguments are used to show particular states are reachable in concatenation automata. We prove some results that seem to capture the essence of many of these induction arguments. Using these results, reachability proofs in concatenation automata can often be done more simply and without using induction directly.


title: The MUSE Hubble Ultra Deep Field Survey: III. Testing photometric redshifts to 30th magnitude

abstract: We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at z>3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.


title: A Distributed Auction Policy for User Association in Device-to-Device Caching Networks

abstract: We propose a distributed bidding-aided Matern carrier sense multiple access (CSMA) policy for device-to-device (D2D) content distribution. The network is composed of D2D receivers and potential D2D transmitters, i.e., transmitters are turned on or off by the scheduling algorithm. Each D2D receiver determines the value of its request, by bidding on the set of potential transmitters in its communication range. Given a medium access probability, a fraction of the potential transmitters are jointly scheduled, i.e., turned on, determined jointly by the auction policy and the power control scheme. The bidding-aided scheduling algorithm exploits (i) the local demand distribution, (ii) spatial distribution of D2D node locations, and (iii) the cache configurations of the potential transmitters. We contrast the performance of the bidding-aided CSMA policy with other well-known CSMA schemes that do not take into account (i)-(iii), demonstrate that our algorithm achieves a higher spectral efficiency in terms of the number of bits transmitted per unit time per unit bandwidth per user. The gain becomes even more visible under randomized configurations and requests rather than more skewed placement configurations and deterministic demand distributions.


title: Unusual negative formation enthalpies and atomic ordering in isovalent alloys of transition metal dichalcogenide monolayers

abstract: Common substitutional isovalent semiconductor alloys usually form disordered metastable phases with positive excess formation enthalpies ({\Delta}H). In contrast, monolayer alloys of transition metal dichalcogenides (TMDs) MX2 (M = Mo, W; X = S, Se) always have negative {\Delta}H, suggesting atomic ordering, which is, however, not yet experimentally observed. Using first-principles calculations, we find that the negative {\Delta}H of cation-mixed TMD alloys results from the charge transfer from weak Mo-X to nearest strong W-X bonds and the negative {\Delta}H of anion-mixed TMD alloys comes from the larger energy gain due to the charge transfer from Se to nearest S atoms than the energy cost due to the lattice mismatch. Consequently, cation-mixed and anion-mixed alloys should energetically prefer to have Mo-X-W and S-M-Se ordering, respectively. The atomic ordering, however, is only locally ordered but disordered in the long range due to the symmetry of TMD monolayers, as demonstrated by many energetically degenerate structures for given alloy compositions. Besides, the local ordering and disordering effects on the macroscopic properties such as bandgaps and optical absorptions are negligible, making the experimental observation of locally ordered TMD alloys challenging. We propose to take the advantage of microscopic properties such as defects which strongly depend on local atomic configurations for experiments to identify the disordering and local ordering in TMD alloys. Finally, quaternary TMD alloys by mixing both cations and anions are studied to have a wide range of bandgaps for optoelectronic applications. Our work is expected to help the formation and utilization of TMD alloys.


title: Portable Devices for In-Vitro Characterization Based on Ultrasound

abstract: We propose a portable cartridge based lab-on-a-chip platform. The main sensing mechanism uses high frequency acoustic waves. Feasibility of two different transducers has been investigated for biological measurements. We used zincoxide (ZnO) based piezoelectric transducers for blood coagulation time measurements. The method is based on the measurement of amplitude of acoustic reflections from a microfluidic channel which is filled with blood. The volume of the microfluidic channel could be as low as 1 micro-litre. A low cost disposable cartridge made of glass has been designed to perform the measurements. The cartridge is composed of two glass substrates and a double sided tape layer where the channel is defined. One of the glass substrates has a micromachined ZnO based transducer that operates at 400 MHz on the outer surface of the cartridge where electrical connections are also provided. The transducer is aligned with the channel to generate acoustic pulses in the fluid that is filling the channel. The reflections coming from the top of the channel propagate through the liquid therefore their amplitude and phase are affected by the fluid properties. In the experiments whole blood was used without any sample preparation. The viscosity of the blood changes during coagulation therefore by monitoring the amplitude of the reflection coming from the top of the channel one can measure the coagulation time. The method which requires only 1 micro-litre of blood has been tested with calcium and activated partial thromboplastin (aPTT) reagents. The proposed method has a potential to be used in a low cost portable coagulation time measurement system for patient self-testing. The second transducer was made of CMUTs (Capacitive Micromachined Ultrasonic Transducers). We introduce this technology for immersion in liquids as a biological sensor for viscosity measurements.


title: Spatial patterns and biodiversity in off-lattice simulations of a cyclic three-species Lotka-Volterra model

abstract: Stochastic simulations of cyclic three-species spatial predator-prey models are usually performed in square lattices with nearest neighbor interactions starting from random initial conditions. In this Letter we describe the results of off-lattice Lotka-Volterra stochastic simulations, showing that the emergence of spiral patterns does occur for sufficiently high values of the (conserved) total density of individuals. We also investigate the dynamics in our simulations, finding an empirical relation characterizing the dependence of the characteristic peak frequency and amplitude on the total density. Finally, we study the impact of the total density on the extinction probability, showing how a low population density may jeopardize biodiversity.


title: Blindfold learning of an accurate neural metric

abstract: The brain has no direct access to physical stimuli, but only to the spiking activity evoked in sensory organs. It is unclear how the brain can structure its representation of the world based on differences between those noisy, correlated responses alone. Here we show how to build a distance map of responses from the structure of the population activity of retinal ganglion cells, allowing for the accurate discrimination of distinct visual stimuli from the retinal response. We introduce the Temporal Restricted Boltzmann Machine to learn the spatiotemporal structure of the population activity, and use this model to define a distance between spike trains. We show that this metric outperforms existing neural distances at discriminating pairs of stimuli that are barely distinguishable. The proposed method provides a generic and biologically plausible way to learn to associate similar stimuli based on their spiking responses, without any other knowledge of these stimuli.


title: Universal behavior of modulationally unstable media

abstract: Evidence is presented of universal behavior in modulationally unstable media. An ensemble of nonlinear evolution equations, including three partial differential equations, an integro-differential equation, a nonlocal system and a differential-difference equation, is studied. Collectively, these systems arise in a variety of applications in the physical and mathematical sciences, including water waves, optics, acoustics, Bose-Einstein condensation, and more. All these models exhibit modulational instability, namely, the property that a constant background is unstable to long-wavelength perturbations. In this work, each of these systems is studied analytically and numerically for a number of different initial perturbations of the constant background, and it is shown that, for all systems and for all initial conditions considered, the dynamics gives rise to a remarkably similar structure comprised of two outer, quiescent sectors separated by a wedge-shaped central region characterized by modulated periodic oscillations. A heuristic criterion that allows one to compute some of the properties of the central oscillation region is also given.


title: Kumaraswamy autoregressive moving average models for double bounded environmental data

abstract: In this paper we introduce the Kumaraswamy autoregressive moving average models (KARMA), which is a dynamic class of models for time series taking values in the double bounded interval $(a,b)$ following the Kumaraswamy distribution. The Kumaraswamy family of distribution is widely applied in many areas, especially hydrology and related fields. Classical examples are time series representing rates and proportions observed over time. In the proposed KARMA model, the median is modeled by a dynamic structure containing autoregressive and moving average terms, time-varying regressors, unknown parameters and a link function. We introduce the new class of models and discuss conditional maximum likelihood estimation, hypothesis testing inference, diagnostic analysis and forecasting. In particular, we provide closed-form expressions for the conditional score vector and conditional Fisher information matrix. An application to environmental real data is presented and discussed.