Showing 1–50 of 53 results for author: Goldman, S

The Llama 3 Herd of Models

Authors: Abhimanyu Dubey, Abhinav Jauhri, Abhinav Pandey, Abhishek Kadian, Ahmad Al-Dahle, Aiesha Letman, Akhil Mathur, Alan Schelten, Amy Yang, Angela Fan, Anirudh Goyal, Anthony Hartshorn, Aobo Yang, Archi Mitra, Archie Sravankumar, Artem Korenev, Arthur Hinsvark, Arun Rao, Aston Zhang, Aurelien Rodriguez, Austen Gregerson, Ava Spataru, Baptiste Roziere, Bethany Biron, Binh Tang , et al. (508 additional authors not shown)

Abstract: Modern artificial intelligence (AI) systems are powered by foundation models. This paper presents a new set of foundation models, called Llama 3. It is a herd of language models that natively support multilinguality, coding, reasoning, and tool usage. Our largest model is a dense Transformer with 405B parameters and a context window of up to 128K tokens. This paper presents an extensive empirical… ▽ More Modern artificial intelligence (AI) systems are powered by foundation models. This paper presents a new set of foundation models, called Llama 3. It is a herd of language models that natively support multilinguality, coding, reasoning, and tool usage. Our largest model is a dense Transformer with 405B parameters and a context window of up to 128K tokens. This paper presents an extensive empirical evaluation of Llama 3. We find that Llama 3 delivers comparable quality to leading language models such as GPT-4 on a plethora of tasks. We publicly release Llama 3, including pre-trained and post-trained versions of the 405B parameter language model and our Llama Guard 3 model for input and output safety. The paper also presents the results of experiments in which we integrate image, video, and speech capabilities into Llama 3 via a compositional approach. We observe this approach performs competitively with the state-of-the-art on image, video, and speech recognition tasks. The resulting models are not yet being broadly released as they are still under development. △ Less

Submitted 31 July, 2024; originally announced July 2024.

A Lyapunov theory demonstrating a fundamental limit on the speed of systems consolidation

Authors: Alireza Alemi, Emre R. F. Aksay, Mark S. Goldman

Abstract: The nervous system reorganizes memories from an early site to a late site, a commonly observed feature of learning and memory systems known as systems consolidation. Previous work has suggested learning rules by which consolidation may occur. Here, we provide conditions under which such rules are guaranteed to lead to stable convergence of learning and consolidation. We use the theory of Lyapunov… ▽ More The nervous system reorganizes memories from an early site to a late site, a commonly observed feature of learning and memory systems known as systems consolidation. Previous work has suggested learning rules by which consolidation may occur. Here, we provide conditions under which such rules are guaranteed to lead to stable convergence of learning and consolidation. We use the theory of Lyapunov functions, which enforces stability by requiring learning rules to decrease an energy-like (Lyapunov) function. We present the theory in the context of a simple circuit architecture motivated by classic models of learning in systems consolidation mediated by the cerebellum. Stability is only guaranteed if the learning rate in the late stage is not faster than the learning rate in the early stage. Further, the slower the learning rate at the late stage, the larger the perturbation the system can tolerate with a guarantee of stability. We provide intuition for this result by mapping the consolidation model to a damped driven oscillator system, and showing that the ratio of early- to late-stage learning rates in the consolidation model can be directly identified with the (square of the) oscillator's damping ratio. This work suggests the power of the Lyapunov approach to provide constraints on nervous system function. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 16 pages, 4 figures

A MultiWavelength Study of the Symbiotic Mira HM Sge with SOFIA and HST

Authors: Steven Goldman, Ravi Sankrit, Edward Montiel, Sean Garner, Nathan Wolthuis, Nicole Karnath

Abstract: We have targeted the dusty symbiotic mira system HM Sge with four instruments from the IR to the UV. We have used these observations along with archival observations to study how the system has been evolving after its 1975 nova-like outburst. We have detected ro-vibrational water emission in a symbiotic system for the first time using new EXES high spectral resolution infrared spectroscopy. The fe… ▽ More We have targeted the dusty symbiotic mira system HM Sge with four instruments from the IR to the UV. We have used these observations along with archival observations to study how the system has been evolving after its 1975 nova-like outburst. We have detected ro-vibrational water emission in a symbiotic system for the first time using new EXES high spectral resolution infrared spectroscopy. The features, detected in emission, have velocities consistent with the systemic velocity but do not show any clear evidence of high velocity outflows. Mid-infrared photometry and grism spectroscopy show that the oxygen-rich Asymptotic Giant Branch (AGB) dust and dust output has shown little to no change over the past 39 years. In the optical/UV, we detect three main [NII] nebular features that were detected 22 years ago. Two of these features show a small amount of movement corresponding to average outflows speeds of 38 kms and 78 kms since they were previously observed; some previously detected [NII] features are no longer visible. New UV spectroscopy has shown that the nebular environment continues to steadily relax after the system's 1975 outburst. The data suggest however, that the hot component has increased in temperature from 200,000 K in 1989 to now greater than 250,000 K. Our new and archival observations suggest that the evolution of the system after its outburst is swift with little to no major changes after a period of a couple years. △ Less

Submitted 15 January, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: Accepted in ApJ

Journal ref: ApJ 961 14 2024

Vast TVB parameter space exploration: A Modular Framework for Accelerating the Multi-Scale Simulation of Human Brain Dynamics

Authors: Michiel van der Vlag, Lionel Kusch, Alain Destexhe, Viktor Jirsa, Sandra Diaz-Pier, Jennifer S. Goldman

Abstract: Global neural dynamics emerge from multi-scale brain structures, with neurons communicating through synapses to form transiently communicating networks. Network activity arises from intercellular communication that depends on the structure of connectome tracts and local connection, intracellular signalling cascades, and the extracellular molecular milieu that regulate cellular properties. Multi-sc… ▽ More Global neural dynamics emerge from multi-scale brain structures, with neurons communicating through synapses to form transiently communicating networks. Network activity arises from intercellular communication that depends on the structure of connectome tracts and local connection, intracellular signalling cascades, and the extracellular molecular milieu that regulate cellular properties. Multi-scale models of brain function have begun to directly link the emergence of global brain dynamics in conscious and unconscious brain states to microscopic changes at the level of cells. In particular, AdEx mean-field models representing statistical properties of local populations of neurons have been connected following human tractography data to represent multi-scale neural phenomena in simulations using The Virtual Brain (TVB). While mean-field models can be run on personal computers for short simulations, or in parallel on high-performance computing (HPC) architectures for longer simulations and parameter scans, the computational burden remains high and vast areas of the parameter space remain unexplored. In this work, we report that our TVB-HPC framework, a modular set of methods used here to implement the TVB-AdEx model for GPU and analyze emergent dynamics, notably accelerates simulations and substantially reduces computational resource requirements. The framework preserves the stability and robustness of the TVB-AdEx model, thus facilitating finer resolution exploration of vast parameter spaces as well as longer simulations previously near impossible to perform. Given that simulation and analysis toolkits are made public as open-source packages, our framework serves as a template onto which other models can be easily scripted and personalized datasets can be used for studies of inter-individual variability of parameters related to functional brain dynamics. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 21 pages, 9 figures

MIST-CF: Chemical formula inference from tandem mass spectra

Authors: Samuel Goldman, Jiayi Xin, Joules Provenzano, Connor W. Coley

Abstract: Chemical formula annotation for tandem mass spectrometry (MS/MS) data is the first step toward structurally elucidating unknown metabolites. While great strides have been made toward solving this problem, the current state-of-the-art method depends on time-intensive, proprietary, and expert-parameterized fragmentation tree construction and scoring. In this work we extend our previous spectrum Tran… ▽ More Chemical formula annotation for tandem mass spectrometry (MS/MS) data is the first step toward structurally elucidating unknown metabolites. While great strides have been made toward solving this problem, the current state-of-the-art method depends on time-intensive, proprietary, and expert-parameterized fragmentation tree construction and scoring. In this work we extend our previous spectrum Transformer methodology into an energy based modeling framework, MIST-CF, for learning to rank chemical formula and adduct assignments given an unannotated MS/MS spectrum. Importantly, MIST-CF learns in a data dependent fashion using a Formula Transformer neural network architecture and circumvents the need for fragmentation tree construction. We train and evaluate our model on a large open-access database, showing an absolute improvement of 10% top 1 accuracy over other neural network architectures. We further validate our approach on the CASMI2022 challenge dataset, achieving nearly equivalent performance to the winning entry within the positive mode category without any manual curation or post-processing of our results. These results demonstrate an exciting strategy to more powerfully leverage MS2 fragment peaks for predicting MS1 precursor chemical formula with data driven learning. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Single-model uncertainty quantification in neural network potentials does not consistently outperform model ensembles

Authors: Aik Rui Tan, Shingo Urata, Samuel Goldman, Johannes C. B. Dietschreit, Rafael Gómez-Bombarelli

Abstract: Neural networks (NNs) often assign high confidence to their predictions, even for points far out-of-distribution, making uncertainty quantification (UQ) a challenge. When they are employed to model interatomic potentials in materials systems, this problem leads to unphysical structures that disrupt simulations, or to biased statistics and dynamics that do not reflect the true physics. Differentiab… ▽ More Neural networks (NNs) often assign high confidence to their predictions, even for points far out-of-distribution, making uncertainty quantification (UQ) a challenge. When they are employed to model interatomic potentials in materials systems, this problem leads to unphysical structures that disrupt simulations, or to biased statistics and dynamics that do not reflect the true physics. Differentiable UQ techniques can find new informative data and drive active learning loops for robust potentials. However, a variety of UQ techniques, including newly developed ones, exist for atomistic simulations and there are no clear guidelines for which are most effective or suitable for a given case. In this work, we examine multiple UQ schemes for improving the robustness of NN interatomic potentials (NNIPs) through active learning. In particular, we compare incumbent ensemble-based methods against strategies that use single, deterministic NNs: mean-variance estimation, deep evidential regression, and Gaussian mixture models. We explore three datasets ranging from in-domain interpolative learning to more extrapolative out-of-domain generalization challenges: rMD17, ammonia inversion, and bulk silica glass. Performance is measured across multiple metrics relating model error to uncertainty. Our experiments show that none of the methods consistently outperformed each other across the various metrics. Ensembling remained better at generalization and for NNIP robustness; MVE only proved effective for in-domain interpolation, while GMM was better out-of-domain; and evidential regression, despite its promise, was not the preferable alternative in any of the cases. More broadly, cost-effective, single deterministic models cannot yet consistently match or outperform ensembling for uncertainty quantification in NNIPs. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: 27 pages, 4 figures, Supporting Information (22 pages)

Generating Molecular Fragmentation Graphs with Autoregressive Neural Networks

Authors: Samuel Goldman, Janet Li, Connor W. Coley

Abstract: The accurate prediction of tandem mass spectra from molecular structures has the potential to unlock new metabolomic discoveries by augmenting the community's libraries of experimental reference standards. Cheminformatic spectrum prediction strategies use a "bond-breaking" framework to iteratively simulate mass spectrum fragmentations, but these methods are (a) slow, due to the need to exhaustivel… ▽ More The accurate prediction of tandem mass spectra from molecular structures has the potential to unlock new metabolomic discoveries by augmenting the community's libraries of experimental reference standards. Cheminformatic spectrum prediction strategies use a "bond-breaking" framework to iteratively simulate mass spectrum fragmentations, but these methods are (a) slow, due to the need to exhaustively and combinatorially break molecules and (b) inaccurate, as they often rely upon heuristics to predict the intensity of each resulting fragment; neural network alternatives mitigate computational cost but are black-box and not inherently more accurate. We introduce a physically-grounded neural approach that learns to predict each breakage event and score the most relevant subset of molecular fragments quickly and accurately. We evaluate our model by predicting spectra from both public and private standard libraries, demonstrating that our hybrid approach offers state of the art prediction accuracy, improved metabolite identification from a database of candidates, and higher interpretability when compared to previous breakage methods and black box neural networks. The grounding of our approach in physical fragmentation events shows especially high promise for elucidating natural product molecules with more complex scaffolds. △ Less

Submitted 8 January, 2024; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: 25 pages, 49 pages with references and appendix, 5 main text figures

Prefix-Tree Decoding for Predicting Mass Spectra from Molecules

Authors: Samuel Goldman, John Bradshaw, Jiayi Xin, Connor W. Coley

Abstract: Computational predictions of mass spectra from molecules have enabled the discovery of clinically relevant metabolites. However, such predictive tools are still limited as they occupy one of two extremes, either operating (a) by fragmenting molecules combinatorially with overly rigid constraints on potential rearrangements and poor time complexity or (b) by decoding lossy and nonphysical discretiz… ▽ More Computational predictions of mass spectra from molecules have enabled the discovery of clinically relevant metabolites. However, such predictive tools are still limited as they occupy one of two extremes, either operating (a) by fragmenting molecules combinatorially with overly rigid constraints on potential rearrangements and poor time complexity or (b) by decoding lossy and nonphysical discretized spectra vectors. In this work, we use a new intermediate strategy for predicting mass spectra from molecules by treating mass spectra as sets of molecular formulae, which are themselves multisets of atoms. After first encoding an input molecular graph, we decode a set of molecular subformulae, each of which specify a predicted peak in the mass spectrum, the intensities of which are predicted by a second model. Our key insight is to overcome the combinatorial possibilities for molecular subformulae by decoding the formula set using a prefix tree structure, atom-type by atom-type, representing a general method for ordered multiset decoding. We show promising empirical results on mass spectra prediction tasks. △ Less

Submitted 3 December, 2023; v1 submitted 11 March, 2023; originally announced March 2023.

Exact Renormalization of Wave Functionals yields Continuous MERA

Authors: Samuel Goldman, Nima Lashkari, Robert G. Leigh, Mudassir Moosa

Abstract: The exact renormalization group (ERG) is a powerful tool for understanding the formal properties of field theories. By adapting generalized ERG schemes to the flow of wavefunctionals, we obtain a large class of continuous unitary networks, a special case of which includes a class of Gaussian continuous Multi-scale Renormalization Ansatzes (cMERAs). The novel feature of these generalized wavefuncti… ▽ More The exact renormalization group (ERG) is a powerful tool for understanding the formal properties of field theories. By adapting generalized ERG schemes to the flow of wavefunctionals, we obtain a large class of continuous unitary networks, a special case of which includes a class of Gaussian continuous Multi-scale Renormalization Ansatzes (cMERAs). The novel feature of these generalized wavefunctional ERG schemes is allowing for modifications of the dispersion relation, which drastically changes the entanglement structure of the ultraviolet states. Through our construction, we demonstrate that cMERA can be derived from a more fundamental "microscopic" principle, which amounts to the usual RG principle of path integral independence, suitably adapted to quantum states of the field theory. The establishment of such a principle may provide a path forward for exploring cMERA beyond the free field regime, and for understanding the nature of entanglement renormalization intrinsically in the continuum. △ Less

Submitted 18 January, 2024; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 12 pages, 2 figures

Journal ref: Phys. Rev. D 108, 085004 (2023)

A Panchromatic Study of Massive Stars in the Extremely Metal-Poor Local Group Dwarf Galaxy Leo A

Authors: Maude Gull, Daniel R. Weisz, Peter Senchyna, Nathan R. Sandford, Yumi Choi, Anna F. McLeod, Kareem El-Badry, Ylva Götberg, Karoline M. Gilbert, Martha Boyer, Julianne J. Dalcanton, Puragra GuhaThakurta, Steven Goldman, Paola Marigo, Kristen B. W. McQuinn, Giada Pastorelli, Daniel P. Stark, Evan Skillman, Yuan-sen Ting, Benjamin F. Williams

Abstract: We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region… ▽ More We characterize massive stars (M>8 M_sun) in the nearby (D~0.8 Mpc) extremely metal-poor (Z~5% Z_sun) galaxy Leo A using Hubble Space Telescope ultra-violet (UV), optical, and near-infrared (NIR) imaging along with Keck/LRIS and MMT/Binospec optical spectroscopy for 18 main sequence OB stars. We find that: (a) 12 of our 18 stars show emission lines, despite not being associated with an H II region, suggestive of stellar activity (e.g., mass loss, accretion, binary star interaction), which is consistent with previous predictions of enhanced activity at low metallicity; (b) 6 are Be stars, which are the first to be spectroscopically studied at such low metallicity -- these Be stars have unusual panchromatic SEDs; (c) for stars well-fit by the TLUSTY non-local thermodynamic equilibrium (non-LTE) models, the photometric and spectroscopic values of T_eff and log(g) agree to within ~0.01 dex and ~0.18 dex, respectively, indicating that NUV/optical/NIR imaging can be used to reliably characterize massive (M ~ 8-30 M_sun) main sequence star properties relative to optical spectroscopy; (d) the properties of the most massive stars in H II regions are consistent with constraints from previous nebular emission line studies; and (e) 13 stars with M>8 M_sun are >40 pc from a known star cluster or H II region. Our sample comprises ~50% of all known massive stars at Z < 10% Z_sun with derived stellar parameters, high-quality optical spectra, and panchromatic photometry. △ Less

Submitted 28 December, 2022; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: 35 pages, 18 figures

Homologous self-assembled superlattices: What causes their periodic polarity switching? Review, model, and experimental test

Authors: Varun Thakur, Dor Benafsha, Yury Turkulets, Almog R. Azulay, Xin Liang, Rachel S. Goldman, Ilan Shalish

Abstract: Quantum semiconductor structures are commonly achieved by bandgap engineering that relies on the ability to switch from one semiconductor to another during their growth. Growth of a superlattice is typically demanding technologically. In contrast, accumulated evidence points to a tendency among a certain class of multiple-cation binary oxides to self-assemble spontaneously as superlattice structur… ▽ More Quantum semiconductor structures are commonly achieved by bandgap engineering that relies on the ability to switch from one semiconductor to another during their growth. Growth of a superlattice is typically demanding technologically. In contrast, accumulated evidence points to a tendency among a certain class of multiple-cation binary oxides to self-assemble spontaneously as superlattice structures. This class has been dubbed the homologous superlattices. For a famous example, when a mixture of indium and zinc is oxidized, the phases of In-O and ZnO separate in an orderly periodic manner, along the ZnO polar axis, with polarity inversion taking place between consecutive ZnO sections. As we review here, the same structure has been observed when the indium was replaced with other metals, and perhaps even in ZnO alone. This peculiar self-assembled structure has been attracting research over the past decade. The purpose of this study is to gain understanding of the physics underlying the formation of this unique structure. Here, we first provide an extensive review of the accumulated literature on these spontaneously-formed structures and then propose an explanation for the long-standing mystery of this intriguing self-assembly in the form of an electrostatic growth phenomenon and test the proposed model on experimental data. △ Less

Submitted 28 August, 2022; originally announced August 2022.

Journal ref: Advanced Physics Research 2023

A Census of Thermally-Pulsing AGB stars in the Andromeda Galaxy and a First Estimate of their Contribution to the Global Dust Budget

Authors: Steven R. Goldman, Martha L. Boyer, Julianne Dalcanton, Iain McDonald, Leo Girardi, Benjamin F. Williams, Sundar Srinivasan, Karl Gordon

