title: Kullback-Leibler Divergence-Based Fuzzy $C$-Means Clustering Incorporating Morphological Reconstruction and Wavelet Frames for Image Segmentation

abstract: Although spatial information of images usually enhance the robustness of the Fuzzy C-Means (FCM) algorithm, it greatly increases the computational costs for image segmentation. To achieve a sound trade-off between the segmentation performance and the speed of clustering, we come up with a Kullback-Leibler (KL) divergence-based FCM algorithm by incorporating a tight wavelet frame transform and a morphological reconstruction operation. To enhance FCM's robustness, an observed image is first filtered by using the morphological reconstruction. A tight wavelet frame system is employed to decompose the observed and filtered images so as to form their feature sets. Considering these feature sets as data of clustering, an modified FCM algorithm is proposed, which introduces a KL divergence term in the partition matrix into its objective function. The KL divergence term aims to make membership degrees of each image pixel closer to those of its neighbors, which brings that the membership partition becomes more suitable and the parameter setting of FCM becomes simplified. On the basis of the obtained partition matrix and prototypes, the segmented feature set is reconstructed by minimizing the inverse process of the modified objective function. To modify abnormal features produced in the reconstruction process, each reconstructed feature is reassigned to the closest prototype. As a result, the segmentation accuracy of KL divergence-based FCM is further improved. What's more, the segmented image is reconstructed by using a tight wavelet frame reconstruction operation. Finally, supporting experiments coping with synthetic, medical and color images are reported. Experimental results exhibit that the proposed algorithm works well and comes with better segmentation performance than other comparative algorithms. Moreover, the proposed algorithm requires less time than most of the FCM-related algorithms.


title: Nondegeneracy of solutions for a critical Hartree equation

abstract: The aim of this paper is to prove the nondegeneracy of the unique positive solutions for the following critical Hartree type equations when $\mu>0$ is close to $0$, $$ -\Delta u=\left(I_{\mu}\ast u^{2^{\ast}_{\mu}}\right)u^{{2}^{\ast}_{\mu}-1},~~x\in\mathbb{R}^{N}, $$ where $ I_{\mu}(x)=\frac{\Gamma(\frac{\mu}{2})}{\Gamma(\frac{{N-\mu}}{2})\pi^{\frac{N}{2}}2^{{N-\mu}}|x|^{\mu}} $ is the Riesz potential and $2^{\ast}_{\mu}=\frac{2{N-\mu}}{N-2}$ is the upper critical exponent due to the Hardy-Littlewood-Sobolev inequality.


title: TFApprox: Towards a Fast Emulation of DNN Approximate Hardware Accelerators on GPU

abstract: Energy efficiency of hardware accelerators of deep neural networks (DNN) can be improved by introducing approximate arithmetic circuits. In order to quantify the error introduced by using these circuits and avoid the expensive hardware prototyping, a software emulator of the DNN accelerator is usually executed on CPU or GPU. However, this emulation is typically two or three orders of magnitude slower than a software DNN implementation running on CPU or GPU and operating with standard floating point arithmetic instructions and common DNN libraries. The reason is that there is no hardware support for approximate arithmetic operations on common CPUs and GPUs and these operations have to be expensively emulated. In order to address this issue, we propose an efficient emulation method for approximate circuits utilized in a given DNN accelerator which is emulated on GPU. All relevant approximate circuits are implemented as look-up tables and accessed through a texture memory mechanism of CUDA capable GPUs. We exploit the fact that the texture memory is optimized for irregular read-only access and in some GPU architectures is even implemented as a dedicated cache. This technique allowed us to reduce the inference time of the emulated DNN accelerator approximately 200 times with respect to an optimized CPU version on complex DNNs such as ResNet. The proposed approach extends the TensorFlow library and is available online at https://github.com/ehw-fit/tf-approximate.


title: Measurements of ZnWO$_4$ anisotropic response to nuclear recoils for the ADAMO project

abstract: Anisotropic scintillators can offer a unique possibility to exploit the so-called directionality approach in order to investigate the presence of those Dark Matter (DM) candidates inducing nuclear recoils. In fact, their use can overcome the difficulty of detecting extremely short nuclear recoil traces. In this paper we present recent measurements performed on the anisotropic response of a ZnWO$_4$ crystal scintillator to nuclear recoils, in the framework of the ADAMO project. The anisotropic features of the ZnWO$_4$ crystal scintillators were initially measured with $\alpha$ particles; those results have been also confirmed by the additional measurements presented here. The experimental nuclear recoil data were obtained by using a neutron generator at ENEA-CASACCIA and neutron detectors to tag the scattered neutrons; in particular, the quenching factor values for nuclear recoils along different crystallographic axes have been determined for three different neutron scattering angles (i.e. nuclear recoils energies). From these measurements, the anisotropy of the light response for nuclear recoils in the ZnWO$_4$ crystal scintillator has been determined at 5.4 standard deviations.


title: A reliable description of the radial oscillations of compact stars

abstract: We develop a numerical algorithm for the solution of the Sturm-Liouville differential equation governing the stationary radial oscillations of nonrotating compact stars. Our method is based on the Numerov's method that turns the Sturm-Liouville differential equation in an eigenvalue problem. In our development we provide a strategy to correctly deal with the star boundaries and the interfaces between layers with different mechanical properties. Assuming that the fluctuations obey the same equation of state of the background, we analyze various different stellar models and we precisely determine hundreds of eigenfrequencies and of eigenmodes. If the equation of state does not present an interface discontinuity, the fundamental radial eigenmode becomes unstable exactly at the critical central energy density corresponding to the largest gravitational mass. However, in the presence of an interface discontinuity, there exist stable configurations with a central density exceeding the critical one and with a smaller gravitational mass.


title: A note on Stein equation for weighted sums of independent $\chi^{2}$ distributions

abstract: This note provides the Stein equation for weighted sums of independent $\chi^{2}$ distributions.


title: The Four Dimensions of Social Network Analysis: An Overview of Research Methods, Applications, and Software Tools

abstract: Social network based applications have experienced exponential growth in recent years. One of the reasons for this rise is that this application domain offers a particularly fertile place to test and develop the most advanced computational techniques to extract valuable information from the Web. The main contribution of this work is three-fold: (1) we provide an up-to-date literature review of the state of the art on social network analysis (SNA);(2) we propose a set of new metrics based on four essential features (or dimensions) in SNA; (3) finally, we provide a quantitative analysis of a set of popular SNA tools and frameworks. We have also performed a scientometric study to detect the most active research areas and application domains in this area. This work proposes the definition of four different dimensions, namely Pattern & Knowledge discovery, Information Fusion & Integration, Scalability, and Visualization, which are used to define a set of new metrics (termed degrees) in order to evaluate the different software tools and frameworks of SNA (a set of 20 SNA-software tools are analyzed and ranked following previous metrics). These dimensions, together with the defined degrees, allow evaluating and measure the maturity of social network technologies, looking for both a quantitative assessment of them, as to shed light to the challenges and future trends in this active area.


title: On shuffle-type functional relations of desingularized multiple zeta-functions

abstract: We treat desingularized multiple zeta-functions introduced by Furusho, Komori, Matsumoto and Tsumura. In this paper, we prove functional relations, which are shuffle type product formulae, between desingularized multiple zeta-functions and desingularized values.


title: Quasi-periodic traveling gravity-capillary waves

abstract: We present a numerical study of spatially quasi-periodic traveling waves on the surface of an ideal fluid of infinite depth. This is a generalization of the classic Wilton ripple problem to the case when the ratio of wave numbers satisfying the dispersion relation is irrational. We propose a conformal mapping formulation of the water wave equations that employs a quasi-periodic variant of the Hilbert transform to compute the normal velocity of the fluid from its velocity potential on the free surface. We develop a Fourier pseudo-spectral discretization of the traveling water wave equations in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on the torus. This leads to an overdetermined nonlinear least squares problem that we solve using a variant of the Levenberg-Marquardt method. We investigate various properties of quasi-periodic traveling waves, including Fourier resonances, time evolution in conformal space on the torus, asymmetric wave crests, capillary wave patterns that change from one gravity wave trough to the next without repeating, and the dependence of wave speed and surface tension on the amplitude parameters that describe a two-parameter family of waves.


title: Optimal Randomized First-Order Methods for Least-Squares Problems

abstract: We provide an exact analysis of a class of randomized algorithms for solving overdetermined least-squares problems. We consider first-order methods, where the gradients are pre-conditioned by an approximation of the Hessian, based on a subspace embedding of the data matrix. This class of algorithms encompasses several randomized methods among the fastest solvers for least-squares problems. We focus on two classical embeddings, namely, Gaussian projections and subsampled randomized Hadamard transforms (SRHT). Our key technical innovation is the derivation of the limiting spectral density of SRHT embeddings. Leveraging this novel result, we derive the family of normalized orthogonal polynomials of the SRHT density and we find the optimal pre-conditioned first-order method along with its rate of convergence. Our analysis of Gaussian embeddings proceeds similarly, and leverages classical random matrix theory results. In particular, we show that for a given sketch size, SRHT embeddings exhibits a faster rate of convergence than Gaussian embeddings. Then, we propose a new algorithm by optimizing the computational complexity over the choice of the sketching dimension. To our knowledge, our resulting algorithm yields the best known complexity for solving least-squares problems with no condition number dependence.


title: Quantifying Interference-Assisted Signal Strength Surveillance of Sound Vibrations

