title: Asymmetric Electric Field Screening in van der Waals Heterostructures

abstract: Electric field screening plays an important role in the physical and chemical properties of materials and their devices. Here, we use a compelling set of theoretical and experimental techniques involving van der Waals (vdW) ab initio density functional theory (DFT) simulations, quantum capacitance-based classical model and electric force microscopy (EFM) to elucidate the intrinsic dielectric screening properties of vdW heterostructures (vdWHs) formed by MoS2 and graphene layers. We experimentally observed an asymmetric electric response in the MoS2/Graphene vdWHs under different directions of the external electric field. That is, when the electric fields are shed towards graphene, a large amount of polarized charges screen the fields, but as the sign of the field was reversed, a strong depolarization field was present, and a partial screening was detected. This effect is thickness-dependent, in particular on the number of the MoS2 layers; whereas increased thickness of graphene showed a small effect on their electrical and screening behavior. Our results indicate that asymmetric dipolar contributions at the interface between graphene and MoS2 are the main cause to the unusual field-effect screening in the vdWHs. This work not only provides new insights on the screening properties of a vast amount of heterojunction fabricated so far, but also uncovers the great potential of controlling a fundamental property, such as screening, for device applications.


title: Exact Solution of a Strongly Coupled Gauge Theory in 0+1 Dimensions

abstract: Gauged tensor models are a class of strongly coupled quantum mechanical theories. We present the exact analytic solution of a specific example of such a theory: namely the smallest colored tensor model due to Gurau and Witten that exhibits non-linearities. We find explicit analytic expressions for the eigenvalues and eigenstates, and the former agree precisely with previous numerical results on (a subset of) eigenvalues of the ungauged theory. The physics of the spectrum, despite the smallness of $N$, exhibits rudimentary signatures of chaos. This Letter is a summary of our main results: the gory details will appear in a companion paper.


title: A Praise for Defensive Programming: Leveraging Uncertainty for Effective Malware Mitigation

abstract: A promising avenue for improving the effectiveness of behavioral-based malware detectors would be to combine fast traditional machine learning detectors with high-accuracy, but time-consuming deep learning models. The main idea would be to place software receiving borderline classifications by traditional machine learning methods in an environment where uncertainty is added, while software is analyzed by more time-consuming deep learning models. The goal of uncertainty would be to rate-limit actions of potential malware during the time consuming deep analysis. In this paper, we present a detailed description of the analysis and implementation of CHAMELEON, a framework for realizing this uncertain environment for Linux. CHAMELEON offers two environments for software: (i) standard - for any software identified as benign by conventional machine learning methods and (ii) uncertain - for software receiving borderline classifications when analyzed by these conventional machine learning methods. The uncertain environment adds obstacles to software execution through random perturbations applied probabilistically on selected system calls. We evaluated CHAMELEON with 113 applications and 100 malware samples for Linux. Our results showed that at threshold 10%, intrusive and non-intrusive strategies caused approximately 65% of malware to fail accomplishing their tasks, while approximately 30% of the analyzed benign software to meet with various levels of disruption. With a dynamic, per-system call threshold, CHAMELEON caused 92% of the malware to fail, and only 10% of the benign software to be disrupted. We also found that I/O-bound software was three times more affected by uncertainty than CPU-bound software. Further, we analyzed the logs of software crashed with non-intrusive strategies, and found that some crashes are due to the software bugs.


title: The revolution in physics of the early Nineteenth century revisited in the context of science-and-society interaction

abstract: The radical changes in the concepts and approach in Physics at the turn of the Nineteenth century were so deep, that is acknowledged as a revolution. However, in 1970 Thomas Kuhn's careful reconstruction of the researches on the black body problem, the concept itself of the revolution seemed to vanish in his diluted discussion of every details. In the present paper, after an examination of the limitations of Kuhn's response to his critics, I put forward the idea, although it is not new, that these changes in Physics cannot be reduced to a point-like event, but happened instead through multiple successive (and even contradictory) changes in the course of decades. Such as the old quantum hypothesis, wave mechanics, orthodox quantum mechanics. In fact, the innovative perspectives started in the 1980s have been considered as a third quantum revolution. My basic argument is that these changes, in order to be really understood, must be interpreted not as mere specific changes in Physics, but framed in the context of the deep social, cultural, and economic changes during those turbulent years. The main steps are outlined.


title: The monodromy of meromorphic projective structures

abstract: We study projective structures on a surface having poles of prescribed orders. We obtain a monodromy map from a complex manifold parameterising such structures to the stack of framed $\mathrm{PGL}_2(\mathbb{C})$ local systems on the associated marked bordered surface. We prove that the image of this map is contained in the union of the domains of the cluster charts. We discuss a number of open questions concerning this monodromy map.


title: Can Local Stress Enhancement Induce Stability in Fracture Processes? Part I: Apparent Stability

abstract: By comparing the evolution of the local and equal load sharing fiber bundle models, we point out the paradoxical result that stresses seem to make the local load sharing model stable when the equal load sharing model is not. We explain this behavior by demonstrating that it is only an apparent stability in the local load sharing model, which originates from a statistical effect due to sample averaging. Even though we use the fiber bundle model to demonstrate the apparent stability, we argue that it is a more general feature of fracture processes.


title: Measuring third party tracker power across web and mobile

abstract: Third-party networks collect vast amounts of data about users via web sites and mobile applications. Consolidations among tracker companies can significantly increase their individual tracking capabilities, prompting scrutiny by competition regulators. Traditional measures of market share, based on revenue or sales, fail to represent the tracking capability of a tracker, especially if it spans both web and mobile. This paper proposes a new approach to measure the concentration of tracking capability, based on the reach of a tracker on popular websites and apps. Our results reveal that tracker prominence and parent-subsidiary relationships have significant impact on accurately measuring concentration.


title: Spectral Phase Control of Interfering Chirped Pulses for High-Energy Narrowband Terahertz Generation

abstract: Highly-efficient optical generation of narrowband terahertz (THz) radiation enables unexplored technologies and sciences from compact electron acceleration to charge manipulation in solids. State-of-the-art conversion efficiencies are currently achieved using difference-frequency generation (DFG) driven by temporal beating of chirped pulses but remain, however, far lower than desired or predicted. Here we show that high-order spectral phase fundamentally limits the efficiency of narrowband DFG using chirped-pulse beating and resolve this limitation by introducing a novel technique based on tuning the relative spectral phase of the pulses. For optical terahertz generation, we demonstrate a 13-fold enhancement in conversion efficiency for 1%-bandwidth, 0.361 THz pulses, yielding a record energy of 0.6 mJ and exceeding previous optically-generated energies by over an order of magnitude. Our results prove the feasibility of millijoule-scale applications like terahertz-based electron accelerators and light sources and solve the long-standing problem of temporal irregularities in the pulse trains generated by interfering chirped pulses.


title: Strong Amplifiers of Natural Selection: Proofs

abstract: We consider the modified Moran process on graphs to study the spread of genetic and cultural mutations on structured populations. An initial mutant arises either spontaneously (aka \emph{uniform initialization}), or during reproduction (aka \emph{temperature initialization}) in a population of $n$ individuals, and has a fixed fitness advantage $r>1$ over the residents of the population. The fixation probability is the probability that the mutant takes over the entire population. Graphs that ensure fixation probability of~1 in the limit of infinite populations are called \emph{strong amplifiers}. Previously, only a few examples of strong amplifiers were known for uniform initialization, whereas no strong amplifiers were known for temperature initialization. In this work, we study necessary and sufficient conditions for strong amplification, and prove negative and positive results. We show that for temperature initialization, graphs that are unweighted and/or self-loop-free have fixation probability upper-bounded by $1-1/f(r)$, where $f(r)$ is a function linear in $r$. Similarly, we show that for uniform initialization, bounded-degree graphs that are unweighted and/or self-loop-free have fixation probability upper-bounded by $1-1/g(r,c)$, where $c$ is the degree bound and $g(r,c)$ a function linear in $r$. Our main positive result complements these negative results, and is as follows: every family of undirected graphs with (i)~self loops and (ii)~diameter bounded by $n^{1-\epsilon}$, for some fixed $\epsilon>0$, can be assigned weights that makes it a strong amplifier, both for uniform and temperature initialization.


title: Control of surface potential at polar domain walls in a nonpolar oxide

abstract: Ferroic domain walls could play an important role in microelectronics, given their nanometric size and often distinct functional properties. Until now, devices and device concepts were mostly based on mobile domain walls in ferromagnetic and ferroelectric materials. A less explored path is to make use of polar domain walls in nonpolar ferroelastic materials. Indeed, while the polar character of ferroelastic domain walls has been demonstrated, polarization control has been elusive. Here, we report evidence for the electrostatic signature of the domain-wall polarization in nonpolar calcium titanate (CaTiO3). Macroscopic mechanical resonances excited by an ac electric field are observed as a signature of a piezoelectric response caused by polar walls. On the microscopic scale, the polarization in domain walls modifies the local surface potential of the sample. Through imaging of surface potential variations, we show that the potential at the domain wall can be controlled by electron injection. This could enable devices based on nondestructive information readout of surface potential.


title: DeepHeart: Semi-Supervised Sequence Learning for Cardiovascular Risk Prediction

abstract: We train and validate a semi-supervised, multi-task LSTM on 57,675 person-weeks of data from off-the-shelf wearable heart rate sensors, showing high accuracy at detecting multiple medical conditions, including diabetes (0.8451), high cholesterol (0.7441), high blood pressure (0.8086), and sleep apnea (0.8298). We compare two semi-supervised train- ing methods, semi-supervised sequence learning and heuristic pretraining, and show they outperform hand-engineered biomarkers from the medical literature. We believe our work suggests a new approach to patient risk stratification based on cardiovascular risk scores derived from popular wearables such as Fitbit, Apple Watch, or Android Wear.


title: Large deviations of reaction fluxes