Abstract: We present a near-complete catalog of the metal-rich population of Thermally-Pulsing Asymptotic Giant Branch stars in the northwest quadrant of M31. This metal-rich sample complements the equally complete metal-poor Magellanic Cloud AGB catalogs produced by the SAGE program. Our catalog includes HST wide-band photometry from the Panchromatic Hubble Andromeda Treasury survey, HST medium-band photom… ▽ More We present a near-complete catalog of the metal-rich population of Thermally-Pulsing Asymptotic Giant Branch stars in the northwest quadrant of M31. This metal-rich sample complements the equally complete metal-poor Magellanic Cloud AGB catalogs produced by the SAGE program. Our catalog includes HST wide-band photometry from the Panchromatic Hubble Andromeda Treasury survey, HST medium-band photometry used to chemically classify a subset of the sample, and Spitzer mid- and far-IR photometry that we have used to isolate dust-producing AGB stars. We have detected 346,623 AGB stars; these include 4,802 AGB candidates producing considerable dust, and 1,356 AGB candidates that lie within clusters with measured ages, and in some cases metallicities. Using the Spitzer data and chemical classifications made with the medium-band data, we have identified both carbon- and oxygen-rich AGB candidates producing significant dust. We have applied color--mass-loss relations based on dusty AGB stars from the LMC to estimate the dust injection by AGB stars in the PHAT footprint. Applying our color relations to a subset of the chemically-classified stars producing the bulk of the dust, we find that ~97.8% of the dust is oxygen-rich. Using several scenarios for the dust lifetime, we have estimated the contribution of AGB stars to the global dust budget of M31 to be 0.9-35.5%, which is in line with previous estimates in the Magellanic Clouds. Follow-up observations of the M31 AGB candidates with the JWST will allow us to further constrain stellar and chemical evolutionary models, and the feedback and dust production of metal-rich evolved stars. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJS. 35 pages, 25 Figures

SOFIA/FORCAST Monitoring of the Dust Emission from R Aqr: Start of the Eclipse

Authors: Ravi Sankrit, Eric Omelian, Uma Gorti, R. Mark Wagner, Steven Goldman, Patricia A. Whitelock

Abstract: We present mid-infrared spectra from our continued monitoring of R Aquarii, the nearest symbiotic Mira, using the Stratospheric Observatory for Infrared Astronomy (SOFIA). New photometric and spectroscipic data were obtained with the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST) in 2018 and 2019 after the system had started its "eclipse", during which it became two magnitudes fain… ▽ More We present mid-infrared spectra from our continued monitoring of R Aquarii, the nearest symbiotic Mira, using the Stratospheric Observatory for Infrared Astronomy (SOFIA). New photometric and spectroscipic data were obtained with the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST) in 2018 and 2019 after the system had started its "eclipse", during which it became two magnitudes fainter in the visual. The mid-IR flux, in particular the 10 mu-m silicate feature, have strengthened compared with the previous cycles. Radiative transfer models for the circumstellar dust emission were calculated for the new spectra, and re-calculated for those previously obtained using more appropriate values of the near-IR magnitudes to constrain the properties of the AGB spectra heating the dust. The modeling shows that the luminosity dependence on pulsation phase is not affected by the onset of the eclipse, and that the increase in the mid-IR flux is due to a higher dust density. The models also confirm our earlier results that micron-size grains are present, and that no changes in the grain composition are required to explain the variations in the spectra. △ Less

Submitted 24 December, 2021; originally announced December 2021.

Comments: 20 pages, 4 figures, accepted for publication in the Astrophysical Journal

Deep Learning for Reaction-Diffusion Glioma Growth Modelling: Towards a Fully Personalised Model?

Authors: Corentin Martens, Antonin Rovai, Daniele Bonatto, Thierry Metens, Olivier Debeir, Christine Decaestecker, Serge Goldman, Gaetan Van Simaeys

Abstract: Reaction-diffusion models have been proposed for decades to capture the growth of gliomas, the most common primary brain tumours. However, severe limitations regarding the estimation of the initial conditions and parameter values of such models have restrained their clinical use as a personalised tool. In this work, we investigate the ability of deep convolutional neural networks (DCNNs) to addres… ▽ More Reaction-diffusion models have been proposed for decades to capture the growth of gliomas, the most common primary brain tumours. However, severe limitations regarding the estimation of the initial conditions and parameter values of such models have restrained their clinical use as a personalised tool. In this work, we investigate the ability of deep convolutional neural networks (DCNNs) to address the pitfalls commonly encountered in the field. Based on 1,200 synthetic tumours grown over real brain geometries derived from magnetic resonance (MR) data of 6 healthy subjects, we demonstrate the ability of DCNNs to reconstruct a whole tumour cell density distribution from only two imaging contours at a single time point. With an additional imaging contour extracted at a prior time point, we also demonstrate the ability of DCNNs to accurately estimate the individual diffusivity and proliferation parameters of the model. From this knowledge, the spatio-temporal evolution of the tumour cell density distribution at later time points can ultimately be precisely captured using the model. We finally show the applicability of our approach to MR data of a real glioblastoma patient. This approach may open the perspective of a clinical application of reaction-diffusion growth models for tumour prognosis and treatment planning. △ Less

Submitted 26 November, 2021; originally announced November 2021.

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

Authors: P. Scicluna, F. Kemper, I. McDonald, S. Srinivasan, A. Trejo, S. H. J. Wallström, J. G. A. Wouterloot, J. Cami, J. Greaves, Jinhua He, D. T. Hoai, Hyosun Kim, O. C. Jones, H. Shinnaga, C. J. R. Clark, T. Dharmawardena, W. Holland, H. Imai, J. Th. van Loon, K. M. Menten, R. Wesson, H. Chawner, S. Feng, S. Goldman, F. C. Liu , et al. (67 additional authors not shown)

Abstract: The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a n… ▽ More The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

Machine learning modeling of family wide enzyme-substrate specificity screens

Authors: Samuel Goldman, Ria Das, Kevin K. Yang, Connor W. Coley

Abstract: Biocatalysis is a promising approach to sustainably synthesize pharmaceuticals, complex natural products, and commodity chemicals at scale. However, the adoption of biocatalysis is limited by our ability to select enzymes that will catalyze their natural chemical transformation on non-natural substrates. While machine learning and in silico directed evolution are well-posed for this predictive mod… ▽ More Biocatalysis is a promising approach to sustainably synthesize pharmaceuticals, complex natural products, and commodity chemicals at scale. However, the adoption of biocatalysis is limited by our ability to select enzymes that will catalyze their natural chemical transformation on non-natural substrates. While machine learning and in silico directed evolution are well-posed for this predictive modeling challenge, efforts to date have primarily aimed to increase activity against a single known substrate, rather than to identify enzymes capable of acting on new substrates of interest. To address this need, we curate 6 different high-quality enzyme family screens from the literature that each measure multiple enzymes against multiple substrates. We compare machine learning-based compound-protein interaction (CPI) modeling approaches from the literature used for predicting drug-target interactions. Surprisingly, comparing these interaction-based models against collections of independent (single task) enzyme-only or substrate-only models reveals that current CPI approaches are incapable of learning interactions between compounds and proteins in the current family level data regime. We further validate this observation by demonstrating that our no-interaction baseline can outperform CPI-based models from the literature used to guide the discovery of kinase inhibitors. Given the high performance of non-interaction based models, we introduce a new structure-based strategy for pooling residue representations across a protein sequence. Altogether, this work motivates a principled path forward in order to build and evaluate meaningful predictive models for biocatalysis and other drug discovery applications. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Strong magnon-photon coupling with chip-integrated YIG in the zero-temperature limit

Authors: Paul G. Baity, Dmytro A. Bozhko, Rair Macêdo, William Smith, Rory C. Holland, Sergey Danilin, Valentino Seferai, João Barbosa, Renju R. Peroor, Sara Goldman, Umberto Nasti, Jharna Paul, Robert H. Hadfield, Stephen McVitie, Martin Weides

Abstract: The cross-integration of spin-wave and superconducting technologies is a promising method for creating novel hybrid devices for future information processing technologies to store, manipulate, or convert data in both classical and quantum regimes. Hybrid magnon-polariton systems have been widely studied using bulk Yttrium Iron Garnet (Y$_{3}$Fe$_{5}$O$_{12}$, YIG) and three-dimensional microwave p… ▽ More The cross-integration of spin-wave and superconducting technologies is a promising method for creating novel hybrid devices for future information processing technologies to store, manipulate, or convert data in both classical and quantum regimes. Hybrid magnon-polariton systems have been widely studied using bulk Yttrium Iron Garnet (Y$_{3}$Fe$_{5}$O$_{12}$, YIG) and three-dimensional microwave photon cavities. However, limitations in YIG growth have thus far prevented its incorporation into CMOS compatible technology such as high quality factor superconducting quantum technology. To overcome this impediment, we have used Plasma Focused Ion Beam (PFIB) technology -- taking advantage of precision placement down to the micron-scale -- to integrate YIG with superconducting microwave devices. Ferromagnetic resonance has been measured at millikelvin temperatures on PFIB-processed YIG samples using planar microwave circuits. Furthermore, we demonstrate strong coupling between superconducting resonator and YIG ferromagnetic resonance modes by maintaining reasonably low loss while reducing the system down to the micron scale. This achievement of strong coupling on-chip is a crucial step toward fabrication of functional hybrid quantum devices that advantage from spin-wave and superconducting components. △ Less

Submitted 14 June, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 10 pages, 6 figures

Infrared variable stars in the compact elliptical galaxy M32

Authors: O. C. Jones, C. Nally, M. J. Sharp, I. McDonald, M. L. Boyer, M. Meixner, F. Kemper, A. M. N. Ferguson, S. R. Goldman, R. M. Rich

