title: A Spectral Method for Stable Bispectrum Inversion with Application to Multireference Alignment

abstract: We focus on an alignment-free method to estimate the underlying signal from a large number of noisy randomly shifted observations. Specifically, we estimate the mean, power spectrum, and bispectrum of the signal from the observations. Since bispectrum contains the phase information of the signal, reliable algorithms for bispectrum inversion is useful in many applications. We propose a new algorithm using spectral decomposition of the normalized bispectrum matrix for this task. For clean signals, we show that the eigenvectors of the normalized bispectrum matrix correspond to the true phases of the signal and its shifted copies. In addition, the spectral method is robust to noise. It can be used as a stable and efficient initialization technique for local non-convex optimization for bispectrum inversion.


title: Global-in-time Stability of 2D MHD boundary Layer in the Prandtl-Hartmann Regime

abstract: In this paper, we prove global existence of solutions with analytic regularity to the 2D MHD boundary layer equations in the mixed Prandtl and Hartmann regime derived by formal multi-scale expansion in \cite{GP}. The analysis shows that the combined effect of the magnetic diffusivity and transveral magnetic field on the boundary leads to a linear damping on the tangential velocity field near the boundary. And this damping effect yields the global in time analytic norm estimate in the tangential space variable on the perturbation of the classical steady Hartmann profile.


title: Pop Music Highlighter: Marking the Emotion Keypoints

abstract: The goal of music highlight extraction is to get a short consecutive segment of a piece of music that provides an effective representation of the whole piece. In a previous work, we introduced an attention-based convolutional recurrent neural network that uses music emotion classification as a surrogate task for music highlight extraction, for Pop songs. The rationale behind that approach is that the highlight of a song is usually the most emotional part. This paper extends our previous work in the following two aspects. First, methodology-wise we experiment with a new architecture that does not need any recurrent layers, making the training process faster. Moreover, we compare a late-fusion variant and an early-fusion variant to study which one better exploits the attention mechanism. Second, we conduct and report an extensive set of experiments comparing the proposed attention-based methods against a heuristic energy-based method, a structural repetition-based method, and a few other simple feature-based methods for this task. Due to the lack of public-domain labeled data for highlight extraction, following our previous work we use the RWC POP 100-song data set to evaluate how the detected highlights overlap with any chorus sections of the songs. The experiments demonstrate the effectiveness of our methods over competing methods. For reproducibility, we open source the code and pre-trained model at https://github.com/remyhuang/pop-music-highlighter/.


title: Epidemiologic analyses with error-prone exposures: Review of current practice and recommendations

abstract: Background: Variables in epidemiological observational studies are commonly subject to measurement error and misclassification, but the impact of such errors is frequently not appreciated or ignored. As part of the STRengthening Analytical Thinking for Observational Studies (STRATOS) Initiative, a Task Group on measurement error and misclassification (TG4) seeks to describe the scope of this problem and the analysis methods currently in use to address measurement error. Methods: TG4 conducted a literature survey of four types of research studies that are typically impacted by exposure measurement error: 1) dietary intake cohort studies, 2) dietary intake population surveys, 3) physical activity cohort studies, and 4) air pollution cohort studies. The survey was conducted to understand current practice for acknowledging and addressing measurement error. Results: The survey revealed that while researchers were generally aware that measurement error affected their studies, very few adjusted their analysis for the error. Most articles provided incomplete discussion of the potential effects of measurement error on their results. Regression calibration was the most widely used method of adjustment. Conclusions: Even in areas of epidemiology where measurement error is a known problem, the dominant current practice is to ignore errors in analyses. Methods to correct for measurement error are available but require additional data to inform the error structure. There is a great need to incorporate such data collection within study designs and improve the analytical approach. Increased efforts by investigators, editors and reviewers are also needed to improve presentation of research when data are subject to error.


title: Learning Discriminative Multilevel Structured Dictionaries for Supervised Image Classification

abstract: Sparse representations using overcomplete dictionaries have proved to be a powerful tool in many signal processing applications such as denoising, super-resolution, inpainting, compression or classification. The sparsity of the representation very much depends on how well the dictionary is adapted to the data at hand. In this paper, we propose a method for learning structured multilevel dictionaries with discriminative constraints to make them well suited for the supervised pixelwise classification of images. A multilevel tree-structured discriminative dictionary is learnt for each class, with a learning objective concerning the reconstruction errors of the image patches around the pixels over each class-representative dictionary. After the initial assignment of the class labels to image pixels based on their sparse representations over the learnt dictionaries, the final classification is achieved by smoothing the label image with a graph cut method and an erosion method. Applied to a common set of texture images, our supervised classification method shows competitive results with the state of the art.


title: Uncovering Multiscale Order in the Prime Numbers via Scattering

abstract: The prime numbers have been a source of fascination for millenia and continue to surprise us. Motivated by the hyperuniformity concept, which has attracted recent attention in physics and materials science, we show that the prime numbers in certain large intervals possess unanticipated order across length scales and represent the first example of a new class of many-particle systems with pure point diffraction patterns, which we call {\it effectively limit-periodic}. In particular, the primes in this regime are hyperuniform. This is shown analytically using the structure factor $S(k)$, proportional to the scattering intensity from a many-particle system. Remarkably, the structure factor for primes is characterized by dense Bragg peaks, like a quasicrystal, but positioned at certain rational wavenumbers, like a limit-periodic point pattern. We identify a transition between ordered and disordered prime regimes that depends on the intervals studied. Our analysis leads to an algorithm that enables one to predict primes with high accuracy. Effective limit-periodicity deserves future investigation in physics, independent of its link to the primes.


title: The initial-boundary value problem for the biharmonic Schr\"odinger equation on the half-line

abstract: We study the local and global wellposedness of the initial-boundary value problem for the biharmonic Schr\"odinger equation on the half-line with inhomogeneous Dirichlet-Neumann boundary data. First, we obtain a representation formula for the solution of the linear nonhomogenenous problem by using the Fokas method (also known as the \emph{unified transform method}). We use this representation formula to prove space and time estimates on the solutions of the linear model in fractional Sobolev spaces by using Fourier analysis. Secondly, we consider the nonlinear model with a power type nonlinearity and prove the local wellposedness by means of a classical contraction argument. We obtain Strichartz estimates to treat the low regularity case by using the oscillatory integral theory directly on the representation formula provided by the Fokas method. Global wellposedness of the defocusing model is established up to cubic nonlinearities by using the multiplier technique and proving hidden trace regularities.


title: Electroweak phase transition in the $\Sigma$SM - I: Dimensional reduction

abstract: In a series of two papers, we make a comparative analysis of the performance of conventional perturbation theory to analyze electroweak phase transition in the real triplet extension of Standard Model ($\Sigma$SM). In Part I (this paper), we derive and present the high-$T$ dimensionally reduced effective theory that is suitable for numerical simulation on the lattice. In the sequel (Part II), we will present results of the numerical simulation and benchmark the performance of conventional perturbation theory. Under the assumption that $\Sigma$ is heavy, the resulting effective theory takes the same form as that derived from the minimal standard model. By recasting the existing non-perturbative results, we map out the phase diagram of the model in the plane of triplet mass $M_\Sigma$ and Higgs portal coupling $a_2$. Contrary to conventional perturbation theory, we find regions of parameter space where the phase transition may be first order, second order, or crossover. We comment on prospects for prospective future colliders to probe the region where the electroweak phase transition is first order by a precise measurement of the $h\rightarrow\gamma\gamma$ partial width.


title: Predictive Uncertainty Estimation via Prior Networks

abstract: Estimating how uncertain an AI system is in its predictions is important to improve the safety of such systems. Uncertainty in predictive can result from uncertainty in model parameters, irreducible data uncertainty and uncertainty due to distributional mismatch between the test and training data distributions. Different actions might be taken depending on the source of the uncertainty so it is important to be able to distinguish between them. Recently, baseline tasks and metrics have been defined and several practical methods to estimate uncertainty developed. These methods, however, attempt to model uncertainty due to distributional mismatch either implicitly through model uncertainty or as data uncertainty. This work proposes a new framework for modeling predictive uncertainty called Prior Networks (PNs) which explicitly models distributional uncertainty. PNs do this by parameterizing a prior distribution over predictive distributions. This work focuses on uncertainty for classification and evaluates PNs on the tasks of identifying out-of-distribution (OOD) samples and detecting misclassification on the MNIST dataset, where they are found to outperform previous methods. Experiments on synthetic and MNIST and CIFAR-10 data show that unlike previous non-Bayesian methods PNs are able to distinguish between data and distributional uncertainty.


title: Coefficient systems on the Bruhat-Tits building and pro-$p$ Iwahori-Hecke modules

abstract: Let $G$ be the group of rational points of a split connected reductive group over a nonarchimedean local field of residue characteristic $p$. Let $I$ be a pro-$p$ Iwahori subgroup of $G$ and let $R$ be a commutative quasi-Frobenius ring. If $H=R[I\backslash G/I]$ denotes the pro-$p$ Iwahori-Hecke algebra of $G$ over $R$ we clarify the relation between the category of $H$-modules and the category of $G$-equivariant coefficient systems on the semisimple Bruhat-Tits building of $G$. If $R$ is a field of characteristic zero this yields alternative proofs of the exactness of the Schneider-Stuhler resolution and of the Zelevinski conjecture for smooth $G$-representations generated by their $I$-invariants. In general, it gives a description of the derived category of $H$-modules in terms of smooth $G$-representations and yields a functor to generalized $(\varphi,\Gamma)$-modules extending the constructions of Colmez, Schneider and Vign\'eras.


title: Anticipation in Human-Robot Cooperation: A Recurrent Neural Network Approach for Multiple Action Sequences Prediction

