As the complexity and volume of datasets have increased along with the capabilities of modular, open-source, easy-to-implement, visualization tools, scientists need for, and appreciation of, data visualization has risen too. Until recently, scientist
s thought of the explanatory graphics created at a research projects conclusion as pretty pictures needed only for journal publication or public outreach. The plots and displays produced during a research project - often intended only for experts - were thought of as a separate category, what we here call exploratory visualization. In this view, discovery comes from exploratory visualization, and explanatory visualization is just for communication. Our aim in this paper is to spark conversation amongst scientists, computer scientists, outreach professionals, educators, and graphics and perception experts about how to foster flexible data visualization practices that can facilitate discovery and communication at the same time. We present an example of a new finding made using the glue visualization environment to demonstrate how the border between explanatory and exploratory visualization is easily traversed. The linked-view principles as well as the actual code in glue are easily adapted to astronomy, medicine, and geographical information science - all fields where combining, visualizing, and analyzing several high-dimensional datasets yields insight. Whether or not scientists can use such a flexible undisciplined environment to its fullest potential without special training remains to be seen. We conclude with suggestions for improving the training of scientists in visualization practices, and of computer scientists in the iterative, non-workflow-like, ways in which modern science is carried out.
In this paper, I present a new set of synthetic spectral energy distributions (SEDs) for young stellar objects (YSOs) spanning a wide range of evolutionary stages, from the youngest deeply embedded protostars to pre-main-sequence stars with few or no
disks. These models include significant improvements on the previous generation of published models: in particular, the new models cover a much wider and more uniform region of parameter space, do not include highly model-dependent parameters, and include a number of improvements that make them more suited to modeling far-infrared and sub-mm observations of forming stars. Rather than all being part of a single monolithic set of models, the new models are split up into sets of varying complexity. The aim of the new set of models is not to provide the most physically realistic models for young stars, but rather to provide deliberately simplified models for initial modeling, which allows a wide range of parameter space to be explored. I present the design of the model set, and show examples of fitting these models to real observations to show how the new grid design can help us better understand what can be determined from limited unresolved observations. The models, as well as a Python-based fitting tool are publicly available to the community.
We analyze near-infrared UKIDSS observations of a sample of 8325 objects taken from a catalog of intrinsically red sources in the Galactic plane selected in the Spitzer-GLIMPSE survey. Given the differences in angular resolution (factor >2 better in
UKIDSS), our aim is to investigate whether there are multiple UKIDSS sources that might all contribute to the GLIMPSE flux, or there is only one dominant UKIDSS counterpart. We then study possible corrections to estimates of the SFR based on counts of GLIMPSE young stellar objects (YSOs). This represents an exploratory work towards the construction of a hierarchical YSO catalog. After performing PSF fitting photometry in the UKIDSS data, we implemented a technique to automatically recognize the dominant UKIDSS sources by evaluating their match with the spectral energy distribution (SED) of the associated GLIMPSE red sources. This is a generic method which could be robustly applied for matching SEDs across gaps at other wavelengths. We found that most (87.0% +- 1.6%) of the candidate YSOs from the GLIMPSE red source catalog have only one dominant UKIDSS counterpart which matches the mid-infrared SED (fainter associated UKIDSS sources might still be present). Though at first sight this could seem surprising, given that YSOs are typically in clustered environments, we argue that within the mass range covered by the GLIMPSE YSO candidates (intermediate to high masses), clustering with objects with comparable mass is unlikely at the GLIMPSE resolution. Indeed, by performing simple clustering experiments based on a population synthesis model of Galactic YSOs, we found that although ~60% of the GLIMPSE YSO enclose at least two UKIDSS sources, in general only one dominates the flux. No significant corrections are needed for estimates of the SFR of the Milky Way based on the assumption that the GLIMPSE YSOs are individual objects. (Abridged)
Through an extensive set of realistic synthetic observations (produced in Paper I), we assess in this part of the paper series (Paper III) how the choice of observational techniques affects the measurement of star-formation rates (SFRs) in star-formi
ng regions. We test the accuracy of commonly used techniques and construct new methods to extract the SFR, so that these findings can be applied to measure the SFR in real regions throughout the Milky Way. We investigate diffuse infrared SFR tracers such as those using 24 {mu}m, 70 {mu}m and total infrared emission, which have been previously calibrated for global galaxy scales. We set up a toy model of a galaxy and show that the infrared emission is consistent with the intrinsic SFR using extra-galactic calibrated laws (although the consistency does not prove their reliability). For local scales, we show that these techniques produce completely unreliable results for single star-forming regions, which are governed by different characteristic timescales. We show how calibration of these techniques can be improved for single star-forming regions by adjusting the characteristic timescale and the scaling factor and give suggestions of new calibrations of the diffuse star-formation tracers. We show that star-forming regions that are dominated by high-mass stellar feedback experience a rapid drop in infrared emission once high-mass stellar feedback is turned on, which implies different characteristic timescales. Moreover, we explore the measured SFRs calculated directly from the observed young stellar population. We find that the measured point sources follow the evolutionary pace of star formation more directly than diffuse star-formation tracers.
We use a large data-set of realistic synthetic observations (PaperI) to assess how observational techniques affect the measurement of physical properties of star-forming regions. In this paper (PaperII), we explore the reliability of the measured tot
al gas mass, dust surface density and dust temperature maps derived from modified blackbody fitting of synthetic Herschel observations. We found from our pixel-by-pixel analysis of the measured dust surface density and dust temperature a worrisome error spread especially close to star-formation sites and low-density regions, where for those contaminated pixels the surface densities can be under/overestimated by up to three orders of magnitude. In light of this, we recommend to treat the pixel-based results from this technique with caution in regions with active star formation. In regions of high background typical in the inner Galactic plane, we are not able to recover reliable surface density maps of individual synthetic regions, since low-mass regions are lost in the FIR background. When measuring the total gas mass of regions in moderate background, we find that modified blackbody fitting works well (absolute error:+9%;-13%) up to 10kpc distance (errors increase with distance). Commonly, the initial images are convolved to the largest common beam-size, which smears contaminated pixels over large areas. The resulting information loss makes this commonly-used technique less verifiable as now chi^2-values cannot be used as a quality indicator of a fitted pixel. Our control measurements of the total gas mass (without the step of convolution to the largest common beam size) produce similar results (absolute error:+20%;-7%) while having much lower median errors especially for the high-mass stellar feedback phase. In upcoming papers (III&IV) we test the reliability of measured star-formation rate with direct and indirect techniques.
We present a catalog of 840 X-ray sources and first results from a 100 ks Chandra X-ray Observatory imaging study of the filamentary infrared dark cloud G014.225$-$00.506, which forms the central regions of a larger cloud complex known as the M17 sou
thwest extension (M17 SWex). In addition to the rich population of protostars and young stellar objects with dusty circumstellar disks revealed by Spitzer Space Telescope archival data, we discover a population of X-ray-emitting, intermediate-mass pre--main-sequence stars (IMPS) that lack infrared excess emission from circumstellar disks. We model the infrared spectral energy distributions of this source population to measure its mass function and place new constraints on the inner dust disk destruction timescales for 2-8 $M_{odot}$ stars. We also place a lower limit on the star formation rate (SFR) and find that it is quite high ($dot{M}ge 0.007~M_{odot}$ yr$^{-1}$), equivalent to several Orion Nebula Clusters in G14.225$-$0.506 alone, and likely accelerating. The cloud complex has not produced a population of massive, O-type stars commensurate with its SFR. This absence of very massive (${ge}20~M_{odot}$) stars suggests that either (1) M17 SWex is an example of a distributed mode of star formation that will produce a large OB association dominated by intermediate-mass stars but relatively few massive clusters, or (2) the massive cores are still in the process of accreting sufficient mass to form massive clusters hosting O stars.
Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region, we assess how the choice of observational techniques affects the measurements of properties which trace star formation. Testing and calibrating observat
ional measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore different techniques for how to map the distributions of densities and temperatures from the particle-based simulations onto a Voronoi mesh suitable for radiative transfer and consequently explore their accuracy. We further test different ways to set up the radiative transfer in order to produce realistic synthetic observations. We give a detailed description of all methods and ultimately recommend techniques. We have found that the flux around 20 microns is strongly overestimated when blindly coupling the dust radiative transfer temperature with the hydrodynamical gas temperature. We find that when instead assuming a constant background dust temperature in addition to the radiative transfer heating, the recovered flux is consistent with actual observations. We present around 5800 realistic synthetic observations for Spitzer and Herschel bands, at different evolutionary time-steps, distances and orientations. In the upcoming papers of this series (Paper II, Paper III and Paper IV), we will test and calibrate measurements of the star-formation rate (SFR), gas mass and the star-formation efficiency (SFE) using our realistic synthetic observations.
We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3 resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm
(source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.
In contrast to most other galaxies, star-formation rates in the Milky Way can be estimated directly from Young Stellar Objects (YSOs). In the Central Molecular Zone (CMZ) the star-formation rate calculated from the number of YSOs with 24 microns emis
sion is up to order of magnitude higher than the value estimated from methods based on diffuse emission (such as free-free emission). Whether this effect is real or whether it indicates problems with either or both star formation rate measures is not currently known. In this paper, we investigate whether estimates based on YSOs could be heavily contaminated by more evolved objects such as main-sequence stars. We present radiative transfer models of YSOs and of main-sequence stars in a constant ambient medium which show that the main-sequence objects can indeed mimic YSOs at 24 microns. However, we show that in some cases the main-sequence models can be marginally resolved at 24 microns, whereas the YSO models are always unresolved. Based on the fraction of resolved MIPS 24 microns sources in the sample of YSOs previously used to compute the star formation rate, we estimate the fraction of misclassified YSOs to be at least 63%, which suggests that the star-formation rate previously determined from YSOs is likely to be at least a factor of three too high.
We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Fle
xible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity
