title: GLRT based Adaptive-Thresholding for CFAR-Detection of Pareto-Target in Pareto-Distributed Clutter

abstract: After Pareto distribution has been validated for sea clutter returns in varied scenarios, some heuristics of adaptive-thresholding appeared in the literature for constant false alarm rate (CFAR) criteria. These schemes used the same adaptive-thresholding form that was originally derived for detecting Swerling-I (exponential) target in exponentially distributed clutter. Statistical procedures obtained under such idealistic assumptions would affect the detection performance when applied to newer target and clutter models, esp. heavy tail distributions like Pareto. Further, in addition to the sea clutter returns, it has also been reported that Generalized Pareto distribution fits best for the measured Radar-cross-section (RCS) data of a SAAB aircraft. Therefore, in Radar application scenarios like Airborne Warning and Control System (AWACS), when both the target and clutter are Pareto distributed, we pose the detection problem as a two-sample, Pareto vs. Pareto composite hypothesis testing problem. We address this problem by corroborating the binary hypothesis framework instead of the conventional way of tweaking the existing adaptive-thresholding CFAR detector. Whereby, for the composite case, considering no knowledge of both scale and shape parameters of Pareto distributed clutter, we derive the new adaptive-thresholding detector based on the generalized likelihood ratio test (GLRT) statistic. We further show that our proposed adaptive-thresholding detector has a CFAR property. We provide extensive simulation results to demonstrate the performance of the proposed detector.


title: On the equipartition theorem and black holes nongaussian entropies

abstract: In this Letter we have shown that, from the standard thermodynamic functions, the mathematical form of an equipartition theorem may be related to the algebraic expression of a particular entropy initially chosen to describe the black hole thermodynamics. Namely, we have different equipartition expressions for distinct statistics. To this end, four different mathematical expressions for the entropy have been selected to demonstrate our objective. Furthermore, a possible phase transition is observed in the heat capacity behavior of the Tsallis and Cirto entropy model.


title: Evolving basal slip under glaciers and ice streams

abstract: Many fast-flowing glaciers and ice streams move over beds consisting of reworked sediments and erosional products, commonly referred to as till. The complex interplay between ice, meltwater, and till at the subglacial bed connects several fundamental problems in glaciology including the debate about rapid mass loss from the ice sheets, the formation and evolution of subglacial landforms, and the storage and transport of subglacial water. In-situ measurements have probed the subglacial bed, but provide surprisingly variable and seemingly inconsistent evidence of the depth where deformation occurs, even at a given field site. These observations suggest that subglacial beds are inherently dynamic. The goal of this paper is to advance our understanding of the physical processes that contribute to the dynamics of subglacial beds as reflected in existing observations of basal deformation. We build on recent advances in modeling dense, granular flows to derive a new numerical model for water-saturated till. Our model demonstrates that changes in the force balance or temporal variations in water pressure can shift slip away from the ice-bed interface and far into the bed, causing episodes of significantly enhanced till transport because the entire till layer above the deep slip interface becomes mobilized. We compare our model results against observations of basal deformation from both mountain glacier and ice-stream settings in the past and present. We also present an analytical solution to assess the variability of till transport for different glacial settings and hydraulic properties.


title: Constraints and correlations of nuclear matter parameters from a density-dependent van der Waals model

abstract: A recently proposed density-dependent van der Waals model, with only $4$ free parameters adjusted to fix binding energy, saturation density, symmetry energy, and incompressibility, is analyzed under symmetric and asymmetric nuclear matter constraints. In a previous paper, it was shown that this model is fully consistent with the constraints related to the binary neutron star merger event named GW170817 and reported by the LIGO and Virgo collaboration. Here, we show that it also describes satisfactorily the low and high-density regions of symmetric nuclear matter, with all the main constraints satisfied. We also found a linear correlation between the incompressibility and the skewness parameter, both at the saturation density and show how it relates to the crossing point presented in the incompressibility as a function of the density. In the asymmetric matter regime, other linear correlations are found, namely, the one between the symmetry energy ($J$) and its slope ($L_0$), and other one establishing the symmetry energy curvature as a function of the combination given by $3J-L_0$.


title: Fluctuation Around the Circular Law for Random Matrices with Real Entries

