UPMASK: unsupervised photometric membership assignment in stellar clusters

We develop a method for membership assignment in stellar clusters using only photometry and positions. The method, UPMASK, is aimed to be unsupervised, data driven, model free, and to rely on as few assumptions as possible. It is based on an iterative process, principal component analysis, clustering algorithm, and kernel density estimations. Moreover, it is able to take into account arbitrary error models. An implementation in R was tested on simulated clusters that covered a broad range of ages, masses, distances, reddenings, and also on real data of cluster fields. Running UPMASK on simulations showed that it effectively separates cluster and field populations. The overall spatial structure and distribution of cluster member stars in the colour-magnitude diagram were recovered under a broad variety of conditions. For a set of 360 simulations, the resulting true positive rates (a measurement of purity) and member recovery rates (a measurement of completeness) at the 90% membership probability level reached high values for a range of open cluster ages ($10^{7.1}-10^{9.5}$ yr), initial masses ($0.5-10\times10^3$M_{\sun}) and heliocentric distances ($0.5-4.0$ kpc). UPMASK was also tested on real data from the fields of the open cluster Haffner~16 and of the closely projected clusters Haffner~10 and Czernik~29. These tests showed that even for moderate variable extinction and cluster superposition, the method yielded useful cluster membership probabilities and provided some insight into their stellar contents. The UPMASK implementation will be available at the CRAN archive.Comment: 12 pages, 13 figures, accepted for publication in Astronomy and Astrophysic

The first analytical expression to estimate photometric redshifts suggested by a machine

We report the first analytical expression purely constructed by a machine to determine photometric redshifts ($z_{\rm phot}$) of galaxies. A simple and reliable functional form is derived using $41,214$ galaxies from the Sloan Digital Sky Survey Data Release 10 (SDSS-DR10) spectroscopic sample. The method automatically dropped the $u$ and $z$ bands, relying only on $g$, $r$ and $i$ for the final solution. Applying this expression to other $1,417,181$ SDSS-DR10 galaxies, with measured spectroscopic redshifts ($z_{\rm spec}$), we achieved a mean $\langle (z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})\rangle\lesssim 0.0086$ and a scatter $\sigma_{(z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})}\lesssim 0.045$ when averaged up to $z \lesssim 1.0$. The method was also applied to the PHAT0 dataset, confirming the competitiveness of our results when faced with other methods from the literature. This is the first use of symbolic regression in cosmology, representing a leap forward in astronomy-data-mining connection.Comment: 6 pages, 4 figures. Accepted for publication in MNRAS Letter

Spatial opinion dynamics and the effects of two types of mixing

Spatially situated opinions that can be held with different degrees of conviction lead to spatiotemporal patterns such as clustering (homophily), polarization, and deadlock. Our goal is to understand how sensitive these patterns are to changes in the local nature of interactions. We introduce two different mixing mechanisms, spatial relocation and nonlocal interaction (“telephoning”), to an earlier fully spatial model (no mixing). Interestingly, the mechanisms that create deadlock in the fully spatial model have the opposite effect when there is a sufficient amount of mixing. With telephoning, not only is polarization and deadlock broken up, but consensus is hastened. The effects of mixing by relocation are even more pronounced. Further insight into these dynamics is obtained for selected parameter regimes via comparison to the mean-field differential equations

His story/her story: A dialogue about including men and masculinities in the women’s studies curriculum

The article discusses the issue of inclusion of men and masculinities in the Women\u27s Studies curriculum. Women\u27s Studies programs were started to compensate for the male domination in the academics. Women\u27s Studies presented a platform where scholarship for women was produced and taken seriously, female students and faculty could find their say or voice, and theoretical investigations required for the advancement of the aims of the women\u27s movement could take place. If the academy as a whole does not sufficiently integrate Women\u27s Studies into the curriculum, integrating Men\u27s Studies into Women\u27s Studies might end up further marginalizing Women\u27s Studies by decreasing the number of classroom hours students spend engaging women\u27s lives and feminist scholarship. Such an integration would presents an another form of male privilege, with men manipulating their way into the only branch of scholarship that has consistently focused on women. On a ground level, feminist scholars are apprehensive that a move from a Women\u27s Studies program to a Gender Studies program will reduce the political aspect of women\u27s programs

NEW EVIDENCES OF THE LOW-P/HIGH-T PRE-ALPINE METAMORPHISM AND MEDIUM-P ALPINE OVERPRINT OF THE PELAGONIAN ZONE DOCUMENTED IN METAPELITES AND ORTHOGNEISSES FROM THE VORAS MASSIF, MACEDONIA, NORTHERN GREECE