Abstract: Variable stars in the compact elliptical galaxy M32 are identified, using three epochs of photometry from the Spitzer Space Telescope at 3.6 and 4.5 $μ$m, separated by 32 to 381 days. We present a high-fidelity catalogue of sources detected in multiple epochs at both 3.6 and 4.5 $μ$m, which we analysed for stellar variability using a joint probability error-weighted flux difference. Of these, 83 s… ▽ More Variable stars in the compact elliptical galaxy M32 are identified, using three epochs of photometry from the Spitzer Space Telescope at 3.6 and 4.5 $μ$m, separated by 32 to 381 days. We present a high-fidelity catalogue of sources detected in multiple epochs at both 3.6 and 4.5 $μ$m, which we analysed for stellar variability using a joint probability error-weighted flux difference. Of these, 83 stars are identified as candidate large-amplitude, long-period variables, with 28 considered high-confidence variables. The majority of the variable stars are classified as asymptotic giant branch star candidates using colour-magnitude diagrams. We find no evidence supporting a younger, infrared-bright stellar population in our M32 field. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 11 pages, 8 figures. Accepted for publication in MNRAS

Initial condition assessment for reaction-diffusion glioma growth models: A translational MRI/histology (in)validation study

Authors: Corentin Martens, Laetitia Lebrun, Christine Decaestecker, Thomas Vandamme, Yves-Rémi Van Eycke, Antonin Rovai, Thierry Metens, Olivier Debeir, Serge Goldman, Isabelle Salmon, Gaetan Van Simaeys

Abstract: Diffuse gliomas are highly infiltrative tumors whose early diagnosis and follow-up usually rely on magnetic resonance imaging (MRI). However, the limited sensitivity of this technique makes it impossible to directly assess the extent of the glioma cell invasion, leading to sub-optimal treatment planing. Reaction-diffusion growth models have been proposed for decades to extrapolate glioma cell infi… ▽ More Diffuse gliomas are highly infiltrative tumors whose early diagnosis and follow-up usually rely on magnetic resonance imaging (MRI). However, the limited sensitivity of this technique makes it impossible to directly assess the extent of the glioma cell invasion, leading to sub-optimal treatment planing. Reaction-diffusion growth models have been proposed for decades to extrapolate glioma cell infiltration beyond margins visible on MRI and predict its spatial-temporal evolution. These models nevertheless require an initial condition, that is the tumor cell density values at every location of the brain at diagnosis time. Several works have proposed to relate the tumor cell density function to abnormality outlines visible on MRI but the underlying assumptions have never been verified so far. In this work we propose to verify these assumptions by stereotactic histological analysis of a non-operated brain with glioblastoma using a tailored 3D-printed slicer. Cell density maps are computed from histological slides using a deep learning approach. The density maps are then registered to a postmortem MR image and related to an MR-derived geodesic distance map to the tumor core. The relation between the edema outlines visible on T2 FLAIR MRI and the distance to the core is also investigated. Our results suggest that (i) the previously suggested exponential decrease of the tumor cell density with the distance to the tumor core is not unreasonable but (ii) the edema outlines may in general not correspond to a cell density iso-contour and (iii) the commonly adopted tumor cell density value at these outlines is likely overestimated. These findings highlight the limitations of using conventional MRI to derive glioma cell density maps and point out the need of validating other methods to initialize reaction-diffusion growth models and make them usable in clinical practice. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Voxelwise principal component analysis of dynamic [S-methyl-11C]methionine PET data in glioma patients

Authors: Corentin Martens, Olivier Debeir, Christine Decaestecker, Thierry Metens, Laetitia Lebrun, Gil Leurquin-Sterk, Nicola Trotta, Serge Goldman, Gaetan Van Simaeys

Abstract: Recent works have demonstrated the added value of dynamic amino acid positron emission tomography (PET) for glioma grading and genotyping, biopsy targeting, and recurrence diagnosis. However, most of these studies are exclusively based on hand-crafted qualitative or semi-quantitative dynamic features extracted from the mean time activity curve (TAC) within predefined volumes. Voxelwise dynamic PET… ▽ More Recent works have demonstrated the added value of dynamic amino acid positron emission tomography (PET) for glioma grading and genotyping, biopsy targeting, and recurrence diagnosis. However, most of these studies are exclusively based on hand-crafted qualitative or semi-quantitative dynamic features extracted from the mean time activity curve (TAC) within predefined volumes. Voxelwise dynamic PET data analysis could instead provide a better insight into intra-tumour heterogeneity of gliomas. In this work, we investigate the ability of the widely used principal component analysis (PCA) method to extract meaningful quantitative dynamic features from high-dimensional motion-corrected dynamic [S-methyl-11C]methionine PET data in a first cohort of 20 glioma patients. By means of realistic numerical simulations, we demonstrate the robustness of our methodology to noise. In a second cohort of 13 glioma patients, we compare the resulting parametric maps to these provided by standard one- and two-tissue compartment pharmacokinetic (PK) models. We show that our PCA model outperforms PK models in the identification of intra-tumour uptake dynamics heterogeneity while being much less computationally expensive. Such parametric maps could be valuable to assess tumour aggressiveness locally with applications in treatment planning as well as in the evaluation of tumour progression and response to treatment. This work also provides further encouraging results on the added value of dynamic over static analysis of [S-methyl-11C]methionine PET data in gliomas, as previously demonstrated for O-(2-[18F]fluoroethyl)-L-tyrosine. △ Less

Submitted 20 November, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

The Dusty Evolved Star Kit (DESK): A Python package for fitting the Spectral Energy Distribution of Evolved Stars

Authors: Steven R. Goldman

Abstract: One of the few ways that we can understand the environment around dusty stars and how much material they contribute back to the Universe, is by fitting their brightness at different wavelengths with models that account for how the energy transfers through the dust. The DESK is a python package designed to compare the best fits of different stellar samples and model grids for a better understanding… ▽ More One of the few ways that we can understand the environment around dusty stars and how much material they contribute back to the Universe, is by fitting their brightness at different wavelengths with models that account for how the energy transfers through the dust. The DESK is a python package designed to compare the best fits of different stellar samples and model grids for a better understanding of the results and their uncertainties. The package fits the Spectral Energy Distribution (SED) of evolved stars, using photometry or spectra, to grids of radiative transfer models using a least-squares method. The package includes newly created grids using a variety of different dust species, and state-of-the-art dust growth grids. A robust method for testing different model grids will be particularly important given the wealth of infrared data to come from the James Webb Space Telescope (JWST). △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: Published in JOSS

Journal ref: Journal of Open Source Software, 5(54), 2554 (2020)

PHAT XX. AGB stars and other cool giants in M31 star clusters

Authors: Leo Girardi, Martha L. Boyer, L. Clifton Johnson, Julianne J. Dalcanton, Philip Rosenfield, Anil C. Seth, Evan D. Skillman, Daniel R. Weisz, Benjamin F. Williams, Antara Raaghavi Bhattacharya, Alessandro Bressan, Nelson Caldwell, Yang Chen, Andrew E. Dolphin, Morgan Fouesneau, Steven Goldman, Puragra Guhathakurta, Paola Marigo, Sagnick Mukherjee, Giada Pastorelli, Amanda Quirk, Monika Soraisam, Michele Trabucchi

Abstract: The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of… ▽ More The presence of AGB stars in clusters provides key constraints for stellar models, as has been demonstrated with historical data from the Magellanic Clouds. In this work, we look for candidate AGB stars in M31 star clusters from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. Our photometric criteria selects stars brighter than the tip of the red giant branch, which includes the bulk of the thermally-pulsing AGB stars as well as early-AGB stars and other luminous cool giants expected in young stellar populations (e.g. massive red supergiants, and intermediate-mass red helium-burning stars). The AGB stars can be differentiated, a posteriori, using the ages already estimated for our cluster sample. 937 candidates are found within the cluster aperture radii, half (450) of which are very likely cluster members. Cross-matching with additional databases reveals two carbon stars and ten secure variables among them. The field-corrected age distribution reveals the presence of young supergiants peaking at ages smaller than 100 Myr, followed by a long tail of AGB stars extending up to the oldest possible ages. This long tail reveals the general decrease in the numbers of AGB stars from initial values of 50e-6/Msun at 100 Myr down to 5e-6/Msun at 10 Gyr. Theoretical models of near-solar metallicity reproduce this general trend, although with localized discrepancies over some age intervals, whose origin is not yet identified. The entire catalogue is released together with finding charts to facilitate follow-up studies. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: To appear in ApJ. Data will be soon on MAST

Measuring shape relations using r-parallel sets

Authors: Hans JT Stephensen, Anne Marie Svane, Carlos Benitez, Steven A. Goldman, Jon Sporring

Abstract: Geometrical measurements of biological objects form the basis of many quantitative analyses. Hausdorff measures such as the volume and the area of objects are simple and popular descriptors of individual objects, however, for most biological processes, the interaction between objects cannot be ignored, and the shape and function of neighboring objects are mutually influential. In this paper, we… ▽ More Geometrical measurements of biological objects form the basis of many quantitative analyses. Hausdorff measures such as the volume and the area of objects are simple and popular descriptors of individual objects, however, for most biological processes, the interaction between objects cannot be ignored, and the shape and function of neighboring objects are mutually influential. In this paper, we present a theory on the geometrical interaction between objects based on the theory of spatial point processes. Our theory is based on the relation between two objects: a reference and an observed object. We generate the $r$-parallel sets of the reference object, we calculate the intersection between the $r$-parallel sets and the observed object, and we define measures on these intersections. Our measures are simple like the volume and area of an object, but describe further details about the shape of individual objects and their pairwise geometrical relation. Finally, we propose a summary statistics for collections of shapes and their interaction. We evaluate these measures on a publicly available FIB-SEM 3D data set of an adult rodent. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 10 pages 5 figures

MSC Class: 51K99 (Primary); 92F05 (Secondary) ACM Class: I.4.7

Improved Conditional Flow Models for Molecule to Image Synthesis

Authors: Karren Yang, Samuel Goldman, Wengong Jin, Alex Lu, Regina Barzilay, Tommi Jaakkola, Caroline Uhler

