title: Solvability of nonlocal systems related to peridynamics

abstract: In this work, we study the Dirichlet problem associated with a strongly coupled system of nonlocal equations. The system of equations comes from a linearization of a model of peridynamics, a nonlocal model of elasticity. It is a nonlocal analogue of the Navier-Lam\'e system of classical elasticity. The leading operator is an integro-differential operator characterized by a distinctive matrix kernel which is used to couple differences of components of a vector field. The paper's main contributions are proving well-posedness of the system of equations and demonstrating optimal local Sobolev regularity of solutions. We apply Hilbert space techniques for well-posedness. The result holds for systems associated with kernels that give rise to non-symmetric bilinear forms. The regularity result holds for systems with symmetric kernels that may be supported only on a cone. For some specific kernels associated energy spaces are shown to coincide with standard fractional Sobolev spaces.


title: Phantom Inflated Planets in Occurrence Rate Based Samples

abstract: The recently launched Transiting Exoplanet Survey Satellite (TESS) is expected to produce many new exoplanet discoveries which will be especially amenable to follow-up study. Assessments of the planet discovery yield of TESS, such as Sullivan et al. (2015) and Barclay et al. (2018), will be important for planning follow-up work. Analyzing these predicted planet samples, however, we find that giant planet radii derived from the current bulk transiting planet sample have been used at all potential orbits without accounting for the temperature dependence of radius inflation. These phantom inflated planets (PIPs) make up just over 1.1% of the Sullivan et al. (2015) predicted population and about 8% of the Barclay et al. (2018) sample. Similar population predictions for direct imaging studies should likewise take care not to include inflated planets at large separations as such planets will appear larger and, depending on assumed albedo, brighter and more easily detectable than physically possible. Despite their high false positive rate, giant planets are the testbeds for atmospheric characterization techniques and care should be taken to understand and account for potential contaminating factors in this population.


title: Learning what and where to attend

abstract: Most recent gains in visual recognition have originated from the inclusion of attention mechanisms in deep convolutional networks (DCNs). Because these networks are optimized for object recognition, they learn where to attend using only a weak form of supervision derived from image class labels. Here, we demonstrate the benefit of using stronger supervisory signals by teaching DCNs to attend to image regions that humans deem important for object recognition. We first describe a large-scale online experiment (ClickMe) used to supplement ImageNet with nearly half a million human-derived "top-down" attention maps. Using human psychophysics, we confirm that the identified top-down features from ClickMe are more diagnostic than "bottom-up" saliency features for rapid image categorization. As a proof of concept, we extend a state-of-the-art attention network and demonstrate that adding ClickMe supervision significantly improves its accuracy and yields visual features that are more interpretable and more similar to those used by human observers.


title: Kepler-503b: An Object at the Hydrogen Burning Mass Limit Orbiting a Subgiant Star

abstract: Using spectroscopic radial velocities with the APOGEE instrument and Gaia distance estimates, we demonstrate that Kepler-503b, currently considered a validated Kepler planet, is in fact a brown-dwarf/low-mass star in a nearly circular 7.2-day orbit around a subgiant star. Using a mass estimate for the primary star derived from stellar models, we derive a companion mass and radius of $0.075\pm0.003 \ M_{\odot}$ ($78.6\pm3.1 \ M_{Jup}$) and $0.099^{+0.006}_{-0.004}\ R_{\odot}$ ($0.96^{+0.06}_{-0.04}\ R_{Jup}$), respectively. Assuming the system is coeval, the evolutionary state of the primary indicates the age is $\sim6.7$ Gyr. Kepler-503b sits right at the hydrogen burning mass limit, straddling the boundary between brown dwarfs and very low-mass stars. More precise radial velocities and secondary eclipse spectroscopy with James Webb Space Telescope will provide improved measurements of the physical parameters and age of this important system to better constrain and understand the physics of these objects and their spectra. This system emphasizes the value of radial velocity observations to distinguish a genuine planet from astrophysical false positives, and is the first result from the SDSS-IV monitoring of Kepler planet candidates with the multi-object APOGEE instrument.


title: Caratheodory convergence and harmonic measure

abstract: We give several new characterizations of Caratheodory convergence of simply connected domains. We then investigate how different definitions of convergence generalize to the multiply-connected case.


title: Estimates for functionals of solutions to higher-order heat-type equations with random initial conditions

abstract: In the present paper we continue the investigation of solutions to higher-order heat-type equations with random initial conditions, which play the important role in many applied areas. We consider the random initial conditions given by harmonizable $\varphi$-sub-Gaussian processes. The main results are the bounds for the distributions of the suprema over bounded and unbounded domains for solutions of such equations. The results obtained in the paper hold, in particular, for the case of Gaussian initial condition.


title: Pulsed excitation dynamics in quantum dot-cavity systems: limits to optimizing the fidelity of on-demand single photon sources

abstract: A quantum dot coupled to an optical cavity has recently proven to be an excellent source of on-demand single photons. Typically, applications require simultaneous high efficiency of the source and quantum indistinguishability of the extracted photons. While much progress has been made both in suppressing background sources of decoherence and utilizing cavity-quantum electrodynamics to overcome fundamental limitations set by the intrinsic exciton-phonon scattering inherent in the solid-state platform, the role of the excitation pulse has been often neglected. We investigate quantitatively the factors associated with pulsed excitation that can limit simultaneous efficiency and indistinguishability, including excitation of multiple excitons, multi-photons, and pump-induced dephasing, and find for realistic single photon sources that these effects cause degradation of the source figures-of-merit comparable to that of phonon scattering. We also develop rigorous open quantum system polaron master equation models of quantum dot dynamics under a time-dependent drive which incorporate non-Markovian effects of both photon and phonon reservoirs, and explicitly show how coupling to a high Q-factor cavity suppresses multi-photon emission in a way not predicted by commonly employed models. We then use our findings to summarize the criteria that can be used for single photon source optimization.


title: Double hysteresis loops in proper uniaxial ferroelectrics

abstract: For the first time in a bulk proper uniaxial ferroelectrics, double antiferroelectric-like hysteresis loops have been observed in the case of Sn$_2$P$_2$S$_6$ crystal. The quantum anharmonic oscillator model was proposed for description of such polarization switching process. This phenomenon is related to three-well local potential of spontaneous polarization fluctuations at peculiar negative ratio of coupling constants which correspond to inter-site interaction in given sublattice and interaction between two sublattices of Sn$_2$P$_2$S$_6$ modeled crystal structure. Obtained data can be used for development of triple-level cell type memory technology.


title: Chirality in Gravitational and Electromagnetic Interactions with Matter

abstract: It has been suggested that single and double jets observed emanating from certain astrophysical objects may have a purely gravitational origin. We discuss new classes of pulsed gravitational wave solutions to the equation for perturbations of Ricci-flat spacetimes around Minkowski metrics, as models for the genesis of such phenomena. We discuss how these solutions are motivated by the analytic structure of spatially compact finite energy pulse solutions of the source-free Maxwell equations generated from complex chiral eigenmodes of a chirality operator. Complex gravitational pulse solutions to the linearised source-free Einstein equations are classified in terms of their chirality and generate a family of non-stationary real spacetime metrics. Particular members of these families are used as backgrounds in analysing time-like solutions to the geodesic equation for test particles. They are found numerically to exhibit both single and double jet-like features with dimensionless aspect ratios suggesting that it may be profitable to include such backgrounds in simulations of astrophysical jet dynamics from rotating accretion discs involving electromagnetic fields.


title: Rapid seismic domain transfer: Seismic velocity inversion and modeling using deep generative neural networks

abstract: Traditional physics-based approaches to infer sub-surface properties such as full-waveform inversion or reflectivity inversion are time-consuming and computationally expensive. We present a deep-learning technique that eliminates the need for these computationally complex methods by posing the problem as one of domain transfer. Our solution is based on a deep convolutional generative adversarial network and dramatically reduces computation time. Training based on two different types of synthetic data produced a neural network that generates realistic velocity models when applied to a real dataset. The system's ability to generalize means it is robust against the inherent occurrence of velocity errors and artifacts in both training and test datasets.


title: Frequency linewidth and decay length of spin waves in curved magnetic membranes