abstract: We study a system of interacting particles that randomly react to form new particles. The reaction flux is the rescaled number of reactions that take place in a time interval. We prove a dynamic large-deviation principle for the reaction fluxes under general assumptions that include mass-action kinetics. This result immediately implies the dynamic large deviations for the empirical concentration.


title: Fra\"iss\'e structures and a conjecture of Furstenberg

abstract: We study problems concerning the Samuel compactification of the automorphism group of a countable first-order structure. A key motivating question is a problem of Furstenberg and a counter-conjecture by Pestov regarding the difference between $S(G)$, the Samuel compactification, and $E(M(G))$, the enveloping semigroup of the universal minimal flow. We resolve Furstenberg's problem for several automorphism groups and give a detailed study in the case of $G = S_\infty$, leading us to define and investigate several new types of ultrafilter on a countable set.


title: B\"uchi-Kamp Theorems for 1-clock ATA

abstract: This paper investigates Kamp-like and B\"uchi-like theorems for 1-clock Alternating Timed Automata (1-ATA) and its natural subclasses. A notion of 1-ATA with loop-free-resets is defined. This automaton class is shown to be expressively equivalent to the temporal logic $\regmtl$ which is $\mathsf{MTL[F_I]}$ extended with a regular expression guarded modality. Moreover, a subclass of future timed MSO with k-variable-connectivity property is introduced as logic $\qkmso$. In a Kamp-like result, it is shown that $\regmtl$ is expressively equivalent to $\qkmso$. As our second result, we define a notion of conjunctive-disjunctive 1-clock ATA ($\wf$ 1-ATA). We show that $\wf$ 1-ATA with loop-free-resets are expressively equivalent to the sublogic $\F\regmtl$ of $\regmtl$. Moreover $\F\regmtl$ is expressively equivalent to $\qtwomso$, the two-variable connected fragment of $\qkmso$. The full class of 1-ATA is shown to be expressively equivalent to $\regmtl$ extended with fixed point operators.


title: Design study of random spectrometers for applications at optical frequencies

abstract: Compact spectrometers based on disordered planar waveguides exhibit a rather high resolution with a relatively small footprint as compared to conventional spectrometers. This is achieved by multiple scattering of light which - if properly engineered - significantly enhances the effective optical path length. Here a design study of random spectrometers for TE- and TM-polarized light is presented that combines the results of Mie theory, multiple-scattering theory and full electromagnetic simulations. It is shown that the performance of such random spectrometers depends on single scattering quantities, notably on the overall scattering efficiency and the asymmetry parameter. Further, the study shows that a well-developed diffusive regime is not required in practice and that a standard integrated-optical layout is sufficient to obtain efficient devices even for rather weakly scattering systems consisting of low index inclusions in high-index matrices such as pores in planar silicon-nitride based waveguides. This allows for both significant reductions in footprint with acceptable losses in resolution and for device operation in the visible and near-infrared frequency range.


title: Towards Understanding the Structure, Dynamics and Bio-activity of Diabetic Drug Metformin

abstract: Small molecules are often found to exhibit extraordinarily diverse biological activities. Metformin is one of them. It is widely used as anti-diabetic drug for type-two diabetes. In addition to that, metformin hydrochloride shows anti-tumour activities and increases the survival rate of patients suffering from certain types of cancer namely colorectal, breast, pancreas and prostate cancer. However, theoretical studies of structure and dynamics of metformin have not yet been fully explored. In this work, we investigate the characteristic structural and dynamical features of three mono-protonated forms of metformin hydrochloride with the help of experiments, quantum chemical calculations and atomistic molecular dynamics simulations. We validate our force field by comparing simulation results to that of the experimental findings. Nevertheless, we discover that the non-planar tautomeric form is the most stable. Metformin forms strong hydrogen bonds with surrounding water molecules and its solvation dynamics show unique features. Because of an extended positive charge distribution, metformin possesses features of being a permanent cationic partner toward several targets. We study its interaction and binding ability with DNA using UV spectroscopy, circular dichroism, fluorimetry and metadynamics simulation. We find a non-intercalating mode of interaction. Metformin feasibly forms a minor/major groove-bound state within a few tens of nanoseconds, preferably with AT rich domains. A significant decrease in the free-energy of binding is observed when it binds to a minor groove of DNA.


title: Blueprint and Evaluation Instruments for a Course on Software Engineering for Sustainability

abstract: We report on a summer school course on Software Engineering for Sustainability (SE4S). We provide a detailed blueprint of the contents taught and its evaluation with the instruments that were used.


title: Disconnection by level sets of the discrete Gaussian free field and entropic repulsion

abstract: We derive asymptotic upper and lower bounds on the large deviation probability that the level set of the Gaussian free field on $Z^d$, d bigger or equal to three, below a given level disconnects the discrete blow-up of a compact set A from the boundary of the discrete blow-up of a box that contains A, when the level set of the Gaussian free field above this level is in a strongly percolative regime. These bounds substantially strengthen the results of arXiv:1412.3960, where A was a box and the convexity of A played an important role in the proof. We also derive an asymptotic upper bound on the probability that the average of the Gaussian free field well inside the discrete blow-up of A is above a certain level when disconnection occurs. The derivation of the upper bounds uses the solidification estimates for porous interfaces that were derived in the work arXiv:1706.07229 of A.-S. Sznitman and the author to treat a similar disconnection problem for the vacant set of random interlacements. If certain critical levels for the Gaussian free field coincide, an open question at the moment, the asymptotic upper and lower bounds that we obtain for the disconnection probability match in principal order, and conditioning on disconnection lowers the average of the Gaussian free field well inside the discrete blow-up of A, which can be understood as entropic repulsion.


title: The Ising model coupled to 2d orders

abstract: In this article we make first steps in coupling matter to causal set theory in the path integral. We explore the case of the Ising model coupled to the 2d discrete Einstein Hilbert action, restricted to the 2d orders. We probe the phase diagram in terms of the Wick rotation parameter \b{eta} and the Ising coupling j and find that the matter and the causal sets together give rise to an interesting phase structure. The couplings give rise to five different phases. The causal sets take on random or crystalline characteristics as described in [1] and the Ising model can be correlated or uncorrelated on the random ordes and correlated, uncorrelated or anti-correlated on the crystalline orders. We find that at least one new phase transition arises, in which the Ising spins push the causal set into the crystalline phase.


title: Combinatorial views on persistent characters in phylogenetics

abstract: The so-called binary perfect phylogeny with persistent characters has recently been thoroughly studied in computational biology as it is less restrictive than the well known binary perfect phylogeny. Here, we focus on the notion of (binary) persistent characters, i.e. characters that can be realized on a phylogenetic tree by at most one $0 \rightarrow 1$ transition followed by at most one $1 \rightarrow 0$ transition in the tree, and analyze these characters under different aspects. First, we illustrate the connection between persistent characters and Maximum Parsimony, where we characterize persistent characters in terms of the first phase of the famous Fitch algorithm. Afterwards we focus on the number of persistent characters for a given phylogenetic tree. We show that this number solely depends on the balance of the tree. To be precise, we develop a formula for counting the number of persistent characters for a given phylogenetic tree based on an index of tree balance, namely the Sackin index. Lastly, we consider the question of how many (carefully chosen) binary characters together with their persistence status are needed to uniquely determine a phylogenetic tree and provide an upper bound for the number of characters needed.


title: Infinite series in cohomology: attractors and Conley index

abstract: In this paper we study the cohomological Conley index of arbitrary isolated invariant continua for continuous maps $f \colon U \subseteq \mathbb{R}^d \to \mathbb{R}^d$ by analyzing the topological structure of their unstable manifold. We provide a simple dynamical interpretation for the first cohomological Conley index, describing it completely, and relate it to the cohomological Conley index in higher degrees. A number of consequences are derived, including new computations of the fixed point indices of isolated invariant continua in dimensions 2 and 3. Our approach exploits certain attractor-repeller decomposition of the unstable manifold, reducing the study of the cohomological Conley index to the relation between the cohomology of an attractor and its basin of attraction. This is a classical problem that, in the present case, is particularly difficult because the dynamics is discrete and the topology of the unstable manifold can be very complicated. To address it we develop a new method that may be of independent interest and involves the summation of power series in cohomology: if $Z$ is a metric space and $K \subseteq Z$ is a compact, global attractor for a continuous map $g \colon Z \to Z$, we show how to interpret series of the form $\sum_{j \ge 0} a_j (g^*)^j$ as endomorphisms of the cohomology group of the pair $(Z,K)$.


title: Joint Attention in Driver-Pedestrian Interaction: from Theory to Practice

abstract: Today, one of the major challenges that autonomous vehicles are facing is the ability to drive in urban environments. Such a task requires communication between autonomous vehicles and other road users in order to resolve various traffic ambiguities. The interaction between road users is a form of negotiation in which the parties involved have to share their attention regarding a common objective or a goal (e.g. crossing an intersection), and coordinate their actions in order to accomplish it. In this literature review we aim to address the interaction problem between pedestrians and drivers (or vehicles) from joint attention point of view. More specifically, we will discuss the theoretical background behind joint attention, its application to traffic interaction and practical approaches to implementing joint attention for autonomous vehicles.


title: Plasma Brain Dynamics (PBD): A Mechanism for EEG Waves Under Human Consciousness

abstract: EEG signals are records of nonlinear solitary waves in human brains. The waves have several types (e.g., a, b, g, q, d) in response to different levels of consciousness. They are classified into two groups: Group-1 consists of complex storm-like waves (a, b, and g); Group-2 is composed of simple quasilinear waves (q and d). In order to elucidate the mechanism of EEG wave formation and propagation, this paper extends the Vlasov-Maxwell equations of Plasma Brain Dynamics (PBD) to a set of two-fluid, self-similar, nonlinear solitary wave equations. Numerical simulations are performed for different EEG signals. Main results include: (1) The excitation and propagation of the EEG wave packets are dependent of electric and magnetic fields, brain aqua-ions, electron and ion temperatures, masses, and their initial fluid speeds; (2) Group-1 complex waves contain three ingredients: the high-frequency ion-acoustic (IA) mode, the intermediate-frequency lower-hybrid (LH) mode, and the low-frequency ion-cyclotron (IC) mode; (3) Group-2 simple waves fall within the IA band, featured by one or a combination of the three envelopes: sinusoidal, sawtooth, and spiky/bipolar. The study proposes an alternative model to Quantum Brain Dynamics (QBD) by suggesting that the formation and propagation of the nonlinear solitary EEG waves in the brain have the same mechanism as that of the waves in space plasmas