Abstract: In this paper, we aim to synthesize cell microscopy images under different molecular interventions, motivated by practical applications to drug development. Building on the recent success of graph neural networks for learning molecular embeddings and flow-based models for image generation, we propose Mol2Image: a flow-based generative model for molecule to cell image synthesis. To generate cell fe… ▽ More In this paper, we aim to synthesize cell microscopy images under different molecular interventions, motivated by practical applications to drug development. Building on the recent success of graph neural networks for learning molecular embeddings and flow-based models for image generation, we propose Mol2Image: a flow-based generative model for molecule to cell image synthesis. To generate cell features at different resolutions and scale to high-resolution images, we develop a novel multi-scale flow architecture based on a Haar wavelet image pyramid. To maximize the mutual information between the generated images and the molecular interventions, we devise a training strategy based on contrastive learning. To evaluate our model, we propose a new set of metrics for biological image generation that are robust, interpretable, and relevant to practitioners. We show quantitatively that our method learns a meaningful embedding of the molecular intervention, which is translated into an image representation reflecting the biological effects of the intervention. △ Less

Submitted 15 June, 2020; originally announced June 2020.

MSC Class: 92-08

TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Motivations and protocol version 1.0

Authors: Thomas J. Fauchez, Martin Turbet, Eric T. Wolf, Ian Boutle, Michael J. Way, Anthony D. Del Genio, Nathan J. Mayne, Konstantinos Tsigaridis, Ravi K. Kopparapu, Jun Yang, Francois Forget, Avi Mandell, Shawn D. Domagal Goldman

Abstract: Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Tim… ▽ More Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: Accepted and published in Geoscientific Model Development (GMD)

Journal ref: Geosci. Model Dev., 13, 707-716, 2020

Studies of Evolved Stars in the Next Decade: EAO Submillimetre Futures White Paper Series, 2019

Authors: Peter Scicluna, Hiroko Shinnaga, Jonathan Marshall, Jan Wouterloot, Iain McDonald, Steven Goldman, Sofia Wallström, Jacco Th. van Loon, Thavisha Dharmawardena, Lapo Fanciullo, Sundar Srinivasan

Abstract: This white paper discusses recent progress in the field of evolved stars, primarily highlighting the contributions of the James Clerk Maxwell Telescope. It discusses the ongoing project, the \emph{Nearby Evolved Stars Survey} (NESS), and the potential to build upon NESS in the next decade. It then outlines a number of science cases which may become feasible with the proposed 850\,$μ$m camera which… ▽ More This white paper discusses recent progress in the field of evolved stars, primarily highlighting the contributions of the James Clerk Maxwell Telescope. It discusses the ongoing project, the \emph{Nearby Evolved Stars Survey} (NESS), and the potential to build upon NESS in the next decade. It then outlines a number of science cases which may become feasible with the proposed 850\,$μ$m camera which is due to become available at the JCMT in late 2022. These include mapping the extended envelopes of evolved stars, including in polarisation, and time-domain monitoring of their variations. The improved sensitivity of the proposed instrument will facilitate statistical studies that put the morphology, polarisation properties and sub-mm variability in perspective with a relatively modest commitment of time that would be impossible with current instrumentation. We also consider the role that could be played by other continuum wavelengths, heterodyne instruments or other facilities in contributing towards these objectives. △ Less

Submitted 8 February, 2020; originally announced February 2020.

Comments: 11 pages, 2 figures; Part of the 2019 EAO Submillimetre Futures Paper Series, https://www.eaobservatory.org//jcmt/wp-content/uploads/sites/2/2019/11/EAOSubmmFuturesWPs.pdf

Mechanisms of GaN quantum dot formation during nitridation of Ga droplets

Authors: H. Lu, C. Reese, S. Jeon, A. Sundar, Y. Fan, E. Rizzi, Y. Zhuo, L. Qi, R. S. Goldman

Abstract: We have examined the formation mechanisms of GaN quantum dots (QDs) via annealing of Ga droplets in a nitrogen flux. We consider the temperature and substrate dependence of the size distributions of droplets and QDs, as well as the relative roles of Ga/N diffusivity and GaN nucleation rates on QD formation. We report on two competing mechanisms mediated by Ga surface diffusion, namely QD formation… ▽ More We have examined the formation mechanisms of GaN quantum dots (QDs) via annealing of Ga droplets in a nitrogen flux. We consider the temperature and substrate dependence of the size distributions of droplets and QDs, as well as the relative roles of Ga/N diffusivity and GaN nucleation rates on QD formation. We report on two competing mechanisms mediated by Ga surface diffusion, namely QD formation at or away from pre-existing Ga droplets. We discuss the relative roles of nucleation and coarsening dominant growth, as well as the polytype selection, on various substrates. The new insights provide an opportunity for tailoring QD size and polytype distributions for a wide range of III-N semiconductor QDs. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Asymptotic Giant Branch Stars in the Nearby Dwarf Galaxy Leo P

Authors: Steven R. Goldman, Martha L. Boyer, Kristen B. McQuinn, Greg C. Sloan, Iain McDonald, Jacco Th. van Loon, Albert A. Zijlstra, Alec S. Hirschauer, Evan D. Skillman, Sundar Srinivasan

Abstract: We have conducted a highly sensitive census of the evolved-star population in the metal-poor dwarf galaxy Leo P and detected four asymptotic giant branch (AGB) star candidates. Leo P is one of the best examples of a nearby analog of high-redshift galaxies because of its primitive metal content (2% of the solar value), proximity, and isolated nature, ensuring a less complicated history. Using mediu… ▽ More We have conducted a highly sensitive census of the evolved-star population in the metal-poor dwarf galaxy Leo P and detected four asymptotic giant branch (AGB) star candidates. Leo P is one of the best examples of a nearby analog of high-redshift galaxies because of its primitive metal content (2% of the solar value), proximity, and isolated nature, ensuring a less complicated history. Using medium-band optical photometry from the Hubble Space Telescope (HST), we have classified the AGB candidates by their chemical type. We have identified one oxygen-rich source which appears to be dusty in both the HST and Spitzer observations. Its brightness, however, suggests it may be a planetary nebula or post-AGB object. We have also identified three carbon-rich candidates, one of which may be dusty. Follow-up observations are needed to confirm the nature of these sources and to study the composition of any dust that they produce. If dust is confirmed, these stars would likely be among the most metal-poor examples of dust-producing stars known and will provide valuable insight into our understanding of dust formation at high redshift. △ Less

Submitted 6 October, 2019; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: Accepted for publication in ApJ

The Nearby Evolved Stars Survey: I. JCMT/SCUBA-2 Sub-millimetre detection of the detached shell of U Antliae

Authors: Thavisha E. Dharmawardena, Francisca Kemper, Sundar Srinivasan, Peter Scicluna, Jonathan P. Marshall, Jan G. A. Wouterloot, Jane Greaves, Steven R. Goldman, Jacco Th. van Loon, Mikako Matsuura, Iain McDonald, Jinhua He, Albert A. Zijlstra, Jesús A. Toalá, Sofia H. J. Wallström, Hyosun Kim, Alfonso Trejo, Paolo Ventura, Eric Lagadec, Martha L. Boyer, Tie Liu, Gioia Rau, Hideyuki Izumiura, Jan Cami, Wayne Holland , et al. (2 additional authors not shown)

Abstract: We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The… ▽ More We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3σ$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS

The mass-loss, expansion velocities and dust production rates of carbon stars in the Magellanic Clouds

Authors: Ambra Nanni, Martin A. T. Groenewegen, Bernhard Aringer, Stefano Rubele, Alessandro Bressan, Jacco Th. van Loon, Steven R. Goldman, Martha L. Boyer

Abstract: The properties of carbon stars in the Magellanic Clouds (MCs) and their total dust production rates are predicted by fitting their spectral energy distributions (SED) over pre-computed grids of spectra reprocessed by dust. The grids are calculated as a function of the stellar parameters by consistently following the growth for several dust species in their circumstellar envelopes, coupled with a s… ▽ More The properties of carbon stars in the Magellanic Clouds (MCs) and their total dust production rates are predicted by fitting their spectral energy distributions (SED) over pre-computed grids of spectra reprocessed by dust. The grids are calculated as a function of the stellar parameters by consistently following the growth for several dust species in their circumstellar envelopes, coupled with a stationary wind. Dust radiative transfer is computed taking as input the results of the dust growth calculations. The optical constants for amorphous carbon are selected in order to reproduce different observations in the infrared and optical bands of \textit{Gaia} Data Release 2. We find a tail of extreme mass-losing carbon stars in the Large Magellanic Cloud (LMC) with low gas-to-dust ratios that is not present in the Small Magellanic Cloud (SMC). Typical gas-to-dust ratios are around $700$ for the extreme stars, but they can be down to $\sim160$--$200$ and $\sim100$ for a few sources in the SMC and in the LMC, respectively. The total dust production rate for the carbon star population is $\sim 1.77\pm 0.45\times10^{-5}$~M$_\odot$~yr$^{-1}$, for the LMC, and $\sim 2.52\pm 0.96 \times 10^{-6}$~M$_\odot$~yr$^{-1}$, for the SMC. The extreme carbon stars observed with the Atacama Large Millimeter Array and their wind speed are studied in detail. For the most dust-obscured star in this sample the estimated mass-loss rate is $\sim 6.3 \times 10^{-5}$~M$_\odot$~yr$^{-1}$. The grids of spectra are available at: https://ambrananni085.wixsite.com/ambrananni/online-data-1 and included in the SED-fitting python package for fitting evolved stars https://github.com/s-goldman/Dusty-Evolved-Star-Kit . △ Less

Submitted 25 May, 2019; v1 submitted 14 April, 2019; originally announced April 2019.

