title: Improving the Detection of Burnt Areas in Remote Sensing using Hyper-features Evolved by M3GP

abstract: One problem found when working with satellite images is the radiometric variations across the image and different images. Intending to improve remote sensing models for the classification of burnt areas, we set two objectives. The first is to understand the relationship between feature spaces and the predictive ability of the models, allowing us to explain the differences between learning and generalization when training and testing in different datasets. We find that training on datasets built from more than one image provides models that generalize better. These results are explained by visualizing the dispersion of values on the feature space. The second objective is to evolve hyper-features that improve the performance of different classifiers on a variety of test sets. We find the hyper-features to be beneficial, and obtain the best models with XGBoost, even if the hyper-features are optimized for a different method.


title: Implementing a neural network interatomic model with performance portability for emerging exascale architectures

abstract: The two main thrusts of computational science are more accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved on-node scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution results in improved performance in many cases, with thread-scalability enabling simulations up to 21 million atoms on a single CPU node and 2 million atoms on a single GPU. We also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~25% and ~10% improvements in performance on a GPU simply by choosing the right level of parallelism and data layout, respectively.


title: A Variational Quantum Algorithm for Preparing Quantum Gibbs States

abstract: Preparation of Gibbs distributions is an important task for quantum computation. It is a necessary first step in some types of quantum simulations and further is essential for quantum algorithms such as quantum Boltzmann training. Despite this, most methods for preparing thermal states are impractical to implement on near-term quantum computers because of the memory overheads required. Here we present a variational approach to preparing Gibbs states that is based on minimizing the free energy of a quantum system. The key insight that makes this practical is the use of Fourier series approximations to the logarithm that allows the entropy component of the free-energy to be estimated through a sequence of simpler measurements that can be combined together using classical post processing. We further show that this approach is efficient for generating high-temperature Gibbs states, within constant error, if the initial guess for the variational parameters for the programmable quantum circuit are sufficiently close to a global optima. Finally, we examine the procedure numerically and show the viability of our approach for five-qubit Hamiltonians using Trotterized adiabatic state preparation as an ansatz.


title: Sub-Doppler Cooling and Compressed Trapping of YO Molecules at $\mu$K Temperatures

abstract: Complex molecular structure demands customized solutions to laser cooling by extending its general set of principles and practices. Yttrium monoxide (YO) has unique intramolecular interactions. The Fermi-contact interaction dominates over the spin-rotation coupling, resulting in two manifolds of closely spaced states, with one of them possessing a negligible Land\'e g-factor. This unique energy level structure favors dual-frequency DC magneto-optical trapping (MOT) and gray molasses cooling (GMC). We report exceptionally robust cooling of YO at 4 $\mu$K over a wide range of laser intensity, detunings (one and two-photon), and magnetic field. The magnetic insensitivity enables the spatial compression of the molecular cloud by alternating GMC and MOT under the continuous operation of the quadrupole magnetic field. A combination of these techniques produces a laser-cooled molecular sample with the highest phase space density in free space.


title: Last Iterate is Slower than Averaged Iterate in Smooth Convex-Concave Saddle Point Problems

abstract: In this paper we study the smooth convex-concave saddle point problem. Specifically, we analyze the last iterate convergence properties of the Extragradient (EG) algorithm. It is well known that the ergodic (averaged) iterates of EG converge at a rate of $O(1/T)$ (Nemirovski, 2004). In this paper, we show that the last iterate of EG converges at a rate of $O(1/\sqrt{T})$. To the best of our knowledge, this is the first paper to provide a convergence rate guarantee for the last iterate of EG for the smooth convex-concave saddle point problem. Moreover, we show that this rate is tight by proving a lower bound of $\Omega(1/\sqrt{T})$ for the last iterate. This lower bound therefore shows a quadratic separation of the convergence rates of ergodic and last iterates in smooth convex-concave saddle point problems.


title: Event triggering control for dynamical systems with designable minimum inter-event time

abstract: This paper presents a class of event-triggering rules for dynamical control systems with guaranteed positive minimum inter-event time (MIET). We first propose an event-based function design with guaranteed control performance under a clock-like variable for general nonlinear systems, and later specify them to general linear systems. Compared to the existing static and dynamic triggering mechanisms, the proposed triggering rules hold the robust global event-separation property, and can be easily implemented on practical digital platform. Namely, it is shown that the minimum inter-event time can be flexibly adapted to the various hardware limitations. Finally, several numerical simulations are given to illustrate the theoretical results.


title: Optimizing Loss Functions Through Multivariate Taylor Polynomial Parameterization

abstract: Metalearning of deep neural network (DNN) architectures and hyperparameters has become an increasingly important area of research. Loss functions are a type of metaknowledge that is crucial to effective training of DNNs, however, their potential role in metalearning has not yet been fully explored. Whereas early work focused on genetic programming (GP) on tree representations, this paper proposes continuous CMA-ES optimization of multivariate Taylor polynomial parameterizations. This approach, TaylorGLO, makes it possible to represent and search useful loss functions more effectively. In MNIST, CIFAR-10, and SVHN benchmark tasks, TaylorGLO finds new loss functions that outperform functions previously discovered through GP, as well as the standard cross-entropy loss, in fewer generations. These functions serve to regularize the learning task by discouraging overfitting to the labels, which is particularly useful in tasks where limited training data is available. The results thus demonstrate that loss function optimization is a productive new avenue for metalearning.


title: Stochastic Extensible Bin Packing

abstract: We consider the stochastic extensible bin packing problem (SEBP) in which $n$ items of stochastic size are packed into $m$ bins of unit capacity. In contrast to the classical bin packing problem, the number of bins is fixed and they can be extended at extra cost. This problem plays an important role in stochastic environments such as in surgery scheduling: Patients must be assigned to operating rooms beforehand, such that the regular capacity is fully utilized while the amount of overtime is as small as possible. This paper focuses on essential ratios between different classes of policies: First, we consider the price of non-splittability, in which we compare the optimal non-anticipatory policy against the optimal fractional assignment policy. We show that this ratio has a tight upper bound of $2$. Moreover, we develop an analysis of a fixed assignment variant of the LEPT rule yielding a tight approximation ratio of $(1+e^{-1}) \approx 1.368$. Furthermore, we prove that the price of fixed assignments, related to the benefit of adaptivity, which describes the loss when restricting to fixed assignment policies, is within the same factor. This shows that in some sense, LEPT is the best fixed assignment policy we can hope for. We also provide a lower bound on the performance of this policy comparing against an optimal fixed assignment policy. Finally, we obtain improved bounds for the case where the processing times are drawn from a particular family of distributions, with either a bounded Pietra index or when the familly is stochastically dominated at the second order.


title: Virial Relations in Density Embedding

abstract: The accuracy of charge-transfer excitation energies, solvatochromic shifts and other environmental effects calculated via various density embedding techniques depend critically on the approximations employed for the non-additive non-interacting kinetic energy functional, $T_{\scriptscriptstyle\rm s}^{\scriptscriptstyle\rm nad}[n]$. Approximating this functional remains an important challenge in electronic structure theory. To assist in the development and testing of approximations for $T_{\scriptscriptstyle\rm s}^{\scriptscriptstyle\rm nad}[n]$, we derive two virial relations for fragments in molecules. These establish separate connections between the non-additive kinetic energies of the non-interacting and interacting systems of electrons, and quantities such as the electron-nuclear attraction forces, the partition (or embedding) energy and potential, and the Kohn-Sham potentials of the system and its parts. We numerically verify both relations on diatomic molecules.


title: Isometric Embeddings of Finite Metric Trees into $(\mathbb{R}^n,d_{1})$ and $(\mathbb{R}^n,d_{\infty})$

abstract: We investigate isometric embeddings of finite metric trees into $(\mathbb{R}^n,d_{1})$ and $( \mathbb{R}^n, d_{\infty})$. We prove that a finite metric tree can be isometrically embedded into $(\mathbb{R}^n,d_{1})$ if and only if the number of its leaves is at most $2n$. We show that a finite star tree with at most $2^n$ leaves can be isometrically embedded into $(\mathbb{R}^{n}, d_{\infty})$ and a finite metric tree with more than $2^n$ leaves cannot be isometrically embedded into $(\mathbb{R}^{n}, d_{\infty})$. We conjecture that an arbitrary finite metric tree with at most $2^n$ leaves can be isometrically embedded into $(\mathbb{R}^{n}, d_{\infty})$.


title: Readout of Superconducting Nanowire Single Photon Detectors through Forward Biased Optical Modulators

abstract: Scalable, high speed data transfer between cryogenic (0.1-4 K) and room temperature environments is instrumental in a broad range of fields including quantum computing, superconducting electronics, single photon imaging and space-based communications. A promising approach to overcome the limitations of conventional wire-based readout is the use of optical fiber communication. Optical fiber presents a 100-1,000x lower heat load than conventional electrical wiring, relaxing the requirements for thermal anchoring, and is also immune to electromagnetic interference, which allows routing of sensitive signals with improved robustness to noise and crosstalk. Most importantly, optical fibers allow for very high bandwidth densities (in the Tbps/mm$^2$ range) by carrying multiple signals through the same physical fiber (Wavelength Division Multiplexing, WDM). Here, we demonstrate for the first time optical readout of a superconducting nanowire single-photon detector (SNSPD) directly coupled to a CMOS photonic modulator, without the need for an interfacing device. By operating the modulator in the forward bias regime at a temperature of 3.6 K, we achieve very high modulation efficiency (1000-10,000 pm/V) and a low input impedance of 500 $\Omega$ with a low power dissipation of 40 $\mu$W. This allows us to obtain optical modulation with the low, millivolt-level signal generated by the SNSPD. Optical communication is becoming the preferred I/O solution in modern, room temperature high performance systems, and this work demonstrates its suitability for scalable cryogenic readout, which could help achieve the full potential of superconducting technologies.


