title: Search for low-mass Dark Matter with the CRESST Experiment

abstract: CRESST is a multi-stage experiment directly searching for dark matter (DM) using cryogenic $\mathrm{CaWO_4}$ crystals. Previous stages established leading limits for the spin-independent DM-nucleon cross section down to DM-particle masses $m_\mathrm{DM}$ below $1\,\mathrm{GeV/c^2}$. Furthermore, CRESST performed a dedicated search for dark photons (DP) which excludes new parameter space between DP masses $m_\mathrm{DP}$ of $300\,\mathrm{eV/c^2}$ to $700\,\mathrm{eV/c^2}$. In this contribution we will discuss the latest results based on the previous CRESST-II phase 2 and we will report on the status of the current CRESST-III phase 1: in this stage we have been operating 10 upgraded detectors with $24,\mathrm{g}$ target mass each and enhanced detector performance since summer 2016. The improved detector design in terms of background suppression and reduction of the detection threshold will be discussed with respect to the previous stage. We will conclude with an outlook on the potential of the next stage, CRESST-III phase 2.


title: Dynamical Imaging with Interferometry

abstract: By linking widely separated radio dishes, the technique of very long baseline interferometry (VLBI) can greatly enhance angular resolution in radio astronomy. However, at any given moment, a VLBI array only sparsely samples the information necessary to form an image. Conventional imaging techniques partially overcome this limitation by making the assumption that the observed cosmic source structure does not evolve over the duration of an observation, which enables VLBI networks to accumulate information as the Earth rotates and changes the projected array geometry. Although this assumption is appropriate for nearly all VLBI, it is almost certainly violated for submillimeter observations of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), which has a gravitational timescale of only ~20 seconds and exhibits intra-hour variability. To address this challenge, we develop several techniques to reconstruct dynamical images ("movies") from interferometric data. Our techniques are applicable to both single-epoch and multi-epoch variability studies, and they are suitable for exploring many different physical processes including flaring regions, stable images with small time-dependent perturbations, steady accretion dynamics, or kinematics of relativistic jets. Moreover, dynamical imaging can be used to estimate time-averaged images from time-variable data, eliminating many spurious image artifacts that arise when using standard imaging methods. We demonstrate the effectiveness of our techniques using synthetic observations of simulated black hole systems and 7mm Very Long Baseline Array observations of M87, and we show that dynamical imaging is feasible for Event Horizon Telescope observations of Sgr A*.


title: Discovering More Precise Process Models from Event Logs by Filtering Out Chaotic Activities

abstract: Process Discovery is concerned with the automatic generation of a process model that describes a business process from execution data of that business process. Real life event logs can contain chaotic activities. These activities are independent of the state of the process and can, therefore, happen at rather arbitrary points in time. We show that the presence of such chaotic activities in an event log heavily impacts the quality of the process models that can be discovered with process discovery techniques. The current modus operandi for filtering activities from event logs is to simply filter out infrequent activities. We show that frequency-based filtering of activities does not solve the problems that are caused by chaotic activities. Moreover, we propose a novel technique to filter out chaotic activities from event logs. We evaluate this technique on a collection of seventeen real-life event logs that originate from both the business process management domain and the smart home environment domain. As demonstrated, the developed activity filtering methods enable the discovery of process models that are more behaviorally specific compared to process models that are discovered using standard frequency-based filtering.


title: Finding purifications with minimal entanglement

abstract: Purification is a tool that allows to represent mixed quantum states as pure states on enlarged Hilbert spaces. A purification of a given state is not unique and its entanglement strongly depends on the particular choice made. Moreover, in one-dimensional systems, the amount of entanglement is linked to how efficiently the purified state can be represented using matrix-product states (MPS). We introduce an MPS based method that allows to find the minimally entangled representation by iteratively minimizing the second Renyi entropy. First, we consider the thermofield double purification and show that its entanglement can be strongly reduced especially at low temperatures. Second, we show that a slowdown of the entanglement growth following a quench of an infinite temperature state is possible.


title: Geometric theory of the natural optical activity in noncentrosymmetric metals

abstract: This is a PhD dissertation. Ignited by the chiral anomaly of recently discovered Weyl (semi-)metals, we study the chiral magnetic effect and the natural optical activity of noncentrosymmetric metals. Both phenomena are related to the linear-in-$\mathbf{q}$ spatial dispersion of the optical conductivity tensor, and can be calculated within the formalism of the semiclassical kinetic equation. Therefore, we calculate the dispersion of optical conductivity up to the linear order of the wave vector, in the low frequency regime, with both the semiclassical Boltzmann equation and the Kubo formula. The two different methods of calculation provide us the same result. In this result, the static and dynamic chiral magnetic effects are revealed to have different origin: one comes from topology, related to Berry monopoles, and the other has a geometric origin, which is determined by the orbital magnetic moment. The Faraday rotation of the polarization of light transmitted through a slab of the sample provides us the most direct way to measure the magnitude of dynamic chiral magnetic effect. We develop an effective medium theory for electromagnetic wave propagating through gapless nonuniform systems, to calculate macroscopic sample inhomogeneities induced corrections to the chiral magnetic conductivity. We show that, in metals with low carrier density, the way in which macroscopic fluctuations of the local conductivity affect the frequency dependent of the measured optical polarization rotation angle: by creating a sharp feature near the plasma edge, which can be detected by experiments. Then we narrow down our research a bit further to the current induced magnetization.


title: Transport of internetwork magnetic flux elements in the solar photosphere

abstract: The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfven wave excitation spectrum for coronal heating models. We examine the motions of internetwork flux elements in a 24 hour long Hinode/NFI magnetogram sequence with 90 second cadence, and study both the scaling of their mean squared displacement and the shape of their displacement probability distribution as a function of time. We find that the mean squared displacement scales super-diffusively with a slope of about 1.48. Super-diffusive scaling has been observed in other studies for temporal increments as small as 5 seconds, increments over which ballistic scaling would be expected. Using high-cadence MURaM simulations, we show that the observed super-diffusive scaling at short temporal increments is a consequence of random changes in the barycenter positions due to flux evolution. We also find that for long temporal increments, beyond granular lifetimes, the observed displacement distribution deviates from that expected for a diffusive process, evolving from Rayleigh to Gaussian. This change in the distribution can be modeled analytically by accounting for supergranular advection along with motions due to granulation. These results complicate the interpretation of magnetic element motions as strictly advective or diffusive on short and long timescales and suggest that measurements of magnetic element motions must be used with caution in turbulent diffusion or wave excitation models. We propose that passive trace motions in measured photospheric flows may yield more robust transport statistics.


title: Degenerate and chiral states in the extended Heisenberg model in the kagome lattice

abstract: We present a study of the low temperature phases of the antiferromagnetic extended classical Heisenberg model in the kagome lattice, up to third nearest neighbors. First, we focus on the degenerate lines in the boundaries of the well-known staggered chiral phases. These boundaries have either semi-extensive or extensive degeneracy, and we discuss the partial selection of states by thermal fluctuations. Then, we study the model under an external magnetic field on these lines and in the staggered chiral phases. We pay particular attention to the highly frustrated point, where the three exchange couplings are equal. We show that this point can me mapped to a model with spin liquid behavior and non-zero chirality. Finally, we explore the effect of Dzyaloshinskii-Moriya (DM) interactions in two ways: an homogeneous and a staggered DM interaction. In both cases, there is a rich low temperature phase diagram, with different spontaneously broken symmetries and non trivial chiral phases.


title: General Physical Properties of CGRaBS Blazars

abstract: We present the results of a multi-frequency, time-averaged analysis of blazars included in the Candidate Gamma-ray Blazar Survey catalog. Our sample consists of 324 $\gamma$-ray detected ($\gamma$-ray loud) and 191 non $\gamma$-ray detected ($\gamma$-ray quiet) blazars, and we consider all the data up to 2016 April 1. We find that both the $\gamma$-ray loud and the $\gamma$-ray quiet blazar populations occupy similar regions in the WISE color-color diagram, and in the radio and X-ray bands $\gamma$-ray loud sources are brighter. A simple one-zone synchrotron inverse-Compton emission model is applied to derive the physical properties of both populations. We find that the central black hole mass and the accretion disk luminosity ($L_{\rm disk}$) computed from the modeling of the optical-UV emission with a Shakura-Sunyaev disk reasonably matches with that estimated from the optical spectroscopic emission-line information. A significantly larger Doppler boosting in the $\gamma$-ray loud blazars is noted, and their jets are more radiatively efficient. On the other hand, the $\gamma$-ray quiet objects are more MeV-peaked, thus could be potential targets for next-generation MeV missions. Our results confirm the earlier findings about the accretion-jet connection in blazars; however, many of the $\gamma$-ray quiet blazars tend to deviate from the recent claim that the jet power exceeds $L_{\rm disk}$ in blazars. A broadband study, considering a larger set of $\gamma$-ray quiet objects and also including BL Lacs, will be needed to confirm/reject this hypothesis and also to verify the evolution of the powerful high-redshift blazars into their low-power nearby counterparts.


title: Localization of Multiple Targets with Identical Radar Signatures in Multipath Environments with Correlated Blocking

abstract: This paper addresses the problem of localizing an unknown number of targets, all having the same radar signature, by a distributed MIMO radar consisting of single antenna transmitters and receivers that cannot determine directions of departure and arrival. Furthermore, we consider the presence of multipath propagation, and the possible (correlated) blocking of the direct paths (going from the transmitter and reflecting off a target to the receiver). In its most general form, this problem can be cast as a Bayesian estimation problem where every multipath component is accounted for. However, when the environment map is unknown, this problem is ill-posed and hence, a tractable approximation is derived where only direct paths are accounted for. In particular, we take into account the correlated blocking by scatterers in the environment which appears as a prior term in the Bayesian estimation framework. A sub-optimal polynomial-time algorithm to solve the Bayesian multi-target localization problem with correlated blocking is proposed and its performance is evaluated using simulations. We found that when correlated blocking is severe, assuming the blocking events to be independent and having constant probability (as was done in previous papers) resulted in poor detection performance, with false alarms more likely to occur than detections.