abstract: The curvature of a magnetic membrane was presented as a means of inducing nonreciprocities in the spin-wave (SW) dispersion relation (see [Ot\'alora et al. Phys. Rev. Lett., 2016 117, 227203] and [Ot\'alora et al. Phys. Rev. B., 2017 95, 184415]), thereby expanding the toolbox for controlling SWs. In this paper, we further complement this toolbox by analytically showing that the membrane curvature is also manifested in the absorption of SWs, leading to a difference in the frequency linewidth (or lifetime) of counterpropagating magnons. Herein, we studied the nanotubular case, predicting changes of approximately greater than 10% and up to 20% in the frequency linewidth of counterpropagating SWs for a wide range of nanotube radii ranging from 30 nm to 260 nm and with a thickness of 10 nm. These percentages are comparable to those that can be extracted from experiments on heavy metal/magnetic metal sandwiches, wherein linewidth asymmetry results from an interfacial Dzyaloshinskii-Moriya interaction (DMI). We also show that the interplay between the frequency linewidth and group velocity leads to asymmetries in the SW decay length, presenting changes between 10\% and 22\% for counterpropagating SWs in the frequency range of 2 - 10 GHz. For the case of the SW dispersion relation, the predicted effects are identified as the classical dipole-dipole interaction, and the analytical expression of the frequency linewidth has the same mathematical form as in thin films with the DMI. Furthermore, we present limiting cases of a tubular geometry with negligible curvature such that our analytical model converges to the case of a planar thin film known from the literature. Our findings represent a step forward toward the realization of three-dimensional curvilinear magnonic devices.


title: An algebraic model for rational toral G-spectra

abstract: For G a compact Lie group, toral G-spectra are those rational G-spectra whose geometric isotropy consists of subgroups of a maximal torus of G. The homotopy category of rational toral G-spectra is a retract of the category of all rational G-spectra. In this paper we show that the abelian category defined by the second author in previous work gives an algebraic model for the toral part of rational G-spectra. This is a major step in establishing an algebraic model for all rational G-spectra, for any compact Lie group G.


title: Aperiodic points in $\mathbb Z^2$-subshifts

abstract: We consider the structure of aperiodic points in $\mathbb Z^2$-subshifts, and in particular the positions at which they fail to be periodic. We prove that if a $\mathbb Z^2$-subshift contains points whose smallest period is arbitrarily large, then it contains an aperiodic point. This lets us characterise the computational difficulty of deciding if an $\mathbb Z^2$-subshift of finite type contains an aperiodic point. Another consequence is that $\mathbb Z^2$-subshifts with no aperiodic point have a very strong dynamical structure and are almost topologically conjugate to some $\mathbb Z$-subshift. Finally, we use this result to characterize sets of possible slopes of periodicity for $\mathbb Z^3$-subshifts of finite type.


title: Stein operators for variables form the third and fourth Wiener chaoses

abstract: Let $Z$ be a standard normal random variable and let $H_n$ denote the $n$-th Hermite polynomial. In this note, we obtain Stein equations for the random variables $H_3(Z)$ and $H_4(Z)$, which represents a first step towards developing Stein's method for distributional limits from the third and fourth Wiener chaoses. Perhaps surprisingly, these Stein equations are fifth and third order linear ordinary differential equations, respectively. As a warm up, we obtain a Stein equation for the random variable $aZ^2+bZ+c$, $a,b,c\in\mathbb{R}$, which leads us to a Stein equation for the non-central chi-square distribution. We also provide a discussion as to why obtaining Stein equations for $H_n(Z)$, $n\geq5$, is more challenging.


title: One machine, one minute, three billion tetrahedra

abstract: This paper presents a new scalable parallelization scheme to generate the 3D Delaunay triangulation of a given set of points. Our first contribution is an efficient serial implementation of the incremental Delaunay insertion algorithm. A simple dedicated data structure, an efficient sorting of the points and the optimization of the insertion algorithm have permitted to accelerate reference implementations by a factor three. Our second contribution is a multi-threaded version of the Delaunay kernel that is able to concurrently insert vertices. Moore curve coordinates are used to partition the point set, avoiding heavy synchronization overheads. Conflicts are managed by modifying the partitions with a simple rescaling of the space-filling curve. The performances of our implementation have been measured on three different processors, an Intel core-i7, an Intel Xeon Phi and an AMD EPYC, on which we have been able to compute 3 billion tetrahedra in 53 seconds. This corresponds to a generation rate of over 55 million tetrahedra per second. We finally show how this very efficient parallel Delaunay triangulation can be integrated in a Delaunay refinement mesh generator which takes as input the triangulated surface boundary of the volume to mesh.


title: The Impact of Uncle Rewards on Selfish Mining in Ethereum

abstract: Many of today's crypto currencies use blockchains as decentralized ledgers and secure them with proof of work. In case of a fork of the chain, Bitcoin's rule for achieving consensus is selecting the longest chain and discarding the other chain as stale. It has been demonstrated that this consensus rule has a weakness against selfish mining in which the selfish miner exploits the variance in block generation by partially withholding blocks. In Ethereum, however, under certain conditions stale blocks don't have to be discarded but can be referenced from the main chain as uncle blocks yielding a partial reward. This concept limits the impact of network delays on the expected revenue for miners. But the concept also reduces the risk for a selfish miner to gain no rewards from withholding a freshly minted block. This paper uses a Monte Carlo simulation to quantify the effect of uncle blocks both to the profitability of selfish mining and the blockchain's security in Ethereum (ETH). A brief outlook about a recent Ethereum Classic (ETC) improvement proposal that weighs uncle blocks during the selection of the main chain will be given.


title: Deep Barcodes for Fast Retrieval of Histopathology Scans

abstract: We investigate the concept of deep barcodes and propose two methods to generate them in order to expedite the process of classification and retrieval of histopathology images. Since binary search is computationally less expensive, in terms of both speed and storage, deep barcodes could be useful when dealing with big data retrieval. Our experiments use the dataset Kimia Path24 to test three pre-trained networks for image retrieval. The dataset consists of 27,055 training images in 24 different classes with large variability, and 1,325 test images for testing. Apart from the high-speed and efficiency, results show a surprising retrieval accuracy of 71.62% for deep barcodes, as compared to 68.91% for deep features and 68.53% for compressed deep features.


title: Asymptotic symmetries of three-dimensional Chern-Simons gravity for the Maxwell algebra

abstract: We study a three-dimensional Chern-Simons gravity theory based on the Maxwell algebra. We find that the boundary dynamics is described by an enlargement and deformation of the $\mathfrak{bms}_3$ algebra with three independent central charges. This symmetry arises from a gravity action invariant under the local Maxwell group and is characterized by presence of Abelian generators which modify the commutation relations of the super-translations in the standard $\mathfrak{bms}_3$ algebra. Our analysis is based on the charge algebra of the theory in the BMS gauge, which includes the known solutions of standard asymptotically flat case. The field content of the theory is different than the one of General Relativity, but it includes all its geometries as particular solutions. In this line, we also study the stationary solutions of the theory in ADM form and we show that the vacuum energy and the vacuum angular momentum of the stationary configuration are influenced by the presence of the gravitational Maxwell field.


title: Revealing Color Forces with Transverse Polarized Electron Scattering

abstract: The Spin Asymmetries of the Nucleon Experiment (SANE) measured two double spin asymmetries using a polarized proton target and polarized electron beam at two beam energies, 4.7 GeV and 5.9 GeV. A large-acceptance open-configuration detector package identified scattered electrons at 40$^{\circ}$ and covered a wide range in Bjorken $x$ ($0.3 < x < 0.8$). Proportional to an average color Lorentz force, the twist-3 matrix element, $\tilde{d}_2^p$, was extracted from the measured asymmetries at $Q^2$ values ranging from 2.0 to 6.0 GeV$^2$. The data display the opposite sign compared to most quark models, including the lattice QCD result, and an apparently unexpected scale dependence. Furthermore when combined with the neutron data in the same $Q^2$ range the results suggest a flavor independent average color Lorentz force.


title: Nonparametric Density Estimation under Adversarial Losses

abstract: We study minimax convergence rates of nonparametric density estimation under a large class of loss functions called "adversarial losses", which, besides classical $\mathcal{L}^p$ losses, includes maximum mean discrepancy (MMD), Wasserstein distance, and total variation distance. These losses are closely related to the losses encoded by discriminator networks in generative adversarial networks (GANs). In a general framework, we study how the choice of loss and the assumed smoothness of the underlying density together determine the minimax rate. We also discuss implications for training GANs based on deep ReLU networks, and more general connections to learning implicit generative models in a minimax statistical sense.


title: Quantum classification of the MNIST dataset with Slow Feature Analysis

abstract: Quantum machine learning carries the promise to revolutionize information and communication technologies. While a number of quantum algorithms with potential exponential speedups have been proposed already, it is quite difficult to provide convincing evidence that quantum computers with quantum memories will be in fact useful to solve real-world problems. Our work makes considerable progress towards this goal. We design quantum techniques for Dimensionality Reduction and for Classification, and combine them to provide an efficient and high accuracy quantum classifier that we test on the MNIST dataset. More precisely, we propose a quantum version of Slow Feature Analysis (QSFA), a dimensionality reduction technique that maps the dataset in a lower dimensional space where we can apply a novel quantum classification procedure, the Quantum Frobenius Distance (QFD). We simulate the quantum classifier (including errors) and show that it can provide classification of the MNIST handwritten digit dataset, a widely used dataset for benchmarking classification algorithms, with $98.5\%$ accuracy, similar to the classical case. The running time of the quantum classifier is polylogarithmic in the dimension and number of data points. We also provide evidence that the other parameters on which the running time depends (condition number, Frobenius norm, error threshold, etc.) scale favorably in practice, thus ascertaining the efficiency of our algorithm.


title: Clustering - What Both Theoreticians and Practitioners are Doing Wrong

abstract: Unsupervised learning is widely recognized as one of the most important challenges facing machine learning nowa- days. However, in spite of hundreds of papers on the topic being published every year, current theoretical understanding and practical implementations of such tasks, in particular of clustering, is very rudimentary. This note focuses on clustering. I claim that the most signif- icant challenge for clustering is model selection. In contrast with other common computational tasks, for clustering, dif- ferent algorithms often yield drastically different outcomes. Therefore, the choice of a clustering algorithm, and their pa- rameters (like the number of clusters) may play a crucial role in the usefulness of an output clustering solution. However, currently there exists no methodical guidance for clustering tool-selection for a given clustering task. Practitioners pick the algorithms they use without awareness to the implications of their choices and the vast majority of theory of clustering papers focus on providing savings to the resources needed to solve optimization problems that arise from picking some concrete clustering objective. Saving that pale in com- parison to the costs of mismatch between those objectives and the intended use of clustering results. I argue the severity of this problem and describe some recent proposals aiming to address this crucial lacuna.


title: Codimension-two defects and Argyres-Douglas theories from outer-automorphism twist in 6d $(2,0)$ theories

abstract: The 6d $(2,0)$ theory has codimension-one symmetry defects associated to the outer-automorphism group of the underlying ADE Lie algebra. These symmetry defects give rise to twisted sectors of codimension-two defects that are either regular or irregular corresponding to simple or higher order poles of the Higgs field. In this paper, we perform a systematic study of twisted irregular codimension-two defects generalizing our earlier work in the untwisted case. In a class S setup, such twisted defects engineer 4d ${\mathcal N}=2$ superconformal field theories of the Argyres-Douglas type whose flavor symmetries are (subgroups of) non-simply-laced Lie groups. We propose formulae for the conformal and flavor central charges of these twisted theories, accompanied by nontrivial consistency checks. We also identify the 2d chiral algebra (vertex operator algebra) of a subclass of these theories and determine their Higgs branch moduli space from the associated variety of the chiral algebra.


title: Tiling the plane with equilateral triangles

abstract: Let $\cal T$ be a tiling of the plane with equilateral triangles no two of which share a side. We prove that if the side lengths of the triangles are bounded from below by a positive constant, then $\cal T$ is periodic and it consists of translates of only at most three different triangles. As a corollary, we prove a theorem of Scherer and answer a question of Nandakumar. The same result has been obtained independently by Richter and Wirth.


title: Distribution Matching Losses Can Hallucinate Features in Medical Image Translation

abstract: This paper discusses how distribution matching losses, such as those used in CycleGAN, when used to synthesize medical images can lead to mis-diagnosis of medical conditions. It seems appealing to use these new image synthesis methods for translating images from a source to a target domain because they can produce high quality images and some even do not require paired data. However, the basis of how these image translation models work is through matching the translation output to the distribution of the target domain. This can cause an issue when the data provided in the target domain has an over or under representation of some classes (e.g. healthy or sick). When the output of an algorithm is a transformed image there are uncertainties whether all known and unknown class labels have been preserved or changed. Therefore, we recommend that these translated images should not be used for direct interpretation (e.g. by doctors) because they may lead to misdiagnosis of patients based on hallucinated image features by an algorithm that matches a distribution. However there are many recent papers that seem as though this is the goal.


title: Compiling with Continuations and LLVM

abstract: LLVM is an infrastructure for code generation and low-level optimizations, which has been gaining popularity as a backend for both research and industrial compilers, including many compilers for functional languages. While LLVM provides a relatively easy path to high-quality native code, its design is based on a traditional runtime model which is not well suited to alternative compilation strategies used in high-level language compilers, such as the use of heap-allocated continuation closures. This paper describes a new LLVM-based backend that supports heap-allocated continuation closures, which enables constant-time callcc and very-lightweight multithreading. The backend has been implemented in the Parallel ML compiler, which is part of the Manticore system, but the results should be useful for other compilers, such as Standard ML of New Jersey, that use heap-allocated continuation closures.


title: Time Dependent Radiation Driven Winds

abstract: We study temporal variability of radiation driven winds using one dimensional, time dependent simulations and an extension of the classic theory of line driven winds developed by Castor Abbott and Klein. We drive the wind with a sinusoidally varying radiation field and find that after a relaxation time, determined by the propagation time for waves to move out of the acceleration zone of the wind, the solution settles into a periodic state. Winds driven at frequencies much higher than the dynamical frequency behave like stationary winds with time averaged radiation flux whereas winds driven at much lower frequencies oscillate between the high and low flux stationary states. Most interestingly, we find a resonance frequency near the dynamical frequency which results in velocity being enhanced or suppressed by a factor comparable to the amplitude of the flux variation. Whether the velocity is enhanced or suppressed depends on the relative phase between the radiation and the dynamical variables. These results suggest that a time-varying radiation source can induce density and velocity perturbations in the acceleration zones of line driven winds.


title: Stellar Companions of Exoplanet Host Stars in K2

abstract: It is well established that roughly half of all nearby solar-type stars have at least one companion. Stellar companions can have significant implications for the detection and characterization of exoplanets, including triggering false positives and masking the true radii of planets. Determining the fraction of exoplanet host stars that are also binaries allows us to better determine planetary characteristics as well as establish the relationship between binarity and planet formation. Using high angular resolution speckle imaging, we detect stellar companions within ~1 arcsec of K2 planet candidate host stars. Comparing our detected companion rate to TRILEGAL star count simulations and known detection limits of speckle imaging we estimate the binary fraction of K2 planet host stars to be 40 - 50%, similar to that of Kepler exoplanet hosts and field stars.


title: Counterfactual Mean Embeddings

abstract: Counterfactual inference has become a ubiquitous tool in online advertisement, recommendation systems, medical diagnosis, and econometrics. Accurate modeling of outcome distributions associated with different interventions -- known as counterfactual distributions -- is crucial for the success of these applications. In this work, we propose to model counterfactual distributions using a novel Hilbert space representation called counterfactual mean embedding (CME). The CME embeds the associated counterfactual distribution into a reproducing kernel Hilbert space (RKHS) endowed with a positive definite kernel, which allows us to perform causal inference over the entire landscape of the counterfactual distribution. Based on this representation, we propose a distributional treatment effect (DTE) that can quantify the causal effect over entire outcome distributions. Our approach is nonparametric as the CME can be estimated under the unconfoundedness assumption from observational data without requiring any parametric assumption about the underlying distributions. We also establish a rate of convergence of the proposed estimator which depends on the smoothness of the conditional mean and the Radon-Nikodym derivative of the underlying marginal distributions. Furthermore, our framework allows for more complex outcomes such as images, sequences, and graphs. Our experimental results on synthetic data and off-policy evaluation tasks demonstrate the advantages of the proposed estimator.


title: CUDACLAW: A high-performance programmable GPU framework for the solution of hyperbolic PDEs

abstract: We present cudaclaw, a CUDA-based high performance data-parallel framework for the solution of multidimensional hyperbolic partial differential equation (PDE) systems, equations describing wave motion. cudaclaw allows computational scientists to solve such systems on GPUs without being burdened by the need to write CUDA code, worry about thread and block details, data layout, and data movement between the different levels of the memory hierarchy. The user defines the set of PDEs to be solved via a CUDA- independent serial Riemann solver and the framework takes care of orchestrating the computations and data transfers to maximize arithmetic throughput. cudaclaw treats the different spatial dimensions separately to allow suitable block sizes and dimensions to be used in the different directions, and includes a number of optimizations to minimize access to global memory.


title: The Ridge Effect, Azimuthal Correlations, and other Novel Features of Gluonic String Collisions in High Energy Photon-Mediated Reactions

abstract: One of the remarkable features of high-multiplicity hadronic events in proton-proton collisions at the LHC is the fact that the produced particles appear as two "ridges", opposite in azimuthal angle $\phi$, with approximately flat rapidity distributions. This phenomena can be identified with the inelastic collision of gluonic flux tubes associated with the QCD interactions responsible for quark confinement in hadrons. In this paper we analyze the ridge phenomena when the collision involves a flux tube connecting the quark and antiquark of a high energy real or virtual photon. We discuss gluonic tube string collisions in the context of two examples: electron-proton scattering at a future electron-ion collider or the peripheral scattering of protons accessible at the LHC. A striking prediction of our analysis is that the azimuthal angle of the produced ridges will be correlated with the scattering plane of the electron or proton producing the virtual photon. In the case of $ep \to eX$, the final state $X$ is expected to exhibit maximal multiplicity when the elliptic flow in $X$ is aligned with the electron scattering plane. In the $pp \to ppX$ example, the multiplicity and elliptic flow in $X$ are estimated to exhibit correlated oscillations as functions of the azimuthal angle $\Phi$ between the proton scattering planes. In the minimum-bias event samples, the amplitude of oscillations is expected to be on the order of 2 to 4 percent of the mean values. In the events with highest multiplicity, the oscillations can be three times larger than in the minimum-bias event samples.


title: Isogeometric Simulation of Thermal Expansion for Twin Screw Compressors

abstract: Isogeometric Analysis (IGA) is a recently introduced computational approach intended to breach the gap between the Finite Element Analysis and the Computer Aided Design worlds. In this work, we apply it to numerically simulate thermal expansion of oil-free twin screw compressors in operation. High global smoothness of IGA leads to a more accurate representation of the compressor geometry. We utilize standard tri-variate B-splines to parametrize the rotors, while the casing is modeled exactly by using NURBS. We employ the Galerkin version of IGA to solve the thermal expansion problem in the stationary case. The results allow to estimate the contraction of the clearance space between the casing and the rotors. The implementation is based on the open source C++ library G+Smo. This work is supported by the European Union within the Project MOTOR: Multi-ObjecTive design Optimization of fluid eneRgy machines.


title: Optimal design of deterministic lateral displacement device for viscosity contrast based cell sorting

abstract: We solve a design optimization problem for deterministic lateral displacement (DLD) device to sort same-size biological cells by their deformability, in particular to sort red blood cells (RBCs) by their viscosity contrast between the fluid in the interior and the exterior of the cells. A DLD device optimized for efficient cell sorting enables rapid medical diagnoses of several diseases such as malaria since infected cells are stiffer than their healthy counterparts. The device consists of pillar arrays in which pillar rows are tilted and hence are not orthogonal to the columns. This arrangement leads cells to have different final vertical displacements depending on their deformability, therefore, it vertically separates the cells. Pillar cross section, tilt angle of the pillar rows and center-to-center distances between pillars define a unique device. For a given pair of viscosity contrast values of the cells we seek optimal DLD designs by fixing the tilt angle and the center-to-center distances. So the only design parameter is the pillar cross section which we parameterize with uniform 5th order B-splines. We propose an objective function to try to capture efficient cell sorting. The objective function is evaluated by simulating the cell flows through a device using our 2D model based on a boundary integral method (Kabacaoglu et al. Journal of Computational Physics, 357:43-77, 2018). We solve the optimization problem using the covariance matrix adaptation evolution strategy (CMA-ES), which is a stochastic, derivative-free algorithm. We present several scenarios where solving the optimization problem finds designs that can separate cells with similar viscosity contrast values. To the best of our knowledge, this is the first study which poses designing a DLD device as a constrained optimization problem and shows that solving this problem systematically discovers optimal designs.


title: Crossover from the Ultracold to the Quasiclassical Regime in State-Selected Photodissociation

abstract: Processes that break molecular bonds are typically observed with molecules occupying a mixture of quantum states and successfully described with quasiclassical models, while a few studies have explored the distinctly quantum mechanical low-energy regime. Here we use photodissociation of diatomic strontium molecules to demonstrate the crossover from the ultracold, quantum regime where the photofragment angular distributions strongly depend on the kinetic energy, to the energy-independent quasiclassical regime. Using time-of-flight velocity map imaging for photodissociation channels with millikelvin reaction barriers, we explore photofragment energies in the 0.1-300 mK range experimentally and up to 3 K theoretically, and discuss the energy scale at which the crossover occurs. Furthermore, we find that the effects of quantum statistics can unexpectedly persist to high photodissociation energies.


title: Non-invariance of weak approximation properties under extension of the ground field

abstract: For rational points on algebraic varieties defined over a number field $K$, we study the behavior of the property of weak approximation with Brauer-Manin obstruction under extension of the ground field. We construct K-varieties accompanied with a quadratic extension $L/K$ such that the property holds over $K$ (conditional on a conjecture) while fails over $L$. The result is unconditional when $K = \mathbb{Q}$ or $K$ is one of several quadratic number fields. Over $\mathbb{Q}$, we give an explicit example.


title: Torques on low-mass bodies in retrograde orbit in gaseous disks

abstract: We evaluate the torque acting on a gravitational perturber on a retrograde circular orbit in the midplane of a gaseous disk. We assume that the mass of this satellite is so low it weakly disturbs the disk (type I migration). The perturber may represent the companion of a binary system with a small mass ratio. We compare the results of hydrodynamical simulations with analytic predictions. Our two-dimensional (2D) simulations indicate that the torque acting on a perturber with softening radius $R_{\rm soft}$ can be accounted for by a scattering approach if $R_{\rm soft}<0.3H$, where $H$ is defined as the ratio between the sound speed and the angular velocity at the orbital radius of the perturber. For $R_{\rm soft}>0.3H$, the torque may present large and persistent oscillations, but the resultant time-averaged torque decreases rapidly with increasing $R_{\rm soft}/H$, in agreement with previous analytical studies. We then focus on the torque acting on small-size perturbers embedded in full three-dimensional (3D) disks and argue that the density waves propagating at distances $\lesssim H$ from the perturber contribute significantly to the torque because they transport angular momentum. We find a good agreement between the torque found in 3D simulations and analytical estimates based on ballistic orbits. We compare the radial migration timescales of prograde versus retrograde perturbers. For a certain range of the perturber's mass and aspect ratio of the disk, the radial migration timescale in the retrograde case may be appreciably shorter than in the prograde case. We also provide the smoothing length required in 2D simulations in order to account for 3D effects.


title: Analysis of the energy dissipation laws in multi-component phase field models

abstract: In this paper, two approaches for modeling three-component fluid flows using diffusive interface method are discussed. Thermodynamic consistency of the proposed models is preserved when using an energetic variational framework to derive the coupled systems of partial differential equations that comprise the resulting models. The issue of algebraic and dynamic consistency is investigated. In addition, the two approaches that are presented are compared analytically and numerically.


title: Advection and diffusion in a chemically induced compressible flow

abstract: We study analytically the joint dispersion of Gaussian patches of salt and colloids in linear flows, and how salt gradients accelerate or delay colloid spreading by diffusiophoretic effects. Because these flows have constant gradients in space, the problem can be solved almost entirely for any set of parameters, leading to predictions of how the mixing time and the Batchelor scale are modified by diffusiophoresis. We observe that the evolution of global concentrations, defined as the inverse of the patches areas, are very similar to those obtained experimentally in chaotic advection. They are quantitatively explained by examining the area dilatation factor, in which diffusive and diffusiophoretic effects are shown to be additive and appear as the divergence of a diffusive contribution or of a drift velocity. An analysis based on compressibility is developed in the salt-attracting case, for which colloids are first compressed before dispersion, to predict the maximal colloid concentration as a function of the parameters. This maximum is found not to depend on the flow stretching rate nor on its topology (strain or shear flow), but only on the characteristics of salt and colloids (diffusion coefficients and diffusiophoretic constant) and the initial size of the patches.


title: Rate-Optimal Denoising with Deep Neural Networks

abstract: Deep neural networks provide state-of-the-art performance for image denoising, where the goal is to recover a near noise-free image from a noisy observation. The underlying principle is that neural networks trained on large datasets have empirically been shown to be able to generate natural images well from a low-dimensional latent representation of the image. Given such a generator network, a noisy image can be denoised by i) finding the closest image in the range of the generator or by ii) passing it through an encoder-generator architecture (known as an autoencoder). However, there is little theory to justify this success, let alone to predict the denoising performance as a function of the network parameters. In this paper we consider the problem of denoising an image from additive Gaussian noise using the two generator based approaches. In both cases, we assume the image is well described by a deep neural network with ReLU activations functions, mapping a $k$-dimensional code to an $n$-dimensional image. In the case of the autoencoder, we show that the feedforward network reduces noise energy by a factor of $O(k/n)$. In the case of optimizing over the range of a generative model, we state and analyze a simple gradient algorithm that minimizes a non-convex loss function, and provably reduces noise energy by a factor of $O(k/n)$. We also demonstrate in numerical experiments that this denoising performance is, indeed, achieved by generative priors learned from data.


