title: Modeling National Latent Socioeconomic Health and Examination of Policy Effects via Causal Inference

abstract: This research develops a socioeconomic health index for nations through a model-based approach which incorporates spatial dependence and examines the impact of a policy through a causal modeling framework. As the gross domestic product (GDP) has been regarded as a dated measure and tool for benchmarking a nation's economic performance, there has been a growing consensus for an alternative measure---such as a composite `wellbeing' index---to holistically capture a country's socioeconomic health performance. Many conventional ways of constructing wellbeing/health indices involve combining different observable metrics, such as life expectancy and education level, to form an index. However, health is inherently latent with metrics actually being observable indicators of health. In contrast to the GDP or other conventional health indices, our approach provides a holistic quantification of the overall `health' of a nation. We build upon the latent health factor index (LHFI) approach that has been used to assess the unobservable ecological/ecosystem health. This framework integratively models the relationship between metrics, the latent health, and the covariates that drive the notion of health. In this paper, the LHFI structure is integrated with spatial modeling and statistical causal modeling, so as to evaluate the impact of a policy variable (mandatory maternity leave days) on a nation's socioeconomic health, while formally accounting for spatial dependency among the nations. We apply our model to countries around the world using data on various metrics and potential covariates pertaining to different aspects of societal health. The approach is structured in a Bayesian hierarchical framework and results are obtained by Markov chain Monte Carlo techniques.


title: An analysis of the induced linear operators associated to divide and color models

abstract: We study the natural linear operators associated to divide and color (DC) models. The degree of nonuniqueness of the random partition yielding a DC model is directly related to the dimension of the kernel of these linear operators. We determine exactly the dimension of these kernels as well as analyze a permutation-invariant version. We also obtain properties of the solution set for certain parameter values which will be important in (1) showing that large threshold discrete Gaussian free fields are DC models and in (2) analyzing when the Ising model with a positive external field is a DC model, both in future work. However, even here, we give an application to the Ising model on a triangle.


title: The spherical multipole resonance probe: kinetic damping in its spectrum

abstract: The multipole resonance probe is one of the recently developed measurement devices to measure plasma parameter like electron density and temperature based on the concept of active plasma resonance spectroscopy. The dynamical interaction between the probe and the plasma in electrostatic, kinetic description can be modeled in an abstract notation based on functional analytic methods. These methods provide the opportunity to derive a general solution, which is given as the response function of the probe-plasma system. It is defined by the matrix elements of the resolvent of an appropriate dynamical operator. Based on the general solution a residual damping for vanishing pressure can be predicted and can only be explained by kinetic effects. Within this manuscript an explicit response function of the multipole resonance probe is derived. Therefore, the resolvent is determined by its algebraic representation based on an expansion in orthogonal basis functions. This allows to compute an approximated response function and its corresponding spectra, which show additional damping due to kinetic effects.


title: The reliability of a deep learning model in clinical out-of-distribution MRI data: a multicohort study

abstract: Deep learning (DL) methods have in recent years yielded impressive results in medical imaging, with the potential to function as clinical aid to radiologists. However, DL models in medical imaging are often trained on public research cohorts with images acquired with a single scanner or with strict protocol harmonization, which is not representative of a clinical setting. The aim of this study was to investigate how well a DL model performs in unseen clinical data sets---collected with different scanners, protocols and disease populations---and whether more heterogeneous training data improves generalization. In total, 3117 MRI scans of brains from multiple dementia research cohorts and memory clinics, that had been visually rated by a neuroradiologist according to Scheltens' scale of medial temporal atrophy (MTA), were included in this study. By training multiple versions of a convolutional neural network on different subsets of this data to predict MTA ratings, we assessed the impact of including images from a wider distribution during training had on performance in external memory clinic data. Our results showed that our model generalized well to data sets acquired with similar protocols as the training data, but substantially worse in clinical cohorts with visibly different tissue contrasts in the images. This implies that future DL studies investigating performance in out-of-distribution (OOD) MRI data need to assess multiple external cohorts for reliable results. Further, by including data from a wider range of scanners and protocols the performance improved in OOD data, which suggests that more heterogeneous training data makes the model generalize better. To conclude, this is the most comprehensive study to date investigating the domain shift in deep learning on MRI data, and we advocate rigorous evaluation of DL models on clinical data prior to being certified for deployment.

categories: physics.med-ph cs.CV cs.LG eess.IV q-bio.QM stat.ML


title: MaxSAT Evaluation 2019 -- Benchmark: Identifying Security-Critical Cyber-Physical Components in Weighted AND/OR Graphs

abstract: This paper presents a MaxSAT benchmark focused on identifying critical nodes in AND/OR graphs. We use AND/OR graphs to model Industrial Control Systems (ICS) as they are able to semantically grasp intricate logical interdependencies among ICS components. However, identifying critical nodes in AND/OR graphs is an NP-complete problem. We address this problem by efficiently transforming the input AND/OR graph-based model into a weighted logical formula that is then used to build and solve a Weighted Partial MaxSAT problem. The benchmark includes 80 cases with AND/OR graphs of different size and composition as well as the optimal cost and solution for each case.


title: The halo mass function of late-type galaxies from HI kinematics

abstract: We present an empirical method to measure the halo mass function (HMF) of galaxies. We determine the relation between the \hi\ line-width from single-dish observations and the dark matter halo mass ($M_{200}$) inferred from rotation curve fits in the SPARC database, then we apply this relation to galaxies from the \hi\ Parkes All Sky Survey (HIPASS) to derive the HMF. This empirical HMF is well fit by a Schecther function, and matches that expected in $\Lambda$CDM over the range $10^{10.5} < M_{200} < 10^{12}\;\mathrm{M}_{\odot}$. More massive halos must be poor in neutral gas to maintain consistency with the power law predicted by $\Lambda$CDM. We detect no discrepancy at low masses. The lowest halo mass probed by HIPASS, however, is just greater than the mass scale where the Local Group missing satellite problem sets in. The integrated mass density associated with the dark matter halos of \hi-detected galaxies sums to $\Omega_{\rm m,gal} \approx 0.03$ over the probed mass range.


title: Spiral arms in the proto-planetary disc HD100453 detected with ALMA: evidence for binary-disc interaction and a vertical temperature gradient

abstract: Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.


title: Bulk Reconstruction Beyond the Entanglement Wedge

abstract: We study the portion of an asymptotically Anti de Sitter geometry's bulk where the metric can be reconstructed, given the areas of minimal 2-surfaces anchored to a fixed boundary subregion. We exhibit situations in which this region can reach parametrically far outside of the entanglement wedge. If the setting is furthermore holographic, so that the bulk geometry is dual to a state in a conformal field theory (CFT), these minimal 2-surface areas can be deduced from the expectation values of operators localized within the boundary subregion. This presents us with an alternative: Either the reduced CFT state encodes significant information about the bulk beyond the entanglement wedge, challenging conventional intuition about holographic subregion duality; or the reduced CFT state fails to contain information about operators whose expectation values give the areas of minimal 2-surfaces anchored within that subregion, challenging conventional intuition about the holographic dictionary.


title: Multi-planet disc interactions in binary systems

abstract: We investigate the evolution of a multi--planet--disc system orbiting one component of a binary star system. The planet--disc system is initially coplanar but misaligned to the binary orbital plane. The planets are assumed to be giants that open gaps in the disc. We first study the role of the disc in shaping the mutual evolution of the two planets using a secular model for low initial tilt. In general we find that the planets and the disc do not remain coplanar, in agreement with Lubow & Martin (2016) for the single planet case. Instead, the planets and the disc undergo tilt oscillations. A high mass disc between the two planets causes the planets and the disc to nodally precess at the same average rate but they are generally misaligned. The amplitude of the tilt oscillations between the planets is larger while the disc is present. We then consider higher initial tilts using hydrodynamical simulations and explore the possibility of the formation of eccentric Kozai-Lidov (KL) planets. We find that the inner planet's orbit undergoes eccentricity growth for a large range of disc masses and initial misalignments. For a low disc mass and large initial misalignment both planets and the disc can undergo KL oscillations. Furthermore, we find that sufficiently massive discs can cause the inner planet to increase its inclination beyond $90^{\circ}$ and therefore to orbit the binary in a retrograde fashion. The results have important implications for the explanation of very eccentric planets and retrograde planets observed in multi-planet systems.


title: Projected velocity statistics of interstellar turbulence

abstract: Velocity statistics is a direct probe of the dynamics of interstellar turbulence. Its observational measurements are very challenging due to the convolution between density and velocity and projection effects. We introduce the projected velocity structure function, which can be generally applied to statistical studies of both sub- and super-sonic turbulence in different interstellar phases. It recovers the turbulent velocity spectrum from the projected velocity field in different regimes, and when the thickness of a cloud is less than the driving scale of turbulence, it can also be used to determine the cloud thickness and the turbulence driving scale. By applying it to the existing core velocity dispersion measurements of the Taurus cloud, we find a transition from the Kolmogorov to the Burgers scaling of turbulent velocities with decreasing length scales, corresponding to the large-scale solenoidal motions and small-scale compressive motions, respectively. The latter occupy a small fraction of the volume and can be selectively sampled by clusters of cores with the typical cluster size indicated by the transition scale.


title: Bootstrapping Mixed Correlators in Three-Dimensional Cubic Theories II