abstract: Close human-robot cooperation is a key enabler for new developments in advanced manufacturing and assistive applications. Close cooperation require robots that can predict human actions and intent, and understand human non-verbal cues. Recent approaches based on neural networks have led to encouraging results in the human action prediction problem both in continuous and discrete spaces. Our approach extends the research in this direction. Our contributions are three-fold. First, we validate the use of gaze and body pose cues as a means of predicting human action through a feature selection method. Next, we address two shortcomings of existing literature: predicting multiple and variable-length action sequences. This is achieved by introducing an encoder-decoder recurrent neural network topology in the discrete action prediction problem. In addition, we theoretically demonstrate the importance of predicting multiple action sequences as a means of estimating the stochastic reward in a human robot cooperation scenario. Finally, we show the ability to effectively train the prediction model on a action prediction dataset, involving human motion data, and explore the influence of the model's parameters on its performance. Source code repository: https://github.com/pschydlo/ActionAnticipation


title: An abelian subfield of the dyadic division field of a hyperelliptic Jacobian

abstract: Given a field $k$ of characteristic different from $2$ and an integer $d \geq 3$, let $J$ be the Jacobian of the "generic" hyperelliptic curve given by $y^2 = \prod_{i = 1}^d (x - \alpha_i)$, where the $\alpha_i$'s are transcendental and independent over $k$; it is defined over the transcendental extension $K / k$ generated by the symmetric functions of the $\alpha_i$'s. We investigate certain subfields of the field $K_{\infty}$ obtained by adjoining all points of $2$-power order of $J(\bar{K})$. In particular, we explicitly describe the maximal abelian subextension of $K_{\infty} / K(J[2])$ and show that it is contained in $K(J[8])$ (resp. $K(J[16])$) if $g \geq 2$ (resp. if $g = 1$). On the way we obtain an explicit description of the abelian subextension $K(J[4])$, and we describe the action of a particular automorphism in $\mathrm{Gal}(K_{\infty} / K)$ on these subfields.


title: Vibrational anomalies in AFe$_\mathbf{2}$As$_\mathbf{2}$ (A$\,=\,$Ca, Sr, and Ba) single crystals

abstract: The detailed behavior of the in-plane infrared-active vibrational modes has been determined in AFe$_2$As$_2$ (A$\,=\,$Ca, Sr, and Ba) above and below the structural and magnetic transition at $T_N=$172, 195 and 138 K, respectively. Above $T_N$, two infrared-active $E_u$ modes are observed. In all three compounds, below $T_N$ the low-frequency $E_u$ mode is observed to split into upper and lower branches; with the exception of the Ba material, the oscillator strength across the transition is conserved. In the Ca and Sr materials, the high-frequency $E_u$ mode splits into an upper and a lower branch; however, the oscillator strengths are quite different. Surprisingly, in both the Sr and Ba materials, below $T_N$ the upper branch appears be either very weak or totally absent, while the lower branch displays an anomalous increase in strength. The frequencies and atomic characters of the lattice modes at the center of the Brillouin zone have been calculated for the high-temperature phase for each of these materials. The high-frequency $E_u$ mode does not change in position or character across this series of compounds. Below $T_N$, the $E_u$ modes are predicted to split into features of roughly equal strength. We discuss the possibility that the anomalous increase in the strength of the lower branch of the high-frequency mode below $T_N$ in the Sr and Ba compounds, and the weak (silent) upper branch, may be related to the orbital ordering and a change in the bonding between the Fe and As atoms in the magnetically-ordered state.


title: Behavior of thin disk crystalline morphology in the presence of corrections to ideal magnetohydrodynamics

abstract: We analyze an axisymmetric magnetohydrodynamics configuration, describing the morphology of a purely differentially rotating thin plasma disk, in which linear and non-linear perturbations are triggered associated with microscopic magnetic structures. We study the evolution of the non-stationary correction in the limit in which the co-rotation condition (i.e., the dependence of the disk angular velocity on the magnetic flux function) is preserved and the poloidal velocity components are neglected. The main feature we address here is the influence of ideal (finite electron inertia) and collisional (resistivity, viscosity, and thermal conductivity) effects on the behavior of the flux function perturbation and of the associated small-scale modifications in the disk. We analyze two different regimes in which resistivity or viscosity dominates and study the corresponding linear and non-linear behaviors of the perturbation evolution, i.e., when the backreaction magnetic field is negligible or comparable to the background one, respectively. We demonstrate that when resistivity dominates, a radial oscillating morphology (crystalline structure) emerges and it turns out to be damped in time, in both the linear and non-linear regimes, but in such a way that the resulting transient can be implemented in the description of relevant astrophysical processes, for instance, associated with jet formation or cataclysmic variables. When the viscosity effect dominates the dynamics, only the non-linear regime is available and a very fast instability is triggered.


title: Unveiling the superconducting mechanism of Ba$_{0.51}$K$_{0.49}$BiO$_3$

abstract: Bismuthates were the first family of oxide high-temperature superconductors, exhibiting superconducting transition temperatures (Tc) up to 32K, but the superconducting mechanism remains under debate despite more than 30 years of extensive research. Our angle-resolved photoemission spectroscopy studies on Ba$_{0.51}$K$_{0.49}$BiO$_3$ reveal an unexpectedly 34% larger bandwidth than in conventional density functional theory calculations. This can be reproduced by calculations that fully account for long-range Coulomb interactions --- the first direct demonstration of bandwidth expansion due to the Fock exchange term, a long-accepted and yet uncorroborated fundamental effect in many body physics. Furthermore, we observe an isotropic superconducting gap with 2\Delta$_0$/k$_B$ T$_c$ = 3.51 $\pm$ 0.05, and strong electron-phonon interactions with a coupling constant \lambda$\sim$ 1.3 $\pm$ 0.2. These findings solve a long-standing mystery --- Ba$_{0.51}$K$_{0.49}$BiO$_3$ is an extraordinary Bardeen-Cooper-Schrieffer (BCS) superconductor, where long-range Coulomb interactions expand the bandwidth, enhance electron-phonon coupling, and generate the high Tc. Such effects will also be critical for finding new superconductors.


title: Brain Tumor Segmentation and Radiomics Survival Prediction: Contribution to the BRATS 2017 Challenge

abstract: Quantitative analysis of brain tumors is critical for clinical decision making. While manual segmentation is tedious, time consuming and subjective, this task is at the same time very challenging to solve for automatic segmentation methods. In this paper we present our most recent effort on developing a robust segmentation algorithm in the form of a convolutional neural network. Our network architecture was inspired by the popular U-Net and has been carefully modified to maximize brain tumor segmentation performance. We use a dice loss function to cope with class imbalances and use extensive data augmentation to successfully prevent overfitting. Our method beats the current state of the art on BraTS 2015, is one of the leading methods on the BraTS 2017 validation set (dice scores of 0.896, 0.797 and 0.732 for whole tumor, tumor core and enhancing tumor, respectively) and achieves very good Dice scores on the test set (0.858 for whole, 0.775 for core and 0.647 for enhancing tumor). We furthermore take part in the survival prediction subchallenge by training an ensemble of a random forest regressor and multilayer perceptrons on shape features describing the tumor subregions. Our approach achieves 52.6% accuracy, a Spearman correlation coefficient of 0.496 and a mean square error of 209607 on the test set.


title: Are there monopoles in the quark-gluon plasma?

abstract: Monopole-like objects have been identified in multiple lattice studies, and there is now a significant amount of literature on their importance in phenomenology. Some analytic indications of their role, however, are still missing. The 't Hooft-Polyakov monopoles, originally derived in the Georgi-Glashow model, are an important dynamical ingredient in theories with extended supersymmetry ${\cal N} = 2,\,4$, and help explain the issues related with electric-magnetic duality. There is no such solution in QCD-like theories without scalar fields. However, all of these theories have instantons and their finite-$T$ constituents known as instanton-dyons (or instanton-monopoles). The latter leads to semiclassical partition functions, which for ${\cal N} = 2,\,4$ theories were shown to be identical ("Poisson dual") to the partition function for monopoles. We show how, in a pure gauge theory, the semiclassical instanton-based partition function can also be Poisson-transformed into a partition function, interpreted as the one of moving and rotating monopoles.


title: Automated design of collective variables using supervised machine learning

abstract: Selection of appropriate collective variables for enhancing sampling of molecular simulations remains an unsolved problem in computational biophysics. In particular, picking initial collective variables (CVs) is particularly challenging in higher dimensions. Which atomic coordinates or transforms there of from a list of thousands should one pick for enhanced sampling runs? How does a modeler even begin to pick starting coordinates for investigation? This remains true even in the case of simple two state systems and only increases in difficulty for multi-state systems. In this work, we solve the initial CV problem using a data-driven approach inspired by the filed of supervised machine learning. In particular, we show how the decision functions in supervised machine learning (SML) algorithms can be used as initial CVs (SML_cv) for accelerated sampling. Using solvated alanine dipeptide and Chignolin mini-protein as our test cases, we illustrate how the distance to the Support Vector Machines' decision hyperplane, the output probability estimates from Logistic Regression, the outputs from deep neural network classifiers, and other classifiers may be used to reversibly sample slow structural transitions. We discuss the utility of other SML algorithms that might be useful for identifying CVs for accelerating molecular simulations.


title: Sidon set systems

abstract: A family ${\mathcal A}$ of $k$-subsets of $\{1,2,\dots, N\}$ is a Sidon system if the sumsets $A+B$, $A,B\in \mathcal{A}$ are pairwise distinct. We show that the largest cardinality $F_k(N)$ of a Sidon system of $k$-subsets of $[N]$ satisfies $F_k(N)\le {N-1\choose k-1}+N-k$ and the asymptotic lower bound $F_k(N)=\Omega_k(N^{k-1})$. More precise bounds on $F_k(N)$ are obtained for $k\le 3$. We also obtain the threshold probability for a random system to be Sidon for $k\ge 2$.


title: Engineered Swift Equilibration for Brownian objects: from underdamped to overdamped dynamics

abstract: We propose a general framework to study transformations that drive an underdamped Brownian particle in contact with a thermal bath from an equilibrium state to a new one in an arbitrarily short time. To this end, we make use of a time and space-dependent potential, that plays a dual role: confine the particle, and manipulate the system. In the special case of an isothermal compression or decompression of a harmonically trapped particle, we derive explicit protocols that perform this quick transformation, following an inverse engineering method. We focus on the properties of these protocols, which crucially depend on two key dimensionless numbers that characterize the relative values of the three timescales of the problem, associated with friction, oscillations in the confinement and duration of the protocol. In particular, we show that our protocols encompass the known overdamped version of this problem and extend it to any friction for decompression and to a large range of frictions for compression.