title: The global oscillations of the sun from tsi, soho and sdo data

abstract: The global oscillations of the Sun are investigated on the base of three independent sets of data: 1. Records of Total Solar Irradiance (TSI). 2. The average brightness of the photosphere in MDI (SOHO) and HMI (SDO) images. 3. Records of the time variations of sunlight reflected from the planets during their passage in the field of view of the LASCO C3 coronagraph. It is found that in the low-frequency spectrum of solar oscillations as a star there are modes of 8-10, 36-38 hours and the main stable 9-day mode. The last mode can be regarded as a reflection of oscillatory processes near the fast rotating core of the Sun according to Fossat et al. (2017).


title: Thin Position through the lens of trisections of 4-manifolds

abstract: Motivated by M. Scharlemann and A. Thompson's definition of thin position of 3-manifolds, we define the width of a handle decomposition a 4-manifold and introduce the notion of thin position of a compact smooth 4-manifold. We determine all manifolds having width equal to $\{1,\dots, 1\}$, and give a relation between the width of $M$ and its double $M\cup_{id_\partial} \overline M$. In particular, we describe how to obtain genus $2g+2$ and $g+2$ trisection diagrams for sphere bundles over orientable and non-orientable surfaces of genus $g$, respectively. By last, we study the problem of describing relative handlebodies as cyclic covers of 4-space branched along knotted surfaces from the width perspective.