abstract: Conformal field theories (CFTs) with cubic global symmetry in 3D are relevant in a variety of condensed matter systems and have been studied extensively with the use of perturbative methods like the $\varepsilon$ expansion. In an earlier work, we used the nonperturbative numerical conformal bootstrap to provide evidence for the existence of a previously unknown 3D CFT with cubic symmetry, dubbed "Platonic CFT". In this work, we make further use of the numerical conformal bootstrap to perform a three-dimensional scan in the space of scaling dimensions of three low-lying operators. We find a three-dimensional isolated allowed region in parameter space, which includes both the 3D (decoupled) Ising model and the Platonic CFT. The essential assumptions on the spectrum of operators used to provide the isolated allowed region include the existence of a stress-energy tensor and the irrelevance of certain operators (in the renormalization group sense).


title: What Gets Echoed? Understanding the "Pointers" in Explanations of Persuasive Arguments

abstract: Explanations are central to everyday life, and are a topic of growing interest in the AI community. To investigate the process of providing natural language explanations, we leverage the dynamics of the /r/ChangeMyView subreddit to build a dataset with 36K naturally occurring explanations of why an argument is persuasive. We propose a novel word-level prediction task to investigate how explanations selectively reuse, or echo, information from what is being explained (henceforth, explanandum). We develop features to capture the properties of a word in the explanandum, and show that our proposed features not only have relatively strong predictive power on the echoing of a word in an explanation, but also enhance neural methods of generating explanations. In particular, while the non-contextual properties of a word itself are more valuable for stopwords, the interaction between the constituent parts of an explanandum is crucial in predicting the echoing of content words. We also find intriguing patterns of a word being echoed. For example, although nouns are generally less likely to be echoed, subjects and objects can, depending on their source, be more likely to be echoed in the explanations.


title: Contributions to the Formalization and Extraction of Generic Bases of Association Rules

abstract: In this thesis, a detailed study shows that closed itemsets and minimal generators play a key role for concisely representing both frequent itemsets and association rules. These itemsets structure the search space into equivalence classes such that each class gathers the itemsets appearing in the same subset aka objects or transactions of the given data. In this respect, we proposed lossless reductions of the minimal generator set thanks to a new substitution-based process. Our theoretical results are extended to the association rule framework in order to reduce as much as possible the number of retained rules without information loss. We then give a thorough formal study of the related inference mechanism allowing to derive all redundant association rules, starting from the retained ones. We also lead a thorough exploration of the disjunctive search space, where itemsets are characterized by their respective disjunctive supports, instead of the conjunctive ones. This exploration is motivated by the fact that, in some applications, such information brings richer knowledge to the end-users. To obtain a redundancy free representation of the disjunctive search space, an interesting solution consists in selecting a unique element to represent itemsets covering the same set of data. Two itemsets are equivalent if their respective items cover the same set of data. In this regard, we introduced a new operator dedicated to this task. In each induced equivalence class, minimal elements are called essential itemsets, while the largest one is called disjunctive closed itemset. The introduced operator is then at the roots of new concise representations of frequent itemsets. We also exploit the disjunctive search space to derive generalized association rules. These latter rules generalize classic ones to also offer disjunction and negation connectors between items, in addition to the conjunctive one.


title: MCTDH-X: The multiconfigurational time-dependent Hartree method for indistinguishable particles software

abstract: We introduce and describe the multiconfigurational time-depenent Hartree for indistinguishable particles (MCTDH-X) software. This powerful tool allows the investigation of ground state properties and dynamics of interacting quantum many-body systems in different spatial dimensions. The MCTDH-X software is a set of programs and scripts to compute, analyze, and visualize solutions for the time-dependent and time-independent many-body Schr\"{o}dinger equation for indistinguishable quantum particles. As the MCTDH-X software represents a general solver for the Schr\"{o}dinger equation, it is applicable to a wide range of problems in the fields of atomic, optical, molecular physics as well as condensed matter systems. In particular, it can be used to study light-matter interactions, correlated dynamics of electrons, as well as some aspects related to quantum information and computing. The MCTDH-X software solves a set of non-linear coupled working equations based on the application of the variational principle to the Schr\"{o}dinger equation. These equations are obtained by using an ansatz for the many-body wavefunction that is a time-dependent expansion in a set of time-dependent many-body basis states. The time-dependence of the basis set enables MCTDH-X to deal with quantum dynamics at a superior accuracy as compared to, for instance, exact diagonalization approaches. Herein, we give an introduction to the MCTDH-X software via an easy-to-follow tutorial with a focus on accessibility. We use the double well to illustrate the fermionization of bosonic particles, the crystallization of fermionic particles, characteristics of the superfluid and Mott-insulator quantum phases in Hubbard models, and even dynamical quantum phase transitions. Our tutorial guides the potential user to apply the MCTDH-X software also to more complex systems.


title: Automated Assignment of Backbone Resonances Using Residual Dipolar Couplings Acquired from a Protein with Known Structure

abstract: Resonance assignment is a critical first step in the investigation of protein structures using NMR spectroscopy. The development of assignment methods that require less experimental data is possible with prior knowledge of the macromolecular structure. Automated methods of performing the task of resonance assignment can significantly reduce the financial cost and time requirement for protein structure determination. Such methods can also be beneficial in validating a protein's solution state structure. Here we present a new approach to the assignment problem. Our approach uses only RDC data to assign backbone resonances. It provides simultaneous order tensor estimation and assignment. Our approach compares independent order tensor estimates to determine when the correct order tensor has been found. We demonstrate the algorithm's viability using simulated data from the protein domain 1A1Z.


title: Memory Requirement Reduction of Deep Neural Networks Using Low-bit Quantization of Parameters

abstract: Effective employment of deep neural networks (DNNs) in mobile devices and embedded systems is hampered by requirements for memory and computational power. This paper presents a non-uniform quantization approach which allows for dynamic quantization of DNN parameters for different layers and within the same layer. A virtual bit shift (VBS) scheme is also proposed to improve the accuracy of the proposed scheme. Our method reduces the memory requirements, preserving the performance of the network. The performance of our method is validated in a speech enhancement application, where a fully connected DNN is used to predict the clean speech spectrum from the input noisy speech spectrum. A DNN is optimized and its memory footprint and performance are evaluated using the short-time objective intelligibility, STOI, metric. The application of the low-bit quantization allows a 50% reduction of the DNN memory footprint while the STOI performance drops only by 2.7%.


title: Mass spectrum of scalar and gauge fields in the cosmological multiply warped brane-world

abstract: In this paper, we compute the massive modes of scalar and gauge fields in a multiply warped brane-world scenario in six-dimensional spacetime with a Freedman-Robertson-Walker (FRW) metric. The massive Kaluza-Klein (KK) modes of the fields are split in two mass towers due to the presence of the warped in both dimensions. This does not occurs in the standard Randall-Sundrum (RS) model. We investigate the dependence of the mass spectrum on Hubble constant. Lighter massive modes of the fields are found, whereby it is possible to infer that scalar and bosonic matter can directly influence in universe expansion.


title: On the logarithmic coefficient of the entanglement entropy of a Maxwell field

abstract: We elucidate the mismatch between the $A$-anomaly coefficient and the coefficient of the logarithmic term in the entanglement entropy of a Maxwell field. In contrast to the usual assumptions about the protection of renormalization group charges at the infrared, the logarithmic term is different for a free Maxwell field and a Maxwell field interacting with heavy charges. This is possible because of the presence of superselection sectors in the IR theory. However, the correction due to the coupling with charged vacuum fluctuations, that restores the anomaly coefficient, is independent of the precise UV dynamics. The problem is invariant under electromagnetic duality, and the solution requires both the existence of electric charges and magnetic monopoles. We use a real-time operator approach but we also show how the results for the free and interacting fields are translated into an effective correction to the four-sphere partition function.


title: Atomistic $T$-matrix theory of disordered 2D materials: Bound states, spectral properties, quasiparticle scattering, and transport

abstract: In this work, we present an atomistic first-principles framework for modeling the low-temperature electronic and transport properties of disordered two-dimensional (2D) materials with randomly distributed point defects (impurities). The method is based on the $T$-matrix formalism in combination with realistic density-functional theory (DFT) descriptions of the defects and their scattering matrix elements. From the $T$-matrix approximations to the disorder-averaged Green's function (GF) and the collision integral in the Boltzmann transport equation, the method allows calculations of, e.g., the density of states (DOS) including contributions from bound defect states, the quasiparticle spectrum and the spectral linewidth (scattering rate), and the conductivity/mobility of disordered 2D materials. We demonstrate the method by examining these quantities in monolayers of the archetypal 2D materials graphene and transition metal dichalcogenides (TMDs) contaminated with vacancy defects and substitutional impurity atoms. By comparing the Born and $T$-matrix approximations, we also demonstrate a strong breakdown of the Born approximation for defects in 2D materials manifested in a pronounced renormalization of, e.g., the scattering rate by the higher-order $T$-matrix method. As the $T$-matrix approximation is essentially exact for dilute disorder, i.e., low defect concentrations ($c_\text{dis} \ll 1$) or density ($n_\text{dis}\ll A_\text{cell}^{-1}$ where $A_\text{cell}$ is the unit cell area), our first-principles method provides an excellent framework for modeling the properties of disordered 2D materials with defect concentrations relevant for devices.


title: Radio-relic and the diffuse emission trail discovered in a low mass galaxy cluster Abell 1697

