title: Electronic ferroelectricity and magneto-electric effect in phase separated magnetic oxidesElectronic ferroelectricity and magneto-electric effect in phase separated magnetic oxides

abstract: We consider phase separated states in magnetic oxides (MO) thin films. We show that these states have a non-zero electric polarization. Moreover, the polarization is intimately related to a spatial distribution of magnetization in the film. Polarized states with opposite polarization and opposite magnetic configuration are degenerate. An external electric field removes the degeneracy and allows to switch between the two states. So, one can control electric polarization and magnetic configuration of the phase separated MO thin film with the external electric field.


title: Nonlinear Multiview Analysis: Identifiability and Neural Network-assisted Implementation

abstract: Multiview analysis aims at extracting shared latent components from data samples that are acquired in different domains, e.g., image, text, and audio. Classic multiview analysis, e.g., canonical correlation analysis (CCA), tackles this problem via matching the linearly transformed views in a certain latent domain. More recently, powerful nonlinear learning tools such as kernel methods and neural networks are utilized for enhancing the classic CCA. However, unlike linear CCA whose theoretical aspects are clearly understood, nonlinear CCA approaches are largely intuition-driven. In particular, it is unclear under what conditions the shared latent components across the views can be identified---while identifiability plays an essential role in many applications. In this work, we revisit nonlinear multiview analysis and address both the theoretical and computational aspects. Our work leverages a useful nonlinear model, namely, the post-nonlinear model, from the nonlinear mixture separation literature. Combining with multiview data, we take a nonlinear multiview mixture learning viewpoint, which is a natural extension of the classic generative models for linear CCA. From there, we derive a learning criterion. We show that minimizing this criterion leads to identification of the latent shared components up to certain ambiguities, under reasonable conditions. Our derivation and formulation also offer new insights and interpretations to existing deep neural network-based CCA formulations. On the computation side, we propose an effective algorithm with simple and scalable update rules. A series of simulations and real-data experiments corroborate our theoretical analysis.


title: Monochromatic Components in Edge-Coloured Graphs with Large Minimum Degree

abstract: For every $n\in\mathbb{N}$ and $k\geq2$, it is known that every $k$-edge-colouring of the complete graph on $n$ vertices contains a monochromatic connected component of order at least $\frac{n}{k-1}$. For $k\geq3$, it is known that the complete graph can be replaced by a graph $G$ with $\delta(G)\geq(1-\varepsilon_k)n$ for some constant $\varepsilon_k$. In this paper, we show that the maximum possible value of $\varepsilon_3$ is $\frac16$. This disproves a conjecture of Gy\'{a}rfas and S\'{a}rk\"{o}zy.


title: Kelvin-Helmholtz versus Tearing Instability: What Drives Turbulence in Stochastic Reconnection?

abstract: Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode or Kelvin-Helmholtz, among other plasma instabilities, could generate turbulence from irregular current sheets. We address the nature of driving mechanism for this process and consider a relative role of tearing and Kelvin-Helmholtz instabilities for the process of turbulence generation. In particular, we analyze the conditions for development of these two instabilities within three-dimensional reconnection regions. We show that both instabilities can excite turbulence fluctuations in reconnection regions. However, tearing mode has relatively slow growth rate, and at later times it becomes partially suppressed by transverse to the current sheet component of magnetic field, generated during the growth of turbulent fluctuations. On the contrary, the Kelvin-Helmholtz instability establishes quickly in the outflow region, and at later times it dominates the turbulence generation comparing to the contribution from tearing mode. Our results demonstrate that the tearing instability is subdominant compared to the the Kelvin-Helmholtz instability in terms of generation of turbulence in the 3D reconnection layers and therefore the self-driven reconnection is turbulent reconnection with tearing instability being important only at the initial stage of the reconnection.

categories: astro-ph.HE astro-ph.GA astro-ph.SR physics.plasm-ph physics.space-ph


title: Theory of nonequilibrium free energy transduction by molecular machines

abstract: Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared to more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

categories: cond-mat.stat-mech cond-mat.soft physics.bio-ph physics.chem-ph q-bio.SC


title: Dictionary Learning for Channel Estimation in Hybrid Frequency-Selective mmWave MIMO Systems

abstract: Exploiting channel sparsity at millimeter wave (mmWave) frequencies reduces the high training overhead associated with the channel estimation stage. Compressive sensing (CS) channel estimation techniques usually adopt the (overcomplete) wavelet/Fourier transform matrix as a sparsifying dictionary. This may not be the best choice when considering non-uniform arrays, antenna gain/phase errors, mutual coupling effects, etc. We propose two dictionary learning (DL) algorithms to learn the best sparsifying dictionaries for channel matrices from observations obtained with hybrid frequency-selective mmWave multiple-input-multiple-output (MIMO) systems. First, we optimize the combined dictionary, i.e., the Kronecker product of transmit and receive dictionaries, as it is used in practice to sparsify the channel matrix. Second, considering the different array structures at the transmitter and receiver, we exploit separable DL to find the best transmit and receive dictionaries. Once the channel is expressed in terms of the optimized dictionaries, various CS-based sparse recovery techniques can be applied for low overhead channel estimation. The proposed DL algorithms perform well under low SNR conditions inherent to any mmWave communication systems before the precoders/combiners can be optimized. The effectiveness of the proposed DL algorithms has been corroborated via numerical simulations with different system configurations, array geometries and hardware impairments.


title: Synthesizing carbon nanotubes in space

abstract: Context. As the 4th most abundant element in the universe, carbon (C) is widespread in the interstellar medium (ISM) in various allotropic forms (e.g., fullerenes have been identified unambiguously in many astronomical environments, the presence of polycyclic aromatic hydrocarbon molecules in space has been commonly admitted, and presolar graphite as well as nanodiamonds have been identified in meteorites). As stable allotropes of these species, whether carbon nanotubes (CNTs) and their hydrogenated counterparts are also present in the ISM or not is unknown. Aims. We explore the possible routes for the formation of CNTs in the ISM and calculate their fingerprint vibrational spectral features in the infrared (IR). Methods. We study the hydrogen-abstraction/acetylene-addition (HACA) mechanism and investigate the synthesis of nanotubes using density functional theory (DFT). The IR vibrational spectra of CNTs and hydrogenated nanotubes (HNTs), as well as their cations, have also been obtained with DFT. Results. We find that CNTs could be synthesized in space through a feasible formation pathway. CNTs and cationic CNTs, as well as their hydrogenated counterparts, exhibit intense vibrational transitions in the IR. Their possible presence in the ISM could be investigated by comparing the calculated vibrational spectra with astronomical observations made by the Infrared Space Observatory, Spitzer Space Telescope, and particularly the upcoming James Webb Space Telescope.


title: Hybrid Inexact BCD for Coupled Structured Matrix Factorization in Hyperspectral Super-Resolution

abstract: This paper develops a first-order optimization method for coupled structured matrix factorization (CoSMF) problems that arise in the context of hyperspectral super-resolution (HSR) in remote sensing. To best leverage the problem structures for computational efficiency, we introduce a hybrid inexact block coordinate descent (HiBCD) scheme wherein one coordinate is updated via the fast proximal gradient (FPG) method, while another via the Frank-Wolfe (FW) method. The FPG-type methods are known to take less number of iterations to converge, by numerical experience, while the FW-type methods can offer lower per-iteration complexity in certain cases; and we wish to take the best of both. We show that the limit points of this HiBCD scheme are stationary. Our proof treats HiBCD as an optimization framework for a class of multi-block structured optimization problems, and our stationarity claim is applicable not only to CoSMF but also to many other problems. Previous optimization research showed the same stationarity result for inexact block coordinate descent with either FPG or FW updates only. Numerical results indicate that the proposed HiBCD scheme is computationally much more efficient than the state-of-the-art CoSMF schemes in HSR.


title: The Gauss map of polyhedral vertex stars

abstract: In discrete differential geometry, it is widely believed that the discrete Gaussian curvature of a polyhedral vertex star equals the algebraic area of its Gauss image. However, no complete proof has yet been described. We present an elementary proof in which we compare, for a particular normal vector, its winding numbers around the Gauss image and its antipode with its critical point index. This index is closely related to the normal degree of the Gauss map. We deduce how the number of inflection faces is related to the numbers of positive and negative components in the Gauss image. The resulting formula significantly limits the possible shapes of the Gauss image of a polyhedral vertex star. For example, if a triangulated vertex star is in general position and its Gauss image has no self-intersections, then it is either a convex spherical polygon if the curvature is positive, or a spherical pseudo-quadrilateral if the curvature is negative.


title: Oxygen Vacancy in ZnO-$w$ Phase: Pseudohybrid Hubbard Density Functional Study

abstract: The study of zinc oxide, within the homogeneous electron gas approximation, results in overhybridization of zinc $3d$ shell with oxygen $2p$ shell, a problem shown for most transition metal chalcogenides. This problem can be partially overcome by using LDA+$U$ (or, GGA+$U$) methodology. However, in contrast to the zinc $3d$ orbital, Hubbard type correction is typically excluded for the oxygen $2p$ orbital. In this work, we provide results of electronic structure calculations of an oxygen vacancy in ZnO supercell from ab initio perspective, with two Hubbard type corrections, $U_{\mathrm{Zn}-3d}$ and $U_{\mathrm{O}-2p}$. The results of our numerical simulations clearly reveal that the account of $U_{\mathrm{O}-2p}$ has a significant impact on the properties of bulk ZnO, in particular the relaxed lattice constants, effective mass of charge carriers as well as the bandgap. For a set of validated values of $U_{\mathrm{Zn}-3d}$ and $U_{\mathrm{O}-2p}$ we demonstrate the appearance of a localized state associated with the oxygen vacancy positioned in the bandgap of the ZnO supercell. Our numerical findings suggest that the defect state is characterized by the highest overlap with the conduction band states as obtained in the calculations with no Hubbard-type correction included. We argue that the electronic density of the defect state is primarily determined by Zn atoms closest to the vacancy.