title: Bosonic dark matter halos: excited states and relaxation in the potential of the ground state

abstract: An ultra-light axion field with mass $\sim 10^{-22}\ {\rm eV}$, also known as wave or fuzzy dark matter, has been proposed as a component of the dark matter in the Universe. We study the evolution of the axion dark matter distribution in the central region of a halo, assuming the mass is dominated by this axion field, and that gravity is the only important interaction. We calculate the excited axion states in the spherical gravitational potential generated by the self-gravitating ground-state, also known as soliton. These excited states are similar to the states of the hydrogen atom with quantum numbers $(n,l,m)$, here designating oscillation modes of a classical wave. At fixed $n$, the modes with highest $l$ have the lowest energy because of the extended mass distribution generating the potential. We use an approximate analytical treatment to derive the distribution of mass in these states when a steady-state is reached by dynamical relaxation, and find that a corona with a mass density profile $\rho \propto r^{-5/3}$ should be set up around the central soliton, analogous to the Bahcall-Wolf cusp predicted for the stellar distribution around a central black hole. The central soliton accretes dark matter from the corona as dynamical relaxation proceeds and negative orbital energy flows out. This density profile should remain valid out to the radius where the mass in the corona is comparable to the mass of the central soliton; further than that, the gravitational potential depends on the initial distribution of dark matter and the relaxation time increases rapidly with radius.


title: Toll Caps in Privatized Road Networks

abstract: We consider a network pricing game on a parallel network with congestion effects in which link owners set tolls for travel so as to maximize profit. A central authority is able to regulate this competition by means of a (uniform) price cap. The first question we want to answer is how such a cap should be designed in order to minimize the total congestion. We provide an algorithm that finds an optimal price cap for networks with affine latency functions and a full support Wardrop equilibrium. Second, we consider the induced network performance at an optimal price cap. We show that for two link networks with affine latency functions, the congestion costs at the optimal price cap are at most $8/7$ times the optimal congestion costs. For more general latency functions, this bound goes up to 2 under the assumption that an uncapped Nash equilibrium exists. However, in general such an equilibrium need not exist and this can be used to show that optimal price caps can induce arbitrarily inefficient flows.


title: Stochastic Dynamic Programming Heuristics for Influence Maximization-Revenue Optimization

abstract: The well-known Influence Maximization (IM) problem has been actively studied by researchers over the past decade, with emphasis on marketing and social networks. Existing research have obtained solutions to the IM problem by obtaining the influence spread and utilizing the property of submodularity. This paper is based on a novel approach to the IM problem geared towards optimizing clicks and consequently revenue within anOnline Social Network (OSN). Our approach diverts from existing approaches by adopting a novel, decision-making perspective through implementing Stochastic Dynamic Programming (SDP). Thus, we define a new problem Influence Maximization-Revenue Optimization (IM-RO) and propose SDP as a method in which this problem can be solved. The SDP method has lucrative gains for an advertiser in terms of optimizing clicks and generating revenue however, one drawback to the method is its associated "curse of dimensionality" particularly for problems involving a large state space. Thus, we introduce the Lawrence Degree Heuristic (LDH), Adaptive Hill-Climbing (AHC) and Multistage Particle Swarm Optimization (MPSO) heuristics as methods which are orders of magnitude faster than the SDP method whilst achieving near-optimal results. Through a comparative analysis on various synthetic and real-world networks we present the AHC and LDH as heuristics well suited to to the IM-RO problem in terms of their accuracy, running times and scalability under ideal model parameters. In this paper we present a compelling survey on the SDP method as a practical and lucrative method for spreading information and optimizing revenue within the context of OSNs.


title: $L^2$ estimates and vanishing theorems for holomorphic vector bundles equipped with singular Hermitian metrics

abstract: We investigate singular Hermitian metrics on vector bundles, especially strictly Griffiths positive ones. $L^2$ esitimates and vanishing theorems usually require an assumption that vector bundles are Nakano positive. However there is no general definition of the Nakano positivity in singular settings. In this paper, we show various $L^2$ estimates and vanishing theromes by assuming that the vector bundle is strictly Griffiths positive and the base manifold is projective.


title: Experimental investigation of heat transport in homogeneous bubbly flow

abstract: We present results on the global and local characterisation of heat transport in homogeneous bubbly flow. Experimental measurements were performed with and without the injection of $\sim 2.5$ mm diameter bubbles (corresponding to $Re_b \approx 600$) in a rectangular water column heated from one side and cooled from the other. The gas volume fraction $\alpha$ was varied in the range $0\% - 5\%$, and the Rayleigh number $Ra_H$ in the range $4.0 \times 10^9 - 1.2 \times 10^{11}$. We find that the global heat transfer is enhanced up to 20 times due to bubble injection. Interestingly, for bubbly flow, for our lowest concentration $\alpha = 0.5\% $ onwards, the Nusselt number $\overline{Nu}$ is nearly independent of $Ra_H$, and depends solely on the gas volume fraction~$\alpha$. We observe the scaling $\overline{Nu} \propto \alpha^{0.45}$, which is suggestive of a diffusive transport mechanism. Through local temperature measurements, we show that the bubbles induce a huge increase in the strength of liquid temperature fluctuations, e.g. by a factor of 200 for $\alpha = 0.9\%$. Further, we compare the power spectra of the temperature fluctuations for the single- and two-phase cases. In the single-phase cases, most of the spectral power of the temperature fluctuations is concentrated in the large-scale rolls. However, with the injection of bubbles, we observe intense fluctuations over a wide range of scales, extending up to very high frequencies. Thus, while in the single-phase flow the thermal boundary layers control the heat transport, once the bubbles are injected, the bubble-induced liquid agitation governs the process from a very small bubble concentration onwards.


title: Mapping mesoscopic phase evolution during e-beam induced transformations via deep learning of atomically resolved images

abstract: Understanding transformations under electron beam irradiation requires mapping the structural phases and their evolution in real time. To date, this has mostly been a manual endeavor comprising of difficult frame-by-frame analysis that is simultaneously tedious and prone to error. Here, we turn towards the use of deep convolutional neural networks (DCNN) to automatically determine the Bravais lattice symmetry present in atomically-resolved images. A DCNN is trained to identify the Bravais lattice class given a 2D fast Fourier transform of the input image. Monte-Carlo dropout is used for determining the prediction probability, and results are shown for both simulated and real atomically-resolved images from scanning tunneling microscopy and scanning transmission electron microscopy. A reduced representation of the final layer output allows to visualize the separation of classes in the DCNN and agrees with physical intuition. We then apply the trained network to electron beam-induced transformations in WS2, which allows tracking and determination of growth rate of voids. These results are novel in two ways: (1) It shows that DCNNs can be trained to recognize diffraction patterns, which is markedly different from the typical "real image" cases, and (2) it provides a method with in-built uncertainty quantification, allowing the real-time analysis of phases present in atomically resolved images.


title: On the Lie bracket approximation approach to distributed optimization: Extensions and limitations

abstract: We consider the problem of solving a smooth convex optimization problem with equality and inequality constraints in a distributed fashion. Assuming that we have a group of agents available capable of communicating over a communication network described by a time-invariant directed graph, we derive distributed continuous-time agent dynamics that ensure convergence to a neighborhood of the optimal solution of the optimization problem. Following the ideas introduced in our previous work, we combine saddle-point dynamics with Lie bracket approximation techniques. While the methodology was previously limited to linear constraints and objective functions given by a sum of strictly convex separable functions, we extend these result here and show that it applies to a very general class of optimization problems under mild assumptions on the communication topology.


title: Glass stability (GS) of chemically complex (natural) sub-alkaline glasses

abstract: Glass stability (GS) indicates the glass reluctance or ability to crystallise upon heating; it can be characterised by several methods and parameters and is frequently used to retrieve glass-forming ability (GFA) of corresponding liquids as the case with which such liquids can be made crystal free via melt-quenching. Here, GS has been determined for the first time on six sub-alkaline glasses having complex (natural) compositions, the most widespread and abundant on Earth. KT, KH, KW, KLL and w2 GS parameters increase linearly and monotonically as a function of SiO2, with very high correlations. Moreover, Tx values and GS parameters highly correlate with GFA via Rc (critical cooling rate), previously determined with ex-situ cooling-induced experiments. Therefore, GS scales with GFA for natural silicate compositions. In addition, the in-situ Rc value of B100 measured with DSC results > 45 {\deg}C/min (> 2700 {\deg}C/h), broadly corroborating the Rc of about 150 {\deg}C/min (9000 {\deg}C/h) determined ex-situ. In turn, relevant solidification parameters on heating or cooling can be obtained by DSC investigations also for chemically complex (natural) systems, similar to simple silicate systems. These outcomes are relevant for lavas or magmas that re-heat glass-bearing volcanic rocks, as well as for fabricate glass-ceramic materials with desirable texture and composition of phases starting from abundant and very cheap raw volcanic rocks.


title: More than five-twelfths of the zeros of $\zeta$ are on the critical line