abstract: We extend our recent result [Cipolloni, Erd\H{o}s, Schr\"oder 2019] on the central limit theorem for the linear eigenvalue statistics of non-Hermitian matrices $X$ with independent, identically distributed complex entries to the real symmetry class. We find that the expectation and variance substantially differ from their complex counterparts, reflecting (i) the special spectral symmetry of real matrices onto the real axis; and (ii) the fact that real i.i.d. matrices have many real eigenvalues. Our result generalizes the previously known special cases where either the test function is analytic [O'Rourke, Renfrew 2016] or the first four moments of the matrix elements match the real Gaussian [Tao, Vu 2015; Kopel 2015]. The key element of the proof is the analysis of several weakly dependent Dyson Brownian motions (DBMs). The conceptual novelty of the real case compared with [Cipolloni, Erd\H{o}s, Schr\"oder 2019] is that the correlation structure of the stochastic differentials in each individual DBM is non-trivial, potentially even jeopardising its well-posedness.


title: Exploring Hubble Constant Data in an Introductory Course

abstract: Popular accounts of exciting discoveries often draw students to physics and astronomy, but at the introductory level it is challenging to connect with these in a meaningful way. The use of real astronomical data in the classroom can help bridge this gap and build valuable quantitative and scientific reasoning skills. This paper presents a strategy for studying Hubble's Law and the accelerating expansion of the universe using actual data. Along with understanding the physical concepts, an explicit goal is to develop skills for analyzing data in terms of a model.


title: A locality-based approach for coded computation

abstract: Modern distributed computation infrastructures are often plagued by unavailabilities such as failing or slow servers. These unavailabilities adversely affect the tail latency of computation in distributed infrastructures. The simple solution of replicating computation entails significant resource overhead. Coded computation has emerged as a resource-efficient alternative, wherein multiple units of data are encoded to create parity units and the function to be computed is applied to each of these units on distinct servers. A decoder can use the available function outputs to decode the unavailable ones. Existing coded computation approaches are resource efficient only for simple variants of linear functions such as multilinear, with even the class of low degree polynomials requiring the same multiplicative overhead as replication for practically relevant straggler tolerance. In this paper, we present a new approach to model coded computation via the lens of locality of codes. We introduce a generalized notion of locality, denoted computational locality, building upon the locality of an appropriately defined code. We show that computational locality is equivalent to the required number of workers for coded computation and leverage results from the well-studied locality of codes to design coded computation schemes. We show that recent results on coded computation of multivariate polynomials can be derived using local recovering schemes for Reed-Muller codes. We present coded computation schemes for multivariate polynomials that adaptively exploit locality properties of input data-- an inadmissible technique under existing frameworks. These schemes require fewer workers than the lower bound under existing coded computation frameworks, showing that the existing multiplicative overhead on the number of servers is not fundamental for coded computation of nonlinear functions.


title: Poincar\'e compactification for non-polynomial vector fields

abstract: In this work a theorical framework to apply the Poincar\'e compactification technique to locally Lipschitz continuous vector fields is developed. It is proved that these vectors fields are compactifiable in the n-dimensional sphere, though the compactified vector field can be identically null in the equator. Moreover, for a fixed projection to the hemisphere, all the compactifications of a vector field, which are not identically null on the equator are equivalent. Also, the conditions determining the invariance of the equator for the compactified vector field are obtained. Up to the knowledge of the authors, this is the first time that the Poincar\'e compactification of locally Lipschitz continuous vector fields is studied. These results are illustrated applying them to some families of vector fields, like polynomial vector fields, vector fields defined as a sum of homogeneous functions and vector fields defined by piecewise linear functions.


title: The Eisenstein ideal for weight k and a Bloch-Kato conjecture for tame families

abstract: We study the Eisenstein ideal for modular forms of even weight $k>2$ and prime level $N$. We pay special attention to the phenomenon of $\mathit{extra \ reducibility}$: the Eisenstein ideal is strictly larger than the ideal cutting out reducible Galois representations. We prove a modularity theorem for these extra reducible representations. As consequences, we relate the derivative of a Mazur-Tate $L$-function to the rank of the Hecke algebra, generalizing a theorem of Merel, and give a new proof of a special case of an equivariant main conjecture of Kato. In the second half of the paper, we recall Kato's formulation of this main conjecture in the case of a family of motives given by twists by characters of conductor $N$ and $p$-power order and its relation to other formulations of the equivariant main conjecture.


title: A homotopy approach to coherent quantum LQG control synthesis using discounted performance criteria

abstract: This paper is concerned with linear-quadratic-Gaussian (LQG) control for a field-mediated feedback connection of a plant and a coherent (measurement-free) controller. Both the plant and the controller are multimode open quantum harmonic oscillators governed by linear quantum stochastic differential equations. The control objective is to make the closed-loop system internally stable and to minimize the infinite-horizon quadratic cost involving the plant variables and the controller output subject to quantum physical realizability (PR) constraints. This coherent quantum LQG (CQLQG) control problem, which has been of active research interest for over ten years, does not admit a solution in the form of separation principle and independent Riccati equations known for its classical counterpart. We apply variational techniques to a family of discounted CQLQG control problems parameterized by an effective time horizon. This gives rise to a homotopy algorithm, which is initialized with a PR (but not necessarily stabilizing) controller and aims at a locally optimal stabilizing controller for the original problem in the limit.


title: Liquid-solid slip on charged walls: dramatic impact of charge distribution

abstract: Nanofluidic systems show great promises for applications in energy conversion, where their performance can be enhanced by nanoscale liquid-solid slip. However, efficiency is also controlled by surface charge, which is known to reduce slip. Combining molecular dynamics simulations and analytical developments, we show the dramatic impact of surface charge distribution on the slip-charge coupling. Homogeneously charged graphene exhibits a very favorable slip-charge relation (rationalized with a new theoretical model correcting some weaknesses of the existing ones), leading to giant electrokinetic energy conversion. In contrast, slip is strongly affected on heterogeneously charged surfaces, due to the viscous drag induced by counter-ions trapped on the surface. In that case slip should depend on the detailed physical chemistry of the interface controlling the fraction of bound ions. Our numerical results and theoretical models provide new fundamental insight on the molecular mechanisms of liquid-solid slip, and practical guidelines for searching new functional interfaces with optimal energy conversion properties, e.g. for blue energy or waste heat harvesting.


title: ViWi Vision-Aided mmWave Beam Tracking: Dataset, Task, and Baseline Solutions

abstract: Vision-aided wireless communication is motivated by the recent advances in deep learning and computer vision as well as the increasing dependence on line-of-sight links in millimeter wave (mmWave) and terahertz systems. By leveraging vision, this new research direction enables an interesting set of new capabilities such as vision-aided mmWave beam and blockage prediction, proactive hand-off, and resource allocation among others. These capabilities have the potential of reliably supporting highly-mobile applications such as vehicular/drone communications and wireless virtual/augmented reality in mmWave and terahertz systems. Investigating these interesting applications, however, requires the development of special dataset and machine learning tasks. Based on the Vision-Wireless (ViWi) dataset generation framework [1], this paper develops an advanced and realistic scenario/dataset that features multiple base stations, mobile users, and rich dynamics. Enabled by this dataset, the paper defines the vision-wireless mmWave beam tracking task (ViWi-BT) and proposes a baseline solution that can provide an initial benchmark for the future ViWi-BT algorithms.


title: Dislocation screening in crystals with spherical topology

abstract: Whereas disclination defects are energetically prohibitive in two-dimensional flat crystals, their existence is necessary in crystals with spherical topology, such as viral capsids, colloidosomes or fullerenes. Such a geometrical frustration gives rise to large elastic stresses, which render the crystal unstable when its size is significantly larger than the typical lattice spacing. Depending on the compliance of the crystal with respect to stretching and bending deformations, these stresses are alleviated by either a local increase of the intrinsic curvature in proximity of the disclinations or by the proliferation of excess dislocations, often organized in the form of one-dimensional chains known as "scars". The associated strain field of the scars is such to counterbalance the one resulting from the isolated disclinations. Here, we develop a continuum theory of dislocation screening in two-dimensional closed crystals with genus one. Upon modeling the flux of scars emanating from a given disclination as an independent scalar field, we demonstrate that the elastic energy of closed two-dimensional crystals with various degrees of asphericity can be expressed as a simple quadratic function of the screened topological charge of the disclinations, both at zero and finite temperature. This allows us to predict the optimal density of the excess dislocations as well as the minimal stretching energy attained by the crystal.


title: Computing the norm of nonnegative matrices and the log-Sobolev constant of Markov chains

abstract: We analyze the global convergence of the power iterates for the computation of a general mixed-subordinate matrix norm. We prove a new global convergence theorem for a class of entrywise nonnegative matrices that generalizes and improves a well-known results for mixed-subordinate $\ell^p$ matrix norms. In particular, exploiting the Birkoff--Hopf contraction ratio of nonnegative matrices, we obtain novel and explicit global convergence guarantees for a range of matrix norms whose computation has been recently proven to be NP-hard in the general case, including the case of mixed-subordinate norms induced by the vector norms made by the sum of different $\ell^p$-norms of subsets of entries. Finally, we use the new results combined with hypercontractive inequalities to prove a new lower bound on the logarithmic Sobolev constant of a Markov chain.


title: The Holographic Nature of Null Infinity

abstract: We argue that, in a theory of quantum gravity in a four dimensional asymptotically flat spacetime, all information about massless excitations can be obtained from an infinitesimal neighbourhood of the past boundary of future null infinity and does not require observations over all of future null infinity. Moreover, all information about the state that can be obtained through observations near a cut of future null infinity can also be obtained from observations near any earlier cut although the converse is not true. We provide independent arguments for these two assertions. Similar statements hold for past null infinity. These statements have immediate implications for the information paradox since they suggest that the fine-grained von Neumann entropy of the state defined on a segment $(-\infty,u)$ of future null infinity is independent of u. This is very different from the oft-discussed Page curve that this entropy is sometimes expected to obey. We contrast our results with recent discussions of the Page curve in the context of black hole evaporation, and also discuss the relation of our results to other proposals for holography in flat space.


title: The OTELO survey. II. The faint-end of the H$\alpha$ luminosity function at z $\sim$ 0.40

abstract: We take advantage of the capability of the OTELO survey to obtain the H$\alpha$ luminosity function (LF) at ${\rm z}\sim0.40$. Because of the deepest coverage of OTELO, we are able to determine the faint end of the LF, and thus better constrain the star formation rate and the number of galaxies at low luminosities. The AGN contribution to this LF is estimated as well. We make use of the multi-wavelength catalogue of objects in the field compiled by the OTELO survey, which is unique in terms of minimum flux and equivalent width. We also take advantage of the pseudo-spectra built for each source, which allow the identification of emission lines and the discrimination of different types of objects. The H$\alpha$ luminosity function at $z\sim0.40$ is obtained, which extends the current faint end by almost 1 dex, reaching minimal luminosities of $\log_{10}L_{\rm lim}=38.5$ erg s$^{-1}$ (or $\sim0.002\, \text{M}_\odot\text{ yr}^{-1})$. The AGN contribution to the total H$\alpha$ luminosity is estimated. We find that no AGN should be expected below a luminosity of $\log_{10}L=38.6$ erg s$^{-1}$. From the sample of non-AGN (presumably, pure SFG) at $z\sim0.40$ we estimated a star formation rate density of $\rho_{\rm SFR}=0.012\pm0.005\ {\rm \text{M}_{\odot}\ yr^{-1}\ Mpc^{-3}}$.


title: Goal-Oriented Multi-Task BERT-Based Dialogue State Tracker

abstract: Dialogue State Tracking (DST) is a core component of virtual assistants such as Alexa or Siri. To accomplish various tasks, these assistants need to support an increasing number of services and APIs. The Schema-Guided State Tracking track of the 8th Dialogue System Technology Challenge highlighted the DST problem for unseen services. The organizers introduced the Schema-Guided Dialogue (SGD) dataset with multi-domain conversations and released a zero-shot dialogue state tracking model. In this work, we propose a GOaL-Oriented Multi-task BERT-based dialogue state tracker (GOLOMB) inspired by architectures for reading comprehension question answering systems. The model "queries" dialogue history with descriptions of slots and services as well as possible values of slots. This allows to transfer slot values in multi-domain dialogues and have a capability to scale to unseen slot types. Our model achieves a joint goal accuracy of 53.97% on the SGD dataset, outperforming the baseline model.


title: Federated Orchestration for Network Slicing of Bandwidth and Computational Resource

abstract: Network slicing has been considered as one of the key enablers for 5G to support diversified IoT services and application scenarios. This paper studies the distributed network slicing for a massive scale IoT network supported by 5G with fog computing. Multiple services with various requirements need to be supported by both spectrum resource offered by 5G network and computational resourc of the fog computing network. We propose a novel distributed framework based on a new control plane entity, federated-orchestrator , which can coordinate the spectrum and computational resources without requiring any exchange of the local data and resource information from BSs. We propose a distributed resource allocation algorithm based on Alternating Direction Method of Multipliers with Partial Variable Splitting . We prove DistADMM-PVS minimizes the average service response time of the entire network with guaranteed worst-case performance for all supported types of services when the coordination between the F-orchestrator and BSs is perfectly synchronized. Motivated by the observation that coordination synchronization may result in high coordination delay that can be intolerable when the network is large in scale, we propose a novel asynchronized ADMM algorithm. We prove that AsynADMM can converge to the global optimal solution with improved scalability and negligible coordination delay. We evaluate the performance of our proposed framework using two-month of traffic data collected in a in-campus smart transportation system supported by a 5G network. Extensive simulation has been conducted for both pedestrian and vehicular-related services during peak and non-peak hours. Our results show that the proposed framework offers significant reduction on service response time for both supported services, especially compared to network slicing with only a single resource.


title: Thermoelectric power factor of pure and doped ZnSb via DFT based defect calculations

abstract: The power factor of pure p-type ZnSb has been calculated via ab initio simulations assuming that the carrier concentrations are due to the doping effect of intrinsic zinc vacancies. With a vacancy concentration close to the experimental solubility limit we were able to perfectly reproduce the Power Factor measured in polycrystalline ZnSb samples. The methodology has then been successfully extended for predicting the effect of extrinsic doping elements on the thermoelectric properties of ZnSb. Germanium and tin seem to be promising p-type doping elements. In addition, we give, for the first time, an explanation of why it is difficult to synthesize polycrystalline n-type ZnSb samples. Indeed, compensative effects between intrinsic defects (zinc vacancies) and doping elements (Ga, or In) explain the existence of an optimal (and relatively high) dopant concentration necessary to convert ZnSb into an n-type semiconductor.


title: Modeling Engagement in Long-Term, In-Home Socially Assistive Robot Interventions for Children with Autism Spectrum Disorders

abstract: Socially assistive robotics (SAR) has great potential to provide accessible, affordable, and personalized therapeutic interventions for children with autism spectrum disorders (ASD). However, human-robot interaction (HRI) methods are still limited in their ability to autonomously recognize and respond to behavioral cues, especially in atypical users and everyday settings. This work applies supervised machine learning algorithms to model user engagement in the context of long-term, in-home SAR interventions for children with ASD. Specifically, we present two types of engagement models for each user: (i) generalized models trained on data from different users; and (ii) individualized models trained on an early subset of the user's data. The models achieved approximately 90% accuracy (AUROC) for post hoc binary classification of engagement, despite the high variance in data observed across users, sessions, and engagement states. Moreover, temporal patterns in model predictions could be used to reliably initiate re-engagement actions at appropriate times. These results validate the feasibility and challenges of recognition and response to user disengagement in long-term, real-world HRI settings. The contributions of this work also inform the design of engaging and personalized HRI, especially for the ASD community.


title: Spatially-resolved UV diagnostics of AGN feedback: radiation pressure dominates in a prototypical quasar-driven superwind

abstract: Galactic-scale winds driven by active galactic nuclei (AGN) are often invoked to suppress star formation in galaxy evolution models, but the mechanisms driving these outflows are hotly debated. Two key AGN feedback models are (1) radiation pressure accelerating cool gas and (2) a hot outflowing wind entraining the ISM. Highly ionized emission-line diagnostics represent a powerful means of differentiating these scenarios because of their sensitivity to the expected compression of the ISM clouds by the hot wind. Here, we report the first spatially resolved UV emission spectroscopy of a prototypical (radio-quiet) quasar-driven superwind around the obscured quasar SDSSJ1356+1026 at z=0.123. We observe ratios of OVI/CIV, NV/CIV, and CIV/HeII that are remarkably similar for outflowing gas clouds <100 pc and ~10 kpc from the nucleus. Such similarity is expected for clouds with AGN radiation pressure dominated dynamics. Comparing the observed line emission to models of clouds in balance with radiation pressure and/or a hot wind, we rule out the presence of a dynamically important hot wind and constrain the ratio of hot gas pressure to radiation pressure to P_hot/P_rad<0.25 both at <100 pc and ~10 kpc from the nucleus. Moreover, the predictions of the radiation pressure confined cloud models that best fit observed UV line ratios are consistent with the observed diffuse X-ray spectrum. These results indicate that this AGN superwind is driven by radiation pressure or was driven by a hot wind that has since dissipated despite on-going AGN activity.


title: Planetary evolution with atmospheric photoevaporation I. Analytical derivation and numerical study of the evaporation valley and transition from super-Earths to sub-Neptunes

abstract: Observations have revealed in the Kepler data a depleted region separating smaller super-Earths from larger sub-Neptunes. This can be explained as an evaporation valley between planets with and without H/He that is caused by atmospheric escape. First, we conduct numerical simulations of the evolution of close-in low-mass planets with H/He undergoing escape. Second, we develop an analytical model for the valley locus. We find that the bottom of the valley quantified by the radius of the largest stripped core $R_{\rm b}$ at a given orbital distance depends only weakly on post-formation H/He mass. The reason is that a high initial H/He mass means that there is more gas to evaporate, but also that the planet density is lower, increasing loss. Regarding stellar $L_{\rm XUV}$, $R_{\rm b}$ scales as $L_{\rm XUV}^{0.135}$. The same weak dependency applies to the efficiency factor $\varepsilon$ of energy-limited evaporation. As found numerically and analytically, $R_{\rm b}$ varies as function of orbital period $P$ for a constant $\varepsilon$ as $P^{-2 p_{\rm c}/3}\approx P^{-0.18}$ where $M \propto R^{p_{\rm c}}$ is the mass-radius relation of solid cores. $R_{\rm b}$ is about 1.7 $R_{\oplus}$ at a 10-day orbit for an Earth-like composition, increasing linearly with ice mass fraction. The numerical results are explained very well with the analytical model where complete evaporation occurs if the temporal integral over the stellar XUV irradiation absorbed by the planet is larger than binding energy of the envelope in the gravitational potential of the core. The weak dependency on primordial H/He mass, $L_{\rm XUV}$ and $\varepsilon$ explains why observationally the valley is visible, and why theoretically models find similar results. At the same time, given the large observed spread of $L_{\rm XUV}$, the dependency on it is still strong enough to explain why the valley is not completely empty.


title: Towards a Swampland Global Symmetry Conjecture using Weak Gravity

abstract: It is widely believed and in part established that exact global symmetries are inconsistent with quantum gravity. One then expects that approximate global symmetries can be quantitatively constrained by quantum gravity or swampland arguments. We provide such a bound for an important class of global symmetries: Those arising from a gauged $U(1)$ with the vector made massive via Higgsing with an axion. The latter necessarily couples to instantons, and their action can be constrained, using both the electric and magnetic version of the axionic weak gravity conjecture, in terms of the cutoff of the theory. As a result, instanton-induced symmetry breaking operators with a suppression factor not smaller than $\exp(-M_{\rm P}^2/\Lambda^2)$ are present, where $\Lambda$ is a cutoff of the 4d effective theory. We provide a general argument and clarify the meaning of $\Lambda$. Simple 4d and 5d models are presented to illustrate this, and we recall that this is the standard way in which things work out in string compactifications with brane instantons. The relation of our constraint to bounds that can be derived from wormholes or gravitational instantons and to those motivated by black-hole effects at finite temperature are discussed, and we present a generalization of the Giddings-Strominger wormhole solution to the case of a gauge-derived $U(1)$ global symmetry. Finally, we discuss potential loopholes to our arguments.


title: Detection of PAH Absorption and Determination of the Mid-Infrared Diffuse Interstellar Extinction Curve from the Sightline Toward Cyg OB2-12

abstract: The sightline toward the luminous blue hypergiant Cyg OB2-12 is widely used in studying interstellar dust on account of its large extinction ($A_V \simeq 10$ mag) and the fact that this extinction appears to be dominated by dust typical of the diffuse interstellar medium. We present a new analysis of archival ISO-SWS and Spitzer IRS observations of Cyg OB2-12 using a model of the emission from the star and its stellar wind to determine the total extinction $A_\lambda$ from 2.4--37 $\mu$m. In addition to the prominent 9.7 and 18 $\mu$m silicate features, we robustly detect absorption features associated with polycyclic aromatic hydrocarbons (PAHs), including the first identification of the 7.7 $\mu$m feature in absorption. The 3.3 $\mu$m aromatic feature is found to be much broader in absorption than is typically seen in emission. The 3.4 and 6.85 $\mu$m aliphatic hydrocarbon features are observed with relative strengths consistent with observation of these features on sightlines toward the Galactic Center. We identify and characterize more than sixty spectral lines in this wavelength range, which may be useful in constraining models of the star and its stellar wind. Based on this analysis, we present an extinction curve $A_\lambda/A_{2.2 \mu m}$ that extrapolates smoothly to determinations of the mean Galactic extinction curve at shorter wavelengths and to dust opacities inferred from emission at longer wavelengths, providing a new constraint on models of interstellar dust in the mid-infrared.


title: General Quantum Resource Theories: Distillation, Formation and Consistent Resource Measures

abstract: Quantum resource theories (QRTs) provide a unified theoretical framework for understanding inherent quantum-mechanical properties that serve as resources in quantum information processing, but resources motivated by physics may possess intractable mathematical structure to analyze, such as non-uniqueness of maximally resourceful states, lack of convexity, and infinite dimension. We investigate state conversion and resource measures in general QRTs under minimal assumptions to figure out universal properties of physically motivated quantum resources that may have such intractable mathematical structure. In the general setting, we prove the existence of maximally resourceful states in one-shot state conversion. Also analyzing asymptotic state conversion, we discover catalytic replication of quantum resources, where a resource state is infinitely replicable by free operations. In QRTs without assuming uniqueness of maximally resourceful states, we formulate the tasks of distillation and formation of quantum resources, and introduce distillable resource and resource cost based on the distillation and the formation, respectively. Furthermore, we introduce consistent resource measures that quantify the amount of quantum resources without contradicting the rate of state conversion even in QRTs with non-unique maximally resourceful states. Progressing beyond the previous work showing a uniqueness theorem for additive resource measures, we prove the corresponding uniqueness inequality for the consistent resource measures; that is, consistent resource measures of a quantum state take values between the distillable resource and the resource cost of the state. These formulations and results establish a foundation of QRTs applicable to mathematically intractable but physically motivated quantum resources in a unified way.


title: Universality in the classical limit of massless gravitational scattering

abstract: We demonstrate the universality of the gravitational classical deflection angle of massless particles through O(G^3) by studying the high-energy limit of full two-loop four-graviton scattering amplitudes in pure Einstein gravity as well as N >= 4 supergravity. As a by-product, our first-principles calculation provides a direct confirmation of the massless deflection angle in Einstein gravity proposed long ago by Amati, Ciafaloni and Veneziano, and is inconsistent with a recently proposed alternative.


title: Intelligent Arxiv: Sort daily papers by learning users topics preference

abstract: Current daily paper releases are becoming increasingly large and areas of research are growing in diversity. This makes it harder for scientists to keep up to date with current state of the art and identify relevant work within their lines of interest. The goal of this article is to address this problem using Machine Learning techniques. We model a scientific paper to be built as a combination of different scientific knowledge from diverse topics into a new problem. In light of this, we implement the unsupervised Machine Learning technique of Latent Dirichlet Allocation (LDA) on the corpus of papers in a given field to: i) define and extract underlying topics in the corpus; ii) get the topics weight vector for each paper in the corpus; and iii) get the topics weight vector for new papers. By registering papers preferred by a user, we build a user vector of weights using the information of the vectors of the selected papers. Hence, by performing an inner product between the user vector and each paper in the daily Arxiv release, we can sort the papers according to the user preference on the underlying topics. We have created the website IArxiv.org where users can read sorted daily Arxiv releases (and more) while the algorithm learns each users preference, yielding a more accurate sorting every day. Current IArxiv.org version runs on Arxiv categories astro-ph, gr-qc, hep-ph and hep-th and we plan to extend to others. We propose several new useful and relevant implementations to be additionally developed as well as new Machine Learning techniques beyond LDA to further improve the accuracy of this new tool.

