Showing 1–50 of 52 results for author: Douglas, K

Visual Episodic Memory-based Exploration

Authors: Jack Vice, Natalie Ruiz-Sanchez, Pamela K. Douglas, Gita Sukthankar

Abstract: In humans, intrinsic motivation is an important mechanism for open-ended cognitive development; in robots, it has been shown to be valuable for exploration. An important aspect of human cognitive development is $\textit{episodic memory}$ which enables both the recollection of events from the past and the projection of subjective future. This paper explores the use of visual episodic memory as a so… ▽ More In humans, intrinsic motivation is an important mechanism for open-ended cognitive development; in robots, it has been shown to be valuable for exploration. An important aspect of human cognitive development is $\textit{episodic memory}$ which enables both the recollection of events from the past and the projection of subjective future. This paper explores the use of visual episodic memory as a source of intrinsic motivation for robotic exploration problems. Using a convolutional recurrent neural network autoencoder, the agent learns an efficient representation for spatiotemporal features such that accurate sequence prediction can only happen once spatiotemporal features have been learned. Structural similarity between ground truth and autoencoder generated images is used as an intrinsic motivation signal to guide exploration. Our proposed episodic memory model also implicitly accounts for the agent's actions, motivating the robot to seek new interactive experiences rather than just areas that are visually dissimilar. When guiding robotic exploration, our proposed method outperforms the Curiosity-driven Variational Autoencoder (CVAE) at finding dynamic anomalies. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: FLAIRS 2023, 7 pages, 11 figures

MSC Class: I.2.9

Journal ref: The International FLAIRS Conference Proceedings. Vol. 36. 2023

The impact of generative artificial intelligence on socioeconomic inequalities and policy making

Authors: Valerio Capraro, Austin Lentsch, Daron Acemoglu, Selin Akgun, Aisel Akhmedova, Ennio Bilancini, Jean-François Bonnefon, Pablo Brañas-Garza, Luigi Butera, Karen M. Douglas, Jim A. C. Everett, Gerd Gigerenzer, Christine Greenhow, Daniel A. Hashimoto, Julianne Holt-Lunstad, Jolanda Jetten, Simon Johnson, Chiara Longoni, Pete Lunn, Simone Natale, Iyad Rahwan, Neil Selwyn, Vivek Singh, Siddharth Suri, Jennifer Sutcliffe , et al. (6 additional authors not shown)

Abstract: Generative artificial intelligence has the potential to both exacerbate and ameliorate existing socioeconomic inequalities. In this article, we provide a state-of-the-art interdisciplinary overview of the potential impacts of generative AI on (mis)information and three information-intensive domains: work, education, and healthcare. Our goal is to highlight how generative AI could worsen existing i… ▽ More Generative artificial intelligence has the potential to both exacerbate and ameliorate existing socioeconomic inequalities. In this article, we provide a state-of-the-art interdisciplinary overview of the potential impacts of generative AI on (mis)information and three information-intensive domains: work, education, and healthcare. Our goal is to highlight how generative AI could worsen existing inequalities while illuminating how AI may help mitigate pervasive social problems. In the information domain, generative AI can democratize content creation and access, but may dramatically expand the production and proliferation of misinformation. In the workplace, it can boost productivity and create new jobs, but the benefits will likely be distributed unevenly. In education, it offers personalized learning, but may widen the digital divide. In healthcare, it might improve diagnostics and accessibility, but could deepen pre-existing inequalities. In each section we cover a specific topic, evaluate existing research, identify critical gaps, and recommend research directions, including explicit trade-offs that complicate the derivation of a priori hypotheses. We conclude with a section highlighting the role of policymaking to maximize generative AI's potential to reduce inequalities while mitigating its harmful effects. We discuss strengths and weaknesses of existing policy frameworks in the European Union, the United States, and the United Kingdom, observing that each fails to fully confront the socioeconomic challenges we have identified. We propose several concrete policies that could promote shared prosperity through the advancement of generative AI. This article emphasizes the need for interdisciplinary collaborations to understand and address the complex challenges of generative AI. △ Less

Submitted 6 May, 2024; v1 submitted 16 December, 2023; originally announced January 2024.

Comments: PNAS Nexus, in press

Predicting Learning Interactions in Social Learning Networks: A Deep Learning Enabled Approach

Authors: Rajeev Sahay, Serena Nicoll, Minjun Zhang, Tsung-Yen Yang, Carlee Joe-Wong, Kerrie A. Douglas, Christopher G Brinton

Abstract: We consider the problem of predicting link formation in Social Learning Networks (SLN), a type of social network that forms when people learn from one another through structured interactions. While link prediction has been studied for general types of social networks, the evolution of SLNs over their lifetimes coupled with their dependence on which topics are being discussed presents new challenge… ▽ More We consider the problem of predicting link formation in Social Learning Networks (SLN), a type of social network that forms when people learn from one another through structured interactions. While link prediction has been studied for general types of social networks, the evolution of SLNs over their lifetimes coupled with their dependence on which topics are being discussed presents new challenges for this type of network. To address these challenges, we develop a series of autonomous link prediction methodologies that utilize spatial and time-evolving network architectures to pass network state between space and time periods, and that models over three types of SLN features updated in each period: neighborhood-based (e.g., resource allocation), path-based (e.g., shortest path), and post-based (e.g., topic similarity). Through evaluation on six real-world datasets from Massive Open Online Course (MOOC) discussion forums and from Purdue University, we find that our method obtains substantial improvements over Bayesian models, linear classifiers, and graph neural networks, with AUCs typically above 0.91 and reaching 0.99 depending on the dataset. Our feature importance analysis shows that while neighborhood and path-based features contribute the most to the results, post-based features add additional information that may not always be relevant for link prediction. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: This work was published in the IEEE/ACM Transactions on Networking

Multi-Layer Personalized Federated Learning for Mitigating Biases in Student Predictive Analytics

Authors: Yun-Wei Chu, Seyyedali Hosseinalipour, Elizabeth Tenorio, Laura Cruz, Kerrie Douglas, Andrew Lan, Christopher Brinton

Abstract: Conventional methods for student modeling, which involve predicting grades based on measured activities, struggle to provide accurate results for minority/underrepresented student groups due to data availability biases. In this paper, we propose a Multi-Layer Personalized Federated Learning (MLPFL) methodology that optimizes inference accuracy over different layers of student grouping criteria, su… ▽ More Conventional methods for student modeling, which involve predicting grades based on measured activities, struggle to provide accurate results for minority/underrepresented student groups due to data availability biases. In this paper, we propose a Multi-Layer Personalized Federated Learning (MLPFL) methodology that optimizes inference accuracy over different layers of student grouping criteria, such as by course and by demographic subgroups within each course. In our approach, personalized models for individual student subgroups are derived from a global model, which is trained in a distributed fashion via meta-gradient updates that account for subgroup heterogeneity while preserving modeling commonalities that exist across the full dataset. The evaluation of the proposed methodology considers case studies of two popular downstream student modeling tasks, knowledge tracing and outcome prediction, which leverage multiple modalities of student behavior (e.g., visits to lecture videos and participation on forums) in model training. Experiments on three real-world online course datasets show significant improvements achieved by our approach over existing student modeling benchmarks, as evidenced by an increased average prediction quality and decreased variance across different student subgroups. Visual analysis of the resulting students' knowledge state embeddings confirm that our personalization methodology extracts activity patterns clustered into different student subgroups, consistent with the performance enhancements we obtain over the baselines. △ Less

Submitted 28 May, 2024; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: IEEE Transactions on Emerging Topics in Computing, 2024

Triadic Temporal Exponential Random Graph Models (TTERGM)

Authors: Yifan Huang, Clayton Barham, Eric Page, Pamela K Douglas

Abstract: Temporal exponential random graph models (TERGM) are powerful statistical models that can be used to infer the temporal pattern of edge formation and elimination in complex networks (e.g., social networks). TERGMs can also be used in a generative capacity to predict longitudinal time series data in these evolving graphs. However, parameter estimation within this framework fails to capture many rea… ▽ More Temporal exponential random graph models (TERGM) are powerful statistical models that can be used to infer the temporal pattern of edge formation and elimination in complex networks (e.g., social networks). TERGMs can also be used in a generative capacity to predict longitudinal time series data in these evolving graphs. However, parameter estimation within this framework fails to capture many real-world properties of social networks, including: triadic relationships, small world characteristics, and social learning theories which could be used to constrain the probabilistic estimation of dyadic covariates. Here, we propose triadic temporal exponential random graph models (TTERGM) to fill this void, which includes these hierarchical network relationships within the graph model. We represent social network learning theory as an additional probability distribution that optimizes Markov chains in the graph vector space. The new parameters are then approximated via Monte Carlo maximum likelihood estimation. We show that our TTERGM model achieves improved fidelity and more accurate predictions compared to several benchmark methods on GitHub network data. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Journal ref: 36th Conference on Neural Information Processing Systems (NeurIPS) 2022 Temporal Graph Learning Workshop

Structure in the Magnetic Field of the Milky Way Disk and Halo traced by Faraday Rotation

Authors: John M. Dickey, Jennifer West, Alec J. M. Thomson, T. L. Landecker, A. Bracco, E. Carretti, J. L. Han, A. S. Hill, Y. K. Ma, S. A. Mao, A. Ordog, Jo-Anne C. Brown, K. A. Douglas, A. Erceg, V. Jelic, R. Kothes, M. Wolleben

Abstract: Magnetic fields in the ionized medium of the disk and halo of the Milky Way impose Faraday rotation on linearly polarized radio emission. We compare two surveys mapping the Galactic Faraday rotation, one showing the rotation measures of extragalactic sources seen through the Galaxy (from Hutschenreuter et al 2022), and one showing the Faraday depth of the diffuse Galactic synchrotron emission from… ▽ More Magnetic fields in the ionized medium of the disk and halo of the Milky Way impose Faraday rotation on linearly polarized radio emission. We compare two surveys mapping the Galactic Faraday rotation, one showing the rotation measures of extragalactic sources seen through the Galaxy (from Hutschenreuter et al 2022), and one showing the Faraday depth of the diffuse Galactic synchrotron emission from the Global Magneto-Ionic Medium Survey. Comparing the two data sets in 5deg x 10deg bins shows good agreement at intermediate latitudes, 10 < |b| < 50 deg, and little correlation between them at lower and higher latitudes. Where they agree, both tracers show clear patterns as a function of Galactic longitude: in the Northern Hemisphere a strong sin(2 x longitude) pattern, and in the Southern hemisphere a sin(longitude + pi) pattern. Pulsars with height above or below the plane |z| > 300 pc show similar longitude dependence in their rotation measures. Nearby non-thermal structures show rotation measure shadows as does the Orion-Eridanus superbubble. We describe families of dynamo models that could explain the observed patterns in the two hemispheres. We suggest that a field reversal, known to cross the plane a few hundred pc inside the solar circle, could shift to positive z with increasing Galactic radius to explain the sin(2xlongitude) pattern in the Northern Hemisphere. Correlation shows that rotation measures from extragalactic sources are one to two times the corresponding rotation measure of the diffuse emission, implying Faraday complexity along some lines of sight, especially in the Southern hemisphere. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 37 pages, 26 figures, Ap. J. accepted