abstract: The second moment of the Riemann zeta-function twisted by a normalized Dirichlet polynomial with coefficients of the form $(\mu \star \Lambda_1^{\star k_1} \star \Lambda_2^{\star k_2} \star \cdots \star \Lambda_d^{\star k_d})$ is computed unconditionally by means of the autocorrelation of ratios of $\zeta$ techniques from Conrey, Farmer, Keating, Rubinstein and Snaith (2005), Conrey, Farmer and Zirnbauer (2008) as well as Conrey and Snaith (2007). This in turn allows us to describe the combinatorial process behind the mollification of \[ \zeta(s) + \lambda_1 \frac{\zeta'(s)}{\log T} + \lambda_2 \frac{\zeta''(s)}{\log^2 T} + \cdots + \lambda_d \frac{\zeta^{(d)}(s)}{\log^d T}, \] where $\zeta^{(k)}$ stands for the $k$th derivative of the Riemann zeta-function and $\{\lambda_k\}_{k=1}^d$ are real numbers. Improving on recent results on long mollifiers and sums of Kloosterman sums due to Pratt and Robles (2017), as an application, we increase the current lower bound of critical zeros of the Riemann zeta-function to slightly over five-twelfths.


title: Universal two-time correlations, out-of-time-ordered correlators and Leggett-Garg inequality violation by edge Majorana fermion qubits

abstract: In the present work we propose that two-time correlations of Majorana edge localized fermions constitute a novel and versatile toolbox for assessing the topological phases of 1D open lattices. Using analytical and numerical calculations on the Kitaev model, we uncover universal relationships between the decay of the short-time correlations and a particular family of out-of-time-ordered correlators, which provide direct experimental alternatives to the quantitative analysis of the system regime, either normal or topological. Furthermore we show that the saturation of two-time correlations possesses features of an order parameter. Finally, we find that violations of Leggett-Garg inequalities can indicate the topological-normal phase transition by looking at different qubits formed by pairing local and non-local edge Majorana fermions.


title: Is Private Browsing in Modern Web Browsers Really Private?

abstract: Web browsers are the most common tool to perform various activities over the internet. Along with normal mode, all modern browsers have private browsing mode. The name of the mode varies from browser to browser but the purpose of the private mode remains same in every browser. In normal browsing mode, the browser keeps track of users' activity and related data such as browsing histories, cookies, auto-filled fields, temporary internet files, etc. In private mode, it is said that no information is stored while browsing or all information is destroyed after closing the current private session. However, some researchers have already disproved this claim by performing various tests in most popular browsers. I have also some personal experience where private mode browsing fails to keep all browsing information as private. In this position paper, I take the position against private browsing. By examining various facts, it is proved that the private browsing mode is not really private as it is claimed; it does not keep everything private. In following sections, I will present the proof to justify my argument. Along with some other already performed research work, I will show my personal case studies and experimental data as well.


title: Quantizing quantum Ricci curvature

abstract: Quantum Ricci curvature has been introduced recently as a new, geometric observable characterizing the curvature properties of metric spaces, without the need for a smooth structure. Besides coordinate invariance, its key features are scalability, computability and robustness. We demonstrate that these properties continue to hold in the context of nonperturbative quantum gravity, by evaluating the quantum Ricci curvature numerically in two-dimensional Euclidean quantum gravity, defined in terms of dynamical triangulations. Despite the well-known, highly nonclassical properties of the underlying quantum geometry, its Ricci curvature can be matched well to that of a five-dimensional round sphere.


title: Searching for supernovae in the multiply-imaged galaxies behind the gravitational telescope A370

abstract: Strong lensing by massive galaxy clusters can provide magnification of the flux and even multiple images of the galaxies that lie behind them. This phenomenon facilitates observations of high-redshift supernovae (SNe), that would otherwise remain undetected. Type Ia supernovae (SNe Ia) detections are of particular interest because of their standard brightness, since they can be used to improve either cluster lensing models or cosmological parameter measurements. We present a ground-based, near-infrared search for lensed SNe behind the galaxy cluster Abell 370. Our survey was based on 15 epochs of J-band observations with the HAWK-I instrument on the Very Large Telescope (VLT). We use Hubble Space Telescope (HST) photometry to infer the global properties of the multiply-imaged galaxies. Using a recently published lensing model of Abell 370, we also present the predicted magnifications and time delays between the images. In our survey, we did not discover any live SNe from the 13 lensed galaxies with 47 multiple images behind Abell 370. This is consistent with the expectation of $0.09\pm0.02$ SNe calculated based on the measured star formation rate. We compare the expectations of discovering strongly lensed SNe in our survey and that performed with HST during the Hubble Frontier Fields (HFF) programme. We also show the expectations of search campaigns that can be conducted with future facilities, such as the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). We show that the NIRCam instrument aboard the JWST will be sensitive to most SN multiple images in the strongly lensed galaxies and thus will be able to measure their time delays if observations are scheduled accordingly.


title: Application of R\'enyi and Tsallis Entropies to Topic Modeling Optimization

abstract: This is full length article (draft version) where problem number of topics in Topic Modeling is discussed. We proposed idea that Renyi and Tsallis entropy can be used for identification of optimal number in large textual collections. We also report results of numerical experiments of Semantic stability for 4 topic models, which shows that semantic stability play very important role in problem topic number. The calculation of Renyi and Tsallis entropy based on thermodynamics approach.


title: Approaching near-perfect state discrimination of photonic Bell states through the use of unentangled ancilla photons

abstract: Despite well-established no-go theorems on a perfect linear optical Bell state analyzer, we find a numerical trend that appears to approach a near-perfect measurement if we incorporate eight or more un-entangled ancilla photons into our device. Following this trend, we begin a promising inductive approach to building an ideal optical Bell measurement device. In the process, we determine that any Bell state analyzer that (even occasionally) bunches all photons into only two of the output modes cannot perform an ideal measurement and we find a set of conditions on our linear optical circuit that prevent this outcome.


title: Dimensional Analysis in Economics: A Study of the Neoclassical Economic Growth Model

abstract: The fundamental purpose of the present research article is to introduce the basic principles of Dimensional Analysis in the context of the neoclassical economic theory, in order to apply such principles to the fundamental relations that underlay most models of economic growth. In particular, basic instruments from Dimensional Analysis are used to evaluate the analytical consistency of the Neoclassical economic growth model. The analysis shows that an adjustment to the model is required in such a way that the principle of dimensional homogeneity is satisfied.


title: Maximum likelihood estimation of a finite mixture of logistic regression models in a continuous data stream

abstract: In marketing we are often confronted with a continuous stream of responses to marketing messages. Such streaming data provide invaluable information regarding message effectiveness and segmentation. However, streaming data are hard to analyze using conventional methods: their high volume and the fact that they are continuously augmented means that it takes considerable time to analyze them. We propose a method for estimating a finite mixture of logistic regression models which can be used to cluster customers based on a continuous stream of responses. This method, which we coin oFMLR, allows segments to be identified in data streams or extremely large static datasets. Contrary to black box algorithms, oFMLR provides model estimates that are directly interpretable. We first introduce oFMLR, explaining in passing general topics such as online estimation and the EM algorithm, making this paper a high level overview of possible methods of dealing with large data streams in marketing practice. Next, we discuss model convergence, identifiability, and relations to alternative, Bayesian, methods; we also identify more general issues that arise from dealing with continuously augmented data sets. Finally, we introduce the oFMLR [R] package and evaluate the method by numerical simulation and by analyzing a large customer clickstream dataset.


title: Orion+: Automated Problem Diagnosis in Computing Systems by Mining Metric Data

abstract: This work presents the suspicious code at a finer granularity of call stack rather than code region, which was being returned by Orion. Call stack based comparison returns call stacks that are most impacted by the bug and save developer time to debug from scratch. This solution has polynomial complexity and hence can be implemented practically.


title: Satellite ruling polynomials, DGA representations, and the colored HOMFLY-PT polynomial

abstract: We establish relationships between two classes of invariants of Legendrian knots in $\mathbb{R}^3$: Representation numbers of the Chekanov-Eliashberg DGA and satellite ruling polynomials. For positive permutation braids, $\beta \subset J^1S^1$, we give a precise formula in terms of representation numbers for the $m$-graded ruling polynomial $R^m_{S(K,\beta)}(z)$ of the satellite of $K$ with $\beta$ specialized at $z=q^{1/2}-q^{-1/2}$ with $q$ a prime power, and we use this formula to prove that arbitrary $m$-graded satellite ruling polynomials, $R^m_{S(K,L)}$, are determined by the Chekanov-Eliashberg DGA of $K$. Conversely, for $m\neq 1$, we introduce an $n$-colored $m$-graded ruling polynomial, $R^m_{n,K}(q)$, in strict analogy with the $n$-colored HOMFLY-PT polynomial, and show that the total $n$-dimensional $m$-graded representation number of $K$ to $\mathbb{F}_q^n$, $\mbox{Rep}_m(K,\mathbb{F}_q^n)$, is exactly equal to $R^m_{n,K}(q)$. In the case of $2$-graded representations, we show that $R^2_{n,K}=\mbox{Rep}_2(K, \mathbb{F}_q^n)$ arises as a specialization of the $n$-colored HOMFLY-PT polynomial.


title: Higgs and $Z$ Assisted Stop Searches at Hadron Colliders

abstract: Current searches for the light top squark (stop) mostly focus on the decay channels of $\tilde{t} \rightarrow t \chi_1^0$ or $\tilde{t} \rightarrow b \chi_1^\pm \rightarrow bW \chi_1^0$, leading to $t\bar{t}/bbWW+\met$ final states for stop pair productions at the LHC. However, in supersymmetric scenarios with light neutralinos and charginos other than the neutralino lightest supersymmetric particle (LSP), more than one decay mode of the stop could be dominant. While those new decay modes could significantly weaken the current stop search limits at the LHC, they also offer alternative discovery channels for stop searches. In this paper, we studied the scenario with light Higgsino next-to-LSPs (NLSPs) and Bino LSP. The light stop decays primarily via $\tilde t_1 \to t \chi_2^0/\chi_3^0$, with the neutralinos subsequent decaying to a $Z$ boson or a Higgs boson: $\chi_2^0/\chi^0_3 \to \chi_1^0 h/Z$. Pair production of light stops at the LHC leads to final states of $t \bar t hh\met$, $t \bar t hZ\met$ or $t \bar t ZZ\met$. We consider three signal regions: one charged lepton (1$\ell$), two opposite sign charged leptons (2 OS $\ell$) and at least three charged leptons ($ \ge 3 \ell$). We found that the 1$\ell$ signal region of channel $t \bar t hZ\met$ has the best reach sensitivity for light stop searches. For 14 TeV LHC with 300 ${\rm fb}^{-1}$ integrated luminosity, a stop mass up to 900 GeV can be discovered at 5$\sigma$ significance, or up to 1050 GeV can be excluded at 95\% C.L. Combining all three decay channels for $1 \ell$ signal region extends the reach for about 100$-$150 GeV. We also studied the stop reach at the 100 TeV $pp$ collider with 3 ${\rm ab}^{-1}$ luminosity, with discovery and exclusion reach being 6 TeV and 7 TeV, respectively.


title: Modern Theory for the Orbital Moment in a Superconductor

abstract: The chiral p-wave superconducting state is comprised of spin triplet Cooper pairs carrying a finite orbital angular momentum. For the case of a periodic lattice, calculating the net magnetisation arising from this orbital component presents a challenge as the circulation operator $\hat{\bf{r}} \times \hat{\bf{p}}$ is not well defined in the Bloch representation. This difficulty has been overcome in the normal state, for which a modern theory is firmly established. Here, we derive the extension of this normal state approach, generating a theory which is valid for a general superconducting state, and go on to perform model calculations for a chiral p-wave state in Sr$_2$RuO$_4$. The results suggest that the magnitude of the elusive edge current in Sr$_2$RuO$_4$ is finite, but lies below experimental resolution. This provides a possible solution to the long-standing controversy concerning the gap symmetry of the superconducting state in this material.


title: Voigt effect-based wide-field magneto-optical microscope integrated in a pump-probe experimental setup

abstract: In this work we describe an experimental setup for spatially-resolved pump-probe experiment with integrated wide-field magneto-optical (MO) microscope. The MO microscope can be used to study ferromagnetic materials with both perpendicular-to-plane and in-plane magnetic anisotropy via polar Kerr and Voigt effects, respectively. The functionality of the Voigt effect-based microscope was tested using an in-plane magnetized ferromagnetic semiconductor (Ga,Mn)As. It was revealed that the presence of mechanical defects in the (Ga,Mn)As epilayer alters significantly the magnetic anisotropy in their proximity. The importance of MO experiments with simultaneous temporal and spatial resolutions was demonstrated using (Ga,Mn)As sample attached to a piezoelectric actuator, which produces a voltage-controlled strain. We observed a considerably different behavior in different parts of the sample that enabled us to identify sample parts where the epilayer magnetic anisotropy was significantly modified by a presence of the piezostressor and where it was not. Finally, we discuss the possible applicability of our experimental setup for the research of compensated antiferromagnets, where only MO effects even in magnetic moments are present.


title: Subsums of conditionally convergent series in finite dimensional spaces

abstract: An achievement set of a series is a set of all its subsums. We study the properties of achievement sets of conditionally convergent series in finite dimensional spaces. The purpose of the paper is to answer some of the open problems formulated in \cite{GM}. We obtain general results for series with harmonic-like coordinates, that is $A((-1)^{n+1}n^{-\alpha_1},\dots,(-1)^{n+1}n^{-\alpha_d})=\mathbb{R}^d$ for pairwise distinct numbers $\alpha_1,\dots,\alpha_d\in(0,1]$. For $d=2$, $\alpha_1=1, \alpha_2=\frac{1}{2}$ it was stated as an open problem in \cite{GM}, that is $A(\frac{(-1)^n}{n},\frac{(-1)^n}{\sqrt{n}})=\mathbb{R}^2$.


title: Photonic glass for high contrast structural color

abstract: Non-iridescent structural colors based on disordered arrangement of monodisperse spherical particles, also called photonic glass, show low color saturation due to gradual transition in reflectivity. No significant improvement is usually expected from particles optimization, as the Mie resonances are broad for small dielectric particles with moderate refractive index. Moreover, the short range order of a photonic glass alone is also insufficient to cause sharp spectral features. We show here, that the combination of a well-chosen particle geometry with the short range order of a photonic glass has strong synergetic effects. We demonstrate how core-shell particles can be used to obtain a sharp transition in the reflection spectrum of photonic glass which is essential to achieve a strong color saturation. The Fourier transform required for a highly saturated color can be achieved by shifting the first zero position of the motif Fourier transform to smaller wave numbers in respect to the peak of the lattice Fourier transform. We show that this can be obtained by choosing a non-monotonous refractive index distribution from the center of the particle through the shell and into the background material. The first-order theoretical predictions are confirmed by numerical simulations.


title: Thermal Conformal Blocks

abstract: We study conformal blocks for thermal one-point-functions on the sphere in conformal field theories of general dimension. These thermal conformal blocks satisfy second order Casimir differential equations and have integral representations related to AdS Witten diagrams. We give an analytic formula for the scalar conformal block in terms of generalized hypergeometric functions. As an application, we deduce an asymptotic formula for the three-point coeffcients of primary operators in the limit where two of the operators are heavy.


title: Multiple stellar populations in Magellanic Cloud clusters. VI. A survey of multiple sequences and Be stars in young clusters

abstract: The split main sequences (MSs) and extended MS turnoffs (eMSTOs) detected in a few young clusters have demonstrated that these stellar systems host multiple populations differing in a number of properties such as rotation and, possibly, age.We analyze Hubble Space Telescope photometry for thirteen clusters with ages between ~40 and ~1000 Myrs and of different masses. Our goal is to investigate for the first time the occurrence of multiple populations in a large sample of young clusters. We find that all the clusters exhibit the eMSTO phenomenon and that MS stars more massive than ~1.6 solar masses define a blue and red MS, with the latter hosting the majority of MS stars. The comparison between the observations and isochrones suggests that the blue MSs are made of slow-rotating stars, while the red MSs host stars with rotational velocities close to the breakup value. About half of the bright MS stars in the youngest clusters are H-alpha emitters. These Be stars populate the red MS and the reddest part of the eMSTO thus supporting the idea that the red MS is made of fast rotators. We conclude that the split MS and the eMSTO are a common feature of young clusters in both Magellanic Clouds. The phenomena of a split MS and an eMSTO occur for stars that are more massive than a specific threshold which is independent of the host-cluster mass. As a by-product, we report the serendipitous discovery of a young SMC cluster, GALFOR1.


title: From synchronous to one-time delayed dynamics in coupled maps

abstract: We study the completely synchronized states (CSSs) of a system of coupled logistic maps as a function of three parameters: interaction strength ($\varepsilon$), range of the interaction ($\alpha$), that can vary from first-neighbors to global coupling, and a parameter ($\beta$) that allows to scan continuously from non-delayed to one-time delayed dynamics. % We identify in the plane $\alpha$-$\varepsilon$ periodic orbits, limit cycles and chaotic trajectories, and describe how these structures change with the delay. These features can be explained by studying the bifurcation diagrams of a two-dimensional non-delayed map. This allows us to understand the effects of one-time delays on CSSs, e.g, regularization of chaotic orbits and synchronization of short-range coupled maps, observed when the dynamics is moderately delayed. Finally, we substitute the logistic map by cubic and logarithmic maps, in order to test the robustness of our findings.


title: Thresholds of Braided Convolutional Codes on the AWGN Channel

abstract: In this paper, we perform a threshold analysis of braided convolutional codes (BCCs) on the additive white Gaussian noise (AWGN) channel. The decoding thresholds are estimated by Monte-Carlo density evolution (MC-DE) techniques and compared with approximate thresholds from an erasure channel prediction. The results show that, with spatial coupling, the predicted thresholds are very accurate and quickly approach capacity if the coupling memory is increased. For uncoupled ensembles with random puncturing, the prediction can be improved with help of the AWGN threshold of the unpunctured ensemble.


title: Conductance signatures of odd-frequency superconductivity in quantum spin Hall systems using a quantum point contact

abstract: Topological superconductors give rise to unconventional superconductivity, which is mainly characterized by the symmetry of the superconducting pairing amplitude. However, since the symmetry of the superconducting pairing amplitude is not directly observable, its experimental identification is rather difficult. In our work, we propose a system, composed of a quantum point contact and proximity induced s-wave superconductivity at the helical edge of a two dimensional topological insulator, for which we demonstrate the presence of odd-frequency pairing and its intimate connection to unambiguous transport signatures. Notably, our proposal requires no time-reversal symmetry breaking terms. We discover the domination of crossed Andreev reflection over electron cotunneling in a wide range of parameter space, which is a quite unusual transport regime.


title: Memory-based Parameter Adaptation

abstract: Deep neural networks have excelled on a wide range of problems, from vision to language and game playing. Neural networks very gradually incorporate information into weights as they process data, requiring very low learning rates. If the training distribution shifts, the network is slow to adapt, and when it does adapt, it typically performs badly on the training distribution before the shift. Our method, Memory-based Parameter Adaptation, stores examples in memory and then uses a context-based lookup to directly modify the weights of a neural network. Much higher learning rates can be used for this local adaptation, reneging the need for many iterations over similar data before good predictions can be made. As our method is memory-based, it alleviates several shortcomings of neural networks, such as catastrophic forgetting, fast, stable acquisition of new knowledge, learning with an imbalanced class labels, and fast learning during evaluation. We demonstrate this on a range of supervised tasks: large-scale image classification and language modelling.


title: A Frequent Itemset Hiding Toolbox

abstract: Advances in data collection and data storage technologies have given way to the establishment of transactional databases among companies and organizations, as they allow enormous amounts of data to be stored efficiently. Useful knowledge can be mined from these data, which can be used in several ways depending on the nature of the data. Quite often companies and organizations are willing to share data for the sake of mutual benefit. However, the sharing of such data comes with risks, as problems with privacy may arise. Sensitive data, along with sensitive knowledge inferred from this data, must be protected from unintentional exposure to unauthorized parties. One form of the inferred knowledge is frequent patterns mined in the form of frequent itemsets from transactional databases. The problem of protecting such patterns is known as the frequent itemset hiding problem. In this paper we present a toolbox, which provides several implementations of frequent itemset hiding algorithms. Firstly, we summarize the most important aspects of each algorithm. We then introduce the architecture of the toolbox and its novel features. Finally, we provide experimental results on real world datasets, demonstrating the efficiency of the toolbox and the convenience it offers in comparing different algorithms.


title: A Fractional Variational Approach for Modelling Dissipative Mechanical Systems: Continuous and Discrete Settings

abstract: Employing a phase space which includes the (Riemann-Liouville) fractional derivative of curves evolving on real space, we develop a restricted variational principle for Lagrangian systems yielding the so-called restricted fractional Euler-Lagrange equations (both in the continuous and discrete settings), which, as we show, are invariant under linear change of variables. This principle relies on a particular restriction upon the admissible variation of the curves. In the case of the half-derivative and mechanical Lagrangians, i.e. kinetic minus potential energy, the restricted fractional Euler-Lagrange equations model a dissipative system in both directions of time, summing up to a set of equations that is invariant under time reversal. Finally, we show that the discrete equations are a meaningful discretisation of the continuous ones.


title: A spectral collocation method for nonlocal diffusion equations

abstract: Nonlocal diffusion model provides an appropriate description of the diffusion process of solute in the complex medium, which cannot be described properly by classical theory of PDE. However, the operators in the nonlocal diffusion models are nonlocal, so the resulting numerical methods generate dense or full stiffness matrices. This imposes significant computational and memory challenge for a nonlocal diffusion model. In this paper, we develop a spectral collocation method for the nonlocal diffusion model and provide a rigorous error analysis which theoretically justifies the spectral rate of convergence provided that the kernel functions and the source functions are sufficiently smooth. Compared to finite difference methods and finite element methods, because of the high order convergence rates, the numerical cost of spectral collocation methods will be greatly decreased. Numerical results confirm the exponential rate of convergence.


title: Computational Theories of Curiosity-Driven Learning

abstract: What are the functions of curiosity? What are the mechanisms of curiosity-driven learning? We approach these questions about the living using concepts and tools from machine learning and developmental robotics. We argue that curiosity-driven learning enables organisms to make discoveries to solve complex problems with rare or deceptive rewards. By fostering exploration and discovery of a diversity of behavioural skills, and ignoring these rewards, curiosity can be efficient to bootstrap learning when there is no information, or deceptive information, about local improvement towards these problems. We also explain the key role of curiosity for efficient learning of world models. We review both normative and heuristic computational frameworks used to understand the mechanisms of curiosity in humans, conceptualizing the child as a sense-making organism. These frameworks enable us to discuss the bi-directional causal links between curiosity and learning, and to provide new hypotheses about the fundamental role of curiosity in self-organizing developmental structures through curriculum learning. We present various developmental robotics experiments that study these mechanisms in action, both supporting these hypotheses to understand better curiosity in humans and opening new research avenues in machine learning and artificial intelligence. Finally, we discuss challenges for the design of experimental paradigms for studying curiosity in psychology and cognitive neuroscience. Keywords: Curiosity, intrinsic motivation, lifelong learning, predictions, world model, rewards, free-energy principle, learning progress, machine learning, AI, developmental robotics, development, curriculum learning, self-organization.


title: Dynamic Pricing with Variable Order Sizes for a Model with Constant Demand Elasticity

abstract: In this paper we investigate a dynamic pricing model for constant demand elasticity where customers have a probability distribution on the number of items they order. This is a generalization from standard models which restrict customers to buy only one item at a time. For the generalized model, we first obtain a closed form expression for the optimal expected revenue and optimal pricing strategy. This expression involves a recursively defined term for which we investigate the behavior. We call comparable models those which have the same demand, which is the customer arrival rate times the average order size. In fact, the average order size plays an important role for results for the generalized model. An important result we show is that comparable models have the same asymptotic pricing behavior. Numerical results also show that comparable models are relatively close even for low inventory levels. Lastly, we prove that the relative difference between comparable models is governed not by the customer arrival rate, but solely by their order size distributions.


title: Using Deep Learning for Segmentation and Counting within Microscopy Data

abstract: Cell counting is a ubiquitous, yet tedious task that would greatly benefit from automation. From basic biological questions to clinical trials, cell counts provide key quantitative feedback that drive research. Unfortunately, cell counting is most commonly a manual task and can be time-intensive. The task is made even more difficult due to overlapping cells, existence of multiple focal planes, and poor imaging quality, among other factors. Here, we describe a convolutional neural network approach, using a recently described feature pyramid network combined with a VGG-style neural network, for segmenting and subsequent counting of cells in a given microscopy image.


title: Automatic topography of high-dimensional data sets by non-parametric Density Peak clustering

abstract: Data analysis in high-dimensional spaces aims at obtaining a synthetic description of a data set, revealing its main structure and its salient features. We here introduce an approach providing this description in the form of a topography of the data, namely a human-readable chart of the probability density from which the data are harvested. The approach is based on an unsupervised extension of Density Peak clustering and a non-parametric density estimator that measures the probability density in the manifold containing the data. This allows finding automatically the number and the height of the peaks of the probability density, and the depth of the "valleys" separating them. Importantly, the density estimator provides a measure of the error, which allows distinguishing genuine density peaks from density fluctuations due to finite sampling. The approach thus provides robust and visual information about the density peaks' height, their statistical reliability, and their hierarchical organization, offering a conceptually powerful extension of the standard clustering partitions. We show that this framework is particularly useful in the analysis of complex data sets.


title: Exclusive double quarkonium production and generalized TMDs of gluons

abstract: Being the "mother distributions" of all types of two-parton correlation functions, generalized transverse momentum dependent parton distributions (GTMDs) have attracted a lot of attention over the last years. We argue that exclusive double production of pseudoscalar quarkonia ($\eta_c$ or $\eta_b$) in nucleon-nucleon collisions gives access to GTMDs of gluons.


title: A Variational Inequality Perspective on Generative Adversarial Networks

abstract: Generative adversarial networks (GANs) form a generative modeling approach known for producing appealing samples, but they are notably difficult to train. One common way to tackle this issue has been to propose new formulations of the GAN objective. Yet, surprisingly few studies have looked at optimization methods designed for this adversarial training. In this work, we cast GAN optimization problems in the general variational inequality framework. Tapping into the mathematical programming literature, we counter some common misconceptions about the difficulties of saddle point optimization and propose to extend techniques designed for variational inequalities to the training of GANs. We apply averaging, extrapolation and a computationally cheaper variant that we call extrapolation from the past to the stochastic gradient method (SGD) and Adam.


title: Equi-coverage Contours in Cellular Networks

abstract: In this letter, we introduce a general cellular network model where i) users and BSs are distributed as two general point processes that may be coupled, ii) pathloss is assumed to follow a multi-slope power-law pathloss model, and iii) fading (power) is assumed to be independent across all wireless links. For this setup, we first obtain a set of contours representing the same meta distribution of SIR, which is the distribution of the conditional coverage probability given the point process, for different values of the parameters of the pathloss function and BS and user point processes. This general result is then specialized to 3GPP-inspired user and BS configurations obtained by combining Poisson point process (PPP) and Poisson cluster process (PCP).