title: Revisit Policy Optimization in Matrix Form

abstract: In tabular case, when the reward and environment dynamics are known, policy evaluation can be written as $\bm{V}_{\bm{\pi}} = (I - \gamma P_{\bm{\pi}})^{-1} \bm{r}_{\bm{\pi}}$, where $P_{\bm{\pi}}$ is the state transition matrix given policy ${\bm{\pi}}$ and $\bm{r}_{\bm{\pi}}$ is the reward signal given ${\bm{\pi}}$. What annoys us is that $P_{\bm{\pi}}$ and $\bm{r}_{\bm{\pi}}$ are both mixed with ${\bm{\pi}}$, which means every time when we update ${\bm{\pi}}$, they will change together. In this paper, we leverage the notation from \cite{wang2007dual} to disentangle ${\bm{\pi}}$ and environment dynamics which makes optimization over policy more straightforward. We show that policy gradient theorem \cite{sutton2018reinforcement} and TRPO \cite{schulman2015trust} can be put into a more general framework and such notation has good potential to be extended to model-based reinforcement learning.


title: Hausdorff dimension of non-conical limit sets

abstract: Geometrically infinite Kleinain groups have nonconical limit sets with the cardinality of the continuum. In this paper, we construct a geometrically infinite Fuchsian group such that the Hausdorff dimension of the nonconical limit set equals zero. For finitely generated, geometrically infinite Kleinian groups, we prove that the Hausdorff dimension of the nonconical limit set is positive.


title: Colloquium: Spintronics in graphene and other two-dimensional materials

abstract: After the first unequivocal demonstration of spin transport in graphene (Tombros et al., 2007), surprisingly at room temperature, it was quickly realized that this novel material was relevant for both fundamental spintronics and future applications. Over the decade since, exciting results have made the field of graphene spintronics blossom, and a second generation of studies has extended to new two-dimensional (2D) compounds. This Colloquium reviews recent theoretical and experimental advances on electronic spin transport in graphene and related 2D materials, focusing on emergent phenomena in van der Waals heterostructures and the new perspectives provided by them. These phenomena include proximity-enabled spin-orbit effects, the coupling of electronic spin to light, electrical tunability, and 2D magnetism.


title: Influence of magnetic ordering on the spectral properties of binary transition metal oxides

abstract: Using the ab initio embedded DMFT (eDMFT) approach, we study the effect of long-range magnetic ordering on the spectral properties in the binary transition metal oxides, and find that the most significant changes appear in the momentum resolved spectral functions, which sharpen into quite well-defined bands in the antiferromagnetic (AFM) phase. The strongest change across the transition is found at the topmost valence band edge (VBE), which is commonly associated with the Zhang-Rice bound state. This VBE state strengthens in the AFM phase, but only for the minority spin component, which is subject to stronger fluctuations. A similar hybridized VBE state also appears in the DFT single-particle description of the AFM phase, but gets much stronger and acquires a well-defined energy in the eDMFT description.


title: Role of Poloidal $\mathbf{E}\times\mathbf{B}$ Drift in Divertor Heat Transport in DIII-D

abstract: Simulations for DIII-D high confinement mode plasmas with the multifluid code UEDGE show a strong role of poloidal $\mathbf{E}\times\mathbf{B}$ drifts on divertor heat transport, challenging the paradigm of conduction limited scrape-off layer (SOL) transport. While simulations with reduced drift magnitude are well aligned with the assumption that electron heat conduction dominates the SOL heat transport, simulations with drifts predict that the poloidal convective $\mathbf{E}\times\mathbf{B}$ heat transport dominates over electron heat conduction in both attached and detached conditions. Since poloidal $\mathbf{E}\times\mathbf{B}$ flow propagates across magnetic field lines, poloidal transport with shallow magnetic pitch angles can reach values that are of the same order as would be provided by sonic flows parallel to the field lines. These flows can lead to strongly convection dominated divertor heat transport, increasing the poloidal volume of radiative power front, consistent with previous measurements at DIII-D. Due to these convective flows, the Lengyel integral approach, assuming zero convective fraction, is expected to provide a pessimistic estimate for radiative capability of impurities in the divertor. For the DIII-D simulations shown here, the Lengyel integral approach underestimates the radiated power by a factor of 6, indicating that for reliable DIII-D divertor power exhaust predictions, full 2D calculations, including drifts, would be necessary.


title: Ultra high energy cosmic rays from super-heavy dark matter in the context of large exposure observatories

abstract: The origin of the ultra high energy cosmic rays (UHECRs, $E>10^{18}$ eV) is still uncertain. However, great progress has been achieved due to the data taken by The Pierre Auger and Telescope Array observatories. The UHECR flux presents two main features, a hardening of the spectrum known as the ankle and a suppression at higher energies. The experimental data suggest that above the ankle the UHECRs flux is dominated by an extragalactic component of astrophysical origin. However, a minority component of exotic origin that dominates the flux beyond the suppression is still compatible with current data. Therefore, there exist the possibility that part of the UHECR flux originates from the decay of super-heavy dark matter particles clustered in the halos of the galaxies. In these scenarios the main contribution comes from the halo of our galaxy. In this article the possibility of identifying these scenarios in the context of the future very large exposure cosmic rays observatories is studied. It is worth mentioning that the contribution of the extragalactic halos located in the nearby universe is also included in these studies.


title: Learning Sparse Mixture of Experts for Visual Question Answering

abstract: There has been a rapid progress in the task of Visual Question Answering with improved model architectures. Unfortunately, these models are usually computationally intensive due to their sheer size which poses a serious challenge for deployment. We aim to tackle this issue for the specific task of Visual Question Answering (VQA). A Convolutional Neural Network (CNN) is an integral part of the visual processing pipeline of a VQA model (assuming the CNN is trained along with entire VQA model). In this project, we propose an efficient and modular neural architecture for the VQA task with focus on the CNN module. Our experiments demonstrate that a sparsely activated CNN based VQA model achieves comparable performance to a standard CNN based VQA model architecture.


title: Running of the top quark mass from proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: The running of the top quark mass is experimentally investigated for the first time. The mass of the top quark in the modified minimal subtraction ($\mathrm{\overline{MS}}$) renormalization scheme is extracted from a comparison of the differential top quark-antiquark ($\mathrm{t\overline{t}}$) cross section as a function of the invariant mass of the $\mathrm{t\overline{t}}$ system to next-to-leading-order theoretical predictions. The differential cross section is determined at the parton level by means of a maximum-likelihood fit to distributions of final-state observables. The analysis is performed using $\mathrm{t\overline{t}}$ candidate events in the e$^\pm\mu^\mp$ channel in proton-proton collision data at a centre-of-mass energy of 13 TeV recorded by the CMS detector at the CERN LHC in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The extracted running is found to be compatible with the scale dependence predicted by the corresponding renormalization group equation. In this analysis, the running is probed up to a scale of the order of 1 TeV.


title: Probing new physics using Rydberg states of atomic hydrogen

abstract: We consider the role of high-lying Rydberg states of simple atomic systems such as $^1$H in setting constraints on physics beyond the Standard Model. We obtain highly accurate bound states energies for a hydrogen atom in the presence of an additional force carrier (the energy levels of the Hellmann potential). These results show that varying the size and shape of the Rydberg state by varying the quantum numbers provides a way to probe the range of new forces. By combining these results with the current state-of-the-art QED corrections, we determine a robust global constraint on new physics that includes all current spectroscopic data in hydrogen. Lastly we show that improved measurements that fully exploit modern cooling and trapping methods as well as higher-lying states could lead to a strong, statistically robust global constraint on new physics based on laboratory measurements only.


title: Real-time Monitoring of Cellular Cultures with Electrolyte-gated Carbon Nanotube Transistors

abstract: Cell-based biosensors constitute a fundamental tool in biotechnology, and their relevance has greatly increased in recent years as a result of a surging demand for reduced animal testing and for high-throughput and cost-effective in vitro screening platforms dedicated to environmental and biomedical diagnostics, drug development and toxicology. In this context, electrochemical/electronic cell-based biosensors represent a promising class of devices that enable long-term and real-time monitoring of cell physiology in a non-invasive and label-free fashion, with a remarkable potential for process automation and parallelization. Common limitations of this class of devices at large include the need for substrate surface modification strategies to ensure cell adhesion and immobilization, limited compatibility with complementary optical cell-probing techniques, and need for frequency-dependent measurements, which rely on elaborated equivalent electrical circuit models for data analysis and interpretation. We hereby demonstrate the monitoring of cell adhesion and detachment through the time-dependent variations in the quasi-static characteristic current curves of a highly stable electrolyte-gated transistor, based on an optically transparent network of printable polymer-wrapped semiconducting carbon-nanotubes.


title: Strategic Scientific Plan for Gemini Observatory