The gas-phase reaction of NH2 with formaldehyde (CH2O) is not a source of formamide (NH2CHO) in interstellar environments

Authors: Kevin M. Douglas, Daniel Lucas, Catherine Walsh, Niclas A. West, Mark A. Blitz, Dwayne E. Heard

Abstract: The first experimental study of the low-temperature kinetics of the gas-phase reaction of NH2 with formaldehyde (CH2O) has been performed. This reaction has previously been suggested as a source of formamide (NH2CHO) in interstellar environments. A pulsed Laval nozzle equipped with laser-flash photolysis and laser-induced fluorescence spectroscopy was used to create and monitor the temporal decay… ▽ More The first experimental study of the low-temperature kinetics of the gas-phase reaction of NH2 with formaldehyde (CH2O) has been performed. This reaction has previously been suggested as a source of formamide (NH2CHO) in interstellar environments. A pulsed Laval nozzle equipped with laser-flash photolysis and laser-induced fluorescence spectroscopy was used to create and monitor the temporal decay of NH2 in the presence of CH2O. No loss of NH2 could be observed via reaction with CH2O and we place an upper-limit on the rate coefficient of <6x10-12 cm3 molecule-1 s-1 at 34K. Ab initio calculations of the potential energy surface were combined with RRKM calculations to predict a rate coefficient of 6.2x10-14 cm3 molecule-1 s-1 at 35K, consistent with the experimental results. The presence of a significant barrier, 18 kJ mol-1, for the formation of formamide as a product, means that only the H-abstraction channel producing NH3 + CHO, in which the transfer of an H-atom can occur by quantum mechanical tunnelling through a 23 kJ mol-1 barrier, is open at low temperatures. These results are in contrast with a recent theoretical study which suggested that the reaction could proceed without a barrier and was therefore a viable route to gas-phase formamide formation. The calculated rate coefficients were used in an astrochemical model which demonstrated that this reaction produces only negligible amounts of gas-phase formamide under interstellar and circumstellar conditions. The reaction of NH2 with CH2O is therefore not an important source of formamide at low temperatures in interstellar environments. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: Manuscript, 14 pages, 4 figures. Supporting Information, 8 pages, 2 figures. Accepted for publication in The Astrophysical Journal Letters

Mitigating Biases in Student Performance Prediction via Attention-Based Personalized Federated Learning

Authors: Yun-Wei Chu, Seyyedali Hosseinalipour, Elizabeth Tenorio, Laura Cruz, Kerrie Douglas, Andrew Lan, Christopher Brinton

Abstract: Traditional learning-based approaches to student modeling generalize poorly to underrepresented student groups due to biases in data availability. In this paper, we propose a methodology for predicting student performance from their online learning activities that optimizes inference accuracy over different demographic groups such as race and gender. Building upon recent foundations in federated l… ▽ More Traditional learning-based approaches to student modeling generalize poorly to underrepresented student groups due to biases in data availability. In this paper, we propose a methodology for predicting student performance from their online learning activities that optimizes inference accuracy over different demographic groups such as race and gender. Building upon recent foundations in federated learning, in our approach, personalized models for individual student subgroups are derived from a global model aggregated across all student models via meta-gradient updates that account for subgroup heterogeneity. To learn better representations of student activity, we augment our approach with a self-supervised behavioral pretraining methodology that leverages multiple modalities of student behavior (e.g., visits to lecture videos and participation on forums), and include a neural network attention mechanism in the model aggregation stage. Through experiments on three real-world datasets from online courses, we demonstrate that our approach obtains substantial improvements over existing student modeling baselines in predicting student learning outcomes for all subgroups. Visual analysis of the resulting student embeddings confirm that our personalization methodology indeed identifies different activity patterns within different subgroups, consistent with its stronger inference ability compared with the baselines. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 10 pages, CIKM 2022

Experimental study of the removal of excited state phosphorus atoms by H$_{2}$O and H$_{2}$: implications for the formation of PO in stellar winds

Authors: Kevin M. Douglas, David Gobrecht, John M. C. Plane

Abstract: The reactions of the low-lying metastable states of atomic phosphorus, P($^2$D) and P($^2$P), with H$_{2}$O and H$_{2}$ were studied by the pulsed laser photolysis at 248 nm of PCl$_{3}$ , combined with laser induced fluorescence detection of P($^2$D), P($^2$P) and PO. Rate coefficients between 291 and 740 K were measured, along with a yield for the production of PO from P($^2$D or $^2$P) + H… ▽ More The reactions of the low-lying metastable states of atomic phosphorus, P($^2$D) and P($^2$P), with H$_{2}$O and H$_{2}$ were studied by the pulsed laser photolysis at 248 nm of PCl$_{3}$ , combined with laser induced fluorescence detection of P($^2$D), P($^2$P) and PO. Rate coefficients between 291 and 740 K were measured, along with a yield for the production of PO from P($^2$D or $^2$P) + H$_{2}$O of (35$\pm$15)%. H$_{2}$ reacts with both excited P states relatively efficiently; physical (i.e. collisional) quenching, rather than chemical reaction to produced PH + H, is shown to be the more likely pathway. A comprehensive phosphorus chemistry network is then developed using a combination of electronic structure theory calculations and a Master Equation treatment of reactions taking place over complex potential energy surfaces. The resulting model shows that at the high temperatures within two stellar radii of a MIRA variable AGB star in oxygen-rich conditions, collisional excitation of ground-state P($^4$S) to P($^2$D), followed by reaction with H$_{2}$O, is a significant pathway for producing PO (in addition to the reaction between P($^4$S) and OH). The model also demonstrates that the PN fractional abundance in a steady (non-pulsating) outflow is under-predicted by about 2 orders of magnitude. However, under shocked conditions where sufficient thermal dissociation of N$_2$ occurs at temperatures above 4000 K, the resulting N atoms convert a substantial fraction of PO to PN. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 24 pages, 7 figures, accepted for publication in MNRAS

Leveraging Evolutionary Algorithms for Feasible Hexapod Locomotion Across Uneven Terrain

Authors: Jack Vice, Gita Sukthankar, Pamela K. Douglas

Abstract: Optimizing gait stability for legged robots is a difficult problem. Even on level surfaces, effectively traversing across different textures (e.g., carpet) rests on dynamically tuning parameters in multidimensional space. Inspired by biology, evolutionary algorithms (EA) remain an attractive solution for feasibly implementing robotic locomotion with both energetic economy and rapid parameter conve… ▽ More Optimizing gait stability for legged robots is a difficult problem. Even on level surfaces, effectively traversing across different textures (e.g., carpet) rests on dynamically tuning parameters in multidimensional space. Inspired by biology, evolutionary algorithms (EA) remain an attractive solution for feasibly implementing robotic locomotion with both energetic economy and rapid parameter convergence. Here, we leveraged this class of algorithms to evolve a stable hexapod gait controller capable of traversing uneven terrain and obstacles. Gait parameters were evolved in a rigid body dynamics simulation on an 8 x 3 meter obstacle course comprised of random step field, linear obstacles and inclined surfaces. Using a fitness function that jointly optimized locomotion velocity and stability, we found that multiple successful gait parameter evolutions yielded specialized functionality for each leg. Specific gait parameters were identified as critical to developing a rough terrain gait. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: To be published in FLAIRS 2022

ACM Class: I.2.9

Reconsidering Spatial Priors In EEG Source Estimation: Does White Matter Contribute to EEG Rhythms?

Authors: PK Douglas, DB Douglas

Abstract: Electroencephalogram (EEG) has been a core tool used in functional neuroimaging in humans for nearly a hundred years. Because it is inexpensive, easy to implement, and noninvasive, it also represents an excellent candidate modality for use in the BCI setting. Nonetheless, a complete understanding of how EEG measurements (voltage fluctuations) relate to information processing in the brain remains s… ▽ More Electroencephalogram (EEG) has been a core tool used in functional neuroimaging in humans for nearly a hundred years. Because it is inexpensive, easy to implement, and noninvasive, it also represents an excellent candidate modality for use in the BCI setting. Nonetheless, a complete understanding of how EEG measurements (voltage fluctuations) relate to information processing in the brain remains somewhat elusive. A deeper understanding of the neuroanatomical underpinnings of the EEG signal may help explain inter-individual variability in evoked and induced potentials, which may improve BCI therapies targeted to the individual. According to classic biophysical models, EEG fluctuations are primarily a reflection of locally synchronized neuronal oscillations within the gray matter oriented approximately orthogonal to the scalp. In contrast, global models ignore local signals due to dendritic processing, and suggest that propagation delays due to white matter architecture are responsible for the EEG signal, and are capable of explaining the coherence between numerous rhythms (e.g., alpha) at spatially distinct areas of the scalp. Recently, combined local-global models suggest that the EEG signal may reflect a superposition of local processing along with global contributors including transduction along white matter tracts in the brain. Incorporating both local and global (e.g., white matter) priors into EEG source models may therefore improve source estimates. These models may also help disentangle which aspects of the EEG signal are predicted to colocalize spatially with measurements from functional MRI (fMRI). Here, we explore the possibility that white matter conductivity contributes to EEG measurements via a generative model based on classic axonal transduction models, and discuss its potential implications for source estimation. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Journal ref: 2019 IEEE 7th International Winter Conference on Brain-Computer Interface (BCI), 2019, pp. 1-12,

Click-Based Student Performance Prediction: A Clustering Guided Meta-Learning Approach

Authors: Yun-Wei Chu, Elizabeth Tenorio, Laura Cruz, Kerrie Douglas, Andrew S. Lan, Christopher G. Brinton