categories: cs.LG astro-ph.HE gr-qc hep-ph hep-th stat.ML


title: Convection with Misaligned Gravity and Rotation: Simulations and Rotating Mixing Length Theory

abstract: We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.


title: Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

abstract: We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.


title: Relaxing the Higgs mass and its vacuum energy by living at the top of the potential

abstract: We consider an ultra-light scalar coupled to the Higgs in the presence of heavier new physics. In the electroweak broken phase the Higgs gives a tree-level contribution to the light-scalar potential, while new physics contributes at loop level. Thereby, the theory has a cosmologically meta-stable phase where the light scalar is around the top of its potential, and the Higgs is a loop factor lighter than new physics. Such regions with precarious naturalness are anthropically and environmentally selected, as regions with heavier Higgs crunch quickly. We expect observable effects of rolling in the dark-energy equation of state. Furthermore, vacuum energies up to the weak scale can be canceled down to anthropically small values.


title: Cosmology with quasars: predictions for eROSITA from a quasar Hubble diagram

abstract: Our group has developed a technique that makes use of the observed non-linear relation between the ultraviolet and the X-ray luminosity in quasars to provide an independent measurement of their distances, thus turning quasars into standardizable candles. This technique, at present, it is mostly based upon quasar samples with data from public catalogues both in the X-rays and in the optical/ultraviolet and extends the Hubble diagram of supernovae to a redshift range still poorly explored (z>2). From the X-ray perspective, we are now on the eve of a major change, as the upcoming mission eROSITA is going to provide us with up to ~3 millions of active galactic nuclei across the entire sky. Here we present predictions for constraining cosmological parameters such as the amount of dark matter ($\Omega_{\rm m}$), dark energy ($\Omega_\Lambda$) and the evolution of the equation of state of dark energy ($w$) through the Hubble diagram of quasars, based on the 4-year eROSITA all-sky survey. Our simulations show that the eROSITA quasars, complemented by redshift and broad-band photometric information, will supply the largest quasar sample at z<2, but with very few objects available for cosmology at higher redshifts that survives the cut for the Malmquist bias, as eROSITA will sample the brighter end of the X-ray luminosity function. The power of the quasar Hubble diagram for precision cosmology lies in the high-redshift regime, where quasars can be observed up to redshift ~7.5, essential to discriminate amongst different model extrapolations. Therefore, to be competitive for cosmology, the eROSITA quasar Hubble diagram must be complemented with the already available quasar samples and dedicated (deep) large programmes at redshift z>3.