title: Topological phase transition in a two-species fermion system: Effects of a rotating trap potential or a synthetic gauge field

abstract: We numerically investigate the quantum phases and phase transition in a system made of two species of fermionic atoms that interact with each other via $s$-wave Feshbach resonance, and are subject to rotation or a synthetic gauge field that puts the fermions at Landau level filling factor $\nu_f = 2$. We show that the system undergoes a continuous quantum phase transition from a $\nu_f = 2$ fermionic integer quantum Hall state formed by atoms, to a $\nu_b = 1/2$ bosonic fractional quantum Hall state formed by bosonic diatomic molecules. In the disk geometry we use, these two different topological phases are distinguished by their different gapless edge excitation spectra, and the quantum phase transition between them is signaled by the closing of the energy gap in the bulk. Comparisons will be made with field theoretical predictions, and the case of $p$-wave pairing.


title: Retrieval and Registration of Long-Range Overlapping Frames for Scalable Mosaicking of In Vivo Fetoscopy

abstract: Purpose: The standard clinical treatment of Twin-to-Twin Transfusion Syndrome consists in the photo-coagulation of undesired anastomoses located on the placenta which are responsible to a blood transfer between the two twins. While being the standard of care procedure, fetoscopy suffers from a limited field-of-view of the placenta resulting in missed anastomoses. To facilitate the task of the clinician, building a global map of the placenta providing a larger overview of the vascular network is highly desired. Methods: To overcome the challenging visual conditions inherent to in vivo sequences (low contrast, obstructions or presence of artifacts, among others), we propose the following contributions: (i) robust pairwise registration is achieved by aligning the orientation of the image gradients, and (ii) difficulties regarding long-range consistency (e.g. due to the presence of outliers) is tackled via a bag-of-word strategy, which identifies overlapping frames of the sequence to be registered regardless of their respective location in time. Results: In addition to visual difficulties, in vivo sequences are characterised by the intrinsic absence of gold standard. We present mosaics motivating qualitatively our methodological choices and demonstrating their promising aspect. We also demonstrate semi-quantitatively, via visual inspection of registration results, the efficacy of our registration approach in comparison to two standard baselines. Conclusion: This paper proposes the first approach for the construction of mosaics of placenta in in vivo fetoscopy sequences. Robustness to visual challenges during registration and long-range temporal consistency are proposed, offering first positive results on in vivo data for which standard mosaicking techniques are not applicable.