title: A Novel Exoplanetary Habitability Score via Particle Swarm Optimization of CES Production Functions

abstract: The search for life has two goals essentially: looking for planets with Earth-like conditions (Earth similarity) and looking for the possibility of life in some form (habitability). Determining habitability from exoplanet data requires that determining parameters are collectively considered before coming up with a conclusion as no single factor alone contributes to it. Our proposed models, would serve as an indicator while looking for new habitable worlds, if computed with precision and efficiency. The models are of the type constrained optimization, multivariate, convex but may suffer from curvature violation and premature convergence impacting desired habitability scores. We mitigate the problem by proposing modified Particle Swarm Optimization (PSO) to tackle constraints and ensuring global optima. In the process, a python library to tackle such problems has been created.


title: Quantum field theory with indefinite causal structure

abstract: Quantum field theory (QFT) in classical spacetime has revealed interesting and puzzling aspects about gravitational systems, in particular black hole thermodynamics and its information processing. Although quantum gravitational effects may be relevant for a better understanding of these topics, a commonly accepted framework for studying QFT with quantum gravitational effects is missing. We present a theory for studying quantum fields in the presence of quantum indefinite causal structure. This theory incorporates quantum properties of spacetime causal structure in a model independent way, and exposes universal features of quantum spacetime which are independent of the details about its microscopic degrees of freedom (strings, loops, causal sets...).


title: Statistical parameter inference of bacterial swimming strategies

abstract: We provide a detailed stochastic description of the swimming motion of an E.coli bacterium in two dimension, where we resolve tumble events in time. For this purpose, we set up two Langevin equations for the orientation angle and speed dynamics. Calculating moments, distribution and autocorrelation functions from both Langevin equations and matching them to the same quantities determined from data recorded in experiments, we infer the swimming parameters of E.coli . They are the tumble rate ${\lambda}$, the tumble time $r^{-1}$ , the swimming speed $v_0$ , the strength of speed fluctuations ${\sigma}$, the relative height of speed jumps ${\eta}$, the thermal value for the rotational diffusion coefficient $D_0$ , and the enhanced rotational diffusivity during tumbling $D_T$ . Conditioning the observables on the swimming direction relative to the gradient of a chemoattractant, we infer the chemotaxis strategies of E.coli . We confirm the classical strategy of a lower tumble rate for swimming up the gradient but also a smaller mean tumble angle (angle bias). The latter is realized by shorter tumbles as well as a slower diffusive reorientation. We also find that speed fluctuations are increased by about 30% when swimming up the gradient compared to the reversed direction.


title: Thermal Creep and the Skyrmion Hall Angle in Driven Skyrmion Crystals

abstract: We numerically examine thermal effects on the skyrmion Hall angle for driven skyrmions interacting with quenched disorder. We identify a creep regime in which motion occurs via intermittent jumps between pinned and flowing states. Here the skyrmion Hall angle is zero since the skyrmions have time to relax into equilibrium positions in the pinning sites, eliminating the side-jump motion induced by the Magnus force. At higher drives we find a crossover to a viscous flow regime where the skyrmion Hall angle is finite and increases with increasing drive or temperature. Our results are in agreement with recent experiments which also show a regime of finite skyrmion velocity with zero skyrmion Hall angle crossing over to a viscous flow regime with a skyrmion Hall angle that increases with drive.


title: Quantum probe of an on-chip broadband interferometer for quantum microwave photonics