Pelitic rocks of the pre-Alpine "gneiss series" from the Voras Massif record a polymetamorphic history of three metamorphic events; a first HT event indicated by migmatitization phenomena, a subsequent LP-HT event (andalusite-sillimanite series) at 2.5-3 Kbar and 610-640°C, and a third medium-P event at ~11 Kbar and ~550°C. The LP-HT metamorphism is associated with intrusion of Variscan granitoid plutois in a magmatic arc setting. The medium-P metamorphism records an Alpine event (probably Early Cretaceous) indicated by kyanite, chlohtoid, garnet and staurolite formation replacing andalusite and cordierite. The Variscan granitoids and the overlying lithologies of the "schist series" are affected only by the medium-P metamorphism. Granitoids are transformed into phengite orthogneisses having the mineral assemblage phen (Si=3.35-3.47 a.p.f.u.)±Bt-Ab-Kfs- Qtz±Czo±Grt and high-alumina pelites into gamet-chloritoid schists with the mineral assemblage Grt-Ctd-Chl-Ms-Qtz-Rt±llm. The "schist series" which underlies the Late Triassic-Early Jurassic marbles and overlies lithologies of the "gneiss series" and metagranitoids represents protoliths of the Permo-Triassic volcanosedimentary series deposited at the eastern margin of the Pelagonian continental bloc

Building an Interdisciplinary Partnership to Improve Math Skills in the Science Classroom

This presentation is intended for University STEM faculty interested in developing interdisciplinary collaborations to improve the ability of students to do math in their science courses. The development of a partnership between faculty in chemistry and mathematics will be discussed. The presenters’ partnership explored ways to improve rates of chemistry student success by identifying key mathematical stumbling blocks encountered by chemistry students. This collaboration resulted in the formation of a student learning community comprised of a General Chemistry course and Precalculus course. Modifications were made to both courses including the addition of material to the mathematics course and changes in pedagogy in the chemistry course. Attendees will be presented with the framework that was used to develop this partnership and will engage in interdisciplinary discussions with other participants. We hope to promote the development of interdisciplinary teams at other colleges and universities

Detecting stars, galaxies, and asteroids with Gaia

(Abridged) Gaia aims to make a 3-dimensional map of 1,000 million stars in our Milky Way to unravel its kinematical, dynamical, and chemical structure and evolution. Gaia's on-board detection software discriminates stars from spurious objects like cosmic rays and Solar protons. For this, parametrised point-spread-function-shape criteria are used. This study aims to provide an optimum set of parameters for these filters. We developed an emulation of the on-board detection software, which has 20 free, so-called rejection parameters which govern the boundaries between stars on the one hand and sharp or extended events on the other hand. We evaluate the detection and rejection performance of the algorithm using catalogues of simulated single stars, double stars, cosmic rays, Solar protons, unresolved galaxies, and asteroids. We optimised the rejection parameters, improving - with respect to the functional baseline - the detection performance of single and double stars, while, at the same time, improving the rejection performance of cosmic rays and of Solar protons. We find that the minimum separation to resolve a close, equal-brightness double star is 0.23 arcsec in the along-scan and 0.70 arcsec in the across-scan direction, independent of the brightness of the primary. We find that, whereas the optimised rejection parameters have no significant impact on the detectability of de Vaucouleurs profiles, they do significantly improve the detection of exponential-disk profiles. We also find that the optimised rejection parameters provide detection gains for asteroids fainter than 20 mag and for fast-moving near-Earth objects fainter than 18 mag, albeit this gain comes at the expense of a modest detection-probability loss for bright, fast-moving near-Earth objects. The major side effect of the optimised parameters is that spurious ghosts in the wings of bright stars essentially pass unfiltered.Comment: Accepted for publication in A&

Genealogies of rapidly adapting populations

The genetic diversity of a species is shaped by its recent evolutionary history and can be used to infer demographic events or selective sweeps. Most inference methods are based on the null hypothesis that natural selection is a weak or infrequent evolutionary force. However, many species, particularly pathogens, are under continuous pressure to adapt in response to changing environments. A statistical framework for inference from diversity data of such populations is currently lacking. Toward this goal, we explore the properties of genealogies in a model of continual adaptation in asexual populations. We show that lineages trace back to a small pool of highly fit ancestors, in which almost simultaneous coalescence of more than two lineages frequently occurs. While such multiple mergers are unlikely under the neutral coalescent, they create a unique genetic footprint in adapting populations. The site frequency spectrum of derived neutral alleles, for example, is non-monotonic and has a peak at high frequencies, whereas Tajima's D becomes more and more negative with increasing sample size. Since multiple merger coalescents emerge in many models of rapid adaptation, we argue that they should be considered as a null-model for adapting populations.Comment: to appear in PNA

Using gamma regression for photometric redshifts of survey galaxies

Machine learning techniques offer a plethora of opportunities in tackling big data within the astronomical community. We present the set of Generalized Linear Models as a fast alternative for determining photometric redshifts of galaxies, a set of tools not commonly applied within astronomy, despite being widely used in other professions. With this technique, we achieve catastrophic outlier rates of the order of ~1%, that can be achieved in a matter of seconds on large datasets of size ~1,000,000. To make these techniques easily accessible to the astronomical community, we developed a set of libraries and tools that are publicly available.Comment: Refereed Proceeding of "The Universe of Digital Sky Surveys" conference held at the INAF - Observatory of Capodimonte, Naples, on 25th-28th November 2014, to be published in the Astrophysics and Space Science Proceedings, edited by Longo, Napolitano, Marconi, Paolillo, Iodice, 6 pages, and 1 figur