title: Inverses of Cartan Matrices of Lie Algebras and Lie Superalgebra

abstract: We derive explicit formulas for the inverses of the Cartan matrices of the simple Lie algebras and the basic classical Lie superalgebras, as well as for their infinite generalizations.


title: Optimal Data Acquisition for Statistical Estimation

abstract: We consider a data analyst's problem of purchasing data from strategic agents to compute an unbiased estimate of a statistic of interest. Agents incur private costs to reveal their data and the costs can be arbitrarily correlated with their data. Once revealed, data are verifiable. This paper focuses on linear unbiased estimators. We design an individually rational and incentive compatible mechanism that optimizes the worst-case mean-squared error of the estimation, where the worst-case is over the unknown correlation between costs and data, subject to a budget constraint in expectation. We characterize the form of the optimal mechanism in closed-form. We further extend our results to acquiring data for estimating a parameter in regression analysis, where private costs can correlate with the values of the dependent variable but not with the values of the independent variables.


title: The classical double copy in maximally symmetric spacetimes

abstract: The classical double copy procedure relates classical asymptotically-flat gravitational field solutions to Yang-Mills and scalar field solutions living in Minkowski space. In this paper we extend this correspondence to maximally symmetric curved spacetimes. We consider asymptotically (A)dS spacetimes in Kerr-Schild form and construct the corresponding single and zeroth copies. In order to clarify the interpretation of these copies, we study several examples including (A)dS-Schwarzschild, (A)dS-Kerr, black strings, black branes, and waves, paying particular attention to the source terms. We find that the single and zeroth copies of stationary solutions satisfy different equations than those of wave solutions. We also consider how to obtain Einstein-Maxwell solutions using this procedure. Finally, we derive the classical single and zeroth copy of the BTZ black hole.


title: Implicit Weight Uncertainty in Neural Networks

abstract: Modern neural networks tend to be overconfident on unseen, noisy or incorrectly labelled data and do not produce meaningful uncertainty measures. Bayesian deep learning aims to address this shortcoming with variational approximations (such as Bayes by Backprop or Multiplicative Normalising Flows). However, current approaches have limitations regarding flexibility and scalability. We introduce Bayes by Hypernet (BbH), a new method of variational approximation that interprets hypernetworks as implicit distributions. It naturally uses neural networks to model arbitrarily complex distributions and scales to modern deep learning architectures. In our experiments, we demonstrate that our method achieves competitive accuracies and predictive uncertainties on MNIST and a CIFAR5 task, while being the most robust against adversarial attacks.


title: About ionization zones remaining after supernovae explosions

abstract: It is shown that investigation of ionization zones left after supernova explosions can yield information on the mechanism active in those explosions.


title: BoostClean: Automated Error Detection and Repair for Machine Learning

abstract: Predictive models based on machine learning can be highly sensitive to data error. Training data are often combined with a variety of different sources, each susceptible to different types of inconsistencies, and new data streams during prediction time, the model may encounter previously unseen inconsistencies. An important class of such inconsistencies is domain value violations that occur when an attribute value is outside of an allowed domain. We explore automatically detecting and repairing such violations by leveraging the often available clean test labels to determine whether a given detection and repair combination will improve model accuracy. We present BoostClean which automatically selects an ensemble of error detection and repair combinations using statistical boosting. BoostClean selects this ensemble from an extensible library that is pre-populated general detection functions, including a novel detector based on the Word2Vec deep learning model, which detects errors across a diverse set of domains. Our evaluation on a collection of 12 datasets from Kaggle, the UCI repository, real-world data analyses, and production datasets that show that Boost- Clean can increase absolute prediction accuracy by up to 9% over the best non-ensembled alternatives. Our optimizations including parallelism, materialization, and indexing techniques show a 22.2x end-to-end speedup on a 16-core machine.


title: Origin of Interstellar Object A/2017 U1 in a Nearby Young Stellar Association?

abstract: We suggest that the kinematics and properties of interstellar object A/2017 U1 point towards its formation in a protoplanetary disk in the ~45 Myr-old Carina or Columba young stellar associations, and subsequent ejection with a low peculiar velocity (1-2 km/sec) during or soon after planet formation inside the ice line. This scenario predicts that the Solar System will encounter more such objects with radiants similar to that of A/2017 U1.


title: Accurate Solutions of Polynomial Eigenvalue Problems

abstract: Quadratic eigenvalue problems (QEP) and more generally polynomial eigenvalue problems (PEP) are among the most common types of nonlinear eigenvalue problems. Both problems, especially the QEP, have extensive applications. A typical approach to solve QEP and PEP is to use a linearization method to reformulate the problem as a higher dimensional linear eigenvalue problem. In this article, we use homotopy continuation to solve these nonlinear eigenvalue problems without passing to higher dimensions. Our main contribution is to show that our method produces substantially more accurate results, and finds all eigenvalues with a certificate of correctness via Smale's $\alpha$-theory. To explain the superior accuracy, we show that the nonlinear eigenvalue problem we solve is better conditioned than its reformulated linear eigenvalue problem, and our homotopy continuation algorithm is more stable than QZ algorithm - theoretical findings that are borne out by our numerical experiments. Our studies provide yet another illustration of the dictum in numerical analysis that, for reasons of conditioning and stability, it is sometimes better to solve a nonlinear problem directly even when it could be transformed into a linear problem with the same solution mathematically.


title: A bifurcation and symmetry discussion of the Sommerfeld effect

abstract: The Arnold Sommerfeld effect is an intriguing resonance capture and release series of events originally demonstrated in 1902. A single event is studied using a two degree of freedom mathematical model of a motor with imbalance mounted to laterally restricted spring connected cart. For a certain power supplied, in general the motor rotates at a speed consistent with a motor on a rigid base. However at speeds close to the natural frequency of the cart, it seemingly takes on extra oscillations where for a single rotation it both speeds up and then slows down. Therefore in a standard experimental demonstration of the effect, as the supplied torque force is increased or decreased, this may give the illusion that the stable operation of the motor is losing and gaining stability. This is not strictly the case, instead small oscillations always present in the system solution are amplified near the resonant frequency. The imbalance in the motor causes a single resonance curve to fold back on itself forming two fold bifurcations which leads to hysteresis and an asymmetry between increasing and decreasing the motor speed. Although as outlined the basic mechanism is due the interplay between two stable and one unstable limit cycles, a more complicated bifurcation scenario is observed for higher imbalances in the motor. The presence of a Z2 phase space symmetry tempers the dynamics and bifurcation picture.


title: On global minimizers of quadratic functions with cubic regularization

abstract: In this paper, we analyze some theoretical properties of the problem of minimizing a quadratic function with a cubic regularization term, arising in many methods for unconstrained and constrained optimization that have been proposed in the last years. First we show that, given any stationary point that is not a global solution, it is possible to compute, in closed form, a new point with a smaller objective function value. Then, we prove that a global minimizer can be obtained by computing a finite number of stationary points. Finally, we extend these results to the case where stationary conditions are approximately satisfied, discussing some possible algorithmic applications.


title: Using ultrafilters to prove Ramsey-type theorems

abstract: Ultrafilters are a tool, originating in mathematical logic and general topology, that has steadily found more and more uses in multiple areas of mathematics, such as combinatorics, dynamics, and algebra, among others. The purpose of this article is to introduce ultrafilters in a friendly manner and present some applications to the branch of combinatorics known as Ramsey theory, culminating with a new ultrafilter-based proof of van der Waerden's theorem.


title: Operation of a superconducting nanowire in two detection modes: KID and SPD

abstract: We present the performance of a superconducting nanowire that can be operated in two detection modes: i) as a kinetic inductance detector (KID) or ii) as a single-photon detector (SPD). Two superconducting nanowires developed for use as single-photon detectors (SNSPDs) are embedded as the inductive (L) component in resonant inductor/capacitor (LC) circuits coupled to a microwave transmission line. The capacitors are low loss commercial chip capacitors and limit the internal quality factor of the resonators to approximately $Q_i = 170$. The resonator quality factor, $Q_r \simeq 23$, is dominated by the coupling to the feedline and limits the detection bandwidth to on the order of 1MHz. When operated in KID mode, the detectors are AC biased with tones at their resonant frequencies of 45.85 and 91.81MHz. In the low-bias, standard KID mode, a single photon produces a hot spot that does not turn an entire section of the line normal but only increases the kinetic inductance. In the high-bias, critical KID mode, a photon event turns a section of the line normal and the resonance is destroyed until the normal region is dissipated. When operated as an SPD in Geiger mode, the resonators are DC biased through cryogenic bias tees and each photon produces a sharp voltage step followed by a ringdown signal at the resonant frequency of the detector which is converted to a standard pulse with an envelop detector. We show that AC biasing in the critical KID mode is inferior to the sensitivity achieved in DC-biased SPD mode due to the small fraction of time spent near the critical current with an AC bias.


title: Deep Learning-Based Dynamic Watermarking for Secure Signal Authentication in the Internet of Things

abstract: Securing the Internet of Things (IoT) is a necessary milestone toward expediting the deployment of its applications and services. In particular, the functionality of the IoT devices is extremely dependent on the reliability of their message transmission. Cyber attacks such as data injection, eavesdropping, and man-in-the-middle threats can lead to security challenges. Securing IoT devices against such attacks requires accounting for their stringent computational power and need for low-latency operations. In this paper, a novel deep learning method is proposed for dynamic watermarking of IoT signals to detect cyber attacks. The proposed learning framework, based on a long short-term memory (LSTM) structure, enables the IoT devices to extract a set of stochastic features from their generated signal and dynamically watermark these features into the signal. This method enables the IoT's cloud center, which collects signals from the IoT devices, to effectively authenticate the reliability of the signals. Furthermore, the proposed method prevents complicated attack scenarios such as eavesdropping in which the cyber attacker collects the data from the IoT devices and aims to break the watermarking algorithm. Simulation results show that, with an attack detection delay of under 1 second the messages can be transmitted from IoT devices with an almost 100% reliability.