title: The ELM Survey. VIII. 98 Double White Dwarf Binaries

abstract: We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15<g_0<20 mag. The measurements are consistent with 100% of ELM WDs being 0.0089 < P < 1.5 d double WD binaries, 35% of which belong to the Galactic halo. We infer these are mostly He+CO WD binaries given the measurement constraints. The merger rate of the observed He+CO WD binaries exceeds the formation rate of stable mass transfer AM CVn binaries by a factor of 25, and so the majority of He+CO WD binaries must experience unstable mass transfer and merge. The shortest-period systems like J0651+2844 are signature {\it LISA} verification binaries that can be studied with gravitational waves and light.


title: Analysis of Gender Inequality In Face Recognition Accuracy

abstract: We present a comprehensive analysis of how and why face recognition accuracy differs between men and women. We show that accuracy is lower for women due to the combination of (1) the impostor distribution for women having a skew toward higher similarity scores, and (2) the genuine distribution for women having a skew toward lower similarity scores. We show that this phenomenon of the impostor and genuine distributions for women shifting closer towards each other is general across datasets of African-American, Caucasian, and Asian faces. We show that the distribution of facial expressions may differ between male/female, but that the accuracy difference persists for image subsets rated confidently as neutral expression. The accuracy difference also persists for image subsets rated as close to zero pitch angle. Even when removing images with forehead partially occluded by hair/hat, the same impostor/genuine accuracy difference persists. We show that the female genuine distribution improves when only female images without facial cosmetics are used, but that the female impostor distribution also degrades at the same time. Lastly, we show that the accuracy difference persists even if a state-of-the-art deep learning method is trained from scratch using training data explicitly balanced between male and female images and subjects.


title: Simultaneous Skull Conductivity and Focal Source Imaging from EEG Recordings with the help of Bayesian Uncertainty Modelling

abstract: The electroencephalography (EEG) source imaging problem is very sensitive to the electrical modelling of the skull of the patient under examination. Unfortunately, the currently available EEG devices and their embedded software do not take this into account; instead, it is common to use a literature-based skull conductivity parameter. In this paper, we propose a statistical method based on the Bayesian approximation error approach to compensate for source imaging errors due to the unknown skull conductivity and, simultaneously, to compute a low-order estimate for the actual skull conductivity value. By using simulated EEG data that corresponds to focal source activity, we demonstrate the potential of the method to reconstruct the underlying focal sources and low-order errors induced by the unknown skull conductivity. Subsequently, the estimated errors are used to approximate the skull conductivity. The results indicate clear improvements in the source localization accuracy and feasible skull conductivity estimates.


title: Local vibrational modes of Si vacancy spin qubits in SiC

abstract: Silicon carbide is a very promising platform for quantum applications because of extraordinary spin and optical properties of point defects in this technologically-friendly material. These properties are strongly influenced by crystal vibrations, but the exact relationship between them and the behavior of spin qubits is not fully investigated. We uncover the local vibrational modes of the Si vacancy spin qubits in as-grown 4H-SiC. We apply the resonant microwave field to isolate the contribution from one particular type of defects, the so-called V2 center, and observe the zero-phonon line together with seven equally-separated phonon replicas. Furthermore, we present first-principles calculations of the photoluminescence lineshape, which are in excellent agreement with our experimental data. To boost up the calculation accuracy and decrease the computation time, we extract the force constants using machine learning algorithms. This allows us to identify dominant modes in the lattice vibrations coupled to an excited electron during optical emission in the Si vacancy. The resonance phonon energy of 36 meV and the Debye-Waller factor of about 6% are obtained. We establish experimentally that the activation energy of the optically-induced spin polarization is given by the local vibrational energy. Our findings give insight into the coupling of electronic states to vibrational modes in SiC spin qubits, which is essential to predict their spin, optical, mechanical and thermal properties. The approach described can be applied to a large variety of spin defects with spectrally overlapped contributions in SiC as well as in other 3D and 2D materials.


title: DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various Density Functional Theory packages

abstract: Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known correlated materials, namely LaNiO3, SrVO3, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.


title: Battery draining attacks against edge computing nodes in IoT networks

abstract: Many IoT devices, especially those deployed at the network edge have limited power resources. A number of attacks aim to exhaust these resources and drain the batteries of such edge nodes. In this work, we study the effects of a variety of battery draining attacks against edge nodes. Through simulation, we clarify the extent to which such attacks are able to increase the usage and hence waste the power resources of edge nodes. Specifically, we implement hello flooding, packet flooding, selective forwarding, rank attack, and versioning attack in ContikiOS and simulate them in the Cooja simulator, and measure and report a number of time and power resource usage metrics including CPU time, low power mode time, TX/RX time, and battery consumption. Besides, we test the stretch attack with three different batteries as an extreme scenario. Our extensive measurements enable us to compare the effectiveness of these attacks. Our results show that Versioning attack is the most severe attack in terms of draining the power resources of the network, followed by Packet Flooding and Hello Flood attacks. Furthermore, we confirm that Selective Forwarding and Rank attacks are not able to considerably increase the power resource usage in our scenarios. By quantifying the effects of these attacks, we demonstrate that under specific scenarios, Versioning attack can be three to four times as effective as Packet Flooding and Hello Flood attacks in wasting network resources, while Packet Flooding is generally comparable to Hello Flood in CPU and TX time usage increase but twice as powerful in draining device batteries.


title: Analysis and optimal control of a malaria mathematical model under resistance and population movement

abstract: In this work, two mathematical models for malaria under resistance are presented. More precisely, the first model shows the interaction between humans and mosquitoes inside a patch under infection of malaria when the human population is resistant to antimalarial drug and mosquitoes population is resistant to insecticides. For the second model, human-mosquitoes population movements in two patches is analyzed under the same malaria transmission dynamic established in one patch. For a single patch, existence and stability conditions for the equilibrium solutions in terms of the local basic reproductive number are developed. These results reveal the existence of a forward bifurcation and the global stability of disease-free equilibrium. In the case of two patches, a theoretical and numerical framework on sensitivity analysis of parameters is presented. After that, the use of antimalarial drugs and insecticides are incorporated as control strategies and an optimal control problem is formulated. Numerical experiments are carried out in both models to show the feasibility of our theoretical results.


title: Rigorous Guarantees for Tyler's M-estimator via quantum expansion

abstract: Estimating the shape of an elliptical distribution is a fundamental problem in statistics. One estimator for the shape matrix, Tyler's M-estimator, has been shown to have many appealing asymptotic properties. It performs well in numerical experiments and can be quickly computed in practice by a simple iterative procedure. Despite the many years the estimator has been studied in the statistics community, there was neither a tight non-asymptotic bound on the rate of the estimator nor a proof that the iterative procedure converges in polynomially many steps. Here we observe a surprising connection between Tyler's M-estimator and operator scaling, which has been intensively studied in recent years in part because of its connections to the Brascamp-Lieb inequality in analysis. We use this connection, together with novel results on quantum expanders, to show that Tyler's M-estimator has the optimal rate up to factors logarithmic in the dimension, and that in the generative model the iterative procedure has a linear convergence rate even without regularization.


title: Data Augmentation for Histopathological Images Based on Gaussian-Laplacian Pyramid Blending

abstract: Data imbalance is a major problem that affects several machine learning (ML) algorithms. Such a problem is troublesome because most of the ML algorithms attempt to optimize a loss function that does not take into account the data imbalance. Accordingly, the ML algorithm simply generates a trivial model that is biased toward predicting the most frequent class in the training data. In the case of histopathologic images (HIs), both low-level and high-level data augmentation (DA) techniques still present performance issues when applied in the presence of inter-patient variability; whence the model tends to learn color representations, which is related to the staining process. In this paper, we propose a novel approach capable of not only augmenting HI dataset but also distributing the inter-patient variability by means of image blending using the Gaussian-Laplacian pyramid. The proposed approach consists of finding the Gaussian pyramids of two images of different patients and finding the Laplacian pyramids thereof. Afterwards, the left-half side and the right-half side of different HIs are joined in each level of the Laplacian pyramid, and from the joint pyramids, the original image is reconstructed. This composition combines the stain variation of two patients, avoiding that color differences mislead the learning process. Experimental results on the BreakHis dataset have shown promising gains vis-a-vis the majority of DA techniques presented in the literature.


title: Constrained Deep Reinforcement Learning for Energy Sustainable Multi-UAV based Random Access IoT Networks with NOMA