abstract: Quantum microwave photonics aims at generating, routing, and manipulating propagating quantum microwave fields in the spirit of optical photonics. To this end, the strong nonlinearities of superconducting quantum circuits can be used to either improve or move beyond the implementation of concepts from the optical domain. In this context, the design of a well-controlled broadband environment for the superconducting quantum circuits is a central task. In this work, we place a superconducting transmon qubit in one arm of an on-chip Mach-Zehnder interferometer composed of two superconducting microwave beam splitters. By measuring its relaxation and dephasing rates we use the qubit as a sensitive spectrometer at the quantum level to probe the broadband electromagnetic environment. At high frequencies, this environment can be well described by an ensemble of harmonic oscillators coupled to the transmon qubit. At low frequencies, we find experimental evidence for colored quasi-static Gaussian noise with a high spectral weight, as it is typical for ensembles of two-level fluctuators. Our work paves the way towards possible applications of propagating microwave photons, such as emulating quantum impurity models or a novel architecture for quantum information processing.


title: Langevin Markov Chain Monte Carlo with stochastic gradients

abstract: Monte Carlo sampling techniques have broad applications in machine learning, Bayesian posterior inference, and parameter estimation. Often the target distribution takes the form of a product distribution over a dataset with a large number of entries. For sampling schemes utilizing gradient information it is cheaper for the derivative to be approximated using a random small subset of the data, introducing extra noise into the system. We present a new discretization scheme for underdamped Langevin dynamics when utilizing a stochastic (noisy) gradient. This scheme is shown to bias computed averages to second order in the stepsize while giving exact results in the special case of sampling a Gaussian distribution with a normally distributed stochastic gradient.


title: Fully discrete DPG methods for the Kirchhoff-Love plate bending model

abstract: We extend the analysis and discretization of the Kirchhoff-Love plate bending problem from [T. F\"uhrer, N. Heuer, A.H. Niemi, An ultraweak formulation of the Kirchhoff-Love plate bending model and DPG approximation, arXiv:1805.07835, 2018] in two aspects. First, we present a well-posed formulation and quasi-optimal DPG discretization that includes the gradient of the deflection. Second, we construct Fortin operators that prove the well-posedness and quasi-optimal convergence of lowest-order discrete schemes with approximated test functions for both formulations. Our results apply to the case of non-convex polygonal plates where shear forces can be less than $L_2$-regular. Numerical results illustrate expected convergence orders.


title: Speaker Recognition using Deep Belief Networks

abstract: Short time spectral features such as mel frequency cepstral coefficients(MFCCs) have been previously deployed in state of the art speaker recognition systems, however lesser heed has been paid to short term spectral features that can be learned by generative learning models from speech signals. Higher dimensional encoders such as deep belief networks (DBNs) could improve performance in speaker recognition tasks by better modelling the statistical structure of sound waves. In this paper, we use short term spectral features learnt from the DBN augmented with MFCC features to perform the task of speaker recognition. Using our features, we achieved a recognition accuracy of 0.95 as compared to 0.90 when using standalone MFCC features on the ELSDSR dataset.


title: Author Obfuscation Using Generalised Differential Privacy

abstract: The problem of obfuscating the authorship of a text document has received little attention in the literature to date. Current approaches are ad-hoc and rely on assumptions about an adversary's auxiliary knowledge which makes it difficult to reason about the privacy properties of these methods. Differential privacy is a well-known and robust privacy approach, but its reliance on the notion of adjacency between datasets has prevented its application to text document privacy. However, generalised differential privacy permits the application of differential privacy to arbitrary datasets endowed with a metric and has been demonstrated on problems involving the release of individual data points. In this paper we show how to apply generalised differential privacy to author obfuscation by utilising existing tools and methods from the stylometry and natural language processing literature.


title: Superconducting receiver arrays for magnetic resonance imaging

abstract: Superconducting QUantum-Interference Devices (SQUIDs) make magnetic resonance imaging (MRI) possible in ultra-low microtesla-range magnetic fields. In this work, we investigate the design parameters affecting the signal and noise performance of SQUID-based sensors and multichannel magnetometers for MRI of the brain. Besides sensor intrinsics, various noise sources along with the size, geometry and number of superconducting detector coils are important factors affecting the image quality. We derive figures of merit based on optimal combination of multichannel data, analyze different sensor array designs, and provide tools for understanding the signal detection and the different noise mechanisms. The work forms a guide to making design decisions for both imagingand sensor-oriented readers.

categories: physics.ins-det cond-mat.supr-con eess.IV eess.SP physics.app-ph


title: Logchain: Blockchain-assisted Log Storage

abstract: During the normal operation of a Cloud solution, no one usually pays attention to the logs except technical department, which may periodically check them to ensure that the performance of the platform conforms to the Service Level Agreements. However, the moment the status of a component changes from acceptable to unacceptable, or a customer complains about accessibility or performance of a platform, the importance of logs increases significantly. Depending on the scope of the issue, all departments, including management, customer support, and even the actual customer, may turn to logs to find out what has happened, how it has happened, and who is responsible for the issue. The party at fault may be motivated to tamper the logs to hide their fault. Given the number of logs that are generated by the Cloud solutions, there are many tampering possibilities. While tamper detection solution can be used to detect any changes in the logs, we argue that critical nature of logs calls for immutability. In this work, we propose a blockchain-based log system, called Logchain, that collects the logs from different providers and avoids log tampering by sealing the logs cryptographically and adding them to a hierarchical ledger, hence, providing an immutable platform for log storage.


title: On Abhyankar's lemma about ramification indices

abstract: We provide a simple proof of the fact that the ramification index of the compositum of two finite extensions of local fields is equal to the least common multiple of the ramification indices when at least one of the extensions is tamely ramified.


title: NuSTAR Detection of a Hard X-ray Source in the Supernova Remnant - Molecular Cloud Interaction Site of IC 443

abstract: We report on a broadband study of a complex X-ray source (1SAX J0618.0+2227) associated with the interaction site of the supernova remnant (SNR) IC 443 and ambient molecular cloud (MC) using NuSTAR, XMM_Newton, and Chandra observations. Its X-ray spectrum is composed of both thermal and non-thermal components. The thermal component can be equally well represented by either a thin plasma model with kT=0.19 keV or a blackbody model with kT=0.11 keV. The non-thermal component can be fit with either a power-law with Gamma~1.7 or a cutoff power-law with Gamma~1.5 and a cutoff energy at E_cut~18 keV. Using the newly obtained NuSTAR dataset, we test three possible scenarios for isolated X-ray sources in the SNR-MC interaction site: 1) pulsar wind nebula (PWN); 2) SNR ejecta fragment; 3) shocked molecular clump. We conclude that this source is most likely composed of a SNR ejecta (or a PWN) and surrounding shocked molecular clumps. The nature of this hard X-ray source in the SNR-MC interaction site of IC 443 may shed light on unidentified X-ray sources with hard X-ray spectra in rich environments for star forming regions, such as the Galactic center.


title: Single beam detection of optically driven spin dynamics in CdTe/(Cd,Mg)Te quantum wells

abstract: We study optical pumping of resident electron spins under resonant excitation of trions in n-type CdTe/(Cd,Mg)Te quantum wells subject to a transverse magnetic field. In contrast to the comprehensively used time-resolved pump-probe techniques with polarimetric detection, we exploit here a single beam configuration in which the time-integrated intensity of the excitation laser light transmitted through the quantum wells is detected. The transmitted intensity reflects the bleaching of light absorption due to optical pumping of the resident electron spins and can be used to evaluate the Larmor precession frequency of the optically oriented carriers and their spin relaxation time. Application of the magnetic field leads to depolarization of the electron spin ensemble so that the Hanle effect is observed. Excitation with a periodic sequence of laser pulses leads to optical pumping in the rotating frame if the Larmor precession frequency is synchronized with the pulse repetition rate. This is manifested by the appearance of Hanle curves every 3.36 or 44.2 mT for pulse repetition rates of 75.8 or 999 MHz, respectively. From the experimental data we evaluate the g factor of |g|=1.61 and the spin relaxation time of 14 ns for the optically pumped resident electrons, in agreement with previous time-resolved pump-probe studies.


title: A simple derivation of the two force laws for elliptic orbits from Proposition 6 in Newton's Principia

abstract: Based on Propostion 6 of his Principia, Newton's geometrical derivation in Propositions 10 and 11 for the radial dependence of the two central forces that lead to elliptical orbits is notoriously difficult. An alternate and more transparent derivation is obtained by applying the affine transformation of a circle into an ellipse.


title: Rigorous Analysis of a Randomised Number Field Sieve

abstract: Factorisation of integers $n$ is of number theoretic and cryptographic significance. The Number Field Sieve (NFS) introduced circa 1990, is still the state of the art algorithm, but no rigorous proof that it halts or generates relationships is known. We propose and analyse an explicitly randomised variant. For each $n$, we show that these randomised variants of the NFS and Coppersmith's multiple polynomial sieve find congruences of squares in expected times matching the best-known heuristic estimates.


title: Unsupervised Domain Adaptation using Regularized Hyper-graph Matching

abstract: Domain adaptation (DA) addresses the real-world image classification problem of discrepancy between training (source) and testing (target) data distributions. We propose an unsupervised DA method that considers the presence of only unlabelled data in the target domain. Our approach centers on finding matches between samples of the source and target domains. The matches are obtained by treating the source and target domains as hyper-graphs and carrying out a class-regularized hyper-graph matching using first-, second- and third-order similarities between the graphs. We have also developed a computationally efficient algorithm by initially selecting a subset of the samples to construct a graph and then developing a customized optimization routine for graph-matching based on Conditional Gradient and Alternating Direction Multiplier Method. This allows the proposed method to be used widely. We also performed a set of experiments on standard object recognition datasets to validate the effectiveness of our framework over state-of-the-art approaches.


title: Measurement of Transverse Single Spin Asymmetries in $\pi^0$ Production from $p^{\uparrow}+p$ and $p^{\uparrow}+A$ Collisions at STAR