abstract: We report the discovery of a putative radio relic, 830 kpc in length and found toward the outskirts of galaxy cluster Abell 1697 ($z=0.181$), using the LOFAR Two Meter Sky Survey (LoTSS) at 144 MHz. With an X-ray-inferred mass of $M_{500}^{X-ray}=2.9^{+0.8}_{-0.7}\times10^{14}~\rm{M_{\odot}}$, this places Abell 1697 among the least massive relic hosts. The relic is also detected at 325 MHz in the Westerbork Northern Sky Survey (WENSS) and at 1.4 GHz in the NRAO VLA Sky Survey (NVSS) with an average spectral index of $\alpha(144,325,1400~\rm{MHz})=-0.98\pm0.01$ and magnetic field of $B_{eq}\sim0.6~\mu$G. This relic, located in the northeast periphery of the cluster, is 300 kpc wide, exhibits a gradual spectral steepening across the width ($\alpha_{144 \rm{MHz}}^{1.4\rm{GHz}}(inj)=-0.70\pm0.11$ to $\alpha_{144 \rm{MHz}}^{1.4\rm{GHz}}(edge)=-1.19\pm0.15$), as well as indications of a co-spatial X-ray (ROSAT) shock and the radio relic emission. The radio power of the relic is $P_{1.4GHz}=8.5\pm1.1\times10^{23}~\rm{W\;Hz^{-1}}$, which is found to be in good agreement with the expected empirical correlation between the radio power and Largest Linear Size (LLS) of relics. The relic is trailed by extended ($790\times550$ kpc) diffuse radio emission towards the cluster center, that is likely an ultra-steep spectrum ($\alpha_{144 \rm{MHz}}^{1.4\rm{GHz}}<-1.84$) radio source. This structure is also found to be older by at least 190 Myrs, has a very low surface brightness of $0.3~\mu$Jy arcsec$^{-2}$ and magnetic field $B_{eq}\sim0.8~\mu$G, similar to that of a radio phoenix. Finally, we discuss the possible mechanisms responsible for the relic and the trailing diffuse radio emission, invoking re-acceleration due to wake turbulence, as well as the revival of fossil electrons from an old AGN activity by the cluster merger shocks.


title: Axion-photon-dark photon oscillation and its implication for 21 cm observation

abstract: We examine the resonant conversion of axion-like particle (ALP) or dark photon to the electromagnetic photon in the early Universe, which takes place due to the ALP-photon-dark photon oscillations in background dark photon gauge fields. It is noted that the corresponding conversion probability can have an unusual spectral feature which allows strong conversion at low frequency domain, but has negligible conversion at high frequencies above certain critical frequency which is determined by the ALP coupling to dark photon and the strength of background dark photon gauge field. We apply this scheme to heat up the 21 cm photons without affecting the Cosmic Microwave Background, which can explain the tentative absorption signal of 21 cm photons detected recently by the EDGES experiment.


title: Traces of Newton-Sobolev, Hajlasz-Sobolev and BV functions on metric spaces

abstract: We study the boundary traces of Newton-Sobolev, Hajlasz-Sobolev, and BV (bounded variation) functions. Assuming less regularity of the domain than is usually done in the literature, we show that all of these function classes achieve the same "boundary values", which in particular implies that the trace spaces coincide provided that they exist. Many of our results seem to be new even in Euclidean spaces but we work in a more general complete metric space equipped with a doubling measure and supporting a Poincare inequality.


title: (the struggle) Towards an open source policy

abstract: Public availability and tracability of results from publically-funded work is a topic that gets more and more attention from funding agencies and scientific policy makers. However, most policies focus on data as the output of research. In this contribution, we focus on research software and we introduce the ASTRON Open Source Policy. Apart from the license used (Apache 2.0), the policy is written as a manual that explains how to license software, when to assign a Digital Object Identifier (DOI), and defines that all code should be put in an ASTRON managed repository. The policy has been made publically available, a DOI has been assigned to it and it has been put in a repository to stimulate the ADASS community to start a conversation on how to make our code publically accessible and citable.


title: Think-aloud interviews: A tool for exploring student statistical reasoning

abstract: Think-aloud interviews have been a valuable but underused tool in statistics education research. Think-alouds, in which students narrate their reasoning in real time while solving problems, differ in important ways from other types of cognitive interviews and related education research methods. Beyond the uses already found in the statistics literature -- mostly validating the wording of statistical concept inventory questions and studying student misconceptions -- we suggest other possible use cases for think-alouds and summarize best-practice guidelines for designing think-aloud interview studies. Using examples from our own experiences studying the local student body for our introductory statistics courses, we illustrate how research goals should inform study-design decisions and what kinds of insights think-alouds can provide. We hope that our overview of think-alouds encourages more statistics educators and researchers to begin using this method.


title: DialoGPT: Large-Scale Generative Pre-training for Conversational Response Generation

abstract: We present a large, tunable neural conversational response generation model, DialoGPT (dialogue generative pre-trained transformer). Trained on 147M conversation-like exchanges extracted from Reddit comment chains over a period spanning from 2005 through 2017, DialoGPT extends the Hugging Face PyTorch transformer to attain a performance close to human both in terms of automatic and human evaluation in single-turn dialogue settings. We show that conversational systems that leverage DialoGPT generate more relevant, contentful and context-consistent responses than strong baseline systems. The pre-trained model and training pipeline are publicly released to facilitate research into neural response generation and the development of more intelligent open-domain dialogue systems.


title: Prompt X-ray emission from Fast Radio Bursts -- Upper limits with AstroSat

abstract: Fast Radio Bursts (FRBs) are short lived ($\sim$ msec), energetic transients (having a peak flux density of $\sim$ Jy) with no known prompt emission in other energy bands. We present results of a search for prompt X-ray emissions from 41 FRBs using the Cadmium Zinc Telluride Imager (CZTI) on AstroSat which continuously monitors $\sim70\%$ of the sky. Our searches on various timescales in the 20-200 keV range, did not yield any counterparts in this hard X-ray band. We calculate upper limits on hard X-ray flux, in the same energy range and convert them to upper bounds for $\eta$: the ratio X-ray to radio fluence of FRBs. We find $\eta \leq 10^{8-10}$ for hard X-ray emission. Our results will help constrain the theoretical models of FRBs as the models become more quantitative and nearer, brighter FRBs are discovered.


title: Optimality of Spectral Clustering in the Gaussian Mixture Model

abstract: Spectral clustering is one of the most popular algorithms to group high dimensional data. It is easy to implement and computationally efficient. Despite its popularity and successful applications, its theoretical properties have not been fully understood. In this paper, we show that spectral clustering is minimax optimal in the Gaussian Mixture Model with isotropic covariance matrix, when the number of clusters is fixed and the signal-to-noise ratio is large enough. Spectral gap conditions are widely assumed in the literature to analyze spectral clustering. On the contrary, these conditions are not needed to establish optimality of spectral clustering in this paper.


title: Spatial structure of lasing modes in wave-chaotic semiconductor microcavities

abstract: We present experimental and numerical studies of broad-area semiconductor lasers with chaotic ray dynamics. The emission intensity distributions at the cavity boundaries are measured and compared to ray tracing simulations and numerical calculations of the passive cavity modes. We study two different cavity geometries, a D-cavity and a stadium, both of which feature fully chaotic ray dynamics. While the far-field distributions exhibit fairly homogeneous emission in all directions, the emission intensity distributions at the cavity boundary are highly inhomogeneous, reflecting the non-uniform intensity distributions inside the cavities. The excellent agreement between experiments and simulations demonstrates that the intensity distributions of wave-chaotic semiconductor lasers are primarily determined by the cavity geometry. This is in contrast to conventional Fabry-Perot broad-area lasers for which the intensity distributions are to a large degree determined by the nonlinear interaction of the lasing modes with the semiconductor gain medium.


title: Three-dimensional general relativistic Poynting-Robertson effect II: Radiation field from a rigidly rotating spherical source

abstract: We investigate the three-dimensional, general relativistic Poynting-Robertson effect in the case of rigidly rotating spherical source which emits radiation radially in the local comoving frame. Such radiation field is meant to approximate the field produced by the surface of a rotating neutron star, or by the central radiating hot corona of accreting black holes; it extends the purely radial radiation field that we considered in a previous study. Its angular momentum is expressed in terms of the rotation frequency and radius of the emitting source. For the background we adopt a Kerr spacetime geometry. We derive the equations of motion for test particles influenced by such radiation field, recovering the classical and weak-field approximation for slow rotation. We concentrate on solutions consisting of particles orbiting along circular orbits off and parallel to the equatorial plane, which are stabilized by the balance between gravitational attraction, radiation force and PR drag. Such solutions are found to lie on a critical hypersurface, whose shape may morph from prolate to oblate depending on the Kerr spin parameter and the luminosity, rotation and radius of the radiating sphere. For selected parameter ranges, the critical hypersurface intersects the radiating sphere giving rise to a bulging equatorial region or, alternatively, two lobes above the poles. We calculate the trajectories of test particles in the close vicinity of the critical hypersurface for a selected set of initial parameters and analyze the spatial and angular velocity of test particles captured on the critical hypersurface.


title: Topical Collection: Testing the Kerr spacetime with gravitational-wave and electromagnetic observations

abstract: One of the Holy Grails of observational astronomy is to confirm the prediction that black holes in the Universe are described by the Kerr solution of Einstein's field equations of general relativity. This Topical Collection provides a status report of theoretical and experimental progress towards confirming the "Kerr paradigm" through X-ray astronomy, gravitational lensing, stellar tidal disruption events, superradiance, and gravitational-wave observations of black hole binary mergers.


title: Sharp regularity for degenerate obstacle type problems: a geometric approach