abstract: In this paper, we apply the Non-Orthogonal Multiple Access (NOMA) technique to improve the massive channel access of a wireless IoT network where solar-powered Unmanned Aerial Vehicles (UAVs) relay data from IoT devices to remote servers. Specifically, IoT devices contend for accessing the shared wireless channel using an adaptive $p$-persistent slotted Aloha protocol; and the solar-powered UAVs adopt Successive Interference Cancellation (SIC) to decode multiple received data from IoT devices to improve access efficiency. To enable an energy-sustainable capacity-optimal network, we study the joint problem of dynamic multi-UAV altitude control and multi-cell wireless channel access management of IoT devices as a stochastic control problem with multiple energy constraints. To learn an optimal control policy, we first formulate this problem as a Constrained Markov Decision Process (CMDP), and propose an online model-free Constrained Deep Reinforcement Learning (CDRL) algorithm based on Lagrangian primal-dual policy optimization to solve the CMDP. Extensive simulations demonstrate that our proposed algorithm learns a cooperative policy among UAVs in which the altitude of UAVs and channel access probability of IoT devices are dynamically and jointly controlled to attain the maximal long-term network capacity while maintaining energy sustainability of UAVs. The proposed algorithm outperforms Deep RL based solutions with reward shaping to account for energy costs, and achieves a temporal average system capacity which is $82.4\%$ higher than that of a feasible DRL based solution, and only $6.47\%$ lower compared to that of the energy-constraint-free system.


title: An embedded variable step IMEX scheme for the incompressible Navier-Stokes equations

abstract: This report presents a series of implicit-explicit (IMEX) variable timestep algorithms for the incompressible Navier-Stokes equations (NSE). With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first-second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable step IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.


title: Nonlinear realization of chiral symmetry breaking in holographic soft wall models

abstract: We investigate nonlinear extensions of the holographic soft wall model proposed by Karch, Katz, Son and Stephanov \cite{Karch:2006pv} with a positive quadratic dilaton. We consider a Higgs potential for the tachyonic field that brings a more natural realisation of chiral symmetry breaking in the infrared regime. Utilising the AdS/CFT dictionary and holographic renormalisation we find the chiral condensate as a function of the quark mass. The nonlinearity of the Higgs potential leads to a nonlinear relation between the chiral condensate and the quark mass. Solving the effective Schr\"odinger equations for the field perturbations we estimate meson masses and decay constants and evaluate their dependence on the quark mass. In the axial and pseudoscalar sector we find an interesting behaviour for the decay constants as the quark mass increases. We also investigate the effect of a 5d running mass for the tachyonic field. We conclude that nonlinear soft wall models with a Higgs potential for the tachyon and a positive quadratic dilaton do not provide spontaneous symmetry breaking in the chiral limit.


title: Symmetry and compact embeddings for critical exponents in metric-measure spaces

abstract: We obtain a compact Sobolev embedding for $H$-invariant functions in compact metric-measure spaces, where $H$ is a subgroup of the measure preserving bijections. In Riemannian manifolds, $H$ is a subgroup of the volume preserving diffeomorphisms: a compact embedding for the critical exponents follows. The results can be viewed as an extension of Sobolev embeddings of functions invariant under isometries in compact manifolds.


title: Topological magnons in CrI$_3$ monolayers: an itinerant fermion description

abstract: Magnons dominate the magnetic response of the recently discovered insulating ferromagnetic two dimensional crystals such as CrI$_3$. Because of the arrangement of the Cr spins in a honeycomb lattice, magnons in CrI$_3$ bear a strong resemblance with electronic quasiparticles in graphene. Neutron scattering experiments carried out in bulk CrI$_3$ show the existence of a gap at the Dirac points, that has been conjectured to have a topological nature. Here we propose a theory for magnons in ferromagnetic CrI$_3$ monolayers based on an itinerant fermion picture, with a Hamiltonian derived from first principles. We obtain the magnon dispersion for 2D CrI$_3$ with a gap at the Dirac points with the same Berry curvature in both valleys. For CrI$_3$ ribbons, we find chiral in-gap edge states. Analysis of the magnon wave functions in momentum space further confirms their topological nature. Importantly, our approach does not require to define a spin Hamiltonian, and can be applied to both insulating and conducting 2D materials with any type of magnetic order.


title: Effects of Roots of Maximal Multiplicity on the Stability of Some Classes of Delay Differential-Algebraic Systems: The Lossless Propagation Case

abstract: It has been observed in several recent works that, for some classes of linear time-delay systems, spectral values of maximal multiplicity are dominant, a property known as multiplicity-induced-dominancy (MID). This paper starts the investigation of whether MID holds for delay differential-algebraic systems by considering a single-delay system composed of two scalar equations. After motivating this problem and recalling some recent results for retarded delay differential equations, we prove that the MID property holds for the delay differential-algebraic system of interest and present some applications.


title: Boosting Algorithms for Estimating Optimal Individualized Treatment Rules

abstract: We present nonparametric algorithms for estimating optimal individualized treatment rules. The proposed algorithms are based on the XGBoost algorithm, which is known as one of the most powerful algorithms in the machine learning literature. Our main idea is to model the conditional mean of clinical outcome or the decision rule via additive regression trees, and use the boosting technique to estimate each single tree iteratively. Our approaches overcome the challenge of correct model specification, which is required in current parametric methods. The major contribution of our proposed algorithms is providing efficient and accurate estimation of the highly nonlinear and complex optimal individualized treatment rules that often arise in practice. Finally, we illustrate the superior performance of our algorithms by extensive simulation studies and conclude with an application to the real data from a diabetes Phase III trial.


title: Convergence rate analysis and improved iterations for numerical radius computation

abstract: The main two algorithms for computing the numerical radius are the level-set method of Mengi and Overton and the cutting-plane method of Uhlig. Via new analyses, we explain why the cutting-plane approach is sometimes much faster or much slower than the level-set one and then propose a new hybrid algorithm that remains efficient in all cases. For matrices whose fields of values are a circular disk centered at the origin, we show that the cost of Uhlig's method blows up with respect to the desired relative accuracy. More generally, we also analyze the local behavior of Uhlig's cutting procedure at outermost points in the field of values, showing that it often has a fast Q-linear rate of convergence and is Q-superlinear at corners. Finally, we identify and address inefficiencies in both the level-set and cutting-plane approaches and propose refined versions of these techniques.


title: Schubert polynomial analogues for degenerate involutions

abstract: We survey the recent study of involution Schubert polynomials and a modest generalization that we call degenerate involution Schubert polynomials. We cite several conditions when (degenerate) involution Schubert polynomials have simple factorization formulae. Such polynomials can be computed by traversing through chains in certain weak order posets, and we provide explicit descriptions of such chains in weak order for involutions and degenerate involutions. As an application, we give several examples of how certain multiplicity-free sums of Schubert polynomials factor completely into very simple linear factors.


title: Regret Minimization in Partially Observable Linear Quadratic Control

abstract: We study the problem of regret minimization in partially observable linear quadratic control systems when the model dynamics are unknown a priori. We propose ExpCommit, an explore-then-commit algorithm that learns the model Markov parameters and then follows the principle of optimism in the face of uncertainty to design a controller. We propose a novel way to decompose the regret and provide an end-to-end sublinear regret upper bound for partially observable linear quadratic control. Finally, we provide stability guarantees and establish a regret upper bound of $\tilde{\mathcal{O}}(T^{2/3})$ for ExpCommit, where $T$ is the time horizon of the problem.


title: The cone theorem and the vanishing of Chow cohomology

abstract: We show that a cone theorem for ${\mathbbA}^1-homotopy invariant contravariant functors implies the vanishing of the positive degree part of the operational Chow cohomology rings of a large class of affine varieties. We also discuss how this vanishing relates to a number of questions about representing Chow cohomology classes of GIT quotients in terms of equivariant cycles.


title: Approximate Summaries for Why and Why-not Provenance (Extended Version)

abstract: Why and why-not provenance have been studied extensively in recent years. However, why-not provenance, and to a lesser degree why provenance, can be very large resulting in severe scalability and usability challenges. In this paper, we introduce a novel approximate summarization technique for provenance which overcomes these challenges. Our approach uses patterns to encode (why-not) provenance concisely. We develop techniques for efficiently computing provenance summaries balancing informativeness, conciseness, and completeness. To achieve scalability, we integrate sampling techniques into provenance capture and summarization. Our approach is the first to scale to large datasets and to generate comprehensive and meaningful summaries.


title: PCA for Implied Volatility Surfaces

abstract: Principal component analysis (PCA) is a useful tool when trying to construct factor models from historical asset returns. For the implied volatilities of U.S. equities there is a PCA-based model with a principal eigenportfolio whose return time series lies close to that of an overarching market factor. The authors show that this market factor is the index resulting from the daily compounding of a weighted average of implied-volatility returns, with weights based on the options' open interest (OI) and Vega. The authors also analyze the singular vectors derived from the tensor structure of the implied volatilities of S&P500 constituents, and find evidence indicating that some type of OI and Vega-weighted index should be one of at least two significant factors in this market.


title: Ultrafast functional magnetic resonance imaging reveals neuroplasticity-driven timing modulations