title: Subleading Regge limit from a soft anomalous dimension

abstract: Wilson lines capture important features of scattering amplitudes, for example soft effects relevant for infrared divergences, and the Regge limit. Beyond the leading power approximation, corrections to the eikonal picture have to be taken into account. In this paper, we study such corrections in a model of massive scattering amplitudes in N = 4 super Yang-Mills, in the planar limit, where the mass is generated through a Higgs mechanism. Using known three-loop analytic expressions for the scattering amplitude, we find that the first power suppressed term has a very simple form, equal to a single power law. We propose that its exponent is governed by the anomalous dimension of a Wilson loop with a scalar inserted at the cusp, and we provide perturbative evidence for this proposal. We also analyze other limits of the amplitude and conjecture an exact formula for a total cross-section at high energies.


title: Spin-charge separation and many-body localization

abstract: We study many-body localization for a disordered chain of spin 1/2 fermions. In [Phys. Rev. B \textbf{94}, 241104 (2016)], when both down and up components are exposed to the same strong disorder, the authors observe a power law growth of the entanglement entropy that suggests that many-body localization is not complete; the density (charge) degree of freedom is localized, while the spin degree of freedom is apparently delocalized. We show that this power-like behavior is only a transient effect and that, for longer times, the growth is logarithmic in time suggesting that the spin degree of freedom is also localized, so that the system follows the standard many-body localization scenario. We also study the experimentally relevant case of quasiperiodic disorder.


title: Detecting false correlations: Uncovering a faked Bell-inequality violation

abstract: It is possible for two parties, Alice and Bob, to establish a secure communication link by sharing an ensemble of entangled particles, and then using these particles to generate a secret key. One way to establish that the particles are indeed entangled is to verify that they violate a Bell inequality. However, it might be the case that Bob is not trustworthy and wishes Alice to believe that their communications are secure, when in fact they are not. He can do this by managing to have prior knowledge of Alice's measurement device settings and then modifying his own settings based upon this information. In this case it is possible for shared particle states that must satisfy a Bell inequality to appear to violate this inequality, which would also make the system appear secure. When Bob modifies his measurement settings, however, he produces false correlations. Here we demonstrate experimentally that Alice can detect these false correlations, and uncover Bob's trickery, by using loop-state-preparation-and-measurement (SPAM) tomography. More generally, we demonstrate that loop SPAM tomography can detect false correlations (correlated errors) in a two-qubit system without needing to know anything about the prepared states or the measurements, other than the dimensions of the operators that describe them.


title: Polariton Hall Effect in Transition-Metal Dichalcogenides

abstract: We analyze the properties of strongly coupled excitons and photons in systems made of semiconducting two-dimensional transition-metal dichalcogenides embedded in optical cavities. Through a detailed microscopic analysis of the coupling we unveil novel, highly tunable features of the spectrum, that result in polariton splitting and a breaking of light-matter selection rules. The dynamics of the composite polaritons is influenced by the Berry phase arising both from their constituents and from the confinement-enhanced coupling. We find that light-matter coupling emerges as a mechanism that enhances the Berry phase of polaritons well beyond that of its elementary constituents, paving the way to achieve a polariton Hall effect.


title: Classification of Things in DBpedia using Deep Neural Networks

abstract: The Semantic Web aims at representing knowledge about the real world at web scale - things, their attributes and relationships among them can be represented as nodes and edges in an inter-linked semantic graph. In the presence of noisy data, as is typical of data on the Semantic Web, a software Agent needs to be able to robustly infer one or more associated actionable classes for the individuals in order to act automatically on it. We model this problem as a multi-label classification task where we want to robustly identify types of the individuals in a semantic graph such as DBpedia, which we use as an exemplary dataset on the Semantic Web. Our approach first extracts multiple features for the individuals using random walks and then performs multi-label classification using fully-connected Neural Networks. Through systematic exploration and experimentation, we identify the effect of hyper-parameters of the feature extraction and the fully-connected Neural Network structure on the classification performance. Our final results show that our method performs better than state-of-the-art inferencing systems like SDtype and SLCN, from which we can conclude that random-walk-based feature extraction of individuals and their multi-label classification using Deep Neural Networks is a promising alternative to these systems for type classification of individuals on the Semantic Web. The main contribution of our work is to introduce a novel approach that allows us to use Deep Neural Networks to identify types of individuals in a noisy semantic graph by extracting features using random walks


title: Comment on 'Simulation of ultra-relativistic electrons and positrons channeling in crystals with MBN Explorer'

abstract: The snapshot model of crystal atoms was implemented in the Monte Carlo code ChaS (Channeling Simulator) and is being successfully used for simulation of ultrarelativistic particle channeling. The model was criticized by Sushko et al. (J. Comp. Phys. 252 (2013) 404-418) who claim that it overestimates the mean scattering angle in a single projectile-atom collision. As a matter of fact, no evidence that would support this claim can be found in the mentioned publication. Moreover, the snapshot model and the model suggested by Sushko et al. yield essentially the same value of the mean scattering angle. Contrary to the claim of Sushko et al., the target electrons can be considered as fixed-position scatterer, corrections due to their finite mass and nonzero initial velocity have a negligible impact on the channeling of light projectiles (electrons and positrons). In contrast to the snapshot model, the model preferred by Sushko et al. does not take into account incoherent scattering of the projectile by crystal electrons. This explains why the two models predict different values of the dechanneling length. The claim that the snapshot model underestimates the dechanneling length is unfounded. In actual fact, this model is in good agreement with experimental data.


title: Magnetic Isotope Effect in the Uranium Isotopes Separation

abstract: When uranyl nitrate is photolysed in the water solution, light isotope ${}^{235}{\rm U}$ is separated from isotope ${}^{238}{\rm U}$ ; an enrichment factor $K = 1.04$ . The initial samples were depleted by a light isotope of uranium ${}^{235}{\rm U}$ (degree $\beta = 0.00455$). Regenerated uranyl nitrate was enriched by an isotope of uranium ${}^{235}{\rm U}$ (degree $\alpha = 0.00464$). Uranium tetrafluoride was depleted by an isotope of uranium ${}^{235}{\rm U}$ (degree $\beta = 0.00446$). Uranium isotopes are separated by the magnetic isotope effect (MIE) in uranyl nitrate photoreduction.


title: Fair comparison of skin detection approaches on publicly available datasets

abstract: Skin detection is the process of discriminating skin and non-skin regions in a digital image and it is widely used in several applications ranging from hand gesture analysis to track body parts and face detection. Skin detection is a challenging problem which has drawn extensive attention from the research community, nevertheless a fair comparison among approaches is very difficult due to the lack of a common benchmark and a unified testing protocol. In this work, we investigate the most recent researches in this field and we propose a fair comparison among approaches using several different datasets. The major contributions of this work are an exhaustive literature review of skin color detection approaches, a framework to evaluate and combine different skin detector approaches, whose source code is made freely available for future research, and an extensive experimental comparison among several recent methods which have also been used to define an ensemble that works well in many different problems. Experiments are carried out in 10 different datasets including more than 10000 labelled images: experimental results confirm that the best method here proposed obtains a very good performance with respect to other stand-alone approaches, without requiring ad hoc parameter tuning. A MATLAB version of the framework for testing and of the methods proposed in this paper will be freely available from https://github.com/LorisNanni


title: A Spatial Mapping Algorithm with Applications in Deep Learning-Based Structure Classification

abstract: Convolutional Neural Network (CNN)-based machine learning systems have made breakthroughs in feature extraction and image recognition tasks in two dimensions (2D). Although there is significant ongoing work to apply CNN technology to domains involving complex 3D data, the success of such efforts has been constrained, in part, by limitations in data representation techniques. Most current approaches rely upon low-resolution 3D models, strategic limitation of scope in the 3D space, or the application of lossy projection techniques to allow for the use of 2D CNNs. To address this issue, we present a mapping algorithm that converts 3D structures to 2D and 1D data grids by mapping a traversal of a 3D space-filling curve to the traversal of corresponding 2D and 1D curves. We explore the performance of 2D and 1D CNNs trained on data encoded with our method versus comparable volumetric CNNs operating upon raw 3D data from a popular benchmarking dataset. Our experiments demonstrate that both 2D and 1D representations of 3D data generated via our method preserve a significant proportion of the 3D data's features in forms learnable by CNNs. Furthermore, we demonstrate that our method of encoding 3D data into lower-dimensional representations allows for decreased CNN training time cost, increased original 3D model rendering resolutions, and supports increased numbers of data channels when compared to purely volumetric approaches. This demonstration is accomplished in the context of a structural biology classification task wherein we train 3D, 2D, and 1D CNNs on examples of two homologous branches within the Ras protein family. The essential contribution of this paper is the introduction of a dimensionality-reduction method that may ease the application of powerful deep learning tools to domains characterized by complex structural data.


title: Factors of generalised polynomials and automatic sequences

abstract: The aim of this short note is to generalise the result of Rampersad--Shallit saying that an automatic sequence and a Sturmian sequence cannot have arbitrarily long common factors. We show that the same result holds if a Sturmian sequence is replaced by an arbitrary sequence whose terms are given by a generalised polynomial (i.e., an expression involving algebraic operations and the floor function) that is not periodic except for a set of density zero.


title: FixaTons: A collection of Human Fixations Datasets and Metrics for Scanpath Similarity