Abstract: We study the problem of predicting student knowledge acquisition in online courses from clickstream behavior. Motivated by the proliferation of eLearning lecture delivery, we specifically focus on student in-video activity in lectures videos, which consist of content and in-video quizzes. Our methodology for predicting in-video quiz performance is based on three key ideas we develop. First, we mod… ▽ More We study the problem of predicting student knowledge acquisition in online courses from clickstream behavior. Motivated by the proliferation of eLearning lecture delivery, we specifically focus on student in-video activity in lectures videos, which consist of content and in-video quizzes. Our methodology for predicting in-video quiz performance is based on three key ideas we develop. First, we model students' clicking behavior via time-series learning architectures operating on raw event data, rather than defining hand-crafted features as in existing approaches that may lose important information embedded within the click sequences. Second, we develop a self-supervised clickstream pre-training to learn informative representations of clickstream events that can initialize the prediction model effectively. Third, we propose a clustering guided meta-learning-based training that optimizes the prediction model to exploit clusters of frequent patterns in student clickstream sequences. Through experiments on three real-world datasets, we demonstrate that our method obtains substantial improvements over two baseline models in predicting students' in-video quiz performance. Further, we validate the importance of the pre-training and meta-learning components of our framework through ablation studies. Finally, we show how our methodology reveals insights on video-watching behavior associated with knowledge acquisition for useful learning analytics. △ Less

Submitted 15 November, 2021; v1 submitted 28 October, 2021; originally announced November 2021.

Comments: 10 pages, IEEE BigData 2021

Diurnal variations of resting-state fMRI data: A graph-based analysis

Authors: Farzad V. Farahani, Waldemar Karwowski, Mark D Esposito, Richard F. Betzel, Magdalena Fafrowicz, Bartosz Bohaterewicz, Tadeusz Marek, Pamela K. Douglas

Abstract: Circadian rhythms synchronize a variety of physiological processes ranging from neural activity and hormone secretion to sleep cycles and feeding habits. Despite significant diurnal variation, time-of-day (TOD) is rarely recorded or analyzed in human brain research. Moreover, sleep-wake patterns, diurnal preferences, and daytime alertness vary across individuals, known as sleep chronotypes. Here,… ▽ More Circadian rhythms synchronize a variety of physiological processes ranging from neural activity and hormone secretion to sleep cycles and feeding habits. Despite significant diurnal variation, time-of-day (TOD) is rarely recorded or analyzed in human brain research. Moreover, sleep-wake patterns, diurnal preferences, and daytime alertness vary across individuals, known as sleep chronotypes. Here, we performed graph-theory network analysis on resting-state functional MRI (rs-fMRI) data to explore topological differences in whole-brain functional networks between morning and evening sessions (TOD effect), and between extreme morning-type and evening-type chronotypes. To that end, 62 individuals (31 extreme morning, 31 evening-type) underwent two fMRI sessions: about 1 hour after the wake-up time (morning), and 10 hours thereafter, scheduled in accord with their declared habitual sleep-wake pattern. TOD significantly altered functional connectivity (FC) patterns, but there was no significant difference in chronotypic categories. Compared to the morning session, we found relatively increased small-worldness, modularity, assortativity, and synchronization in the evening session, indicating more efficient functional topology. Local measures were changed during the day predominantly across the areas involved in somatomotor, ventral attention, as well as default mode networks. Also, connectivity and hub analyses showed that the somatomotor, ventral attention, and visual networks are the most densely-connected brain areas in both sessions, respectively, with the first being more active in the evening session and the two latter in the morning session. Collectively, these findings suggest TOD can impact classic analyses used in human neuroimaging such as functional connectivity, and should be recorded and included in the analysis of functional neuroimaging data. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 39 pages

CCAT-prime Collaboration: Science Goals and Forecasts with Prime-Cam on the Fred Young Submillimeter Telescope

Authors: CCAT-Prime collaboration, M. Aravena, J. E. Austermann, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, F. Bigiel, J. R. Bond, P. C. Breysse, C. Broughton, R. Bustos, S. C. Chapman, M. Charmetant, S. K. Choi, D. T. Chung, S. E. Clark, N. F. Cothard, A. T. Crites, A. Dev, K. Douglas, C. J. Duell, R. Dunner, H. Ebina, J. Erler , et al. (62 additional authors not shown)

Abstract: We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Corn… ▽ More We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Cornell University and sited at more than 5600 meters on Cerro Chajnantor in northern Chile. Prime-Cam is one of two instruments planned for FYST and will provide unprecedented spectroscopic and broadband measurement capabilities to address important astrophysical questions ranging from Big Bang cosmology through reionization and the formation of the first galaxies to star formation within our own Milky Way galaxy. Prime-Cam on the FYST will have a mapping speed that is over ten times greater than existing and near-term facilities for high-redshift science and broadband polarimetric imaging at frequencies above 300 GHz. We describe details of the science program enabled by this system and our preliminary survey strategies. △ Less

Submitted 8 August, 2022; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: 61 pages, 16 figures. Resubmitted to ApJSS July 11, 2022

The Global Magneto-Ionic Medium Survey: A Faraday Depth Survey of the Northern Sky Covering 1280-1750 MHz

Authors: M. Wolleben, T. L. Landecker, K. A. Douglas, A. D. Gray, A. Ordog, J. M. Dickey, A. S. Hill, E. Carretti, J. C. Brown, B. M. Gaensler, J. L. Han, M. Haverkorn, R. Kothes, J. P. Leahy, N. McClure-Griffiths, D. McConnell, W. Reich, A. R. Taylor, A. J. M. Thomson, J. L. West

Abstract: The Galactic interstellar medium hosts a significant magnetic field, which can be probed through the synchrotron emission produced from its interaction with relativistic electrons. Linearly polarized synchrotron emission is generated throughout the Galaxy, and at longer wavelengths, modified along nearly every path by Faraday rotation in the intervening magneto-ionic medium. Full characterization… ▽ More The Galactic interstellar medium hosts a significant magnetic field, which can be probed through the synchrotron emission produced from its interaction with relativistic electrons. Linearly polarized synchrotron emission is generated throughout the Galaxy, and at longer wavelengths, modified along nearly every path by Faraday rotation in the intervening magneto-ionic medium. Full characterization of the polarized emission requires wideband observations with many frequency channels. We have surveyed polarized radio emission from the Northern sky over the the range 1280-1750 MHz, with channel width 236.8 kHz, using the John A. Galt Telescope (diameter 25.6 m) at the Dominion Radio Astrophysical Observatory, as part of the Global Magneto-Ionic Medium Survey. The survey covered 72% of the sky, declinations -30 to +87 degrees at all right ascensions. The intensity scale was absolutely calibrated, based on the flux density and spectral index of Cygnus A. Polarization angle was calibrated using the extended polarized emission of the Fan Region. Data are presented as brightness temperatures with angular resolution 40'. Sensitivity in Stokes Q and U is 45 mK rms in a 1.18 MHz band. We have applied rotation measure synthesis to the data to obtain a Faraday depth cube of resolution 150 radians per square metre and sensitivity 3 mK rms of polarized intensity. Features in Faraday depth up to a width of 110 radians per square metre are represented. The maximum detectable Faraday depth is +/- 20,000 radians per square metre. The survey data are available at the Canadian Astronomy Data Centre. △ Less

Submitted 24 July, 2021; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the Astronomical Journal

Structural Brain Asymmetries in Youths with Combined and Inattentive Presentations of Attention Deficit Hyperactivity Disorder

Authors: Cintya Nirvana Dutta, Pamela K. Douglas, Hernando Ombao

Abstract: Alterations in structural brain laterality are reported in attention-deficit/hyperactivity disorder (ADHD). However, few studies examined differences within presentations of ADHD. We investigate asymmetry index (AI) across 13 subcortical and 33 cortical regions from anatomical metrics of volume, surface area, and thickness. Structural T1-weighted MRI data were obtained from youths with inattentive… ▽ More Alterations in structural brain laterality are reported in attention-deficit/hyperactivity disorder (ADHD). However, few studies examined differences within presentations of ADHD. We investigate asymmetry index (AI) across 13 subcortical and 33 cortical regions from anatomical metrics of volume, surface area, and thickness. Structural T1-weighted MRI data were obtained from youths with inattentive (n = 64) and combined (n = 51) presentations, and aged-matched controls (n = 298). We used a linear mixed effect model that accounts for data site heterogeneity, while studying associations between AI and covariates of presentation and age. Our paper contributes to the functional results seen among ADHD presentations evidencing disrupted connectivity in motor networks from ADHD-C and cingulo-frontal networks from ADHD-I, as well as new findings in the temporal cortex and default mode subnetworks. Age patterns of structural asymmetries vary with presentation type. Linear mixed effects model is a practical tool for characterizing associations between brain asymmetries, diagnosis, and neurodevelopment. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: 5 pages, 3 figures, 1 table, submitted to ISBI conference

Convalescent Blood Treatment for COVID-19: Are Local Donors Enough?

Authors: Pamela K. Douglas, Farzad V. Farahani, David B. Douglas, Susan Bookheimer

Abstract: COVID-19 is now a global pandemic, and an effective vaccine may be many months away. Over 100 years ago, Spanish flu fatalities were attenuated when doctors began treating patients with blood plasma donated by recovered (or convalesced) survivors. Passive immunity transfer via administration of convalesced blood product (CBP) appears to represent a readily available and promising avenue for mitiga… ▽ More COVID-19 is now a global pandemic, and an effective vaccine may be many months away. Over 100 years ago, Spanish flu fatalities were attenuated when doctors began treating patients with blood plasma donated by recovered (or convalesced) survivors. Passive immunity transfer via administration of convalesced blood product (CBP) appears to represent a readily available and promising avenue for mitigating mortalities, expediting recovery time, and even prophylaxis against the SARS-CoV-2 virus. Here, we review challenges to CBP efficacy, and present a graph theoretical model of transmission dynamics that identifies evolving hubs of COVID-19 cases. Importantly, this model suggests that CBP efficacy may rest on an efficient and distributed global sampling scheme as opposed to CBP pooled from local donors alone. △ Less

Submitted 14 March, 2022; v1 submitted 27 September, 2020; originally announced September 2020.

Comments: This content is no longer relevant

On the Similarity of Deep Learning Representations Across Didactic and Adversarial Examples

Authors: Pamela K. Douglas, Farzad Vasheghani Farahani