abstract: Functional Magnetic Resonance Imaging (fMRI) is predominantly harnessed for spatially mapping activation foci along distributed pathways. However, resolving dynamic information on activation sequence remains elusive. Here, we show an ultra-fast fMRI (ufMRI) approach - a facilitating non-invasive methodology for mapping Blood-Oxygenation-Level-Dependent (BOLD) response timings in distributed pathways with high spatiotemporal sensitivity and resolution. The mouse visual pathway was investigated under both normal and dark reared conditions. Results show that BOLD responses of normal reared mice preserve the neural input order from onset to peak times. However, modulatory effects in cortical responses due to dark rearing are only measurable at early response timings while a general delay in responses is measured at later timings. Our findings highlight the importance of robustly measuring early BOLD timings and pave the way for a better understanding and interpretation of functional MRI in healthy and aberrant conditions.


title: Exact and Robust Reconstructions of Integer Vectors Based on Multidimensional Chinese Remainder Theorem (MD-CRT)

abstract: The robust Chinese remainder theorem (CRT) has been recently proposed for robustly reconstructing a large nonnegative integer from erroneous remainders. It has found many applications in signal processing, including phase unwrapping and frequency estimation under sub-Nyquist sampling. Motivated by the applications in multidimensional (MD) signal processing, in this paper we propose the MD-CRT and robust MD-CRT for integer vectors. Specifically, by rephrasing the abstract CRT for rings in number-theoretic terms, we first derive the MD-CRT for integer vectors with respect to a general set of integer matrix moduli, which provides an algorithm to uniquely reconstruct an integer vector from its remainders, if it is in the fundamental parallelepiped of the lattice generated by a least common right multiple of all the moduli. For some special forms of moduli, we present explicit reconstruction formulae. Moreover, we derive the robust MD-CRT for integer vectors when the remaining integer matrices of all the moduli left divided by their greatest common left divisor (gcld) are pairwise commutative and coprime. Two different reconstruction algorithms are proposed, and accordingly, two different conditions on the remainder error bound for the reconstruction robustness are obtained, which are related to a quarter of the minimum distance of the lattice generated by the gcld of all the moduli or the Smith normal form of the gcld.


title: PatchworkWave: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe/Multimethod Simulations at Arbitrary Order

abstract: We present an extension of the PatchworkMHD code [1], itself an MHD-capable extension of the Patchwork code [2], for which several algorithms presented here were co-developed. Its purpose is to create a multipatch scheme compatible with numerical simulations of arbitrary equations of motion at any discretization order in space and time. In the Patchwork framework, the global simulation is comprised of an arbitrary number of moving, local meshes, or patches, which are free to employ their own resolution, coordinate system/topology, physics equations, reference frame, and in our new approach, numerical method. Each local patch exchanges boundary data with a single global patch on which all other patches reside through a client-router-server parallelization model. In generalizing Patchwork to be compatible with arbitrary order time integration, PatchworkMHD and PatchworkWave have significantly improved the interpatch interpolation accuracy by removing an interpolation of interpolated data feedback present in the original Patchwork code. Furthermore, we extend Patchwork to be multimethod by allowing multiple state vectors to be updated simultaneously, with each state vector providing its own interpatch interpolation and transformation procedures. As such, our scheme is compatible with nearly any set of hyperbolic partial differential equations. We demonstrate our changes through the implementation of a scalar wave toy-model that is evolved on arbitrary, time dependent patch configurations at 4th order accuracy.


title: Harmonizing Child Mortality Data at Disparate Geographic Levels

abstract: There is an increasing focus on reducing inequalities in health outcomes in developing countries. Subnational variation is of particular interest, with geographic data used to understand the spatial risk of detrimental outcomes and to identify who is at greatest risk. While some health surveys provide observations with associated geographic coordinates, many others provide data that have their locations masked and instead only report the strata within which the data resides. How to harmonize these data sources for spatial analysis has seen previously considered though no method has been agreed upon and comparison of the validity of methods are lacking. In this paper, we present a new method for analyzing masked survey data alongside traditional geolocated data, using a method that is consistent with the data generating process. In addition, we critique two proposed approaches to analyzing masked data and illustrate that they are fundamentally flawed methodologically. To validate our method, we compare our approach with previously formulated solutions in several realistic simulation environments in which the underlying structure of the risk field is known. We simulate samples from spatial fields in a way that mimics the sampling frame implemented in the most common health surveys in low and middle income countries, the DHS and MICS. In simulations, the newly proposed approach outperforms previously proposed approaches in terms of minimizing error while increasing the precision of estimates. The approaches are subsequently compared using child mortality data from the Dominican Republic where our findings are reinforced. Accurately increasing precision of child mortality estimates, and health estimates in general, by leveraging various types of data improves our ability to implement precision public health initiatives and better understand the landscape of geographic health inequalities.


title: Mathematical Modeling, Laboratory Experiments, and Sensitivity Analysis of Bioplug Technology at Darcy Scale

abstract: In this paper we study a Darcy-scale mathematical model for biofilm formation in porous media. The pores in the core are divided into three phases: water, oil, and biofilm. The water and oil flow are modeled by an extended version of Darcy's law and the substrate is transported by diffusion and convection in the water phase. Initially there is biofilm on the pore walls. The biofilm consumes substrate for production of biomass and modifies the pore space which changes the rock permeability. The model includes detachment of biomass due to water flux and death of bacteria, and is implemented in MRST. We discuss the capability of the numerical simulator to capture results from laboratory experiments. We perform a novel sensitivity analysis based on sparse-grid interpolation and multi-wavelet expansion to identify the critical model parameters. Numerical experiments using diverse injection strategies are performed to study the impact of different porosity-permeability relations in a core saturated with water and oil.


title: Using Inaudible Audio to Improve Indoor-Localization- and Proximity-Aware Intelligent Applications

abstract: While it is often critical for indoor-location- and proximity-aware applications to know whether a user is in a space or not (e.g., a specific room or office), a key challenge is that the difference between standing on one side or another of a doorway or wall is well within the error range of most RF-based approaches. In this work, we address this challenge by augmenting RF-based localization and proximity detection with active ultrasonic sensing, taking advantage of the limited propagation of sound waves. This simple and cost-effective approach can allow, for example, a Bluetooth smart-lock to discern whether a user is inside or outside their home in order to lock or unlock doors automatically. We describe a configurable architecture for our solution and present experiments that validate this approach but also demonstrate that different user behavior and application needs can impact system configuration decisions. Finally, we describe applications that could benefit from our solution and address privacy concerns.


title: Hybrid Graph Neural Networks for Crowd Counting

abstract: Crowd counting is an important yet challenging task due to the large scale and density variation. Recent investigations have shown that distilling rich relations among multi-scale features and exploiting useful information from the auxiliary task, i.e., localization, are vital for this task. Nevertheless, how to comprehensively leverage these relations within a unified network architecture is still a challenging problem. In this paper, we present a novel network structure called Hybrid Graph Neural Network (HyGnn) which targets to relieve the problem by interweaving the multi-scale features for crowd density as well as its auxiliary task (localization) together and performing joint reasoning over a graph. Specifically, HyGnn integrates a hybrid graph to jointly represent the task-specific feature maps of different scales as nodes, and two types of relations as edges:(i) multi-scale relations for capturing the feature dependencies across scales and (ii) mutual beneficial relations building bridges for the cooperation between counting and localization. Thus, through message passing, HyGnn can distill rich relations between the nodes to obtain more powerful representations, leading to robust and accurate results. Our HyGnn performs significantly well on four challenging datasets: ShanghaiTech Part A, ShanghaiTech Part B, UCF_CC_50 and UCF_QNRF, outperforming the state-of-the-art approaches by a large margin.


title: Time-smoothing for parabolic variational problems in metric measure spaces

abstract: In 2013, Masson and Siljander determined a method to prove that the $p$-minimal upper gradient $g_{f_\varepsilon}$ for the time mollification $f_\varepsilon$, $\varepsilon>0$, of a parabolic Newton-Sobolev function $f\in L^p_\mathrm{loc}(0,\tau;N^{1,p}_\mathrm{loc}(\Omega))$, with $\tau>0$ and $\Omega$ open domain in a doubling metric measure space $(\mathbb{X},d,\mu)$ supporting a weak $(1,p)$-Poincar\'e inequality, $p\in(1,\infty)$, is such that $g_{f-f_\varepsilon}\to0$ as $\varepsilon\to0$ in $L^p_\mathrm{loc}(\Omega_\tau)$, $\Omega_\tau$ being the parabolic cylinder $\Omega_{\tau}:=\Omega\times(0,\tau)$. Their approach involved the use of Cheeger's differential structure, and therefore exhibited some limitations; here, we shall see that the definition and the formal properties of the parabolic Sobolev spaces themselves allow to find a more direct method to show such convergence, which relies on $p$-weak upper gradients only and which is valid regardless of structural assumptions on the ambient space, also in the limiting case when $p=1$.


title: Relations between large scale brain connectivity and effects of regional stimulation depend on collective dynamical state

abstract: At the macroscale, the brain operates as a network of interconnected neuronal populations, which display rhythmic dynamics that support interareal communication. Understanding how stimulation of a particular brain area impacts such concerted activity is important for gaining basic insights into brain function and for developing neuromodulation as a therapeutic tool. However, it remains difficult to predict the downstream effects of focal stimulation. Specifically, little is known about how the collective oscillatory regime of network activity may affect the outcomes of regional perturbations on cooperative dynamics. Here, we combine connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control the collective dynamics of the network, we identify distinct states of simulated brain activity, and investigate how the distributed effects of stimulation manifest in different states. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, nearby regions develop phase locked activity in the excited frequency band. Importantly, we find that focal stimulation also causes more distributed modifications to network coherence at regions' baseline oscillation frequencies, and that these effects are better predicted by functional rather than structural connectivity. In contrast, when the network operates in a regime of stronger endogenous oscillations, stimulation causes only slight shifts in power and frequency, and network averaged changes in coherence are more homogenous across the choice of the stimulated area. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and highlights that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network wide responses to local perturbations.