title: Design of two-dimensional particle assemblies using isotropic pair interactions with an attractive well

abstract: Using ground-state and relative-entropy based inverse design strategies, isotropic interactions with an attractive well are determined to stabilize and promote as- sembly of particles into two-dimensional square, honeycomb, and kagome lattices. The design rules inferred from these results are discussed and validated in the dis- covery of interactions that favor assembly of the highly open truncated-square and truncated-hexagonal lattices.


title: On graphs of bounded semilattices

abstract: In this paper, we introduce the graph $G(S)$ of a bounded semilattice $S$, which is a generalization of the intersection graph of the substructures of an algebraic structure. We prove some general theorems about these graphs; as an example, we show that if $S$ is a product of three or more chains, then $G(S)$ is Eulerian if and only if either the length of every chain is even or all the chains are of length one. We also show that if $G(S)$ contains a cycle, then $girth(G(S)) = 3$. Finally, we show that if $(S,+,\cdot,0,1)$ is a dually atomic bounded distributive lattice whose set of dual atoms is nonempty, and the graph $G(S)$ of $S$ has no isolated vertex, then $G(S)$ is connected with $diam(G(S))\leq 4$.


title: VIP2 at Gran Sasso - Test of the validity of the spin statistics theorem for electrons with X-ray spectroscopy

abstract: In the VIP2 (VIolation of the Pauli Exlusion Principle) experiment at the Gran Sasso underground laboratory (LNGS) we are searching for possible violations of standard quantum mechanics predictions. With high precision we investigate the Pauli Exclusion Principle and the collapse of the wave function (collapse models). We will present our experimental method of searching for possible small violations of the Pauli Exclusion Principle for electrons, via the search for "anomalous" X-ray transitions in copper atoms, produced by "new" electrons (brought inside a copper bar by circulating current) which could have the probability to undergo Pauli-forbidden transition to the ground state (1 s level) already occupied by two electrons. We will describe the concept of the VIP2 experiment taking data at LNGS presently. The goal of VIP2 is to test the PEP for electrons with unprecedented accuracy, down to a limit in the probability that PEP is violated at the level of 10$^{-31}$. We will show preliminary experimental results obtained at LNGS and discuss implications of a possible violation.


title: Assessing numerical methods for molecular and particle simulation

abstract: We discuss the design of state-of-the-art numerical methods for molecular dynamics, focusing on the demands of soft matter simulation, where the purposes include sampling and dynamics calculations both in and out of equilibrium. We discuss the characteristics of different algorithms, including their essential conservation properties, the convergence of averages, and the accuracy of numerical discretizations. Formulations of the equations of motion which are suited to both equilibrium and nonequilibrium simulation include Langevin dynamics, dissipative particle dynamics (DPD), and the more recently proposed "pairwise adaptive Langevin" (PAdL) method, which, like DPD but unlike Langevin dynamics, conserves momentum and better matches the relaxation rate of orientational degrees of freedom. PAdL is easy to code and suitable for a variety of problems in nonequilibrium soft matter modeling, our simulations of polymer melts indicate that this method can also provide dramatic improvements in computational efficiency. Moreover we show that PAdL gives excellent control of the relaxation rate to equilibrium. In the nonequilibrium setting, we further demonstrate that while PAdL allows the recovery of accurate shear viscosities at higher shear rates than are possible using the DPD method at identical timestep, it also outperforms Langevin dynamics in terms of stability and accuracy at higher shear rates.


title: Open clusters: III. Fundamental parameters of B stars in NGC 6087, NGC 6250, NGC 6383 and NGC 6530. B type stars with circumstellar envelopes

abstract: Context. Stellar physical properties of star clusters are poorly known. Aims. Our goals are to perform a spectrophotometric study of the B star population in open clusters to derive accurate stellar parameters, search for the presence of circumstellar envelopes, and discuss the characteristics of these stars. Methods. The BCD spectrophotometric system is a powerful tool to obtain fundamental parameters and infer the main properties of open clusters: distance modulus, color excess, and age. We inspected the Balmer discontinuity to seek circumstellar disks and identify Be-star candidates. High resolution spectra in the H$\alpha$ region are used to confirm the Be nature. Results. We provide Teff , log g, Mv , Mbol and spectral types for a sample of 68 stars in the field of the open clusters NGC 6087, NGC 6250, NGC 6383, and NGC 6530, as well as the cluster distances, ages and reddening. Then, based on a sample of 230 B stars in the direction of the 11 open clusters studied along this series of three papers, we report 6 new Be stars, 4 blue straggler candidates, and 15 B-type stars (called Bdd) with a double Balmer discontinuity. We also find that the majority of the Be stars are dwarfs and present a maximum at the spectral type B2-B4 in young and intermediate-age open clusters. Another maximum of Be stars is observed at the spectral type B6-B8 in open clusters older than 40 Myr, where the population of Bdd stars also becomes relevant. Conclusions. Our results support previous statements that the Be phenomenon is present along the whole main sequence band and occurs in very different evolutionary states. We find clear evidence of an increase of stars with circumstellar envelopes with cluster age. The Be phenomenon reaches its maximum in clusters of intermediate age and the number of B stars with circumstellar envelopes (Be+Bdd stars) is also high for the older clusters.


title: NeuralFDR: Learning Discovery Thresholds from Hypothesis Features

abstract: As datasets grow richer, an important challenge is to leverage the full features in the data to maximize the number of useful discoveries while controlling for false positives. We address this problem in the context of multiple hypotheses testing, where for each hypothesis, we observe a p-value along with a set of features specific to that hypothesis. For example, in genetic association studies, each hypothesis tests the correlation between a variant and the trait. We have a rich set of features for each variant (e.g. its location, conservation, epigenetics etc.) which could inform how likely the variant is to have a true association. However popular testing approaches, such as Benjamini-Hochberg's procedure (BH) and independent hypothesis weighting (IHW), either ignore these features or assume that the features are categorical or uni-variate. We propose a new algorithm, NeuralFDR, which automatically learns a discovery threshold as a function of all the hypothesis features. We parametrize the discovery threshold as a neural network, which enables flexible handling of multi-dimensional discrete and continuous features as well as efficient end-to-end optimization. We prove that NeuralFDR has strong false discovery rate (FDR) guarantees, and show that it makes substantially more discoveries in synthetic and real datasets. Moreover, we demonstrate that the learned discovery threshold is directly interpretable.


title: On causality violation on a Kerr-de Sitter spacetime

abstract: The causal properties of the family of Kerr-de Sitter spacetimes are analyzed and compared to those of the Kerr family. At first we show that a Kerr-de Sitter spacetime can be viewed as an assembly of Carter's blocks i.e. four dimensional spacetime regions contained within Killing horizons or a Killing horizon and the asymptotic de Sitter region. From this perspective and leaving aside topological identifications, the causal properties of a Kerr de Sitter spacetime are determined by the causal properties of the individual Carter's blocks viewed as spacetimes in their own right. We show that any Carter's block is stably causal except for the blocks that contain the ring singularity. The latter are vicious sets: any two events within such block can be connected by a future (respectively past) directed timelike curve. This behavior is identical to the causal behavior of the Boyer-Lindquist blocks that contains the Kerr ring singularity which are also vicious sets. On the other hand, while for the case of a naked Kerr singularity the entire spacetime is a vicious set and thus closed timelike curves pass through any event including events lying in the asymptotic region, for the case of a Kerr-de Sitter spacetime the cosmological horizons imply that the asymptotic de Sitter region is causally well behaved. In that regard a positive cosmological constant appears to improve the causal behavior of the underlying spacetime.


title: Intensity-Coupled-Polarization in Instruments with a Continuously Rotating Half-Wave Plate

abstract: We discuss a systematic effect associated with measuring polarization with a continuously rotating half-wave plate. The effect was identified with the data from the E and B Experiment (EBEX), which was a balloon-borne instrument designed to measure the polarization of the CMB as well as that from Galactic dust. The data show polarization fraction larger than 10\% while less than 3\% were expected from instrumental polarization. We give evidence that the excess polarization is due to detector non-linearity in the presence of a continuously rotating HWP. The non-linearity couples intensity signals into polarization. We develop a map-based method to remove the excess polarization. Applying this method for the 150 (250) GHz bands data we find that 81\% (92\%) of the excess polarization was removed. Characterization and mitigation of this effect is important for future experiments aiming to measure the CMB B-modes with a continuously rotating HWP.


title: Theory of magnetotransport in artificial kagome spin ice

abstract: Magnetic nanoarray with special geometries exhibits nontrivial collective behaviors similar to those observed in the spin ice materials. Here we present a novel circuit model to describe the complex magnetotransport phenomena in artificial kagome spin ice. In this picture, the system can be viewed as a resistor network driven by voltage sources that are located at vertices of the honeycomb array. The differential voltages across different terminals of these sources are related to the ice-rules that govern the local magnetization ordering. The circuit model relates the transverse Hall voltage of kagome ice to the underlying spin correlations. Treating the magnetic nanoarray as metamaterials, we present a mesoscopic constitutive equation relating the Hall resistance to magnetization components of the system. We further show that the Hall signal is significantly enhanced when the kagome ice undergoes a magnetic-charge ordering transition. Our analysis can be readily generalized to other lattice geometry, providing a quantitative method for the design of magnetoresistance devices based on artificial spin ices.