Comments: 22 pages, 16 figures, accepted for publication in MNRAS

An Infrared Census Of Dust In Nearby Galaxies With Spitzer (DUSTiNGS): V. The Period-luminosity Relation For Dusty Metal-Poor AGB Stars

Authors: Steven R. Goldman, Martha L. Boyer, Kristen B. McQuinn, Patricia A. Whitelock, Iain McDonald, Jacco Th. van Loon, Evan D. Skillman, Robert D. Gehrz, Atefeh Javadi, Gregory C. Sloan, Olivia C. Jones, Martin A. T. Groenewegen, John W. Menzies

Abstract: The survey for DUST In Nearby Galaxies with Spitzer (DUSTiNGS) has identified hundreds of candidate dust-producing Asymptotic Giant Branch (AGB) stars in several nearby metal-poor galaxies. We have obtained multi-epoch follow-up observations for these candidates with the Spitzer Space Telescope and measured their infrared (IR) lightcurves. This has allowed us to confirm their AGB nature and invest… ▽ More The survey for DUST In Nearby Galaxies with Spitzer (DUSTiNGS) has identified hundreds of candidate dust-producing Asymptotic Giant Branch (AGB) stars in several nearby metal-poor galaxies. We have obtained multi-epoch follow-up observations for these candidates with the Spitzer Space Telescope and measured their infrared (IR) lightcurves. This has allowed us to confirm their AGB nature and investigate pulsation behavior at very low metallicity. We have obtained high-confidence pulsation periods for 88 sources in seven galaxies. We have confirmed DUSTiNGS variable star candidates with a 20% success rate, and determined the pulsation properties of 19 sources already identified as Thermally-Pulsing AGB (TP-AGB) stars. We find that the AGB pulsation properties are similar in all galaxies surveyed here, with no discernible difference between the DUSTiNGS galaxies (down to 1.4% solar metallicity; [Fe/H]=-1.85) and the far more metal-rich Magellanic Clouds (up to 50% solar metallicity; [Fe/H]=-0.38). These results strengthen the link between dust production and pulsation in AGB stars and establish the IR Period-Luminosity (P - L) relation as a reliable tool (+/- 4%) for determining distances to galaxies, regardless of metallicity. △ Less

Submitted 22 February, 2019; v1 submitted 19 February, 2019; originally announced February 2019.

Comments: accepted for publication in ApJ, 26 pages, 16 figures, 7 tables

SPIRITS Catalog of Infrared Variables: Identification of Extremely Luminous Long Period Variables

Authors: V. R. Karambelkar, S. M. Adams, P. A. Whitelock, M. M. Kasliwal, J. E. Jencson, M. L. Boyer, S. R. Goldman, F. Masci, A. M. Cody, J. Bally, H. E. Bond, R. D. Gehrz, M. Parthasarathy, R. M. Lau

Abstract: We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud.… ▽ More We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $μ$m photometric catalogs for all sources in these 20 galaxies. △ Less

Submitted 6 May, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 18 pages, 25 figures

Explicit cloud representation in the Atmos 1D climate model for Earth and rocky planet applications

Authors: Thomas Fauchez, Giada Arney, Ravi Kumar Kopparapu, Shawn Domagal Goldman

Abstract: 1D climate models are less sophisticated than 3D global circulation models (GCMs), however their computational time is much less expensive, allowing a large number of runs in a short period of time to explore a wide parameter space. Exploring parameter space is particularly important for predicting the observable properties of exoplanets, for which few parameters are known with certainty. Therefor… ▽ More 1D climate models are less sophisticated than 3D global circulation models (GCMs), however their computational time is much less expensive, allowing a large number of runs in a short period of time to explore a wide parameter space. Exploring parameter space is particularly important for predicting the observable properties of exoplanets, for which few parameters are known with certainty. Therefore, 1D climate models are still very useful tools for planetary studies. In most of these 1D models, clouds are not physically represented in the atmosphere, despite having a well-known, significant impact on a planetary radiative budget. This impact is simulated by artificially raising surface albedo, in order to reproduce the observed-averaged surface temperature (i.e. 288 K for modern Earth) and a radiative balance at the top of the atmosphere. This non-physical representation of clouds, causes atmospheric longwave and shortwaves fluxes to not match observational data. Additionally, this technique represents a parameter that is highly-tuned to modern Earth's climate, and may not be appropriate for planets that deviate from modern Earth's climate conditions. In this paper, we present an update to the climate model within the Atmos 1D atmospheric modeling package with a physical representation of clouds. We show that this physical representation of clouds in the atmosphere allows both longwave and shortwave fluxes to match observational data. This improvement will allow us to study the energy fluxes for a variety of cloudy rocky planets, and increase our confidence in future simulations of temperature profile and net energy balance. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: Accepted in AIMS Geosciences

The planetary nebula IC 5148 and its ionized halo

Authors: D. Barría, S. Kimeswenger, W. Kausch, D. S. Goldman

Abstract: Many round or nearly roundish Planetary Nebulae (PNe) show multiple shells and halo structures during their evolutionary stage near the maximum temperature of their central star. Controversial debate is currently ongoing if these structures are recombination halos, as suggested by hydrodynamic modelling efforts, or ionized material. Recently we discovered a halo with even somewhat unusual structur… ▽ More Many round or nearly roundish Planetary Nebulae (PNe) show multiple shells and halo structures during their evolutionary stage near the maximum temperature of their central star. Controversial debate is currently ongoing if these structures are recombination halos, as suggested by hydrodynamic modelling efforts, or ionized material. Recently we discovered a halo with even somewhat unusual structures around the sparsely studied PN IC~5148 and present for the first time spectroscopy going out to the halo of such a PN.} resolution spectroscopy is used to derive dust chemistry and mineralogy. We investigate the spatial distribution of material and its ionization state from the center of the nebula up to the very outskirts of the halo. We obtained long-slit low resolution spectroscopy (FORS2@VLT) of the nebula in two position angles, which we used to investigate the nebular structure and its halo in the optical range from 450 to 880\,nm. In addition we used medium resolution spectra taken with X-SHOOTER@VLT ranging from 320 nm to 2.4 mu to derive atmospheric parameters for the central star. We obtained the distance and position in the Galaxy from various methods combined with GAIA DR2 data. We also applied Cloudy models to the nebula in order to derive physical parameters of the various regions. We obtained spatially resolved structures and detailed descriptions of the outrunning shock front and a set of unusual halo structures denoted to further shock. The halo structures appears clearly as hot ionized material. Furthermore we derived a reliable photometric value for the central star at a GAIA distance of D=1.3kpc. Considering the large distance $z=1.0$\,kpc from the galactic plane together to its non-circular motion in the galaxy and, a metallicity only slightly below that of typical disk PNe, most likely IC 5148 originates from a thick disk population star. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 12 pages, 17 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 620, A84 (2018)

Cholinergic switch between two types of slow waves in cerebral cortex

Authors: Trang-Anh E. Nghiem, Núria Tort-Colet, Tomasz Górski, Ulisse Ferrari, Shayan Moghimyfiroozabad, Jennifer S. Goldman, Bartosz Teleńczuk, Cristiano Capone, Thierry Bal, Matteo di Volo, Alain Destexhe

Abstract: Sleep slow waves are known to participate in memory consolidation, yet slow waves occurring under anesthesia present no positive effects on memory. Here, we shed light onto this paradox, based on a combination of extracellular recordings in vivo, in vitro, and computational models. We find two types of slow waves, based on analyzing the temporal patterns of successive slow-wave events. The first t… ▽ More Sleep slow waves are known to participate in memory consolidation, yet slow waves occurring under anesthesia present no positive effects on memory. Here, we shed light onto this paradox, based on a combination of extracellular recordings in vivo, in vitro, and computational models. We find two types of slow waves, based on analyzing the temporal patterns of successive slow-wave events. The first type is consistently observed in natural slow-wave sleep, while the second is shown to be ubiquitous under anesthesia. Network models of spiking neurons predict that the two slow wave types emerge due to a different gain on inhibitory vs excitatory cells and that different levels of spike-frequency adaptation in excitatory cells can account for dynamical distinctions between the two types. This prediction was tested in vitro by varying adaptation strength using an agonist of acetylcholine receptors, which demonstrated a neuromodulatory switch between the two types of slow waves. Finally, we show that the first type of slow-wave dynamics is more sensitive to external stimuli, which can explain how slow waves in sleep and anesthesia differentially affect memory consolidation, as well as provide a link between slow-wave dynamics and memory diseases. △ Less

Submitted 28 November, 2019; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 37 pages, 5 main figures, 4 supplementary figures

Journal ref: Cerebral Cortex 2019

Why do nanowires grow with their c-axis vertically-aligned in the absence of epitaxy?

Authors: Almog R. Azulay, Yury Turkulets, Davide Del Gaudio, Rachel S. Goldman, Ilan Shalish

Abstract: Images of uniform and upright nanowires are fascinating, but often, they are quite puzzling, when epitaxial templating from the substrate is clearly absent. Here, we reveal the physics underlying one such hidden growth guidance mechanism through a specific example - the case of ZnO nanowires grown on silicon oxide and glass. We show how electric fields exerted by the insulating substrate may be ma… ▽ More Images of uniform and upright nanowires are fascinating, but often, they are quite puzzling, when epitaxial templating from the substrate is clearly absent. Here, we reveal the physics underlying one such hidden growth guidance mechanism through a specific example - the case of ZnO nanowires grown on silicon oxide and glass. We show how electric fields exerted by the insulating substrate may be manipulated through the surface charge to define the orientation and polarity of the nanowires. Surface charge is ubiquitous on the surfaces of semiconductors and insulators, and as a result, substrate electric fields need always be considered. Our results suggest a new concept, according to which the growth of wurtzite semiconductors may often be described as a process of electric-charge-induced self assembly, wherein the internal built-in field in the polar material tends to align in parallel to an external field exerted by the substrate to minimize the interfacial energy of the system. △ Less