title: Variation of fundamental constants and white dwarfs

abstract: Theories that attempt to unify the four fundamental interactions and alternative theories of gravity predict time and/or spatial variation of the fundamental constants of nature. Different versions of these theories predict different behaviours for these variations. In consequence, experimental and observational bounds are an important tool to check the validity of such proposals. In this paper, we review constraints on the possible variation of the fundamental constants from astronomical observations and geophysical experiments designed to test the constancy of the fundamental constants of nature over different timescales. We also focus on the limits that can be obtained from white dwarfs, which can constrain the variation of the constants with the gravitational potential.


title: Simultaneous nonvanishing of the Products of L-functions associated to elliptic cusp forms

abstract: A generalized Riemann hypothesis states that all zeros of the completed Hecke $L$-function $L^*(f,s)$ of a normalized Hecke eigenform $f$ on the full modular group should lie on the vertical line $Re(s)=\frac{k}{2}.$ It was shown by Kohnen that there exists a Hecke eigenform $f$ of weight $k$ such that $L^*(f,s) \neq 0$ for sufficiently large $k$ and any point on the line segments $Im(s)=t_0, \frac{k-1}{2} < Re(s) < \frac{k}{2}-\epsilon, \frac{k }{2}+\epsilon < Re(s) < \frac{k+1}{2},$ for any given real number $t_0$ and a positive real number $\epsilon.$ This paper concerns the non-vanishing of the product $L^*(f,s)L^*(f,w)$ $(s,w\in \mathbb{C})$ on average.


title: Physics-Guided Deep Neural Networks for Power Flow Analysis

abstract: Solving power flow (PF) equations is the basis of power flow analysis, which is important in determining the best operation of existing systems, performing security analysis, etc. However, PF equations can be out-of-date or even unavailable due to system dynamics and uncertainties, making traditional numerical approaches infeasible. To address these concerns, researchers have proposed data-driven approaches to solve the PF problem by learning the mapping rules from historical system operation data. Nevertheless, prior data-driven approaches suffer from poor performance and generalizability, due to overly simplified assumptions of the PF problem or ignorance of physical laws governing power systems. In this paper, we propose a physics-guided neural network to solve the PF problem, with an auxiliary task to rebuild the PF model. By encoding different granularity of Kirchhoff's laws and system topology into the rebuilt PF model, our neural-network based PF solver is regularized by the auxiliary task and constrained by the physical laws. The simulation results show that our physics-guided neural network methods achieve better performance and generalizability compared to existing unconstrained data-driven approaches. Furthermore, we demonstrate that the weight matrices of our physics-guided neural networks embody power system physics by showing their similarities with the bus admittance matrices.


title: Lifetime of almost strong edge-mode operators in one dimensional, interacting, symmetry protected topological phases

abstract: Almost strong edge-mode operators arising at the boundaries of certain interacting 1D symmetry protected topological phases with \(Z_2\) symmetry have infinite temperature lifetimes that are non-perturbatively long in the integrability breaking terms, making them promising as bits for quantum information processing. We extract the lifetime of these edge-mode operators for small system sizes as well as in the thermodynamic limit. For the latter, a Lanczos scheme is employed to map the operator dynamics to a one dimensional tight-binding model of a single particle in Krylov space. We find this model to be that of a spatially inhomogeneous Su-Schrieffer-Heeger model with a hopping amplitude that increases away from the boundary, and a dimerization that decreases away from the boundary. We associate this dimerized or staggered structure with the existence of the almost strong mode. Thus the short time dynamics of the almost strong mode is that of the edge-mode of the Su-Schrieffer-Heeger model, while the long time dynamics involves decay due to tunneling out of that mode, followed by chaotic operator spreading. We also show that competing scattering processes can lead to interference effects that can significantly enhance the lifetime.


title: Dynamic and Distributed Online Convex Optimization for Demand Response of Commercial Buildings

abstract: We extend the regret analysis of the online distributed weighted dual averaging (DWDA) algorithm [1] to the dynamic setting and provide the tightest dynamic regret bound known to date with respect to the time horizon for a distributed online convex optimization (OCO) algorithm. Our bound is linear in the cumulative difference between consecutive optima and does not depend explicitly on the time horizon. We use dynamic-online DWDA (D-ODWDA) and formulate a performance-guaranteed distributed online demand response approach for heating, ventilation, and air-conditioning (HVAC) systems of commercial buildings. We show the performance of our approach for fast timescale demand response in numerical simulations and obtain demand response decisions that closely reproduce the centralized optimal ones.


title: Scalable bundling via dense product embeddings

abstract: Bundling, the practice of jointly selling two or more products at a discount, is a widely used strategy in industry and a well examined concept in academia. Historically, the focus has been on theoretical studies in the context of monopolistic firms and assumed product relationships, e.g., complementarity in usage. We develop a new machine-learning-driven methodology for designing bundles in a large-scale, cross-category retail setting. We leverage historical purchases and consideration sets created from clickstream data to generate dense continuous representations of products called embeddings. We then put minimal structure on these embeddings and develop heuristics for complementarity and substitutability among products. Subsequently, we use the heuristics to create multiple bundles for each product and test their performance using a field experiment with a large retailer. We combine the results from the experiment with product embeddings using a hierarchical model that maps bundle features to their purchase likelihood, as measured by the add-to-cart rate. We find that our embeddings-based heuristics are strong predictors of bundle success, robust across product categories, and generalize well to the retailer's entire assortment.


title: Effects of the transverse coherence length in relativistic collisions

abstract: Effects of the quantum interference in collisions of particles have a twofold nature: they arise because of the auto-correlation of a complex scattering amplitude and due to spatial coherence of the incoming wave packets. Both these effects are neglected in a conventional scattering theory dealing with the delocalized plane waves, although they sometimes must be taken into account in particle and atomic physics. Here, we study the role of a transverse coherence length of the packets, putting special emphasis on the case in which one of the particles is twisted, that is, it carries an orbital angular momentum $\ell\hbar$. In $ee, ep$, and $pp$ collisions the interference results in corrections to the plane-wave cross sections, usually negligible at the energies $\sqrt{s} \gg 1$ GeV but noticeable for smaller ones, especially if there is a twisted hadron with $|\ell| > 10^3$ in initial state. Beyond the perturbative QCD, these corrections become only moderately attenuated allowing one to probe a phase of the hadronic amplitude as a function of $s$ and $t$. In this regime, the coherence effects can compete with the loop corrections in QED and facilitate testing the phenomenological models of the strong interaction at intermediate and low energies.


title: Edge-based sequential graph generation with recurrent neural networks

abstract: Graph generation with Machine Learning is an open problem with applications in various research fields. In this work, we propose to cast the generative process of a graph into a sequential one, relying on a node ordering procedure. We use this sequential process to design a novel generative model composed of two recurrent neural networks that learn to predict the edges of graphs: the first network generates one endpoint of each edge, while the second network generates the other endpoint conditioned on the state of the first. We test our approach extensively on five different datasets, comparing with two well-known baselines coming from graph literature, and two recurrent approaches, one of which holds state of the art performances. Evaluation is conducted considering quantitative and qualitative characteristics of the generated samples. Results show that our approach is able to yield novel, and unique graphs originating from very different distributions, while retaining structural properties very similar to those in the training sample. Under the proposed evaluation framework, our approach is able to reach performances comparable to the current state of the art on the graph generation task.


title: Estimating Welfare Effects in a Nonparametric Choice Model: The Case of School Vouchers

abstract: We develop new robust discrete choice tools to learn about the average willingness to pay for a price subsidy and its effects on demand given exogenous, discrete variation in prices. Our starting point is a nonparametric, nonseparable model of choice. We exploit the insight that our welfare parameters in this model can be expressed as functions of demand for the different alternatives. However, while the variation in the data reveals the value of demand at the observed prices, the parameters generally depend on its values beyond these prices. We show how to sharply characterize what we can learn when demand is specified to be entirely nonparametric or to be parameterized in a flexible manner, both of which imply that the parameters are not necessarily point identified. We use our tools to analyze the welfare effects of price subsidies provided by school vouchers in the DC Opportunity Scholarship Program. We robustly find that the provision of the status quo voucher and a wide range of counterfactual vouchers of different amounts have positive benefits net of costs. This positive effect can be explained by the popularity of low-tuition schools in the program; removing them from the program can result in a negative net benefit. Relative to our bounds, we also find that comparable logit estimates potentially understate the benefits for certain voucher amounts, and provide a misleading sense of robustness for alternative amounts.


title: Post-Training Piecewise Linear Quantization for Deep Neural Networks