abstract: A malicious attacker could, by taking control of internet-of-things devices, use them to capture received signal strength (RSS) measurements and perform surveillance on a person's vital signs, activities, and sound in their environment. This article considers an attacker who looks for subtle changes in the RSS in order to eavesdrop sound vibrations. The challenge to the adversary is that sound vibrations cause very low amplitude changes in RSS, and RSS is typically quantized with a significantly larger step size. This article contributes a lower bound on an attacker's monitoring performance as a function of the RSS step size and sampling frequency so that a designer can understand their relationship. Our bound considers the little-known and counter-intuitive fact that an adversary can improve their sinusoidal parameter estimates by making some devices transmit to add interference power into the RSS measurements. We demonstrate this capability experimentally. As we show, for typical transceivers, the RSS surveillance attacker can monitor sound vibrations with remarkable accuracy. New mitigation strategies will be required to prevent RSS surveillance attacks.


title: A Possible Distance Bias for Type Ia Supernovae with Different Ejecta Velocities

abstract: After correcting for their light-curve shape and color, Type Ia supernovae (SNe Ia) are precise cosmological distance indicators. However, there remains a non-zero intrinsic scatter in the differences between measured distance and that inferred from a cosmological model (i.e., Hubble residuals or HRs), indicating that SN Ia distances can potentially be further improved. We use the open-source relational database kaepora to generate composite spectra with desired average properties of phase, light-curve shape, and HR. At many phases, the composite spectra from two subsamples with positive and negative average HRs are significantly different. In particular, in all spectra from 9 days before to 15 days after peak brightness, we find that SNe with negative HRs have, on average, higher ejecta velocities (as seen in nearly every optical spectral feature) than SNe with positive HRs. At +4 days relative to B-band maximum, using a sample of 62 SNe Ia, we measure a 0.091 +/- 0.035 mag HR step between SNe with Si II 6355 line velocities higher/lower than -11,000 km/s (the median velocity). After light-curve shape and color correction, SNe with higher velocities tend to have underestimated distance moduli relative to a cosmological model. The intrinsic scatter in our sample reduces from 0.094 mag to 0.082 mag after making this correction. Using the Si II 6355 velocity evolution of 115 SNe Ia, we estimate that a velocity difference > 500 km/s exists at each epoch between the positive-HR and negative-HR samples with 99.4% confidence. Finally at epochs later than +37 days, we observe that negative-HR composite spectra tend to have weaker spectral features in comparison to positive-HR composite spectra.


title: Gravitational wave population inference with deep flow-based generative network

abstract: We combine hierarchical Bayesian modeling with a flow-based deep generative network, in order to demonstrate that one can efficiently constraint numerical gravitational wave (GW) population models at a previously intractable complexity. Existing techniques for comparing data to simulation,such as discrete model selection and Gaussian process regression, can only be applied efficiently to moderate-dimension data. This limits the number of observable (e.g. chirp mass, spins.) and hyper-parameters (e.g. common envelope efficiency) one can use in a population inference. In this study, we train a network to emulate a phenomenological model with 6 observables and 4 hyper-parameters, use it to infer the properties of a simulated catalogue and compare the results to the phenomenological model. We find that a 10-layer network can emulate the phenomenological model accurately and efficiently. Our machine enables simulation-based GW population inferences to take on data at a new complexity level.


title: Constructing d-log integrands and computing master integrals for three-loop four-particle scattering

abstract: We compute all master integrals for massless three-loop four-particle scattering amplitudes required for processes like di-jet or di-photon production at the LHC. We present our result in terms of a Laurent expansion of the integrals in the dimensional regulator up to 8$^{\text{th}}$ power, with coefficients expressed in terms of harmonic polylogarithms. As a basis of master integrals we choose integrals with integrands that only have logarithmic poles - called $d$log forms. This choice greatly facilitates the subsequent computation via the method of differential equations. We detail how this basis is obtained via an improved algorithm originally developed by one of the authors. We provide a public implementation of this algorithm. We explain how the algorithm is naturally applied in the context of unitarity. In addition, we classify our $d$log forms according to their soft and collinear properties.


title: Orbital evidences for dark-matter-free Milky Way dwarf spheroidal galaxies

abstract: The nature of Milky Way dwarf spheroidals (MW dSphs) has been questioned, in particular whether they are dominated by dark matter (DM). Here we investigate an alternative scenario, for which tidal shocks are exerted by the MW to DM-free dSphs after a first infall of their gas-rich progenitors, and for which theoretical calculations have been verified by pure N-body simulations. Whether or not the dSphs are on their first infall cannot be resolved on the sole basis of their star formation history. In fact, gas removal may cause complex gravitational instabilities and near-pericenter passages can give rise to tidal disruptive processes. Advanced precision with the Gaia satellite in determining both their past orbital motions and the MW velocity curve is, however, providing crucial results. First, tidal shocks explain why DM-free dSphs are found preferentially near their pericenter, where they are in a destructive process, while their chance to be long-lived satellites is associated with a very low probability P~ 2 10^-7, which is at odds with the current DM-dominated dSph scenario. Second, most dSph binding energies are consistent with a first infall. Third, the MW tidal shocks that predict the observed dSph velocity dispersions are themselves predicted in amplitude by the most accurate MW velocity curve. Fourth, tidal shocks accurately predict the forces or accelerations exerted at half-light radius of dSphs, including the MW and the Magellanic System gravitational attractions. The above is suggestive of dSphs that are DM-free and tidally shocked near their pericenters, which may provoke a significant quake in our understanding of near-field cosmology.


title: Understanding the large inferred Einstein radii of observed low-mass galaxy clusters

abstract: We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe.


title: The Stars in M15 Were Born with the r-process

abstract: High-resolution spectroscopy of stars on the red giant branch (RGB) of the globular cluster M15 has revealed a large (~1 dex) dispersion in the abundances of r-process elements, like Ba and Eu. Neutron star mergers (NSMs) have been proposed as a major source of the r-process. However, most NSM models predict a delay time longer than the timescale for cluster formation. One possibility is that a NSM polluted the surfaces of stars in M15 long after the cluster finished forming. In this case, the abundances of the polluting elements would decrease in the first dredge-up as stars turn on to the RGB. We present Keck/DEIMOS abundances of Ba in 66 stars along the entire RGB and the top of the main sequence. The Ba abundances have no trend with stellar luminosity (evolutionary phase). Therefore, the stars were born with the Ba they have today, and Ba did not originate in a source with a delay time longer than the timescale for cluster formation. In particular, if the source of Ba was a neutron star merger, it would have had a very short delay time. Alternatively, if Ba enrichment took place before the formation of the cluster, an inhomogeneity of a factor of 30 in Ba abundance needs to be able to persist over the length scale of the gas cloud that formed M15, which is unlikely.


title: From birth associations to field stars: mapping the small-scale orbit distribution in the Galactic disc

abstract: Stars born at the same time in the same place should have formed from gas of the same element composition. But most stars subsequently disperse from their birth siblings, in orbit and orbital phase, becoming 'field stars'. Here we explore and provide direct observational evidence for this process in the Milky Way disc, by quantifying the probability that orbit-similarity among stars implies indistinguishable metallicity. We define the orbit similarity among stars through their distance in action-angle space, $\Delta (J,\theta)$, and their abundance similarity simply by $\Delta$[Fe/H]. Analyzing a sample of main sequence stars from Gaia DR2 and LAMOST, we find an excess of pairs with the same metallicity ($\Delta\mathrm{[Fe/H]}<0.1$) that extends to remarkably large separations in $\Delta (J,\theta)$ that correspond to nearly 1 kpc distances. We assess the significance of this effect through a mock sample, drawn from a smooth and phase-mixed orbit distribution. Through grouping such star pairs into associations with a friend-of-friends algorithm linked by $\Delta (J,\theta)$, we find 100s of mono-abundance groups with $\ge 3$ (to $\gtrsim 20$) members; these groups -- some clusters, some spread across the sky -- are over an order-of-magnitude more abundant than expected for a smooth phase-space distribution, suggesting that we are witnessing the 'dissolution' of stellar birth associations into the field.


title: The Edge of the Galaxy

abstract: We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc.


title: ALMA Observations of the Molecular Clouds in NGC 625