abstract: We present the Strategic Scientific Plan (SSP) for the direction and activities of the Gemini Observatory in the 2020s. The overarching goal is to ensure that Gemini best uses the available resources to serve the needs of its international user community throughout the coming decade. The actionable items fall into three general categories: (1) preserving Gemini's current facilities and strengths; (2) developing instrumentation and software systems, including data pipelines, to enable new scientific capabilities that build on those strengths; (3) strategizing how visiting instruments can deliver additional valuable capabilities. We provide a high-level timeline (schematically illustrated in one figure) for the main developments discussed in this SSP. The schedule is ambitious, but in light of the recent Gemini in the Era of Multi-Messenger Astronomy (GEMMA) award from the NSF, the plan becomes achievable. Lists of milestones are given for gauging progress. As these milestones are reached and new instruments become available, some current instruments will need to be retired; we make recommendations in this regard. The final section concludes by reemphasizing the importance of a strong partnership committed to the needs of all members.


title: Perovskite PV-powered RFID: enabling low-cost self-powered IoT sensors

abstract: Photovoltaic (PV) cells have the potential to serve as on-board power sources for low-power IoT devices. Here, we explore the use of perovskite solar cells to power Radio Frequency (RF) backscatter-based IoT devices with a few {\mu}W power demand. Perovskites are suitable for low-cost, high-performance, low-temperature processing, and flexible light energy harvesting that hold the possibility to significantly extend the range and lifetime of current backscatter techniques such as Radio Frequency Identification (RFID). For these reasons, perovskite solar cells are prominent candidates for future low-power wireless applications. We report on realizing a functional perovskite-powered wireless temperature sensor with 4 m communication range. We use a 10.1% efficient perovskite PV module generating an output voltage of 4.3 V with an active area of 1.06 cm2 under 1 sun illumination, with AM 1.5G spectrum, to power a commercial off-the-shelf RFID IC, requiring 10 - 45 {\mu}W of power. Having an on-board energy harvester provides extra-energy to boost the range of the sensor (5x) in addition to providing energy to carry out high-volume sensor measurements (hundreds of measurements per min). Our evaluation of the prototype suggests that perovskite photovoltaic cells are able to meet the energy needs to enable fully autonomous low-power RF backscatter applications of the future. We conclude with an outlook into a range of applications that we envision to leverage the synergies offered by combining perovskite photovoltaics and RFID.


title: Methods, Models, and the Evolution of Moral Psychology