abstract: In the last three decades, human visual attention has been a topic of great interest in various disciplines. In computer vision, many models have been proposed to predict the distribution of human fixations on a visual stimulus. Recently, thanks to the creation of large collections of data, machine learning algorithms have obtained state-of-the-art performance on the task of saliency map estimation. On the other hand, computational models of scanpath are much less studied. Works are often only descriptive or task specific. This is due to the fact that the scanpath is harder to model because it must include the description of a dynamic. General purpose computational models are present in the literature, but are then evaluated in tasks of saliency prediction, losing therefore information about the dynamics and the behaviour. In addition, two technical reasons have limited the research. The first reason is the lack of robust and uniformly used set of metrics to compare the similarity between scanpath. The second reason is the lack of sufficiently large and varied scanpath datasets. In this report we want to help in both directions. We present FixaTons, a large collection of datasets human scanpaths (temporally ordered sequences of fixations) and saliency maps. It comes along with a software library for easy data usage, statistics calculation and implementation of metrics for scanpath and saliency prediction evaluation.


title: Directly and Efficiently Optimizing Prediction Error and AUC of Linear Classifiers

abstract: The predictive quality of machine learning models is typically measured in terms of their (approximate) expected prediction error or the so-called Area Under the Curve (AUC) for a particular data distribution. However, when the models are constructed by the means of empirical risk minimization, surrogate functions such as the logistic loss are optimized instead. This is done because the empirical approximations of the expected error and AUC functions are nonconvex and nonsmooth, and more importantly have zero derivative almost everywhere. In this work, we show that in the case of linear predictors, and under the assumption that the data has normal distribution, the expected error and the expected AUC are not only smooth, but have closed form expressions, which depend on the first and second moments of the normal distribution. Hence, we derive derivatives of these two functions and use these derivatives in an optimization algorithm to directly optimize the expected error and the AUC. In the case of real data sets, the derivatives can be approximated using empirical moments. We show that even when data is not normally distributed, computed derivatives are sufficiently useful to render an efficient optimization method and high quality solutions. Thus, we propose a gradient-based optimization method for direct optimization of the prediction error and AUC. Moreover, the per-iteration complexity of the proposed algorithm has no dependence on the size of the data set, unlike those for optimizing logistic regression and all other well known empirical risk minimization problems.


title: Low-Frequency Noise and Sliding of the Charge Density Waves in Two-Dimensional Materials

abstract: There has been a recent renewal of interest in charge-density-wave (CDW) phenomena, primarily driven by the emergence of two-dimensional (2D) layered CDW materials, such as 1T-TaS2, characterized by very high transition temperatures to CDW phases. In the extensively studied classical bulk CDW materials with quasi-1D crystal structure, the charge carrier transport exhibits intriguing sliding behavior, which reveals itself in the frequency domain as "narrowband" and "broadband" noise. Despite the increasing attention on physics of 2D CDWs, there have been few reports of CDW sliding, specifically in quasi-2D rare-earth tritellurides and none on the noise in any of 2D CDW systems. Here we report the results of low-frequency noise (LFN) measurements on 1T-TaS2 thin films - archetypal 2D CDW systems, as they are driven from the nearly commensurate (NC) to incommensurate (IC) CDW phases by voltage and temperature stimuli. We have found that noise in 1T-TaS2 devices has two pronounced maxima at the bias voltages, which correspond to the onset of CDW sliding and the NC-to-IC phase transition. We observed unusual Lorentzian noise features and exceptionally strong noise dependence on electric bias and temperature. We argue that LFN in 2D CDW systems has unique physical origin, different from known fundamental noise types. The specifics of LFN in 2D CDW materials can be explained by invoking the concept of interacting discrete fluctuators in the NC-CDW phase. Noise spectroscopy can serve as a useful tool for understanding electronic transport phenomena in 2D CDW materials characterized by coexistence of different phases and strong CDW pinning.


title: Minimally modified theories of gravity: a playground for testing the uniqueness of general relativity

abstract: In a recent paper [1], it was introduced a new class of gravitational theories with two local degrees of freedom. The existence of these theories apparently challenges the distinctive role of general relativity as the unique non-linear theory of massless spin-2 particles. Here we perform a comprehensive analysis of these theories with the aim of (i) understanding whether or not these are actually equivalent to general relativity, and (ii) finding the root of the variance in case these are not. We have found that a broad set of seemingly different theories actually pass all the possible tests of equivalence to general relativity (in vacuum) that we were able to devise, including the analysis of scattering amplitudes using on-shell techniques. These results are complemented with the observation that the only examples which are manifestly not equivalent to general relativity either do not contain gravitons in their spectrum, or are not guaranteed to include only two local degrees of freedom once radiative corrections are taken into account. Coupling to matter is also considered: we show that coupling these theories to matter in a consistent way is not as straightforward as one could expect. Minimal coupling, as well as the most straightforward non-minimal couplings, cannot be used. Therefore, before being able to address any issues in the presence of matter, it would be necessary to find a consistent (and in any case rather peculiar) coupling scheme.


title: Yes, but Did It Work?: Evaluating Variational Inference

abstract: While it's always possible to compute a variational approximation to a posterior distribution, it can be difficult to discover problems with this approximation. We propose two diagnostic algorithms to alleviate this problem. The Pareto-smoothed importance sampling (PSIS) diagnostic gives a goodness of fit measurement for joint distributions, while simultaneously improving the error in the estimate. The variational simulation-based calibration (VSBC) assesses the average performance of point estimates.


title: Multiple positive bound states for critical Schr\"odinger-Poisson systems

abstract: Using variational methods we prove some results about existence and multiplicity of positive bound states of to the following Schr\"odinger-Poisson system: $$ \left\{ \begin{array}{l} \vspace{2mm} -\Delta u+V(x)u+K(x)\phi(x)u=u^5 -\Delta \phi =K(x)u^2\qquad x\in\R^3 \end{array}\right.\quad\quad (SP) $$ We remark that $(SP)$ exhibits a "double" lack of compactness because of the unboundedness of $\R^3$ and the critical growth of the nonlinear term and that in our assumptions ground state solutions of $(SP)$ do not exist.


title: Molecular Regulation of Histamine Synthesis

abstract: Histamine is a critical mediator of IgE/ cell-mediated anaphylaxis, a neurotransmitter and a regulator of gastric acid secretion. Histamine is a monoamine synthesized from the amino acid histidine through a reaction catalyzed by the enzyme histidine decarboxylase (HDC), which removes carboxyl group from histidine. Despite the importance of histamine, transcriptional regulation of HDC gene expression in mammals is still poorly understood. In this Review, we focus on discussing advances in the understanding of molecular regulation of mammalian histamine synthesis.


title: A Morphing Continuum Analysis of Energy Transfer in Compressible Turbulence

abstract: A shock-preserving finite volume solver with the generalized Lax-Friedrichs splitting flux for Morphing Continuum Theory (MCT) is presented and verified. The numerical MCT solver is showcased in a supersonic turbulent flow with Mach 2.93 over an $8^{\circ}$ compression ramp. The simulation results validated MCT with experiments as an alternative for modeling compressible turbulence. The required size of the smallest mesh cell for the MCT simulation is shown to be almost an order larger than that in a similar DNS study. The comparison shows MCT is a much more computationally friendly theory than the classical NS equations. The dynamics of energy cascade at the length-scale of individual eddies is illuminated through the subscale rotation introduced by MCT. In this regard, MCT provides a statistical averaging procedure for capturing energy transfer in compressible turbulence, not found in classical fluid theories. Analysis of the MCT results show the existence of a statistical coupling of the internal and translational kinetic energy fluctuations with the corresponding eddy rotational energy fluctuations, indicating a multiscale transfer of energy. In conclusion, MCT gives a new characterization of the energy cascade within compressible turbulence without the use of excessive computational resources.


title: Disordered fermions, extra dimensions and a solvable Yang-Mills theory

abstract: Generalizing disorder couplings of the SYK model by means of SU(N) matrices we formulate a lattice model of fermions in d+1 dimensions. Integration of fermions yields an effective theory of Yang-Mills fields in d dimensions, the latter approaching the standard Yang-Mills theory in the case of heavy fermions and the classical limit of vanishing coupling constant of the theory. Quantum mechanically, the theory is solved using large N approximation of the dual effective theory of Hermitian matrices in d dimensions. The theory is asymptotically free and confines the color. In case of massless fermions the emerging theory is an asymptotic safe QCD theory. We discuss also the relationship of this theory to the SYK model.


title: Self-stabilizing processes

abstract: We construct `self-stabilizing' processes {Z(t), t $\in [t_0,t_1)$}. These are random processes which when `localized', that is scaled around t to a fine limit, have the distribution of an $\alpha$(Z(t))-stable process, where $\alpha$ is some given function on R. Thus the stability index at t depends on the value of the process at t. Here we address the case where $\alpha$: R $\to$ (0,1). We first construct deterministic functions which satisfy a kind of autoregressive property involving sums over a plane point set $\Pi$. Taking $\Pi$ to be a Poisson point process then defines a random pure jump process, which we show has the desired localized distributions.


title: A Dynamic Programming Approach to Evaluating Multivariate Gaussian Probabilities

abstract: We propose a method of approximating multivariate Gaussian probabilities using dynamic programming. We show that solving the optimization problem associated with a class of discrete-time finite horizon Markov decision processes with non-Lipschitz cost functions is equivalent to integrating a Gaussian functions over polytopes. An approximation scheme for this class of MDPs is proposed and explicit error bounds under the supremum norm for the optimal cost to go functions are derived.


title: In aqua electrochemistry probed by XPEEM: experimental setup, examples, and challenges