abstract: We prove sharp regularity estimates for solutions of obstacle type problems driven by a class of degenerate fully nonlinear operators; more specifically, we consider viscosity solutions of \[ |D u|^\gamma F(x, D^2u) = f(x)\chi_{\{u>\phi\}} \textrm{ in } B_1 \] with $\gamma>0$, $\phi \in C^{1, \alpha}(B_1)$ for some $\alpha\in(0,1]$ and $f\in L^\infty(B_1)$ constrained to satisfy \[ u\geq \phi\textrm{ in } B_1 \] and prove that they are $C^{1,\beta}(B_{1/2})$ (and in particular along free boundary points) where $\beta=\min\left\{\alpha, \frac{1}{\gamma+1}\right\}$. Moreover, we achieve such a feature by using a recently developed geometric approach which is a novelty for these kind of free boundary problems. Further, under a natural non-degeneracy assumption on the obstacle, we prove that the free boundary $\partial\{u>\phi\}$ has zero Lebesgue measure. Our results are new even for seemingly simple model as follows \[ |Du|^\gamma \Delta u=\chi_{\{u>\phi\}} \quad \text{with}\quad \gamma>0. \]


title: The contribution from stars stripped in binaries to cosmic reionization of hydrogen and helium

abstract: Massive stars are often found in binary systems and it has been argued that binary products boost the ionizing radiation of stellar populations. Accurate predictions for binary products are needed to understand and quantify their contribution to Cosmic Reionization. We investigate the contribution of stars stripped in binaries since (1) they are, arguably, the best-understood products of binary evolution, (2) we recently produced the first non-LTE radiative transfer calculations for the atmospheres of these stripped stars that predict their ionizing spectra, and (3) they are very promising sources since they boost the ionizing emission of stellar populations at late times. This allows stellar feedback to clear the surroundings such that a higher fraction of their photons can escape and ionize the intergalactic medium. Combining our detailed predictions for the ionizing spectra with a simple cosmic reionization model, we estimate that stripped stars contributed tens of percent of the photons that caused cosmic reionization of hydrogen, depending on the assumed escape fractions. More importantly, stripped stars harden the ionizing emission. We estimate that the spectral index for the ionizing part of the spectrum can increase to -1 compared to <-2 for single stars. At high redshift, stripped stars and massive single stars combined dominate the HeII-ionizing emission, but we expect active galactic nuclei drive cosmic helium reionization. Further observational consequences we expect are (1) high ionization states for the intergalactic gas surrounding stellar systems, such as CIV and SiIV and (2) additional heating of the intergalactic medium of up to a few thousand Kelvin. Quantifying these warrants the inclusion of accurate models for stripped stars and other binary products in full cosmological simulations.


title: Dark Matter Cores and Cusps in Spiral Galaxies and their Explanations

abstract: We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.


title: Multi-octave, CEP-stable source for high-energy field synthesis

abstract: The development of high-energy, high-power, multi-octave light-transients is currently subject of intense research driven by emerging applications in attosecond spectroscopy and coherent control. We report on a phase-stable, multi-octave source based on a Yb:YAG amplifier for light-transient generation. We demonstrate the amplification of a two-octave spectrum to 25\,$\mu$J of energy in two broadband amplification channels and their temporal compression to 6\,fs and 18\,fs at 1\,$\mu$m and 2\,$\mu$m, respectively. In this scheme due to the intrinsic temporal synchronization between the pump and seed pulses, the temporal jitter is restricted to long-term drift. We show that the intrinsic stability of the synthesizer allows for sub-cycle detection of an electric field at 0.15\,PHz. The complex electric field of the 0.15\,PHz pulses and their free-induction decay after interaction with water molecules are resolved by electro-optic sampling over 2\,ps. The scheme is scalable in peak- and average-power.


title: Stern-Gerlach: conceptually clean or acceptably vague?

abstract: This paper develops a number of quantum mechanical characterisations of Stern-Gerlach. It discusses areas of vagueness in their formulation. Philosophers criticise quantum mechanics for unacceptable vagueness in connection with the measurement problem. The quantum formulation problems identified by this paper go beyond the locus of philosophical criticism. It concludes with an open question, are some areas of vagueness in quantum mechanics more acceptable philosophically than others and, if so, why?


title: Uncover Sexual Harassment Patterns from Personal Stories by Joint Key Element Extraction and Categorization

abstract: The number of personal stories about sexual harassment shared online has increased exponentially in recent years. This is in part inspired by the \#MeToo and \#TimesUp movements. Safecity is an online forum for people who experienced or witnessed sexual harassment to share their personal experiences. It has collected \textgreater 10,000 stories so far. Sexual harassment occurred in a variety of situations, and categorization of the stories and extraction of their key elements will provide great help for the related parties to understand and address sexual harassment. In this study, we manually annotated those stories with labels in the dimensions of location, time, and harassers' characteristics, and marked the key elements related to these dimensions. Furthermore, we applied natural language processing technologies with joint learning schemes to automatically categorize these stories in those dimensions and extract key elements at the same time. We also uncovered significant patterns from the categorized sexual harassment stories. We believe our annotated data set, proposed algorithms, and analysis will help people who have been harassed, authorities, researchers and other related parties in various ways, such as automatically filling reports, enlightening the public in order to prevent future harassment, and enabling more effective, faster action to be taken.


title: A Framework to Explore Workload-Specific Performance and Lifetime Trade-offs in Neuromorphic Computing

abstract: Neuromorphic hardware with non-volatile memory (NVM) can implement machine learning workload in an energy-efficient manner. Unfortunately, certain NVMs such as phase change memory (PCM) require high voltages for correct operation. These voltages are supplied from an on-chip charge pump. If the charge pump is activated too frequently, its internal CMOS devices do not recover from stress, accelerating their aging and leading to negative bias temperature instability (NBTI) generated defects. Forcefully discharging the stressed charge pump can lower the aging rate of its CMOS devices, but makes the neuromorphic hardware unavailable to perform computations while its charge pump is being discharged. This negatively impacts performance such as latency and accuracy of the machine learning workload being executed. In this paper, we propose a novel framework to exploit workload-specific performance and lifetime trade-offs in neuromorphic computing. Our framework first extracts the precise times at which a charge pump in the hardware is activated to support neural computations within a workload. This timing information is then used with a characterized NBTI reliability model to estimate the charge pump's aging during the workload execution. We use our framework to evaluate workload-specific performance and reliability impacts of using 1) different SNN mapping strategies and 2) different charge pump discharge strategies. We show that our framework can be used by system designers to explore performance and reliability trade-offs early in the design of neuromorphic hardware such that appropriate reliability-oriented design margins can be set.


title: Enhanced Seeing Mode at the LBT: A Method to Significantly Improve Angular Resolution over a 4' x 4' Field of View

abstract: Since 2014, the LBT's First Light Adaptive Optics (FLAO) system has also included a seldom used capability, known as Enhanced Seeing Mode (ESM), that can improve the angular resolution over a 4' x 4' field of view (FOV). In full AO operation, FLAO provides diffraction limited (DL) capabilities over a small (30" x 30") FOV. By comparison, ESM can achieve significantly enhanced resolution, over natural seeing, across a far larger FOV. This improves operational efficiency over standard seeing limited (SL) observations and is applicable across a broader range of scientific targets. ESM uses 11 modes of correction (including tip and tilt) to remove residual aberrations and jitter which significantly improves angular resolution over the full FOV. While this mode does not reach the DL, it can achieve uniform angular resolutions as good as 0".22 over the FOV. Furthermore, it allows for the use of multi-object spectroscopy with R~10,000 or imaging with angular resolution similar to that achieved by the Wide-Field Camera 3 infrared channel on the Hubble Space Telescope, but powered by 11.6 meters of effective aperture in binocular mode. As part of the on-going characterization of ESM, we have demonstrated that even in poor seeing conditions (1".5-2") the image quality delivered to the focal station is improved by factors of 2-3. Here, we present the first results of the characterization of ESM, including systematic tests of the delivered PSF across the FOV as a function of the brightness of, and distance from, the AO Reference Star. We present a range of galactic and extra-galactic targets showing the improvements obtained over a broad range of seeing conditions and propose ESM as a standard observational mode for near-Infrared observations.


title: Decoding of visual-related information from the human EEG using an end-to-end deep learning approach

abstract: There is increasing interest in using deep learning approach for EEG analysis as there are still rooms for the improvement of EEG analysis in its accuracy. Convolutional long short-term (CNNLSTM) has been successfully applied in time series data with spatial structure through end-to-end learning. Here, we proposed a CNNLSTM based neural network architecture termed EEG_CNNLSTMNet for the classification of EEG signals in response to grating stimuli with different spatial frequencies. EEG_CNNLSTMNet comprises two convolutional layers and one bidirectional long short-term memory (LSTM) layer. The convolutional layers capture local temporal characteristics of the EEG signal at each channel as well as global spatial characteristics across channels, while the LSTM layer extracts long-term temporal dependency of EEG signals. Our experiment showed that EEG_CNNLSTMNet performed much better at EEG classification than a traditional machine learning approach, i.e. a support vector machine (SVM) with features. Additionally, EEG_CNNLSTMNet outperformed EEGNet, a state-of-art neural network architecture for the intra-subject case. We infer that the underperformance when using an LSTM layer in the inter-subject case is due to long-term dependency characteristics in the EEG signal that vary greatly across subjects. Moreover, the inter-subject fine-tuned classification model using very little data of the new subject achieved much higher accuracy than that trained only on the data from the other subjects. Our study suggests that the fine-tuned inter-subject model can be a potential end-to-end EEG analysis method considering both the accuracy and the required training data of the new subject.


title: A remark on the well-posedness of the modified KdV equation in the Fourier-Lebesgue spaces