Submitted 4 November, 2019; v1 submitted 6 September, 2018; originally announced September 2018.

Comments: Revised title, abstract, and discussion

Journal ref: Scientific Reports, 10, 6554 (2020)

State equation from the spectral structure of human brain activity

Authors: Trang-Anh Nghiem, Jean-Marc Lina, Matteo di Volo, Cristiano Capone, Alan C. Evans, Alain Destexhe, Jennifer S. Goldman

Abstract: Neural electromagnetic (EM) signals recorded non-invasively from individual human subjects vary in complexity and magnitude. Nonetheless, variation in neural activity has been difficult to quantify and interpret, due to complex, broad-band features in the frequency domain. Studying signals recorded with magnetoencephalography (MEG) from healthy young adult subjects while in resting and active stat… ▽ More Neural electromagnetic (EM) signals recorded non-invasively from individual human subjects vary in complexity and magnitude. Nonetheless, variation in neural activity has been difficult to quantify and interpret, due to complex, broad-band features in the frequency domain. Studying signals recorded with magnetoencephalography (MEG) from healthy young adult subjects while in resting and active states, a systematic framework inspired by thermodynamics is applied to neural EM signals. Despite considerable inter-subject variation in terms of spectral entropy and energy across time epochs, data support the existence of a robust and linear relationship defining an effective state equation, with higher energy and lower entropy in the resting state compared to active, consistently across subjects. Mechanisms underlying the emergence of relationships between empirically measured effective state functions are further investigated using a model network of coupled oscillators, suggesting an interplay between noise and coupling strength can account for coherent variation of empirically observed quantities. Taken together, the results show macroscopic neural observables follow a robust, non-trivial conservation rule for energy modulation and information generation. △ Less

Submitted 3 July, 2018; v1 submitted 19 June, 2018; originally announced June 2018.

Comments: 6 pages, 4 figures

Synchronization of electrically coupled resonate-and-fire neurons

Authors: Thomas Chartrand, Mark S. Goldman, Timothy J. Lewis

Abstract: Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with strong electrical coupling have been shown to exhibit resonant properties, and the subthreshold fluctuations arising from resonance are transmitted through elec… ▽ More Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with strong electrical coupling have been shown to exhibit resonant properties, and the subthreshold fluctuations arising from resonance are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with linear subthreshold dynamics and discrete post-spike reset. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuity in the dynamics. We find that both spikes and resonant subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs. △ Less

Submitted 17 January, 2018; originally announced January 2018.

MSC Class: 92B25; 34D06; 34C15; 34C20; 92C20; 34A38; 37G15; 37N25

Emission line galaxies behind the planetary nebula IC 5148: Potential for a serendipity survey with archival data

Authors: S. Kimeswenger, D. Barria, W. Kausch, D. S. Goldman

Abstract: During the start of a survey program using FORS2 long slit spectroscopy on planetary nebulae (PN) and their haloes, we serendipitously discovered six background emission line galaxies (ELG) with redshifts of z = 0.2057, 0.3137, 0.37281, 0.4939, 0.7424 and 0.8668. Thus they clearly do not belong to a common cluster structure. We derived the major physical properties of the targets. Since the used l… ▽ More During the start of a survey program using FORS2 long slit spectroscopy on planetary nebulae (PN) and their haloes, we serendipitously discovered six background emission line galaxies (ELG) with redshifts of z = 0.2057, 0.3137, 0.37281, 0.4939, 0.7424 and 0.8668. Thus they clearly do not belong to a common cluster structure. We derived the major physical properties of the targets. Since the used long slit covers a sky area of only 570 arcsec^2, we discuss further potential of serendipitous discoveries in archival data, beside the deep systematic work of the ongoing and upcoming big surveys. We conclude that archival data provide a decent potential for extending the overall data on ELGs without any selection bias. △ Less

Submitted 9 October, 2017; originally announced October 2017.

Comments: 9 pages, 10 figures, accepted for publication in Revista Mexicana de Astronomía y Astrofísica

A dearth of OH/IR stars in the Small Magellanic Cloud

Authors: Steven R. Goldman, Jacco Th. van Loon, José F. Gómez, James A. Green, Albert A. Zijlstra, Ambra Nanni, Hiroshi Imai, Patricia A. Whitelock, Martin A. T. Groenewegen, Joana M. Oliveira

Abstract: We present the results of targeted observations and a survey of 1612-, 1665-, and 1667-MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Small Magellanic Cloud (SMC), using the Parkes and Australia Telescope Compact Array radio telescopes. No clear OH maser emission has been detected in any of our observations targeting luminous, long-pe… ▽ More We present the results of targeted observations and a survey of 1612-, 1665-, and 1667-MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Small Magellanic Cloud (SMC), using the Parkes and Australia Telescope Compact Array radio telescopes. No clear OH maser emission has been detected in any of our observations targeting luminous, long-period, large-amplitude variable stars, which have been confirmed spectroscopically and photometrically to be mid- to late-M spectral type. These observations have probed 3 - 4 times deeper than any OH maser survey in the SMC. Using a bootstrapping method with LMC and Galactic OH/IR star samples and our SMC observation upper limits, we have calculated the likelihood of not detecting maser emission in any of the two sources considered to be the top maser candidates to be less than 0.05%, assuming a similar pumping mechanism as the LMC and Galactic OH/IR sources. We have performed a population comparison of the Magellanic Clouds and used Spitzer IRAC and MIPS photometry to confirm that we have observed all high luminosity SMC sources that are expected to exhibit maser emission. We suspect that, compared to the OH/IR stars in the Galaxy and LMC, the reduction in metallicity may curtail the dusty wind phase at the end of the evolution of the most massive cool stars. We also suspect that the conditions in the circumstellar envelope change beyond a simple scaling of abundances and wind speed with metallicity. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Fast and Precise Type Checking for JavaScript

Authors: Avik Chaudhuri, Panagiotis Vekris, Sam Goldman, Marshall Roch, Gabriel Levi

Abstract: In this paper we present the design and implementation of Flow, a fast and precise type checker for JavaScript that is used by thousands of developers on millions of lines of code at Facebook every day. Flow uses sophisticated type inference to understand common JavaScript idioms precisely. This helps it find non-trivial bugs in code and provide code intelligence to editors without requiring signi… ▽ More In this paper we present the design and implementation of Flow, a fast and precise type checker for JavaScript that is used by thousands of developers on millions of lines of code at Facebook every day. Flow uses sophisticated type inference to understand common JavaScript idioms precisely. This helps it find non-trivial bugs in code and provide code intelligence to editors without requiring significant rewriting or annotations from the developer. We formalize an important fragment of Flow's analysis and prove its soundness. Furthermore, Flow uses aggressive parallelization and incrementalization to deliver near-instantaneous response times. This helps it avoid introducing any latency in the usual edit-refresh cycle of rapid JavaScript development. We describe the algorithms and systems infrastructure that we built to scale Flow's analysis. △ Less

Submitted 30 August, 2017; v1 submitted 26 August, 2017; originally announced August 2017.

TAPAS: Two-pass Approximate Adaptive Sampling for Softmax

Authors: Yu Bai, Sally Goldman, Li Zhang

Abstract: TAPAS is a novel adaptive sampling method for the softmax model. It uses a two pass sampling strategy where the examples used to approximate the gradient of the partition function are first sampled according to a squashed population distribution and then resampled adaptively using the context and current model. We describe an efficient distributed implementation of TAPAS. We show, on both syntheti… ▽ More TAPAS is a novel adaptive sampling method for the softmax model. It uses a two pass sampling strategy where the examples used to approximate the gradient of the partition function are first sampled according to a squashed population distribution and then resampled adaptively using the context and current model. We describe an efficient distributed implementation of TAPAS. We show, on both synthetic data and a large real dataset, that TAPAS has low computational overhead and works well for minimizing the rank loss for multi-class classification problems with a very large label space. △ Less

Submitted 14 July, 2017; v1 submitted 10 July, 2017; originally announced July 2017.

Current-induced spin polarization in InGaAs and GaAs epilayers with varying doping densities

Authors: M. Luengo-Kovac, S. Huang, D. Del Gaudio, J. Occena, R. S. Goldman, R. Raimondi, V. Sih

Abstract: The current-induced spin polarization and momentum-dependent spin-orbit field were measured in In$_{x}$Ga$_{1-x}$As epilayers with varying indium concentrations and silicon doping densities. Samples with higher indium concentrations and carrier concentrations and lower mobilities were found to have larger electrical spin generation efficiencies. Furthermore, current-induced spin polarization was d… ▽ More The current-induced spin polarization and momentum-dependent spin-orbit field were measured in In$_{x}$Ga$_{1-x}$As epilayers with varying indium concentrations and silicon doping densities. Samples with higher indium concentrations and carrier concentrations and lower mobilities were found to have larger electrical spin generation efficiencies. Furthermore, current-induced spin polarization was detected in GaAs epilayers despite the absence of measurable spin-orbit fields, indicating that the extrinsic contributions to the spin polarization mechanism must be considered. Theoretical calculations based on a model that includes extrinsic contributions to the spin dephasing and the spin Hall effect, in addition to the intrinsic Rashba and Dresselhaus spin-orbit coupling, are found to qualitatively agree with the experimental results. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Comments: 16 pages, 8 figures

Journal ref: Phys. Rev. B 96, 195206 (2017)

Astrometry of OH/IR stars using 1612 MHz hydroxyl masers. I. Annual parallaxes of WX Psc and OH138.0+7.2

Authors: G. Orosz, H. Imai, R. Dodson, M. J. Rioja, S. Frey, R. A. Burns, S. Etoka, A. Nakagawa, H. Nakanishi, Y. Asaki, S. R. Goldman, D. Tafoya