abstract: Quantization plays an important role in the energy-efficient deployment of deep neural networks on resource-limited devices. Post-training quantization is highly desirable since it does not require retraining or access to the full training dataset. The well-established uniform scheme for post-training quantization achieves satisfactory results by converting neural networks from full-precision to 8-bit fixed-point integers. However, it suffers from significant performance degradation when quantizing to lower bit-widths. In this paper, we propose a piecewise linear quantization (PWLQ) scheme to enable accurate approximation for tensor values that have bell-shaped distributions with long tails. Our approach breaks the entire quantization range into non-overlapping regions for each tensor, with each region being assigned an equal number of quantization levels. Optimal breakpoints that divide the entire range are found by minimizing the quantization error. Compared to state-of-the-art post-training quantization methods, experimental results show that our proposed method achieves superior performance on image classification, semantic segmentation, and object detection with minor overhead.


title: General upper bound on the game domination number

abstract: It is conjectured that the game domination number is at most $3n/5$ for every $n$-vertex graph which does not contain isolated vertices. It was proved in the recent years that the conjecture holds for several graph classes, including the class of forests and that of graphs with minimum degree at least two. Here we prove that the slightly bigger upper bound $5n/8$ is valid for every isolate-free graph.


title: Orbital dynamics of highly probable but rare Orionid outbursts possibly observed by the ancient Maya

abstract: Using orbital integrations of particles ejected from Comet Halley's passages between 1404 BC and 240 BC, the authors investigate possible outbursts of the Orionids (twin shower of the Eta Aquariids) that may have been observed in the western hemisphere. In an earlier orbital integration study the authors determined there was a high probability linking probable outbursts of the Eta Aquariid meteor shower with certain events recorded in inscriptions during the Maya Classic Period, AD 250-900. This prior examination was the first scientific inquiry of its kind into ancient meteor outbursts possibly recorded in the western hemisphere where previously no pre-Columbian observations had existed. In the current paper the aim is to describe orbital dynamics of rare but probable Orionid outbursts that would have occurred on or near applicable dates recorded in the Classic Maya inscriptions. Specifically, significant probable outbursts are found in AD 417 and 585 out of 30 possible target years. The driving mechanisms for outbursts in those two years are Jovian 1:6 and 1:7 mean motion resonances acting to maintain compact structures within the Orionid stream for over 1 kyr. Furthermore, an Orionid outburst in AD 585 recorded by China is confirmed.


title: Generative Modeling with Denoising Auto-Encoders and Langevin Sampling

abstract: We study convergence of a generative modeling method that first estimates the score function of the distribution using Denoising Auto-Encoders (DAE) or Denoising Score Matching (DSM) and then employs Langevin diffusion for sampling. We show that both DAE and DSM provide estimates of the score of the Gaussian smoothed population density, allowing us to apply the machinery of Empirical Processes. We overcome the challenge of relying only on $L^2$ bounds on the score estimation error and provide finite-sample bounds in the Wasserstein distance between the law of the population distribution and the law of this sampling scheme. We then apply our results to the homotopy method of arXiv:1907.05600 and provide theoretical justification for its empirical success.


title: Measurement of the Spectral Shape of the beta-decay of 137Xe to the Ground State of 137Cs in EXO-200 and Comparison with Theory

abstract: We report on a comparison between the theoretically predicted and experimentally measured spectra of the first-forbidden non-unique $\beta$-decay transition $^{137}\textrm{Xe}(7/2^-)\to\,^{137}\textrm{Cs}(7/2^+)$. The experimental data were acquired by the EXO-200 experiment during a deployment of an AmBe neutron source. The ultra-low background environment of EXO-200, together with dedicated source deployment and analysis procedures, allowed for collection of a pure sample of the decays, with an estimated signal-to-background ratio of more than 99-to-1 in the energy range from 1075 to 4175 keV. In addition to providing a rare and accurate measurement of the first-forbidden non-unique $\beta$-decay shape, this work constitutes a novel test of the calculated electron spectral shapes in the context of the reactor antineutrino anomaly and spectral bump.


title: Pierce stalks in preprimal varieties

abstract: An algebra $\mathbf{P}$ is called \textit{preprimal} if $\mathbf{P}$ is finite and $\func{Clo}(\mathbf{P})$ is a maximal clone. A \textit{preprimal variety} is a variety generated by a preprimal algebra. After Rosenberg's classification of maximal clones \cite{ro}; we have that a finite algebra is preprimal if and only if its term operations are exactly the functions preserving a relation of one of the following seven types: 1. Permutations with cycles all the same prime length, 2. Proper subsets, 3 Prime-affine relations, 4. Bounded partial orders, 5. $h$-adic relations, 6. Central relations $h\geq 2$, 7. Proper, non-trivial equivalence relations. In \cite{kn} Knoebel studies the Pierce sheaf of the different preprimal varieties and he asks for a description of the Pierce stalks. He solves this problem for the cases 1.,2. and 3. and left open the remaining cases. In this paper, using central element theory we succeeded in describing the Pierce stalks of the cases 6. and 7..


title: Lightcone Commutator and Stress-Tensor Exchange in $d>2$ CFTs

abstract: Motivated by developing a field-theoretic algebraic approach to the universal part of the stress-tensor sector of a scalar four-point function in a class of higher-dimensional CFTs, we construct a mode operator, ${\cal L}_m$, near the lightcone in $d=4$ CFTs and show that it leads to a Virasoro-like commutator, including a regularized central-term. As an example, we describe how to reproduce the $d=4$ single-stress tensor exchange contribution in the lightcone limit by a mode summation. A general-$d$ extension is included. We comment on possible generalizations.


title: Modeling supercoiled DNA interacting with an anchored cluster of proteins: towards a quantitative estimation of chromosomal DNA supercoiling

abstract: We investigate the measurement of DNA supercoiling density ($\sigma$) along chromosomes using interaction frequencies between DNA and DNA-anchored clusters of proteins. Specifically, we show how the physics of DNA supercoiling leads, in bacteria, to the quantitative modeling of binding properties of ParB proteins around their centromere-like site, {\it parS}. Using this framework, we provide an upper bound for $\sigma$ in the {\it Escherichia coli} chromosome, consistent with plasmid values, and offer a proof of concept for a high accuracy measurement. To reach these conclusions, we revisit the problem of the formation of ParB clusters. We predict, in particular, that they result from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., they behave like liquid-like protein condensates with unconventional ``leaky'' boundaries.


title: Pseudo-multipliers and smooth molecules on Hermite Besov and Hermite Triebel--Lizorkin spaces

abstract: We obtain new molecular decompositions and molecular synthesis estimates for Hermite Besov and Hermite Triebel--Lizorkin spaces and use such tools to prove boundedness properties of Hermite pseudo-multipliers on those spaces. The notion of molecule we develop leads to boundedness of pseudo-multipliers associated to symbols of H\"ormander-type adapted to the Hermite setting on spaces with non-positive smoothness; in particular, we obtain continuity results for such operators on Lebesgue and Hermite local Hardy spaces. As a byproduct of our results on boundedness properties of pseudo-multipliers, we show that Hermite Besov spaces and Hermite Triebel--Lizorkin spaces are closed under non-linearities.


title: Better Compression with Deep Pre-Editing

abstract: Could we compress images via standard codecs while avoiding visible artifacts? The answer is obvious -- this is doable as long as the bit budget is generous enough. What if the allocated bit-rate for compression is insufficient? Then unfortunately, artifacts are a fact of life. Many attempts were made over the years to fight this phenomenon, with various degrees of success. In this work we aim to break the unholy connection between bit-rate and image quality, and propose a way to circumvent compression artifacts by pre-editing the incoming image and modifying its content to fit the given bits. We design this editing operation as a learned convolutional neural network, and formulate an optimization problem for its training. Our loss takes into account a proximity between the original image and the edited one, a bit-budget penalty over the proposed image, and a no-reference image quality measure for forcing the outcome to be visually pleasing. The proposed approach is demonstrated on the popular JPEG compression, showing savings in bits and/or improvements in visual quality, obtained with intricate editing effects.


title: Discontinuous Galerkin Approach to Large Bending Deformation of a Bilayer Plate with Isometry Constraint

abstract: We present a computational model of thin elastic bilayers that undergo large bending isometric deformations when actuated by non-mechanical stimuli. We propose a discontinuous Galerkin approximation of the variational formulation discussed in Bartels et al. (2018). We showcase its advantages and good computational performance with configurations of interest in both engineering and medicine and either Dirichlet or free boundary conditions.


title: Dissecting the Workload of a Major Adult Video Portal

abstract: Adult content constitutes a major source of Internet traffic. As with many other platforms, these sites are incentivized to engage users and maintain them on the site. This engagement (e.g., through recommendations) shapes the journeys taken through such sites. Using data from a large content delivery network, we explore session journeys within an adult website. We take two perspectives. We first inspect the corpus available on these platforms. Following this, we investigate the session access patterns. We make a number of observations that could be exploited for optimizing delivery, e.g., that users often skip within video streams.


title: Hybridized Summation-By-Parts Finite Difference Methods

abstract: We present a hybridization technique for summation-by-parts finite difference methods with weak enforcement of interface and boundary conditions for second order, linear elliptic partial differential equations. The method is based on techniques from the hybridized discontinuous Galerkin literature where local and global problems are defined for the volume and trace grid points, respectively. By using a Schur complement technique the volume points can be eliminated, which drastically reduces the system size. We derive both the local and global problems, and show that the linear systems that must be solved are symmetric positive definite. The theoretical stability results are confirmed with numerical experiments as is the accuracy of the method.