title: Quantum vs classical dynamics in a spin-boson system: manifestations of spectral correlations and scarring

abstract: We compare the entire classical and quantum evolutions of the Dicke model in its regular and chaotic domains. This is a paradigmatic interacting spin-boson model of great experimental interest. By studying the classical and quantum survival probabilities of initial coherent states, we identify features of the long-time dynamics that are purely quantum and discuss their impact on the equilibration times. We show that the ratio between the quantum and classical asymptotic values of the survival probability serves as a metric to determine the proximity to a separatrix in the regular regime and to distinguish between two manifestations of quantum chaos: scarring and ergodicity. In the case of maximal quantum ergodicity, our results are analytical and show that quantum equilibration takes longer than classical equilibration.


title: Cross-correlating galaxy catalogs and gravitational waves: a tomographic approach

abstract: Unveiling the origin of the coalescing binaries detected via gravitational waves (GW) is challenging, notably if no multi-wavelength counterpart is detected. One important diagnostic tool is the coalescing binary distribution with respect to the large scale structures (LSS) of the universe, which we quantify via the cross-correlation of galaxy catalogs with GW ones. By using both existing and forthcoming galaxy catalogs and using realistic Monte Carlo simulations of GW events, we find that the cross-correlation signal should be marginally detectable in a 10-year data taking of advanced LIGO-Virgo detectors at design sensitivity, at least for binary neutron star mergers. The expected addition of KAGRA and LIGO-India to the GW detector network would allow for a firmer detection of this signal and, in combination with future cosmological surveys, would also permit the detection of cross-correlation for coalescing black holes. Such a measurement may unveil, for instance, a primordial origin of coalescing black holes. To attain this goal, we find that it is crucial to adopt a tomographic approach and to reach a sufficiently accurate localization of GW events. The depth of forthcoming surveys will be fully exploited by third generation GW detectors such as the Einstein Telescope or the Cosmic Explorer, which will allow one to perform precision studies of the coalescing black hole LSS distribution and attain rather advanced model discrimination capabilities.


title: Hydrodynamic Response of the Intergalactic Medium to Reionization

abstract: The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ $M_\odot$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $\Delta t = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $\Delta t = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.


title: Long-term NIR Variability in the UKIDSS Ultra Deep Survey: a new probe of AGN activity at high redshift

abstract: We present the first attempt to select AGN using long-term NIR variability. By analysing the K-band light curves of all the galaxies in the UKIDSS Ultra Deep Survey, the deepest NIR survey over ~1 sq degree, we have isolated 393 variable AGN candidates. A comparison to other selection techniques shows that only half of the variable sources are also selected using either deep Chandra X-ray imaging or IRAC colour selection, suggesting that using NIR variability can locate AGN that are missed by more standard selection techniques. In particular, we find that long-term NIR variability identifies AGN at low luminosities and in host galaxies with low stellar masses, many of which appear relatively X-ray quiet.


title: Prevalence of Complex Organic Molecules in Starless and Prestellar Cores within the Taurus Molecular Cloud

abstract: The detection of complex organic molecules (COMs) toward dense, collapsing prestellar cores has sparked interest in the fields of astrochemistry and astrobiology, yet the mechanisms for COM formation are still debated. It was originally believed that COMs initially form in ices which are then irradiated by UV radiation from the surrounding interstellar radiation field as well as forming protostars and subsequently photodesorbed into the gas-phase. However, starless and prestellar cores do not have internal protostars to heat-up and sublimate the ices. Alternative models using chemical energy have been developed to explain the desorption of COMs, yet in order to test these models robust measurements of COM abundances are needed toward representative samples of cores. We've conducted a large-sample survey of 31 starless and prestellar cores in the Taurus Molecular Cloud, detecting methanol (CH$_3$OH) in 100$\%$ of the cores targeted and acetaldehyde (CH$_3$CHO) in 70$\%$. At least two transition lines of each molecule were measured, allowing us to place tight constraints on excitation temperature, column density and abundance. Additional mapping of methanol revealed extended emission, detected down to A$_\mathrm{V}$ as low as $\sim$ 3 mag. We find complex organic molecules are detectable in the gas-phase and are being formed early, at least hundreds of thousands of years prior to star and planet formation. The precursor molecule, CH$_3$OH, may be chemically linked to the more complex CH$_3$CHO, however higher spatial resolution maps are needed to further test chemical models.


title: Melting holographic mesons by cooling a magnetized quark gluon plasma

abstract: We extend our analysis of holographic meson dissociation in the presence of an intense magnetic field. In addition to the previously known critical temperature above which the mesons melt, we found that for certain magnetic field intensities there exists a second lower critical temperature, below which stable mesons cease to exist. While we showed before that there is a range of high temperatures for which mesons can be melted by changing the magnetic field intensity, here we show that, as a consequence of the second critical point, there is also a range of low temperatures for which this phenomenom, which we term Magnetic Meson Melting (MMM), can be triggered. Additionaly, we also show that the magnetic field decreases the mass gap of the meson spectrum along with their masses. We are able to observe this by constructing a configuration that makes it possible to apply gauge/gravity methods to study fundamental degrees of freedom in a quark-gluon plasma subject to a magnetic field as intense as that expected in high energy collisions. This is achieved by the confection of a ten-dimensional background that is dual to the magnetized plasma and nonetheless permits the embedding of D7-branes in it. The main difference with previous approaches, which in consequence gives the novel results, is that the magnetic field retroacts in the geometry itself, as opposed to be confined to the world volume of the probe D7-branes.


title: Relative homological representations of framed mapping class groups