Abstract: The increasing use of deep neural networks (DNNs) has motivated a parallel endeavor: the design of adversaries that profit from successful misclassifications. However, not all adversarial examples are crafted for malicious purposes. For example, real world systems often contain physical, temporal, and sampling variability across instrumentation. Adversarial examples in the wild may inadvertently p… ▽ More The increasing use of deep neural networks (DNNs) has motivated a parallel endeavor: the design of adversaries that profit from successful misclassifications. However, not all adversarial examples are crafted for malicious purposes. For example, real world systems often contain physical, temporal, and sampling variability across instrumentation. Adversarial examples in the wild may inadvertently prove deleterious for accurate predictive modeling. Conversely, naturally occurring covariance of image features may serve didactic purposes. Here, we studied the stability of deep learning representations for neuroimaging classification across didactic and adversarial conditions characteristic of MRI acquisition variability. We show that representational similarity and performance vary according to the frequency of adversarial examples in the input space. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 2 figures

ACM Class: I.2.4

Journal ref: Med NeurIPS 2019

Thermodynamic Computing

Authors: Tom Conte, Erik DeBenedictis, Natesh Ganesh, Todd Hylton, John Paul Strachan, R. Stanley Williams, Alexander Alemi, Lee Altenberg, Gavin Crooks, James Crutchfield, Lidia del Rio, Josh Deutsch, Michael DeWeese, Khari Douglas, Massimiliano Esposito, Michael Frank, Robert Fry, Peter Harsha, Mark Hill, Christopher Kello, Jeff Krichmar, Suhas Kumar, Shih-Chii Liu, Seth Lloyd, Matteo Marsili , et al. (14 additional authors not shown)

Abstract: The hardware and software foundations laid in the first half of the 20th Century enabled the computing technologies that have transformed the world, but these foundations are now under siege. The current computing paradigm, which is the foundation of much of the current standards of living that we now enjoy, faces fundamental limitations that are evident from several perspectives. In terms of hard… ▽ More The hardware and software foundations laid in the first half of the 20th Century enabled the computing technologies that have transformed the world, but these foundations are now under siege. The current computing paradigm, which is the foundation of much of the current standards of living that we now enjoy, faces fundamental limitations that are evident from several perspectives. In terms of hardware, devices have become so small that we are struggling to eliminate the effects of thermodynamic fluctuations, which are unavoidable at the nanometer scale. In terms of software, our ability to imagine and program effective computational abstractions and implementations are clearly challenged in complex domains. In terms of systems, currently five percent of the power generated in the US is used to run computing systems - this astonishing figure is neither ecologically sustainable nor economically scalable. Economically, the cost of building next-generation semiconductor fabrication plants has soared past $10 billion. All of these difficulties - device scaling, software complexity, adaptability, energy consumption, and fabrication economics - indicate that the current computing paradigm has matured and that continued improvements along this path will be limited. If technological progress is to continue and corresponding social and economic benefits are to continue to accrue, computing must become much more capable, energy efficient, and affordable. We propose that progress in computing can continue under a united, physically grounded, computational paradigm centered on thermodynamics. Herein we propose a research agenda to extend these thermodynamic foundations into complex, non-equilibrium, self-organizing systems and apply them holistically to future computing systems that will harness nature's innate computational capacity. We call this type of computing "Thermodynamic Computing" or TC. △ Less

Submitted 14 November, 2019; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: A Computing Community Consortium (CCC) workshop report, 36 pages

Report number: ccc2019report_6

On the black hole content and initial mass function of 47 Tuc

Authors: Vincent Hénault-Brunet, Mark Gieles, Jay Strader, Miklos Peuten, Eduardo Balbinot, Kaela E. K. Douglas

Abstract: The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfu… ▽ More The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfully matches the observables and correctly predicts the radial distribution of millisecond pulsars and their gravitational accelerations inferred from long-term timing observations. The data favours a population of BHs with a total mass of $430^{+386}_{-301}$ $M_{\odot}$, but the most likely model has very few BHs. Since our models do not include a central IMBH and accurately reproduce the observations, we conclude that there is currently no need to invoke the presence of an IMBH in 47 Tuc. The global present-day mass function inferred is significantly depleted in low-mass stars (power-law slope $α=-0.52^{+0.17}_{-0.16}$). Given the orbit and predicted mass-loss history of this massive GC, the dearth of low-mass stars is difficult to explain with a standard initial mass function (IMF) followed by long-term preferential escape of low-mass stars driven by two-body relaxation, and instead suggests that 47 Tuc may have formed with a bottom-light IMF. We discuss alternative evolutionary origins for the flat mass function and ways to reconcile this with the low BH retention fraction. Finally, by capturing the effect of dark remnants, our method offers a new way to probe the IMF in a GC above the current main-sequence turn-off mass, for which we find a slope of $-2.49\pm0.08$. △ Less

Submitted 25 October, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 16 pages, 7 figures, accepted to MNRAS after minor revision

The Galactic Magneto-Ionic Medium Survey: Moments of the Faraday Spectra

Authors: John M. Dickey, T. L. Landecker, A. J. M. Thomson, M. Wolleben, X. Sun, E. Carretti, K. Douglas, A. Fletcher, B. M. Gaensler, A. Gray, M. Haverkorn, A. S. Hill, S. A. Mao, N. M. McClure-Griffiths

Abstract: Faraday rotation occurs along every line of sight in the Galaxy; Rotation Measure (RM) synthesis allows a three-dimensional representation of the interstellar magnetic field. This study uses data from the Global Magneto-Ionic Medium Survey, a combination of single-antenna spectro-polarimetric studies, including northern sky data from the DRAO 26-m Telescope (1270-1750 MHz) and southern sky data fr… ▽ More Faraday rotation occurs along every line of sight in the Galaxy; Rotation Measure (RM) synthesis allows a three-dimensional representation of the interstellar magnetic field. This study uses data from the Global Magneto-Ionic Medium Survey, a combination of single-antenna spectro-polarimetric studies, including northern sky data from the DRAO 26-m Telescope (1270-1750 MHz) and southern sky data from the Parkes 64-m Telescope (300-480 MHz). From the synthesized Faraday spectral cubes we compute the zeroth, first, and second moments to find the total polarized emission, mean and RM-width of the polarized emission. From DRAO first moments we find a weak vertical field directed from Galactic North to South, but Parkes data reveal fields directed towards the Sun at high latitudes in both hemispheres: the two surveys clearly sample different volumes. DRAO second moments show feature widths in Faraday spectra increasing with decreasing positive latitudes, implying that longer lines of sight encounter more Faraday rotating medium, but this is not seen at negative latitudes. Parkes data show the opposite: at positive latitudes the second moment decreases with decreasing latitude, but not at negative latitudes. Comparing first moments with RMs of pulsars and extragalactic sources and a study of depolarization together confirm that the DRAO survey samples to larger distances than the Parkes data. Emission regions in the DRAO survey are typically 700 to 1000 pc away, slightly beyond the scale-height of the magneto-ionic medium; emission detected in the Parkes survey is entirely within the magneto-ionic disk, less than 500 pc away. △ Less

Submitted 13 December, 2018; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: Astrophysical Journal in press

Interpreting Encoding and Decoding Models

Authors: Nikolaus Kriegeskorte, Pamela K. Douglas

Abstract: Encoding and decoding models are widely used in systems, cognitive, and computational neuroscience to make sense of brain-activity data. However, the interpretation of their results requires care. Decoding models can help reveal whether particular information is present in a brain region in a format the decoder can exploit. Encoding models make comprehensive predictions about representational spac… ▽ More Encoding and decoding models are widely used in systems, cognitive, and computational neuroscience to make sense of brain-activity data. However, the interpretation of their results requires care. Decoding models can help reveal whether particular information is present in a brain region in a format the decoder can exploit. Encoding models make comprehensive predictions about representational spaces. In the context of sensory systems, encoding models enable us to test and compare brain-computational models, and thus directly constrain computational theory. Encoding and decoding models typically include fitted linear-model components. Sometimes the weights of the fitted linear combinations are interpreted as reflecting, in an encoding model, the contribution of different sensory features to the representation or, in a decoding model, the contribution of different measured brain responses to a decoded feature. Such interpretations can be problematic when the predictor variables or their noise components are correlated and when priors (or penalties) are used to regularize the fit. Encoding and decoding models are evaluated in terms of their generalization performance. The correct interpretation depends on the level of generalization a model achieves (e.g. to new response measurements for the same stimuli, to new stimuli from the same population, or to stimuli from a different population). Significant decoding or encoding performance of a single model (at whatever level of generality) does not provide strong constraints for theory. Many models must be tested and inferentially compared for analyses to drive theoretical progress. △ Less

Submitted 26 April, 2019; v1 submitted 1 December, 2018; originally announced December 2018.

Comments: 19 pages, 2 figures, author preprint

Cognitive computational neuroscience

Authors: Nikolaus Kriegeskorte, Pamela K. Douglas

Abstract: To learn how cognition is implemented in the brain, we must build computational models that can perform cognitive tasks, and test such models with brain and behavioral experiments. Cognitive science has developed computational models of human cognition, decomposing task performance into computational components. However, its algorithms still fall short of human intelligence and are not grounded in… ▽ More To learn how cognition is implemented in the brain, we must build computational models that can perform cognitive tasks, and test such models with brain and behavioral experiments. Cognitive science has developed computational models of human cognition, decomposing task performance into computational components. However, its algorithms still fall short of human intelligence and are not grounded in neurobiology. Computational neuroscience has investigated how interacting neurons can implement component functions of brain computation. However, it has yet to explain how those components interact to explain human cognition and behavior. Modern technologies enable us to measure and manipulate brain activity in unprecedentedly rich ways in animals and humans. However, experiments will yield theoretical insight only when employed to test brain-computational models. It is time to assemble the pieces of the puzzle of brain computation. Here we review recent work in the intersection of cognitive science, computational neuroscience, and artificial intelligence. Computational models that mimic brain information processing during perceptual, cognitive, and control tasks are beginning to be developed and tested with brain and behavioral data. △ Less

Submitted 31 July, 2018; originally announced July 2018.

Comments: 31 pages, 4 figures

Hemispheric Preference and Cyclic Variation of Solar Filament Chirality from 2000 to 2016

Authors: Soumitra Hazra, Sushant S. Mahajan, William Keith Douglas Jr., Petrus C. H. Martens