abstract: In 2015 the first collisions between polarized protons and nuclei occurred at the Relativistic Heavy Ion Collider (RHIC), at a center-of-mass energy of $\sqrt{s_{NN}}=200$ GeV. Comparisons between spin asymmetries and cross-sections in $p+p$ production to those in $p+A$ production provide insight into nuclear structure, namely nuclear modification factors, nuclear dependence of spin asymmetries, and comparison to models with saturation effects. The transverse single-spin asymmetry, $A_{N}$, has been measured in $\pi^{0}$ production in the STAR Forward Meson Spectrometer (FMS), an electromagnetic calorimeter covering a forward psuedorapidity range of $2.6<\eta<4$. Within this kinematic range, STAR has previously reported the persistence of large $\pi^0$ asymmetries with unexpected dependences on $p_T$ and event topology in $p+p$ collisions. This talk will compare these dependences to those in $p+A$ production.


title: Soteria: Automated IoT Safety and Security Analysis

abstract: Broadly defined as the Internet of Things (IoT), the growth of commodity devices that integrate physical processes with digital systems have changed the way we live, play and work. Yet existing IoT platforms cannot evaluate whether an IoT app or environment is safe, secure, and operates correctly. In this paper, we present Soteria, a static analysis system for validating whether an IoT app or IoT environment (collection of apps working in concert) adheres to identified safety, security, and functional properties. Soteria operates in three phases; (a) translation of platform-specific IoT source code into an intermediate representation (IR), (b) extracting a state model from the IR, (c) applying model checking to verify desired properties. We evaluate Soteria on 65 SmartThings market apps through 35 properties and find nine (14%) individual apps violate ten (29%) properties. Further, our study of combined app environments uncovered eleven property violations not exhibited in the isolated apps. Lastly, we demonstrate Soteria on MalIoT, a novel open-source test suite containing 17 apps with 20 unique violations.


title: Adversarial Label Learning

abstract: We consider the task of training classifiers without labels. We propose a weakly supervised method---adversarial label learning---that trains classifiers to perform well against an adversary that chooses labels for training data. The weak supervision constrains what labels the adversary can choose. The method therefore minimizes an upper bound of the classifier's error rate using projected primal-dual subgradient descent. Minimizing this bound protects against bias and dependencies in the weak supervision. Experiments on three real datasets show that our method can train without labels and outperforms other approaches for weakly supervised learning.


title: ARiA: Utilizing Richard's Curve for Controlling the Non-monotonicity of the Activation Function in Deep Neural Nets

abstract: This work introduces a novel activation unit that can be efficiently employed in deep neural nets (DNNs) and performs significantly better than the traditional Rectified Linear Units (ReLU). The function developed is a two parameter version of the specialized Richard's Curve and we call it Adaptive Richard's Curve weighted Activation (ARiA). This function is non-monotonous, analogous to the newly introduced Swish, however allows a precise control over its non-monotonous convexity by varying the hyper-parameters. We first demonstrate the mathematical significance of the two parameter ARiA followed by its application to benchmark problems such as MNIST, CIFAR-10 and CIFAR-100, where we compare the performance with ReLU and Swish units. Our results illustrate a significantly superior performance on all these datasets, making ARiA a potential replacement for ReLU and other activations in DNNs.


title: Globule-like conformation and enhanced diffusion of active polymers

abstract: We study the dynamics and conformation of polymers composed by active monomers. By means of Brownian dynamics simulations we show that when the direction of the self-propulsion of each monomer is aligned with the backbone, the polymer undergoes a coil-to-globule-like transition, highlighted by a marked change of the scaling exponent of the gyration radius. Concurrently, the diffusion coefficient of the center of mass of the polymer becomes essentially independent of the polymer size for sufficiently long polymers or large magnitudes of the self-propulsion. These effects are reduced when the self-propulsion of the monomers is not bound to be tangent to the backbone of the polymer. Our results, rationalized by a minimal stochastic model, open new routes for activity-controlled polymer and, possibly, for a new generation of polymer-based drug carriers.


title: Benchmarking van der Waals-treated DFT: The case of hexagonal boron nitride and graphene on Ir(111)

abstract: There is enormous recent interest in weak, van der Waals-type (vdW) interactions due to their fundamental relevance for two-dimensional materials and the so-called vdW heterostructures. Tackling this problem using computer simulation is very challenging due to the non-trivial, non-local nature of these interactions. We benchmark different treatments of London dispersion forces within the density functional theory (DFT) framework on hexagonal boron nitride or graphene monolayers on Ir(111) by comparing the calculated geometries to a comprehensive set of experimental data. The geometry of these systems crucially depends on the interplay between vdW interactions and wave function hybridisation, making them excellent test cases for vdW-treated DFT. Our results show strong variations in the calculated atomic geometry. While some of the approximations reproduce the experimental structure, this is rather based on \textit{a posteriori} comparison with the ``target results''. General predictive power in vdW-treated DFT is not achieved yet and might require new approaches.


title: Double Helicity Asymmetries of Forward Neutral Pions from $\sqrt{s}=510$ GeV $pp$ Collisions at STAR

abstract: Longitudinally polarized $p+p$ scattering experiments provide access to gluon polarization via measurement of the double helicity asymmetry, $A_{LL}$. At the completion of the 2013 RHIC running period, a significant dataset of $\pi^{0}$s corresponding to an integrated luminosity of 46 $\text{pb}^{-1}$ (2012) and 8 $\text{pb}^{-1}$ (2013) produced from polarized $p+p$ scattering at $\sqrt{s}=510$ GeV with an average beam polarization of approximately $50\%$ was acquired. The $\pi^{0}$ kinematics were measured via isolation cones by the STAR Forward Meson Spectrometer, an electromagnetic calorimeter covering a forward pseudorapidity range of $2.6 < \eta <4$. The asymmetric $qg \to qg$ subprocess becomes more dominant in this forward region than in the midrapidity region; furthermore, asymmetry measurements in the forward region are sensitive to low-$x$ gluons. Progress on $A_{LL}$ determined from forward $\pi^{0}$ events, complementing previous midrapidity measurements, are presented.


title: Multi-task Maximum Entropy Inverse Reinforcement Learning

abstract: Multi-task Inverse Reinforcement Learning (IRL) is the problem of inferring multiple reward functions from expert demonstrations. Prior work, built on Bayesian IRL, is unable to scale to complex environments due to computational constraints. This paper contributes a formulation of multi-task IRL in the more computationally efficient Maximum Causal Entropy (MCE) IRL framework. Experiments show our approach can perform one-shot imitation learning in a gridworld environment that single-task IRL algorithms need hundreds of demonstrations to solve. We outline preliminary work using meta-learning to extend our method to the function approximator setting of modern MCE IRL algorithms. Evaluating on multi-task variants of common simulated robotics benchmarks, we discover serious limitations of these IRL algorithms, and conclude with suggestions for further work.


title: Sensitivity of Regular Estimators

abstract: This paper studies local asymptotic relationship between two scalar estimates. We define sensitivity of a target estimate to a control estimate to be the directional derivative of the target functional with respect to the gradient direction of the control functional. Sensitivity according to the information metric on the model manifold is the asymptotic covariance of regular efficient estimators. Sensitivity according to a general policy metric on the model manifold can be obtained from influence functions of regular efficient estimators. Policy sensitivity has a local counterfactual interpretation, where the ceteris paribus change to a counterfactual distribution is specified by the combination of a control parameter and a Riemannian metric on the model manifold.


title: The strong Frobenius numbers for cyclic defect blocks are equal to one

abstract: The purpose of this note is to provide a reference for the fact that the strong Frobenius number, in the sense of Eaton and Livesey, of a block of a finite group with a cyclic defect group is equal to one. This answers a question of Farrell and Kessar.


title: The primitive spectrum and category O for the periplectic Lie superalgebra

abstract: We solve two problems in representation theory for the periplectic Lie superalgebra pe(n), namely the description of the primitive spectrum in terms of functorial realisations of the braid group and the decomposition of category O into indecomposable blocks. To solve the first problem we establish a new type of equivalence between category O for all (not just simple or basic) classical Lie superalgebras and a category of Harish-Chandra bimodules. The latter bimodules have a left action of the Lie superalgebra but a right action of the underlying Lie algebra. To solve the second problem we establish a BGG reciprocity result for the periplectic Lie superalgebra.


title: Oberwolfach Report: Spacetime Intrinsic Flat Convergence

abstract: This Oberwolfach report describes published work with Carlos Vega introducing the null distance on a spacetime, and announces unpublished applications of this work jointly with Carlos Vega and with Anna Sakovich applying the null distance to define a spacetime intrinsic flat convergence (SF convergence). The announced work with Vega concerns convergence of sequences of big bang spacetimes and the announced work with Anna Sakovich concerns convergence of future maximal developments of initial data sets. Both were originally announced at the 2018 JMM. This report lists some open questions that will be discussed at the Oberwolfach meeting this coming August.


title: Geodesics on a Kerr-Newman-(anti-)de Sitter instanton

abstract: We study geodesics along a noncompact Kerr-Newman instanton, where the asymptotic geometry is either de Sitter or anti-de Sitter. We use first integrals for the Hamilton-Jacobi equation to characterize trajectories both near and away from horizons. We study the interaction of geodesics with special features of the metric, particularly regions of angular degeneracy or "theta horizons" in the de Sitter case. Finally, we characterize a number of stable equilibrium orbits.


title: KMT-2016-BLG-1107: A New Hollywood-Planet Close/Wide Degeneracy

abstract: We show that microlensing event KMT-2016-BLG-1107 displays a new type of degeneracy between wide-binary and close-binary Hollywood events in which a giant-star source envelops the planetary caustic. The planetary anomaly takes the form of a smooth, two-day "bump" far out on the falling wing of the light curve, which can be interpreted either as the source completely enveloping a minor-image caustic due to a close companion with mass ratio $q=0.036$, or partially enveloping a major-image caustic due to a wide companion with $q=0.004$. The best estimates of the companion masses are both in the planetary regime ($3.3^{+3.5}_{-1.8}\,M_{\rm jup}$ and $0.090^{+0.096}_{-0.037}\,M_{\rm jup}$) but differ by an even larger factor than the mass ratios due to different inferred host masses. We show that the two solutions can be distinguished by high-resolution imaging at first light on next-generation ("30m") telescopes. We provide analytic guidance to understand the conditions under which this new type of degeneracy can appear.