title: Diagnosis inspired by the thermodynamic geometry for different thermodynamic schemes of the charged BTZ black hole

abstract: Due to the asymptotic structure of the black hole solution, there are two different thermodynamic schemes for the charged Banados-Teitelboim-Zanelli (BTZ) black hole. In one scheme, the charged BTZ black hole is super-entropic, while in the other, it is not (the reverse isoperimetric inequality is saturated). In this paper, we investigate the thermodynamic curvature of the charged BTZ black hole in different coordinate spaces. We find that in both schemes, the thermodynamic curvature is always positive, which may be related to the information of repulsive interaction between black hole molecules for the charged BTZ black hole if we accept an empirical relationship between the thermodynamic curvature and interaction of a system. More importantly, we provide a diagnosis for the discrimination of the two schemes from the point of view of the thermodynamics geometry. For the charged BTZ black hole, when the reverse isoperimetric inequality is saturated, the thermodynamic curvature of an extreme black hole tends to be infinity, while when the reverse isoperimetric inequality is violated, the thermodynamic curvature of the extreme black hole goes to a finite value.


title: AdvectiveNet: An Eulerian-Lagrangian Fluidic reservoir for Point Cloud Processing

abstract: This paper presents a novel physics-inspired deep learning approach for point cloud processing motivated by the natural flow phenomena in fluid mechanics. Our learning architecture jointly defines data in an Eulerian world space, using a static background grid, and a Lagrangian material space, using moving particles. By introducing this Eulerian-Lagrangian representation, we are able to naturally evolve and accumulate particle features using flow velocities generated from a generalized, high-dimensional force field. We demonstrate the efficacy of this system by solving various point cloud classification and segmentation problems with state-of-the-art performance. The entire geometric reservoir and data flow mimics the pipeline of the classic PIC/FLIP scheme in modeling natural flow, bridging the disciplines of geometric machine learning and physical simulation.


title: Improving Domain-Adapted Sentiment Classification by Deep Adversarial Mutual Learning

abstract: Domain-adapted sentiment classification refers to training on a labeled source domain to well infer document-level sentiment on an unlabeled target domain. Most existing relevant models involve a feature extractor and a sentiment classifier, where the feature extractor works towards learning domain-invariant features from both domains, and the sentiment classifier is trained only on the source domain to guide the feature extractor. As such, they lack a mechanism to use sentiment polarity lying in the target domain. To improve domain-adapted sentiment classification by learning sentiment from the target domain as well, we devise a novel deep adversarial mutual learning approach involving two groups of feature extractors, domain discriminators, sentiment classifiers, and label probers. The domain discriminators enable the feature extractors to obtain domain-invariant features. Meanwhile, the label prober in each group explores document sentiment polarity of the target domain through the sentiment prediction generated by the classifier in the peer group, and guides the learning of the feature extractor in its own group. The proposed approach achieves the mutual learning of the two groups in an end-to-end manner. Experiments on multiple public datasets indicate our method obtains the state-of-the-art performance, validating the effectiveness of mutual learning through label probers.


title: On the Consistency of Optimal Bayesian Feature Selection in the Presence of Correlations

abstract: Optimal Bayesian feature selection (OBFS) is a multivariate supervised screening method designed from the ground up for biomarker discovery. In this work, we prove that Gaussian OBFS is strongly consistent under mild conditions, and provide rates of convergence for key posteriors in the framework. These results are of enormous importance, since they identify precisely what features are selected by OBFS asymptotically, characterize the relative rates of convergence for posteriors on different types of features, provide conditions that guarantee convergence, justify the use of OBFS when its internal assumptions are invalid, and set the stage for understanding the asymptotic behavior of other algorithms based on the OBFS framework.


title: 28 GHz mmWave Channel Measurements: A Comparison of Horn and Phased Array Antennas and Coverage Enhancement Using Passive and Active Repeaters

abstract: Propagation channel characteristics are substantially different in sub-6 GHz and millimeter-wave (mmWave) bands. A typical mmWave link experiences more than an order-of-magnitude larger path loss and is more susceptible to blockages than a traditional sub-6 GHz link. Therefore, extensive measurement campaigns with efficient channel sounders are needed for a complete understanding of the mmWave propagation channel characteristics. In addition, solutions should be sought to overcome the outage problems associated with the mmWave bands and to achieve acceptable communication performance. In this paper, our 28 GHz channel sounder is introduced that can be used with horn antennas as well as phased array antennas. The two antenna types are then compared in terms of speed they switch from one direction to another and extracted multipath component (MPC) properties, such as path gain. We then propose the use of passive reflectors and an active repeater for coverage enhancement by improving the received signal strength. In particular, ECHO reflector and TURBO repeater are introduced, and measurement results using them are presented.


title: Multi-channel Acoustic Modeling using Mixed Bitrate OPUS Compression

abstract: Recent literature has shown that a learned front end with multi-channel audio input can outperform traditional beam-forming algorithms for automatic speech recognition (ASR). In this paper, we present our study on multi-channel acoustic modeling using OPUS compression with different bitrates for the different channels. We analyze the degradation in word error rate (WER) as a function of the audio encoding bitrate and show that the WER degrades by 12.6% relative with 16kpbs as compared to uncompressed audio. We show that its always preferable to have a multi-channel audio input over a single channel audio input given limited bandwidth. Our results show that for the best WER, when one of the two channels can be encoded with a bitrate higher than 32kbps, its optimal to encode the other channel with the highest bitrate possible. For bitrates lower than that, its preferable to distribute the bitrate equally between the two channels. We further show that by training the acoustic model on mixed bitrate input, up to 50% of the degradation can be recovered using a single model.


title: Model Extraction Attacks against Recurrent Neural Networks

abstract: Model extraction attacks are a kind of attacks in which an adversary obtains a new model, whose performance is equivalent to that of a target model, via query access to the target model efficiently, i.e., fewer datasets and computational resources than those of the target model. Existing works have dealt with only simple deep neural networks (DNNs), e.g., only three layers, as targets of model extraction attacks, and hence are not aware of the effectiveness of recurrent neural networks (RNNs) in dealing with time-series data. In this work, we shed light on the threats of model extraction attacks against RNNs. We discuss whether a model with a higher accuracy can be extracted with a simple RNN from a long short-term memory (LSTM), which is a more complicated and powerful RNN. Specifically, we tackle the following problems. First, in a case of a classification problem, such as image recognition, extraction of an RNN model without final outputs from an LSTM model is presented by utilizing outputs halfway through the sequence. Next, in a case of a regression problem. such as in weather forecasting, a new attack by newly configuring a loss function is presented. We conduct experiments on our model extraction attacks against an RNN and an LSTM trained with publicly available academic datasets. We then show that a model with a higher accuracy can be extracted efficiently, especially through configuring a loss function and a more complex architecture different from the target model.


title: Efficient q-Integer Linear Decomposition of Multivariate Polynomials

abstract: We present two new algorithms for the computation of the q-integer linear decomposition of a multivariate polynomial. Such a decomposition is essential for the treatment of q-hypergeometric symbolic summation via creative telescoping and for describing the q-counterpart of Ore-Sato theory. Both of our algorithms require only basic integer and polynomial arithmetic and work for any unique factorization domain containing the ring of integers. Complete complexity analyses are conducted for both our algorithms and two previous algorithms in the case of multivariate integer polynomials, showing that our algorithms have better theoretical performances. A Maple implementation is also included which suggests that our algorithms are also much faster in practice than previous algorithms.


title: Fully Learnable Front-End for Multi-Channel Acoustic Modeling using Semi-Supervised Learning

abstract: In this work, we investigated the teacher-student training paradigm to train a fully learnable multi-channel acoustic model for far-field automatic speech recognition (ASR). Using a large offline teacher model trained on beamformed audio, we trained a simpler multi-channel student acoustic model used in the speech recognition system. For the student, both multi-channel feature extraction layers and the higher classification layers were jointly trained using the logits from the teacher model. In our experiments, compared to a baseline model trained on about 600 hours of transcribed data, a relative word-error rate (WER) reduction of about 27.3% was achieved when using an additional 1800 hours of untranscribed data. We also investigated the benefit of pre-training the multi-channel front end to output the beamformed log-mel filter bank energies (LFBE) using L2 loss. We find that pre-training improves the word error rate by 10.7% when compared to a multi-channel model directly initialized with a beamformer and mel-filter bank coefficients for the front end. Finally, combining pre-training and teacher-student training produces a WER reduction of 31% compared to our baseline.


title: Instant single-pixel imaging: on-chip real-time implementation based on instant ghost imaging algorithm

abstract: Single-pixel imaging (SPI) uses a single-pixel detector to create an image of an object. SPI relies on a computer to construct an image, thus increasing both the size and cost of SPI and limiting its application. We developed instant single-pixel imaging (ISPI), an on-chip SPI system that implements real-time imaging at a rate of 25 fps. ISPI uses the instant ghost imaging algorithm we proposed which leverages signal differences for image creation. It does not require a computer, which greatly reduces its cost and size. The reconstruct time of ISPI for image creation is almost zero because little processing is required after signal detection. ISPI paves the way for the practical application of SPI.


title: Probabilistic direction-dependent ionospheric calibration for LOFAR-HBA