abstract: Recent developments in environmental and liquid cells equipped with electron transparent graphene windows have enabled traditional surface science spectromicroscopy tools, such as X-ray photoelectron spectroscopy (XPS), photoemission electron microscopy (PEEM), and scanning electron microscopy (SEM) to be applied to study solid-liquid and liquid-gas interfaces. Here, we focus on the experimental implementation of PEEM to probe electrified graphene-liquid interfaces using electrolyte-filled microchannel arrays as a new sample platform. We demonstrate the important methodological advantage of these multi-sample arrays: they enable the combination of the wide field of view hyperspectral imaging capabilities from PEEM with the use of powerful data mining algorithms to reveal spectroscopic and temporal behaviors at the level of the individual microsample or the entire array ensemble


title: A mathematical model of the metabolism of a cell. Self-organization and chaos

abstract: Using the classical tools of nonlinear dynamics, we study the process of self-organization and the appearance of the chaos in the metabolic process in a cell with the help of a mathematical model of the transformation of steroids by a cell Arthrobacter globiformis. We constructed the phase-parametric diagrams obtained under a variation of the dissipation of the kinetic membrane potential. The oscillatory modes obtained are classified as regular and strange attractors. We calculated the bifurcations, by which the self-organization and the chaos occur in the system, and the transitions "chaos-order", "order-chaos", "order-order", and "chaos-chaos" arise. Feigenbaum's scenarios and the intermittences are found. For some selected modes, the projections of the phase portraits of attractors, Poincar\'e sections, and Poincar\'e maps are constructed. The total spectra of Lyapunov indices for the modes under study are calculated. The structural stability of the attractors is demonstrated. A general scenario of the formation of regular and strange attractors in the given metabolic process in a cell is found. The physical nature of their appearance in the metabolic process is studied.


title: Learning One Convolutional Layer with Overlapping Patches

abstract: We give the first provably efficient algorithm for learning a one hidden layer convolutional network with respect to a general class of (potentially overlapping) patches. Additionally, our algorithm requires only mild conditions on the underlying distribution. We prove that our framework captures commonly used schemes from computer vision, including one-dimensional and two-dimensional "patch and stride" convolutions. Our algorithm-- $Convotron$ -- is inspired by recent work applying isotonic regression to learning neural networks. Convotron uses a simple, iterative update rule that is stochastic in nature and tolerant to noise (requires only that the conditional mean function is a one layer convolutional network, as opposed to the realizable setting). In contrast to gradient descent, Convotron requires no special initialization or learning-rate tuning to converge to the global optimum. We also point out that learning one hidden convolutional layer with respect to a Gaussian distribution and just $one$ disjoint patch $P$ (the other patches may be arbitrary) is $easy$ in the following sense: Convotron can efficiently recover the hidden weight vector by updating $only$ in the direction of $P$.


title: Predicting Hurricane Trajectories using a Recurrent Neural Network

abstract: Hurricanes are cyclones circulating about a defined center whose closed wind speeds exceed 75 mph originating over tropical and subtropical waters. At landfall, hurricanes can result in severe disasters. The accuracy of predicting their trajectory paths is critical to reduce economic loss and save human lives. Given the complexity and nonlinearity of weather data, a recurrent neural network (RNN) could be beneficial in modeling hurricane behavior. We propose the application of a fully connected RNN to predict the trajectory of hurricanes. We employed the RNN over a fine grid to reduce typical truncation errors. We utilized their latitude, longitude, wind speed, and pressure publicly provided by the National Hurricane Center (NHC) to predict the trajectory of a hurricane at 6-hour intervals. Results show that this proposed technique is competitive to methods currently employed by the NHC and can predict up to approximately 120 hours of hurricane path.


title: Maurer-Cartan moduli and theorems of Riemann-Hilbert type

abstract: We study Maurer-Cartan moduli spaces of dg algebras and associated dg categories and show that, while not quasi-isomorphism invariants, they are invariants of strong homotopy type, a natural notion that has not been studied before. We prove, in several different contexts, Schlessinger-Stasheff type theorems comparing the notions of homotopy and gauge equivalence for Maurer-Cartan elements as well as their categorified versions. As an application, we re-prove and generalize Block-Smith's higher Riemann-Hilbert correspondence, and develop its analogue for simplicial complexes and topological spaces.


title: Semi-Amortized Variational Autoencoders

abstract: Amortized variational inference (AVI) replaces instance-specific local inference with a global inference network. While AVI has enabled efficient training of deep generative models such as variational autoencoders (VAE), recent empirical work suggests that inference networks can produce suboptimal variational parameters. We propose a hybrid approach, to use AVI to initialize the variational parameters and run stochastic variational inference (SVI) to refine them. Crucially, the local SVI procedure is itself differentiable, so the inference network and generative model can be trained end-to-end with gradient-based optimization. This semi-amortized approach enables the use of rich generative models without experiencing the posterior-collapse phenomenon common in training VAEs for problems like text generation. Experiments show this approach outperforms strong autoregressive and variational baselines on standard text and image datasets.


title: A general existence result for stationary solutions to the Keller-Segel system

abstract: We consider a Liouville-type PDE on a smooth bounded planar domain, which is related to stationary solutions of the Keller-Segel's model for chemotaxis. We prove existence of solutions under some algebraic conditions on the parameters. In particular, if the domain is not simply connected, then we can find solution for a generic choice of the parameters. We use variational and Morse-theoretical methods.


title: Cracking the Conundrum of F-Supergiant Coronae

abstract: Chandra X-ray and HST far-ultraviolet (FUV) observations of three early-F supergiants have shed new light on a previous puzzle involving a prominent member of the class: Alpha Persei (HD20902: F5Ib). The warm supergiant is a moderately strong, hard coronal (T~10MK) X-ray source, but has ten times weaker "sub-coronal" Si IV 139 nm (T~80,000 K) emissions than early-G supergiants of similar high-energy properties. The Alpha Per X-ray excess speculatively was ascribed to a close-in hyperactive G-dwarf companion, which could have escaped previous notice, lost in the glare of the bright star. However, a subsequent dedicated multi-wavelength imaging campaign failed to find any evidence for a resolved secondary. The origin of the Alpha Per high-energy dichotomy then devolved to: (1) an unresolved companion; or (2) intrinsic coronal behavior. Exploring the second possibility, the present program has found that early-F supergiants do appear to belong to a distinct coronal class, characterized by elevated X-ray/FUV ratios, although sharing some similarities with Cepheid variables in their transitory X-ray "high states." Remarkably, the early-F supergiants now are seen to align with the low-activity end of the X-ray/FUV sequence defined by late-type dwarfs, suggesting that the disjoint behavior relative to the G supergiants might be attributed to thinner outer atmospheres on the F types, as in dwarfs, but in this case perhaps caused by a weakened "ionization valve" effect due to overly warm photospheres.


title: Wrinkles, folds and plasticity in granular rafts

abstract: We investigate the mechanical response of a compressed monolayer of large and dense particles at a liquid-fluid interface: a granular raft. Upon compression, rafts first wrinkle; then, as the confinement increases, the deformation localizes in a unique fold. This characteristic buckling pattern is usually associated to floating elastic sheets and as a result, particle laden interfaces are often modeled as such. Here, we push this analogy to its limits by comparing the first quantitative measurements of the raft morphology to a theoretical continuous elastic model of the interface. We show that although powerful to describe the wrinkle wavelength, the wrinkle-to-fold transition and the fold shape, this elastic description does not capture the finer details of the experiment. We describe an unpredicted secondary wavelength, a compression discrepancy with the model and a hysteretic behavior during compression cycles, all of which are a signature of the intrinsic discrete and frictional nature of granular rafts. It suggests also that these composite materials exhibit both plastic transition and jamming dynamics.


title: Local Convergence Properties of SAGA/Prox-SVRG and Acceleration

abstract: Over the past ten years, driven by large scale optimisation problems arising from machine learning, the development of stochastic optimisation methods have witnessed a tremendous growth. However, despite their popularity, the theoretical understandings of these methods are quite limited in contrast to the deterministic optimisation methods. In this paper, we present a local convergence analysis for a typical type of stochastic optimisation methods: proximal variance reduced stochastic gradient methods, and mainly focus on the SAGA [12] and Prox-SVRG [43] algorithms. Under the assumption that the non-smooth component of the optimisation problem is partly smooth relative to a smooth manifold, we present a unified framework for the local convergence analysis of the SAGA/Prox-SVRG algorithms: (i) the sequences generated by the SAGA/Prox-SVRG are able to identify the smooth manifold in a finite number of iterations; (ii) then the sequence enters a local linear convergence regime. Beyond local convergence analysis, we also discuss various possibilities for accelerating these algorithms, including adapting to better local parameters, and applying higher-order deterministic/stochastic optimisation methods which can achieve super-linear convergence. Concrete examples arising from machine learning are considered to verify the obtained results.


title: Correlations of occupation numbers in the canonical ensemble and application to BEC in a 1D harmonic trap

abstract: We study statistical properties of $N$ non-interacting identical bosons or fermions in the canonical ensemble. We derive several general representations for the $p$-point correlation function of occupation numbers $\overline{n_1\cdots n_p}$. We demonstrate that it can be expressed as a ratio of two $p\times p$ determinants involving the (canonical) mean occupations $\overline{n_1}$, ..., $\overline{n_p}$, which can themselves be conveniently expressed in terms of the $k$-body partition functions (with $k\leq N$). We draw some connection with the theory of symmetric functions, and obtain an expression of the correlation function in terms of Schur functions. Our findings are illustrated by revisiting the problem of Bose-Einstein condensation in a 1D harmonic trap, for which we get analytical results. We get the moments of the occupation numbers and the correlation between ground state and excited state occupancies. In the temperature regime dominated by quantum correlations, the distribution of the ground state occupancy is shown to be a truncated Gumbel law. The Gumbel law, describing extreme value statistics, is obtained when the temperature is much smaller than the Bose-Einstein temperature.


title: Current Flow Group Closeness Centrality for Complex Networks