title: Spiking Linear Dynamical Systems on Neuromorphic Hardware for Low-Power Brain-Machine Interfaces

abstract: Neuromorphic architectures achieve low-power operation by using many simple spiking neurons in lieu of traditional hardware. Here, we develop methods for precise linear computations in spiking neural networks and use these methods to map the evolution of a linear dynamical system (LDS) onto an existing neuromorphic chip: IBM's TrueNorth. We analytically characterize, and numerically validate, the discrepancy between the spiking LDS state sequence and that of its non-spiking counterpart. These analytical results shed light on the multiway tradeoff between time, space, energy, and accuracy in neuromorphic computation. To demonstrate the utility of our work, we implemented a neuromorphic Kalman filter (KF) and used it for offline decoding of human vocal pitch from neural data. The neuromorphic KF could be used for low-power filtering in domains beyond neuroscience, such as navigation or robotics.


title: Step Size Matters in Deep Learning

abstract: Training a neural network with the gradient descent algorithm gives rise to a discrete-time nonlinear dynamical system. Consequently, behaviors that are typically observed in these systems emerge during training, such as convergence to an orbit but not to a fixed point or dependence of convergence on the initialization. Step size of the algorithm plays a critical role in these behaviors: it determines the subset of the local optima that the algorithm can converge to, and it specifies the magnitude of the oscillations if the algorithm converges to an orbit. To elucidate the effects of the step size on training of neural networks, we study the gradient descent algorithm as a discrete-time dynamical system, and by analyzing the Lyapunov stability of different solutions, we show the relationship between the step size of the algorithm and the solutions that can be obtained with this algorithm. The results provide an explanation for several phenomena observed in practice, including the deterioration in the training error with increased depth, the hardness of estimating linear mappings with large singular values, and the distinct performance of deep residual networks.


title: Losing the IR: a Holographic Framework for Area Theorems

abstract: Gravitational area laws are expected to arise as a result of ignorance of "UV gravitational data". In AdS/CFT, the UV/IR correspondence suggests that this data is dual to infrared physics in the CFT. Motivated by these heuristic expectations, we define a precise framework for explaining bulk area laws (in any dimension) by discarding IR CFT data. In (1+1) boundary dimensions, our prescribed mechanism shows explicitly that the boundary dual to these area laws is strong subadditivity of von Neumann entropy. Moreover, such area laws may be of arbitrary (and mixed) signature; thus our framework gives the first entropic explanation of mixed signature area laws (as well as area laws for certain dynamical causal horizons). In general dimension, the framework is easily modified to include bulk quantum corrections, thus giving rise to an infinite family of bulk generalized second laws.


title: Loss of Hall Conductivity Quantization in a Non-Hermitian Quantum Anomalous Hall Insulator

abstract: Recent work has extended topological band theory to open, non-Hermitian Hamiltonians, yet little is understood about how non-Hermiticity alters the topological quantization of associated observables. We address this problem by studying the quantum anomalous Hall effect (QAHE) generated in the Dirac surface states of a 3D time-reversal-invariant topological insulator (TI) that is proximity-coupled to a metallic ferromagnet. By constructing a contact self-energy for the ferromagnet, we show that in addition to generating a mass gap in the surface spectrum, the ferromagnet can introduce a non-Hermitian broadening term, which can obscure the mass gap in the spectral function. We calculate the Hall conductivity for the effective non-Hermitian Hamiltonian describing the heterostructure and show that it is no longer quantized despite being classified as a Chern insulator based on non-Hermitian topological band theory. Our results indicate that the QAHE will be challenging to experimentally observe in ferromagnet-TI heterostructures due to the finite lifetime of quasi-particles at the interface.


title: On-the-fly Vertex Reuse for Massively-Parallel Software Geometry Processing

abstract: Compute-mode rendering is becoming more and more attractive for non-standard rendering applications, due to the high flexibility of compute-mode execution. These newly designed pipelines often include streaming vertex and geometry processing stages. In typical triangle meshes, the same transformed vertex is on average required six times during rendering. To avoid redundant computation, a post-transform cache is traditionally suggested to enable reuse of vertex processing results. However, traditional caching neither scales well as the hardware becomes more parallel, nor can be efficiently implemented in a software design. We investigate alternative strategies to reusing vertex shading results on-the-fly for massively parallel software geometry processing. Forming static and dynamic batching on the data input stream, we analyze the effectiveness of identifying potential local reuse based on sorting, hashing, and efficient intra-thread-group communication. Altogether, we present four vertex reuse strategies, tailored to modern parallel architectures. Our simulations showcase that our batch-based strategies significantly outperform parallel caches in terms of reuse. On actual GPU hardware, our evaluation shows that our strategies not only lead to good reuse of processing results, but also boost performance by $2-3\times$ compared to na\"ively ignoring reuse in a variety of practical applications.


title: Are Computer Science and Engineering Graduates Ready for the Software Industry? Experiences from an Industrial Student Training Program

abstract: It has been 50 years since the term software engineering was coined in 1968 at a NATO conference. The field should be relatively mature by now, with most established universities covering core software engineering topics in their Computer Science programs and others offering specialized degrees. However, still many practitioners lament a lack of skills in new software engineering hires. With the growing demand for software engineers from the industry, this apparent gap becomes more and more pronounced. One corporate strategy to address this gap is for the industry to develop supplementary training programs before the hiring process, which could also help them screen viable candidates. In this paper, we report on our experiences and lessons learned in conducting a summer school program aimed at screening new graduates, introducing them to core skills relevant to the organization and the industry, and assessing their attitudes toward mastering those skills before the hiring process begins. Our experience suggests that such initiatives can be mutually beneficial for new hires and companies. We support this insight with pre- and post-training data collected from the participants during the first edition of such a summer school and a follow-up questionnaire conducted after a year with the graduates, 50% of whom was hired by the company shortly after the summer school.


title: Iterated local cohomology groups and Lyubeznik numbers for determinantal rings

abstract: We give an explicit recipe for determining iterated local cohomology groups with support in ideals of minors of a generic matrix in characteristic zero, expressing them as direct sums of indecomposable D-modules. For non-square matrices these indecomposables are simple, but this is no longer true for square matrices where the relevant indecomposables arise from the pole order filtration associated with the determinant hypersurface. Specializing our results to a single iteration, we determine the Lyubeznik numbers for all generic determinantal rings, thus answering a question of Hochster.


title: Rate-Maximizing OFDM Pilot Patterns for UAV Communications in Nonstationary A2G Channels

abstract: In this paper, we propose and evaluate rate-maximizing pilot configurations for Unmanned Aerial Vehicle (UAV) communications employing OFDM waveforms. OFDM relies on pilot symbols for effective communications. We formulate a rate-maximization problem in which the pilot spacing (in the time-frequency resource grid) and power is varied as a function of the time-varying channel statistics. The receiver solves this rate-maximization problem, and the optimal pilot spacing and power are explicitly fed back to the transmitter to adapt to the time-varying channel statistics in an air-to-ground (A2G) environment. We show the enhanced throughput performance of this scheme for UAV communications in sub-6 GHz bands. These performance gains are achieved at the cost of very low computational complexity and feedback requirements, making it attractive for A2G UAV communications in 5G.


title: Teacher's Perception in the Classroom

abstract: The ability for a teacher to engage all students in active learning processes in classroom constitutes a crucial prerequisite for enhancing students achievement. Teachers' attentional processes provide important insights into teachers' ability to focus their attention on relevant information in the complexity of classroom interaction and distribute their attention across students in order to recognize the relevant needs for learning. In this context, mobile eye tracking is an innovative approach within teaching effectiveness research to capture teachers' attentional processes while teaching. However, analyzing mobile eye-tracking data by hand is time consuming and still limited. In this paper, we introduce a new approach to enhance the impact of mobile eye tracking by connecting it with computer vision. In mobile eye tracking videos from an educational study using a standardized small group situation, we apply a state-ofthe-art face detector, create face tracklets, and introduce a novel method to cluster faces into the number of identity. Subsequently, teachers' attentional focus is calculated per student during a teaching unit by associating eye tracking fixations and face tracklets. To the best of our knowledge, this is the first work to combine computer vision and mobile eye tracking to model teachers' attention while instructing.


title: Energy Sustainable IoT with Individual QoS Constraints Through MISO SWIPT Multicasting

abstract: Enabling technologies for energy sustainable Internet of Things (IoT) are of paramount importance since the proliferation of high data communication demands of low power network devices. In this paper, we consider a Multiple Input Single Output (MISO) multicasting IoT system comprising of a multiantenna Transmitter (TX) simultaneously transferring information and power to low power and data hungry IoT Receivers (RXs). Each IoT device is assumed to be equipped with Power Splitting (PS) hardware that enables Energy Harvesting (EH) and imposes an individual Quality of Service (QoS) constraint to the downlink communication. We study the joint design of TX precoding and IoT PS ratios for the considered MISO Simultaneous Wireless Information and Power Transfer (SWIPT) multicasting IoT system with the objective of maximizing the minimum harvested energy among IoT, while satisfying their individual QoS requirements. In our novel EH fairness maximization formulation, we adopt a generic Radio Frequency (RF) EH model capturing practical rectification operation, and resulting in a nonconvex optimization problem. For this problem, we first present an equivalent semi-definite relaxation formulation and then prove it possesses unique global optimality. We also derive tight upper and lower bounds on the globally optimal solution that are exploited in obtaining low complexity algorithmic implementations for the targeted joint design. Analytical expressions for the optimal TX beamforming directions, power allocation, and IoT PS ratios are also presented. Our representative numerical results including comparisons with benchmark designs corroborate the usefulness of proposed framework and provide useful insights on the interplay of critical system parameters.


title: Echo: Compiler-based GPU Memory Footprint Reduction for LSTM RNN Training