abstract: Why are we good? Why are we bad? Questions regarding the evolution of morality have spurred an astoundingly large interdisciplinary literature. Some significant subset of this body of work addresses questions regarding our moral psychology: how did humans evolve the psychological properties which underpin our systems of ethics and morality? Here I do three things. First, I discuss some methodological issues, and defend particularly effective methods for addressing many research questions in this area. Second, I give an in-depth example, describing how an explanation can be given for the evolution of guilt---one of the core moral emotions---using the methods advocated here. Last, I lay out which sorts of strategic scenarios generally are the ones that our moral psychology evolved to `solve', and thus which models are the most useful in further exploring this evolution.


title: Photonic Higher-Order Topological States Induced by Long Range Interactions

abstract: The discovery of topological phases has recently led to a paradigm shift in condensed matter physics, and facilitated breakthroughs in engineered photonics and acoustic metamaterials. Topological insulators (TIs) enable the generation of electronic, photonic, and acoustic modes exhibiting wave propagation that is resilient to disorder, irrespective of manufacturing precision or unpredictable defects induced by the operational environment, known as topological protection. While originally limited to a dimensionality of the protected states that is one dimension lower than the host TI material, the recent discovery of higher-order topological insulators (HOTIs) provides the potential to overcome this dimensionality limitations by offering topological protection over an extended range of dimensionalities. Here we demonstrate 2D photonic HOTI (PHOTI) with topological states two dimensions lower than the one of the host system. We consider a photonic metacrystal of distorted Kagome lattice geometry that exhibits topological bulk polarization, leading to the emergence of 1D topological edge states and of higher order 0D states confined to the corners of the structure. Interestingly, in addition to corner states due to the nearest neighbour interactions and protected by generalized chiral symmetry 1, we discover and take advantage of a new class of topological corner states sustained by long-range interactions, available in wave-based systems, such as in photonics. Our findings demonstrate that photonic HOTIs possess richer physics compared to their condensed matter counterparts, offering opportunities for engineering novel designer electromagnetic states with unique topological robustness.


title: Polar Codes and Their Quantum-Domain Counterparts

abstract: Arikan's polar codes are capable of achieving the Shannon's capacity at a low encoding and decoding complexity, while inherently supporting rate adaptation. By virtue of these attractive features, polar codes have provided fierce competition to both the turbo as well as the Low Density Parity Check (LDPC) codes, making its way into the 5G New Radio (NR). Realizing the significance of polar codes, in this paper we provide a comprehensive survey of polar codes, highlighting the major milestones achieved in the last decade. Furthermore, we also provide tutorial insights into the polar encoder, decoders as well as the code construction methods. We also extend our discussions to quantum polar codes with an emphasis on the underlying quantum-to-classical isomorphism and the syndrome-based quantum polar codes.


title: A canonical form for pairs consisting of a Hermitian form and a self-adjoint antilinear operator

abstract: Motivated by a problem in local differential geometry of Cauchy--Riemann (CR) structures of hypersurface type, we find a canonical form for pairs consisting of a nondegenerate Hermitian form and a self-adjoint antilinear operator, or, equivalently, consisting of a nondegenerate Hermitian form and a symmetric bilinear form. This generalizes the only previously known results on simultaneous normalization of such pairs, namely, the results of Benedetti and Cragnolini (1984) on simultaneous diagonalization of these pairs in the case where the Hermitian form is positive definite and of Hong and Horn (1986), where a criterion for simultaneous diagonalization is given.


title: Modification of Turbulent Boundary Layer Coherent Structures with Drag Reducing Polymer Solution

abstract: The modification of dominant coherent structures that extend through the log-region of a drag reduced turbulent boundary layer is studied via examination of two-point correlations from time-resolved particle-image-velocimetry. Measurements were acquired in polymer oceans (uniform concentration) at drag reduction levels corresponding to the low drag reduction regime (< 40%), the high drag reduction (HDR) regime (>40%), and at an intermediate level (46%). The mean velocity profiles and two-point correlations were compared with that of water (Newtonian, DR = 0%). These results show that with increasing drag reduction the inclination of these dominant coherent structures decrease, their streamwise extent increases, and the fluctuations in the correlations are suppressed (especially at HDR). These observations are examined in comparison with the coherent structure literature (Newtonian and polymeric).


title: Rate Control for Wireless-Powered Communication Network with Reliability and Delay Constraints

abstract: We consider a two-phase Wireless-Powered Communication Network under Nakagami-m fading, where a wireless energy transfer process first powers a sensor node that then uses such energy to transmit its data in the wireless information transmission phase. We explore a fixed transmit rate scheme designed to cope with the reliability and delay constraints of the system while attaining closed-form approximations for the optimum wireless energy transfer and wireless information transmission blocklength. Then, a more-elaborate rate control strategy exploiting the readily available battery charge information is proposed and the results evidence its outstanding performance when compared with the fixed transmit rate, for which no battery charge information is available. It even reaches an average rate performance close to that of an ideal scheme requiring full Channel State Information at transmitter side. Numerical results show the positive impact of a greater number of antennas at the destination, and evidence that the greater the reliability constraints, the smaller the message sizes on average, and the smaller the optimum information blocklengths. Finally, we corroborate the appropriateness of using the asymptotic blocklength formulation as an approximation of the non-asymptotic finite blocklength results.


title: Mass loss on the red giant branch: plasmoid-driven winds above the RGB bump

abstract: The onset of cool massive winds in evolved giants is correlated with an evolutionary feature on the red giant branch known as the bump. Also at the bump, shear instability in the star leads to magnetic fields that occur preferentially on small length scales. Pneuman (1983) has suggested that the emergence of small scale flux tubes in the Sun can give rise to enhanced acceleration of the solar wind as a result of plasmoid acceleration (the melon seed mechanism). In this paper, we examine the Pneuman formalism to determine if it may shed some light on the process that drives mass loss from stars above the bump. Because we do not currently have detailed information for some of the relevant physical parameters, we are not yet able to derive a detailed model. Instead, our goal in this paper is to explore a proof of concept. Using parameters that are known to be plausible in cool giants, we find that the total mass loss rate from such stars can be replicated. Moreover, we find that the radial profile of the wind speed in such stars can be steep or shallow depending on the fraction of the mass loss which is contained in the plasmoids. This is consistent with empirical data which indicate that the velocity profiles of winds from cool giants range from shallow to steep.


title: Existence of Non-Obvious Divergent Trajectories in homogeneous spaces

abstract: We prove a modified version for a conjecture of Weiss from 2004. Let $G$ be a semisimple real algebraic group defined over $\mathbb{Q}$, $\Gamma$ be an arithmetic subgroup of $G$. A trajectory in $G/\Gamma$ is divergent if eventually it leaves every compact subset, and is obvious divergent if there is a finite collection of algebraic data which cause the divergence. Let $A$ be a diagonalizable subgroup of $G$ of positive dimension. We show that if the projection of $A$ to any $\mathbb{Q}$-factor of $G$ is of small enough dimension (relatively to the $\mathbb{Q}$-rank of the $\mathbb{Q}$-factor), then there are non-obvious divergent trajectories for the action of $A$ on $G/\Gamma$.


title: Periodic Jacobi Operators with Complex Coefficients

abstract: We present certain results on the direct and inverse spectral theory of the Jacobi operator with complex periodic coefficients. For instance, we show that any $N$-th degree polynomial whose leading coefficient is $(-1)^N$ is the Hill discriminant of finitely many discrete $N$-periodic Schr\"{o}dinger operators (Theorem 1). Also, in the case where the spectrum is a closed interval we prove a result (Theorem 5) which is the analog of Borg's Theorem for the non-self-adjoint Jacobi case.


title: An Empirical Mass-Radius Relation for Cool Giant Planets

abstract: Probabilistic relationships between the mass and radius of planets are useful for a variety of purposes, including estimating the yields of planet discovery efforts and the radii of discovered planets given their masses. Previous work on giant planets often do not make the important distinction between cool and (inflated) hot giants. This can lead to needlessly large uncertainties and overestimated radii for the cooler objects. Using hierarchical bayesian modeling and the AIC model selection criterion, we identify an empirical mass-radius relationship for cool giant planets above 15 Earth masses. It takes the form of a quadratic function in log mass, with a power-law dispersion. Parametric uncertainties from MCMC are provided, and we compare our results with the outputs of interior structure models and the radii of the solar system giant planets, which fit comfortably within our relation.


title: Five-dimensional Supergravity in N = 1/2 Superspace

abstract: We construct 5D, N = 1 supergravity in a 4D, N = 1 superspace with an extra bosonic coordinate. This represents four of the supersymmetries and the associated Poincar\'e symmetries manifestly. The remaining four supersymmetries and the rest of the Poincar\'e symmetries are represented linearly but not manifestly. In the linearized approximation, the action reduces to the known superspace result. As an application of the formalism, we construct the $\int A\wedge R\wedge R$ invariant in this superspace.


title: A Human-Centered Data-Driven Planner-Actor-Critic Architecture via Logic Programming

abstract: Recent successes of Reinforcement Learning (RL) allow an agent to learn policies that surpass human experts but suffers from being time-hungry and data-hungry. By contrast, human learning is significantly faster because prior and general knowledge and multiple information resources are utilized. In this paper, we propose a Planner-Actor-Critic architecture for huMAN-centered planning and learning (PACMAN), where an agent uses its prior, high-level, deterministic symbolic knowledge to plan for goal-directed actions, and also integrates the Actor-Critic algorithm of RL to fine-tune its behavior towards both environmental rewards and human feedback. This work is the first unified framework where knowledge-based planning, RL, and human teaching jointly contribute to the policy learning of an agent. Our experiments demonstrate that PACMAN leads to a significant jump-start at the early stage of learning, converges rapidly and with small variance, and is robust to inconsistent, infrequent, and misleading feedback.


title: ELMAG 3.01: A three-dimensional Monte Carlo simulation of electromagnetic cascades on the extragalactic background light and in magnetic fields

abstract: The version 3.01 of ELMAG, a Monte Carlo program for the simulation of electromagnetic cascades initiated by high-energy photons and electrons interacting with extragalactic background light (EBL), is presented. Pair production and inverse Compton scattering on EBL photons as well as synchrotron losses are implemented using weighted sampling of the cascade development. New features include, among others, the implementation of turbulent extragalactic magnetic fields and the calculation of three-dimensional electron and positron trajectories, solving the Lorentz force equation. As final result of the three-dimensional simulations, the program provides two-dimensional source images as function of the energy and the time delay of secondary cascade particles.


title: Border effect corrections for diagonal line based recurrence quantification analysis measures

abstract: Recurrence Quantification Analysis (RQA) defines a number of quantifiers, which base upon diagonal line structures in the recurrence plot (RP). Due to the finite size of an RP, these lines can be cut by the borders of the RP and, thus, bias the length distribution of diagonal lines and, consequently, the line based RQA measures. In this letter we investigate the impact of the mentioned border effects and of the thickening of diagonal lines in an RP (caused by tangential motion) on the estimation of the diagonal line length distribution, quantified by its entropy. Although a relation to the Lyapunov spectrum is theoretically expected, the mentioned entropy yields contradictory results in many studies. Here we summarize correction schemes for both, the border effects and the tangential motion and systematically compare them to methods from the literature. We show that these corrections lead to the expected behavior of the diagonal line length entropy, in particular meaning zero values in case of a regular motion and positive values for chaotic motion. Moreover, we test these methods under noisy conditions, in order to supply practical tools for applied statistical research.


title: RELICS: Spectroscopy of gravitationally-lensed $z\simeq 2$ reionization-era analogs and implications for CIII] detections at $z>6$

abstract: Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$.


title: Experimenting with Constraint Programming on GPU

abstract: The focus of my PhD thesis is on exploring parallel approaches to efficiently solve problems modeled by constraints and presenting a new proposal. Current solvers are very advanced; they are carefully designed to effectively manage the high-level problems' description and include refined strategies to avoid useless work. Despite this, finding a solution can take an unacceptable amount of time. Parallelization can mitigate this problem when the instance of the problem modeled is large, as it happens in real world problems. It is done by propagating constraints in parallel and concurrently exploring different parts of the search space. I am developing on a constraint solver that exploits the many cores available on Graphics Processing Units (GPU) to speed up the search.


title: Asymptotic reduced density matrix of discrete-time quantum walks

abstract: In this article we show that for any quantum walker with \textit{m}-dimensional coin subspace, we have $m^2\times m^2$ specific constant matrix $\mathcal{C}$ where it completely determines the asymptotic reduced density matrix of the walker. We show that for any initial state with $P_0$ projector, reduced density matrix, can be obtained by $Tr_1\left(P_0\otimes I\;\mathcal{C}\right)$ or equivalently $Tr_2\left(I\otimes P_0\;\mathcal{C}\right)$. It is worth to mention that characteristic matrix $\mathcal{C}$ is independent of the initial state and just depends on coin operator, so by finding this matrix for specific type of QW the long-time behavior of it, such as local state of the coin after a long time walking and asymptotic entanglement between coin and position will be completely known for any initial state. We have found the characteristic matrix $\mathcal{C}$ for general coin operator, $U\left(2\right)$, as well as exact form of this matrix for local initial state.


title: Emergent structural length scales in a model binary glass --- the micro-second molecular dynamics time-scale regime

abstract: It is now possible to routinely perform atomistic simulations at the microsecond timescale. In the present work, we exploit this for a model binary Lennard-Jones glass to study structural relaxation at a timescale spanning up to 80 microseconds. It is found that at these longer time-scales, significant mobility and relaxation occurs, demonstrating a reproducible relaxation trajectory towards a class of asymptotic structural states which are strongly heterogeneous. These structures are discussed both in terms of the ideal glass state and a production route towards novel amorphous crystalline nano-composite micro-structures with strong interface stability.


title: Control landscape for ultrafast manipulation by a qubit

abstract: In this work we study extrema of objective functionals for ultrafast manipulation by a qubit. Traps are extrema of the objective functionals which are optimal for manipulation by quantum systems only locally but not globally. Much effort in prior works was devoted to the analysis of traps for quantum systems controlled by long enough laser pulses and, for example, manipulation by a qubit with long control pulses was shown to be trap-free. Ultrafast femtosecond and attosecond control becomes now widely applicable that motivates the necessity for the analysis of traps on the ultrafast time scale. We do such analysis for a qubit and show that ultrafast state transfer in a qubit remains trap-free for a wide range of the initial and final states of the qubit. We prove that for this range the probability of transition between the initial and the final states has a saddle but has no traps.


title: Structured Discrete Shape Approximation: Theoretical Complexity and Practical Algorithm

abstract: We consider the problem of approximating a two-dimensional shape contour (or curve segment) using discrete assembly systems, which allow to build geometric structures based on limited sets of node and edge types subject to edge length and orientation restrictions. We show that already deciding feasibility of such approximation problems is NP-hard, and remains intractable even for very simple setups. We then devise an algorithmic framework that combines shape sampling with exact cardinality-minimization to obtain good approximations using few components. As a particular application and showcase example, we discuss approximating shape contours using the classical Zometool construction kit and provide promising computational results, demonstrating that our algorithm is capable of obtaining good shape representations within reasonable time, in spite of the problem's general intractability. We conclude the paper with an outlook on possible extensions of the developed methodology, in particular regarding 3D shape approximation tasks.


title: Interpretable Discriminative Dimensionality Reduction and Feature Selection on the Manifold

abstract: Dimensionality reduction (DR) on the manifold includes effective methods which project the data from an implicit relational space onto a vectorial space. Regardless of the achievements in this area, these algorithms suffer from the lack of interpretation of the projection dimensions. Therefore, it is often difficult to explain the physical meaning behind the embedding dimensions. In this research, we propose the interpretable kernel DR algorithm (I-KDR) as a new algorithm which maps the data from the feature space to a lower dimensional space where the classes are more condensed with less overlapping. Besides, the algorithm creates the dimensions upon local contributions of the data samples, which makes it easier to interpret them by class labels. Additionally, we efficiently fuse the DR with feature selection task to select the most relevant features of the original space to the discriminative objective. Based on the empirical evidence, I-KDR provides better interpretations for embedding dimensions as well as higher discriminative performance in the embedded space compared to the state-of-the-art and popular DR algorithms.


title: Doppler tomography as a tool for detecting exoplanet atmospheres

abstract: High-resolution Doppler spectroscopy is a powerful tool for identifying molecular species in the atmospheres of both transiting and non-transiting exoplanets. Currently, such data is analysed using cross-correlation techniques to detect the Doppler shifting signal from the orbiting planet. In this paper we demonstrate that, compared to cross-correlation methods currently used, the technique of Doppler tomography has improved sensitivity in detecting the subtle signatures expected from exoplanet atmospheres. This is partly due to the use of a regularizing statistic, which acts to suppress noise, coupled to the fact that all the data is fit simultaneously. In addition, we show that the technique can also effectively suppress contanimating spectral features that may arise due to overlapping lines, repeating line patterns, or the use of incorrect linelists. These issues can confuse conventional cross-correlation approaches, primarily due to aliasing issues inherent in such techniques, whereas Doppler tomography is less susceptible to such effects. In particular, Doppler tomography shows exceptional promise for simultaneously detecting multiple line species (e.g. isotopologues), even when there are high contrasts between such species -- and far outperforms current CCF analyses in this respect. Finally, we demonstrate that Doppler tomography is capable of recovering molecular signals from exoplanets using real data, by confirming the strong detection of CO in the atmosphere of Tau Boo b. We recover a signal with a planetary radial velocity semi-amplitude Kp = 109.6 +/- 2.2 km/s, in excellent agreement with the previously reported value of 110.0 +/- 3.2 km/s.


title: Goal-Embedded Dual Hierarchical Model for Task-Oriented Dialogue Generation

abstract: Hierarchical neural networks are often used to model inherent structures within dialogues. For goal-oriented dialogues, these models miss a mechanism adhering to the goals and neglect the distinct conversational patterns between two interlocutors. In this work, we propose Goal-Embedded Dual Hierarchical Attentional Encoder-Decoder (G-DuHA) able to center around goals and capture interlocutor-level disparity while modeling goal-oriented dialogues. Experiments on dialogue generation, response generation, and human evaluations demonstrate that the proposed model successfully generates higher-quality, more diverse and goal-centric dialogues. Moreover, we apply data augmentation via goal-oriented dialogue generation for task-oriented dialog systems with better performance achieved.


title: Berezin regularity of domains in C^n and the essential norms of Toeplitz operators

abstract: For the open unit disc $\mathbb{D}$ in the complex plane, it is well known that if $\phi \in C(\overline{\mathbb{D}})$ then its Berezin transform $\widetilde{\phi}$ also belongs to $C(\overline{\mathbb{D}})$. We say that $\mathbb{D}$ is BC-regular. In this paper we study BC-regularity of some pseudoconvex domains in $\mathbb{C}^n$ and show that the boundary geometry plays an important role. We also establish a relationship between the essential norm of an operator in a natural Toeplitz subalgebra and its Berezin transform.


title: Large- and Medium-Scale Anisotropies in the Arrival Directions of Cosmic Rays observed with KASCADE-Grande

abstract: We search for anisotropies in the arrival directions of cosmic rays observed by the KASCADE-Grande air shower experiment. The analysis is based on public data of about 23.7 million events with reconstructed primary energies above 1 PeV. We apply a novel maximum-likelihood reconstruction method for the cosmic ray anisotropy, that compensates for spurious anisotropies induced by local detector effects. We find no evidence for a large-scale dipole anisotropy in the data, consistent with official results based on the conventional East-West derivative method. On the other hand, a subset of cosmic rays with median energy of 33 PeV shows strong evidence for a medium-scale feature with an angular diameter of 40 degrees. After accounting for the look-elsewhere effect, the post-trial significance of this medium-scale feature is at the level of 4$\sigma$.


title: InterpretML: A Unified Framework for Machine Learning Interpretability

abstract: InterpretML is an open-source Python package which exposes machine learning interpretability algorithms to practitioners and researchers. InterpretML exposes two types of interpretability - glassbox models, which are machine learning models designed for interpretability (ex: linear models, rule lists, generalized additive models), and blackbox explainability techniques for explaining existing systems (ex: Partial Dependence, LIME). The package enables practitioners to easily compare interpretability algorithms by exposing multiple methods under a unified API, and by having a built-in, extensible visualization platform. InterpretML also includes the first implementation of the Explainable Boosting Machine, a powerful, interpretable, glassbox model that can be as accurate as many blackbox models. The MIT licensed source code can be downloaded from github.com/microsoft/interpret.


title: The Colliding Reciprocal Dance Problem: A Mitigation Strategy with Application to Automotive Active Safety Systems

abstract: A reciprocal dance occurs when two mobile agents attempt to pass each other but incompatible interaction models result in repeated attempts to take mutually blocking actions. Often, such a situation simply results in deadlock. But in systems with significant inertial constraints, it can result in collision. This abstract presents this colliding variant of the reciprocal dance, how it arises, and a mitigation strategy that can improve safety without sacrificing flexibility. A demonstration of the concept is provided in the context of automotive active safety.


title: Gaze Estimation for Assisted Living Environments

abstract: Effective assisted living environments must be able to perform inferences on how their occupants interact with one another as well as with surrounding objects. To accomplish this goal using a vision-based automated approach, multiple tasks such as pose estimation, object segmentation and gaze estimation must be addressed. Gaze direction in particular provides some of the strongest indications of how a person interacts with the environment. In this paper, we propose a simple neural network regressor that estimates the gaze direction of individuals in a multi-camera assisted living scenario, relying only on the relative positions of facial keypoints collected from a single pose estimation model. To handle cases of keypoint occlusion, our model exploits a novel confidence gated unit in its input layer. In addition to the gaze direction, our model also outputs an estimation of its own prediction uncertainty. Experimental results on a public benchmark demonstrate that our approach performs on pair with a complex, dataset-specific baseline, while its uncertainty predictions are highly correlated to the actual angular error of corresponding estimations. Finally, experiments on images from a real assisted living environment demonstrate the higher suitability of our model for its final application.


title: Search for direct production of electroweakinos in final states with one lepton, missing transverse momentum and a Higgs boson decaying into two $b$-jets in $pp$ collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

abstract: The results of a search for electroweakino pair production $pp \rightarrow \tilde\chi^\pm_1 \tilde\chi^0_2$ in which the chargino ($\tilde\chi^\pm_1$) decays into a $W$ boson and the lightest neutralino ($\tilde\chi^0_1$), while the heavier neutralino ($\tilde\chi^0_2$) decays into the Standard Model 125 GeV Higgs boson and a second $\tilde\chi^0_1$ are presented. The signal selection requires a pair of $b$-tagged jets consistent with those from a Higgs boson decay, and either an electron or a muon from the $W$ boson decay, together with missing transverse momentum from the corresponding neutrino and the stable neutralinos. The analysis is based on data corresponding to 139 $\mathrm{fb}^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collisions provided by the Large Hadron Collider and recorded by the ATLAS detector. No statistically significant evidence of an excess of events above the Standard Model expectation is found. Limits are set on the direct production of the electroweakinos in simplified models, assuming pure wino cross-sections. Masses of $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ up to 740 GeV are excluded at 95% confidence level for a massless $\tilde{\chi}^{0}_{1}$.


title: An Introduction to Quaternion-Valued Recurrent Projection Neural Networks

abstract: Hypercomplex-valued neural networks, including quaternion-valued neural networks, can treat multi-dimensional data as a single entity. In this paper, we introduce the quaternion-valued recurrent projection neural networks (QRPNNs). Briefly, QRPNNs are obtained by combining the non-local projection learning with the quaternion-valued recurrent correlation neural network (QRCNNs). We show that QRPNNs overcome the cross-talk problem of QRCNNs. Thus, they are appropriate to implement associative memories. Furthermore, computational experiments reveal that QRPNNs exhibit greater storage capacity and noise tolerance than their corresponding QRCNNs.


title: Non-Bayesian Social Learning with Uncertain Models

abstract: Non-Bayesian social learning theory provides a framework that models distributed inference for a group of agents interacting over a social network. In this framework, each agent iteratively forms and communicates beliefs about an unknown state of the world with their neighbors using a learning rule. Existing approaches assume agents have access to precise statistical models (in the form of likelihoods) for the state of the world. However in many situations, such models must be learned from finite data. We propose a social learning rule that takes into account uncertainty in the statistical models using second-order probabilities. Therefore, beliefs derived from uncertain models are sensitive to the amount of past evidence collected for each hypothesis. We characterize how well the hypotheses can be tested on a social network, as consistent or not with the state of the world. We explicitly show the dependency of the generated beliefs with respect to the amount of prior evidence. Moreover, as the amount of prior evidence goes to infinity, learning occurs and is consistent with traditional social learning theory.


title: On the Mathematical Foundations of Causal Fermion Systems in Minkowski Space

abstract: The emergence of the concept of a causal fermion system is revisited and further investigated for the vacuum Dirac equation in Minkowski space. After a brief recap of the Dirac equation and its solution space, in order to allow for the effects of a possibly nonstandard structure of spacetime at the Planck scale, a regularization by a smooth cutoff in momentum space is introduced, and its properties are discussed. Given an ensemble of solutions, we recall the construction of a local correlation function, which realizes spacetime in terms of operators. It is shown in various situations that the local correlation function maps spacetime points to operators of maximal rank and that it is closed and homeomorphic onto its image. It is inferred that the corresponding causal fermion systems are regular and have a smooth manifold structure. The cases considered include a Dirac sea vacuum and systems involving a finite number of particles and antiparticles.


title: Determining the Composition of Relativistic Jets from Polarization Maps

abstract: We present a stationary, axisymmetric, self-similar semi-analytic model of magnetically dominated jet plasma based on force-free regions of a relativistic magnetohydrodynamic simulation. We use this model to illustrate how the composition of relativistic jet plasma can be determined, with special attention to the example of M87. In particular, we compute synthetic Stokes maps in e-e+p plasmas with various positron-to-proton ratios using synchrotron emission models scaling the partial pressure of electrons and positrons emitting at the observed frequency to the magnetic pressure, taking into account Faraday rotation and conversion. The lepton-dominated models produce bilaterally asymmetric radio intensity profiles with strong linear polarization and Stokes Q and U maps that are bilaterally asymmetric (but strongly up-down correlated) and antisymmetric (and sometimes up-down anticorrelated), respectively. The hadronic models produce more centrally brightened intensity and polarization maps. Circular polarization provides the cleanest observational tool for distinguishing the plasmas, as it increases outward from the jet core and central axis for highly ionic plasma, and is suppressed for pair dominated plasma. We find a measurable degree of circular polarization V/I of O(10e-3) for sub-equipartition hadronic jet plasmas. Our stationary model predicts that the intensity-normalized autocorrelation functions of Q and U increase and decrease with frequency, respectively. On the other hand, the autocorrelation of V is less sensitive to the frequency. Multi-band polarimetric observations could therefore be used as a novel probe of the composition of jet plasma.


title: Chaitin's Omega and an Algorithmic Phase Transition

abstract: We consider the statistical mechanical ensemble of bit string histories that are computed by a universal Turing machine. The role of the energy is played by the program size. We show that this ensemble has a first-order phase transition at a critical temperature, at which the partition function equals Chaitin's halting probability $\Omega$. This phase transition has curious properties: the free energy is continuous near the critical temperature, but almost jumps: it converges more slowly to its finite critical value than any computable function. At the critical temperature, the average size of the bit strings diverges. We define a non-universal Turing machine that approximates this behavior of the partition function in a computable way by a super-logarithmic singularity, and discuss its thermodynamic properties. We also discuss analogies and differences between Chaitin's Omega and the partition functions of a quantum mechanical particle, a spin model, random surfaces, and quantum Turing machines. For universal Turing machines, we conjecture that the ensemble of bit string histories at the critical temperature has a continuum formulation in terms of a string theory.

categories: cs.CC cs.IT hep-th math-ph math.IT math.LO math.MP


title: The arithmetic Kuznetsov formula on $GL(3)$, II: The general case

abstract: We obtain the last of the standard Kuznetsov formulas for $SL(3,\Bbb{Z})$. In the previous paper, we were able to exploit the relationship between the positive-sign Bessel function and the Whittaker function to apply Wallach's Whittaker expansion; now we demonstrate the expansion of functions into Bessel functions for all four signs, generalizing Wallach's theorem for $SL(3)$. As applications, we again consider the Kloosterman zeta functions and smooth sums of Kloosterman sums. The new Kloosterman zeta functions pose the same difficulties as we saw with the positive-sign case, but for the negative-sign case, we obtain some analytic continuation of the unweighted zeta function and give a sort of reflection formula that exactly demonstrates the obstruction when the moduli are far apart. The completion of the remaining sign cases means this work now both supersedes the author's thesis and completes the work started in the original paper of Bump, Friedberg and Goldfeld.


title: Robust Humanoid Contact Planning with Learned Zero- and One-Step Capturability Prediction

abstract: Humanoid robots maintain balance and navigate by controlling the contact wrenches applied to the environment. While it is possible to plan dynamically-feasible motion that applies appropriate wrenches using existing methods, a humanoid may also be affected by external disturbances. Existing systems typically rely on controllers to reactively recover from disturbances. However, such controllers may fail when the robot cannot reach contacts capable of rejecting a given disturbance. In this paper, we propose a search-based footstep planner which aims to maximize the probability of the robot successfully reaching the goal without falling as a result of a disturbance. The planner considers not only the poses of the planned contact sequence, but also alternative contacts near the planned contact sequence that can be used to recover from external disturbances. Although this additional consideration significantly increases the computation load, we train neural networks to efficiently predict multi-contact zero-step and one-step capturability, which allows the planner to generate robust contact sequences efficiently. Our results show that our approach generates footstep sequences that are more robust to external disturbances than a conventional footstep planner in four challenging scenarios.


title: Dark CP-violation through the $Z$-portal

abstract: Despite great agreement with experiment, the Standard Model (SM) of particle physics lacks a viable Dark Matter (DM) candidate and sufficient amount of CP-violation to account for the observed baryon excess in the universe. Non-minimal Higgs frameworks are economic extensions of the SM which could remedy these shortcomings. Within the framework of a three Higgs doublet model, we introduce an extended dark sector which accommodates both DM and CP-violation. Such dark sources of CP-violation do not contribute to the Electric Dipole Moments and are therefore unconstrained. We present a novel mechanism in which the CP-violating dark particles only interact with the SM through the gauge bosons, primarily the $Z$ boson. Such $Z$-portal dark CP violation is realised in the regions of the parameter space where Higgs-mediated (co)annihilation processes are sub-dominant and have negligible contributions to the DM relic density. We show that such $Z$-portal CP violating DM can still thermalise and satisfy all experimental and observational bounds and discuss the implications of such phenomena for electroweak baryogenesis.


title: On the Impact of Ground Sound

abstract: Rigid-body impact sound synthesis methods often omit the ground sound. In this paper we analyze an idealized ground-sound model based on an elastodynamic halfspace, and use it to identify scenarios wherein ground sound is perceptually relevant versus when it is masked by the impacting object's modal sound or transient acceleration noise. Our analytical model gives a smooth, closed-form expression for ground surface acceleration, which we can then use in the Rayleigh integral or in an "acoustic shader" for a finite-difference time-domain wave simulation. We find that when modal sound is inaudible, ground sound is audible in scenarios where a dense object impacts a soft ground and scenarios where the impact point has a low elevation angle to the listening point.


title: The influence of cut vertices and eigenvalues on character graphs of solvable groups

abstract: Given a finite group $G$, the character graph, denoted by $\Delta(G)$, for its irreducible character degrees is a graph with vertex set $\rho(G)$ which is the set of prime numbers that divide the irreducible character degrees of $G$, and with $\{p,q\}$ being an edge if there exist a non-linear $\chi\in {\rm Irr}(G)$ whose degree is divisible by $pq$. In this paper, we discuss the influences of cut vertices and eigenvalues of $\Delta(G)$ on the group structure of $G$. Recently, Lewis and Meng proved the character graph of each solvable group has at most one cut vertex. Now, we determine the structure of character graphs of solvable groups with a cut vertex and diameter $3$. Furthermore, we study solvable groups whose character graphs have at most two distinct eigenvalues. Moreover, we investigate the solvable groups whose character graphs are regular with three distinct eigenvalues. In addition, we give some lower bounds for the number of edges of $\Delta(G)$.


title: Improved Variational Neural Machine Translation by Promoting Mutual Information

abstract: Posterior collapse plagues VAEs for text, especially for conditional text generation with strong autoregressive decoders. In this work, we address this problem in variational neural machine translation by explicitly promoting mutual information between the latent variables and the data. Our model extends the conditional variational autoencoder (CVAE) with two new ingredients: first, we propose a modified evidence lower bound (ELBO) objective which explicitly promotes mutual information; second, we regularize the probabilities of the decoder by mixing an auxiliary factorized distribution which is directly predicted by the latent variables. We present empirical results on the Transformer architecture and show the proposed model effectively addressed posterior collapse: latent variables are no longer ignored in the presence of powerful decoder. As a result, the proposed model yields improved translation quality while demonstrating superior performance in terms of data efficiency and robustness.


title: A necessary and sufficient condition for radial property of positive entire solutions of $\Delta^2 u+u^{-q}=0$ in $\mathbf{R}^3$

abstract: In this article, we are concerned with the following geometric equation \begin{equation}\label{MainEq} \Delta^2 u = -u^{-q} \qquad \text{in } \mathbf{R}^3 \end{equation} for $q>0$. Recently in \cite{GWZ18}, Guo, Wei and Zhou have established the relationship between the radial symmetry and the exact growth rate at infinity of a positive entire solution of that equation as $1<q<3$. The aim of this paper is to obtain the similar result in the case $q>3$ thanks to the method of moving plane.


title: Systemic Cascades On Inhomogeneous Random Financial Networks

abstract: This systemic risk paper introduces inhomogeneous random financial networks (IRFNs). Such models are intended to describe parts, or the entirety, of a highly heterogeneous network of banks and their interconnections, in the global financial system. Both the balance sheets and the stylized crisis behaviour of banks are ingredients of the network model. A systemic crisis is pictured as triggered by a shock to banks' balance sheets, which then leads to the propagation of damaging shocks and the potential for amplification of the crisis, ending with the system in a cascade equilibrium. Under some conditions the model has ``locally tree-like independence (LTI)'', where a general percolation theoretic argument leads to an analytic fixed point equation describing the cascade equilibrium when the number of banks $N$ in the system is taken to infinity. This paper focusses on mathematical properties of the framework in the context of Eisenberg-Noe solvency cascades generalized to account for fractional bankruptcy charges. New results including a definition and proof of the ``LTI property'' of the Eisenberg-Noe solvency cascade mechanism lead to explicit $N=\infty$ fixed point equations that arise under very general model specifications. The essential formulas are shown to be implementable via well-defined approximation schemes, but numerical exploration of some of the wide range of potential applications of the method is left for future work.


title: Experimental Demonstration of Probabilistic Spin Logic by Magnetic Tunnel Junctions

abstract: The recently proposed probabilistic spin logic presents promising solutions to novel computing applications. Multiple cases of implementations, including invertible logic gate, have been studied numerically by simulations. Here we report an experimental demonstration of a magnetic tunnel junction-based hardware implementation of probabilistic spin logic.


title: Nguyen's Tridents and the Classification of Semigraphical Translators for Mean Curvature Flow

abstract: We construct a one-parameter family of singly periodic translating solutions to mean curvature flow that converge as the period tends to $0$ to the union of a grim reaper surface and a plane that bisects it lengthwise. The surfaces are semigraphical: they are properly embedded, and, after removing a discrete collection of vertical lines, they are graphs. We also provide a nearly complete classification of semigraphical translators.


title: How large a union of proper subgroups of a finite group can be?

abstract: Let $k$ be a positive integer and $G$ be a finite group that cannot be written as the union of $k$ proper subgroups. In this short note, we study the existence of a constant $c_k\in (0,1)$ such that $|\cup_{i=1}^k H_i| \leq c_k|G|$, for all proper subgroups $H_1$, ..., $H_k$ of $G$.


title: Generalized numerical radius and related inequalities

abstract: They proved several properties and introduced some inequalities. We continue with the study of this generalized numerical radius and we develop diverse inequalities involving w_N. We also study particular cases with a fixed N(.), for instance the p-Schatten norms. In ["A generalization of the numerical radius". Linear Algebra Appl. 569 (2019)], Abu Omar and Kittaneh defined a new generalization of the numerical radius. That is, given a norm $N(\cdot)$ on $\bh$, the space of bounded linear operators over a Hilbert space H, and A in B(H) w_N(A)=sup_{\theta\in \R}N(Re(e^{i\theta}A)). They proved several properties and introduced some inequalities. We continue with the study of this generalized numerical radius and we develop diverse inequalities involving $w_N$. We also study particular cases when N(.) is the p- Schatten norm with p>1.


title: Inversion of lattice models from the observations of microscopic degrees of freedom: parameter estimation with uncertainty quantification

abstract: Experimental advances in condensed matter physics and material science have enabled ready access to atomic-resolution images, with resolution of modern tools often sufficient to extract minute details of symmetry-breaking distortions such as polarization, octahedra tilts, or other structure-coupled order parameters. The patterns of observed distortions in turn contain the information on microscopic driving forces defining the development of materials microstructure and associated thermodynamics. However, the analysis of underpinning physical models from experimentally observed microscopic degrees of freedom remains a largely unresolved issue. Here, we explore such an approach using the paradigmatic Ising model on a square lattice. We show that the microscopic parameters of the Ising model both for ferromagnetic and antiferromagnetic case can be extracted from the spin configurations for temperatures an order of magnitude higher than the phase transition and perform uncertainty analysis for such reconstructions. This suggests that microscopic observations of materials with sufficiently high precision can provide information on generative physics at temperatures well above corresponding phase transition, opening new horizons for scientific exploration via high-resolution imaging.


title: Annular Rasmussen invariants: Properties and 3-braid classification

abstract: We prove that for a fixed braid index there are only finitely many possible shapes of the annular Rasmussen $d_t$ invariant of braid closures. Applying the same perspective to the knot Floer invariant $\Upsilon_K(t)$, we show that for a fixed concordance genus of $K$ there are only finitely many possibilities for $\Upsilon_K(t)$. Focusing on the case of 3-braids, we compute the Rasmussen $s$ invariant and the annular Rasmussen $d_t$ invariant of all 3-braid closures. As a corollary, we show that the vanishing/non-vanishing of the $\psi$ invariant is entirely determined by the $s$ invariant and the self-linking number.


title: Machine Learning for Clinical Predictive Analytics

abstract: In this chapter, we provide a brief overview of applying machine learning techniques for clinical prediction tasks. We begin with a quick introduction to the concepts of machine learning and outline some of the most common machine learning algorithms. Next, we demonstrate how to apply the algorithms with appropriate toolkits to conduct machine learning experiments for clinical prediction tasks. The objectives of this chapter are to (1) understand the basics of machine learning techniques and the reasons behind why they are useful for solving clinical prediction problems, (2) understand the intuition behind some machine learning models, including regression, decision trees, and support vector machines, and (3) understand how to apply these models to clinical prediction problems using publicly available datasets via case studies.


title: Mobility edge of the two dimensional Bose-Hubbard model

abstract: We analyze the disorder driven localization of the two dimensional Bose-Hubbard model by evaluating the full low energy quasiparticle spectrum via a recently developed fluctuation operator expansion method. For any strength of the local interaction we find a mobility edge that displays an approximately exponential decay with increasing disorder strength. We determine the finite-size scaling collapse and exponents at this critical line finding that the localization of excitations is characterized by weak multi-fractality and a thermal-like critical gap ratio. A direct comparison to a recent experiment yields an excellent match of the predicted finite-size transition point and scaling of single particle correlations.


title: Representation Learning for Electronic Health Records

abstract: Information in electronic health records (EHR), such as clinical narratives, examination reports, lab measurements, demographics, and other patient encounter entries, can be transformed into appropriate data representations that can be used for downstream clinical machine learning tasks using representation learning. Learning better representations is critical to improve the performance of downstream tasks. Due to the advances in machine learning, we now can learn better and meaningful representations from EHR through disentangling the underlying factors inside data and distilling large amounts of information and knowledge from heterogeneous EHR sources. In this chapter, we first introduce the background of learning representations and reasons why we need good EHR representations in machine learning for medicine and healthcare in Section 1. Next, we explain the commonly-used machine learning and evaluation methods for representation learning using a deep learning approach in Section 2. Following that, we review recent related studies of learning patient state representation from EHR for clinical machine learning tasks in Section 3. Finally, in Section 4 we discuss more techniques, studies, and challenges for learning natural language representations when free texts, such as clinical notes, examination reports, or biomedical literature are used. We also discuss challenges and opportunities in these rapidly growing research fields.


title: A consensus-based global optimization method for high dimensional machine learning problems

abstract: We improve recently introduced consensus-based optimization method, proposed in [R. Pinnau, C. Totzeck, O. Tse and S. Martin, Math. Models Methods Appl. Sci., 27(01):183--204, 2017], which is a gradient-free optimization method for general non-convex functions. We first replace the isotropic geometric Brownian motion by the component-wise one, thus removing the dimensionality dependence of the drift rate, making the method more competitive for high dimensional optimization problems. Secondly, we utilize the random mini-batch ideas to reduce the computational cost of calculating the weighted average which the individual particles tend to relax toward. For its mean-field limit--a nonlinear Fokker-Planck equation--we prove, in both time continuous and semi-discrete settings, that the convergence of the method, which is exponential in time, is guaranteed with parameter constraints {\it independent} of the dimensionality. We also conduct numerical tests to high dimensional problems to check the success rate of the method.


title: Distribution function of the blow up time of the solution of an anticipating random fatigue equation

abstract: In this paper, we study the distribution function of the time of explosion of a stochastic differential equation modeling the length of the dominant crack due to fatigue. The main novelty is that initial condition is regarded as an anticipating random variable and the stochastic integral is in the forward sense. Under suitable conditions, we use the substitution formula from Russo and Vallois to find the local solution of this equation. Then, we find the law of blow up time by proving some results on barrier crossing probabilities of Brownian bridge.


title: AllenNLP Interpret: A Framework for Explaining Predictions of NLP Models

abstract: Neural NLP models are increasingly accurate but are imperfect and opaque---they break in counterintuitive ways and leave end users puzzled at their behavior. Model interpretation methods ameliorate this opacity by providing explanations for specific model predictions. Unfortunately, existing interpretation codebases make it difficult to apply these methods to new models and tasks, which hinders adoption for practitioners and burdens interpretability researchers. We introduce AllenNLP Interpret, a flexible framework for interpreting NLP models. The toolkit provides interpretation primitives (e.g., input gradients) for any AllenNLP model and task, a suite of built-in interpretation methods, and a library of front-end visualization components. We demonstrate the toolkit's flexibility and utility by implementing live demos for five interpretation methods (e.g., saliency maps and adversarial attacks) on a variety of models and tasks (e.g., masked language modeling using BERT and reading comprehension using BiDAF). These demos, alongside our code and tutorials, are available at https://allennlp.org/interpret .


title: HyperLearn: A Distributed Approach for Representation Learning in Datasets With Many Modalities

abstract: Multimodal datasets contain an enormous amount of relational information, which grows exponentially with the introduction of new modalities. Learning representations in such a scenario is inherently complex due to the presence of multiple heterogeneous information channels. These channels can encode both (a) inter-relations between the items of different modalities and (b) intra-relations between the items of the same modality. Encoding multimedia items into a continuous low-dimensional semantic space such that both types of relations are captured and preserved is extremely challenging, especially if the goal is a unified end-to-end learning framework. The two key challenges that need to be addressed are: 1) the framework must be able to merge complex intra and inter relations without losing any valuable information and 2) the learning model should be invariant to the addition of new and potentially very different modalities. In this paper, we propose a flexible framework which can scale to data streams from many modalities. To that end we introduce a hypergraph-based model for data representation and deploy Graph Convolutional Networks to fuse relational information within and across modalities. Our approach provides an efficient solution for distributing otherwise extremely computationally expensive or even unfeasible training processes across multiple-GPUs, without any sacrifices in accuracy. Moreover, adding new modalities to our model requires only an additional GPU unit keeping the computational time unchanged, which brings representation learning to truly multimodal datasets. We demonstrate the feasibility of our approach in the experiments on multimedia datasets featuring second, third and fourth order relations.


title: What's Missing: A Knowledge Gap Guided Approach for Multi-hop Question Answering

abstract: Multi-hop textual question answering requires combining information from multiple sentences. We focus on a natural setting where, unlike typical reading comprehension, only partial information is provided with each question. The model must retrieve and use additional knowledge to correctly answer the question. To tackle this challenge, we develop a novel approach that explicitly identifies the knowledge gap between a key span in the provided knowledge and the answer choices. The model, GapQA, learns to fill this gap by determining the relationship between the span and an answer choice, based on retrieved knowledge targeting this gap. We propose jointly training a model to simultaneously fill this knowledge gap and compose it with the provided partial knowledge. On the OpenBookQA dataset, given partial knowledge, explicitly identifying what's missing substantially outperforms previous approaches.


title: Microwave filtering using forward Brillouin scattering in photonic-phononic emit-receive devices

abstract: Microwave photonic systems are compelling for their ability to process signals at high frequencies and over extremely wide bandwidths as a basis for next generation communication and radar technologies. However, many applications also require narrow-band $(\sim\text{MHz})$ filtering operations that are challenging to implement using optical filtering techniques, as this requires reliable integration of ultra-high quality factor $(\sim 10^8)$ optical resonators. One way to address this challenge is to utilize long-lived acoustic resonances, taking advantage of their narrow-band frequency response to filter microwave signals. In this paper, we examine new strategies to harness a narrow-band acoustic response within a microwave-photonic signal processing platform through use of light-sound coupling. Our signal processing scheme is based on a recently demonstrated photon-phonon emitter-receiver device, which transfers information between the optical and acoustic domains using Brillouin interactions, and produces narrow-band filtering of a microwave signal. To understand the best way to use this device technology, we study the properties of a microwave-photonic link using this filtering scheme. We theoretically analyze the noise characteristics of this microwave-photonic link, and explore the parameter space for the design and optimization of such systems.


title: CP violation for four generations of quarks

abstract: We discuss the generalization of the Jarlskog condition of CP conservation for the case of the Standard Model with 4~quark generations. We express this condition in terms of the 3 Jarlskog invariants of the CKM matrix. Next we present the test for the existence of the 4-th quark generation in terms of the Jarlskog invariants involving only known particles.


title: Triplet-Aware Scene Graph Embeddings

abstract: Scene graphs have become an important form of structured knowledge for tasks such as for image generation, visual relation detection, visual question answering, and image retrieval. While visualizing and interpreting word embeddings is well understood, scene graph embeddings have not been fully explored. In this work, we train scene graph embeddings in a layout generation task with different forms of supervision, specifically introducing triplet super-vision and data augmentation. We see a significant performance increase in both metrics that measure the goodness of layout prediction, mean intersection-over-union (mIoU)(52.3% vs. 49.2%) and relation score (61.7% vs. 54.1%),after the addition of triplet supervision and data augmentation. To understand how these different methods affect the scene graph representation, we apply several new visualization and evaluation methods to explore the evolution of the scene graph embedding. We find that triplet supervision significantly improves the embedding separability, which is highly correlated with the performance of the layout prediction model.


title: How to design a derivatives market?

abstract: We consider the problem of designing a derivatives exchange aiming at addressing clients needs in terms of listed options and providing suitable liquidity. We proceed into two steps. First we use a quantization method to select the options that should be displayed by the exchange. Then, using a principal-agent approach, we design a make take fees contract between the exchange and the market maker. The role of this contract is to provide incentives to the market maker so that he offers small spreads for the whole range of listed options, hence attracting transactions and meeting the commercial requirements of the exchange.


title: Why phonon scattering in glasses is universally small at low temperatures

abstract: We present a novel view of the standard model of tunneling two level systems (TLS) to explain the puzzling universal value of a quantity, $C\sim 3\times 10^{-4}$, that characterizes phonon scattering in glasses below 1 K as reflected in thermal conductivity, ultrasonic attenuation, internal friction, and the change in sound velocity. Physical considerations lead to a broad distribution of phonon-TLS couplings that (1) exponentially renormalize tunneling matrix elements, and (2) reduce the TLS density of states through TLS-TLS interactions. We find good agreement between theory and experiment for a variety of individual glasses.


title: How to put quantum particles on magic bullet trajectories that can hit two targets without a clear line-of-sight

abstract: Quantum particles move in strange ways, even when they propagate freely in space. As a result of the uncertainty principle, it is not possible to control the initial conditions of particle emission in such a way that the particle will definitely pass through two precisely defined positions along its path, even if it is possible to line up the two positions with the emitter. However, there is also an upside to the quantum mechanical laws of motion: constructive quantum interferences can actually raise probabilities to values higher than those permitted by classical causality. Here, it is shown that conventional interferometric methods can be used to prepare photons in a quantum state in which a non-vanishing fraction of particles will hit both of two possible targets, even though the direct line-of-sight connecting the two targets is blocked at the source. The demonstration of the effect is complicated by the uncertainty principle because the physical detection of a particle at one target disturbs the motion of the particle, making it impossible to determine whether the initial state of motion would have allowed the particle to hit the other target or not. It is nonetheless possible to determine the minimal fraction of "magic bullet" particles that must have hit both targets by showing that the number of particles hitting target A is larger than the number of particles missing target B. Quantum interference effects can thus be used to optimize the path of particles in free space beyond the classical limit of motion along a straight line.


title: Generating few-cycle pulses with integrated nonlinear photonics

abstract: Ultrashort laser pulses that last only a few optical cycles have been transformative tools for studying and manipulating light--matter interactions. Few-cycle pulses are typically produced from high-peak-power lasers, either directly from a laser oscillator, or through nonlinear effects in bulk or fiber materials. Now, an opportunity exists to explore the few-cycle regime with the emergence of fully integrated nonlinear photonics. Here, we experimentally and numerically demonstrate how lithographically patterned waveguides can be used to generate few-cycle laser pulses from an input seed pulse. Moreover, our work explores a design principle in which lithographically varying the group-velocity dispersion in a waveguide enables the creation of highly constant-intensity supercontinuum spectra across an octave of bandwidth. An integrated source of few-cycle pulses could broaden the range of applications for ultrafast light sources, including supporting new lab-on-a-chip systems in a scalable form factor.


title: Posterior Contraction Rate of Sparse Latent Feature Models with Application to Proteomics

abstract: The Indian buffet process (IBP) and phylogenetic Indian buffet process (pIBP) can be used as prior models to infer latent features in a data set. The theoretical properties of these models are under-explored, however, especially in high dimensional settings. In this paper, we show that under mild sparsity condition, the posterior distribution of the latent feature matrix, generated via IBP or pIBP priors, converges to the true latent feature matrix asymptotically. We derive the posterior convergence rate, referred to as the contraction rate. We show that the convergence holds even when the dimensionality of the latent feature matrix increases with the sample size, therefore making the posterior inference valid in high dimensional setting. We demonstrate the theoretical results using computer simulation, in which the parallel-tempering Markov chain Monte Carlo method is applied to overcome computational hurdles. The practical utility of the derived properties is demonstrated by inferring the latent features in a reverse phase protein arrays (RPPA) dataset under the IBP prior model. Software and dataset reported in the manuscript are provided at http://www.compgenome.org/IBP.


title: Extremal rays of the embedded subgroup saturation cone

abstract: We examine the extremal rays of the cone of dominant weights $(\mu, \widehat\mu)$ for groups $G\subseteq \widehat G$ for which there exists $N \gg0$ such that $$ \left(V(N\mu)\otimes V(N\widehat \mu)\right)^G\ne (0). $$ We exhibit formulas for a class of rays ("type I") on any regular face of the cone. These rays are identified thanks to a generalization of Fulton's conjecture, which we prove along the way. We verify that the remaining rays ("type II") on the face are the images of extremal rays for a smaller cone under a certain map, whose formula is given. A procedure is given for finding the rays of the cone not on any regular face. This is a generalization of the work of Belkale and Kiers on extremal rays for the saturated tensor cone; the specialization is given by $\widehat G = G\times G$ with the diagonal embedding of $G$. We include several examples to illustrate the formulas.


title: Propagated Perturbation of Adversarial Attack for well-known CNNs: Empirical Study and its Explanation

abstract: Deep Neural Network based classifiers are known to be vulnerable to perturbations of inputs constructed by an adversarial attack to force misclassification. Most studies have focused on how to make vulnerable noise by gradient based attack methods or to defense model from adversarial attack. The use of the denoiser model is one of a well-known solution to reduce the adversarial noise although classification performance had not significantly improved. In this study, we aim to analyze the propagation of adversarial attack as an explainable AI(XAI) point of view. Specifically, we examine the trend of adversarial perturbations through the CNN architectures. To analyze the propagated perturbation, we measured normalized Euclidean Distance and cosine distance in each CNN layer between the feature map of the perturbed image passed through denoiser and the non-perturbed original image. We used five well-known CNN based classifiers and three gradient-based adversarial attacks. From the experimental results, we observed that in most cases, Euclidean Distance explosively increases in the final fully connected layer while cosine distance fluctuated and disappeared at the last layer. This means that the use of denoiser can decrease the amount of noise. However, it failed to defense accuracy degradation.


title: Comparing distributions: $\ell_1$ geometry improves kernel two-sample testing

abstract: Are two sets of observations drawn from the same distribution? This problem is a two-sample test. Kernel methods lead to many appealing properties. Indeed state-of-the-art approaches use the $L^2$ distance between kernel-based distribution representatives to derive their test statistics. Here, we show that $L^p$ distances (with $p\geq 1$) between these distribution representatives give metrics on the space of distributions that are well-behaved to detect differences between distributions as they metrize the weak convergence. Moreover, for analytic kernels, we show that the $L^1$ geometry gives improved testing power for scalable computational procedures. Specifically, we derive a finite dimensional approximation of the metric given as the $\ell_1$ norm of a vector which captures differences of expectations of analytic functions evaluated at spatial locations or frequencies (i.e, features). The features can be chosen to maximize the differences of the distributions and give interpretable indications of how they differs. Using an $\ell_1$ norm gives better detection because differences between representatives are dense as we use analytic kernels (non-zero almost everywhere). The tests are consistent, while much faster than state-of-the-art quadratic-time kernel-based tests. Experiments on artificial and real-world problems demonstrate improved power/time tradeoff than the state of the art, based on $\ell_2$ norms, and in some cases, better outright power than even the most expensive quadratic-time tests.