abstract: Current flow closeness centrality (CFCC) has a better discriminating ability than the ordinary closeness centrality based on shortest paths. In this paper, we extend this notion to a group of vertices in a weighted graph, and then study the problem of finding a subset $S$ of $k$ vertices to maximize its CFCC $C(S)$, both theoretically and experimentally. We show that the problem is NP-hard, but propose two greedy algorithms for minimizing the reciprocal of $C(S)$ with provable guarantees using the monotoncity and supermodularity. The first is a deterministic algorithm with an approximation factor $(1-\frac{k}{k-1}\cdot\frac{1}{e})$ and cubic running time; while the second is a randomized algorithm with a $(1-\frac{k}{k-1}\cdot\frac{1}{e}-\epsilon)$-approximation and nearly-linear running time for any $\epsilon > 0$. Extensive experiments on model and real networks demonstrate that our algorithms are effective and efficient, with the second algorithm being scalable to massive networks with more than a million vertices.


title: Neyman-Pearson classification: parametrics and sample size requirement

abstract: The Neyman-Pearson (NP) paradigm in binary classification seeks classifiers that achieve a minimal type II error while enforcing the prioritized type I error controlled under some user-specified level $\alpha$. This paradigm serves naturally in applications such as severe disease diagnosis and spam detection, where people have clear priorities among the two error types. Recently, Tong, Feng and Li (2018) proposed a nonparametric umbrella algorithm that adapts all scoring-type classification methods (e.g., logistic regression, support vector machines, random forest) to respect the given type I error upper bound $\alpha$ with high probability, without specific distributional assumptions on the features and the responses. Universal the umbrella algorithm is, it demands an explicit minimum sample size requirement on class $0$, which is often the more scarce class, such as in rare disease diagnosis applications. In this work, we employ the parametric linear discriminant analysis (LDA) model and propose a new parametric thresholding algorithm, which does not need the minimum sample size requirements on class $0$ observations and thus is suitable for small sample applications such as rare disease diagnosis. Leveraging both the existing nonparametric and the newly proposed parametric thresholding rules, we propose four LDA-based NP classifiers, for both low- and high-dimensional settings. On the theoretical front, we prove NP oracle inequalities for one proposed classifier, where the rate for excess type II error benefits from the explicit parametric model assumption. Furthermore, as NP classifiers involve a sample splitting step of class $0$ observations, we construct a new adaptive sample splitting scheme that can be applied universally to NP classifiers, and this adaptive strategy reduces the type II error of these classifiers.


title: Intentional Control of Type I Error over Unconscious Data Distortion: a Neyman-Pearson Approach to Text Classification

abstract: This paper addresses the challenges in classifying textual data obtained from open online platforms, which are vulnerable to distortion. Most existing classification methods minimize the overall classification error and may yield an undesirably large type I error (relevant textual messages are classified as irrelevant), particularly when available data exhibit an asymmetry between relevant and irrelevant information. Data distortion exacerbates this situation and often leads to fallacious prediction. To deal with inestimable data distortion, we propose the use of the Neyman-Pearson (NP) classification paradigm, which minimizes type II error under a user-specified type I error constraint. Theoretically, we show that the NP oracle is unaffected by data distortion when the class conditional distributions remain the same. Empirically, we study a case of classifying posts about worker strikes obtained from a leading Chinese microblogging platform, which are frequently prone to extensive, unpredictable and inestimable censorship. We demonstrate that, even though the training and test data are susceptible to different distortion and therefore potentially follow different distributions, our proposed NP methods control the type I error on test data at the targeted level. The methods and implementation pipeline proposed in our case study are applicable to many other problems involving data distortion.


title: Mechanical Stresses Estimation in Silicon and Glass Bonded at Elevated Temperature

abstract: During electrostatic bonding, also known as anodic bonding, silicon is bonded to glass by applying an external voltage and simultaneous heating to temperatures of 200...450 $\deg$C. While cooling to working temperature after bonding happened pieces are mutually deformed. Due to linear thermal expansion coefficients mismatch of anodically bonded glass and silicon samples an internal stress state is generated. Such stresses are called thermal mismatch stresses. The aim of this paper is a determination of technological and design solutions to achieve minimal thermal mismatch stresses in resulting bond. The nonlinear dependence of linear thermal expansion coefficients of bonded samples' materials on temperature makes it difficult to minimize thermal mismatch stresses by chosing materials with close average thermal expansion coefficients in particular temperature range. To assess means of lowering thermal mismatch stress in this paper two different ways to describe assembly are used: two thin bonded layers and multilayered composite material. Based on properties of two brands of glass (LK5, Borofloat 33) and silicon used with described mathematical models thermal mismatch stresses at temperature $T_w$ in samples bonded at several different temperatures $T_b$ are evaluated. Bonded silicon surface stress dependence of glass to silicon wafer thickness ratio is evaluated. Based on such evaluations one can say that by varying thickness of glass bonded to silicon one can obtain zero thermal mismatch stress at a particular depth of material or obtain stress of some defined value at this depth. Models of assembly description used in this paper can be used to optimize anodic bonding process parameters. Such usage aimed to minimize thermal mismatch stresses at device working temperatures is presented in this paper.


title: Cubic Preferences and the Character Admissibility Problem

abstract: In multiple-question referendum elections, the separability problem occurs when a voter's preferences on some questions or proposals depend on the predicted outcomes of others. The notion of separability formalizes the study of interdependence in multidimensional preferences, and the character admissibility problem deals with the construction of voter preferences with given separability structures. In this paper, we develop a graph theoretic approach to the character admissibilty problem, using Hamiltonian paths to generate voter preferences. We apply this method specifically to the hypercube graph, defining the class of cubic preferences. We then explore how the algebraic structure of the group of symmetries of the hypercube impacts the separability structures exhibited by cubic preferences. We prove that the characters of cubic preferences satisfy set theoretic properties distinct from those produced by previous methods, and we define two functions to construct cubic preferences. Our results have potential applications to experimental work involving election simulation.


title: Polisis: Automated Analysis and Presentation of Privacy Policies Using Deep Learning

abstract: Privacy policies are the primary channel through which companies inform users about their data collection and sharing practices. These policies are often long and difficult to comprehend. Short notices based on information extracted from privacy policies have been shown to be useful but face a significant scalability hurdle, given the number of policies and their evolution over time. Companies, users, researchers, and regulators still lack usable and scalable tools to cope with the breadth and depth of privacy policies. To address these hurdles, we propose an automated framework for privacy policy analysis (Polisis). It enables scalable, dynamic, and multi-dimensional queries on natural language privacy policies. At the core of Polisis is a privacy-centric language model, built with 130K privacy policies, and a novel hierarchy of neural-network classifiers that accounts for both high-level aspects and fine-grained details of privacy practices. We demonstrate Polisis' modularity and utility with two applications supporting structured and free-form querying. The structured querying application is the automated assignment of privacy icons from privacy policies. With Polisis, we can achieve an accuracy of 88.4% on this task. The second application, PriBot, is the first freeform question-answering system for privacy policies. We show that PriBot can produce a correct answer among its top-3 results for 82% of the test questions. Using an MTurk user study with 700 participants, we show that at least one of PriBot's top-3 answers is relevant to users for 89% of the test questions.


title: Fair-by-design matching

abstract: Matching algorithms are used routinely to match donors to recipients for solid organs transplantation, for the assignment of medical residents to hospitals, record linkage in databases, scheduling jobs on machines, network switching, online advertising, and image recognition, among others. Although many optimal solutions may exist to a given matching problem, when the elements that shall or not be included in a solution correspond to individuals, it becomes of paramount importance that the solution be selected fairly. In this paper we study individual fairness in matching problems. Given that many maximum matchings may exist, each one satisfying a different set of individuals, the only way to guarantee fairness is through randomization. Hence we introduce the distributional maxmin fairness framework which provides, for any given input instance, the strongest guarantee possible simultaneously for all individuals in terms of satisfaction probability (the probability of being matched in the solution). Specifically, a probability distribution over feasible solutions is maxmin-fair if it is not possible to improve the satisfaction probability of any individual without decreasing it for some other individual which is no better off. In the special case of matchings in bipartite graphs, our framework is equivalent to the egalitarian mechanism of Bogomolnaia and Mouline. Our main contribution is a polynomial-time algorithm for fair matching building on techniques from minimum cuts, and edge-coloring algorithms for regular bipartite graphs, and transversal theory. For bipartite graphs, our algorithm runs in $O((|V|^2 + |E||V|^{2/3}) \cdot (\log |V|)^2)$ expected time and scales to graphs with tens of millions of vertices and hundreds of millions of edges. To the best of our knowledge, this provides the first large-scale implementation of the egalitarian mechanism.


title: A global linear and local superlinear/quadratic inexact non-interior continuation method for variational inequalities

abstract: We use the concept of barrier-based smoothing approximations introduced in [ C. B. Chua and Z. Li, A barrier-based smoothing proximal point algorithm for NCPs over closed convex cones, SIOPT 23(2), 2010] to extend the non-interior continuation method proposed in [B. Chen and N. Xiu, A global linear and local quadratic noninterior continuation method for nonlinear complementarity problems based on Chen-Mangasarian smoothing functions, SIOPT 9(3), 1999] to an inexact non-interior continuation method for variational inequalities over general closed convex sets. Newton equations involved in the method are solved inexactly to deal with high dimension problems. The method is proved to have global linear and local superlinear/quadratic convergence under suitable assumptions. We apply the method to non-negative orthants, positive semidefinite cones, polyhedral sets, epigraphs of matrix operator norm cone and epigraphs of matrix nuclear norm cone.


title: Numerical Semigroups Generated by Concatenation of Arithmetic Sequences

abstract: We introduce the notion of numerical semigroups generated by concatenation of arithmetic sequences and show that this class of numerical semigroups exhibit multiple interesting behaviours.


title: Applying Cooperative Machine Learning to Speed Up the Annotation of Social Signals in Large Multi-modal Corpora