title: Learning Stable and Energetically Economical Walking with RAMone

abstract: In this paper, we optimize over the control parameter space of our planar-bipedal robot, RAMone, for stable and energetically economical walking at various speeds. We formulate this task as an episodic reinforcement learning problem and use Covariance Matrix Adaptation. The parameters we are interested in modifying include gains from our Hybrid Zero Dynamics style controller and from RAMone's low-level motor controllers.


title: Small-scale equidistribution for random spherical harmonics

abstract: We study random spherical harmonics at shrinking scales. We compare the mass assigned to a small spherical cap with its area, and find the smallest possible scale at which, with high probability, the discrepancy between them is small simultaneously at every point on the sphere.


title: Improving Exoplanet Detection Power: Multivariate Gaussian Process Models for Stellar Activity

abstract: The radial velocity method is one of the most successful techniques for detecting exoplanets. It works by detecting the velocity of a host star induced by the gravitational effect of an orbiting planet, specifically the velocity along our line of sight, which is called the radial velocity of the star. Low-mass planets typically cause their host star to move with radial velocities of 1 m/s or less. By analyzing a time series of stellar spectra from a host star, modern astronomical instruments can in theory detect such planets. However, in practice, intrinsic stellar variability (e.g., star spots, convective motion, pulsations) affects the spectra and often mimics a radial velocity signal. This signal contamination makes it difficult to reliably detect low-mass planets. A principled approach to recovering planet radial velocity signals in the presence of stellar activity was proposed by Rajpaul et al. (2015). It uses a multivariate Gaussian process model to jointly capture time series of the apparent radial velocity and multiple indicators of stellar activity. We build on this work in two ways: (i) we propose using dimension reduction techniques to construct new high-information stellar activity indicators; and (ii) we extend the Rajpaul et al. (2015) model to a larger class of models and use a power-based model comparison procedure to select the best model. Despite significant interest in exoplanets, previous efforts have not performed large-scale stellar activity model selection or attempted to evaluate models based on planet detection power. In the case of main sequence G2V stars, we find that our method substantially improves planet detection power compared to previous state-of-the-art approaches.


title: Optimal entangling operations between deterministic blocks of qubits encoded into single photons

abstract: Here, we numerically simulate probabilistic elementary entangling operations between rail-encoded photons for the purpose of scalable universal quantum computation or communication. We propose grouping logical qubits into single-photon blocks wherein single-qubit rotations and the CNOT gate are fully deterministic and simple to implement. Inter-block communication is then allowed through said probabilistic entangling operations. We find a promising trend in the increasing probability of successful inter-block communication as we increase the number of optical modes operated on by our elementary entangling operations.


title: RF Wireless Power Transfer: Regreening Future Networks

abstract: Green radio communication is an emerging topic since the overall footprint of information and communication technology (ICT) services is predicted to triple between 2007 and 2020. Given this research line, energy harvesting (EH) and wireless power transfer (WPT) networks can be evaluated as promising approaches. In this paper, an overview of recent trends for future green networks on the platforms of EH and WPT is provided. By rethinking the application of radio frequency (RF)-WPT, a new concept, namely green RF-WTP, is introduced. Accordingly, opening challenges and promising combinations among current technologies, such as small-cell, millimeter (mm)-wave, and Internet of Things (IoT) networks, are discussed in details to seek solutions for the question with how to re-green the future networks?


title: Traveling foreshocks and transient foreshock phenomena

abstract: We use the multi-spacecraft capabilities of the Cluster and THEMIS missions to show that two types of foreshock may be detected in spacecraft data. One is the global foreshock that appears upstream of the Earth's quasi-parallel bow-shock under steady or variable interplanetary magnetic field. Another type is a traveling foreshock that is bounded by two rotational discontinuities in the interplanetary magnetic field and propagates along the bow-shock. Foreshock compressional boundaries are found at the edges of both types of foreshock. We show that isolated foreshock cavities are a subset of the traveling foreshock that form when two bounding rotational discontinuities are so close that the ultra-low frequency waves do not develop in the region between them. We also report observations of a spontaneous hot flow anomaly inside a traveling foreshock. This means that other phenomena, such as foreshock cavitons, may also exist inside this type of foreshock. In the second part of this work we present statistical properties of phenomena related to the foreshock, namely foreshock cavities, cavitons, spontaneous hot flow anomalies and foreshock compressional boundaries. We show that spontaneous hot flow anomalies are the most depleted transient structures in terms of the B-field and plasma density inside them and that the foreshock compressional boundaries and foreshock cavities are closely related structures.


title: A unified continuum and variational multiscale formulation for fluids, solids, and fluid-structure interaction

abstract: We develop a unified continuum modeling framework for viscous fluids and hyperelastic solids using the Gibbs free energy as the thermodynamic potential. This framework naturally leads to a pressure primitive variable formulation for the continuum body, which is well-behaved in both compressible and incompressible regimes. Our derivation also provides a rational justification of the isochoric-volumetric additive split of free energies in nonlinear continuum mechanics. The variational multiscale analysis is performed for the continuum model to construct a foundation for numerical discretization. We first consider the continuum body instantiated as a hyperelastic material and develop a variational multiscale formulation for the hyper-elastodynamic problem. The generalized-alpha method is applied for temporal discretization. A segregated algorithm for the nonlinear solver is designed and carefully analyzed. Second, we apply the new formulation to construct a novel unified formulation for fluid-solid coupled problems. The variational multiscale formulation is utilized for spatial discretization in both fluid and solid subdomains. The generalized-alpha method is applied for the whole continuum body, and optimal high-frequency dissipation is achieved in both fluid and solid subproblems. A new predictor multi-corrector algorithm is developed based on the segregated algorithm to attain a good balance between robustness and efficiency. The efficacy of the new formulations is examined in several benchmark problems. The results indicate that the proposed modeling and numerical methodologies constitute a promising technology for biomedical and engineering applications, particularly those necessitating incompressible models.

categories: physics.comp-ph cs.NA math.NA physics.flu-dyn


title: Constant Approximation for $k$-Median and $k$-Means with Outliers via Iterative Rounding

abstract: In this paper, we present a new iterative rounding framework for many clustering problems. Using this, we obtain an $(\alpha_1 + \epsilon \leq 7.081 + \epsilon)$-approximation algorithm for $k$-median with outliers, greatly improving upon the large implicit constant approximation ratio of Chen [Chen, SODA 2018]. For $k$-means with outliers, we give an $(\alpha_2+\epsilon \leq 53.002 + \epsilon)$-approximation, which is the first $O(1)$-approximation for this problem. The iterative algorithm framework is very versatile; we show how it can be used to give $\alpha_1$- and $(\alpha_1 + \epsilon)$-approximation algorithms for matroid and knapsack median problems respectively, improving upon the previous best approximations ratios of $8$ [Swamy, ACM Trans. Algorithms] and $17.46$ [Byrka et al, ESA 2015]. The natural LP relaxation for the $k$-median/$k$-means with outliers problem has an unbounded integrality gap. In spite of this negative result, our iterative rounding framework shows that we can round an LP solution to an almost-integral solution of small cost, in which we have at most two fractionally open facilities. Thus, the LP integrality gap arises due to the gap between almost-integral and fully-integral solutions. Then, using a pre-processing procedure, we show how to convert an almost-integral solution to a fully-integral solution losing only a constant-factor in the approximation ratio. By further using a sparsification technique, the additive factor loss incurred by the conversion can be reduced to any $\epsilon > 0$.