abstract: The Long-Short-Term-Memory Recurrent Neural Networks (LSTM RNNs) are a popular class of machine learning models for analyzing sequential data. Their training on modern GPUs, however, is limited by the GPU memory capacity. Our profiling results of the LSTM RNN-based Neural Machine Translation (NMT) model reveal that feature maps of the attention and RNN layers form the memory bottleneck and runtime is unevenly distributed across different layers when training on GPUs. Based on these two observations, we propose to recompute the feature maps rather than stashing them persistently in the GPU memory. While the idea of feature map recomputation has been considered before, existing solutions fail to deliver satisfactory footprint reduction, as they do not address two key challenges. For each feature map recomputation to be effective and efficient, its effect on (1) the total memory footprint, and (2) the total execution time has to be carefully estimated. To this end, we propose *Echo*, a new compiler-based optimization scheme that addresses the first challenge with a practical mechanism that estimates the memory benefits of recomputation over the entire computation graph, and the second challenge by non-conservatively estimating the recomputation overhead leveraging layer specifics. *Echo* reduces the GPU memory footprint automatically and transparently without any changes required to the training source code, and is effective for models beyond LSTM RNNs. We evaluate *Echo* on numerous state-of-the-art machine learning workloads on real systems with modern GPUs and observe footprint reduction ratios of 1.89X on average and 3.13X maximum. Such reduction can be converted into faster training with a larger batch size, savings in GPU energy consumption (e.g., training with one GPU as fast as with four), and/or an increase in the maximum number of layers under the same GPU memory budget.


title: A new bound on Erd\H{o}s distinct distances problem in the plane over prime fields

abstract: In this paper we obtain a new lower bound on the Erd\H{o}s distinct distances problem in the plane over prime fields. More precisely, we show that for any set $A\subset \mathbb{F}_p^2$ with $|A|\le p^{7/6}$, the number of distinct distances determined by pairs of points in $A$ satisfies $$ |\Delta(A)| \gg |A|^{\frac{1}{2}+\frac{149}{4214}}.$$ Our result gives a new lower bound of $|\Delta{(A)}|$ in the range $|A|\le p^{1+\frac{149}{4065}}$. The main tools we employ are the energy of a set on a paraboloid due to Rudnev and Shkredov, a point-line incidence bound given by Stevens and de Zeeuw, and a lower bound on the number of distinct distances between a line and a set in $\mathbb{F}_p^2$. The latter is the new feature that allows us to improve the previous bound due Stevens and de Zeeuw.


title: Community Detection with Side Information: Exact Recovery under the Stochastic Block Model

abstract: The community detection problem involves making inferences about node labels in a graph, based on observing the graph edges. This paper studies the effect of additional, non-graphical side information on the phase transition of exact recovery in the binary stochastic block model (SBM) with $n$ nodes. When side information consists of noisy labels with error probability $\alpha$, it is shown that phase transition is improved if and only if $\log(\frac{1-\alpha}{\alpha})=\Omega(\log(n))$. When side information consists of revealing a fraction $1-\epsilon$ of the labels, it is shown that phase transition is improved if and only if $\log(1/\epsilon)=\Omega(\log(n))$. For a more general side information consisting of $K$ features, two scenarios are studied: (1)~$K$ is fixed while the likelihood of each feature with respect to corresponding node label evolves with $n$, and (2)~The number of features $K$ varies with $n$ but the likelihood of each feature is fixed. In each case, we find when side information improves the exact recovery phase transition and by how much. The calculated necessary and sufficient conditions for exact recovery are tight except for one special case. In the process of deriving inner bounds, a variation of an efficient algorithm is proposed for community detection with side information that uses a partial recovery algorithm combined with a local improvement procedure.


title: On Picard groups of blocks of finite groups

abstract: We show that the subgroup of the Picard group of a $p$-block of a finite group given by bimodules with endopermutation sources modulo the automorphism group of a source algebra is determined locally in terms of the fusion system on a defect group. We show that the Picard group of a block over the a complete discrete valuation ring ${\mathcal O}$ of characteristic zero with an algebraic closure $k$ of ${\mathbb F}_p$ as residue field is a colimit of finite Picard groups of blocks over $p$-adic subrings of ${\mathcal O}$. We apply the results to blocks with an abelian defect group and Frobenius inertial quotient, and specialise this further to blocks with cyclic or Klein four defect groups.


title: Compressive sensing lattice dynamics. II. Efficient phonon calculations and long-range interactions

abstract: We apply the compressive sensing lattice dynamics (CSLD) method to calculate phonon dispersion for crystalline solids. While existing methods such as frozen phonon, small displacement, and linear response are routinely applied for phonon calculations, they are considerable more expensive or cumbersome to apply to certain solids, including structures with large unit cells or low symmetry, systems that require more expensive electronic structure treatment, and polar semiconductors/insulators. In the latter case, we propose an approach based on a corrected long-range force constant model with proper treatment of the acoustic sum rule and the symmetric on-site force constant matrix. Our approach is demonstrated to be accurate and efficient for these systems through case studies of NaCl, CeO$_2$, Y$_3$Al$_5$O$_{12}$ and La$_2$Fe$_{14}$B.


title: Compressive sensing lattice dynamics. I. General formalism

abstract: {\it Ab initio\} calculations have been successfully used for evaluating lattice dynamical properties of solids within the (quasi-)harmonic approximation (i.e., assuming non-interacting phonons with infinite lifetimes), but it remains difficult to treat anharmonicity in all but the simplest compounds. We detail a systematic information theory based approach to deriving {\it ab initio\} anharmonic force constants: compressive sensing lattice dynamics (CSLD). The non-negligible terms that are necessary to reproduce the first-principles calculated interatomic forces are automatically selected by minimizing the $\ell_1$ norm (sum of absolute values) of the scaled force constants. By using efficient sampling of the configuration space using a modest number of atomic configurations with quasi-random displacements, CSLD is well suited for deriving accurate anharmonic potentials for complex multicomponent compounds with large unit cells. We demonstrate the power and generality of CSLD by calculating the phonon lifetimes and thermal transport properties of Type-I Si clathrates.


title: AffinityNet: semi-supervised few-shot learning for disease type prediction

abstract: While deep learning has achieved great success in computer vision and many other fields, currently it does not work very well on patient genomic data with the "big p, small N" problem (i.e., a relatively small number of samples with high-dimensional features). In order to make deep learning work with a small amount of training data, we have to design new models that facilitate few-shot learning. Here we present the Affinity Network Model (AffinityNet), a data efficient deep learning model that can learn from a limited number of training examples and generalize well. The backbone of the AffinityNet model consists of stacked k-Nearest-Neighbor (kNN) attention pooling layers. The kNN attention pooling layer is a generalization of the Graph Attention Model (GAM), and can be applied to not only graphs but also any set of objects regardless of whether a graph is given or not. As a new deep learning module, kNN attention pooling layers can be plugged into any neural network model just like convolutional layers. As a simple special case of kNN attention pooling layer, feature attention layer can directly select important features that are useful for classification tasks. Experiments on both synthetic data and cancer genomic data from TCGA projects show that our AffinityNet model has better generalization power than conventional neural network models with little training data. The code is freely available at https://github.com/BeautyOfWeb/AffinityNet .


title: Joint Optimal Design for Outage Minimization in DF Relay-assisted Underwater Acoustic Networks

abstract: This letter minimizes outage probability in a single decode-and-forward (DF) relay-assisted underwater acoustic network (UAN) without direct source-to-destination link availability. Specifically, a joint global-optimal design for relay positioning and allocating power to source and relay is proposed. For analytical insights, a novel low-complexity tight approximation method is also presented. Selected numerical results validate the analysis and quantify the comparative gains achieved using optimal power allocation (PA) and relay placement (RP) strategies.


title: Spatial analysis of airborne laser scanning point clouds for predicting forest variables

abstract: With recent developments in remote sensing technologies, plot-level forest resources can be predicted utilizing airborne laser scanning (ALS). The prediction is often assisted by mostly vertical summaries of the ALS point clouds. We present a spatial analysis of the point cloud by studying the horizontal distribution of the pulse returns through canopy height models thresholded at different height levels. The resulting patterns of patches of vegetation and gabs on each layer are summarized to spatial ALS features. We propose new features based on the Euler number, which is the number of patches minus the number of gaps, and the empty-space function, which is a spatial summary function of the gab space. The empty-space function is also used to describe differences in the gab structure between two different layers. We illustrate usefulness of the proposed spatial features for predicting different forest variables that summarize the spatial structure of forests or their breast height diameter distribution. We employ the proposed spatial features, in addition to commonly used features from literature, in the well-known k-nn estimation method to predict the forest variables. We present the methodology on the example of a study site in Central Finland.


title: Unitary evolution with fluctuations and dissipation

abstract: We outline an extension of the classical Langevin equation to a quantum formulation of the treatment of dissipation and fluctuations of all collective degrees of freedom with unitary evolution of a many-fermion system within an extension of the time-dependent density functional theory. We illustrate the method by computing the distribution of fission fragment yields for $^{258}$Fm in a quantum hydrodynamic approach and a typical trajectory with full unrestricted density functional theory augmented with dissipation and fluctuations.


title: Magneto-elastic induced vibronic bound state in the spin ice pyrochlore Ho$_2$Ti$_2$O$_7$

abstract: The single ion physics of Ho$_2$Ti$_2$O$_7$ is well-understood to produce strong Ising anisotropy, which is an essential ingredient to its low-temperature spin ice state. We present inelastic neutron scattering measurements on Ho$_2$Ti$_2$O$_7$ that reveal a clear inconsistency with its established single ion Hamiltonian. Specifically, we show that a crystal field doublet near 60~meV is split by approximately 3~meV. Furthermore, this crystal field splitting is not isolated to Ho$_2$Ti$_2$O$_7$ but can also be found in its chemical pressure analogs, Ho$_2$Ge$_2$O$_7$ and Ho$_2$Sn$_2$O$_7$. We demonstrate that the origin of this effect is a vibronic bound state, resulting from the entanglement of a phonon and crystal field excitation. We derive the microscopic Hamiltonian that describes the magneto-elastic coupling and provides a quantitative description of the inelastic neutron spectra.