abstract: Let $\Sigma$ be a surface with either boundary or marked points, equipped with an arbitrary framing. In this note we determine the action of the associated "framed mapping class group" on the homology of $\Sigma$ relative to its boundary (respectively marked points), describing the image as the kernel of a certain crossed homomorphism related to classical spin structures. Applying recent work of the authors, we use this to describe the monodromy action of the orbifold fundamental group of a stratum of abelian differentials on the relative periods.


title: Framed mapping class groups and the monodromy of strata of Abelian differentials

abstract: This paper investigates the relationship between strata of abelian differentials and various mapping class groups afforded by means of the topological monodromy representation. Building off of prior work of the authors, we show that the fundamental group of a stratum surjects onto the subgroup of the mapping class group which preserves a fixed framing of the underlying Riemann surface, thereby giving a complete characterization of the monodromy group. In the course of our proof we also show that these "framed mapping class groups" are finitely generated (even though they are of infinite index) and give explicit generating sets.


title: A Novel Survey for Young Substellar Objects with the W-band filter II. The Coolest and Lowest Mass Members of the Serpens-South Star-forming Region

abstract: Given its relative proximity (~430 pc), compact size (< 20'), young age (~ 0.5 Myr) and rich number of young stellar objects, the Serpens-South star forming region is a promising site for studying young sub-stellar objects, yet the low-mass members of this region remain largely undiscovered. In this paper we report on a deep photometric survey using a custom 1.45 um filter (W-band), as well as standard J and H near-IR filters, in order to identify candidate low-mass young brown-dwarfs in the Serpens-South region. We constructed a reddening-insensitive index (Q) by combining J, H and W-band photometry for survey objects, in order to identify candidate low-mass members of Serpens based on the strength of the water absorption feature at 1.45 um in the atmospheres of mid-M and later objects. We then conducted spectroscopic follow up to confirm youth and spectral type for our candidates. This is the first survey to identify the very low-mass and coolest members of Serpens-South. We identify 4 low-mass candidate Serpens members, which all display IR excess emission, indicating the likely presence of circumstellar disks around them. One of the four candidate low-mass members in our list, SERP182918-020245, exhibits Pa_beta and Br_gamma emission features, confirming its youth and ongoing magnetospheric accretion. Our new candidate members have spectral types >M4 and are the coolest and lowest mass candidate members yet identified in Serpens-South.


title: Design of a Fully Actuated Robotic Hand With Multiple Gelsight Tactile Sensors

abstract: This work details the design of a novel two finger robot gripper with multiple Gelsight based optical-tactile sensors covering the inner surface of the hand. The multiple Gelsight sensors can gather the surface topology of the object from multiple views simultaneously as well as can track the shear and tensile stress. In addition, other sensing modalities enable the hand to gather the thermal, acoustic and vibration information from the object being grasped. The force controlled gripper is fully actuated so that it can be used for various grasp configurations and can also be used for in-hand manipulation tasks. Here we present the design of such a gripper.


title: Monitoring strong coupling in nonlocal plasmonics with electron spectroscopies

abstract: Plasmon-exciton polaritons provide exciting possibilities to control light-matter interactions at the nanoscale by enabling closer investigation of quantum optical effects and facilitating novel technologies based, for instance, on Bose-Einstein condensation and polaritonic lasing. Nevertheless, observing and visualising polaritons is challenging, and traditional optical microscopy techniques often lead to ambiguities regarding the emergence and strength of the plasmon-exciton coupling. Electron microscopy offers a more robust means to study and verify the nature of plexcitons, but is still hindered by instrument limitations and resolution. A simple theoretical description of electron beam-excited plexcitons is therefore vital to complement ongoing experimental efforts. Here we apply analytic solutions for the electron-loss and photon-emission probabilities to evaluate plasmon-exciton coupling studied either with the recently adopted technique of electron energy-loss spectroscopy, or with the so-far unexplored in this context cathodoluminescence spectroscopy. Foreseeing the necessity to account for quantum corrections in the plasmonic response, we extend these solutions within the framework of general nonlocal hydrodynamic descriptions. As a specific example we study core-shell spherical emitter-molecule hybrids, going beyond the standard local-response approximation through the hydrodynamic Drude model for screening and the generalised nonlocal optical response theory for nonlocal damping. We show that electron microscopies are extremely powerful in describing the interaction of emitters with the otherwise weakly excited by optical means higher-order plasmonic multipoles, a response that survives when quantum-informed models are considered.


title: The Sum Composition Problem

abstract: In this paper, we study the "sum composition problem" between two lists $A$ and $B$ of positive integers. We start by saying that $B$ is "sum composition" of $A$ when there exists an ordered $m$-partition $[A_1,\ldots,A_m]$ of $A$ where $m$ is the length of $B$ and the sum of each part $A_k$ is equal to the corresponding part of $B$. Then, we consider the following two problems: $i)$ the "exhaustive problem", consisting in the generation of all partitions of $A$ for which $B$ is sum composition of $A$, and $ii)$ the "existential problem", consisting in the verification of the existence of a partition of $A$ for which $B$ is sum composition of $A$. Starting from some general properties of the sum compositions, we present a first algorithm solving the exhaustive problem and then a second algorithm solving the existential problem. We also provide proofs of correctness and experimental analysis for assessing the quality of the proposed solutions along with a comparison with related works.


title: Interaction Networks from Discrete Event Data by Poisson Multivariate Mutual Information Estimation and Information Flow with Applications from Gene Expression Data

abstract: In this work, we introduce a new methodology for inferring the interaction structure of discrete valued time series which are Poisson distributed. While most related methods are premised on continuous state stochastic processes, in fact, discrete and counting event oriented stochastic process are natural and common, so called time-point processes (TPP). An important application that we focus on here is gene expression. Nonparameteric methods such as the popular k-nearest neighbors (KNN) are slow converging for discrete processes, and thus data hungry. Now, with the new multi-variate Poisson estimator developed here as the core computational engine, the causation entropy (CSE) principle, together with the associated greedy search algorithm optimal CSE (oCSE) allows us to efficiently infer the true network structure for this class of stochastic processes that were previously not practical. We illustrate the power of our method, first in benchmarking with synthetic datum, and then by inferring the genetic factors network from a breast cancer micro-RNA (miRNA) sequence count data set. We show the Poisson oCSE gives the best performance among the tested methods anfmatlabd discovers previously known interactions on the breast cancer data set.


title: Homogenization of periodic parabolic systems in the $L_2(\mathbb{R}^d)$-norm with the corrector taken into account

abstract: In $L_2(\mathbb{R}^d;\mathbb{C}^n)$, consider a self-adjoint matrix second order elliptic differential operator $\mathcal{B}_\varepsilon$, $0<\varepsilon \leqslant 1$. The principal part of the operator is given in a factorised form, the operator contains first and zero order terms. The operator $\mathcal{B}_\varepsilon$ is positive definite, its coefficients are periodic and depend on $\mathbf{x}/\varepsilon$. We study the behaviour in the small period limit of the operator exponential $e^{-\mathcal{B}_\varepsilon t}$, $t\geqslant 0$. The approximation in the $(L_2\rightarrow L_2)$-operator norm with error estimate of order $O(\varepsilon ^2)$ is obtained. The corrector is taken into account in this approximation. The results are applied to homogenization of the solutions for the Cauchy problem for parabolic systems.


title: Fractional derivative order determination from harmonic oscillator damping factor

abstract: This article analysis differential equations which represents damped and fractional oscillators. First, it is shown that prior to using physical quantities in fractional calculus, it is imperative that they are turned dimensionless. Afterwards, approximated expressions that relate the two equations parameters for the case that the fractional order is close to an integer number are presented. Following, a numerical regression is made using power series expansion, and, also from fractional calculus, the fact that both equations cannot be equivalent is concluded. In the end, from the numerical regression data, the analytical approximated expressions that relate the two equations' parameters are refined.


title: Extremal triangle-free and odd-cycle-free colourings of uncountable graphs

abstract: The optimality of the Erd\H{o}s-Rado theorem for pairs is witnessed by the colouring $\Delta_\kappa : [2^\kappa]^2 \rightarrow \kappa$ recording the least point of disagreement between two functions. This colouring has no monochromatic triangles or, more generally, odd cycles. We investigate a number of questions investigating the extent to which $\Delta_\kappa$ is an \emph{extremal} such triangle-free or odd-cycle-free colouring. We begin by introducing the notion of $\Delta$-regressive and almost $\Delta$-regressive colourings and studying the structures that must appear as monochromatic subgraphs for such colourings. We also consider the question as to whether $\Delta_\kappa$ has the minimal cardinality of any \emph{maximal} triangle-free or odd-cycle-free colouring into $\kappa$. We resolve the question positively for odd-cycle-free colourings.


title: Sensitivity Analysis in the Dupire Local Volatility Model with Tensorflow

abstract: In a recent paper, we have demonstrated how the affinity between TPUs and multi-dimensional financial simulation resulted in fast Monte Carlo simulations that could be setup in a few lines of python Tensorflow code. We also presented a major benefit from writing high performance simulations in an automated differentiation language such as Tensorflow: a single line of code enabled us to estimate sensitivities, i.e. the rate of change in price of financial instrument with respect to another input such as the interest rate, the current price of the underlying, or volatility. Such sensitivities (otherwise known as the famous financial "Greeks") are fundamental for risk assessment and risk mitigation. In the present follow-up short paper, we extend the developments exposed in our previous work about the use of Tensor Processing Units and Tensorflow for TPUs.


title: An Asymptotic Determinant Bound on 0-1 Matrices with Fixed Row Sums