title: Growing a `Cosmic Beast': Observations and Simulations of MACS J0717.5+3745

abstract: We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.


title: GMC Collisions as Triggers of Star Formation. VI. Collision-Induced Turbulence

abstract: We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent, GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for injection of turbulent momentum into high density gas at high rates.


title: Tachyonic quench in a free bosonic field theory

abstract: We present a characterization of a bosonic field theory driven by a free (Gaussian) tachyonic Hamiltonian. This regime is obtained from a theory describing two coupled bosonic fields after a regular quench. Relevant physical quantities such as simple correlators, entanglement entropies, and the mutual information of disconnected subregions are computed. We show that the causal structure resembles a critical (massless) quench. For short times, physical quantities also resemble critical quenches. However, exponential divergences end up dominating the dynamics in a very characteristic way. This is related to the fact that the low-frequency modes do not equilibrate. Some applications and extensions are outlined.


title: Phototactic Supersmarticles

abstract: Smarticles, or smart active particles, are small robots equipped with only basic movement and sensing abilities that are incapable of rotating or displacing individually. We study the ensemble behavior of smarticles, i.e., the behavior a collective of these very simple computational elements can achieve, and how such behavior can be implemented using minimal programming. We show that an ensemble of smarticles constrained to remain close to one another (which we call a supersmarticle), achieves directed locomotion toward or away from a light source, a phenomenon known as phototaxing. We present experimental and theoretical models of phototactic supersmarticles that collectively move with a directed displacement in response to light. The motion of the supersmarticle is approximately Brownian, and is a result of chaotic interactions among smarticles. The system can be directed by introducing asymmetries among the individual smarticle's behavior, in our case by varying activity levels in response to light, resulting in supersmarticle biased motion.


title: An homotopy method for $\ell_p$ regression provably beyond self-concordance and in input-sparsity time

abstract: We consider the problem of linear regression where the $\ell_2^n$ norm loss (i.e., the usual least squares loss) is replaced by the $\ell_p^n$ norm. We show how to solve such problems up to machine precision in $O^*(n^{|1/2 - 1/p|})$ (dense) matrix-vector products and $O^*(1)$ matrix inversions, or alternatively in $O^*(n^{|1/2 - 1/p|})$ calls to a (sparse) linear system solver. This improves the state of the art for any $p\not\in \{1,2,+\infty\}$. Furthermore we also propose a randomized algorithm solving such problems in {\em input sparsity time}, i.e., $O^*(Z + \mathrm{poly}(d))$ where $Z$ is the size of the input and $d$ is the number of variables. Such a result was only known for $p=2$. Finally we prove that these results lie outside the scope of the Nesterov-Nemirovski's theory of interior point methods by showing that any symmetric self-concordant barrier on the $\ell_p^n$ unit ball has self-concordance parameter $\tilde{\Omega}(n)$.


title: Fast Path Localization on Graphs via Multiscale Viterbi Decoding

abstract: We consider a problem of localizing a path-signal that evolves over time on a graph. A path-signal can be viewed as the trajectory of a moving agent on a graph in several consecutive time points. Combining dynamic programming and graph partitioning, we propose a path-localization algorithm with significantly reduced computational complexity. We analyze the localization error for the proposed approach both in the Hamming distance and the destination's distance between the path estimate and the true path using numerical bounds. Unlike usual theoretical bounds that only apply to restricted graph models, the obtained numerical bounds apply to all graphs and all non-overlapping graph-partitioning schemes. In random geometric graphs, we are able to derive a closed-form expression for the localization error bound, and a tradeoff between localization error and the computational complexity. Finally, we compare the proposed technique with the maximum likelihood estimate under the path constraint in terms of computational complexity and localization error, and show significant speedup (100 times) with comparable localization error (4 times) on a graph from real data. Variants of the proposed technique can be applied to tracking, road congestion monitoring, and brain signal processing.


title: Asymptotics of extreme statistics of escape time in 1,2 and 3-dimensional diffusions

abstract: The first of $N$ identical independently distributed (i.i.d.) Brownian trajectories that arrives to a small target, sets the time scale of activation, which in general is much faster than the arrival to the target of only a single trajectory. Analytical asymptotic expressions for the minimal time is notoriously difficult to compute in general geometries. We derive here asymptotic laws for the probability density function of the first and second arrival times of a large number of i.i.d. Brownian trajectories to a small target in 1,2, and 3 dimensions and study their range of validity by stochastic simulations. The results are applied to activation of biochemical pathways in cellular transduction.


title: Colloidal particle electrorotation in a non-uniform electric field

abstract: A model to study the dynamics of colloidal particles in nonuniform electric fields is proposed. For an isolated sphere, the conditions and threshold for sustained (Quincke) rotation in a linear direct current (dc) field are determined. Particle dynamics becomes more complex with increasing electric field strength, changing from steady spinning around the particle center to time-dependent orbiting motion around the minimum field location. Pairs of particles exhibit intricate trajectories, which are combination of translation, due to dielectrophoresis, and rotation, due to the Quincke effect. Our model provides a basis to study the collective dynamics of many particles in a general electric field.


title: Molecular Hydrodynamics: Vortex Formation and Sound Wave Propagation

abstract: In the present study, quantitative feasibility tests of the hydrodynamic description of a two-dimensional fluid at the molecular level are performed, both with respect to length and time scales. Using high-resolution fluid velocity data obtained from extensive molecular dynamics simulations, we computed the transverse and longitudinal components of the velocity field by the Helmholtz decomposition and compared them with those obtained from the linearized Navier-Stokes (LNS) equations with time-dependent transport coefficients. By investigating the vortex dynamics and the sound wave propagation in terms of these field components, we confirm the validity of the LNS description for times comparable to or larger than several mean collision times. The LNS description still reproduces the transverse velocity field accurately at smaller times, but it fails to predict characteristic patterns of molecular origin visible in the longitudinal velocity field. Based on these observations, we validate the main assumptions of the mode-coupling approach. The assumption that the velocity autocorrelation function can be expressed in terms of the fluid velocity field and the tagged particle distribution is found to be remarkably accurate even for times comparable to or smaller than the mean collision time. This suggests that the hydrodynamic-mode description remains valid down to the molecular scale.


title: Learning to Bid Without Knowing your Value

abstract: We address online learning in complex auction settings, such as sponsored search auctions, where the value of the bidder is unknown to her, evolving in an arbitrary manner and observed only if the bidder wins an allocation. We leverage the structure of the utility of the bidder and the partial feedback that bidders typically receive in auctions, in order to provide algorithms with regret rates against the best fixed bid in hindsight, that are exponentially faster in convergence in terms of dependence on the action space, than what would have been derived by applying a generic bandit algorithm and almost equivalent to what would have been achieved in the full information setting. Our results are enabled by analyzing a new online learning setting with outcome-based feedback, which generalizes learning with feedback graphs. We provide an online learning algorithm for this setting, of independent interest, with regret that grows only logarithmically with the number of actions and linearly only in the number of potential outcomes (the latter being very small in most auction settings). Last but not least, we show that our algorithm outperforms the bandit approach experimentally and that this performance is robust to dropping some of our theoretical assumptions or introducing noise in the feedback that the bidder receives.


title: Two Error Bounds of Imperfect Binary Search

abstract: Suppose we know that an object is in a sorted table and we want to determine the index of that object. To achieve this goal we could perform a binary search. However, suppose it is time-consuming to determine the relative position of that object to any other objects in the table. In this scenario, we might want to resort to an incomplete solution: we could device an algorithm that quickly predicts the result of comparing two objects, and replace the actual comparison with this algorithm during a binary search. The question then is how far away are the results yielded by the imperfect binary search from the correct answers. We present two quick lemmas that answer this question.


title: Differentially Private ANOVA Testing

abstract: Modern society generates an incredible amount of data about individuals, and releasing summary statistics about this data in a manner that provably protects individual privacy would offer a valuable resource for researchers in many fields. We present the first algorithm for analysis of variance (ANOVA) that preserves differential privacy, allowing this important statistical test to be conducted (and the results released) on databases of sensitive information. In addition to our private algorithm for the F test statistic, we show a rigorous way to compute p-values that accounts for the added noise needed to preserve privacy. Finally, we present experimental results quantifying the statistical power of this differentially private version of the test, finding that a sample of several thousand observations is frequently enough to detect variation between groups. The differentially private ANOVA algorithm is a promising approach for releasing a common test statistic that is valuable in fields in the sciences and social sciences.


title: A Study of Optimizing Heterogeneous Resources for Open IoT

abstract: Recently, IoT technologies have been progressed, and many sensors and actuators are connected to networks. Previously, IoT services were developed by vertical integration style. But now Open IoT concept has attracted attentions which achieves various IoT services by integrating horizontal separated devices and services. For Open IoT era, we have proposed the Tacit Computing technology to discover the devices with necessary data for users on demand and use them dynamically. We also implemented elemental technologies of Tacit Computing. In this paper, we propose three layers optimizations to reduce operation cost and improve performance of Tacit computing service, in order to make as a continuous service of discovered devices by Tacit Computing. In optimization process, appropriate function allocations are calculated for device, network and cloud layer before full-scale operation.


title: On Spaces with the Maximal Number of Conformal Killing Vectors

abstract: It is natural to expect and simple to prove that every conformally flat space possess the maximal number of conformal Killing vector fields (CKVs). On the other hand, it is interesting to ask whether the converse is true. Is conformal flatness a necessary condition for the existence of the maximal number of CKVs? In this review article it is proven that the answer is yes, a space admits the maximal number of CKVs if, and only if, it is conformally flat.


title: Probing the structure of Kepler ZZ Ceti stars with full evolutionary models-based asteroseismology

abstract: We present an asteroseismological analysis of four ZZ Ceti stars observed with \emph{Kepler}: GD 1212, SDSS J113655.17+040952.6, KIC 11911480 and KIC 4552982, based on a grid of full evolutionary models of DA white dwarf stars. We employ a grid of carbon-oxygen core white dwarfs models, characterized by a detailed and consistent chemical inner profile for the core and the envelope. In addition to the observed periods, we take into account other information from the observational data, as amplitudes, rotational splittings and period spacing, as well as photometry and spectroscopy. For each star, we present an asteroseismological model that closely reproduce their observed properties. The asteroseismological stellar mass and effective temperature of the target stars are (0.632 +/- 0.027 Msun, 10737 +/- 73 K) for GD 1212, (0.745 +/- 0.007 Msun, 11110 +/- 69 K) for KIC 4552982, (0.5480 +/- 0.01 Msun, 12721 +/- 228 K) for KIC1191480 and (0.570 +/- 0.01 Msun, 12060 +/- 300 K) for SDSS J113655.17+040952.6. In general, the asteroseismological values are in good agreement with the spectroscopy. For KIC 11911480 and SDSS J113655.17+040952.6 we derive a similar seismological mass, but the hydrogen envelope is an order of magnitude thinner for SDSS J113655.17+040952.6, that is part of a binary system and went through a common envelope phase.


title: Binary Linear Codes with Optimal Scaling: Polar Codes with Large Kernels

abstract: We prove that, for the binary erasure channel (BEC), the polar-coding paradigm gives rise to codes that not only approach the Shannon limit but do so under the best possible scaling of their block length as a~function of the gap to capacity. This result exhibits the first known family of binary codes that attain both optimal scaling and quasi-linear complexity of encoding and decoding. Our proof is based on the construction and analysis of binary polar codes with large kernels. When communicating reliably at rates within $\varepsilon > 0$ of capacity, the code length $n$ often scales as $O(1/\varepsilon^{\mu})$, where the constant $\mu$ is called the scaling exponent. It is known that the optimal scaling exponent is $\mu=2$, and it is achieved by random linear codes. The scaling exponent of conventional polar codes (based on the $2\times 2$ kernel) on the BEC is $\mu=3.63$. This falls far short of the optimal scaling guaranteed by random codes. Our main contribution is a rigorous proof of the following result: for the BEC, there exist $\ell\times\ell$ binary kernels, such that polar codes constructed from these kernels achieve scaling exponent $\mu(\ell)$ that tends to the optimal value of $2$ as $\ell$ grows. We furthermore characterize precisely how large $\ell$ needs to be as a function of the gap between $\mu(\ell)$ and $2$. The resulting binary codes maintain the recursive structure of conventional polar codes, and thereby achieve construction complexity $O(n)$ and encoding/decoding complexity $O(n\log n)$.


title: Algebras of diagonal operators of the form scalar-plus-compact are Calkin algebras

abstract: For every Banach space $X$ with a Schauder basis consider the Banach algebra $\mathbb{R} I\oplus\mathcal{K}_\mathrm{diag}(X)$ of all diagonal operators that are of the form $\lambda I + K$. We prove that $\mathbb{R} I\oplus\mathcal{K}_\mathrm{diag}(X)$ is a Calkin algbra i.e., there exists a Banach space $\mathcal{Y}_X$ so that the Calkin algebra of $\mathcal{Y}_X$ is isomorphic as a Banach algebra to $\mathbb{R} I\oplus\mathcal{K}_\mathrm{diag}(X)$. Among other applications of this theorem we obtain that certain hereditarily indecomposable spaces and the James spaces $J_p$ and their duals endowed with natural multiplications are Calkin algebras, that all non-reflexive Banach spaces with unconditional bases are isomorphic as Banach spaces to Calkin algebras, and that sums of reflexive spaces with unconditional bases with certain James-Tsirelson type spaces are isomorphic as Banach spaces to Calkin algebras.


title: Generalized Linear Model Regression under Distance-to-set Penalties

abstract: Estimation in generalized linear models (GLM) is complicated by the presence of constraints. One can handle constraints by maximizing a penalized log-likelihood. Penalties such as the lasso are effective in high dimensions, but often lead to unwanted shrinkage. This paper explores instead penalizing the squared distance to constraint sets. Distance penalties are more flexible than algebraic and regularization penalties, and avoid the drawback of shrinkage. To optimize distance penalized objectives, we make use of the majorization-minimization principle. Resulting algorithms constructed within this framework are amenable to acceleration and come with global convergence guarantees. Applications to shape constraints, sparse regression, and rank-restricted matrix regression on synthetic and real data showcase strong empirical performance, even under non-convex constraints.


title: Solar System Science with ESA Euclid

abstract: The ESA Euclid mission has been designed to map the geometry of the dark Universe. Scheduled for launch in 2020, it will conduct a six-years visible and NIR imaging and spectroscopic survey over 15,000 deg 2 down to mag~24.5. Although the survey will avoid low ecliptic latitudes, the survey pattern in repeated sequences of four broad-band filters seems well-adapted to Solar System objects (SSOs) detection and characterization. We aim at evaluating Euclid capability to discover SSOs, and measure their position, apparent magnitude, and SED. Also, we investigate how these measurements can lead to the determination of their orbits, morphology, physical properties, and surface composition. We use current census of SSOs to estimate the number of SSOs detectable by Euclid. Then we estimate how Euclid will constrain the SSOs dynamical, physical, and compositional properties. With current survey design, about 150,000 SSOs, mainly from the asteroid main-belt, should be observed by Euclid. These objects will all have high inclination. There is a potential for discovery of several 10,000 SSOs, in particular KBOs at high declination. Euclid observations will refine the spectral classification of SSOs by extending the spectral coverage provided by, e.g. Gaia and the LSST to 2 microns. The time-resolved photometry, combined with sparse photometry will contribute to the determination of SSO rotation period, spin orientation, and shape model. The sharp and stable point-spread function of Euclid will also allow to resolve KBO binary systems and detect activity around Centaurs. The depth of Euclid survey, its spectral coverage, and observation cadence has great potential for Solar System research. A dedicated processing for SSOs is being set in place to produce catalogs of astrometry, multi-color and time-resolved photometry, and spectral classification of some 10$^5$ SSOs, delivered as Legacy Science.


title: Accelerating Training of Deep Neural Networks via Sparse Edge Processing

abstract: We propose a reconfigurable hardware architecture for deep neural networks (DNNs) capable of online training and inference, which uses algorithmically pre-determined, structured sparsity to significantly lower memory and computational requirements. This novel architecture introduces the notion of edge-processing to provide flexibility and combines junction pipelining and operational parallelization to speed up training. The overall effect is to reduce network complexity by factors up to 30x and training time by up to 35x relative to GPUs, while maintaining high fidelity of inference results. This has the potential to enable extensive parameter searches and development of the largely unexplored theoretical foundation of DNNs. The architecture automatically adapts itself to different network sizes given available hardware resources. As proof of concept, we show results obtained for different bit widths.


title: Implications for planetary system formation from interstellar object 1I/2017 U1 (`Oumuamua)