abstract: We study the complex-valued modified Korteweg-de Vries equation (mKdV) on the circle. We first consider the real-valued setting and show global well-posedness of the (usual) renormalized mKdV equation in the Fourier-Lebesgue spaces. In the complex-valued setting, we observe that the momentum plays an important role in the well-posedness theory. In particular, we prove that the complex-valued mKdV equation is ill-posed in the sense of non-existence of solutions when the momentum is infinite, in the spirit of the work on the nonlinear Schr\"odinger equation by Guo-Oh (2018). This non-existence result motivates the introduction of the second renormalized mKdV equation, which we propose as the correct model in the complex-valued setting outside of $H^\frac12(\mathbb{T})$. Furthermore, imposing a new notion of finite momentum for the initial data, at low regularity, we show existence of solutions to the complex-valued mKdV equation. In particular, we require an energy estimate, from which conservation of momentum follows.


title: On Magnetic Forces and Work

abstract: We address a long-standing debate over whether classical magnetic forces can do work, ultimately answering the question in the affirmative. In detail, we couple a classical particle with intrinsic spin and elementary dipole moments to the electromagnetic field, derive the appropriate generalization of the Lorentz force law, show that the particle's dipole moments must be collinear with its spin axis, and argue that the magnetic field does mechanical work on the particle's elementary magnetic dipole moment. As consistency checks, we calculate the overall system's energy-momentum and angular momentum, and show that their local conservation equations lead to the same force law and therefore the same conclusions about magnetic forces and work. We also compute the system's Belinfante-Rosenfeld energy-momentum tensor.


title: A tunable High-Q millimeter wave cavity for hybrid circuit and cavity QED experiments

abstract: The millimeter wave (mm-wave) frequency band provides exciting prospects for quantum science and devices, since many high-fidelity quantum emitters, including Rydberg atoms, molecules and silicon vacancies, exhibit resonances near 100 GHz. High-Q resonators at these frequencies would give access to strong interactions between emitters and single photons, leading to rich and unexplored quantum phenomena at temperatures above 1K. We report a 3D mm-wave cavity with a measured single-photon internal quality factor of $3 \times 10^{7}$ and mode volume of $0.14 \times \lambda^3$ at $98.2$ GHz, sufficient to reach strong coupling in a Rydberg cavity QED system. An in-situ piezo tunability of $18$ MHz facilitates coupling to specific atomic transitions. Our unique, seamless and optically accessible resonator design is enabled by the realization that intersections of 3D waveguides support tightly confined bound states below the waveguide cutoff frequency. Harnessing the features of our cavity design, we realize a hybrid mm-wave and optical cavity, designed for interconversion and entanglement of mm-wave and optical photons using Rydberg atoms.


title: Penalized robust estimators in logistic regression with applications to sparse models

abstract: Sparse covariates are frequent in classification and regression problems and in these settings the task of variable selection is usually of interest. As it is well known, sparse statistical models correspond to situations where there are only a small number of non--zero parameters and for that reason, they are much easier to interpret than dense ones. In this paper, we focus on the logistic regression model and our aim is to address robust and penalized estimation for the regression parameter. We introduce a family of penalized weighted $M-$type estimators for the logistic regression parameter that are stable against atypical data. We explore different penalizations functions and we introduce the so--called Sign penalization. This new penalty has the advantage that it depends only on one penalty parameter, avoiding arbitrary tuning constants. We discuss the variable selection capability of the given proposals as well as their asymptotic behaviour. Through a numerical study, we compare the finite sample performance of the proposal corresponding to different penalized estimators either robust or classical, under different scenarios. A robust cross--validation criterion is also presented. The analysis of two real data sets enables to investigate the stability of the penalized estimators to the presence of outliers.


title: The Power Graph of a Torsion-Free Group of Nilpotency Class $2$

abstract: The directed power graph $\mathcal G(\mathbf G)$ of a group $\mathbf G$ is the simple digraph with vertex set $G$ in which $x\rightarrow y$ if $y$ is a power of $x$, the power graph is the underlying simple graph, and the enhanced power graph of $\mathbf G$ is the simple graph with the same vertex set such that two vertices are adjacent if they are powers of some element of $\mathbf G$. In this paper three versions of the definition of the power graphs are discussed, and it is proved that the power graph by any of the three versions of the definitions determines the other two up to isomorphism. It is also proved that, if $\mathbf G$ is a torsion-free group of nilpotency class $2$ and if $\mathbf H$ is a group such that $\mathcal G(\mathbf H)\cong\mathcal G(\mathbf G)$, then $\mathbf G$ and $\mathbf H$ have isomorphic directed power graphs, which was an open problem proposed by Cameron, Guerra and Jurina.


title: Viability of bouncing cosmology in energy-momentum-squared gravity

abstract: We analyze the early-time isotropic cosmology in the so-called energy-momentum-squared gravity (EMSG). In this theory, a $T_{\mu\nu}T^{\mu\nu}$ term is added to the Einstein-Hilbert action, which has been shown to replace the initial singularity by a regular bounce. We show that this is not the case, and the bouncing solution obtained does not describe our Universe since it belongs to a different solution branch. The solution branch that corresponds to our Universe, while nonsingular, is geodesically incomplete. We analyze the conditions for having viable regular-bouncing solutions in a general class of theories that modify gravity by adding higher order matter terms. Applying these conditions on generalizations of EMSG that add a $\left(T_{\mu\nu}T^{\mu\nu}\right)^{n}$ term to the action, we show that the case of $n=5/8$ is the only one that can give a viable bouncing solution, while the $n>5/8$ cases suffer from the same problem as EMSG, i.e. they give nonsingular, geodesically incomplete solutions. Furthermore, we show that the $1/2<n<5/8$ cases can provide a nonsingular initially de Sitter solution. Finally, the expanding, geodesically incomplete branch of EMSG or its generalizations can be combined with its contracting counterpart using junction conditions to provide a (weakly) singular bouncing solution. We outline the junction conditions needed for this extension and provide the extended solution explicitly for EMSG. In this sense, EMSG replaces the standard early-time singularity by a singular bounce instead of a regular one.


title: The full replica symmetry breaking solution in mean-field spin glass models

abstract: This thesis focus on the extension of the Parisi full replica symmetry breaking solution to the Ising spin glass on a random regular graph. We propose a new martingale approach, that overcomes the limits of the Parisi-M\'ezard cavity method, providing a well-defined formulation of the full replica symmetry breaking problem in random regular graphs. We obtain a variational free energy functional, defined by the sum of two variational functionals (auxiliary variational functionals), that are an extension of the Parisi functional of the Sherrington-Kirkpatrick model. We study the properties of the two variational functionals in detailed, providing representation through the solution of a proper backward stochastic differential equation, that generalize the Parisi partial differential equation. Finally, we define the order parameters of the system and get a set of self-consistency equations for the order parameters and free energy.


title: Prediction Modeling and Analysis for Telecom Customer Churn in Two Months

abstract: A practical churn customer prediction model is critical to retain customers for telecom companies in the saturated and competitive market. Previous studies focus on predicting churn customers in current or next month, in which telecom companies don't have enough time to develop and carry out churn management strategies. In this paper, we propose a new T+2 churn customer prediction model, in which the churn customers in two months are recognized and the one-month window T+1 is reserved to carry out churn management strategies. However, the predictions for churn customers in two months are much more difficult than in current or next month because of the weaker correlation between the customer information and churn states. Two characteristics of telecom dataset, the discrimination between churn and non-churn customers is complicated and the class imbalance problem is serious, are observed. To discriminate the churn customers accurately, random forest (RF) classifier is chosen because RF solves the nonlinear separable problem with low bias and low variance and handles high feature spaces and large number of training examples. To overcome the imbalance problem, synthetic minority over-sampling with borderline or tomek link, in which the distribution of the samples remains and the number of the training examples becomes larger, is applied. Overall, a precision ratio of about 50% with a recall ratio of about 50% is achieved in the T+2 churn prediction. The proposed prediction model provides an accurate and operable churn customer prediction model for telecom companies.


title: Effect of color reconnection and rope formation on resonance production in p$-$p collisions in Pythia 8

abstract: The resonance production in proton$-$proton collisions at $\sqrt{s}$ = 7 TeV and 13 TeV have been investigated using Pythia 8 event generator within the framework of microscopic processes like color reconnection and rope hadronization. Specifically, the observable effects of different modes of color reconnections on the ratio of yields of mesonic and baryonic resonances with respect to their stable counterpart have been explored as a function of event activity. A suppression in the ratio is observed as a function of number of multi-partonic interactions for mesonic resonances. The $\mathrm{\phi/K}$ and $\mathrm{\phi/\pi}$ ratios show an enhancement for high multiplicity events due to enhanced production of strange quarks via the microscopic process of rope hadronization in the partonic phase. The mechanism of the hadronization of color ropes together with the QCD-based color reconnection of partons predicted an enhancement in the ratio for baryonic resonances to non-resonance baryons having similar quark content. The yield ratios of resonances are found to be independent of the collision energy and strongly dependent on event activity.


title: Unveiling the merger dynamics of the most massive MaDCoWS cluster at $z = 1.2$ from a multi-wavelength mapping of its intracluster medium properties

abstract: The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)


title: Twin-Finder: Integrated Reasoning Engine for Pointer-related Code Clone Detection

abstract: Detecting code clones is crucial in various software engineering tasks. In particular, code clone detection can have significant uses in the context of analyzing and fixing bugs in large scale applications. However, prior works, such as machine learning-based clone detection, may cause a considerable amount of false positives. In this paper, we propose Twin-Finder, a novel, closed-loop approach for pointer-related code clone detection that integrates machine learning and symbolic execution techniques to achieve precision. Twin-Finder introduces a clone verification mechanism to formally verify if two clone samples are indeed clones and a feedback loop to automatically generated formal rules to tune machine learning algorithm and further reduce the false positives. Our experimental results show that Twin-Finder can swiftly identify up 9X more code clones comparing to a tree-based clone detector, Deckard and remove an average 91.69% false positives.


title: Anderson transition in three-dimensional systems with non-Hermitian disorder