abstract: This paper improves previously known bounds on the determinant of 0-1 matrices where each row has fixed support size. This uses a method based on Scheinerman's, with new analyses to improve upon his conjectures.


title: Poorly separated infinite normal products

abstract: A product of compact normal spaces is normal; the product of a countably infinite collection of non-trivial spaces is normal if and only if it is countably paracompact and each of its finite sub-products is normal; if all powers of a space X are normal then X is compact: provided in each case that the spaces involved are T1. Here I examine the situation for infinite products not required to be T1 (or regular), extending or generalizing each of these facts. In addition, I prove some related results, give a number of examples, explore some alternative proofs, and close with some speculation regarding potential applications of these findings to category theory and lattice theory.


title: Conservative numerical schemes with optimal dispersive wave relations -- Part II. Numerical evaluations

abstract: A new energy and enstrophy conserving scheme is evaluated using a suite of test cases over the global spherical domain or bounded domains. The evaluation is organized around a set of pre-defined properties: accuracy of individual opeartors, accuracy of the whole scheme, conservation, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. The scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that the new scheme can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

categories: math.NA cs.NA physics.ao-ph physics.geo-ph


title: Rigidity theorems for higher rank lattice actions

abstract: Let $\Gamma$ be a weakly irreducible higher rank lattice. In this paper, we will prove various rigidity results for the $\Gamma$-action following a philosophy of the Zimmer program. We provide new rigidity results including local and global rigidity of the $\Gamma$-action when $\Gamma$ does not have Property (T). The new ingredient is a dynamical cocycle super-rigidity theorem. It can be thought of as the generalization of Zimmer's cocycle super-rigidity theorem since it provides almost the same consequences dynamically and algebraically. This allows us to derive various our rigidity results using dynamical super-rigidity instead of Zimmer's cocycle super-rigidity theorem.


title: Calibrating GONG Magnetograms with End-to-end Instrument Simulation III: Comparison, Calibration, and Results

abstract: This is the last of three papers describing an `absolute' calibration of the GONG magnetograph using and end-to-end simulation of its measurement process. The simulation begins with a MURaM 3D MHD datacube and ends with a `synthetic magnetogram' of the corresponding magnetic field values as they would be observed by GONG. We determine a calibration by comparing the synthetic magnetic field measurements with the MURaM magnetic field values that produced them. The previous two papers have described the GONG measurement process (both instrument and data processing), our simulation of it, and the theory of magnetogram comparison and calibration. In this paper, we address some final points on calibration, combine all of this work into a set of calibration curves, and consider the results. We also review the results of the previous two papers for locality of reference. Our calibration indicates that GONG magnetograms underestimate weak flux by a factor of over 2 near disk center, but that factor decreases to about 1 as the line-of-sight approaches the limb. A preliminary investigation of the generalizability of these results suggests other instruments will be affected in a similar way. We also find that some differences in previous magnetograph comparisons are artifacts of instrumental resolution which do not reflect an intrinsic calibration difference, and the measurements are more similar than sometimes thought. These results are directly applicable to question of solar wind prediction model accuracies, particularly in the search for the cause of the common discrepancy between predicted solar wind magnetic flux at 1 AU and values measured in situ by current satellite missions.


title: Efficient Algorithms for Generating Provably Near-Optimal Cluster Descriptors for Explainability

abstract: Improving the explainability of the results from machine learning methods has become an important research goal. Here, we study the problem of making clusters more interpretable by extending a recent approach of [Davidson et al., NeurIPS 2018] for constructing succinct representations for clusters. Given a set of objects $S$, a partition $\pi$ of $S$ (into clusters), and a universe $T$ of tags such that each element in $S$ is associated with a subset of tags, the goal is to find a representative set of tags for each cluster such that those sets are pairwise-disjoint and the total size of all the representatives is minimized. Since this problem is NP-hard in general, we develop approximation algorithms with provable performance guarantees for the problem. We also show applications to explain clusters from datasets, including clusters of genomic sequences that represent different threat levels.


title: Cusps and $q$-expansion principles for modular curves at infinite level

abstract: We develop an analytic theory of cusps for Scholze's $p$-adic modular curves at infinite level in terms of perfectoid parameter spaces for Tate curves. As an application, we describe a canonical tilting isomorphism between an anticanonical overconvergent neighbourhood of the ordinary locus of the modular curve at level $\Gamma_1(p^\infty)$ and the analogous locus of an infinite level perfected Igusa variety. We also prove various $q$-expansion principles for functions on modular curves at infinite level, namely that the properties of extending to the cusps, vanishing, coming from finite level, and being bounded, can all be detected on $q$-expansions.


title: Calibrating GONG Magnetograms with End-to-end Instrument Simulation I: Background, the GONG Instrument, and End-to-end Simulation