title: Continuity of Functional Transducers: A Profinite Study of Rational Functions

abstract: A word-to-word function is continuous for a class of languages~$\mathcal{V}$ if its inverse maps $\mathcal{V}$_languages to~$\mathcal{V}$. This notion provides a basis for an algebraic study of transducers, and was integral to the characterization of the sequential transducers computable in some circuit complexity classes. Here, we report on the decidability of continuity for functional transducers and some standard classes of regular languages. To this end, we develop a robust theory rooted in the standard profinite analysis of regular languages. Since previous algebraic studies of transducers have focused on the sole structure of the underlying input automaton, we also compare the two algebraic approaches. We focus on two questions: When are the automaton structure and the continuity properties related, and when does continuity propagate to superclasses?


title: The hierarchy of Poisson brackets for the open Toda lattice and its' spectral curves

abstract: We establish a new representation of the infinite hierarchy of Pois- son brackets (PB) for the open Toda lattice in terms of its spectral curve. For the classical Poisson bracket (PB) we give a representation in the form of a contour integral of some special Abelian differential (meromorphic one-form) on the spectral curve. All higher brackets of the infinite hierarchy are obtained by multiplication of the one-form by a power of the spectral parameter.


title: General-type discrete self-adjoint Dirac systems: explicit solutions of direct and inverse problems, asymptotics of Verblunsky-type coefficients and stability of solving inverse problem