abstract: We study the Anderson transition for three-dimensional (3D) $N \times N \times N$ tightly bound cubic lattices where both real and imaginary parts of onsite energies are independent random variables distributed uniformly between $-W/2$ and $W/2$. Such a non-Hermitian analog of the Anderson model is used to describe random-laser medium with local loss and amplification. We employ eigenvalue statistics to search for the Anderson transition. For 25\% smallest-modulus complex eigenvalues we find the average ratio $r$ of distances to the first and the second nearest neighbor as a function of $W$. For a given $N$ the function $r(W)$ crosses from $0.72$ to 2/3 with a growing $W$ demonstrating a transition from delocalized to localized states. When plotted at different $N$ all $r(W)$ cross at $W_c = 6.0 \pm 0.1$ (in units of nearest neighbor overlap integral) clearly demonstrating the 3D Anderson transition. We find that in the non-Hermitian 2D Anderson model, the transition is replaced by a crossover.


title: Assessing the Privacy Benefits of Domain Name Encryption

abstract: As Internet users have become more savvy about the potential for their Internet communication to be observed, the use of network traffic encryption technologies (e.g., HTTPS/TLS) is on the rise. However, even when encryption is enabled, users leak information about the domains they visit via DNS queries and via the Server Name Indication (SNI) extension of TLS. Two recent proposals to ameliorate this issue are DNS over HTTPS/TLS (DoH/DoT) and Encrypted SNI (ESNI). In this paper we aim to assess the privacy benefits of these proposals by considering the relationship between hostnames and IP addresses, the latter of which are still exposed. We perform DNS queries from nine vantage points around the globe to characterize this relationship. We quantify the privacy gain offered by ESNI for different hosting and CDN providers using two different metrics, the k-anonymity degree due to co-hosting and the dynamics of IP address changes. We find that 20% of the domains studied will not gain any privacy benefit since they have a one-to-one mapping between their hostname and IP address. On the other hand, 30% will gain a significant privacy benefit with a k value greater than 100, since these domains are co-hosted with more than 100 other domains. Domains whose visitors' privacy will meaningfully improve are far less popular, while for popular domains the benefit is not significant. Analyzing the dynamics of IP addresses of long-lived domains, we find that only 7.7% of them change their hosting IP addresses on a daily basis. We conclude by discussing potential approaches for website owners and hosting/CDN providers for maximizing the privacy benefits of ESNI.


title: Infrared dynamics of a light scalar field in de Sitter

abstract: Inertial observers in de Sitter are surrounded by a horizon and see thermal fluctuations. To them, a massless scalar field appears to follow a random motion but any attractive potential, no matter how weak, will eventually stabilize the field. We study this thermalization process in the static patch (the spacetime region accessible to an individual observer) via a truncation to the low frequency spectrum. We focus on the distribution of the field averaged over a subhorizon region. At timescales much longer than the inverse temperature and to leading order in the coupling, we find the evolution to be Markovian, governed by the same Fokker-Planck equation that arises when the theory is studied in the inflationary setup.


title: Dimensionality Reduction and Reduced Order Modeling for Traveling Wave Physics

abstract: We develop an unsupervised machine learning algorithm for the automated discovery and identification of traveling waves in spatio-temporal systems governed by partial differential equations (PDEs). Our method uses sparse regression and subspace clustering to robustly identify translational invariances that can be leveraged to build improved reduced order models (ROMs). Invariances, whether translational or rotational, are well known to compromise the ability of ROMs to produce accurate and/or low-rank representations of the spatio-temporal dynamics. However, by discovering translations in a principled way, data can be shifted into a coordinate systems where quality, low-dimensional ROMs can be constructed. This approach can be used on either numerical or experimental data with or without knowledge of the governing equations. We demonstrate our method on a variety of PDEs of increasing difficulty, taken from the field of fluid dynamics, showing the efficacy and robustness of the proposed approach.


title: Predicting word error rate for reverberant speech

abstract: Reverberation negatively impacts the performance of automatic speech recognition (ASR). Prior work on quantifying the effect of reverberation has shown that clarity (C50), a parameter that can be estimated from the acoustic impulse response, is correlated with ASR performance. In this paper we propose predicting ASR performance in terms of the word error rate (WER) directly from acoustic parameters via a polynomial, sigmoidal, or neural network fit, as well as blindly from reverberant speech samples using a convolutional neural network (CNN). We carry out experiments on two state-of-the-art ASR models and a large set of acoustic impulse responses (AIRs). The results confirm C50 and C80 to be highly correlated with WER, allowing WER to be predicted with the proposed fitting approaches. The proposed non-intrusive CNN model outperforms C50-based WER prediction, indicating that WER can be estimated blindly, i.e., directly from the reverberant speech samples without knowledge of the acoustic parameters.


title: Frequentist Regret Bounds for Randomized Least-Squares Value Iteration

abstract: We consider the exploration-exploitation dilemma in finite-horizon reinforcement learning (RL). When the state space is large or continuous, traditional tabular approaches are unfeasible and some form of function approximation is mandatory. In this paper, we introduce an optimistically-initialized variant of the popular randomized least-squares value iteration (RLSVI), a model-free algorithm where exploration is induced by perturbing the least-squares approximation of the action-value function. Under the assumption that the Markov decision process has low-rank transition dynamics, we prove that the frequentist regret of RLSVI is upper-bounded by $\widetilde O(d^2 H^2 \sqrt{T})$ where $ d $ are the feature dimension, $ H $ is the horizon, and $ T $ is the total number of steps. To the best of our knowledge, this is the first frequentist regret analysis for randomized exploration with function approximation.


title: Goals, Process, and Challenges of Exploratory Data Analysis: An Interview Study

abstract: How do analysis goals and context affect exploratory data analysis (EDA)? To investigate this question, we conducted semi-structured interviews with 18 data analysts. We characterize common exploration goals: profiling (assessing data quality) and discovery (gaining new insights). Though the EDA literature primarily emphasizes discovery, we observe that discovery only reliably occurs in the context of open-ended analyses, whereas all participants engage in profiling across all of their analyses. We describe the process and challenges of EDA highlighted by our interviews. We find that analysts must perform repetitive tasks (e.g., examine numerous variables), yet they may have limited time or lack domain knowledge to explore data. Analysts also often have to consult other stakeholders and oscillate between exploration and other tasks, such as acquiring and wrangling additional data. Based on these observations, we identify design opportunities for exploratory analysis tools, such as augmenting exploration with automation and guidance.


title: Mitigating the Effects of Non-Identifiability on Inference for Bayesian Neural Networks with Latent Variables

abstract: Bayesian Neural Networks with Latent Variables (BNN+LVs) capture predictive uncertainty by explicitly modeling model uncertainty (via priors on network weights) and environmental stochasticity (via a latent input noise variable). In this work, we first show that BNN+LV suffers from a serious form of non-identifiability: explanatory power can be transferred between the model parameters and latent variables while fitting the data equally well. We demonstrate that as a result, in the limit of infinite data, the posterior mode over the network weights and latent variables is asymptotically biased away from the ground-truth. Due to this asymptotic bias, traditional inference methods may in practice yield parameters that generalize poorly and misestimate uncertainty. Next, we develop a novel inference procedure that explicitly mitigates the effects of likelihood non-identifiability during training and yields high-quality predictions as well as uncertainty estimates. We demonstrate that our inference method improves upon benchmark methods across a range of synthetic and real data-sets.


title: MPro: Combining Static and Symbolic Analysis for Scalable Testing of Smart Contract

abstract: Smart contracts are executable programs that enable the building of a programmable trust mechanism between multiple entities without the need of a trusted third-party. Researchers have developed several security scanners in the past couple of years. However, many of these analyzers either do not scale well, or if they do, produce many false positives. This issue is exacerbated when bugs are triggered only after a series of interactions with the functions of the contract-under-test. A depth-n vulnerability, refers to a vulnerability that requires invoking a specific sequence of n functions to trigger. Depth-n vulnerabilities are time-consuming to detect by existing automated analyzers, because of the combinatorial explosion of sequences of functions that could be executed on smart contracts. In this paper, we present a technique to analyze depth-n vulnerabilities in an efficient and scalable way by combining symbolic execution and data dependency analysis. A significant advantage of combining symbolic with static analysis is that it scales much better than symbolic alone and does not have the problem of false positive that static analysis tools typically have. We have implemented our technique in a tool called MPro, a scalable and automated smart contract analyzer based on the existing symbolic analysis tool Mythril-Classic and the static analysis tool Slither. We analyzed 100 randomly chosen smart contracts on MPro and our evaluation shows that MPro is about n-times faster than Mythril-Classic for detecting depth-n vulnerabilities, while preserving all the detection capabilities of Mythril-Classic.


title: Cylindrical Shape Decomposition for 3D Segmentation of Tubular Objects

abstract: We develop a cylindrical shape decomposition (CSD) algorithm to decompose an object, a union of several tubular structures, into its semantic components. We decompose the object using its curve skeleton and restricted translational sweeps. For that, CSD partitions the curve skeleton into maximal-length sub-skeletons over an orientation cost, each sub-skeleton corresponds to a semantic component. To find the intersection of the tubular components, CSD translationally sweeps the object in decomposition intervals to identify critical points at which the shape of the object changes substantially. CSD cuts the object at critical points and assigns the same label to parts along the same sub-skeleton, thereby constructing a semantic component. The proposed method further reconstructs the acquired semantic components at the intersection of object parts using generalized cylinders. We apply CSD for segmenting axons in large 3D electron microscopy images and decomposing vascular networks and synthetic objects. We show that our proposal is robust to severe surface noise and outperforms state-of-the-art decomposition techniques in its applications.