abstract: Scientific disciplines, such as Behavioural Psychology, Anthropology and recently Social Signal Processing are concerned with the systematic exploration of human behaviour. A typical work-flow includes the manual annotation (also called coding) of social signals in multi-modal corpora of considerable size. For the involved annotators this defines an exhausting and time-consuming task. In the article at hand we present a novel method and also provide the tools to speed up the coding procedure. To this end, we suggest and evaluate the use of Cooperative Machine Learning (CML) techniques to reduce manual labelling efforts by combining the power of computational capabilities and human intelligence. The proposed CML strategy starts with a small number of labelled instances and concentrates on predicting local parts first. Afterwards, a session-independent classification model is created to finish the remaining parts of the database. Confidence values are computed to guide the manual inspection and correction of the predictions. To bring the proposed approach into application we introduce NOVA - an open-source tool for collaborative and machine-aided annotations. In particular, it gives labellers immediate access to CML strategies and directly provides visual feedback on the results. Our experiments show that the proposed method has the potential to significantly reduce human labelling efforts.


title: Nash multiplicity sequences and Hironaka's order function

abstract: When $X$ is a $d$-dimensional variety defined over a field $k$ of characteristic zero, a constructive resolution of singularities can be achieved by successively lowering the maximum multiplicity via blow ups at smooth equimultiple centers. This is done by stratifying the maximum multiplicity locus of $X$ by means of the so called {\em resolution functions}. The most important of these functions is what we know as {\em Hironaka's order function in dimension $d$}. Actually, this function can be defined for varieties when the base field is perfect; however if the characteristic of $k$ is positive, the function is, in general, too coarse and does not provide enough information so as to define a resolution. It is very natural to ask what the meaning of this function is in this case, and to try to find refinements that could lead, ultimately, to a resolution. In this paper we show that Hironaka's order function in dimension $d$ can be read in terms of the {\em Nash multiplicity sequences} introduced by Lejeune-Jalabert. Therefore, the function is intrinsic to the variety and has a geometrical meaning in terms of its space of arcs.


title: An improved multi-parametric programming algorithm for flux balance analysis of metabolic networks

abstract: Flux balance analysis has proven an effective tool for analyzing metabolic networks. In flux balance analysis, reaction rates and optimal pathways are ascertained by solving a linear program, in which the growth rate is maximized subject to mass-balance constraints. A variety of cell functions in response to environmental stimuli can be quantified using flux balance analysis by parameterizing the linear program with respect to extracellular conditions. However, for most large, genome-scale metabolic networks of practical interest, the resulting parametric problem has multiple and highly degenerate optimal solutions, which are computationally challenging to handle. An improved multi-parametric programming algorithm based on active-set methods is introduced in this paper to overcome these computational difficulties. Degeneracy and multiplicity are handled, respectively, by introducing generalized inverses and auxiliary objective functions into the formulation of the optimality conditions. These improvements are especially effective for metabolic networks because their stoichiometry matrices are generally sparse; thus, fast and efficient algorithms from sparse linear algebra can be leveraged to compute generalized inverses and null-space bases. We illustrate the application of our algorithm to flux balance analysis of metabolic networks by studying a reduced metabolic model of Corynebacterium glutamicum and a genome-scale model of Escherichia coli. We then demonstrate how the critical regions resulting from these studies can be associated with optimal metabolic modes and discuss the physical relevance of optimal pathways arising from various auxiliary objective functions. Achieving more than five-fold improvement in computational speed over existing multi-parametric programming tools, the proposed algorithm proves promising in handling genome-scale metabolic models.


title: VISER: Visual Self-Regularization

abstract: In this work, we propose the use of large set of unlabeled images as a source of regularization data for learning robust visual representation. Given a visual model trained by a labeled dataset in a supervised fashion, we augment our training samples by incorporating large number of unlabeled data and train a semi-supervised model. We demonstrate that our proposed learning approach leverages an abundance of unlabeled images and boosts the visual recognition performance which alleviates the need to rely on large labeled datasets for learning robust representation. To increment the number of image instances needed to learn robust visual models in our approach, each labeled image propagates its label to its nearest unlabeled image instances. These retrieved unlabeled images serve as local perturbations of each labeled image to perform Visual Self-Regularization (VISER). To retrieve such visual self regularizers, we compute the cosine similarity in a semantic space defined by the penultimate layer in a fully convolutional neural network. We use the publicly available Yahoo Flickr Creative Commons 100M dataset as the source of our unlabeled image set and propose a distributed approximate nearest neighbor algorithm to make retrieval practical at that scale. Using the labeled instances and their regularizer samples we show that we significantly improve object categorization and localization performance on the MS COCO and Visual Genome datasets where objects appear in context.


title: First-order structural transition and pressure-induced lattice/phonon anomalies in Sr$_2$IrO$_4$

abstract: We investigate the crystal structure and lattice vibrations of Sr$_2$IrO$_4$ by a combined phonon Raman scattering and x-ray powder diffraction experiment under pressures up to $66$ GPa and room temperature. Density functional theory (DFT) and $ab$-initio lattice dynamics calculations were also carried out. A first-order structural phase transition associated with an $8$ % collapse of the $c$-axis is observed at high pressures, with phase coexistence being observed between $\sim 40$ and $55$ GPa. At lower pressures, lattice and phonon anomalies were observed, reflecting crossovers between isostructural competing states. A critical pressure of $P_1=17$ GPa is associated with: (i) a reduction of lattice volume compressibility and a change of behavior of the tetragonal $c/a$ ratio take place above $P_1$; (ii) a four-fold symmetry-breaking lattice strain associated with lattice disorder; (iii) disappearance of two Raman active modes (at $\sim 180$ and $\sim 260$ cm$^{-1}$); and (iv) development of an asymmetric Fano lineshape for the $\sim 390$ cm$^{-1}$ mode. DFT indicates that the phase above $P_1$ is most likely non-magnetic. Exploring the similarities between iridate and cuprate physics, we argue that these observations are consistent with the emergence of a rotational symmetry-breaking electronic instability at $P_1$, providing hints for the avoided metallization under pressure and supporting the hypothesis of possible competing orders that are detrimental to superconductivity in this family. Alternative scenarios for the transition at $P_1$ are also suggested and critically discussed. Additional phonon and lattice anomalies in the tetragonal phase are observed at $P_2=30$ and $P_3=40$ GPa, indicating further competing phases that are stabilized at high pressures.


title: Fabrication and characterization of pH responsive nanoprobes based on ion current rectification

abstract: In this study, we investigated the ionic current rectification of glass nanopipettes modified with bovine serum albumin - glutaraldehyde (BSA-GA) artificial membrane using solutions with various pHs. Ionic current rectification is a phenomenon that is observed with nanopores as asymmetric I-V curves, where the ionic currents recorded through a nanopore differ at the same magnitude of applied electrical potentials biased with opposite polarities. The results clearly showed that modifying the tip of a nanopipette results in a pH dependent ionic current behavior. The proposed strategy is a facile method for fabrication of a pH responsive nanoprobe that has a potential for intracellular pH measurement.


title: Bitewing Radiography Semantic Segmentation Base on Conditional Generative Adversarial Nets

abstract: Currently, Segmentation of bitewing radiograpy images is a very challenging task. The focus of the study is to segment it into caries, enamel, dentin, pulp, crowns, restoration and root canal treatments. The main method of semantic segmentation of bitewing radiograpy images at this stage is the U-shaped deep convolution neural network, but its accuracy is low. in order to improve the accuracy of semantic segmentation of bitewing radiograpy images, this paper proposes the use of Conditional Generative Adversarial network (cGAN) combined with U-shaped network structure (U-Net) approach to semantic segmentation of bitewing radiograpy images. The experimental results show that the accuracy of cGAN combined with U-Net is 69.7%, which is 13.3% higher than the accuracy of u-shaped deep convolution neural network of 56.4%.


title: Dynamical analysis on $f(R,\mathcal{G})$ cosmology

abstract: We use a dynamical system approach to study the cosmological viability of $f(R,\mathcal{G})$ gravity theories. The method consists of formulating the evolution equations as an autonomous system of ODEs, using suitable variables. The formalism is applied to a class of models in which $f(R,\mathcal{G})\propto R^{n}\mathcal{G}^{1-n}$ and its solutions and corresponding stability are analysed in detail. New accelerating solutions that can be attractors in the phase space are found. We also find that this class of models does not exhibit a matter-dominated epoch, a solution which is inconsistent with current cosmological observations.


title: Zorua: Enhancing Programming Ease, Portability, and Performance in GPUs by Decoupling Programming Models from Resource Management

abstract: The application resource specification--a static specification of several parameters such as the number of threads and the scratchpad memory usage per thread block--forms a critical component of the existing GPU programming models. This specification determines the performance of the application during execution because the corresponding on-chip hardware resources are allocated and managed purely based on this specification. This tight coupling between the software-provided resource specification and resource management in hardware leads to significant challenges in programming ease, portability, and performance, as we demonstrate in this work. Our goal in this work is to reduce the dependence of performance on the software-provided resource specification to simultaneously alleviate the above challenges. To this end, we introduce Zorua, a new resource virtualization framework, that decouples the programmer-specified resource usage of a GPU application from the actual allocation in the on-chip hardware resources. Zorua enables this decoupling by virtualizing each resource transparently to the programmer. We demonstrate that by providing the illusion of more resources than physically available, Zorua offers several important benefits: (i) Programming Ease: Zorua eases the burden on the programmer to provide code that is tuned to efficiently utilize the physically available on-chip resources. (ii) Portability: Zorua alleviates the necessity of re-tuning an application's resource usage when porting the application across GPU generations. (iii) Performance: By dynamically allocating resources and carefully oversubscribing them when necessary, Zorua improves or retains the performance of applications that are already highly tuned to best utilize the resources. The holistic virtualization provided by Zorua has many other potential uses which we describe in this paper.


title: Transport of the light polarization in the weak gravitational wave background

abstract: The influence of the weak gravitational wave on the light polarization is considered. Oscillations in the direction of the polarization vector is found.


title: New variable Stars from the Photographic Archive: Semi-automated Discoveries, Attempts of Automatic Classification, and the New Field 104 Her

abstract: Using 172 plates taken with the 40-cm astrograph of the Sternberg Astronomical Institute (Lomonosov Moscow University) in 1976-1994 and digitized with the resolution of 2400 dpi, we discovered and studied 275 new variable stars. We present the list of our new variables with all necessary information concerning their brightness variations. As in our earlier studies, the new discoveries show a rather large number of high-amplitude Delta Scuti variables, predicting that many stars of this type remain not detected in the whole sky. We also performed automated classification of the newly discovered variable stars based on the Random Forest algorithm. The results of the automated classification were compared to traditional classification and showed that automated classification was possible even with noisy photographic data. However, further improvement of automated techniques is needed, which is especially important having in mind the very large numbers of new discoveries expected from all-sky surveys.


title: The Chemical Homogeneity of Sun-like Stars in the Solar Neighborhood

abstract: The compositions of stars are a critical diagnostic tool for many topics in astronomy such as the evolution of our Galaxy, the formation of planets, and the uniqueness of the Sun. Previous spectroscopic measurements indicate a large intrinsic variation in the elemental abundance patterns of stars with similar overall metal content. However, systematic errors arising from inaccuracies in stellar models are known to be a limiting factor in such studies, and thus it is uncertain to what extent the observed diversity of stellar abundance patterns is real. Here we report the abundances of 30 elements with precisions of 2% for 79 Sun-like stars within 100 parsecs. Systematic errors are minimized in this study by focusing on solar twin stars and performing a line-by-line differential analysis using high-resolution, high-signal-to-noise spectra. We resolve [X/Fe] abundance trends in galactic chemical evolution at precisions of $10^{-3}$ dex Gyr$^{-1}$ and reveal that stars with similar ages and metallicities have nearly identical abundance patterns. Contrary to previous results, we find that the ratios of carbon-to-oxygen and magnesium-to-silicon in solar metallicity stars are homogeneous to within 10% throughout the solar neighborhood, implying that exoplanets may exhibit much less compositional diversity than previously thought. Finally, we demonstrate that the Sun has a subtle deficiency in refractory material relative to >80% of solar twins (at 2$\sigma$ confidence), suggesting a possible signpost for planetary systems like our own.


title: Supernova Neutrino Neutrino Astronomy

abstract: Modern neutrino facilities will be able to detect a large number of neutrinos from the next Galactic supernova. We investigate the viability of the triangulation method to locate a core-collapse supernova by employing the neutrino arrival time differences at various detectors. We perform detailed numerical fits in order to determine the uncertainties of these time differences for the cases when the core collapses into a neutron star or a black hole. We provide a global picture by combining all the relevant current and future neutrino detectors. Our findings indicate that in the scenario of a neutron star formation, supernova can be located with precision of 1.5 and 3.5 degrees in declination and right ascension, respectively. For the black hole scenario, sub-degree precision can be reached.


title: Signatures of Cosmic Reionization on the 21cm 2- and 3-point Correlation Function I: Quadratic Bias Modeling

abstract: The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local quadratic bias model. We find that fits of bias model predictions for the 2PCF and 3PCF deviate from our measurements by $\sim 20\%$ at scales above the typical size of ionized regions ($\simeq 30$ Mpc) and at early times with global neutral fractions of $\langle x_{\rm HI} \rangle \gtrsim 0.7$. At later times and smaller scales these deviations increase strongly, indicating a break down of the bias model. The 2PCF and 3PCF fits of the linear bias parameter agree at the $10\%$ level for different EoR model configurations. This agreement holds, when adding redshift space distortions to the simulations. The relation between spatial fluctuations in the matter density and the 21cm signal, as predicted by the bias model, is consistent with direct measurements of this relation in simulations for large smoothing scales ($\gtrsim 30$ Mpc). From this latter test we conclude that negative amplitudes of the 21cm 3PCF result from negative bias parameters, which describe the anti-correlation between the matter over-densities and the 21cm signal during the EoR. However, a more detailed interpretation of the bias parameters may require a description of non-local contributions to the bias model.


title: The Planes of Satellite Galaxies Problem, Suggested Solutions, and Open Questions

abstract: Satellite galaxies of the Milky Way and of the Andromeda galaxy have been found to preferentially align in significantly flattened planes of satellite galaxies, and available velocity measurements are indicative of a preference of satellites in those structures to co-orbit. There is increasing evidence that such kinematically correlated satellite planes are also present around more distant hosts. Detailed comparisons show that similarly anisotropic phase-space distributions of sub-halos are exceedingly rare in cosmological simulations based on the $\Lambda$CDM paradigm. Analogs to the observed systems have frequencies of $\leq 0.5$ per cent in such simulations. In contrast to other small-scale problems, the satellite planes issue is not strongly affected by baryonic processes because the distribution of sub-halos on scales of hundreds of kpc is dominated by gravitational effects. This makes the satellite planes one of the most serious small-scale problem for $\Lambda$CDM. This review summarizes the observational evidence for planes of satellite galaxies in the Local Group and beyond, and provides an overview of how they compare to cosmological simulations. It also discusses scenarios which aim at explaining the coherence of satellite positions and orbits, and why they all are currently unable to satisfactorily resolve the issue.


title: Interactions resolve state-dependence in a toy-model of AdS black holes

abstract: We show that the holographic description of a class of AdS black holes with scalar hair involves dual field theories with a double well effective potential. Black hole microstates have significant support around both vacua in the dual, which correspond to perturbative degrees of freedom on opposite sides of the horizon. A solvable toy-model version of this dual is given by a quantum mechanical particle in a double well potential. In this we show explicitly that the interactions replace the state-dependence that is needed to describe black hole microstates in a low energy effective model involving the tensor product of two decoupled harmonic oscillators. A naive number operator signals the presence of a firewall but a careful construction of perturbative states and operators extinguishes this.


title: Self-consistent redshift estimation using correlation functions without a spectroscopic reference sample

abstract: We present a new method to estimate redshift distributions and galaxy-dark matter bias parameters using correlation functions in a fully data driven and self-consistent manner. Unlike other machine learning, template, or correlation redshift methods, this approach does not require a reference sample with known redshifts. By measuring the projected cross- and auto- correlations of different galaxy sub-samples, e.g., as chosen by simple cells in color-magnitude space, we are able to estimate the galaxy-dark matter bias model parameters, and the shape of the redshift distributions of each sub-sample. This method fully marginalises over a flexible parameterisation of the redshift distribution and galaxy-dark matter bias parameters of sub-samples of galaxies, and thus provides a general Bayesian framework to incorporate redshift uncertainty into the cosmological analysis in a data-driven, consistent, and reproducible manner. This result is improved by an order of magnitude by including cross-correlations with the CMB and with galaxy-galaxy lensing. We showcase how this method could be applied to real galaxies. By using idealised data vectors, in which all galaxy-dark matter model parameters and redshift distributions are known, this method is demonstrated to recover unbiased estimates on important quantities, such as the offset $\Delta_z$ between the mean of the true and estimated redshift distribution and the 68\% and 95\% and 99.5\% widths of the redshift distribution to an accuracy required by current and future surveys.


title: OGLE-2017-BLG-1434Lb: Eighth q < 1 * 10^-4 Mass-Ratio Microlens Planet Confirms Turnover in Planet Mass-Ratio Function

abstract: We report the discovery of a cold Super-Earth planet (m_p=4.4 +/- 0.5 M_Earth) orbiting a low-mass (M=0.23 +/- 0.03 M_Sun) M dwarf at projected separation a_perp = 1.18 +/- 0.10 AU, i.e., about 1.9 times the snow line. The system is quite nearby for a microlensing planet, D_Lens = 0.86 +/- 0.09 kpc. Indeed, it was the large lens-source relative parallax pi_rel=1.0 mas (combined with the low mass M) that gave rise to the large, and thus well-measured, "microlens parallax" that enabled these precise measurements. OGLE-2017-BLG-1434Lb is the eighth microlensing planet with planet-host mass ratio q < 1 * 10^-4. We apply a new planet-detection sensitivity method, which is a variant of "V/V_max", to seven of these eight planets to derive the mass-ratio function in this regime. We find dN/d(ln q) ~ q^p, with p = 1.05 (+0.78,-0.68), which confirms the "turnover" in the mass function found by Suzuki et al. relative to the power law of opposite sign n = -0.93 +/- 0.13 at higher mass ratios q >~ 2 * 10^-4. We combine our result with that of Suzuki et al. to obtain p = 0.73 (+0.42,-0.34).


title: The dragonfly nearby galaxies survey. Iv. A giant stellar disk in ngc 2841

abstract: Neutral gas is commonly believed to dominate over stars in the outskirts of galaxies, and investigations of the disk-halo interface are generally considered to be in the domain of radio astronomy. This may simply be a consequence of the fact that deep HI observations typically probe to a lower mass surface density than visible wavelength data. This paper presents low surface brightness optimized visible wavelength observations of the extreme outskirts of the nearby spiral galaxy NGC 2841. We report the discovery of an enormous low-surface brightness stellar disk in this object. When azimuthally averaged, the stellar disk can be traced out to a radius of $\sim$70 kpc (5 $R_{25}$ or 23 inner disk scale lengths). The structure in the stellar disk traces the morphology of HI emission and extended UV emission. Contrary to expectations, the stellar mass surface density does not fall below that of the gas mass surface density at any radius. In fact, at all radii greater than $\sim$20 kpc, the ratio of the stellar to gas mass surface density is a constant 3:1. Beyond $\sim$30 kpc, the low surface brightness stellar disk begins to warp, which may be an indication of a physical connection between the outskirts of the galaxy and infall from the circumgalactic medium. A combination of stellar migration, accretion and in-situ star formation might be responsible for building up the outer stellar disk, but whatever mechanisms formed the outer disk must also explain the constant ratio between stellar and gas mass in the outskirts of this galaxy.