abstract: We consider discrete self-adjoint Dirac systems determined by the potentials (sequences) $\{C_k\}$ such that the matrices $C_k$ are positive definite and $j$-unitary, where $j$ is a diagonal $m\times m$ matrix and has $m_1$ entries $1$ and $m_2$ entries $-1$ ($m_1+m_2=m$) on the main diagonal. We construct systems with rational Weyl functions and explicitly solve inverse problem to recover systems from the contractive rational Weyl functions. Moreover, we study the stability of this procedure. The matrices $C_k$ (in the potentials) are so called Halmos extensions of the Verblunsky-type coefficients $\rho_k$. We show that in the case of the contractive rational Weyl functions the coefficients $\rho_k$ tend to zero and the matrices $C_k$ tend to the indentity matrix $I_m$.


title: Exactly Robust Kernel Principal Component Analysis

abstract: Robust principal component analysis (RPCA) can recover low-rank matrices when they are corrupted by sparse noises. In practice, many matrices are, however, of high-rank and hence cannot be recovered by RPCA. We propose a novel method called robust kernel principal component analysis (RKPCA) to decompose a partially corrupted matrix as a sparse matrix plus a high or full-rank matrix with low latent dimensionality. RKPCA can be applied to many problems such as noise removal and subspace clustering and is still the only unsupervised nonlinear method robust to sparse noises. Our theoretical analysis shows that, with high probability, RKPCA can provide high recovery accuracy. The optimization of RKPCA involves nonconvex and indifferentiable problems. We propose two nonconvex optimization algorithms for RKPCA. They are alternating direction method of multipliers with backtracking line search and proximal linearized minimization with adaptive step size. Comparative studies in noise removal and robust subspace clustering corroborate the effectiveness and superiority of RKPCA.


title: Quantum work for sudden quenches in Gaussian random Hamiltonians

abstract: In the context of nonequilibrium quantum thermodynamics, variables like work behave stochastically. A particular definition of the work probability density function (pdf) for coherent quantum processes allows the verification of the quantum version of the celebrated fluctuation theorems, due to Jarzynski and Crooks, that apply when the system is driven away from an initial equilibrium thermal state. Such a particular pdf depends basically on the details of the initial and final Hamiltonians, on the temperature of the initial thermal state and on how some external parameter is changed during the coherent process. Using random matrix theory we derive a simple analytic expression that describes the general behavior of the work characteristic function $G(u)$, associated with this particular work pdf for sudden quenches, valid for all the traditional Gaussian ensembles of Hamiltonians matrices. This formula well describes the general behavior of $G(u)$ calculated from single draws of the initial and final Hamiltonians in all ranges of temperatures.


title: Novelty Detection with GAN

abstract: The ability of a classifier to recognize unknown inputs is important for many classification-based systems. We discuss the problem of simultaneous classification and novelty detection, i.e. determining whether an input is from the known set of classes and from which specific class, or from an unknown domain and does not belong to any of the known classes. We propose a method based on the Generative Adversarial Networks (GAN) framework. We show that a multi-class discriminator trained with a generator that generates samples from a mixture of nominal and novel data distributions is the optimal novelty detector. We approximate that generator with a mixture generator trained with the Feature Matching loss and empirically show that the proposed method outperforms conventional methods for novelty detection. Our findings demonstrate a simple, yet powerful new application of the GAN framework for the task of novelty detection.


title: Entropy ratio for infinite sequences with positive entropy

abstract: The complexity function of an infinite word $w$ on a finite alphabet $A$ is the sequence counting, for each non-negative $n$, the number of words of length $n$ on the alphabet $A$ that are factors of the infinite word $w$. For any given function $f$ with exponential growth, we introduced in [MM17] the notion of {\it word entropy} $E_W(f)$ associated to $f$ and we described the combinatorial structure of sets of infinite words with a complexity function bounded by $f$. The goal of this work is to give estimates on the word entropy $E_W(f)$ in terms of the limiting lower exponential growth rate of $f$.


title: Splitting of integer polynomials over fields of prime order

abstract: It is well known that a polynomial $\phi(X)\in \mathbb{Z}[X]$ of given degree $d$ factors into at most $d$ factors in $\mathbb{F}_p$ for any prime $p$. We prove in this paper the existence of infinitely many primes $q$ so that the given polynomial $\phi$(X) splits into exactly $d$ linear factors in $\mathbb{F}_q$ by using only elementary results in field theory and some elementary number theory by proving that $\phi$ splits in $\mathbb{F}_q$ iff $P$ has a root in $\mathbb{F}_q$ for all sufficiently large primes $q$, where $P\in \mathbb{Z}[X]$ is any polynomial such that $P$ has a root $\beta \in \mathbb{C}$ for which $\mathbb{Q}(\beta)$ is the splitting field of $\phi$ over $\mathbb{Q}$. Furthermore, we prove that any such $P$ splits in $\mathbb{F}_r$ iff it has a root in $\mathbb{F}_r$, for all sufficiently large primes $r$. Existence of infinitely many such $P$ for any given $\phi$ is also proven.


title: Marstrand type projection theorems for normed spaces

abstract: We consider Marstrand type projection theorems for closest-point projections in the normed space $\mathbb{R}^2$. We prove that if a norm on $\mathbb{R}^2$ is regular enough, then the analogues of the well-known statements from the Euclidean setting hold, while they fail for norms whose unit balls have corners. We establish our results by verifying Peres and Schlag's transversality property and thereby also obtain a Besicovitch-Federer type characterization of purely unrectifiable sets.


title: Evaluation of an Integral

abstract: The Moll-Arias de Reyna integral [1] $$\int_0^{\infty}\frac{dx}{(x^2+1)^{3/2}}\frac{1}{\sqrt{\varphi(x)+\sqrt{\varphi(x)}}}$$ $$\varphi(x)=1+\frac{4}{3}\left(\frac{x}{x^2+1}\right)^2$$ is generalised and several values are given.


title: An illumination effect and an eccentric orbit for the symbiotic binary PU Vul revealed by 32 years of optical spectroscopy

abstract: We analyze $\sim$32 years of optical spectra and photometry for the symbiotic binary PU Vul. Light curves for the He I $\lambda$4471, He II $\lambda$4686 and H$\beta$ $\lambda$4861 emission lines reveal an illumination effect, where the hot white dwarf ionizes the outflowing wind of the red giant, and evidence for an eccentric orbit with e $\approx$ 0.16. Along with the gradual appearance of high ionization emission from [Fe VII] and O VI, the relative fluxes of these lines suggest an increase in the effective temperature of the hot component, from roughly $10^5$ K on JD 2448000 (1990) to roughly $2\times 10^5$ K on JD 2455000 (2009). During this period, the luminosity of the hot component dropped by a factor of 4$-$6 to a current value of roughly 1000 L$_{\odot}$.


title: Characterizing Demand Graphs for (Fixed-Parameter) Shallow-Light Steiner Network

abstract: We consider the Shallow-Light Steiner Network problem from a fixed-parameter perspective. Given a graph $G$, a distance bound $L$, and $p$ pairs of vertices $(s_1,t_1),\cdots,(s_p,t_p)$, the objective is to find a minimum-cost subgraph $G'$ such that $s_i$ and $t_i$ have distance at most $L$ in $G'$ (for every $i \in [p]$). Our main result is on the fixed-parameter tractability of this problem with parameter $p$. We exactly characterize the demand structures that make the problem "easy", and give FPT algorithms for those cases. In all other cases, we show that the problem is W$[1]$-hard. We also extend our results to handle general edge lengths and costs, precisely characterizing which demands allow for good FPT approximation algorithms and which demands remain W$[1]$-hard even to approximate.


title: Learning by Playing - Solving Sparse Reward Tasks from Scratch

abstract: We propose Scheduled Auxiliary Control (SAC-X), a new learning paradigm in the context of Reinforcement Learning (RL). SAC-X enables learning of complex behaviors - from scratch - in the presence of multiple sparse reward signals. To this end, the agent is equipped with a set of general auxiliary tasks, that it attempts to learn simultaneously via off-policy RL. The key idea behind our method is that active (learned) scheduling and execution of auxiliary policies allows the agent to efficiently explore its environment - enabling it to excel at sparse reward RL. Our experiments in several challenging robotic manipulation settings demonstrate the power of our approach.


title: Machine learning and genomics: precision medicine vs. patient privacy

abstract: Machine learning can have major societal impact in computational biology applications. In particular, it plays a central role in the development of precision medicine, whereby treatment is tailored to the clinical or genetic features of the patient. However, these advances require collecting and sharing among researchers large amounts of genomic data, which generates much concern about privacy. Researchers, study participants and governing bodies should be aware of the ways in which the privacy of participants might be compromised, as well as of the large body of research on technical solutions to these issues. We review how breaches in patient privacy can occur, present recent developments in computational data protection, and discuss how they can be combined with legal and ethical perspectives to provide secure frameworks for genomic data sharing.


title: Simultaneously Self-Attending to All Mentions for Full-Abstract Biological Relation Extraction

abstract: Most work in relation extraction forms a prediction by looking at a short span of text within a single sentence containing a single entity pair mention. This approach often does not consider interactions across mentions, requires redundant computation for each mention pair, and ignores relationships expressed across sentence boundaries. These problems are exacerbated by the document- (rather than sentence-) level annotation common in biological text. In response, we propose a model which simultaneously predicts relationships between all mention pairs in a document. We form pairwise predictions over entire paper abstracts using an efficient self-attention encoder. All-pairs mention scores allow us to perform multi-instance learning by aggregating over mentions to form entity pair representations. We further adapt to settings without mention-level annotation by jointly training to predict named entities and adding a corpus of weakly labeled data. In experiments on two Biocreative benchmark datasets, we achieve state of the art performance on the Biocreative V Chemical Disease Relation dataset for models without external KB resources. We also introduce a new dataset an order of magnitude larger than existing human-annotated biological information extraction datasets and more accurate than distantly supervised alternatives.