abstract: The recently discovered minor body 1I/2017 U1 (`Oumuamua) is the first known object in our Solar System that is not bound by the Sun's gravity. Its hyperbolic orbit (eccentricity greater than unity) strongly suggests that it originated outside our Solar System; its red color is consistent with substantial space weathering experienced over a long interstellar journey. We carry out an simple calculation of the probability of detecting such an object. We find that the observed detection rate of 1I-like objects can be satisfied if the average mass of ejected material from nearby stars during the process of planetary formation is ~20 Earth masses, similar to the expected value for our Solar System. The current detection rate of such interstellar interlopers is estimated to be 0.2/year, and the expected number of detections over the past few years is almost exactly one. When the Large Synoptic Survey Telescope begins its wide, fast, deep all-sky survey the detection rate will increase to 1/year. Those expected detections will provide further constraints on nearby planetary system formation through a better estimate of the number and properties of interstellar objects.


title: Computationally efficient cardiac views projection using 3D Convolutional Neural Networks

abstract: 4D Flow is an MRI sequence which allows acquisition of 3D images of the heart. The data is typically acquired volumetrically, so it must be reformatted to generate cardiac long axis and short axis views for diagnostic interpretation. These views may be generated by placing 6 landmarks: the left and right ventricle apex, and the aortic, mitral, pulmonary, and tricuspid valves. In this paper, we propose an automatic method to localize landmarks in order to compute the cardiac views. Our approach consists of first calculating a bounding box that tightly crops the heart, followed by a landmark localization step within this bounded region. Both steps are based on a 3D extension of the recently introduced ENet. We demonstrate that the long and short axis projections computed with our automated method are of equivalent quality to projections created with landmarks placed by an experienced cardiac radiologist, based on a blinded test administered to a different cardiac radiologist.


title: Reflection Spectra of the Black Hole Binary Candidate MAXI J1535-571 in the Hard State Observed by NuSTAR

abstract: We report on a NuSTAR observation of the recently discovered bright black hole candidate MAXI J1535-571. NuSTAR observed the source on MJD 58003 (five days after the outburst was reported). The spectrum is characteristic of a black hole binary in the hard state. We observe clear disk reflection features, including a broad Fe K$\alpha$ line and a Compton hump peaking around 30 keV. Detailed spectral modeling reveals narrow Fe K$\alpha$ line complex centered around 6.5 keV on top of the strong relativistically broadened Fe K$\alpha$ line. The narrow component is consistent with distant reflection from moderately ionized material. The spectral continuum is well described by a combination of cool thermal disk photons and a Comptonized plasma with the electron temperature $kT_{\rm e}=19.7\pm{0.4}$ keV. An adequate fit can be achieved for the disk reflection features with a self-consistent relativistic reflection model that assumes a lamp-post geometry for the coronal illuminating source. The spectral fitting measures a black hole spin $a>0.84$, inner disk radius $R_{\rm in}<2.01~r_{\rm ISCO}$, and a lamp-post height $h=7.2^{+0.8}_{-2.0} r_{\rm g}$ (statistical errors, 90% confidence), indicating no significant disk truncation and a compact corona. Although the distance and mass of this source are not currently known, this suggests the source was likely in the brighter phases of the hard state during this NuSTAR observation.


title: Merging Cluster Collaboration: Optical and Spectroscopic Survey of a Radio-Selected Sample of Twenty Nine Merging Galaxy Clusters

abstract: Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system.


title: Automatic Differentiation for Tensor Algebras

abstract: Kjolstad et. al. proposed a tensor algebra compiler. It takes expressions that define a tensor element-wise, such as $f_{ij}(a,b,c,d) = \exp\left[-\sum_{k=0}^4 \left((a_{ik}+b_{jk})^2\, c_{ii} + d_{i+k}^3 \right) \right]$, and generates the corresponding compute kernel code. For machine learning, especially deep learning, it is often necessary to compute the gradient of a loss function $l(a,b,c,d)=l(f(a,b,c,d))$ with respect to parameters $a,b,c,d$. If tensor compilers are to be applied in this field, it is necessary to derive expressions for the derivatives of element-wise defined tensors, i.e. expressions for $(da)_{ik}=\partial l/\partial a_{ik}$. When the mapping between function indices and argument indices is not 1:1, special attention is required. For the function $f_{ij} (x) = x_i^2$, the derivative of the loss is $(dx)_i=\partial l/\partial x_i=\sum_j (df)_{ij}2x_i$; the sum is necessary because index $j$ does not appear in the indices of $f$. Another example is $f_{i}(x)=x_{ii}^2$, where $x$ is a matrix; here we have $(dx)_{ij}=\delta_{ij}(df)_i2x_{ii}$; the Kronecker delta is necessary because the derivative is zero for off-diagonal elements. Another indexing scheme is used by $f_{ij}(x)=\exp x_{i+j}$; here the correct derivative is $(dx)_{k}=\sum_i (df)_{i,k-i} \exp x_{k}$, where the range of the sum must be chosen appropriately. In this publication we present an algorithm that can handle any case in which the indices of an argument are an arbitrary linear combination of the indices of the function, thus all the above examples can be handled. Sums (and their ranges) and Kronecker deltas are automatically inserted into the derivatives as necessary. Additionally, the indices are transformed, if required (as in the last example). The algorithm outputs a symbolic expression that can be subsequently fed into a tensor algebra compiler. Source code is provided.


title: Non-Wire Alternatives to Capacity Expansion

abstract: Distributed energy resources (DERs) can serve as non-wire alternatives to capacity expansion by managing peak load to avoid or defer traditional expansion projects. In this paper, we study a planning problem that co-optimizes DERs investment and operation (e.g., energy efficiency, energy storage, demand response, solar photovoltaic) and the timing of capacity expansion. We formulate the problem as a large scale (in the order of millions of variables because we model the operation of DERs over a period of decades) non-convex optimization problem. Despite its non-convexities, we find its optimal solution by decomposing it using the Dantzig-Wolfe Decomposition Algorithm and solving a series of small linear problems. Finally, we present a real planning problem at the University of Washington Seattle Campus.


title: HESS J1741-302: a hidden accelerator in the Galactic plane

abstract: The H.E.S.S. collaboration has discovered a new very high energy (VHE, E $>$ 0.1 TeV) $\gamma$-ray source, HESS J1741-302, located in the Galactic plane. Despite several attempts to constrain its nature, no plausible counterpart has been found so far at X-ray and MeV/GeV $\gamma$-ray energies, and the source remains unidentified. An analysis of 145-hour of observations of HESS J1741-302 at VHEs has revealed a steady and relatively weak TeV source ($\sim$1$\%$ of the Crab Nebula flux), with a spectral index of $\Gamma$ = 2.3 $\pm$ 0.2$_{\text{stat}}$ $\pm$ 0.2$_{\text{sys}}$, extending to energies up to 10 TeV without any clear signature of a cut-off. In a hadronic scenario, such a spectrum implies an object with particle acceleration up to energies of several hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size of HESS J1741-302 is compatible with the H.E.S.S. point spread function at VHEs, with an extension constrained to be below 0.068$^{\circ}$ at a 99$\%$ confidence level. The $\gamma$-ray emission detected by H.E.S.S. can be explained both within a hadronic scenario, due to collisions of protons with energies of hundreds of TeV with dense molecular clouds, and in a leptonic scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged (20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio source 1LC 358.266+0.038 found to be spatially coincident with the best fit position of HESS J1741-302, is also envisaged.


title: Uplink Performance Analysis of a Drone Cell in a Random Field of Ground Interferers

abstract: Aerial base stations are a promising technology to increase the capabilities of the existing communication networks. However, the existing analytical frameworks do not sufficiently characterize the impact of ground interferers on the aerial base stations. In order to address this issue, we model the effect of interference coming from the coexisting ground networks on the aerial link, which could be the uplink of an aerial cell served by a drone base station. By considering a Poisson field of ground interferers, we characterize the aggregate interference experienced by the drone. This result includes the effect of the drone antenna pattern, the height-dependent shadowing, and various types of environment. We show that the benefits that a drone obtains from a better line-of-sight (LoS) at high altitudes is counteracted by a high vulnerability to the interference coming from the ground. However, by deriving the link coverage probability and transmission rate we show that a drone base station is still a promising technology if the overall system is properly dimensioned according to the given density and transmission power of the interferers. Particularly, our results illustrate how the benefits of such network is maximized by defining the optimal drone altitude and signal-to- interference (SIR) requirement.


title: Experimentally Quantifying Anion Polarizability at the Air/Water Interface

abstract: The adsorption of anions from aqueous solution on the air/water interface controls important heterogeneous chemistry in the atmosphere and is thought to have similar physics to anion adsorption at hydrophobic interfaces more generally. While much theoretical and experimental work over the last 20 years has made clear that the adsorption of large, polarizable anions is thermodynamically favorable, determining the role of polarizability in adsorption has proven surprisingly challenging. Simple physical models make clear that nonpolarizable anions will not adsorb, but trends in adsorption are difficult to rationalize based on polarizability and in some theoretical approaches adsorption is observed without change in anion polarization. Because there are no experimental studies of interfacial anion polarizability, one possible explanation of these results is that theoretical descriptions suffer from systematic error. In this study we use interface specific nonlinear optical spectroscopy to extract the spectral response of the interfacial ClO$_{4}^{-}$ anion. We find that (i) the interfacial environment induces a break in symmetry of the anion due to its solvation anisotropy (ii) the Raman depolarization ratio, a measure of the ratio of the change in two components of ClO$_{4}^{-}$'s polarizability tensor with change in nuclear positions, is $>2\times$ larger than its value in the adjoining bulk phase and interfacial concentration dependent and (iii) using a simple theoretical description we find that our measured changes in depolarization ratio are consistent with known changes in surface potential and tension with small changes in bulk ClO$_{4}^{-}$ concentration. The notion that interfacial anion polarizability differs from that in bulk, and that this polarizability is coverage dependent, is not accounted for in any current theoretical treatment of ions at the air/water interface.


title: Decentralised firewall for malware detection

abstract: This paper describes the design and development of a decentralized firewall system powered by a novel malware detection engine. The firewall is built using blockchain technology. The detection engine aims to classify Portable Executable (PE) files as malicious or benign. File classification is carried out using a deep belief neural network (DBN) as the detection engine. Our approach is to model the files as grayscale images and use the DBN to classify those images into the aforementioned two classes. An extensive data set of 10,000 files is used to train the DBN. Validation is carried out using 4,000 files previously unexposed to the network. The final result of whether to allow or block a file is obtained by arriving at a proof of work based consensus in the blockchain network.


title: Collective behavior of active topological solitons, knotted streamlines, and transport of cargo in liquid crystals

abstract: Active colloids and liquid crystals are capable of locally converting the macroscopically-supplied energy into directional motion and promise a host of new applications, ranging from drug delivery to cargo transport at the mesoscale. Here we uncover how topological solitons in liquid crystals can locally transform electric energy to translational motion and allow for the transport of cargo along directions dependent on frequency of the applied electric field. By combining polarized optical video microscopy and numerical modeling that reproduces both the equilibrium structures of solitons and their temporal evolution in applied fields, we uncover the physical underpinnings behind this reconfigurable motion and study how it depends on the structure and topology of solitons. We show that, unexpectedly, the directional motion of solitons with and without the cargo arises mainly from the asymmetry in rotational dynamics of molecular ordering in liquid crystal rather than from the asymmetry of fluid flows, as in conventional active soft matter systems.


title: Counting Roots of Polynomials Over Prime Power Rings

abstract: Suppose $p$ is a prime, $t$ is a positive integer, and $f\!\in\!\mathbb{Z}[x]$ is a univariate polynomial of degree $d$ with coefficients of absolute value $<\!p^t$. We show that for any fixed $t$, we can compute the number of roots in $\mathbb{Z}/(p^t)$ of $f$ in deterministic time $(d+\log p)^{O(1)}$. This fixed parameter tractability appears to be new for $t\!\geq\!3$. A consequence for arithmetic geometry is that we can efficiently compute Igusa zeta functions $Z$, for univariate polynomials, assuming the degree of $Z$ is fixed.


title: Dunkl Operators for Arbitrary Finite Groups

abstract: The Dunkl operators associated to a necessarily finite Coxeter group acting on a Euclidean space are generalized to any finite group using the techniques of non-commutative geometry, as introduced by the authors to view the usual Dunkl operators as covariant derivatives in a quantum principal bundle with a quantum connection. The definitions of Dunkl operators and their corresponding Dunkl connections are generalized to quantum principal bundles over quantum spaces which possess a classical finite structure group. We introduce cyclic Dunkl connections and their cyclic Dunkl operators. Then we establish a number of interesting properties of these structures, including the characteristic zero curvature property. Particular attention is given to the example of complex reflection groups, and their naturally generalized siblings called groups of Coxeter type.


title: Reconstructing Video from Interferometric Measurements of Time-Varying Sources

abstract: Very long baseline interferometry (VLBI) makes it possible to recover images of astronomical sources with extremely high angular resolution. Most recently, the Event Horizon Telescope (EHT) has extended VLBI to short millimeter wavelengths with a goal of achieving angular resolution sufficient for imaging the event horizons of nearby supermassive black holes. VLBI provides measurements related to the underlying source image through a sparse set spatial frequencies. An image can then be recovered from these measurements by making assumptions about the underlying image. One of the most important assumptions made by conventional imaging methods is that over the course of a night's observation the image is static. However, for quickly evolving sources, such as the galactic center's supermassive black hole (Sgr A*) targeted by the EHT, this assumption is violated and these conventional imaging approaches fail. In this work we propose a new way to model VLBI measurements that allows us to recover both the appearance and dynamics of an evolving source by reconstructing a video rather than a static image. By modeling VLBI measurements using a Gaussian Markov Model, we are able to propagate information across observations in time to reconstruct a video, while simultaneously learning about the dynamics of the source's emission region. We demonstrate our proposed Expectation-Maximization (EM) algorithm, StarWarps, on realistic synthetic observations of black holes, and show how it substantially improves results compared to conventional imaging algorithms. Additionally, we demonstrate StarWarps on real VLBI data of the M87 Jet from the VLBA.


title: Strengthening Convex Relaxations of 0/1-Sets Using Boolean Formulas

abstract: In convex integer programming, various procedures have been developed to strengthen convex relaxations of sets of integer points. On the one hand, there exist several general-purpose methods that strengthen relaxations without specific knowledge of the set $ S $, such as popular linear programming or semi-definite programming hierarchies. On the other hand, various methods have been designed for obtaining strengthened relaxations for very specific sets that arise in combinatorial optimization. We propose a new efficient method that interpolates between these two approaches. Our procedure strengthens any convex set $ Q \subseteq \mathbb{R}^n $ containing a set $ S \subseteq \{0,1\}^n $ by exploiting certain additional information about $ S $. Namely, the required extra information will be in the form of a Boolean formula $ \phi $ defining the target set $ S $. The aim of this work is to analyze various aspects regarding the strength of our procedure. As one result, interpreting an iterated application of our procedure as a hierarchy, our findings simplify, improve, and extend previous results by Bienstock and Zuckerberg on covering problems.


title: Magellan/PFS Radial Velocities of GJ 9827, a late K dwarf at 30 pc with Three Transiting Super-Earths

abstract: The Kepler mission showed us that planets with sizes between that of Earth and Neptune appear to be the most common type in our Galaxy. These "super-Earths" continue to be of great interest for exoplanet formation, evolution, and composition studies. However, the number of super-Earths with well-constrained mass and radius measurements remains small (40 planets with $\sigma_{\rm{mass}}<$ 25\%), due in part to the faintness of their host stars causing ground-based mass measurements to be challenging. Recently, three transiting super-Earth planets were detected by the K2 mission around the nearby star GJ 9827/HIP 115752, at only 30 pc away. The radii of the planets span the "radius gap"' detected by Fulton et al. (2017), and all orbit within ~6.5 days, easing follow-up observations. Here we report radial velocity (RV) observations of GJ 9827, taken between 2010 and 2016 with the Planet Finder Spectrograph on the Magellan II Telescope. We employ two different RV analysis packages, SYSTEMIC and RadVel, to derive masses and thus densities of the GJ 9827 planets. We also test a Gaussian Process regression analysis, but find the correlated stellar noise is not well constrained by the PFS data, and that the GP tends to over fit the RV semi-amplitudes resulting in a lower K value. Our RV observations are not able to place strong mass constraints on the two outer planets (c & d) but do indicate that planet b, at 1.64 R$_{\oplus}$ and ~8 M$_{\oplus}$, is one of the most massive (and dense) super-Earth planets detected to date.


title: Dynamical Freezing in a Spin Glass System with Logarithmic Correlations

abstract: We consider a continuous time random walk on the two-dimensional discrete torus, whose motion is governed by the discrete Gaussian free field on the corresponding box acting as a potential. More precisely, at any vertex the walk waits an exponentially distributed time with mean given by the exponential of the field and then jumps to one of its neighbors, chosen uniformly at random. We prove that throughout the low-temperature regime and at in-equilibrium timescales, the process admits a scaling limit as a spatial K-process driven by a random trapping landscape, which is explicitly related to the limiting extremal process of the field. Alternatively, the limiting process is a supercritical Liouville Brownian motion with respect to the continuum Gaussian free field on the box. This demonstrates rigorously and for the first time, as far as we know, a dynamical freezing in a spin glass system with logarithmically correlated energy levels.


title: An Optimal Distributed $(\Delta+1)$-Coloring Algorithm?

abstract: Vertex coloring is one of the classic symmetry breaking problems studied in distributed computing. In this paper we present a new algorithm for $(\Delta+1)$-list coloring in the randomized ${\sf LOCAL}$ model running in $O(\mathsf{Det}_{\scriptscriptstyle d}(\text{poly} \log n))$ time, where $\mathsf{Det}_{\scriptscriptstyle d}(n')$ is the deterministic complexity of $(\text{deg}+1)$-list coloring on $n'$-vertex graphs. (In this problem, each $v$ has a palette of size $\text{deg}(v)+1$.) This improves upon a previous randomized algorithm of Harris, Schneider, and Su [STOC'16, JACM'18] with complexity $O(\sqrt{\log \Delta} + \log\log n + \mathsf{Det}_{\scriptscriptstyle d}(\text{poly} \log n))$, and, for some range of $\Delta$, is much faster than the best known deterministic algorithm of Fraigniaud, Heinrich, and Kosowski [FOCS'16] and Barenboim, Elkin, and Goldenberg [PODC'18], with complexity $O(\sqrt{\Delta\log \Delta}\log^\ast \Delta + \log^* n)$. Our algorithm "appears to be" optimal, in view of the $\Omega(\mathsf{Det}(\text{poly} \log n))$ randomized lower bound due to Chang, Kopelowitz, and Pettie [FOCS'16], where $\mathsf{Det}$ is the deterministic complexity of $(\Delta+1)$-list coloring. At present, the best upper bounds on $\mathsf{Det}_{\scriptscriptstyle d}(n')$ and $\mathsf{Det}(n')$ are both $2^{O(\sqrt{\log n'})}$ and use a black box application of network decompositions (Panconesi and Srinivasan [Journal of Algorithms'96]). It is quite possible that the true complexities of both problems are the same, asymptotically, which would imply the randomized optimality of our $(\Delta+1)$-list coloring algorithm.


title: "Attention" for Detecting Unreliable News in the Information Age

abstract: An Unreliable news is any piece of information which is false or misleading, deliberately spread to promote political, ideological and financial agendas. Recently the problem of unreliable news has got a lot of attention as the number instances of using news and social media outlets for propaganda have increased rapidly. This poses a serious threat to society, which calls for technology to automatically and reliably identify unreliable news sources. This paper is an effort made in this direction to build systems for detecting unreliable news articles. In this paper, various NLP algorithms were built and evaluated on Unreliable News Data 2017 dataset. Variants of hierarchical attention networks (HAN) are presented for encoding and classifying news articles which achieve the best results of 0.944 ROC-AUC. Finally, Attention layer weights are visualized to understand and give insight into the decisions made by HANs. The results obtained are very promising and encouraging to deploy and use these systems in the real world to mitigate the problem of unreliable news.


title: The Green function for waves on the $2$-regular Bethe lattice

abstract: In this paper, we compute an explicit analytic expression for the Green function of the wave operator on the $2$-regular lattice called the "Bethe lattice" equipped with its standard metric. In particular, we exhibit a phenomena of abnormal speed of propagation for waves: the effective speed of propagation of energy for large time is $c_*=2\sqrt 2/3 <1$, and there exists a true propagation at any speed $c<c_*$.


title: A Faster Distributed Single-Source Shortest Paths Algorithm

abstract: We devise new algorithms for the single-source shortest paths (SSSP) problem with non-negative edge weights in the CONGEST model of distributed computing. While close-to-optimal solutions, in terms of the number of rounds spent by the algorithm, have recently been developed for computing SSSP approximately, the fastest known exact algorithms are still far away from matching the lower bound of $ \tilde \Omega (\sqrt{n} + D) $ rounds by Peleg and Rubinovich [SIAM Journal on Computing 2000], where $ n $ is the number of nodes in the network and $ D $ is its diameter. The state of the art is Elkin's randomized algorithm [STOC 2017] that performs $ \tilde O(n^{2/3} D^{1/3} + n^{5/6}) $ rounds. We significantly improve upon this upper bound with our two new randomized algorithms for polynomially bounded integer edge weights, the first performing $ \tilde O (\sqrt{n D}) $ rounds and the second performing $ \tilde O (\sqrt{n} D^{1/4} + n^{3/5} + D) $ rounds. Our bounds also compare favorably to the independent result by Ghaffari and Li [STOC 2018]. As side results, we obtain a $ (1 + \epsilon) $-approximation $ \tilde O ((\sqrt{n} D^{1/4} + D) / \epsilon) $-round algorithm for directed SSSP and a new work/depth trade-off for exact SSSP on directed graphs in the PRAM model.


title: A generalized MBO diffusion generated motion for orthogonal matrix-valued fields

abstract: We consider the problem of finding stationary points of the Dirichlet energy for orthogonal matrix-valued fields. Following the Ginzburg-Landau approach, this energy is relaxed by penalizing the matrix-valued field when it does not take orthogonal matrix values. A generalization of the MBO diffusion generated motion is introduced that effectively finds local minimizers of this energy by iterating two steps until convergence. In the first step, as in the original method, the current matrix-valued field is evolved by the diffusion equation. In the second step, the field is pointwise reassigned to the closest orthogonal matrix, which can be computed via the singular value decomposition. We extend the Lyapunov function of Esedoglu and Otto to show that the method is non-increasing on iterates and hence, unconditionally stable. We also prove that spatially discretized iterates converge to a stationary solution in a finite number of iterations. The algorithm is implemented using the closest point method and non-uniform fast Fourier transform. We conclude with several numerical experiments on flat tori and closed surfaces, which, unsurprisingly, exhibit classical behavior from the Allen-Cahn and complex Ginzburg Landau equations, but also new phenomena.


title: Sequential two-fold Pearson chi-squared test and tails of the Bessel process distributions

abstract: We find asymptotic formulas for error probabilities of two-fold Pearson goodness-of-fit test as functions of two critical levels. These results may be reformulated in terms of tails of two-dimensional distributions of the Bessel process. Necessary properties of the Infeld function are obtained.


title: Proximal Alternating Penalty Algorithms for Constrained Convex Optimization

abstract: We develop two new proximal alternating penalty algorithms to solve a wide range class of constrained convex optimization problems. Our approach mainly relies on a novel combination of the classical quadratic penalty, alternating minimization, Nesterov's acceleration, and adaptive strategy for parameters. The first algorithm is designed to solve generic and possibly nonsmooth constrained convex problems without requiring any Lipschitz gradient continuity or strong convexity, while achieving the best-known $O(\frac{1}{k})$-convergence rate in a non-ergodic sense, where $k$ is the iteration counter. The second algorithm is also designed to solve non-strongly convex, but semi-strongly convex problems. This algorithm can achieve the best-known $O(\frac{1}{k^2})$-convergence rate on the primal constrained problem. Such a rate is obtained in two cases: (i)~averaging only on the iterate sequence of the strongly convex term, or (ii) using two proximal operators of this term without averaging. In both algorithms, we allow one to linearize the second subproblem to use the proximal operator of the corresponding objective term. Then, we customize our methods to solve different convex problems, and lead to new variants. As a byproduct, these algorithms preserve the same convergence guarantees as in our main algorithms. We verify our theoretical development via different numerical examples and compare our methods with some existing state-of-the-art algorithms.