abstract: Direction dependent calibration and imaging is a vital part of producing deep, high fidelity, high-dynamic range radio images with a wide-field low-frequency array like LOFAR. Currently, state-of-the-art facet-based direction dependent calibration algorithms rely on the assumption that the isoplanatic-patch size is much larger than the separation between bright in-field calibrators. This assumption is often violated due to the dynamic nature of the ionosphere, and as a result direction dependent errors affect image quality between calibrators. In this paper we propose a probabilistic physics-informed model for inferring ionospheric phase screens, providing a calibration for all sources in the field of view. We apply our method to a randomly selected observation from the LOFAR Two-Metre Sky Survey archive, and show that almost all direction dependent effects between bright calibrators are corrected and that the root-mean-squared residuals around bright sources is reduced by 32\% on average.


title: Lens Flare: Magnified X-Ray Binaries as Passive Beacons in SETI

abstract: Low mass X-ray binaries (LMXBs) containing neutron stars are both extremely luminous and compact, emitting up to ~10^6 L_sun within a kilometer-scale boundary layer. This combination allows for easy modulation, motivating X-ray SETI. When X-ray lenses with radii 100--1,000 km magnify the LMXB boundary layer, it brightens by a factor of several thousand for a fraction of a second. In addition, there should be occultation events where the neutron star is blocked out. Passive X-ray lenses could require little internal power and the LMXB light source itself shines for millions of years, with potential for an effective beacon for interstellar communication. A very large number of lenses would be needed to ensure frequent signals in all directions, however, and gathering material to construct them could be very difficult. Avoiding collisions between lenses, aiming them, and building and maintaining their precise shapes pose additional challenges. "Lens flares" of bright LMXBs are easily detectable in the Galaxy, although they would be rare events, occurring perhaps once per decade. Our sensitive X-ray instruments could detect the eclipses of Galactic LMXBs and possibly intergalactic flares, but it is unlikely they would be observing the LMXB at the right time.


title: Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016

abstract: 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E >100 GeV) gamma-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ~ 3 times the Crab Nebula flux was measured by MAGIC on the 13th, 14th of June and 1st July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to high values of the Doppler factors (>=30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L_{jet}~10^{46} erg/s and under high values of the magnetic field strength (~100 G) and maximum proton energy (~few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016.


title: Higher Criticism Tuned Regression For Weak And Sparse Signals

abstract: Here we propose a novel searching scheme for a tuning parameter in high-dimensional penalized regression methods to address variable selection and modeling when sample sizes are limited compared to the data dimensions. Our method is motivated by high-throughput biological data such as genome-wide association studies (GWAS) and epigenome-wide association studies (EWAS). We propose a new estimate of the regularization parameter $\lambda$ in penalized regression methods based on an estimated lower bound of the proportion of false null hypotheses with confidence $(1 - \alpha)$. The bound is estimated by applying the empirical null distribution of the higher criticism statistic, a second-level significance test constructed by dependent $p$-values using a multi-split regression and aggregation method. A tuning parameter estimate in penalized regression, $\lambda$, corresponds with the lower bound of the proportion of false null hypotheses. Different penalized regression methods with varied signal sparsity and strength are compared in the multi-split method setting. We demonstrate the performance of our method using both simulation experiments and the applications of real data on (1) lipid-trait genetics from the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial and (2) epigenetic analysis evaluating smoking's influence in differential methylation in the Agricultural Lung Health Study. The proposed algorithm is included in the HCTR package, available at https://cran.r-project.org/web/packages/HCTR/index.html.


title: An optimum energy consumption hybrid algorithm for xln strategic design in wsns

abstract: In this paper, X-Layer protocol is originated which executes mobility error prediction (MEP) algorithm to calculate the remaining energy level of each node. This X-Layer protocol structure employs the mobility aware protocol that senses the mobility concerned to each node with the utilization of Ad-hoc On-Demand Distance Vector (AODV), which shares the information or data specific to the distance among individual nodes. With the help of this theory, the neighbour list will be updated only to those nodes which are mobile resulting in less energy consumption when compared to all (static/mobile) other nodes in the network. Apart from the MEP algorithm, clustering head (CH) election algorithm has also been specified to identify the relevant clusters whether they exists within the network region or not. Also clustering multi-hop routing (CMHR) algorithm was implemented in which the node can identify the cluster to which it belongs depending upon the distance from each cluster surrounding the node. Finally comprising the AODV routing protocol with the Two-Ray Ground method, we implement X-Layer protocol structure by considering MAC protocol in accordance to IEEE 802.15.4 to obtain the best results in energy consumption and also by reducing the energy wastage with respect to each node. The effective results had been illustrated through Network Simulator-II platform.


title: Superconducting Quantum Interference Devices Made of Sb-doped Bi2Se3 Topological Insulator Nanoribbons

abstract: We report the fabrication and characterization of superconducting quantum interference devices (SQUIDs) made of Sb-doped Bi2Se3 topological insulator (TI) nanoribbon (NR) contacted with PbIn superconducting electrodes. When an external magnetic field was applied along the NR axis, the TI NR exhibited periodic magneto-conductance oscillations, the so-called Aharonov-Bohm oscillations, owing to one-dimensional subbands. Below the superconducting transition temperature of PbIn electrodes, we observed supercurrent flow through TI NR-based SQUID. The critical current periodically modulates with a magnetic field perpendicular to the SQUID loop, revealing that the periodicity corresponds to the superconducting flux quantum. Our experimental observations can be useful to explore Majorana bound states (MBS) in TI NR, promising for developing topological quantum information devices.


title: Global Texture Enhancement for Fake Face Detection in the Wild

abstract: Generative Adversarial Networks (GANs) can generate realistic fake face images that can easily fool human beings.On the contrary, a common Convolutional Neural Network(CNN) discriminator can achieve more than 99.9% accuracyin discerning fake/real images. In this paper, we conduct an empirical study on fake/real faces, and have two important observations: firstly, the texture of fake faces is substantially different from real ones; secondly, global texture statistics are more robust to image editing and transferable to fake faces from different GANs and datasets. Motivated by the above observations, we propose a new architecture coined as Gram-Net, which leverages global image texture representations for robust fake image detection. Experimental results on several datasets demonstrate that our Gram-Net outperforms existing approaches. Especially, our Gram-Netis more robust to image editings, e.g. down-sampling, JPEG compression, blur, and noise. More importantly, our Gram-Net generalizes significantly better in detecting fake faces from GAN models not seen in the training phase and can perform decently in detecting fake natural images.


title: Rational solutions of the defocusing nonlocal nonlinear Schrodinger equation: Asymptotic analysis and soliton interactions

abstract: In this paper, we obtain the Nth-order rational solutions for the defocusing nonlocal nonlinear Schrodinger equation by the Darboux transformation and some limit technique. Then, via an improved asymptotic analysis method relying on the balance between different algebraic terms, we derive the explicit expressions of all asymptotic solitons of the rational solutions with the order 1<=N<=4. It turns out that the asymptotic solitons are localized in the straight or algebraic curves, and the exact solutions approach the curved asymptotic solitons with a slower rate than the straight ones. Moreover, we find that all the rational solutions exhibit just five different types of soliton interactions, and the interacting solitons are divided into two halves with each having the same amplitudes. Particularly for the curved asymptotic solitons, there may exist a slight difference for their velocities between at t and -t with certain parametric condition. In addition, we reveal that the soliton interactions in the rational solutions with N>=2 are stronger than those in the exponential and exponential-and-rational solutions.


title: Random magnetic anisotropy driven transitions in layered perovskite LaSrCoO$_4$

abstract: Attempts to unravel the nature of magnetic ordering in LaSrCoO$_4$ (Co$^{3+}$), a compound intermediate between antiferromagnetic (AFM) La$_2$CoO$_4$ (Co$^{2+}$) and ferromagnetic (FM) Sr$_2$CoO$_4$ (Co$^{4+}$), have met with a limited success so far. In this report, the results of a thorough investigation of dc magnetization and ac susceptibility (ACS) in single-phase LaSrCoO$_4$ provide clinching evidence for a thermodynamic paramagnetic (PM) - ferromagnetic (FM) phase transition at T$_{c}$ = 220.5 K, followed at lower temperature (T$_{g}$ = 7.7 K) by a transition to the cluster spin glass (CSG) state. Analysis of the low-field Arrott plot isotherms, in the critical region near T$_{c}$, in terms of the Aharony-Pytte scaling equation of state clearly establishes that the PM-FM transition is basically driven by random magnetic anisotropy (RMA). For temperatures below $\approx$ 30 K, large enough RMA destroys long-range FM order by breaking up the infinite FM network into FM clusters of finite size and leads to the formation of a CSG state at temperatures T $\lesssim$ 8 K by promoting freezing of finite FM clusters in random orientations. Increasing strength of the single-ion magnetocrystalline anisotropy (and hence RMA) with decreasing temperature is taken to reflect an increase in the number of low-spin (LS) Co$^{3+}$ ions at the expense of that of high-spin (HS) Co$^{3+}$ ions. At intermediate temperatures (30 K $\lesssim T \lesssim$ 180 K), spin dynamics has contributions from the infinite FM network (fast relaxation governed by a single anisotropy energy barrier) and finite FM clusters (extremely slow stretched exponential relaxation due to hierarchical energy barriers).