Abstract: It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challengi… ▽ More It is well known that solar filaments are features in the solar atmosphere which show a hemispheric preference in their chirality. The hemispheric preference is such that the dextral chirality dominates in the northern hemisphere while the sinistral chirality dominates in the southern. Determining the strength and cyclic variation of the degree of this hemispheric preference, however, is challenging and tedious and thus needs to be automated. In this paper, we follow Dr. Pietro Bernasconi's algorithm (Bernasconi et al. 2005) to detect filament chirality with two parallel channels of application. The algorithm is applied to H-alpha images by "Advanced Automated Filament Detection and Characterization Code" (AAFDCC) (Bernasconi et al. 2005) or the algorithm is explained to a human and the human determines the chirality of the solar filament. We have conducted this exercise on data during the month of August from years 2000 to 2016 and we found that 83 % of our visually determined filaments follow the hemispheric chirality preference, while 58 % of automatically determined filaments follow it. Our visually compiled results have over 90 % agreement with those of Pevtsov et al. (2003), yet the visually determined chiralities of filaments disagree with automated determinations significantly. We find that the hemispheric preference remained the same between solar cycles 23 and 24 but the preference is very difficult to determine during the solar minimum of 2008-2010 primarily due to the absence of filaments. △ Less

Submitted 5 October, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: Accepted in Astrophysical Journal

Journal ref: 2018, ApJ, 865, 108

The Fan Region at 1.5 GHz. I: Polarized synchrotron emission extending beyond the Perseus Arm

Authors: A. S. Hill, T. L. Landecker, E. Carretti, K. Douglas, X. H. Sun, B. M. Gaensler, S. A. Mao, N. M. McClure-Griffiths, W. Reich, M. Wolleben, J. M. Dickey, A. D. Gray, M. Haverkorn, J. P. Leahy, D. H. F. M. Schnitzeler

Abstract: The Fan Region is one of the dominant features in the polarized radio sky, long thought to be a local (distance < 500 pc) synchrotron feature. We present 1.3-1.8 GHz polarized radio continuum observations of the region from the Global Magneto-Ionic Medium Survey (GMIMS) and compare them to maps of Halpha and polarized radio continuum intensity from 0.408-353 GHz. The high-frequency (> 1 GHz) and l… ▽ More The Fan Region is one of the dominant features in the polarized radio sky, long thought to be a local (distance < 500 pc) synchrotron feature. We present 1.3-1.8 GHz polarized radio continuum observations of the region from the Global Magneto-Ionic Medium Survey (GMIMS) and compare them to maps of Halpha and polarized radio continuum intensity from 0.408-353 GHz. The high-frequency (> 1 GHz) and low-frequency (< 600 MHz) emission have different morphologies, suggesting a different physical origin. Portions of the 1.5 GHz Fan Region emission are depolarized by about 30% by ionized gas structures in the Perseus Arm, indicating that this fraction of the emission originates >2 kpc away. We argue for the same conclusion based on the high polarization fraction at 1.5 GHz (about 40%). The Fan Region is offset with respect to the Galactic plane, covering -5° < b < +10°; we attribute this offset to the warp in the outer Galaxy. We discuss origins of the polarized emission, including the spiral Galactic magnetic field. This idea is a plausible contributing factor although no model to date readily reproduces all of the observations. We conclude that models of the Galactic magnetic field should account for the > 1 GHz emission from the Fan Region as a Galactic-scale, not purely local, feature. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: Accepted to MNRAS

Journal ref: MNRAS 467 4631 (2017)

A High-Velocity Cloud Impact Forming a Supershell in the Milky Way

Authors: Geumsook Park, Bon-Chul Koo, Ji-hyun Kang, Steven J. Gibson, J. E. G. Peek, Kevin A. Douglas, Eric J. Korpela, Carl E. Heiles

Abstract: Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HV… ▽ More Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HVCs) falling to the Galactic disk. Here we report the detection of a kiloparsec (kpc)-size supershell in the outskirts of the Milky Way with the compact HVC 040+01$-$282 (hereafter CHVC040) at its geometrical center using the "Inner-Galaxy Arecibo L-band Feed Array" HI 21-cm survey data. The morphological and physical properties of both objects suggest that CHVC040, which is either a fragment of a nearby disrupted galaxy or a cloud originated from an intergalactic accreting flow, collided with the disk $\sim 5$ Myrs ago to form the supershell. Our result shows that some compact HVCs can survive their trip through the Galactic halo and inject energy and momentum into the Milky Way disk. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 10 pages, 5 figures, Accepted for publication in the Astrophysical Journal Letters on July 26, 2016

Gain and Polarization Properties of a Large Radio Telescope from Calculation and Measurement: The John A. Galt Telescope

Authors: Xuan Du, Thomas L. Landecker, Timothy Robishaw, Andrew D. Gray, Kevin A. Douglas, Maik Wolleben

Abstract: Measurement of the brightness temperature of extended radio emission demands knowledge of the gain (or aperture efficiency) of the telescope and measurement of the polarized component of the emission requires correction for the conversion of unpolarized emission from sky and ground to apparently polarized signal. Radiation properties of the John A. Galt Telescope at the Dominion Radio Astrophysica… ▽ More Measurement of the brightness temperature of extended radio emission demands knowledge of the gain (or aperture efficiency) of the telescope and measurement of the polarized component of the emission requires correction for the conversion of unpolarized emission from sky and ground to apparently polarized signal. Radiation properties of the John A. Galt Telescope at the Dominion Radio Astrophysical Observatory were studied through analysis and measurement in order to provide absolute calibration of a survey of polarized emission from the entire northern sky from 1280 to 1750 MHz, and to understand the polarization performance of the telescope. Electromagnetic simulators CST and GRASP-10 were used to compute radiation patterns of the telescope in all Stokes parameters, and aperture efficiency. Aperture efficiency was also evaluated using geometrical optics and was measured using Cyg A. Measured aperture efficiency varied smoothly with frequency between values of 0.49 and 0.54; GRASP-10 yielded values 6.5% higher but with closely similar variation with frequency. Overall error across the frequency band is 3%, but values at any two frequencies are relatively correct to ~1%. Dominant influences on aperture efficiency are illumination taper of the feed radiation pattern and shadowing by the feed-support struts. A model of ground emission was developed based on measurements and on empirical data from remote sensing of the Earth from satellite-borne telescopes. This model was convolved with the computed antenna response to estimate conversion of ground emission into spurious polarized signal. The computed spurious signal is comparable to measured values, but is not accurate enough to be used to correct observations. A simpler model, in which the ground is considered as an unpolarized emitter with a brightness temperature of ~240 K, is shown to have useful accuracy when compared to measurements. △ Less

Submitted 21 July, 2016; originally announced July 2016.

Comments: 14 pages, 12 figures

Decoding the Encoding of Functional Brain Networks: an fMRI Classification Comparison of Non-negative Matrix Factorization (NMF), Independent Component Analysis (ICA), and Sparse Coding Algorithms

Authors: Jianwen Xie, Pamela K. Douglas, Ying Nian Wu, Arthur L. Brody, Ariana E. Anderson