title: Statistical shape analysis in a Bayesian framework for shapes in two and three dimensions

abstract: In this paper, we describe a novel shape classification method which is embedded in the Bayesian paradigm. We discuss the modelling and the resulting shape classification algorithm for two and three dimensional data shapes. We conclude by evaluating the efficiency and efficacy of the proposed algorithm on the Kimia shape database for the two dimensional case.


title: Spatially-correlated Site Occupancy in the Nonstoichiometric Meta-stable {\epsilon}-Al60Sm11 Phase during Devitrification of Al-10.2 at.% Sm Glasses

abstract: A metastable {\epsilon}-Al60Sm11 phase appears during the initial devitrification of as-quenched Al-10.2 at.% Sm glasses. The {\epsilon} phase is nonstoichiometric in nature since Al occupation is observed on the 16f Sm lattice sites. Scanning transmission electron microscopic images reveal profound spatial correlation of Sm content on these sites, which cannot be explained by the "average crystal" description from Rietveld analysis of diffraction data. Thermodynamically favorable configurations, established by Monte Carlo (MC) simulations based on a cluster-expansion model, also give qualitatively different correlation functions from experimental observations. On the other hand, molecular dynamics simulations of the growth of {\epsilon}-Al60Sm11 in undercooled liquid show that when the diffusion range of Sm is limited to ~ 4 {\AA}, the correlation function of the as-grown crystal structure agrees well with that of the STEM images. Our results show that kinetic effects, especially the limited diffusivity of Sm atoms plays the fundamental role in determining the nonstoichiometric site occupancies of the {\epsilon}-Al60Sm11 phase during the crystallization process.


title: Optimal work extraction from quantum states by photo-assisted Cooper pair tunneling

abstract: The theory of quantum thermodynamics predicts fundamental bounds on work extraction from quantum states. As these bounds are derived in a very general and abstract setting, it is unclear how relevant they are in an experimental context, where control is typically limited. Here we address this question by showing that optimal work extraction is possible for a realistic engine. The latter consists of a superconducting circuit, where a LC-resonator is coupled to a Josephson junction. The oscillator state fuels the engine, providing energy absorbed by Cooper pairs, thus producing work in the form of an electrical current against an external voltage bias. We show that this machine can extract the maximal amount of work from all Gaussian and Fock states. Furthermore, we consider work extraction from a continuously stabilized oscillator state. In both scenarios, coherence between energy eigenstates is beneficial, increasing the power output of the machine. This is possible because the phase difference across the Josephson junction provides a phase reference.


title: Modelling the diffuse dust emission around Orion

abstract: We have studied the diffuse radiation in the surroundings of M42 using photometric data from the Galaxy Evolution Explorer (GALEX) in the far-ultraviolet (FUV) and infrared observations of the AKARI space telescope. The main source of the FUV diffuse emission is the starlight from the Trapezium stars scattered by dust in front of the nebula. We initially compare the diffuse FUV with the far-infrared (FIR) observations at the same locations. The FUV-IR correlations enable us to determine the type of dust contributing to this emission. We then use an existing model for studying the FUV dust scattering in Orion to check if it can be extended to regions away from the centre in a 10 deg radius. We obtain an albedo, {\alpha} = 0.7 and scattering phase function asymmetry factor, g = 0.6 as the median values for our dust locations on different sides of the central Orion region. We find a uniform value of optical parameters across our sample of locations with the dust properties varying significantly from those at the centre of the nebula.


title: Sparse Tensor Algebra Optimizations with Workspaces

abstract: This paper shows how to optimize sparse tensor algebraic expressions by introducing temporary tensors, called workspaces, into the resulting loop nests. We develop a new intermediate language for tensor operations called concrete index notation that extends tensor index notation. Concrete index notation expresses when and where sub-computations occur and what tensor they are stored into. We then describe the workspace optimization in this language, and how to compile it to sparse code by building on prior work in the literature. We demonstrate the importance of the optimization on several important sparse tensor kernels, including sparse matrix-matrix multiplication (SpMM), sparse tensor addition (SpAdd), and the matricized tensor times Khatri-Rao product (MTTKRP) used to factorize tensors. Our results show improvements over prior work on tensor algebra compilation and brings the performance of these kernels on par with state-of-the-art hand-optimized implementations. For example, SpMM was not supported by prior tensor algebra compilers, the performance of MTTKRP on the nell-2 data set improves by 35%, and MTTKRP can for the first time have sparse results.


title: Near-Optimal Sample Complexity Bounds for Maximum Likelihood Estimation of Multivariate Log-concave Densities

abstract: We study the problem of learning multivariate log-concave densities with respect to a global loss function. We obtain the first upper bound on the sample complexity of the maximum likelihood estimator (MLE) for a log-concave density on $\mathbb{R}^d$, for all $d \geq 4$. Prior to this work, no finite sample upper bound was known for this estimator in more than $3$ dimensions. In more detail, we prove that for any $d \geq 1$ and $\epsilon>0$, given $\tilde{O}_d((1/\epsilon)^{(d+3)/2})$ samples drawn from an unknown log-concave density $f_0$ on $\mathbb{R}^d$, the MLE outputs a hypothesis $h$ that with high probability is $\epsilon$-close to $f_0$, in squared Hellinger loss. A sample complexity lower bound of $\Omega_d((1/\epsilon)^{(d+1)/2})$ was previously known for any learning algorithm that achieves this guarantee. We thus establish that the sample complexity of the log-concave MLE is near-optimal, up to an $\tilde{O}(1/\epsilon)$ factor.


title: Modeling Activity Tracker Data Using Deep Boltzmann Machines

abstract: Commercial activity trackers are set to become an essential tool in health research, due to increasing availability in the general population. The corresponding vast amounts of mostly unlabeled data pose a challenge to statistical modeling approaches. To investigate the feasibility of deep learning approaches for unsupervised learning with such data, we examine weekly usage patterns of Fitbit activity trackers with deep Boltzmann machines (DBMs). This method is particularly suitable for modeling complex joint distributions via latent variables. We also chose this specific procedure because it is a generative approach, i.e., artificial samples can be generated to explore the learned structure. We describe how the data can be preprocessed to be compatible with binary DBMs. The results reveal two distinct usage patterns in which one group frequently uses trackers on Mondays and Tuesdays, whereas the other uses trackers during the entire week. This exemplary result shows that DBMs are feasible and can be useful for modeling activity tracker data.


title: Constraining Baryon--Dark Matter Scattering with the Cosmic Dawn 21-cm Signal

abstract: The recent detection of an anomalously strong 21-cm signal of neutral hydrogen from Cosmic Dawn by the EDGES Low-Band radio experiment can be explained if cold dark matter particles scattered off the baryons draining excess energy from the gas. In this Letter we explore the expanded range of the 21-cm signal that is opened up by this interaction, varying the astrophysical parameters as well as the properties of dark matter particles in the widest possible range. We identify models consistent with current data by comparing to both the detection in the Low-Band and the upper limits from the EDGES High-Band antenna. We find that consistent models predict a 21-cm fluctuation during Cosmic Dawn that is between 3 and 30 times larger than the largest previously expected without dark matter scattering. The expected power spectrum exhibits strong Baryon Acoustic Oscillations imprinted by the velocity-dependent cross-section. The latter signature is a smoking gun of the velocity-dependent scattering and could be used by interferometers to verify the dark matter explanation of the EDGES detection.


title: Rings of constants of linear derivations on Fermat rings

abstract: In this paper we characterize all the linear $\mathbb{C}$-derivations of the Fermat ring. We show that the Fermat ring has linear $\mathbb{C}$-derivations with trivial ring of constants and construct some examples.


title: Renormalization group analysis for the quasi-1D superconductor BaFe$_2$S$_3$

abstract: Motivated by the discovery of superconductivity in the two-leg, quasi-one dimensional ladder compound, BaFe$_2$S$_3$ we present a renormalization group study of electrons moving on a two leg, two orbital ladder, subjected to Hubbard repulsion $U$ and Hund's coupling $J$. In our calculations, we adopt tightbinding parameters obtained from ab-initio studies on this material. At incommensurate filling, the long wavelength analysis displays four phases as a function of $0 \leq J/U <1$. We show that a fully gapped superconductor is stabilized at sufficiently large Hund's coupling, the relative phases at the three Fermi points are "+,-,-". By contrast, when the system is tuned to half filling, Umklapp scattering gives rise to Mott insulating phases. We discuss the general implications of our study for the broad class of iron-based superconductors.


title: The Coulomb problem in iron based superconductors

abstract: We discuss the role of strong Coulomb interactions in iron-based superconductors (FeSCs). The presumed $s^{\pm}$ character of these superconductors means that the condensate is not symmetry protected against Coulomb repulsion. Remarkably, the transition temperatures and the excitation gap are quite robust across the large family of iron based superconductors, despite drastic changes in Fermi surface geometry. The Coulomb problem is to understand how these superconductors avoid the strong onsite Coulomb interaction at the iron atoms, while maintaining a robust transition temperature. Within the dominant space of $t_{2g}$ orbitals, on-site repulsion in the FeSCs enforces two linearly independent components of the condensate to vanish. This raises the possibility that iron-based superconductors might adapt their condensate to the Coulomb constraints by rotating the pairing state within the large manifold of entangled, extended s-wave gap functions with different orbital and momentum space structure. We examine this "orbital and k-space flexibility" (OKF) mechanism using both Landau theory and microscopic calculations within a multi-orbital t-J model. Based on our results, we conclude that OKF necessitates a large condensate degeneracy. One interesting possibility raised by our results, is that a resolution to the Coulomb problem in FeSC might require a reconsideration of triplet pairing.


title: Orbifolds of Lattice Vertex Operator Algebras at $d=48$ and $d=72$

abstract: Motivated by the notion of extremal vertex operator algebras, we investigate cyclic orbifolds of vertex operator algebras coming from extremal even self-dual lattices in $d=48$ and $d=72$. In this way we construct about one hundred new examples of holomorphic VOAs with a small number of low weight states.