title: Probabilistic Formulation of the Take The Best Heuristic

abstract: The framework of cognitively bounded rationality treats problem solving as fundamentally rational, but emphasises that it is constrained by cognitive architecture and the task environment. This paper investigates a simple decision making heuristic, Take The Best (TTB), within that framework. We formulate TTB as a likelihood-based probabilistic model, where the decision strategy arises by probabilistic inference based on the training data and the model constraints. The strengths of the probabilistic formulation, in addition to providing a bounded rational account of the learning of the heuristic, include natural extensibility with additional cognitively plausible constraints and prior information, and the possibility to embed the heuristic as a subpart of a larger probabilistic model. We extend the model to learn cue discrimination thresholds for continuous-valued cues and experiment with using the model to account for biased preference feedback from a boundedly rational agent in a simulated interactive machine learning task.


title: An Experimental Study of a 1-planarity Testing and Embedding Algorithm

abstract: The definition of $1$-planar graphs naturally extends graph planarity, namely a graph is $1$-planar if it can be drawn in the plane with at most one crossing per edge. Unfortunately, while testing graph planarity is solvable in linear time, deciding whether a graph is $1$-planar is NP-complete, even for restricted classes of graphs. Although several polynomial-time algorithms have been described for recognizing specific subfamilies of $1$-planar graphs, no implementations of general algorithms are available to date. We investigate the feasibility of a $1$-planarity testing and embedding algorithm based on a backtracking strategy. While the experiments show that our approach can be successfully applied to graphs with up to 30 vertices, they also suggest the need of more sophisticated techniques to attack larger graphs. Our contribution provides initial indications that may stimulate further research on the design of practical approaches for the $1$-planarity testing problem.


title: Local regularity result for an optimal transportation problem with rough measures in the plane

abstract: We investigate the properties of convex functions in the plane that satisfy a local inequality which generalizes the notion of sub-solution of Monge-Ampere equation for a Monge-Kantorovich problem with quadratic cost between non-absolutely continuous measures. For each measure, we introduce a discrete scale so that the measure behaves as an absolutely continuous measure up to that scale. Our main theorem then proves that such convex functions cannot exhibit any flat part at a scale larger than the corresponding discrete scales on the measures. This, in turn, implies a $C^1$ regularity result up to the discrete scale for the Legendre transform. Our result applies in particular to any Kantorovich potential associated to an optimal transportation problem between two measures that are (possibly only locally) sums of uniformly distributed Dirac masses. The proof relies on novel explicit estimates directly based on the optimal transportation problem, instead of the Monge-Ampere equation.


title: Dynamic continuous-wave spectroscopy of coherent population trapping at phase-jump modulation

abstract: A method of dynamic continuous-wave spectroscopy of coherent population trapping (CPT) resonances using phase modulation of the jump type is developed. The time evolution of the spectroscopic signal is investigated. A method for the formation of an error signal for frequency stabilization is proposed. We show that our approach has a reduced sensitivity to the lineshape asymmetry of the CPT resonance. The experimental results are in good qualitative agreement with theoretical predictions based on a mathematical model of a three-level $\Lambda$ system in a bichromatic field. This method can be used in atomic frequency standards (including chip-scale atomic clocks).


title: Development of the Preparation Session Observational Tool (PSOT) for Quantitative Analysis by Probing the Underlying Behaviors that Contribute to Various Learning Assistant-Faculty Relationships

abstract: In order to have a successful Learning Assistant Program it is important to create adequate support tools and resources to adopt and implement the model. The Learning Assistant Model was created at the University of Colorado-Boulder and consists of three essential elements (pedagogy, content, and practice). Our research focuses on the content component, or the weekly content preparation session, in which faculty and their Learning Assistants (LAs) meet to discuss course objectives, content, pedagogical approaches, and student learning. For this study, video recordings of weekly content preparation sessions of various LA-faculty assignments at Chicago State University were analyzed. This data and a review of previous research on observation protocols led us to develop the Preparation Session Observation Tool (PSOT). PSOT is used to develop a better understanding of how the weekly content preparation sessions are fostering different types of partnerships. It is expected that PSOT data can be used by new adoptees of the LA Model and experienced members of LA Programs to analyze their weekly content preparation session in order to improve or change LA-faculty working relationships.


title: Comment on "Quantum supremacy using a programmable superconducting processor"

abstract: Even if Google AI's Sycamore processor is efficient for the particular task it has been designed for it fails to deliver universal computational capacity. Furthermore, even classical devices implementing transverse homoclinic orbits realize exponential speedups with respect to universal classical as well as quantum computations. Moreover, relative to the validity of quantum mechanics, there already exist quantum oracles which violate the Church-Turing thesis.


title: 2-D Covariant Affine Integral Quantization(s)

abstract: Covariant affine integral quantization is studied and applied to the motion of a particle in a punctured plane R^2_\ast=R^2\{0}, for which the phase space is R^2_\ast=R^2\{0}X R^2. We examine the consequences of different quantizer operators built from weight functions on this phase space. To illustrate the procedure, we examine two examples of weights. The first one corresponds to 2-D coherent state families, while the second one corresponds to the affine inversion in the punctured plane. The later yields the usual canonical quantization and a quasi-probability distribution (2-D affine Wigner function) which is real, marginal in both position and momentum.


title: Some majorization inequalities induced by Schur products in Euclidean Jordan algebras