Abstract: We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtained a 3-sigma upper limit of <=5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of >~190 pc. The obtained parallax of… ▽ More We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtained a 3-sigma upper limit of <=5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of >~190 pc. The obtained parallax of OH138.0+7.2 is 0.52+/-0.09 mas (+/-18%), corresponding to a distance of 1.9(+0.4,-0.3) kpc, making this the first hydroxyl maser parallax below one milliarcsecond. We also introduce a new method of error analysis for detecting systematic errors in the astrometry. Finally, we compare our trigonometric distances to published phase-lag distances toward these stars and find a good agreement between the two methods. △ Less

Submitted 18 January, 2017; originally announced January 2017.

Comments: Preprint, accepted for publication in The Astronomical Journal (January 17, 2017)

The wind speeds, dust content, and mass-loss rates of evolved AGB and RSG stars at varying metallicity

Authors: Steven R. Goldman, Jacco Th. van Loon, Albert A. Zijlstra, James A. Green, Peter R. Wood, Ambra Nanni, Hiroshi Imai, Patricia A. Whitelock, Mikako Matsuura, Martin A. T. Groenewegen, José F. Gómez

Abstract: We present the results of our survey of 1612 MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Large Magellanic Cloud. We have discovered four new circumstellar maser sources in the LMC, and increased the number of reliable wind speeds from IR stars in the LMC from 5 to 13. Using our new wind speeds, as well as those from Galactic source… ▽ More We present the results of our survey of 1612 MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Large Magellanic Cloud. We have discovered four new circumstellar maser sources in the LMC, and increased the number of reliable wind speeds from IR stars in the LMC from 5 to 13. Using our new wind speeds, as well as those from Galactic sources, we have derived an updated relation for dust driven winds: $v_{exp} \propto Z L^{0.4}$. We compare the sub-solar metallicity LMC OH/IR stars with carefully selected samples of more metal-rich OH/IR stars, also at known distances, in the Galactic Centre and Galactic Bulge. For 8 of the Bulge stars we derive pulsation periods for the first time, using near-IR photometry from the VVV survey. We have modeled our LMC OH/IR stars and developed an empirical method of deriving gas-to-dust ratios and mass loss rates by scaling the models to the results from maser profiles. We have done this also for samples in the Galactic Centre and Bulge and derived a new mass loss prescription that includes luminosity, pulsation period, and gas-to-dust ratio $\dot{M} = 1.06^{+3.5}_{-0.8} \rm{ \cdot }10^{-5}\,(L/10^4\,\rm{L}_\odot)^{0.9\pm0.1}(P/500\,\rm{d})^{0.75\pm0.3} (r_{gd}/200)^{-0.03\pm0.07}\,\rm{M_{\odot}}\, yr^{-1}$. The tightest correlation is found between mass loss rate and luminosity. We find that the gas-to-dust ratio has little effect on the mass loss of oxygen-rich AGB stars and RSGs within the Galaxy and the LMC. This suggests that mass loss of oxygen-rich AGB stars and RSGs is (nearly) independent of metallicity between a half and twice solar. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: 33 pages, 29 figures, accepted for publication in MNRAS

The ALMA detection of CO rotational line emission in AGB stars in the Large Magellanic Cloud

Authors: M. A. T. Groenewegen, W. H. T. Vlemmings, P. Marigo, G. C. Sloan, L. Decin, M. W. Feast, S. R. Goldman, K. Justtanont, F. Kerschbaum, M. Matsuura, I. McDonald, H. Olofsson, R. Sahai, J. Th. van Loon, P. R. Wood, A. A. Zijlstra, J. Bernard-Salas, M. L. Boyer, L. Guzman-Ramirez, O. C. Jones, E. Lagadec, M. Meixner, M. G. Rawlings, S. Srinivasan

Abstract: Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB… ▽ More Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models. △ Less

Submitted 30 September, 2016; originally announced September 2016.

Comments: Accepted for publication by A&A

Journal ref: A&A 596, A50 (2016)

ALMA reveals sunburn: CO dissociation around AGB stars in the globular cluster 47 Tucanae

Authors: Iain McDonald, Albert A. Zijlstra, Eric Lagadec, Gregory C. Sloan, Martha L. Boyer, Mikako Matsuura, Rowan J. Smith, Christina L. Smith, Jeremy A. Yates, Jacco Th. van Loon, Olivia C. Jones, Sofia Ramstedt, Adam Avison, Kay Justtanont, Hans Olofsson, Joris A. D. L. Blommaert, Steven R. Goldman, Martin A. T. Groenewegen

Abstract: ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultra… ▽ More ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultraviolet radiation field from post-AGB stars and white dwarfs in 47 Tuc, we conclude CO should be dissociated abnormally close to the stars. We estimate that the CO envelopes will be truncated at a few hundred stellar radii from their host stars and that the line intensities are about two orders of magnitude below our current detection limits. The truncation of CO envelopes should be important for AGB stars in dense clusters. Observing the CO (3-2) and higher transitions and targeting stars far from the centres of clusters should result in the detections needed to measure the outflow velocities from these stars. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: 14 pages, 4 figures, accepted MNRAS

3D super-resolved in vitro multiphoton microscopy by saturation of excitation

Authors: Anh Dung Nguyen, François Duport, Arno Bouwens, Frédérique Vanholsbeeck, Dominique Egrise, Gaetan Van Simayes, Philippe Emplit, Serge Goldman, Simon-Pierre Gorza

Abstract: We demonstrate a significant resolution enhancement beyond the conventional limit in multiphoton microscopy (MPM) using saturated excitation of fluorescence. Our technique achieves super-resolved imaging by temporally modulating the excitation laser-intensity and demodulating the higher harmonics from the saturated fluorescence signal. The improvement of the lateral and axial resolutions is measur… ▽ More We demonstrate a significant resolution enhancement beyond the conventional limit in multiphoton microscopy (MPM) using saturated excitation of fluorescence. Our technique achieves super-resolved imaging by temporally modulating the excitation laser-intensity and demodulating the higher harmonics from the saturated fluorescence signal. The improvement of the lateral and axial resolutions is measured on a sample of fluorescent microspheres. While the third harmonic already provides an enhanced resolution, we show that a further improvement can be obtained with an appropriate linear combination of the demodulated harmonics. Finally, we present in vitro imaging of fluorescent microspheres incorporated in HeLa cells to show that this technique performs well in biological samples. △ Less

Submitted 10 August, 2015; v1 submitted 6 July, 2015; originally announced July 2015.

g-Factor Modification in a Bulk InGaAs Epilayer by an In-plane Electric Field

Authors: M. Luengo-Kovac, M. Macmahon, S. Huang, R. S. Goldman, V. Sih

Abstract: We report on the modification of the g-factor by an in-plane electric field in an In$_{0.031}$Ga$_{0.969}$As epilayer. We performed external magnetic field scans of the Kerr rotation of the InGaAs film in order to independently determine the g-factor and the spin-orbit fields. The g-factor increases from $-0.4473\pm0.0001$ at 0 V/cm to $-0.4419\pm0.0001$ at 25 V/cm applied along the [1… ▽ More We report on the modification of the g-factor by an in-plane electric field in an In$_{0.031}$Ga$_{0.969}$As epilayer. We performed external magnetic field scans of the Kerr rotation of the InGaAs film in order to independently determine the g-factor and the spin-orbit fields. The g-factor increases from $-0.4473\pm0.0001$ at 0 V/cm to $-0.4419\pm0.0001$ at 25 V/cm applied along the [1$\overline{1}$0] crystal axis. In addition, spatially-resolved spin measurements show a g-factor dependence on diffusive velocity. The change in g-factor with electric field can have a large effect on the determination of the internal spin-orbit and nuclear fields from Larmor precession frequency measurements. △ Less

Submitted 31 March, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 5 pages, 4 figures

Gas bubble dynamics in soft materials

Authors: J. M. Solano-Altamirano, John D. Malcolm, Saul Goldman

Abstract: Epstein and Plesset's seminal work on the rate of gas bubble dissolution and growth in a simple liquid is generalized to render it applicable to a gas bubble embedded in a soft elastic medium. Both the underlying diffusion equation and the expression for the gas bubble pressure were modified to allow for the non-zero shear modulus of the elastic medium. The extension of the diffusion equation resu… ▽ More Epstein and Plesset's seminal work on the rate of gas bubble dissolution and growth in a simple liquid is generalized to render it applicable to a gas bubble embedded in a soft elastic medium. Both the underlying diffusion equation and the expression for the gas bubble pressure were modified to allow for the non-zero shear modulus of the elastic medium. The extension of the diffusion equation results in a trivial shift (by an additive constant) in the value of the diffusion coefficient, and does not change the form of the rate equations. But the use of a Generalized Young-Laplace equation for the bubble pressure resulted in significant differences on the dynamics of bubble dissolution and growth, relative to a simple liquid medium. Depending on whether the salient parameters (solute concentration, initial bubble radius, surface tension, and shear modulus) lead to bubble growth or dissolution, the effect of allowing for a non-zero shear modulus in the Generalized Young-Laplace equation is to speed up the rate of bubble growth, or to reduce the rate of bubble dissolution, respectively. The relation to previous work on visco-elastic materials is discussed, as is the connection of this work to the problem of Decompression Sickness (specifically, "the bends"). Examples of tissues to which our expressions can be applied are provided. Also, a new phenomenon is predicted whereby, for some parameter values, a bubble can be metastable and persist for long times, or it may grow, when embedded in a homogeneous under-saturated soft elastic medium. △ Less

Submitted 13 October, 2014; v1 submitted 10 September, 2014; originally announced September 2014.

Comments: 18 pages, 4 figures, 1 table, (Included also "Supplementary information": 2 pages, 1 table)

Journal ref: Soft Matter, 11 (2015) 202