abstract: This is the first of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. The input to this simulation is a MURaM 3D MHD photospheric simulation and the output is the corresponding set of simulated data numbers which would be recorded by the GONG detectors. These simulated data numbers are then used to produce `synthetic magnetograms' which can be compared with the simulation inputs. This paper describes the GONG instrument, the MURaM datacube, our instrument simulator, and calculation of synthetic magnetograms, setting the stage for the subsequent two papers. These will first lay groundwork for calibration (and magnetogram comparison in general), then apply them to calibration of GONG using the simulation results.


title: Calibrating GONG Magnetograms with End-to-end Instrument Simulation II: Theory of Calibration

abstract: This is the second of three papers describing an `absolute' calibration of the GONG magnetograph using an end-to-end simulation of its measurement process. In the first paper, we described the GONG instrument and our `end-to-end' simulation of its measurement process. In this paper, we consider the theory of calibration, and magnetograph comparison in general, identifying some of the significant issues and pitfalls. The calibration of a magnetograph is a function of whether or not it preserves flux, independent of its spatial resolution. However, we find that the one-dimensional comparison methods most often used for magnetograph calibration and comparison will show dramatic differences between two magnetograms with differing spatial resolution, even if they both preserve flux. Some of the apparent disagreement between magnetograms found in the literature are likely a result of these instrumental resolution differences rather than any intrinsic calibration differences. To avoid them, spatial resolution must be carefully matched prior to comparing magnetograms or making calibration curves. In the third paper, we apply the lessons learned here to absolute calibration of GONG using our `end-to-end' measurement simulation.


title: Quantum Phase Diagrams of Matter-Field Hamiltonians I: Fidelity, Bures Distance, and Entanglement

abstract: A general procedure is established to calculate the quantum phase diagrams for finite matter-field Hamiltonian models. The minimum energy surface associated to the different symmetries of the model is calculated as a function of the matter-field coupling strengths. By means of the ground state wave functions, one looks for minimal fidelity or maximal Bures distance surfaces in terms of the parameters, and from them the critical regions of those surfaces characterize the finite quantum phase transitions. Following this procedure for $N_a=1$ and $N_a=4$ particles, the quantum phase diagrams are calculated for the generalised Tavis-Cummings and Dicke models of 3-level systems interacting dipolarly with $2$ modes of electromagnetic field. For $N_a=1$, the reduced density matrix of the matter allows us to determine the phase regions in a $2$-simplex (associated to a general three dimensional density matrix), on the different $3$-level atomic configurations, together with a measurement of the quantum correlations between the matter and field sectors. As the occupation probabilities can be measured experimentally, the existence of a quantum phase diagram for a finite system can be established.


title: Consistency of a Recurrent Language Model With Respect to Incomplete Decoding

abstract: Despite strong performance on a variety of tasks, neural sequence models trained with maximum likelihood have been shown to exhibit issues such as length bias and degenerate repetition. We study the related issue of receiving infinite-length sequences from a recurrent language model when using common decoding algorithms. To analyze this issue, we first define inconsistency of a decoding algorithm, meaning that the algorithm can yield an infinite-length sequence that has zero probability under the model. We prove that commonly used incomplete decoding algorithms - greedy search, beam search, top-k sampling, and nucleus sampling - are inconsistent, despite the fact that recurrent language models are trained to produce sequences of finite length. Based on these insights, we propose two remedies which address inconsistency: consistent variants of top-k and nucleus sampling, and a self-terminating recurrent language model. Empirical results show that inconsistency occurs in practice, and that the proposed methods prevent inconsistency.


title: On Ridership and Frequency

abstract: Even before the start of the COVID-19 pandemic, bus ridership in the United States had attained its lowest level since 1973. If transit agencies hope to reverse this trend, they must understand how their service allocation policies affect ridership. This paper is among the first to model ridership trends on a hyper-local level over time. A Poisson fixed-effects model is developed to evaluate the ridership elasticity to frequency on weekdays using passenger count data from Portland, Miami, Minneapolis/St-Paul, and Atlanta between 2012 and 2018. In every agency, ridership is found to be elastic to frequency when observing the variation between individual route-segments at one point in time. In other words, the most frequent routes are already the most productive in terms of passengers per vehicle-trip. When observing the variation within each route-segment over time, however, ridership is inelastic; each additional vehicle-trip is expected to generate less ridership than the average bus already on the route. In three of the four agencies, the elasticity is a decreasing function of prior frequency, meaning that low-frequency routes are the most sensitive to changes in frequency. This paper can help transit agencies anticipate the marginal effect of shifting service throughout the network. As the quality and availability of passenger count data improve, this paper can serve as the methodological basis to explore the dynamics of bus ridership.


title: Domain structure dynamics in the ferromagnetic Kagome-lattice Weyl semimetal Co$_3$Sn$_2$S$_2$

abstract: Co$_3$Sn$_2$S$_2$, a Weyl semimetal that consists of layers of Kagome lattices, transitions from a high-temperature paramagnetic phase to a low-temperature ferromagnetic phase below 177 K. The phase transition occurs through an intermediate non-trivial magnetic phase, the so-called "A"-phase just below the Curie temperature. The "A"-phase was earlier linked with a competing anti-ferromagnetic phase, a spin-glass phase and certain indirect measurements indicated the possibility of magnetic Skyrmions in this phase. We have imaged the magnetic domain structure in a single crystal of Co$_3$Sn$_2$S$_2$ at different temperatures, magnetic fields and field-angles by magnetic force microscopy. At low temperatures, we observed stripe domains indicating presence of uniaxial anisotropy. Above 130 K, the domain walls become mobile and they tend to align relatively easily when the magnetic field is increased along the $c$-axis than in the $a-b$ plane. Our detailed study of field-dependent domain dynamics reveal that the anomalous nature of the phase transition just below $T_c$ is dominantly governed by domain wall motion.


title: Search for exotic spin-dependent interactions using polarized helium

abstract: We investigate the sensitivities of searches for exotic spin-dependent interactions between the polarized nuclear spins of $^3$He and the particles of unpolarized or polarized solid-state masses using the frequency method and the resonance method. In the frequency method, the spin-dependent interactions act as an effective static magnetic field, causing the frequency shift to the spin precession of $^{3}$He. In the resonance method, proposed by Arvanitaki and Geraci [Phys. Rev. Lett. 113, 161801 (2014)] for the significant improvement of the experimental sensitivities on the spin-dependent interactions, the mass movement is modulated at the Larmor frequency of $^3$He. This results in the modulating spin-dependent interactions inducing an effective oscillatory magnetic field, which can tilt the $^3$He spins, similarly as an oscillatory magnetic field in nuclear magnetic resonance. We estimate the sensitivities of the searches using a room-temperature $^3$He target for its extremely long relaxation time. New limits on the coupling strengths of the spin-dependent interactions can be set in the interaction length range below $10^{-1}$ m.


title: Implementing Inductive bias for different navigation tasks through diverse RNN attractors

abstract: Navigation is crucial for animal behavior and is assumed to require an internal representation of the external environment, termed a cognitive map. The precise form of this representation is often considered to be a metric representation of space. An internal representation, however, is judged by its contribution to performance on a given task, and may thus vary between different types of navigation tasks. Here we train a recurrent neural network that controls an agent performing several navigation tasks in a simple environment. To focus on internal representations, we split learning into a task-agnostic pre-training stage that modifies internal connectivity and a task-specific Q learning stage that controls the network's output. We show that pre-training shapes the attractor landscape of the networks, leading to either a continuous attractor, discrete attractors or a disordered state. These structures induce bias onto the Q-Learning phase, leading to a performance pattern across the tasks corresponding to metric and topological regularities. By combining two types of networks in a modular structure, we could get better performance for both regularities. Our results show that, in recurrent networks, inductive bias takes the form of attractor landscapes -- which can be shaped by pre-training and analyzed using dynamical systems methods. Furthermore, we demonstrate that non-metric representations are useful for navigation tasks, and their combination with metric representation leads to flexibile multiple-task learning.


title: On the limits of cross-domain generalization in automated X-ray prediction

abstract: This large scale study focuses on quantifying what X-rays diagnostic prediction tasks generalize well across multiple different datasets. We present evidence that the issue of generalization is not due to a shift in the images but instead a shift in the labels. We study the cross-domain performance, agreement between models, and model representations. We find interesting discrepancies between performance and agreement where models which both achieve good performance disagree in their predictions as well as models which agree yet achieve poor performance. We also test for concept similarity by regularizing a network to group tasks across multiple datasets together and observe variation across the tasks. All code is made available online and data is publicly available: https://github.com/mlmed/torchxrayvision


title: The Structure of Tidal Disruption Event Host Galaxies on Scales of Tens to Thousands of Parsecs

abstract: We explore the galaxy structure of four tidal disruption event (TDE) host galaxies on 30 pc to kpc scales using HST WFC3 multi-band imaging. The star formation histories of these hosts are diverse, including one post-starburst galaxy (ASASSN-14li), two hosts with recent weak starbursts (ASASSN-14ae and iPTF15af), and one early type (PTF09ge). Compared to early type galaxies of similar stellar masses, the TDE hosts have higher central surface brightnesses and stellar mass surface densities on 30-100 pc scales. The TDE hosts do not show the large, kpc-scale tidal disruptions seen in some post-starburst galaxies; the hosts have low morphological asymmetries similar to those of early type galaxies. The lack of strong asymmetries are inconsistent with a recent major (~1:1 mass) merger, although minor ($\lesssim$1:3) mergers are possible. Given the time elapsed since the end of the starbursts in the three post-burst TDE hosts and the constraints on the merger mass ratios, it is unlikely that a bound supermassive black hole binary (SMBHB) has had time to coalesce. The TDE hosts have low central (<140 pc) ellipticities compared to early type galaxies. The low central ellipticities disfavor a strong radial anisotropy as the cause for the enhanced TDE rate, although we cannot rule out eccentric disks at the scale of the black hole gravitational radius of influence (~1 pc). These observations suggest that the high central stellar densities are a more important driver than SMBHBs or radial anisotropies in increasing the TDE rate in galaxies with recent starbursts.


title: Optical/near-infrared observations of the Fried Egg Nebula: Multiple shell ejections on a 100 yr timescale from a massive yellow hypergiant

abstract: Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 {\mu}m, 11.85 {\mu}m and 12.81 {\mu}m simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 {\mu}m show that the Br{\gamma} emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Br{\gamma} has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms.


title: Spectral Analysis of the Schrodinger Operator with an Optical Potential

abstract: In this paper we give a complete description of the spectral analysis of the Schrodinger operator L(V) with the optical potentil. First we consider the Bolch eigenvalues and spectrum of L(V). Then using it we investigate spectral singularities and essential spectral singularities (ESS). We prove that the operator L(V) has no ESS and has ESS respectively if and only if V is not a critical point and V is a critical point. Using it we classify the spectral expansion in term of the critical points. Finally we discuss the critical points, formulate some conjectures and describe the changes of the spectrum of L(V) when V changes.


title: Classical Tools for Antipodal Identification in Reissner-Nordstr\"om Spacetime

abstract: We extend the discussion of the antipodal identification of black holes to the Reissner-Nordstr\"om (RN) spacetime by developing the classical tools necessary to define the corresponding quantum field theory (QFT). We solve the massless Klein-Gordon equation in the RN background in terms of scattering coefficients and provide a procedure for constructing a solution for an arbitrary analytic extension of RN. The behavior of the maximally extended solution is highly dependent upon the coefficients of scattering between the inner and outer horizons, so we present the low-frequency behavior of, and numerical solutions for, these quantities. We find that, for low enough frequency, field amplitudes of solutions with purely positive or negative frequency at each horizon will acquire only a phase after passing both the inner and outer horizons, while at higher frequencies the amplitudes will tend to grow exponentially either to the future or to the past, and decay exponentially in the other direction. Regardless, we can always construct a basis of globally antipodal symmetric and antisymmetric solutions for any finite analytic extension of RN. We have characterized this basis in terms of positive and negative frequency solutions for future use in constructing the corresponding QFT.


title: On uniform convergence of the inverse Fourier transform for differential equations and Hamiltonian systems with degenerating weight

abstract: We study pseudospectral and spectral functions for Hamiltonian system $Jy'-B(t)=\lambda\Delta(t)y$ and differential equation $l[y]=\lambda\Delta(t)y$ with matrix-valued coefficients defined on an interval $\mathcal{I}=[a,b)$ with the regular endpoint $a$. It is not assumed that the matrix weight $\Delta(t)\geq 0$ is invertible a.e. on $\mathcal{I}$. In this case a pseudospectral function always exists, but the set of spectral functions may be empty. We obtain a parametrization $\sigma=\sigma_\tau$ of all pseudospectral and spectral functions $\sigma$ by means of a Nevanlinna parameter $\tau$ and single out in terms of $\tau$ and boundary conditions the class of functions $y$ for which the inverse Fourier transform $y(t)=\int\limits_{\mathbb{R}} \varphi(t,s)\, d\sigma (s) \widehat y(s)$ converges uniformly. We also show that for scalar equation $l[y]=\lambda \Delta(t)y$ the set of spectral functions is not empty. This enables us to extend the Kats-Krein and Atkinson results for scalar Sturm - Liouville equation $-(p(t)y')'+q(t)y=\lambda \Delta (t) y$ to such equations with arbitrary coefficients $p(t)$ and $q(t)$ and arbitrary non trivial weight $\Delta (t)\geq 0$.


title: Filtered Lattice Boltzmann Collision Formulation Enforcing Isotropy and Galilean Invariance

abstract: The central framework of a filtered lattice Boltzmann collision operator formulation is to remove hydrodynamic moments that are not supported by the order of isotropy of a given lattice velocity set. Due to the natural moment orthogonality of the Hermite polynomials, the form of a filtered collision operator is obtained directly via truncation of the Hermite expansion. In this paper, we present an extension of the filtered collision operator formulation to enforce Galilean invariance. This is accomplished by representing hydrodynamic moments in the relative reference frame with respect to local fluid velocity. The resulting collision operator has a compact and fully Galilean invariant form, and it can then be exactly expressed in terms of an infinite Hermite expansion. Giving a lattice velocity set of specific order of isotropy, a proper truncation of this expansion can be directly determined. Higher order terms are retained in the truncation if a higher order lattice velocity set is used, so that Galilean invariance is attained asymptotically. The previously known filtered collision operator forms can be seen as a limit of zero fluid velocity.


title: Attosecond spectroscopy reveals alignment dependent core-hole dynamics in the ICl molecule

abstract: The removal of electrons located in the core shells of molecules creates transient states that live between a few femtoseconds to attoseconds. Owing to these short lifetimes, time-resolved studies of these states are challenging and complex molecular dynamics driven solely by electronic correlation are difficult to observe. Here, few-femtosecond core-excited state lifetimes of iodine monochloride are obtained by attosecond transient absorption on iodine 4d-16p transitions around 55 eV. Core-level ligand field splitting allows direct access of excited states aligned along and perpendicular to the ICl molecular axis. Lifetimes of 3.5 fs and 4.3 fs are obtained for core-hole states parallel to the bond and 6.5 fs and 6.9 fs for perpendicular states, while nuclear motion is essentially frozen on this timescale. Theory shows that the dramatic decrease of lifetime for core-vacancies parallel to the covalent bond is a manifestation of non-local interactions with the neighboring Cl atom of ICl.


title: Relativistic coupled-cluster study of BaF in search of $\mathcal{CP}$ violation

abstract: BaF is one of the potential candidates for the experimental search of the electric dipole moment of the electron (eEDM). The NL-eEDM collaboration is building a new experimental set up to measure the eEDM using the BaF molecule [The NL-eEDM collaboration, Eur. Phys. J. D (2018) 72: 197]. To analyze the results of such an experiment, one would require the accurate value of the molecular ${\mathcal{P,T}}$-odd interaction parameters that cannot be measured from any experiment. In this work, we report the precise value of the ${\mathcal{P,T}}$-odd interaction parameters of the BaF molecule obtained from the four-component relativistic coupled-cluster calculations. We also calculate the hyperfine structure (HFS) constants of the same molecule to assess the reliability of the reported molecular parameters. The calculated HFS constants show good agreement with the available experimental values. Further, the systematic effects of electron-correlation along with the roles of inner-core electrons and the virtual energy functions in the calculation of the studied properties of BaF are investigated.


title: Continuous Geodesic Convolutions for Learning on 3D Shapes

abstract: The majority of descriptor-based methods for geometric processing of non-rigid shape rely on hand-crafted descriptors. Recently, learning-based techniques have been shown effective, achieving state-of-the-art results in a variety of tasks. Yet, even though these methods can in principle work directly on raw data, most methods still rely on hand-crafted descriptors at the input layer. In this work, we wish to challenge this practice and use a neural network to learn descriptors directly from the raw mesh. To this end, we introduce two modules into our neural architecture. The first is a local reference frame (LRF) used to explicitly make the features invariant to rigid transformations. The second is continuous convolution kernels that provide robustness to sampling. We show the efficacy of our proposed network in learning on raw meshes using two cornerstone tasks: shape matching, and human body parts segmentation. Our results show superior results over baseline methods that use hand-crafted descriptors.


title: Giant Stark splitting of an exciton in bilayer MoS$_2$

abstract: Transition metal dichalcogenides (TMDs) constitute a versatile platform for atomically thin optoelectronics devices and spin-valley memory applications. In monolayers optical absorption is strong, but the transition energy is not tunable as the neutral exciton has essentially no out-of-plane electric dipole. In contrast, interlayer exciton transitions in heterobilayers are widely tunable in applied electric fields, but their coupling to light is considerably reduced. Here, we show tuning over 120 meV of interlayer excitons with high oscillator strength in bilayer MoS2. These shifts are due to the quantum confined Stark effect, here the electron is localised to one of the layers yet the hole is delocalised across the bilayer. We optically probe the interaction between intra- and interlayer excitons as they are energetically tuned into resonance. This allows studying their mixing supported by beyond standard density functional theory calculations including excitonic effects. In MoS2 trilayers our experiments uncover two types of interlayer excitons with and without in-built electric dipoles, respectively. Highly tunable excitonic transitions with large oscillator strength and in-built dipoles, that lead to considerable exciton-exciton interactions, hold great promise for non-linear optics with polaritons.


title: Differentially Quantized Gradient Methods

abstract: Consider the following distributed optimization scenario. A worker has access to training data that it uses to compute the gradients while a server decides when to stop iterative computation based on its target accuracy or delay constraints. The server receives all its information about the problem instance from the worker via a rate-limited noiseless communication channel. We introduce the principle we call Differential Quantization (DQ) that prescribes compensating the past quantization errors to direct the descent trajectory of a quantized algorithm towards that of its unquantized counterpart. Assuming that the objective function is smooth and strongly convex, we prove that Differentially Quantized Gradient Descent (DQ-GD) attains a linear contraction factor of $\max\{\sigma_{\mathrm{GD}}, \rho_n 2^{-R}\}$, where $\sigma_{\mathrm{GD}}$ is the contraction factor of unquantized gradient descent (GD), $\rho_n \geq 1$ is the covering efficiency of the quantizer, and $R$ is the bitrate per problem dimension $n$. Thus at any $R\geq\log_2 \rho_n /\sigma_{\mathrm{GD}}$ bits, the contraction factor of DQ-GD is the same as that of unquantized GD, i.e., there is no loss due to quantization. We show that no algorithm within a certain class can converge faster than $\max\{\sigma_{\mathrm{GD}}, 2^{-R}\}$. Since quantizers exist with $\rho_n \to 1$ as $n \to \infty$ (Rogers, 1963), this means that DQ-GD is asymptotically optimal. The principle of differential quantization continues to apply to gradient methods with momentum such as Nesterov's accelerated gradient descent, and Polyak's heavy ball method. For these algorithms as well, if the rate is above a certain threshold, there is no loss in contraction factor obtained by the differentially quantized algorithm compared to its unquantized counterpart. Experimental results on least-squares problems validate our theoretical analysis.


title: Scalable Communication Endpoints for MPI+Threads Applications

abstract: Hybrid MPI+threads programming is gaining prominence as an alternative to the traditional "MPI everywhere'" model to better handle the disproportionate increase in the number of cores compared with other on-node resources. Current implementations of these two models represent the two extreme cases of communication resource sharing in modern MPI implementations. In the MPI-everywhere model, each MPI process has a dedicated set of communication resources (also known as endpoints), which is ideal for performance but is resource wasteful. With MPI+threads, current MPI implementations share a single communication endpoint for all threads, which is ideal for resource usage but is hurtful for performance. In this paper, we explore the tradeoff space between performance and communication resource usage in MPI+threads environments. We first demonstrate the two extreme cases---one where all threads share a single communication endpoint and another where each thread gets its own dedicated communication endpoint (similar to the MPI-everywhere model) and showcase the inefficiencies in both these cases. Next, we perform a thorough analysis of the different levels of resource sharing in the context of Mellanox InfiniBand. Using the lessons learned from this analysis, we design an improved resource-sharing model to produce \emph{scalable communication endpoints} that can achieve the same performance as with dedicated communication resources per thread but using just a third of the resources.


title: Isoperimetric estimates in low dimensional Riemannian products

abstract: Let $(T^k,h_k)=(S_{r_1}^1\times S_{r_2}^1 \times ... \times S_{r_k}^1, dt_1^2+dt_2^2+...+dt_k^2)$ be flat tori, $r_k\geq ...\geq r_2\geq r_1>0$ and $(\mathbb R^n,g_E)$ the Euclidean space with the flat metric. We compute the isoperimetric profile of $(T^2\times \mathbb R^n, h_2+g_E)$, $2\leq n\leq 5$, for small and big values of the volume. These computations give explicit lower bounds for the isoperimetric profile of $T^2\times\mathbb R^n$. We also note that similar estimates for $(T^k\times \mathbb R^n, h_k+g_E)$, $2\leq k\leq5$, $2\leq n\leq 7-k$, may be computed, provided estimates for $(T^{k-1}\times \mathbb R^{n+1}, h_{k-1}+g_E)$ exist. We compute this explicitly for $k=3$. We use symmetrization techniques for product manifolds, based on work of A. Ros and F. Morgan.


title: Introducing Aspects of Creativity in Automatic Poetry Generation

abstract: Poetry Generation involves teaching systems to automatically generate text that resembles poetic work. A deep learning system can learn to generate poetry on its own by training on a corpus of poems and modeling the particular style of language. In this paper, we propose taking an approach that fine-tunes GPT-2, a pre-trained language model, to our downstream task of poetry generation. We extend prior work on poetry generation by introducing creative elements. Specifically, we generate poems that express emotion and elicit the same in readers, and poems that use the language of dreams---called dream poetry. We are able to produce poems that correctly elicit the emotions of sadness and joy 87.5 and 85 percent, respectively, of the time. We produce dreamlike poetry by training on a corpus of texts that describe dreams. Poems from this model are shown to capture elements of dream poetry with scores of no less than 3.2 on the Likert scale. We perform crowdsourced human-evaluation for all our poems. We also make use of the Coh-Metrix tool, outlining metrics we use to gauge the quality of text generated.