abstract: We present the highest resolution (1") $^{12}$CO(1-0) observations of molecular gas in the dwarf starburst galaxy NGC 625, to date, obtained with ALMA. Molecular gas is distributed in discrete clouds within an area of $0.4$ kpc$^2$ and does not have well-ordered large-scale motions. We measure a total molecular mass in NGC 625 of $5.3\times 10^6 M_\odot$, assuming a Milky Way CO-to-H$_2$ conversion factor. We use the CPROPS package to identify molecular clouds and measure their properties. The 19 resolved CO clouds have a median radius of 20 pc, a median line width 2.5 km s$^{-1}$, and a median surface density of 169 $M_\odot$pc$^{-2}$. Larson scaling relations suggest that molecular clouds in NGC 625 are mostly in virial equilibrium. Comparison of our high-resolution CO observations with a star formation rate map, inferred from ancillary optical observations observations, suggests that about 40% of the molecular clouds coincide with the brightest HII regions. These bright HII regions have a range of molecular gas depletion timescales, all within a factor of $\sim 3$ of the global depletion time in NGC 625 of 106-134 Myr. The highest surface density molecular clouds towards the southwest of the galaxy, in a region we call the Butterfly, do not show strong star formation activity and suggest a depletion time scale longer than 5 Gyr.


title: The regulation of galaxy growth along the size-mass relation by star-formation, as traced by H${\alpha}$ in KMOS$^{3D}$ galaxies at 0.7 < z < 2.7

abstract: We present half-light sizes measured from H${\alpha}$ emission tracing star-formation in 281 star-forming galaxies from the KMOS3D survey at 0.7 < z < 2.7. Sizes are derived by fitting 2D exponential disk models, with bootstrap errors averaging 20%. H${\alpha}$ sizes are a median (mean) of 1.19 (1.26) times larger than those of the stellar continuum, which due to radial dust gradients places an upper limit on the growth in stellar size via star formation, with just 43% intrinsic scatter. At fixed continuum size the H${\alpha}$ size shows no residual trend with stellar mass, star formation rate, redshift or morphology. The only significant residual trend is with the excess obscuration of H${\alpha}$ by dust, at fixed continuum obscuration. The scatter in continuum size at fixed stellar mass is likely driven by the scatter in halo spin parameters. The stability of the ratio of H${\alpha}$ size to continuum size demonstrates a high degree of stability in halo spin and in the transfer of angular momentum to the disk over a wide range of physical conditions and cosmic time. This may require local regulation by feedback processes. The implication of our results, as we demonstrate using a toy model, is that our upper limit on star-formation driven growth is sufficient only to evolve star-forming galaxies approximately along the observed size-mass relation, consistent with the size growth of galaxies at constant cumulative co-moving number density. To explain the observed evolution of the size-mass relation of star-forming disk galaxies other processes, such as the preferential quenching of compact galaxies or galaxy mergers, may be required.


title: A pre-metric generalization of the Lorentz transformation

abstract: The concept of an observer and their associated rest space is defined in a pre-metric (i.e., projective-geometric) context that relates to time+space decompositions of the tangent bundle to space-time. The transformation from one observer to another when the two are in a state of relative motion is then defined, and its relationship to the Lorentz transformation is discussed. The group of all linear transformations that preserve the observer quadric, which generalizes the proper-time hyperboloid in Minkowski space, is defined and the reductions to some of its subgroups are described, as well as its extension to the group that preserves the fundamental quadric, which generalizes the light cone.


title: Generation of chiral asymmetry via helical magnetic fields

abstract: It is well known that helical magnetic fields undergo a so-called inverse cascade by which their correlation length grows due to the conservation of magnetic helicity in classical ideal magnetohydrodynamics (MHD). At high energies above approximately $10$ MeV, however, classical MHD is necessarily extended to chiral MHD and then the conserved quantity is $\langle\mathcal{H}\rangle + 2 \langle\mu_5\rangle / \lambda$ with $\langle\mathcal{H}\rangle$ being the mean magnetic helicity and $\langle\mu_5\rangle$ being the mean chiral chemical potential of charged fermions. Here, $\lambda$ is a (phenomenological) chiral feedback parameter. In this paper, we study the evolution of the chiral MHD system with the initial condition of nonzero $\langle\mathcal{H}\rangle$ and vanishing $\mu_5$. We present analytic derivations for the time evolution of $\langle\mathcal{H}\rangle$ and $\langle\mu_5\rangle$ that we compare to a series of laminar and turbulent three-dimensional direct numerical simulations. We find that the late-time evolution of $\langle\mathcal{H}\rangle$ depends on the magnetic and kinetic Reynolds numbers ${\rm Re}_{_\mathrm{M}}$ and ${\rm Re}_{_\mathrm{K}}$. For a high ${\rm Re}_{_\mathrm{M}}$ and ${\rm Re}_{_\mathrm{K}}$ where turbulence occurs, $\langle\mathcal{H}\rangle$ eventually evolves in the same way as in classical ideal MHD where the inverse correlation length of the helical magnetic field scales with time $t$ as $k_\mathrm{p} \propto t^{-2/3}$. For a low Reynolds numbers where the velocity field is negligible, the scaling is changed to $k_\mathrm{p} \propto t^{-1/2}\mathrm{ln}\left(t/t_\mathrm{log}\right)$. After being rapidly generated, $\langle\mu_5\rangle$ always decays together with $k_\mathrm{p}$, i.e. $\langle\mu_5\rangle \approx k_\mathrm{p}$, with a time evolution that depends on whether the system is in the limit of low or high Reynolds numbers.

categories: physics.plasm-ph astro-ph.CO astro-ph.HE cond-mat.soft hep-ph


title: Illuminating black hole subsystems in young star clusters

abstract: There is increasing evidence that globular clusters retain sizeable black hole populations at present day. This is supported by dynamical simulations of cluster evolution, which have unveiled the spatial distribution and mass spectrum of black holes in clusters across cosmic age. However, black hole populations of young, high metallicity clusters remain unconstrained. Black holes hosted by these clusters mass segregate early in their evolutionary history, forming central subsystems of hundreds to thousands of black holes. We argue that after supernova feedback has subsided ($\gtrsim 50\,{\rm Myr}$), the host cluster can accumulate gas from its dense surroundings, from which the black hole subsystem accretes at highly enhanced rates. The collective accretion luminosity can be substantial and provides a novel observational constraint for young massive clusters. We test this hypothesis by performing 3D hydrodynamic simulations where we embed discretized potentials, representing our black holes, within the potential of a massive cluster. This system moves supersonically with respect to a gaseous medium from which it accretes. We study the accretion of this black hole subsystem for different subsystem populations and determine the integrated accretion luminosity of the black hole subsystem. We apply our results to the young massive clusters of the Antennae Galaxies and find that a typical subsystem accretion luminosity should be in excess of $\approx 10^{40}\,{\rm ergs\,\,s^{-1}}$. We argue that no strong candidates of this luminous signal have been observed and constrain the subsystem population of a typical cluster in the Antennae Galaxies to $\lesssim10-2\times10^2$ $10\,M_\odot$ black holes, given that feedback doesn't significantly impede accretion and that the gas remains optically thin.


title: Multiple Chemodynamic Stellar Populations of the Ursa Minor Dwarf Spheroidal Galaxy

abstract: We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.


title: Robust Numerical Tracking of One Path of a Polynomial Homotopy on Parallel Shared Memory Computers

abstract: We consider the problem of tracking one solution path defined by a polynomial homotopy on a parallel shared memory computer. Our robust path tracker applies Newton's method on power series to locate the closest singular parameter value. On top of that, it computes singular values of the Hessians of the polynomials in the homotopy to estimate the distance to the nearest different path. Together, these estimates are used to compute an appropriate adaptive stepsize. For n-dimensional problems, the cost overhead of our robust path tracker is O(n), compared to the commonly used predictor-corrector methods. This cost overhead can be reduced by a multithreaded program on a parallel shared memory computer.


title: Estimating Q(s,s') with Deep Deterministic Dynamics Gradients

abstract: In this paper, we introduce a novel form of value function, $Q(s, s')$, that expresses the utility of transitioning from a state $s$ to a neighboring state $s'$ and then acting optimally thereafter. In order to derive an optimal policy, we develop a forward dynamics model that learns to make next-state predictions that maximize this value. This formulation decouples actions from values while still learning off-policy. We highlight the benefits of this approach in terms of value function transfer, learning within redundant action spaces, and learning off-policy from state observations generated by sub-optimal or completely random policies. Code and videos are available at http://sites.google.com/view/qss-paper.


title: Modelling the 3D Climate of Venus with OASIS

abstract: Flexible 3D models to explore the vast diversity of terrestrial planets and interpret observational data are still in their early stages. In this work, we present OASIS: a novel and flexible 3D virtual planet laboratory. With OASIS we envision a platform that couples self-consistently seven individual modules representing the main physical and chemical processes that shape planetary environments. Additionally, OASIS is capable of producing simulated spectra from different instruments and observational techniques. In this work we focus on the benchmark test of coupling four of the physical modules: fluid dynamics, radiation, turbulence and surface/soil. To test the OASIS platform, we produced 3D simulations of the Venus climate and its atmospheric circulation and study how the modeled atmosphere changes with various cloud covers, atmospheric heat capacity, and surface friction. 3D simulations of Venus are challenging because they require long integration times with a computationally expensive radiative transfer code. By comparing OASIS results with observational data, we verify that the new model is able to successfully simulate Venus. With simulated spectra produced directly from the 3D simulations, we explore the capabilities of future missions, like LUVOIR, to observe Venus analogs located at a distance of 10 pc. With OASIS, we have taken the first steps to build a sophisticated and very flexible platform capable of studying the environment of terrestrial planets, which will be an essential tool to characterize observed terrestrial planets and plan future observations.


title: Molecular beam epitaxy preparation and in situ characterization of FeTe thin films

abstract: We have synthesized Fe$_{1+y}$Te thin films by means of molecular beam epitaxy (MBE) under Te-limited growth conditions. We found that epitaxial layer-by-layer growth is possible for a wide range of excess Fe values, wider than expected from what is known on the bulk material. Using x-ray magnetic circular dichroism spectroscopy at the Fe L$_{2,3}$ and Te M$_{4,5}$ edges, we observed that films with high excess Fe contain ferromagnetic clusters while films with lower excess Fe remain nonmagnetic. Moreover, x-ray absorption spectroscopy showed that it is possible to obtain films with very similar electronic structure as that of a high quality bulk single crystal Fe$_{1.14}$Te. Our results suggest that MBE with Te-limited growth may provide an opportunity to synthesize FeTe films with smaller amounts of excess Fe as to come closer to a possible superconducting phase.


title: The Pre-Explosion Mass Distribution of Hydrogen-Poor Superluminous Supernova Progenitors and New Evidence for a Mass-Spin Correlation

abstract: Despite indications that superluminous supernovae (SLSNe) originate from massive progenitors, the lack of a uniformly analyzed statistical sample has so far prevented a detailed view of the progenitor mass distribution. Here we present and analyze the pre-explosion mass distribution of hydrogen-poor SLSN progenitors as determined from uniformly modelled light curves of 62 events. We construct the distribution by summing the ejecta mass posteriors of each event, using magnetar light curve models presented in our previous works (and using a nominal neutron star remnant mass). The resulting distribution spans $3.6-40$ M$_{\odot}$, with a sharp decline at lower masses, and is best fit by a broken power law described by ${\rm d}N/{\rm dlog}M \propto M^{-0.41 \pm 0.06}$ at $3.6-8.6$ M$_{\odot}$ and $\propto M^{-1.26 \pm 0.06}$ at $8.6-40$ M$_{\odot}$. We find that observational selection effects cannot account for the shape of the distribution. Relative to Type Ib/c SNe, the SLSN mass distribution extends to much larger masses and has a different power-law shape, likely indicating that the formation of a magnetar allows more massive stars to explode as some of the rotational energy accelerates the ejecta. Comparing the SLSN distribution with predictions from single and binary star evolution models, we find that binary models for a metallicity of $Z\lesssim 1/3$ Z$_{\odot}$ are best able to reproduce its broad shape, in agreement with the preference of SLSNe for low metallicity environments. Finally, we uncover a correlation between the pre-explosion mass and the magnetar initial spin period, where SLSNe with low masses have slower spins, a trend broadly consistent with the effects of angular momentum transport evident in models of rapidly-rotating carbon-oxygen stars.


title: Gowers norms for automatic sequences

abstract: We show that any automatic sequence can be separated into a structured part and a Gowers uniform part in a way that is considerably more efficient than guaranteed by the Arithmetic Regularity Lemma. For sequences produced by strongly connected and prolongable automata, the structured part is rationally almost periodic, while for general sequences the description is marginally more complicated. In particular, we show that all automatic sequences orthogonal to periodic sequences are Gowers uniform. As an application, we obtain for any $l \geq 2$ and any automatic set $A \subset \mathbb{N}_0$ lower bounds on the number of $l$-term arithmetic progressions - contained in $A$ - with a given difference. The analogous result is false for general subsets of $\mathbb{N}_0$ and progressions of length $\geq 5$.


title: Spin pumping and large field-like torque at room temperature in sputtered amorphous WTe$_{2-x}$ films

abstract: We studied the spin-to-charge and charge-to-spin conversion at room temperature in sputtered WTe2-x (x=0.8)(t)/Co20Fe60B20(6 nm) heterostructures. Spin pumping measurements were used to characterize the spin-to-charge efficiency and the spin efficiency was calculated to be larger than ~0.035. Second harmonic Hall measurements were carried out to estimate the charge-to-spin conversion ratio. We found that the system exhibits a large field-like-torque (spin torque efficiency ~ 0.1) and small damping-like-torque (spin torque efficiency ~0.001) compared to that reported for heavy metals. High-resolution transmission electron microscopy images show that the WTe2-x layer is amorphous, which may enhance the spin swapping effect by inducing large interfacial spin orbit scattering, thus contributing to a large field like torque.


title: Chronofold: a data structure for versioned text

abstract: Chronofold is a replicated data structure for versioned text. It is designed for use in collaborative editors and revision control systems. Past models of this kind either retrofitted local linear orders to a distributed system (the OT approach) or employed distributed data models locally (the CRDT approach). That caused either extreme fragility in a distributed setting or egregious overheads in local use. Overall, that local/distributed impedance mismatch is cognitively taxing and causes lots of complexity. We solve that by using subjective linear orders locally at each replica, while inter-replica communication uses a distributed model. A separate translation layer insulates local data structures from the distributed environment. We modify the Lamport timestamping scheme to make that translation as trivial as possible. We believe our approach has applications beyond the domain of collaborative editing.


title: On The Microscopic Modeling of Vehicular Traffic on General Networks

abstract: We introduce a formalism to deal with the microscopic modeling of vehicular traffic on a road network. Traffic on each road is uni-directional, and the dynamics of each vehicle is described by a Follow-the-Leader model. From a mathematical point of view, this amounts to define a system of ordinary differential equations on an arbitrary network. A general existence and uniqueness result is provided, while priorities at junctions are shown to hinder the stability of solutions. We investigate the occurrence of the Braess paradox in a time-dependent setting within this model. The emergence of Nash equilibria in a non-stationary situation results in the appearance of Braess type paradoxes, and this is supported by numerical simulations.


title: Knowledge Transfer between Buildings for Seismic Damage Diagnosis through Adversarial Learning

abstract: Automated structural damage diagnosis after earthquakes is important for improving the efficiency of disaster response and rehabilitation. In conventional data-driven frameworks which use machine learning or statistical models, structural damage diagnosis models are often constructed using supervised learning. The supervised learning requires historical structural response data and corresponding damage states (i.e., labels) for each building to learn the building-specific damage diagnosis model. However, in post-earthquake scenarios, historical data with labels are often not available for many buildings in the affected area. This makes it difficult to construct a damage diagnosis model. Further, directly using the historical data from other buildings to construct a damage diagnosis model for the target building would lead to inaccurate results. This is because each building has unique physical properties and thus unique data distribution. To this end, we introduce a new framework to transfer the model learned from other buildings to diagnose structural damage states in the target building without any labels. This framework is based on an adversarial domain adaptation approach that extracts domain-invariant feature representations of data from different buildings. The feature extraction function is trained in an adversarial way, which ensures that the extracted feature distributions are robust to changes in structures while being predictive of the damage states. With the extracted domain-invariant feature representations, the data distributions become consistent across different buildings. We evaluate our framework on both numerical simulation and field data collected from multiple building structures, which outperforms the state-of-the-art benchmark methods.


title: A polynomial time algorithm for studying physical observables in chaotic eigenstates

abstract: We introduce an algorithm, the Orthogonal Operator Polynomial Expansion (OOPEX), to approximately compute expectation values in energy eigenstates at finite energy density of non-integrable quantum many-body systems with polynomial effort, whereas exact diagonalization (ED) of the Hamiltonian $H$ is exponentially hard. The OOPEX relies on the eigenstate thermalization hypothesis, which conjectures that eigenstate expectation values of physical observables in such systems vary smoothly with the eigenstate energy (and other macroscopic conserved quantities, if any), and computes them through a series generated by repeated multiplications, rather than diagonalization, of $H$ and whose successive terms oscillate faster with the energy. The hypothesis guarantees that only the first few terms of this series contribute appreciably. We further show that the OOPEX, in a sense, is the most optimum algorithm based on series expansions of $H$ as it avoids computing the many-body density of states which plagues other similar algorithms. Then, we argue non-rigorously that working in the Fock space of operators, rather than that of states as is usually done, yields convergent results with computational resources that scale polynomially with $N$. We demonstrate the polynomial scaling by applying the OOPEX to the non-integrable Ising chain and comparing with ED and high-temperature expansion (HTX) results. The OOPEX provides access to much larger $N$ than ED and HTX do, which facilitates overcoming finite-size effects that plague the other methods to extract correlation lengths in chaotic eigenstates. In addition, access to large systems allows testing a recent conjecture that the Renyi entropy of chaotic eigenstates has positive curvature if the Renyi index $>1$, and we find encouraging supporting evidence.

categories: cond-mat.stat-mech cond-mat.str-el hep-th physics.comp-ph


title: Petrophysically and geologically guided multi-physics inversion using a dynamic Gaussian mixture model

abstract: In a previous paper, we introduced a framework for carrying out petrophysically and geologically guided geophysical inversions. In that framework, petrophysical and geological information is modelled with a Gaussian Mixture Model (GMM). In the inversion, the GMM serves as a prior for the geophysical model. The formulation was confined to problems in which a single physical property model was sought, with a single geophysical dataset. In this paper, we extend that framework to jointly invert multiple geophysical datasets that depend on multiple physical properties. The petrophysical and geological information is used to couple geophysical surveys that, otherwise, rely on independent physics. This requires advancements in two areas. First, an extension from a univariate to a multivariate analysis of the petrophysical data, and their inclusion within the inverse problem, is necessary. Second, we address the practical issues of simultaneously inverting data from multiple surveys and finding a solution that acceptably reproduces each one, along with the petrophysical and geological information. To illustrate the efficacy of our approach and the advantages of carrying out multi-physics inversions, we invert synthetic gravity and magnetic data associated with a kimberlite deposit. The kimberlite pipe contains two distinct facies embedded in a host rock. Inverting the datasets individually leads to a binary geological model: background or kimberlite. A multi-physics inversion, with petrophysical information, differentiates between the two main kimberlite facies of the pipe. Through this example, we also highlight the capabilities of our framework to work with interpretive geologic assumptions when minimal quantitative information is available. In those cases, the dynamic updates of the Gaussian Mixture Model allow us to perform multi-physics inversions by learning a petrophysical model.


title: Minimax-Optimal Off-Policy Evaluation with Linear Function Approximation

abstract: This paper studies the statistical theory of batch data reinforcement learning with function approximation. Consider the off-policy evaluation problem, which is to estimate the cumulative value of a new target policy from logged history generated by unknown behavioral policies. We study a regression-based fitted Q iteration method, and show that it is equivalent to a model-based method that estimates a conditional mean embedding of the transition operator. We prove that this method is information-theoretically optimal and has nearly minimal estimation error. In particular, by leveraging contraction property of Markov processes and martingale concentration, we establish a finite-sample instance-dependent error upper bound and a nearly-matching minimax lower bound. The policy evaluation error depends sharply on a restricted $\chi^2$-divergence over the function class between the long-term distribution of the target policy and the distribution of past data. This restricted $\chi^2$-divergence is both instance-dependent and function-class-dependent. It characterizes the statistical limit of off-policy evaluation. Further, we provide an easily computable confidence bound for the policy evaluator, which may be useful for optimistic planning and safe policy improvement.


title: The Blue Compact Dwarf Galaxy VCC848 Formed by Dwarf-Dwarf Merging

abstract: It has long been speculated that many starburst or compact dwarf galaxies are resulted from dwarf-dwarf galaxy merging, but unequivocal evidence for this possibility has rarely been reported in the literature. We present the first study of deep optical broadband images of a gas-dominated blue compact dwarf galaxy (BCD) VCC848 (Mstar=2e8Msun) which hosts extended stellar shells and thus is confirmed to be a dwarf-dwarf merger. VCC848 is located in the outskirts of the Virgo Cluster. By analyzing the stellar light distribution, we found that VCC848 is the result of a merging between two dwarf galaxies with a primary-to-secondary mass ratio < ~ 5 for the stellar components and < ~ 2 for the presumed dark matter halos. The secondary progenitor galaxy has been almost entirely disrupted. The age-mass distribution of photometrically selected star cluster candidates in VCC848 implies that the cluster formation rate (CFR, proportional to star formation rate) was enhanced by a factor of ~ 7 - 10 during the past 1 Gyr. The merging-induced enhancement of CFR peaked near the galactic center a few hundred Myr ago and has started declining in the last few tens of Myr. The current star formation activities, as traced by the youngest clusters, mainly occur at large galactocentric distances (> ~ 1 kpc). The fact that VCC848 is still (atomic) gas-dominated after the period of most violent collision suggests that gas-rich dwarf galaxy merging can result in BCD-like remnants with extended atomic gas distribution surrounding a blue compact center, in general agreement with previous numerical simulations.


title: Memory-Based Graph Networks

abstract: Graph neural networks (GNNs) are a class of deep models that operate on data with arbitrary topology represented as graphs. We introduce an efficient memory layer for GNNs that can jointly learn node representations and coarsen the graph. We also introduce two new networks based on this layer: memory-based GNN (MemGNN) and graph memory network (GMN) that can learn hierarchical graph representations. The experimental results shows that the proposed models achieve state-of-the-art results in eight out of nine graph classification and regression benchmarks. We also show that the learned representations could correspond to chemical features in the molecule data. Code and reference implementations are released at: https://github.com/amirkhas/GraphMemoryNet


title: Exponential Automatic Amortized Resource Analysis

abstract: Automatic amortized resource analysis (AARA) is a type-based technique for inferring concrete (non-asymptotic) bounds on a program's resource usage. Existing work on AARA has focused on bounds that are polynomial in the sizes of the inputs. This paper presents and extension of AARA to exponential bounds that preserves the benefits of the technique, such as compositionality and efficient type inference based on linear constraint solving. A key idea is the use of the Stirling numbers of the second kind as the basis of potential functions, which play the same role as the binomial coefficients in polynomial AARA. To formalize the similarities with the existing analyses, the paper presents a general methodology for AARA that is instantiated to the polynomial version, the exponential version, and a combined system with potential functions that are formed by products of Stirling numbers and binomial coefficients. The soundness of exponential AARA is proved with respect to an operational cost semantics and the analysis of representative example programs demonstrates the effectiveness of the new analysis.


title: Proper actions of discrete groups of affine transformations

abstract: In the early 1980's Margulis startled the world by showing the existence of proper affine actions of free groups on 3-space, answering a provocative and suggestive question Milnor posed in 1977. In this paper we discuss the historical background motivating this question, recent progress on this subject, and future directions inspired by this discovery.


title: Improved bounds on sizes of generalized caps in $AG(n,q)$

abstract: An $m$-general set in $AG(n,q)$ is a set of points such that any subset of size $m$ is in general position. A $3$-general set is often called a capset. In this paper, we study the maximum size of an $m$-general set in $AG(n,q)$, significantly improving previous results. When $m=4$ and $q=2$ we give a precise estimate, solving a problem raised by Bennett.


title: Resonant Quantum Search with Monitor Qubits

abstract: We present an algorithm for the generalized search problem (searching $k$ marked items among $N$ items) based on a continuous Hamiltonian and exploiting resonance. This resonant algorithm has the same time complexity $O(\sqrt{N/k})$ as the Grover algorithm. A natural extension of the algorithm, incorporating auxiliary "monitor" qubits, can determine $k$ precisely, if it is unknown. The time complexity of our counting algorithm is $O(\sqrt{N})$, similar to the best quantum approximate counting algorithm, or better, given appropriate physical resources.


title: Struct-MMSB: Mixed Membership Stochastic Blockmodels with Interpretable Structured Priors

abstract: The mixed membership stochastic blockmodel (MMSB) is a popular framework for community detection and network generation. It learns a low-rank mixed membership representation for each node across communities by exploiting the underlying graph structure. MMSB assumes that the membership distributions of the nodes are independently drawn from a Dirichlet distribution, which limits its capability to model highly correlated graph structures that exist in real-world networks. In this paper, we present a flexible richly structured MMSB model, \textit{Struct-MMSB}, that uses a recently developed statistical relational learning model, hinge-loss Markov random fields (HL-MRFs), as a structured prior to model complex dependencies among node attributes, multi-relational links, and their relationship with mixed-membership distributions. Our model is specified using a probabilistic programming templating language that uses weighted first-order logic rules, which enhances the model's interpretability. Further, our model is capable of learning latent characteristics in real-world networks via meaningful latent variables encoded as a complex combination of observed features and membership distributions. We present an expectation-maximization based inference algorithm that learns latent variables and parameters iteratively, a scalable stochastic variation of the inference algorithm, and a method to learn the weights of HL-MRF structured priors. We evaluate our model on six datasets across three different types of networks and corresponding modeling scenarios and demonstrate that our models are able to achieve an improvement of 15\% on average in test log-likelihood and faster convergence when compared to state-of-the-art network models.


title: Efficient unitary designs with a system-size independent number of non-Clifford gates

abstract: Many quantum information protocols require the implementation of random unitaries. Because it takes exponential resources to produce Haar-random unitaries drawn from the full $n$-qubit group, one often resorts to $t$-designs. Unitary $t$-designs mimic the Haar-measure up to $t$-th moments. It is known that Clifford operations can implement at most $3$-designs. In this work, we quantify the non-Clifford resources required to break this barrier. We find that it suffices to inject $O(t^{4}\log^{2}(t)\log(1/\varepsilon))$ many non-Clifford gates into a polynomial-depth random Clifford circuit to obtain an $\varepsilon$-approximate $t$-design. Strikingly, the number of non-Clifford gates required is independent of the system size -- asymptotically, the density of non-Clifford gates is allowed to tend to zero. We also derive novel bounds on the convergence time of random Clifford circuits to the $t$-th moment of the uniform distribution on the Clifford group. Our proofs exploit a recently developed variant of Schur-Weyl duality for the Clifford group, as well as bounds on restricted spectral gaps of averaging operators.


title: On the refined Strichartz estimates

abstract: We present a slightly simpler proof of the multilinear refined Strichartz estimate, and prove a slightly more general linear refined Strichartz estimate. Our arguments seek to clarify the connection between these estimates, refined decoupling and tube incidences.


title: Stochastic Subspace Cubic Newton Method

abstract: In this paper, we propose a new randomized second-order optimization algorithm---Stochastic Subspace Cubic Newton (SSCN)---for minimizing a high dimensional convex function $f$. Our method can be seen both as a {\em stochastic} extension of the cubically-regularized Newton method of Nesterov and Polyak (2006), and a {\em second-order} enhancement of stochastic subspace descent of Kozak et al. (2019). We prove that as we vary the minibatch size, the global convergence rate of SSCN interpolates between the rate of stochastic coordinate descent (CD) and the rate of cubic regularized Newton, thus giving new insights into the connection between first and second-order methods. Remarkably, the local convergence rate of SSCN matches the rate of stochastic subspace descent applied to the problem of minimizing the quadratic function $\frac12 (x-x^*)^\top \nabla^2f(x^*)(x-x^*)$, where $x^*$ is the minimizer of $f$, and hence depends on the properties of $f$ at the optimum only. Our numerical experiments show that SSCN outperforms non-accelerated first-order CD algorithms while being competitive to their accelerated variants.


title: Spreading of correlations in Markovian open quantum systems

abstract: Understanding the spreading of quantum correlations in out-of-equilibrium many-body systems is one of the major challenges in physics. For {\it isolated} systems, a hydrodynamic theory explains the origin and spreading of entanglement via the propagation of quasi-particle pairs. However, when systems interact with their surrounding much less has been established. Here we show that the quasi-particle picture remains valid for open quantum systems: while information is still spread by quasiparticles, the environment modifies their correlation and introduces incoherent and mixing effects. For free fermions with gain/loss dissipation we provide formulae fully describing incoherent and quasiparticle contributions in the spreading of entropy and mutual information. Importantly, the latter is not affected by incoherent correlations. The mutual information is exponentially damped at short times and eventually vanishes signalling the onset of a classical limit. The behaviour of the logarithmic negativity is similar and this scenario is common to other dissipations. For weak dissipation, the presence of quasiparticles underlies remarkable scaling behaviors.

categories: cond-mat.stat-mech cond-mat.quant-gas cond-mat.str-el hep-th quant-ph


title: Measuring topological invariants in polaritonic graphene

abstract: Topological materials rely on engineering global properties of their bulk energy bands called topological invariants. These invariants, usually defined over the entire Brillouin zone, are related to the existence of protected edge states. However, for an important class of Hamiltonians corresponding to 2D lattices with time-reversal and chiral symmetry (e.g. graphene), the existence of edge states is linked to invariants that are not defined over the full 2D Brillouin zone, but on reduced 1D sub-spaces. Here, we demonstrate a novel scheme based on a combined real- and momentum-space measurement to directly access these 1D topological invariants in lattices of semiconductor microcavities confining exciton-polaritons. We extract these invariants in arrays emulating the physics of regular and critically compressed graphene sucht that Dirac cones have merged. Our scheme provides a direct evidence of the bulk-edge correspondence in these systems, and opens the door to the exploration of more complex topological effects, for example involving disorder and interactions.


title: The failure of the fractal uncertainty principle for the Walsh--Fourier transform

abstract: We construct $\delta$-regular sets with $\delta\ge \frac12$ for which the analog of the Bourgain--Dyatlov Fractal Uncertainty Principle fails for the Walsh--Fourier transform.


title: Using machine learning to separate hadronic and electromagnetic interactions in the GlueX forward calorimeter

abstract: The GlueX forward calorimeter is an array of 2800 lead glass modules that was constructed to detect photons produced in the decays of hadrons. A background to this process originates from hadronic interactions in the calorimeter, which, in some instances, can be difficult to distinguish from low energy photon interactions. Machine learning techniques were applied to the classification of particle interactions in the GlueX forward calorimeter. The algorithms were trained on data using decays of the $\omega$ meson, which contain both true photons and charged particles that interact with the calorimeter. Algorithms were evaluated on efficiency, rate of false positives, run time, and implementation complexity. An algorithm that utilizes a multi-layer perceptron neural net was deployed in the GlueX software stack and provides a signal efficiency of 85% with a background rejection of 60% for an inclusive $\pi^0$ data sample for an intermediate quality constraint.


title: Solitary wave solutions and global well-posedness for a coupled system of gKdV equations

abstract: In this work we consider the initial-value problem associated with a coupled system of generalized Korteweg-de Vries equations. We present a relationship between the best constant for a Gagliardo-Nirenberg type inequality and a criterion for the existence of global solutions in the energy space. We prove that such a constant is directly related to the existence problem of solitary-wave solutions with minimal mass, the so called ground state solutions. To guarantee the existence of ground states we use a variational method.


title: Deterministic correction of qubit loss

abstract: The loss of qubits - the elementary carriers of quantum information - poses one of the fundamental obstacles towards large-scale and fault-tolerant quantum information processors. In this work, we experimentally demonstrate a complete toolbox and the implementation of a full cycle of qubit loss detection and correction on a minimal instance of a topological surface code. This includes a quantum non-demolition measurement of a qubit loss event that conditionally triggers a restoration procedure, mapping the logical qubit onto a new encoding on the remaining qubits. The demonstrated methods, implemented here in a trapped-ion quantum processor, are applicable to other quantum computing architectures and codes, including leading 2D and 3D topological quantum error correcting codes. These tools complement previously demonstrated techniques to correct computational errors, and in combination constitute essential building blocks for complete and scalable quantum error correction.


title: Real-Time Visualization in Non-Isotropic Geometries

abstract: Non-isotropic geometries are of interest to low-dimensional topologists, physicists and cosmologists. However, they are challenging to comprehend and visualize. We present novel methods of computing real-time native geodesic rendering of non-isotropic geometries. Our methods can be applied not only to visualization, but also are essential for potential applications in machine learning and video games.


title: Hyperbolic Minesweeper is in P

abstract: We show that, while Minesweeper is NP-complete, its hyperbolic variant is in P. Our proof does not rely on the rules of Minesweeper, but is valid for any puzzle based on satisfying local constraints on a graph embedded in the hyperbolic plane.


title: RobustPeriod: Time-Frequency Mining for Robust Multiple Periodicity Detection

abstract: Periodicity detection is a crucial step in time series tasks, including monitoring and forecasting of metrics in many areas, such as IoT applications and self-driving database management system. In many of these applications, multiple periodic components exist and are often interlaced with each other. Such dynamic and complicated periodic patterns make the accurate periodicity detection difficult. In addition, other components in the time series, such as trend, outliers and noises, also pose additional challenges for accurate periodicity detection. In this paper, we propose a robust and general framework for multiple periodicity detection. Our algorithm applies maximal overlap discrete wavelet transform to transform the time series into multiple temporal-frequency scales such that different periodic components can be isolated. We rank them by wavelet variance, and then at each scale detect single periodicity by our proposed Huber-periodogram and Huber-ACF robustly. We rigorously prove the theoretical properties of Huber-periodogram and justify the use of Fisher's test on Huber-periodogram for periodicity detection. To further refine the detected periods, we compute unbiased autocorrelation function based on Wiener-Khinchin theorem from Huber-periodogram for improved robustness and efficiency. Experiments on synthetic and real-world datasets show that our algorithm outperforms other popular ones for both single and multiple periodicity detection.


title: Image to Language Understanding: Captioning approach

abstract: Extracting context from visual representations is of utmost importance in the advancement of Computer Science. Representation of such a format in Natural Language has a huge variety of applications such as helping the visually impaired etc. Such an approach is a combination of Computer Vision and Natural Language techniques which is a hard problem to solve. This project aims to compare different approaches for solving the image captioning problem. In specific, the focus was on comparing two different types of models: Encoder-Decoder approach and a Multi-model approach. In the encoder-decoder approach, inject and merge architectures were compared against a multi-modal image captioning approach based primarily on object detection. These approaches have been compared on the basis on state of the art sentence comparison metrics such as BLEU, GLEU, Meteor, and Rouge on a subset of the Google Conceptual captions dataset which contains 100k images. On the basis of this comparison, we observed that the best model was the Inception injected encoder model. This best approach has been deployed as a web-based system. On uploading an image, such a system will output the best caption associated with the image.


title: Dynamics of Magnus Dominated Particle Clusters, Collisions, Pinning and Ratchets

abstract: Motivated by the recent work in skyrmions and active chiral matter systems, we examine pairs and small clusters of repulsively interacting point particles in the limit where the dynamics is dominated by the Magnus force. We find that particles with the same Magnus force can form stable pairs, triples and higher ordered clusters or exhibit chaotic motion. For mixtures of particles with opposite Magnus force, particle pairs can combine to form translating dipoles. Under an applied drive, particles with the same Magnus force translate; however, particles with different or opposite Magnus force exhibit a drive-dependent decoupling transition. When the particles interact with a repulsive obstacle, they can form localized orbits with depinning or unwinding transitions under an applied drive. We examine the interaction of these particles with clusters or lines of obstacles, and find that the particles can become trapped in orbits that encircle multiple obstacles. Under an ac drive, we observe a series of ratchet effects, including ratchet reversals, for particles interacting with a line of obstacles due to the formation of commensurate orbits. Finally, in assemblies of particles with mixed Magnus forces of the same sign, we find that the particles with the largest Magnus force become localized in the center of the cluster, while for mixtures with opposite Magnus forces, the motion is dominated by transient local pairs or clusters, where the translating pairs can be regarded as a form of active matter.


title: Knot Selection in Sparse Gaussian Processes

abstract: Knot-based, sparse Gaussian processes have enjoyed considerable success as scalable approximations to full Gaussian processes. Problems can occur, however, when knot selection is done by optimizing the marginal likelihood. For example, the marginal likelihood surface is highly multimodal, which can cause suboptimal knot placement where some knots serve practically no function. This is especially a problem when many more knots are used than are necessary, resulting in extra computational cost for little to no gains in accuracy. We propose a one-at-a-time knot selection algorithm to select both the number and placement of knots. Our algorithm uses Bayesian optimization to efficiently propose knots that are likely to be good and largely avoids the pathologies encountered when using the marginal likelihood as the objective function. We provide empirical results showing improved accuracy and speed over the current standard approaches.


title: Overlap Local-SGD: An Algorithmic Approach to Hide Communication Delays in Distributed SGD

abstract: Distributed stochastic gradient descent (SGD) is essential for scaling the machine learning algorithms to a large number of computing nodes. However, the infrastructures variability such as high communication delay or random node slowdown greatly impedes the performance of distributed SGD algorithm, especially in a wireless system or sensor networks. In this paper, we propose an algorithmic approach named Overlap-Local-SGD (and its momentum variant) to overlap the communication and computation so as to speedup the distributed training procedure. The approach can help to mitigate the straggler effects as well. We achieve this by adding an anchor model on each node. After multiple local updates, locally trained models will be pulled back towards the synchronized anchor model rather than communicating with others. Experimental results of training a deep neural network on CIFAR-10 dataset demonstrate the effectiveness of Overlap-Local-SGD. We also provide a convergence guarantee for the proposed algorithm under non-convex objective functions.


title: Aphelion Cloud Belt Phase Function Investigations with Mars Color Imager (MARCI)

abstract: This paper constrains the scattering phase function of water ice clouds (WICs) found within Mars' Aphelion Cloud Belt (ACB), determined from orbit by processing publicly available raw Mars Color Imager (MARCI) data spanning solar longitudes (Ls) 42-170 during Mars Years (MYs) 28 and 29. MARCI visible wavelength data were calibrated and then pipeline-processed to select the pixels most likely to possess clouds. Mean phase function curves for the MARCI blue filter data were derived, and for all seasons investigated, modeled aggregates, plates, solid and hollow columns, bullet rosettes, and droxtals were all found to be plausible habits. Spheres were found to be the least plausible, but still possible. Additionally, this work probed the opposition surge to examine the slope of the linear relationship between column ice water content and cloud opacity on Mars, and found a significant dependence on particle radius. The half-width-half-maxima (HWHM) of the visible 180 degree peak of five MARCI images were found to agree better with modeled HWHM for WICs than with modeled HWHM for dust.


title: Evaluation of Automatic FPGA Offloading for Loop Statements of Applications

abstract: In recent years, with the prediction of Moore's law slowing down, utilization of hardware other than CPU such as FPGA which is energy effective is increasing. However, when using heterogeneous hardware other than CPUs, barriers of technical skills such as OpenCL are high. Based on that, I have proposed environment adaptive software that enables automatic conversion, configuration, and high-performance operation of once written code, according to the hardware to be placed. Partly of the offloading to the GPU was automated previously. In this paper, I propose and evaluate an automatic extraction method of appropriate offload target loop statements of source code as the first step of offloading to FPGA. I evaluate the effectiveness of the proposed method using an existing application.


title: Adaptation in general temporally changing environments

abstract: We analyze a nonlocal PDE model describing the dynamics of adaptation of a phenotypically structured population, under the effects of mutation and selection, in a changing environment. Previous studies have analyzed the large-time behavior of such models, with particular forms of environmental changes, either linearly changing or periodically fluctuating. We use here a completely different mathematical approach, which allows us to consider very general forms of environmental variations and to give an analytic description of the full trajectories of adaptation, including the transient phase, before a stationary behavior is reached. The main idea behind our approach is to study a bivariate distribution of two `fitness components' which contains enough information to describe the distribution of fitness at any time. This distribution solves a degenerate parabolic equation that is dealt with by defining a multidimensional cumulant generating function associated with the distribution, and solving the associated transport equation. We apply our results to several examples, and check their accuracy, using stochastic individual-based simulations as a benchmark. These examples illustrate the importance of being able to describe the transient dynamics of adaptation to understand the development of drug resistance in pathogens.


title: Modelling Latent Skills for Multitask Language Generation

abstract: We present a generative model for multitask conditional language generation. Our guiding hypothesis is that a shared set of latent skills underlies many disparate language generation tasks, and that explicitly modelling these skills in a task embedding space can help with both positive transfer across tasks and with efficient adaptation to new tasks. We instantiate this task embedding space as a latent variable in a latent variable sequence-to-sequence model. We evaluate this hypothesis by curating a series of monolingual text-to-text language generation datasets - covering a broad range of tasks and domains - and comparing the performance of models both in the multitask and few-shot regimes. We show that our latent task variable model outperforms other sequence-to-sequence baselines on average across tasks in the multitask setting. In the few-shot learning setting on an unseen test dataset (i.e., a new task), we demonstrate that model adaptation based on inference in the latent task space is more robust than standard fine-tuning based parameter adaptation and performs comparably in terms of overall performance. Finally, we examine the latent task representations learnt by our model and show that they cluster tasks in a natural way.


title: Quasinormal modes and thermalization in Improved Holographic QCD

abstract: Following a series of similar calculations in simpler non-conformal holographic setups, we determine the quasinormal mode spectrum for an operator dual to a gauge-invariant scalar field within the Improved Holographic QCD framework. At temperatures somewhat above the critical temperature of the deconfinement transition, we find a small number of clearly separated modes followed by a branch-cut-like structure parallel to the real axis, the presence of which is linked to the form of the IHQCD potential employed. The temperature dependence of the lowest nonzero mode is furthermore used to study the thermalization time of the corresponding correlator, which is found to be of the order of the inverse critical temperature near the phase transition and decrease slightly faster than $1/T$ at higher temperatures.


title: RobustTAD: Robust Time Series Anomaly Detection via Decomposition and Convolutional Neural Networks

abstract: The monitoring and management of numerous and diverse time series data at Alibaba Group calls for an effective and scalable time series anomaly detection service. In this paper, we propose RobustTAD, a Robust Time series Anomaly Detection framework by integrating robust seasonal-trend decomposition and convolutional neural network for time series data. The seasonal-trend decomposition can effectively handle complicated patterns in time series, and meanwhile significantly simplifies the architecture of the neural network, which is an encoder-decoder architecture with skip connections. This architecture can effectively capture the multi-scale information from time series, which is very useful in anomaly detection. Due to the limited labeled data in time series anomaly detection, we systematically investigate data augmentation methods in both time and frequency domains. We also introduce label-based weight and value-based weight in the loss function by utilizing the unbalanced nature of the time series anomaly detection problem. Compared with the widely used forecasting-based anomaly detection algorithms, decomposition-based algorithms, traditional statistical algorithms, as well as recent neural network based algorithms, RobustTAD performs significantly better on public benchmark datasets. It is deployed as a public online service and widely adopted in different business scenarios at Alibaba Group.


title: IMDfence: Architecting a Secure Protocol for Implantable Medical Devices

abstract: Over the past decade, focus on the security and privacy aspects of implantable medical devices (IMDs) has intensified, driven by the multitude of cybersecurity vulnerabilities found in various existing devices. However, due to their strict computational, energy and physical constraints, conventional security protocols are not directly applicable to IMDs. Custom-tailored schemes have been proposed instead which, however, fail to cover the full spectrum of security features that modern IMDs and their ecosystems so critically require. In this paper we propose IMDfence, a security protocol for IMD ecosystems that provides a comprehensive yet practical security portfolio, which includes availability, non-repudiation, access control, entity authentication, remote monitoring and system scalability. The protocol also allows emergency access that results in the graceful degradation of offered services without compromising security and patient safety. The performance of the security protocol as well as its feasibility and impact on modern IMDs are extensively analyzed and evaluated. We find that IMDfence achieves the above security requirements at a mere less than 7% increase in total IMD energy consumption, and less than 14 ms and 9 kB increase in system delay and memory footprint, respectively.


title: Stochastic Normalizing Flows

abstract: We introduce stochastic normalizing flows, an extension of continuous normalizing flows for maximum likelihood estimation and variational inference (VI) using stochastic differential equations (SDEs). Using the theory of rough paths, the underlying Brownian motion is treated as a latent variable and approximated, enabling efficient training of neural SDEs as random neural ordinary differential equations. These SDEs can be used for constructing efficient Markov chains to sample from the underlying distribution of a given dataset. Furthermore, by considering families of targeted SDEs with prescribed stationary distribution, we can apply VI to the optimization of hyperparameters in stochastic MCMC.


title: Three-Nucleon Forces: Implementation and Applications to Atomic Nuclei and Dense Matter

abstract: Recent advances in nuclear structure theory have significantly enlarged the accessible part of the nuclear landscape via ab initio many-body calculations. These developments open new ways for microscopic studies of light, medium-mass and heavy nuclei as well as nuclear matter and represent an important step toward a systematic and comprehensive understanding of atomic nuclei across the nuclear chart. While remarkable agreement has been found between different many-body methods for a given nuclear Hamiltonian, the comparison with experiment and the understanding of theoretical uncertainties are still important open questions. The observed discrepancies to experiment indicate deficiencies in presently used nuclear interactions and operators. Chiral effective field theory (EFT) allows to systematically derive contributions to nucleon-nucleon (NN), three-nucleon (3N) and higher-body interactions including estimates of theoretical uncertainties. While the treatment of NN interactions in many-body calculations is well established, the calculation of 3N interactions and their incorporation in ab initio frameworks is still a frontier. This work reviews in detail recent and current developments on the derivation and implementation of improved 3N interactions and provides a comprehensive introduction to fundamental methods for their practical calculation and representation. We further give an overview of novel and established methods that facilitate the inclusion and treatment of 3N interactions in ab initio nuclear structure frameworks and present a selection of the latest calculations of atomic nuclei as well as nuclear matter based on state-of-the-art nuclear NN and 3N interactions derived within chiral EFT. Finally, we discuss ongoing efforts, open questions and future directions.


title: Optimal Signal-Adaptive Trading with Temporary and Transient Price Impact

abstract: We study optimal liquidation in the presence of linear temporary and transient price impact along with taking into account a general price predicting finite-variation signal. We formulate this problem as minimization of a cost-risk functional over a class of absolutely continuous and signal-adaptive strategies. The stochastic control problem is solved by following a probabilistic and convex analytic approach. We show that the optimal trading strategy is given by a system of four coupled forward-backward SDEs, which can be solved explicitly. Our results reveal how the induced transient price distortion provides together with the predictive signal an additional predictor about future price changes. As a consequence, the optimal signal-adaptive trading rate trades off exploiting the predictive signal against incurring the transient displacement of the execution price from its unaffected level. This answers an open question from Lehalle and Neuman [29] as we show how to derive the unique optimal signal-adaptive liquidation strategy when price impact is not only temporary but also transient.


title: Iron Line Tomography of General Relativistic Hydrodynamic Accretion around Kerr Black Holes

abstract: We consider a temporal response of relativistically broadened line spectrum of iron from black hole accretion irradiated by an X-ray echo under strong gravity. The physical condition of accreting gas is numerically calculated in the context of general relativistic hydrodynamics under steady-state, axisymmetry in Kerr geometry. With the onset of a point-like X-ray flare of a short finite duration just above the accretion surface, the gas is assumed to be ionized to produce a neutral iron fluorescent line. Using a fully relativistic ray-tracing approach, the response of line photons due to the X-ray illumination is traced as a function of time and energy for different source configurations around \sw and Kerr black holes. Our calculations show that the X-ray echo on the accretion surface clearly imprints a characteristic time-variability in the line spectral features depending on those parameters. Simulated line profiles, aimed for the future microcalorimeter missions of large collecting area such as {\it Athena}/X-IFU for typical radio-quiet Seyfert galaxies, are presented to demonstrate that state-of-the-art new observations could differentiate various source parameters by such an X-ray tomographic line reverberation.


title: On the Stable Transfer for $\mathrm{Sym}^{n}$ Lifting of $\mathrm{GL}_{2}$

abstract: Following the paradigm of \cite{MR3117742}, we are going to explore the stable transfer factors for $\mathrm{Sym}^{n}$ lifting from $\mathrm{GL}_{2}$ to $\mathrm{GL}_{n+1}$ over any local fields $F$ of characteristic zero with residue characteristic not equal to $2$. When $F=\mathbb{C}$ we construct an explicit stable transfer factor for any $n$. When $n$ is odd, we provide a reduction formula, reducing the question to the construction of the stable transfer factors when the $L$-morphism is the diagonal embedding from $\mathrm{GL}_{2}(\mathbb{C})$ to finitely many copies of $\mathrm{GL}_{2}(\mathbb{C})$ under mild assumptions on the residue characteristic of $F$. With the assumptions on the residue characteristic, the reduction formula works uniformly over any local fields of characteristic zero, except that for $p$-adic situation we need to exclude the twisted Steinberg representations.


title: Visible and near-infrared spectro-interferometric analysis of the edge-on Be star $\omicron$ Aquarii

abstract: We present a detailed visible and near-IR spectro-interferometric analysis of the Be-shell star $\omicron$ Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. We measured the stellar radius of $\omicron$ Aquarii as 4.0 $\pm$ 0.3 $\mathrm{R_{\odot}}$. We constrained the disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The disk sizes in H$\alpha$ and Br$\gamma$ were found to be similar, at $\sim$10-12 $\mathrm{D_{\star}}$, which is uncommon since most results for Be stars show a larger extension in H$\alpha$ than in Br$\gamma$. We found that the inclination angle $i$ derived from H$\alpha$ is significantly lower ($\sim$15 deg) than the one derived from Br$\gamma$. The disk kinematics were found to be near to the Keplerian rotation in Br$\gamma$, but not in H$\alpha$. After analyzing all our data using a grid of HDUST models (BeAtlas), we found a common physical description for the disk in both lines: $\Sigma_{0}$ = 0.12 g cm\textsuperscript{-2} and $m$ = 3.0. The stellar rotational rate was found to be very close ($\sim$96\%) to the critical value. Our analysis of multi-epoch H$\alpha$ profiles and imaging polarimetry indicates that the disk has been stable for at least 20 years. Compared to Br$\gamma$, the data in H$\alpha$ shows a substantially different picture that cannot fully be understood using the current physical models of Be star disks. $\omicron$ Aquarii presents a stable disk, but the measured $m$ is lower than the standard value in the VDD model for steady-state. Such long-term stability can be understood in terms of the high rotational rate for this star, the rate being a main source for the mass injection in the disk. Our results on the stellar rotation and disk stability are consistent with results in the literature showing that late-type Be stars are more likely to be fast rotators and have stable disks.


title: Sequential decomposition of discrete memoryless channel with noisy feedback

abstract: In this paper, we consider a discrete memoryless point to point channel with noisy feedback, where there is a sender with a private message that she wants to communicate to a receiver by sequentially transmitting symbols over a noisy channel. After each transmission, she receives a noisy feedback of the symbol received by the receiver. The goal is to design transmission control strategy of the sender that minimize the average probability of error. This is an instance of decentralized control of information where the two controllers, the sender and the receiver have no common information. There exist no methodology in the literature that provides a notion of "state" and a dynamic program to find optimal policies for this problem In this paper, we show introduce a notion of state, based on which we provide a sequential decomposition methodology that finds optimum policies within the class of Markov strategies with respect to this state (which need not be globally optimum). This allows to decompose the problem across time and reduce the complexity dependence on time from double exponential to linear in time.


title: Particle Filter Based Monocular Human Tracking with a 3D Cardbox Model and a Novel Deterministic Resampling Strategy

abstract: The challenge of markerless human motion tracking is the high dimensionality of the search space. Thus, efficient exploration in the search space is of great significance. In this paper, a motion capturing algorithm is proposed for upper body motion tracking. The proposed system tracks human motion based on monocular silhouette-matching, and it is built on the top of a hierarchical particle filter, within which a novel deterministic resampling strategy (DRS) is applied. The proposed system is evaluated quantitatively with the ground truth data measured by an inertial sensor system. In addition, we compare the DRS with the stratified resampling strategy (SRS). It is shown in experiments that DRS outperforms SRS with the same amount of particles. Moreover, a new 3D articulated human upper body model with the name 3D cardbox model is created and is proven to work successfully for motion tracking. Experiments show that the proposed system can robustly track upper body motion without self-occlusion. Motions towards the camera can also be well tracked.


title: Invariant measures and global well posedness for SQG equation

abstract: We construct an invariant measure $\mu$ for the Surface Quasi-Geostrophic (SQG) equation and show that almost all functions in the support of $\mu$ are initial conditions of global, unique solutions of SQG, that depend continuously on the initial data. In addition, we show that the support of $\mu$ is infinite dimensional, meaning that it is not locally a subset of any compact set with finite Hausdorff dimension. Also, there are global solutions that have arbitrarily large initial condition. The measures $\mu$ is obtained via fluctuation-dissipation method, that is, as a limit of invariant measures for stochastic SQG with a carefully chosen dissipation and random forcing.


title: CALICE SiW ECAL -- Development and performance of a highly compact digital readout system

abstract: A highly granular silicon-tungsten electromagnetic calorimeter (SiW-ECAL) is the reference design of the ECAL of the International Large Detector (ILD) concept, one of the two detector concepts for the detector(s) at the future International Linear Collider. Prototypes for this type of detector are developed within the CALICE Collaboration. During the last year a highly compact digital readout system has been built. The system has been used for the first time in a beam test in Summer 2019 at DESY. This article summarises the main features of the system and report on its performance during the beam test.


title: Thermalization and entropy production of a 1D bipartite fermionic lattice under influence of a dephasing noise

abstract: In this paper we study entropy production and transport of heat and matter in an 1D bipartite fermionic chain. Moreover, we input a phenomenological dephasing noise in order to show its effect on the dynamics of the system. Specially, it is shown the importance of such noise in the emergence of the entropy production. For the fluxes, analytical solutions are obtained using both Boltzmann and Fermi-Dirac distribuition. Finally, we show that a particular fluctuation theorem for energy and matter exchange is obeyed even in the presence of such noise.


title: Self-Supervised Poisson-Gaussian Denoising

abstract: We extend the blindspot model for self-supervised denoising to handle Poisson-Gaussian noise and introduce an improved training scheme that avoids hyperparameters and adapts the denoiser to the test data. Self-supervised models for denoising learn to denoise from only noisy data and do not require corresponding clean images, which are difficult or impossible to acquire in some application areas of interest such as low-light microscopy. We introduce a new training strategy to handle Poisson-Gaussian noise which is the standard noise model for microscope images. Our new strategy eliminates hyperparameters from the loss function, which is important in a self-supervised regime where no ground truth data is available to guide hyperparameter tuning. We show how our denoiser can be adapted to the test data to improve performance. Our evaluations on microscope image denoising benchmarks validate our approach.


title: Temperature dependence of bend elastic constant in oblique helicoidal cholesterics

abstract: Elastic moduli of liquid crystals, known as Frank constants, are of quintessential importance for understanding fundamental properties of these materials and for the design of their applications. Although there are many methods to measure the Frank constants in the nematic phase, little is known about the elastic constants of the chiral version of the nematic, the so-called cholesteric liquid crystal, since the helicoidal structure of the cholesteric renders these methods inadequate. Here we present a technique to measure the bend modulus $K_{33}$ of cholesterics that is based on the electrically tunable reflection of light at an oblique helicoidal $Ch_{OH}$ cholesteric structure. $K_{33}$ is typically smaller than 0.6 pN, showing a non-monotonous temperature dependence with a slight increase near the transition to the twist-bend phase. $K_{33}$ depends strongly on the molecular composition. In particular, chiral mixtures that contain the flexible dimer 1'',7''-bis(4-cyanobiphenyl-4'-yl) heptane (CB7CB) and rod-like molecules such as pentylcyanobiphenyl (5CB) show a $K_{33}$ value that is 5 times smaller than $K_{33}$ of pure CB7CB or of mixtures of CB7CB with chiral dopants. Furthermore, $K_{33}$ in CB11CB doped with a chiral agent is noticeably smaller than $K_{33}$ in a similarly doped CB7CB which is explained by the longer flexible link in CB11CB. The proposed technique allows a direct in-situ determination of how the molecular composition, molecular structure and molecular chirality affect the elastic properties of chiral liquid crystals.


title: Practical Verification of MapReduce Computation Integrity via Partial Re-execution

abstract: Big data processing is often outsourced to powerful, but untrusted cloud service providers that provide agile and scalable computing resources to weaker clients. However, untrusted cloud services do not ensure the integrity of data and computations while clients have no control over the outsourced computation or no means to check the correctness of the execution. Despite a growing interest and recent progress in verifiable computation, the existing techniques are still not practical enough for big data processing due to high verification overhead. In this paper, we present a solution called V-MR (Verifiable MapReduce), which is a framework that verifies the integrity of MapReduce computation outsourced in the untrusted cloud via partial re-execution. V-MR is practically effective and efficient in that (1) it can detect the violation of MapReduce computation integrity and identify the malicious workers involved in the that produced the incorrect computation. (2) it can reduce the overhead of verification via partial re-execution with carefully selected input data and program code using program analysis. The experiment results of a prototype of V-MR show that V-MR can verify the integrity of MapReduce computation effectively with small overhead for partial re-execution.