Abstract: Brain networks in fMRI are typically identified using spatial independent component analysis (ICA), yet mathematical constraints such as sparse coding and positivity both provide alternate biologically-plausible frameworks for generating brain networks. Non-negative Matrix Factorization (NMF) would suppress negative BOLD signal by enforcing positivity. Spatial sparse coding algorithms ($L1$ Regula… ▽ More Brain networks in fMRI are typically identified using spatial independent component analysis (ICA), yet mathematical constraints such as sparse coding and positivity both provide alternate biologically-plausible frameworks for generating brain networks. Non-negative Matrix Factorization (NMF) would suppress negative BOLD signal by enforcing positivity. Spatial sparse coding algorithms ($L1$ Regularized Learning and K-SVD) would impose local specialization and a discouragement of multitasking, where the total observed activity in a single voxel originates from a restricted number of possible brain networks. The assumptions of independence, positivity, and sparsity to encode task-related brain networks are compared; the resulting brain networks for different constraints are used as basis functions to encode the observed functional activity at a given time point. These encodings are decoded using machine learning to compare both the algorithms and their assumptions, using the time series weights to predict whether a subject is viewing a video, listening to an audio cue, or at rest, in 304 fMRI scans from 51 subjects. For classifying cognitive activity, the sparse coding algorithm of $L1$ Regularized Learning consistently outperformed 4 variations of ICA across different numbers of networks and noise levels (p$<$0.001). The NMF algorithms, which suppressed negative BOLD signal, had the poorest accuracy. Within each algorithm, encodings using sparser spatial networks (containing more zero-valued voxels) had higher classification accuracy (p$<$0.001). The success of sparse coding algorithms may suggest that algorithms which enforce sparse coding, discourage multitasking, and promote local specialization may capture better the underlying source processes than those which allow inexhaustible local processes such as ICA. △ Less

Submitted 1 July, 2016; originally announced July 2016.

The Future of Computing Research: Industry-Academic Collaborations

Authors: Nady Boules, Khari Douglas, Stuart Feldman, Limor Fix, Gregory Hager, Brent Hailpern, Martial Hebert, Dan Lopresti, Beth Mynatt, Chris Rossbach, Helen Wright

Abstract: IT-driven innovation is an enormous factor in the worldwide economic leadership of the United States. It is larger than finance, construction, or transportation, and it employs nearly 6% of the US workforce. The top three companies, as measured by market capitalization, are IT companies - Apple, Google (now Alphabet), and Microsoft. Facebook, a relatively recent entry in the top 10 list by market… ▽ More IT-driven innovation is an enormous factor in the worldwide economic leadership of the United States. It is larger than finance, construction, or transportation, and it employs nearly 6% of the US workforce. The top three companies, as measured by market capitalization, are IT companies - Apple, Google (now Alphabet), and Microsoft. Facebook, a relatively recent entry in the top 10 list by market capitalization has surpassed Walmart, the nation's largest retailer, and the largest employer in the world. The net income of just the top three exceeds $80 billion - roughly 100 times the total budget of the NSF CISE directorate which funds 87% of computing research. In short, the direct return on federal research investments in IT research has been enormously profitable to the nation. The IT industry ecosystem is also evolving. The time from conception to market of successful products has been cut from years to months. Product life cycles are increasingly a year or less. This change has pressured companies to focus industrial R&D on a pipeline or portfolio of technologies that bring immediate, or almost immediate, value to the companies. To defeat the competition and stay ahead of the pack, a company must devote resources to realizing gains that are shorter term, and must remain agile to respond quickly to market changes driven by new technologies, new startups, evolving user experience expectations, and the continuous consumer demand for new and exciting products. Amidst this landscape, the Computing Community Consortium convened a round-table of industry and academic participants to better understand the landscape of industry-academic interaction, and to discuss possible actions that might be taken to enhance those interactions. We close with some recommendations for actions that could expand the lively conversation we experienced at the round-table to a national scale. △ Less

Submitted 29 June, 2016; originally announced June 2016.

Comments: A Computing Community Consortium (CCC) white paper, 19 pages

Faraday Tomography of the North Polar Spur: Constraints on the distance to the Spur and on the Magnetic Field of the Galaxy

Authors: X. H. Sun, T. L. Landecker, B. M. Gaensler, E. Carretti, W. Reich, J. P. Leahy, N. M. McClure-Griffiths, R. M. Crocker, M. Wolleben, M. Haverkorn, K. A. Douglas, A. D. Gray

Abstract: We present radio continuum and polarization images of the North Polar Spur (NPS) from the Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. Combining this RM map with a p… ▽ More We present radio continuum and polarization images of the North Polar Spur (NPS) from the Global Magneto-Ionic Medium Survey (GMIMS) conducted with the Dominion Radio Astrophysical Observatory 26-m Telescope. We fit polarization angle versus wavelength squared over 2048 frequency channels from 1280 to 1750 MHz to obtain a Faraday Rotation Measure (RM) map of the NPS. Combining this RM map with a published Faraday depth map of the entire Galaxy in this direction, we derive the Faraday depth introduced by the NPS and the Galactic interstellar medium (ISM) in front of and behind the NPS. The Faraday depth contributed by the NPS is close to zero, indicating that the NPS is an emitting only feature. The Faraday depth caused by the ISM in front of the NPS is consistent with zero at b>50 degree, implying that this part of the NPS is local at a distance of approximately several hundred parsecs. The Faraday depth contributed by the ISM behind the NPS gradually increases with Galactic latitude up to b=44 degree, and decreases at higher Galactic latitudes. This implies that either the part of the NPS at b<44 degree is distant or the NPS is local but there is a sign change of the large-scale magnetic field. If the NPS is local, there is then no evidence for a large-scale anti-symmetry pattern in the Faraday depth of the Milky Way. The Faraday depth introduced by the ISM behind the NPS at latitudes b>50 degree can be explained by including a coherent vertical magnetic field. △ Less

Submitted 16 August, 2015; originally announced August 2015.

Comments: 9 pages, 8 figures, accepted for publication in ApJ. Some figures have been degraded to reduce sizes, for a high resolution version, see http://physics.usyd.edu.au/~xhsun/ms_nps.pdf

Synthetic CO, H2 and HI surveys of the Galactic 2nd Quadrant, and the properties of molecular gas

Authors: A. Duarte-Cabral, D. M. Acreman, C. L. Dobbs, J. C. Mottram, S. J. Gibson, C. M. Brunt, K. A. Douglas

Abstract: We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic dat… ▽ More We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic datacubes to observations of the 2nd Quadrant of the Milky Way to test the ability of the current models to reproduce the basic chemistry of the Galactic ISM, as well as to test how sensitive such galaxy models are to different recipes of chemistry and/or feedback. We find that models which include feedback and self-gravity can reproduce the production of CO with respect to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge impact on the statistical properties of the chemistry in the simulated galaxies, we find that the inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback failed to reproduce all of the observables. Finally, even though the transition from H2 to CO seems to be robust, we find that all models seem to underproduce molecular gas, and have a lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with feedback and self-gravity has shown to be robust in reproducing the statistical properties of the basic molecular gas components of the ISM in our Galaxy. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: Accepted for publication in MNRAS

Exploring GLIMPSE Bubble N107: Multiwavelength Observations and Simulations

Authors: Vojtech Sidorin, Kevin A. Douglas, Jan Palous, Richard Wunsch, Sona Ehlerova

Abstract: Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of… ▽ More Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of numerical simulations, we aim to estimate the bubble age and other parameters which cannot be derived directly from observations. Methods. From the observations of the HI (I-GALFA) and 13CO (GRS) lines we derive the bubble's kinematical distance and masses of the atomic and molecular components. With the algorithm DENDROFIND, we decompose molecular material into individual clumps. From the continuum observations at 1420 MHz (VGPS) and 327 MHz (WSRT), we derive the radio flux density and the spectral index. With the numerical code ring, we simulate the evolution of stellar-blown bubbles similar to N107. Results. The total HI mass associated with N107 is 5.4E3 Msun. The total mass of the molecular component (a mixture of cold gasses of H2, CO, He and heavier elements) is 1.3E5 Msun, from which 4.0E4 Msun is found along the bubble border. We identified 49 molecular clumps distributed along the bubble border, with the slope of the clump mass function of -1.1. The spectral index of -0.30 of a strong radio source located apparently within the bubble indicates nonthermal emission, hence part of the flux likely originates in a supernova remnant, not yet catalogued. The numerical simulations suggest N107 is likely less than 2.25 Myr old. Since first supernovae explode only after 3 Myr or later, no supernova remnant should be present within the bubble. It may be explained if there is a supernova remnant in the direction towards the bubble, however not associated with it. △ Less

Submitted 23 February, 2014; originally announced February 2014.

Comments: 15 pages, 11 figures

Characterizing the Turbulent Properties of the Starless Molecular Cloud MBM16

Authors: Nickolas M. Pingel, Snezana Stanimirovic, J. E. G. Peek, Min-Young Lee, Alex Lazarian, Blakesley Burkhart, Ayesha Begum, Kevin A. Douglas, Carl Heiles, Steven J. Gibson, Jana Grcevich, Eric J. Korpela, Allen Lawrence, Claire Murray, Mary E. Putman, Destry Saul

Abstract: We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a si… ▽ More We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from -3.3 to -3.7 and being consistent over the velocity range of MBM16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to HI intensity fluctuations. By using this variation we estimate the slope of 3D density fluctuations in MBM16 to be -3.7\pm0.2. This is significantly steeper than what has been found for HI in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum. △ Less

Submitted 27 October, 2013; originally announced October 2013.

Comments: 12 Pages; 10 figures

Hi Shells and Supershells in the I-GALFA Hi 21-cm Line Survey: I. Fast-Expanding Hi Shells Associated with Supernova Remnants

Authors: Geumsook Park, Bon-Chul Koo, Steven J. Gibson, Ji-hyun Kang, Daria C. Lane, Kevin A. Douglas, Joshua E. G. Peek, Eric J. Korpela, Carl E. Heiles, Jonathan H. Newton

Abstract: We search for fast-expanding Hi shells associated with Galactic supernova remnants (SNRs) in the longitude range l \approx 32\arcdeg to 77\arcdeg using 21-cm line data from the Inner-Galaxy Arecibo L-band Feed Array (I-GALFA) Hi survey. Among the 39 known Galactic SNRs in this region, we find such Hi shells in four SNRs: W44, G54.4-0.3, W51C, and CTB 80. All four were previously identified in low-… ▽ More We search for fast-expanding Hi shells associated with Galactic supernova remnants (SNRs) in the longitude range l \approx 32\arcdeg to 77\arcdeg using 21-cm line data from the Inner-Galaxy Arecibo L-band Feed Array (I-GALFA) Hi survey. Among the 39 known Galactic SNRs in this region, we find such Hi shells in four SNRs: W44, G54.4-0.3, W51C, and CTB 80. All four were previously identified in low-resolution surveys, and three of those (excluding G54.4-0.3) were previously studied with the Arecibo telescope. A remarkable new result, however, is the detection of Hi emission at both very high positive and negative velocities in W44 from the receding and approaching parts of the Hi expanding shell, respectively. This is the first detection of both sides of an expanding shell associated with an SNR in Hi 21-cm emission. The high-resolution I-GALFA survey data also reveal a prominent expanding Hi shell with high circular symmetry associated with G54.4-0.3. We explore the physical characteristics of four SNRs and discuss what differentiates them from other SNRs in the survey area. We conclude that these four SNRs are likely the remnants of core-collapse supernovae interacting with a relatively dense (>~1 cm^-3) ambient medium, and we discuss the visibility of SNRs in the Hi 21-cm line. △ Less

Submitted 27 June, 2013; originally announced June 2013.

Comments: 36 pages, 11 figures, 3 tables; accepted for publication in ApJ

The GALFA-HI Compact Cloud Catalog

Authors: Destry R. Saul, J. E. G. Peek, J. Grcevich, M. E. Putman, K. A. Douglas, E. J. Korpela, S. Stanimirovic, C. Heiles, S. J. Gibson, M. Lee, A. Begum, A. R. H. Brown, B. Burkhart, E. T. Hamden, N. M. Pingel, S. Tonnesen

Abstract: We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and col… ▽ More We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating. △ Less

Submitted 21 May, 2015; v1 submitted 20 August, 2012; originally announced August 2012.

Comments: 34 Pages, 9 Figures, published in ApJ (2012, ApJ, 758, 44), this version has the corrected fluxes and corresponding flux histogram and masses

GASKAP -- The Galactic ASKAP Survey

Authors: J. M. Dickey, N. McClure-Griffiths, S. J. Gibson, J. F. Gomez, H. Imai, P. Jones, S. Stanimirovic, J. Th. van Loon, A. Walsh, A. Alberdi, G. Anglada, L. Uscanga, H. Arce, M. Bailey, A. Begum, B. Wakker, N. Ben Bekhti, P. Kalberla, B. Winkel, K. Bekki, B. -Q. For, L. Staveley-Smith, T. Westmeier, M. Burton, M. Cunningham , et al. (33 additional authors not shown)

Abstract: A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emis… ▽ More A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream. △ Less

Submitted 3 July, 2012; originally announced July 2012.

Comments: 45 pages, 8 figures, Pub. Astron. Soc. Australia (in press)

Hydroxyl as a Tracer of H2 in the Envelope of MBM40

Authors: David L. Cotten, Loris Magnani, Elizabeth A. Wennerstrom, Kevin A. Douglas, Joseph S. Onello

Abstract: We observed 51 positions in the OH 1667 MHz main line transitions in the translucent, high latitude cloud MBM40. We detected OH emission in 8 out of 8 positions in the molecular core of the cloud and 24 out of 43 in the surrounding, lower extinction envelope and periphery of the cloud. Using a linear relationship between the integrated OH line intensity and E(B-V), we estimate the mass in the core… ▽ More We observed 51 positions in the OH 1667 MHz main line transitions in the translucent, high latitude cloud MBM40. We detected OH emission in 8 out of 8 positions in the molecular core of the cloud and 24 out of 43 in the surrounding, lower extinction envelope and periphery of the cloud. Using a linear relationship between the integrated OH line intensity and E(B-V), we estimate the mass in the core, the envelope, and the periphery of the cloud to be 4, 8, and 5 solar masses. As much as a third of the total cloud mass may be found in the in the periphery (E(B-V) $<$ 0.12 mag) and about a half in the envelope (0.12 $\le$ E(B-V) $\le$ 0.17 mag). If these results are applicable to other translucent clouds the OH 1667 MHz line is an excellent tracer of gas in very low extinction regions and high-sensitivity mapping of the envelopes of molecular clouds may reveal the presence of significant quantities of molecular mass. △ Less

Submitted 24 September, 2012; v1 submitted 3 March, 2012; originally announced March 2012.

Comments: 26 pages, 3 figures, and 5 tables

The structure of HI in galactic disks: Simulations vs observations

Authors: David M. Acreman, Clare L. Dobbs, Christopher M. Brunt, Kevin A. Douglas

Abstract: We generate synthetic HI Galactic plane surveys from spiral galaxy simulations which include stellar feedback processes. Compared to a model without feedback we find an increased scale height of HI emission (in better agreement with observations) and more realistic spatial structure (including supernova blown bubbles). The synthetic data show HI self-absorption with a morphology similar to that se… ▽ More We generate synthetic HI Galactic plane surveys from spiral galaxy simulations which include stellar feedback processes. Compared to a model without feedback we find an increased scale height of HI emission (in better agreement with observations) and more realistic spatial structure (including supernova blown bubbles). The synthetic data show HI self-absorption with a morphology similar to that seen in observations. The density and temperature of the material responsible for HI self-absorption is consistent with observationally determined values, and is found to be only weakly dependent on absorption strength and star formation efficiency. △ Less

Submitted 24 January, 2012; originally announced January 2012.

Comments: 12 pages, 7 figures. Accepted for publication in MNRAS

A High Resolution Study of the HI-H2 Transition across the Perseus Molecular Cloud

Authors: Min-Young Lee, Snezana Stanimirovic, Kevin A. Douglas, Lewis B. G. Knee, James Di Francesco, Steven J. Gibson, Ayesha Begum, Jana Grcevich, Carl Heiles, Eric J. Korpela, Adam K. Leroy, J. E. G. Peek, Nick Pingel, Mary E. Putman, Destry Saul

Abstract: To investigate the fundamental principles of H2 formation in a giant molecular cloud (GMC), we derive the HI and H2 surface density (Sigma_HI and Sigma_H2) images of the Perseus molecular cloud on sub-pc scales (~0.4 pc). We use the far-infrared data from the Improved Reprocessing of the IRAS Survey and the V-band extinction image provided by the COMPLETE Survey to estimate the dust column density… ▽ More To investigate the fundamental principles of H2 formation in a giant molecular cloud (GMC), we derive the HI and H2 surface density (Sigma_HI and Sigma_H2) images of the Perseus molecular cloud on sub-pc scales (~0.4 pc). We use the far-infrared data from the Improved Reprocessing of the IRAS Survey and the V-band extinction image provided by the COMPLETE Survey to estimate the dust column density image of Perseus. In combination with the HI data from the Galactic Arecibo L-band Feed Array HI Survey and an estimate of the local dust-to-gas ratio, we then derive the Sigma_H2 distribution across Perseus. We find a relatively uniform Sigma_HI ~ 6-8 Msun pc^-2 for both dark and star-forming regions, suggesting a minimum HI surface density required to shield H2 against photodissociation. As a result, a remarkably tight and consistent relation is found between Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2. The transition between the HI- and H2-dominated regions occurs at N(HI)+2N(H2) ~ (8-14) x 10^20 cm^-2. Our findings are consistent with predictions for H2 formation in equilibrium, suggesting that turbulence may not be of primary importance for H2 formation. However, the importance of a warm neutral medium for H2 shielding, an internal radiation field, and the timescale of H2 formation still remain as open questions. We also compare H2 and CO distributions and estimate the fraction of "CO-dark" gas, f_DG ~ 0.3. While significant spatial variations of f_DG are found, we do not find a clear correlation with the mean V-band extinction. △ Less

Submitted 11 May, 2012; v1 submitted 12 October, 2011; originally announced October 2011.

Comments: updated to match the final version published in April 2012

Journal ref: Astrophys.J.748:75,2012

The GALFA-HI Survey: Data Release 1

Authors: J. E. G. Peek, Carl Heiles, Kevin A. Douglas, Min-Young Lee, Jana Grcevich, Snezana Stanimirovic, M. E. Putman, Eric J. Korpela, Steven J. Gibson, Ayesha Begum, Destry Saul, Timothy Robishaw, Marko Krco

Abstract: We present the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, and its first full data release (DR1). GALFA-HI is a high resolution (~ 4'), large area (13000 deg^2), high spectral resolution (0.18 km/s), wide band (-700 < v_LSR < +700 km/s) survey of the Galactic interstellar medium in the 21-cm line hyperfine transition of neutral hydrogen conducted at Arecibo Observatory. Typical noise… ▽ More We present the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, and its first full data release (DR1). GALFA-HI is a high resolution (~ 4'), large area (13000 deg^2), high spectral resolution (0.18 km/s), wide band (-700 < v_LSR < +700 km/s) survey of the Galactic interstellar medium in the 21-cm line hyperfine transition of neutral hydrogen conducted at Arecibo Observatory. Typical noise levels are 80 mK RMS in an integrated 1 km/s channel. GALFA-HI is a dramatic step forward in high-resolution, large-area Galactic HI surveys, and we compare GALFA-HI to past, present, and future Galactic HI surveys. We describe in detail new techniques we have developed to reduce these data in the presence of fixed pattern noise, gain variation, and inconsistent beam shapes, and we show how we have largely mitigated these effects. We present our first full data release, covering 7520 square degrees of sky and representing 3046 hours of integration time, and discuss the details of these data. △ Less

Submitted 11 January, 2011; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: Accepted to the ApJS

Galactic Small Scale Structure Revealed by the GALFA-HI Survey

Authors: Ayesha Begum, Snezana Stanimirovic, Joshua E. Peek, Nicholas Ballering, Carl Heiles, Kevin A. Douglas, Mary Putman, Steven Gibson, Jana Grcevich, Eric Korpela, Min-Young Lee, Destry Saul, John S. Gallagher III

Abstract: The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and an angular resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of many isolated, very compact HI clouds at low Galactic velocities… ▽ More The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and an angular resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of many isolated, very compact HI clouds at low Galactic velocities which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ searched so far, we have detected 96 such compact clouds. The detected clouds are cold with kinetic temperature less than 300 K. Moreover, they are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5\times10^{18}$ cm$^{-2}$ in HI column density. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of clouds are related to the Galactic disk, and are within a few kpc distance. We present properties of the compact clouds sample and discuss various possible scenarios for the origin of this clouds population and its role in the Galactic interstellar medium studies. △ Less

Submitted 18 August, 2010; originally announced August 2010.

Comments: 7 pages, 3 figures. To appear in "The Dynamic ISM: A celebration of the Canadian Galactic Plane Survey," ASP Conference Series

Compact HI clouds from the GALFA-HI survey

Authors: Ayesha Begum, Snezana Stanimirovic, Joshua E. Peek, Nicholas P. Ballering, Carl Heiles, Kevin A. Douglas, Mary Putman, Steven J. Gibson, Jana Grcevich, Eric J. Korpela, Min-Young Lee, Destry Saul, John S. Gallagher III

Abstract: The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and a spatial resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of numerous isolated, very compact HI clouds at low Galactic velocit… ▽ More The Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey is mapping the entire Arecibo sky at 21-cm, over a velocity range of -700 to +700 km/s (LSR), at a velocity resolution of 0.18 km/s and a spatial resolution of 3.5 arcmin. The unprecedented resolution and sensitivity of the GALFA-HI survey have resulted in the detection of numerous isolated, very compact HI clouds at low Galactic velocities, which are distinctly separated from the HI disk emission. In the limited area of ~4600 deg$^2$ surveyed so far, we have detected 96 of such compact clouds. The detected clouds are cold with a median T$_{k,max}$ (the kinetic temperature in the case in which there is no non-thermal broadening) of 300 K. Moreover, these clouds are quite compact and faint, with median values of 5 arcmin in angular size, 0.75 K in peak brightness temperature, and $5 \times 10^{18}$ cm$^{-2}$ in HI column density. Most of the clouds deviate from Galactic rotation at the 20-30 km/s level, and a significant fraction show evidence for a multiphase medium and velocity gradients. No counterparts for these clouds were found in other wavebands. From the modeling of spatial and velocity distributions of the whole compact cloud population, we find that the bulk of the compact clouds are related to the Galactic disk, and their distances are likely to be in the range of 0.1 to a few kpc. We discuss various possible scenarios for the formation and maintenance of this cloud population and its significance for Galactic ISM studies. △ Less

Submitted 7 August, 2010; originally announced August 2010.

Comments: Accepted for publication in the Astrophysical Journal

A Synthetic 21-cm Galactic Plane Survey of an SPH Galaxy Simulation

Authors: Kevin A. Douglas, David M. Acreman, Clare L. Dobbs, Christopher M. Brunt

Abstract: We have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90 degrees < l < 180 degrees, using a model spiral galaxy from SPH simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile i… ▽ More We have created synthetic neutral hydrogen (HI) Galactic Plane Survey data cubes covering 90 degrees < l < 180 degrees, using a model spiral galaxy from SPH simulations and the radiative transfer code TORUS. The density, temperature and other physical parameters are fed from the SPH simulation into TORUS, where the HI emissivity and opacity are calculated before the 21-cm line emission profile is determined. Our main focus is the observation of Outer Galaxy `Perseus Arm' HI, with a view to tracing atomic gas as it encounters shock motions as it enters a spiral arm interface, an early step in the formation of molecular clouds. The observation of HI self-absorption features at these shock sites (in both real observations and our synthetic data) allows us to investigate further the connection between cold atomic gas and the onset of molecular cloud formation. △ Less

Submitted 23 April, 2010; originally announced April 2010.

Comments: MNRAS accepted; 11 pages, 12 figures

Journal ref: MNRAS 407 (2010) 405-414

Synthetic HI observations of a simulated spiral galaxy

Authors: David M. Acreman, Kevin A. Douglas, Clare L. Dobbs, Christopher M. Brunt

Abstract: Using the Torus radiative transfer code we produce synthetic observations of the 21 cm neutral hydrogen line from an SPH simulation of a spiral galaxy. The SPH representation of the galaxy is mapped onto an AMR grid, and a ray tracing method is used to calculate 21 cm line emission for lines of sight through the AMR grid in different velocity channels and spatial pixels. The result is a synthetic… ▽ More Using the Torus radiative transfer code we produce synthetic observations of the 21 cm neutral hydrogen line from an SPH simulation of a spiral galaxy. The SPH representation of the galaxy is mapped onto an AMR grid, and a ray tracing method is used to calculate 21 cm line emission for lines of sight through the AMR grid in different velocity channels and spatial pixels. The result is a synthetic spectral cube which can be directly compared to real observations. We compare our synthetic spectral cubes to observations of M31 and M33 and find good agreement, whereby increasing velocity channels trace the main disc of the galaxy. The synthetic data also show kinks in the velocity across the spiral arms, evidence of non-circular velocities. These are still present even when we blur our data to a similar resolution as the observations, but largely absent in M31 and M33, indicating those galaxies do not contain significant spiral shocks. Thus the detailed velocity structure of our maps better represent previous observations of the grand design spiral M81. △ Less

Submitted 19 April, 2010; originally announced April 2010.

Comments: MNRAS accepted; 13 pages; 11 figures, 4 in colour

Tracing the Universe's Most Abundant Atom with the World's Largest Filled-Aperture Telescope

Authors: Kevin A. Douglas

Abstract: Among present-day observatories, the Arecibo Radio Telescope represents an extension of Galileo's vision to its logical extreme. With a diameter of 305 metres and state-of-the-art instrumentation, Arecibo continues to build on its legacy of world-class scientific achievement in radio astronomy. This paper highlights milestones in the remarkable history of this telescope, and also discusses curre… ▽ More Among present-day observatories, the Arecibo Radio Telescope represents an extension of Galileo's vision to its logical extreme. With a diameter of 305 metres and state-of-the-art instrumentation, Arecibo continues to build on its legacy of world-class scientific achievement in radio astronomy. This paper highlights milestones in the remarkable history of this telescope, and also discusses current surveys that are imaging the hydrogen content of our Milky Way Galaxy, and far beyond, in unprecedented detail. △ Less

Submitted 26 January, 2010; originally announced January 2010.

Comments: 12 pages, to appear in the proceedings of the conference "Astronomy and its Instruments Before and After Galileo" held in Venice, IT, 28 September - 2 October 2009

I-GALFA: The Inner-Galaxy ALFA Low-Latitude H I Survey

Authors: Bon-Chul Koo, Steven J. Gibson, Ji-hyun Kang, Kevin A. Douglas, Geumsook Park, Joshua E. G. Peek, Eric J. Korpela, Carl E. Heiles, Thomas M. Bania

Abstract: The I-GALFA survey is mapping all the HI in the inner Galactic disk visible to the Arecibo 305m telescope within 10 degrees of the Galactic plane (longitudes of 32 to 77 deg at b=0 deg). The survey, which will obtain 1.3 million independent spectra, became possible with the installation of the 7-beam Arecibo L-Band Feed Array (ALFA) receiver in 2004. ALFA's 3.4 arcmin resolution and tremendous s… ▽ More The I-GALFA survey is mapping all the HI in the inner Galactic disk visible to the Arecibo 305m telescope within 10 degrees of the Galactic plane (longitudes of 32 to 77 deg at b=0 deg). The survey, which will obtain 1.3 million independent spectra, became possible with the installation of the 7-beam Arecibo L-Band Feed Array (ALFA) receiver in 2004. ALFA's 3.4 arcmin resolution and tremendous sensitivity offer a great opportunity to observe the fine details of HI in the Galaxy. The I-GALFA survey began in 2008 May and will be completed in 2009 September. The data will be made publicly available when the calibrated and gridded cubes are completed. Further information on the I-GALFA project may be found at www.naic.edu/~galfa. △ Less

Submitted 11 December, 2009; originally announced December 2009.

Comments: 1 page, 2 figures, to appear in proceedings of IAU 2009 Special Session "The Galactic Plane - in Depth and Across the Spectrum"

Synthetic Observations of the HI Line in SPH-Simulated Spiral Galaxies

Authors: Kevin A. Douglas, David Acreman, Clare Dobbs, Chris Brunt

Abstract: Using the radiative transfer code Torus, we produce spectral-line cubes of the predicted HI profile from global SPH simulations of spiral galaxies. Torus grids the SPH galaxy using Adaptive Mesh Refinement, then applies a ray-tracing method to infer the HI profile along the line(s) of sight. The gridded galaxy can be observed from any direction, which enables us to model the observed HI profile… ▽ More Using the radiative transfer code Torus, we produce spectral-line cubes of the predicted HI profile from global SPH simulations of spiral galaxies. Torus grids the SPH galaxy using Adaptive Mesh Refinement, then applies a ray-tracing method to infer the HI profile along the line(s) of sight. The gridded galaxy can be observed from any direction, which enables us to model the observed HI profile for galaxies of any orientation. We can also place the observer inside the galaxy, to simulate HI observations taken from the Earth's position in the Milky Way. △ Less

Submitted 1 September, 2009; originally announced September 2009.

Comments: 4 pages, 2 figures, conference proceedings for "Panoramic Radio Astronomy: 1-2 Ghz Research on Galaxy Evolution" June 2-5, 2009 Groningen

A peculiar HI cloud near the distant globular cluster Pal 4

Authors: Jacco Th. van Loon, Snezana Stanimirovic, Mary Putman, Joshua E. G. Peek, Steven J. Gibson, Kevin A. Douglas, Eric J. Korpela

Abstract: We present 21-cm observations of four Galactic globular clusters, as part of the on-going GALFA-HI Survey at Arecibo. We discovered a peculiar HI cloud in the vicinity of the distant (109 kpc) cluster Pal 4, and discuss its properties and likelihood of association with the cluster. We conclude that an association of the HI cloud and Pal 4 is possible, but that a chance coincidence between Pal 4… ▽ More We present 21-cm observations of four Galactic globular clusters, as part of the on-going GALFA-HI Survey at Arecibo. We discovered a peculiar HI cloud in the vicinity of the distant (109 kpc) cluster Pal 4, and discuss its properties and likelihood of association with the cluster. We conclude that an association of the HI cloud and Pal 4 is possible, but that a chance coincidence between Pal 4 and a nearby compact high-velocity cloud cannot be ruled out altogether. New, more stringent upper limits were derived for the other three clusters: M 3, NGC 5466, and Pal 13. We briefly discuss the fate of globular cluster gas and the interaction of compact clouds with the Galactic Halo gas. △ Less

Submitted 13 March, 2009; originally announced March 2009.

Comments: Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.396:1096-1105,2009

The Disruption and Fueling of M33

Authors: M. E. Putman, J. E. G. Peek, A. Muratov, O. Y. Gnedin, W. Hsu, K. A. Douglas, C. Heiles, S. Stanimirovic, E. J. Korpela, S. J. Gibson

Abstract: The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of M33 (M_HI(tot) =1.4 x 10^9 Msun) is found beyond the star forming disk as traced in the far-ultraviolet (FU… ▽ More The disruption of the M33 galaxy is evident from its extended gaseous structure. We present new data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey that show the full extent and detailed spatial and kinematic structure of M33's neutral hydrogen. Over 18% of the HI mass of M33 (M_HI(tot) =1.4 x 10^9 Msun) is found beyond the star forming disk as traced in the far-ultraviolet (FUV). The most distinct features are extended warps, an arc from the northern warp to the disk, diffuse gas surrounding the galaxy, and a southern cloud with a filament back to the galaxy. The features extend out to 22 kpc from the galaxy center (18 kpc from the edge of the FUV disk) and the gas is directly connected to M33's disk. Only five discrete clouds (i.e., gas not directly connected to M33 in position-velocity space) are catalogued in the vicinity of M33, and these clouds show similar properties to Galactic and M31 halo clouds. M33's gaseous features most likely originate from the tidal disruption of M33 by M31 1-3 Gyr ago as shown from an orbit analysis which results in a tidal radius < 15 kpc in the majority of M33's possible orbits. M33 is now beyond the disruptive gravitational influence of M31 and the gas appears to be returning to M33's disk and redistributing its star formation fuel. M33's high mean velocity dispersion in the disk (~18.5 km/s) may also be consistent with the previous interaction and high rate of star formation. M33 will either exhaust its star formation fuel in the next few Gyrs or eventually become star formation fuel for M31. The latter represents the accretion of a large gaseous satellite by a spiral galaxy, similar to the Magellanic Clouds' relationship to the Galaxy. △ Less

Submitted 23 December, 2009; v1 submitted 12 December, 2008; originally announced December 2008.

Comments: published in ApJ (13 figures, several with color), see published version for full resolution figures; data available at http://sites.google.com/site/galfahi/data

Journal ref: Astrophys.J.703:1486-1501,2009

Low-Velocity Halo Clouds

Authors: J. E. G. Peek, Carl Heiles, M. E. Putman, Kevin A. Douglas

Abstract: Models that reproduce the observed high-velocity clouds (HVCs) also predict clouds at lower radial velocities that may easily be confused with Galactic disk (|z| < 1 kpc) gas. We describe the first search for these low-velocity halo clouds (LVHCs) using IRAS data and the initial data from the Galactic Arecibo L-band Feed Array survey in HI (GALFA-HI). The technique is based upon the expectation… ▽ More Models that reproduce the observed high-velocity clouds (HVCs) also predict clouds at lower radial velocities that may easily be confused with Galactic disk (|z| < 1 kpc) gas. We describe the first search for these low-velocity halo clouds (LVHCs) using IRAS data and the initial data from the Galactic Arecibo L-band Feed Array survey in HI (GALFA-HI). The technique is based upon the expectation that such clouds should, like HVCs, have very limited infrared thermal dust emission as compared to their HI column density. We describe our 'displacement-map' technique for robustly determining the dust-to-gas ratio of clouds and the associated errors that takes into account the significant scatter in the infrared flux from the Galactic disk gas. We find that there exist lower-velocity clouds that have extremely low dust-to-gas ratios, consistent with being in the Galactic halo - candidate LVHCs. We also confirm the lack of dust in many HVCs with the notable exception of complex M, which we consider to be the first detection of warm dust in HVCs. We do not confirm the previously reported detection of dust in complex C. In addition, we find that most Intermediate- and Low-Velocity clouds that are part of the Galactic disk have a higher 60 micron/100 micron flux ratio than is typically seen in Galactic HI, which is consistent with a previously proposed picture in which fast-moving Galactic clouds have smaller, hotter dust grains. △ Less

Submitted 22 October, 2008; originally announced October 2008.

Comments: 30 pages, 7 figures. Accepted to the ApJ

Journal ref: Astrophys.J.692:827-838,2009