abstract: In an Euclidean Jordan algebra V of rank n, an element x is said to be majorized by an element y, if the corresponding eigenvalue vector of x is majorized by the eigenvalue vector of y in R^n. In this article, we describe pointwise majorization inequalities of the form `T(x) majorized by S(x)', where T and S are linear transformations induced by Schur products. Specializing, we recover analogs of majorization inequalities of Schur, Hadamard, and Oppenheimer stated in the setting of Euclidean Jordan algebras, as well as majorization inequalities connecting quadratic and Lyapunov transformations on V. We also show how Schur products induced by certain scalar means (such as arithmetic, geometric, harmonic, and logarithmic means) naturally lead to majorization inequalities.


title: Introduction to Univalent Foundations of Mathematics with Agda

abstract: We introduce Voevodsky's univalent foundations and univalent mathematics, and explain how to develop them with the computer system Agda, which is based on Martin-L\"of type theory. Agda allows us to write mathematical definitions, constructions, theorems and proofs, for example in number theory, analysis, group theory, topology, category theory or programming language theory, checking them for logical and mathematical correctness. Agda is a constructive mathematical system by default, which amounts to saying that it can also be considered as a programming language for manipulating mathematical objects. But we can assume the axiom of choice or the principle of excluded middle for pieces of mathematics that require them, at the cost of losing the implicit programming-language character of the system. For a fully constructive development of univalent mathematics in Agda, we would need to use its new cubical flavour, and we hope these notes provide a base for researchers interested in learning cubical type theory and cubical Agda as the next step. Compared to most expositions of the subject, we work with explicit universe levels.


title: SDSS-IV MaNGA: Variations in the N/O -- O/H relation bias metallicity gradient measurements

abstract: In this paper we use strong line calibrations of N/O and O/H in MaNGA spaxel data to explore the systematics introduced by variations in N/O on various strong-line metallicity diagnostics. We find radial variations in N/O at fixed O/H which correlate with total galaxy stellar-mass; and which can induce $\sim 40 \%$ systematic uncertainties in oxygen abundance gradients when nitrogen-dependent abundance calibrations are used. Empirically, we find that these differences are associated with variation in the local star formation efficiency, as predicted by recent chemical evolution models for galaxies, but we cannot rule out other processes such as radial migration and the accretion of passive dwarf galaxies also playing a role.


title: VoteNet+ : An Improved Deep Learning Label Fusion Method for Multi-atlas Segmentation

abstract: In this work, we improve the performance of multi-atlas segmentation (MAS) by integrating the recently proposed VoteNet model with the joint label fusion (JLF) approach. Specifically, we first illustrate that using a deep convolutional neural network to predict atlas probabilities can better distinguish correct atlas labels from incorrect ones than relying on image intensity difference as is typical in JLF. Motivated by this finding, we propose VoteNet+, an improved deep network to locally predict the probability of an atlas label to differs from the label of the target image. Furthermore, we show that JLF is more suitable for the VoteNet framework as a label fusion method than plurality voting. Lastly, we use Platt scaling to calibrate the probabilities of our new model. Results on LPBA40 3D MR brain images show that our proposed method can achieve better performance than VoteNet.


title: Program Sketching with Live Bidirectional Evaluation

abstract: We present a system called Smyth for program sketching in a typed functional language whereby the concrete evaluation of ordinary assertions gives rise to input-output examples, which are then used to guide the search to complete the holes. The key innovation, called live bidirectional evaluation, propagates examples "backward" through partially evaluated sketches. Live bidirectional evaluation enables Smyth to (a) synthesize recursive functions without trace-complete sets of examples and (b) specify and solve interdependent synthesis goals. Eliminating the trace-completeness requirement resolves a significant limitation faced by prior synthesis techniques when given partial specifications in the form of input-output examples. To assess the practical implications of our techniques, we ran several experiments on benchmarks used to evaluate Myth, a state-of-the-art example-based synthesis tool. First, given expert examples (and no partial implementations), we find that Smyth requires on average 66% of the number of expert examples required by Myth. Second, we find that Smyth is robust to randomly-generated examples, synthesizing many tasks with relatively few more random examples than those provided by an expert. Third, we create a suite of small sketching tasks by systematically employing a simple sketching strategy to the Myth benchmarks; we find that user-provided sketches in Smyth often further reduce the total specification burden (i.e. the combination of partial implementations and examples). Lastly, we find that Leon and Synquid, two state-of-the-art logic-based synthesis tools, fail to complete several tasks on which Smyth succeeds.


title: A Perceived Environment Design using a Multi-Modal Variational Autoencoder for learning Active-Sensing

abstract: This contribution comprises the interplay between a multi-modal variational autoencoder and an environment to a perceived environment, on which an agent can act. Furthermore, we conclude our work with a comparison to curiosity-driven learning.


title: A limit theorem for the $1$st Betti number of layer-$1$ subgraphs in random graphs

abstract: We initiate the study of local topology of random graphs. The high level goal is to characterize local "motifs" in graphs. In this paper, we consider what we call the layer-$r$ subgraphs for an input graph $G = (V,E)$: Specifically, the layer-$r$ subgraph at vertex $u \in V$, denoted by $G_{u; r}$, is the induced subgraph of $G$ over vertex set $\Delta_{u}^{r}:= \left\{v \in V: d_G(u,v) = r \right\}$, where $d_G$ is shortest-path distance in $G$. Viewing a graph as a 1-dimensional simplicial complex, we then aim to study the $1$st Betti number of such subgraphs. Our main result is that the $1$st Betti number of layer-$1$ subgraphs in Erd\H{o}s--R\'enyi random graphs $G(n,p)$ satisfies a central limit theorem.


title: CNF Encodings of Cardinality Constraints Based on Comparator Networks

abstract: Boolean Satisfiability Problem (SAT) is one of the core problems in computer science. As one of the fundamental NP-complete problems, it can be used - by known reductions - to represent instances of variety of hard decision problems. Additionally, those representations can be passed to a program for finding satisfying assignments to Boolean formulas, for example, to a program called MiniSat. Those programs (called SAT-solvers) have been intensively developed for many years and - despite their worst-case exponential time complexity - are able to solve a multitude of hard practical instances. A drawback of this approach is that clauses are neither expressive, nor compact, and using them to describe decision problems can pose a big challenge on its own. We can improve this by using high-level constraints as a bridge between a problem at hand and SAT. Such constraints are then automatically translated to eqisatisfiable Boolean formulas. The main theme of this thesis revolves around one type of such constraints, namely Boolean Cardinality Constraints (or simply cardinality constraints). Cardinality constraints state that at most (at least, or exactly) k out of n propositional literals can be true. Such cardinality constraints appear naturally in formulations of different real-world problems including cumulative scheduling, timetabling or formal hardware verification. The goal of this thesis is to propose and analyze new and efficient methods to encode (translate) cardinality constraints into equisatisfiable proposition formulas in CNF, such that the resulting SAT instances are small and that the SAT-solver runtime is as short as possible.


title: Functional transcendence for the unipotent Albanese map

abstract: We prove a certain transcendence property of the unipotent Albanese map of a smooth variety, conditional on the Ax-Schanuel conjecture for variations of mixed Hodge structure. We show that this property allows the Chabauty-Kim method to be generalized to higher-dimensional varieties. In particular, we conditionally generalize several of the main Diophantine finiteness results in Chabauty-Kim theory to arbitrary number fields.


title: Identifying Dehn Functions of Bestvina--Brady Groups From Their Defining Graphs

abstract: Let $\Gamma$ be a finite simplicial graph such that the flag complex on $\Gamma$ is a $2$-dimensional triangulated disk. We show that with some assumptions, the Dehn function of the associated Bestvina--Brady group is either quadratic, cubic, or quartic. Furthermore, we can identify the Dehn function from the defining graph $\Gamma$.


title: Line Operators of Gauge Theories on Non-Spin Manifolds

abstract: We study four-dimensional gauge theories on oriented and non-spin spacetime manifolds. On such manifolds, each line operator arises only either as a boson or a fermion. Based on physical arguments, a method of systematically assigning spin labels to line operators is proposed, and several consistency checks are performed. This is used to classify all possible sets of allowed line operators -- including their spins -- for gauge theories with simple Lie algebras. The Lagrangian descriptions of the theories with these sets of allowed line operators are given. Finally, the one-form symmetries of these theories are studied by coupling to background gauge fields, and their 't Hooft anomalies are computed.


title: Graph inverse semigroups and Leavitt path algebras

abstract: We study two classes of inverse semigroups built from directed graphs, namely graph inverse semigroups and a new class of semigroups that we refer to as Leavitt inverse semigroups. These semigroups are closely related to graph $C^*$-algebras and Leavitt path algebras. We provide a topological characterization of the universal groups of the local submonoids of these inverse semigroups. We study the relationship between the graph inverse semigroups of two graphs when there is a directed immersion between the graphs. We describe the structure of graphs that admit a directed cover or directed immersion into a circle and we provide structural information about graph inverse semigroups of finite graphs that admit a directed cover onto a bouquet of circles. We also find necessary and sufficient conditions for a homomorphic image of a graph inverse semigroup to be another graph inverse semigroup. We find a presentation for the Leavitt inverse semigroup of a graph in terms of generators and relations. We describe the structure of the Leavitt inverse semigroup and the Leavitt path algebra of a graph that admits a directed immersion into a circle. We show that two graphs that have isomorphic Leavitt inverse semigroups have isomorphic Leavitt path algebras and we classify graphs that have isomorphic Leavitt inverse semigroups. As a consequence, we show that Leavitt path algebras are $0$-retracts of certain matrix algebras.


title: Non-standard anomalous heat conduction in harmonic chains with correlated isotopic disorder

abstract: We address the general problem of heat conduction in one dimensional harmonic chain, with correlated isotopic disorder, attached at its ends to white noise or oscillator heat baths. When the low wavelength $\mu$ behavior of the power spectrum $W$ (of the fluctuations of the random masses around their common mean value) scales as $W(\mu)\sim \mu^\beta$, the asymptotic thermal conductivity $\kappa$ scales with the system size $N$ as $\kappa \sim N^{(1+\beta)/(2+\beta)}$ for free boundary conditions, whereas for fixed boundary conditions $\kappa \sim N^{(\beta-1)/(2+\beta)}$; where $\beta>-1$, which is the usual power law scaling for one dimensional systems. Nevertheless, if $W$ does not scale as a power law in the low wavelength limit, the thermal conductivity may not scale in its usual form $\kappa\sim N^{\alpha}$, where the value of $\alpha$ depends on the particular one dimensional model. As an example of the latter statement, if $W(\mu)\sim \exp(-1/\mu)/\mu^2$, $\kappa \sim N/(\log N)^3$ for fixed boundary conditions and $\kappa \sim N/\log(N)$ for free boundary conditions, which represent non-standard scalings of the thermal conductivity.


title: Heat conduction in harmonic chains with Levy-type disorder

abstract: We consider heat transport in one-dimensional harmonic chains attached at its ends to Langevin heat baths. The harmonic chain has mass impurities where the separation $d$ between any two successive impurities is randomly distributed according to a power-law distribution $P(d)\sim 1/d^{\alpha+1}$, being $\alpha>0$. In the regime where the first moment of the distribution is well defined ($1<\alpha<2$) the thermal conductivity $\kappa$ scales with the system size $N$ as $\kappa\sim N^{(\alpha-3)/\alpha}$ for fixed boundary conditions, whereas for free boundary conditions $\kappa\sim N^{(\alpha-1)/\alpha}$ if $N\gg1$. When $\alpha=2$, the inverse localization length $\lambda$ scales with the frequency $\omega$ as $\lambda\sim \omega^2 \ln \omega$ in the low frequency regime, due to the logarithmic correction, the size scaling law of the thermal conductivity acquires a non-closed form. When $\alpha>2$, the thermal conductivity scales as in the uncorrelated disorder case. The situation $\alpha<1$ is only analyzed numerically, where $\lambda(\omega)\sim \omega^{2-\alpha}$ which leads to the following asymptotic thermal conductivity: $\kappa \sim N^{-(\alpha+1)/(2-\alpha)}$ for fixed boundary conditions and $\kappa \sim N^{(1-\alpha)/(2-\alpha)}$ for free boundary conditions.


title: Entropy-stable positivity-preserving DG schemes for Boltzmann-Poisson models of collisional electronic transport along energy bands

abstract: This work is related to developing entropy-stable positivity-preserving Discontinuous Galerkin (DG) methods as a computational scheme for Boltzmann-Poisson systems modeling the probability density of collisional electronic transport along semiconductor energy bands. In momentum coordinates representing spherical / energy-angular variables, we pose the respective Vlasov-Boltzmann equation with a linear collision operator and a singular measure, modeling scatterings as functions of the band structure appropriately for hot electron nanoscale transport. We show stability results of semi-discrete DG schemes under an entropy norm for 1D-position (2D-momentum) and 2D-position (3D-momentum), using dissipative properties of the collisional operator given its entropy inequality. The latter depends on an exponential of the Hamiltonian rather than the Maxwellian associated with only kinetic energy. For the 1D problem, knowing the analytic solution to the Poisson equation and convergence to a constant current is crucial to obtaining full stability (weighted entropy norm decreasing over time). For the 2D problem, specular reflection boundary conditions and periodicity are considered in estimating stability under an entropy norm. Regarding the positivity-preservation proofs in the DG scheme for the 1D problem, inspired by \cite{ZhangShu1}, \cite{ZhangShu2}, and \cite{CGP}, \cite{EECHXM-JCP}, we treat collisions as a source and find convex combinations of the transport and collision terms which guarantee positivity of the cell average of our numerical probability density at the next time. The positivity of the numerical solution to the probability density in the domain is guaranteed by applying the limiters in \cite{ZhangShu1} and \cite{ZhangShu2} that preserve the cell average modifying the slope of the piecewise linear solutions to make the function non-negative.