Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Pandeia is the exposure time calculator (ETC) system developed for the James Webb Space Telescope (JWST) that will be used for creating JWST proposals. It includes a simulation-hybrid Python engine that calculates the two-dimensional pixel-by-pixel signal and noise properties of the JWST instruments. This allows for appropriate handling of realistic point spread functions, MULTIACCUM detector readouts, correlated detector readnoise, and multiple photometric and spectral extraction strategies. Pandeia includes support for all the JWST observing modes, including imaging, slitted/slitless spectroscopy, integral field spectroscopy, and coronagraphy. Its highly modular, data-driven design makes it easily adaptable to other observatories. An implementation for use with WFIRST is also available.

With the prospect of the next generation of ground-based telescopes, the extremely large telescopes (ELTs), increasingly complex and demanding adaptive optics (AO) systems are needed. This is to compensate for image distortion caused by atmospheric turbulence and fully take advantage of mirrors with diameters of 30 to 40 m. This requires a more precise characterization of the turbulence. The PML (Profiler of Moon Limb) was developed within this context. The PML aims to provide high-resolution altitude profiles of the turbulence using differential measurements of the Moon limb position to calculate the transverse spatio-angular covariance of the Angle of Arrival fluctuations. The covariance of differential image motion for different separation angles is sensitive to the altitude distribution of the seeing. The use of the continuous Moon limb provides a large number of separation angles allowing for the high-resolution altitude of the profiles. The method is presented and tested with simulated data. Moreover a PML instrument was deployed at the Sutherland Observatory in South Africa in August 2011. We present here the results of this measurement campaign.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present ultraviolet, optical, and near-infrared observations of SN 2012ap, a broad-lined Type Ic supernova in the galaxy NGC 1729 that produced a relativistic and rapidly decelerating outflow without a gamma-ray burst signature. Photometry and spectroscopy follow the flux evolution from -13 to +272 days past the B-band maximum of -17.4 +/- 0.5 mag. The spectra are dominated by Fe II, O I, and Ca II absorption lines at ejecta velocities of 20,000 km/s that change slowly over time. Other spectral absorption lines are consistent with contributions from photospheric He I, and hydrogen may also be present at higher velocities (> 27,000 km/s). We use these observations to estimate explosion properties and derive a total ejecta mass of 2.7 Msolar, a kinetic energy of 1.0x10^{52} erg, and a 56Ni mass of 0.1-0.2 Msolar. Nebular spectra (t > 200d) exhibit an asymmetric double-peaked [OI] 6300,6364 emission profile that we associate with absorption in the supernova interior, although toroidal ejecta geometry is an alternative explanation. SN 2012ap joins SN 2009bb as another exceptional supernova that shows evidence for a central engine (e.g., black-hole accretion or magnetar) capable of launching a non-negligible portion of ejecta to relativistic velocities without a coincident gamma-ray burst detection. Defining attributes of their progenitor systems may be related to notable properties including above-average environmental metallicities of Z > Zsolar, moderate to high levels of host-galaxy extinction (E(B-V) > 0.4 mag), detection of high-velocity helium at early epochs, and a high relative flux ratio of [Ca II]/[O I] > 1 at nebular epochs. These events support the notion that jet activity at various energy scales may be present in a wide range of supernovae.

The diffuse interstellar bands (DIBs) are absorption features observed in optical and near-infrared spectra that are thought to be associated with carbon-rich polyatomic molecules in interstellar gas. However, because the central wavelengths of these bands do not correspond with electronic transitions of any known atomic or molecular species, their nature has remained uncertain since their discovery almost a century ago. Here we report on unusually strong DIBs in optical spectra of the broad-lined Type Ic supernova SN 2012ap that exhibit changes in equivalent width over short (~30 days) timescales. The 4428 and 6283 Angstrom DIB features get weaker with time, whereas the 5780 Angstrom feature shows a marginal increase. These nonuniform changes suggest that the supernova is interacting with a nearby source of the DIBs and that the DIB carriers possess high ionization potentials, such as small cations or charged fullerenes. We conclude that moderate-resolution spectra of supernovae with DIB absorptions obtained within weeks of outburst could reveal unique information about the mass-loss environment of their progenitor systems and provide new constraints on the properties of DIB carriers.

Observations spanning a large wavelength range, from X-ray to radio, of the Type IIb supernova 2011hs are presented, covering its evolution during the first year after explosion. The optical light curve presents a narrower shape and a fainter luminosity at peak than previously observed for Type IIb SNe. High expansion velocities are measured from the broad absorption H I and He I lines. From the comparison of the bolometric light curve and the time evolution of the photospheric velocities with hydrodynamical models, we found that SN 2011hs is consistent with the explosion of a 3-4 Msun He-core progenitor star, corresponding to a main sequence mass of 12-15 Msun, that ejected a mass of 56Ni of about 0.04 Msun, with an energy of E= 8.5 x 10^50 erg. Such a low-mass progenitor scenario is in full agreement with the modelling of the nebular spectrum taken at $\sim$215 days from maximum. From the modelling of the adiabatic cooling phase, we infer a progenitor radius of $\approx$500-600 Rsun, clearly pointing to an extended progenitor star. The radio light curve of SN 2011hs yields a peak luminosity similar to that of SN 1993J, but with a higher mass loss rate and a wind density possibly more similar to that of SN 2001ig. Although no significant deviations from a smooth decline have been found in the radio light curves, we cannot rule out the presence of a binary companion star.

We present results from the first year of a campaign to characterize and monitor the optical turbulence profile at the SAAO Sutherland observing station in South Africa. A MASS-DIMM (MultiAperture Scintillation Sensor - Differential Image Motion Monitor) was commissioned in March 2010 to provide continuous monitoring of the seeing conditions. Over the first month of the campaign, a SLODAR (SLOpe Detection And Ranging) from Durham University was also installed allowing an independent verification of the performance of the MASS-DIMM device. After the first year of data collection, the overall median seeing value is found to be 1.32" as measured at ground level. The ground layer which includes all layers below 1 km accounts for 84% of the turbulence, while the free atmosphere above 1 km accounts for 16% with a median value of 0.41". The median isoplanatic angle value is 1.92", which is similar to other major astronomical sites. The median coherence time, calculated from corrected MASS measurements, is 2.85 ms. The seeing conditions at the site do show a strong correlation with wind direction with bad seeing conditions being associated with winds from the South-East.

We present 65 optical spectra of the Type Ia supernova SN 2012fr, of which 33 were obtained before maximum light. At early times SN 2012fr shows clear evidence of a high-velocity feature (HVF) in the Si II 6355 line which can be cleanly decoupled from the lower velocity "photospheric" component. This Si II 6355 HVF fades by phase -5; subsequently, the photospheric component exhibits a very narrow velocity width and remains at a nearly constant velocity of v~12,000 km/s until at least 5 weeks after maximum brightness. The Ca II infrared (IR) triplet exhibits similar evidence for both a photospheric component at v~12,000 km/s with narrow line width and long velocity plateau, as well as a high-velocity component beginning at v~31,000 km/s two weeks before maximum. SN 2012fr resides on the border between the "shallow silicon" and "core-normal" subclasses in the Branch et al. (2009) classification scheme, and on the border between normal and "high-velocity" SNe Ia in the Wang et al. (2009a) system. Though it is a clear member of the "low velocity gradient" (LVG; Benetii et al., 2005) group of SNe Ia and exhibits a very slow light-curve decline, it shows key dissimilarities with the overluminous SN 1991T or SN 1999aa subclasses of SNe Ia. SN 2012fr represents a well-observed SN Ia at the luminous end of the normal SN Ia distribution, and a key transitional event between nominal spectroscopic subclasses of SNe Ia.

WISE J104915.57$-$531906.1 is a L/T brown dwarf binary located 2pc from the Sun. The pair contains the closest known brown dwarfs and is the third closest known system, stellar or sub-stellar. We report comprehensive follow-up observations of this newly uncovered system. We have determined the spectral types of both components (L8+/-1, for the primary, agreeing with the discovery paper; T1.5+/-2 for the secondary, which was lacking spectroscopic type determination in the discovery paper) and, for the first time, their radial velocities (V_rad~23.1, 19.5 km/s) using optical spectra obtained at the Southern African Large Telescope (SALT) and other facilities located at the South African Astronomical Observatory (SAAO). The relative radial velocity of the two components is smaller than the range of orbital velocities for theoretically predicted masses, implying that they form a gravitationally bound system. We report resolved near-infrared JHK_S photometry from the IRSF telescope at the SAAO which yields colors consistent with the spectroscopically derived spectral types. The available kinematic and photometric information excludes the possibility that the object belongs to any of the known nearby young moving groups or associations. Simultaneous optical polarimetry observations taken at the SAAO 1.9-m give a non-detection with an upper limit of 0.07%. For the given spectral types and absolute magnitudes, 1Gyr theoretical models predict masses of 0.04--0.05 M_odot for the primary, and 0.03--0.05 M_odot for the secondary.

We present results of a multi-site photometric campaign on the high-amplitude $\delta$\,Scuti star KIC\,6382916 in the {\it Kepler} field. The star was observed over a 85-d interval at five different sites in North America and Europe during 2011. {\it Kepler} photometry and ground-based multicolour light curves of KIC\,6382916 are used to investigate the pulsational content and to identify the principal modes. High-dispersion spectroscopy was also obtained in order to derive the stellar parameters and projected rotational velocity. From an analysis of the {\it Kepler} time series, three independent frequencies and a few hundred combination frequencies are found. The light curve is dominated by two modes with frequencies $f_{1}$= 4.9107 and $f_{2}$= 6.4314\,d$^{-1}$. The third mode with $f_{3}$= 8.0350\,d$^{-1}$ has a much lower amplitude. We attempt mode identification by examining the amplitude ratios and phase differences in different wavebands from multicolour photometry and comparing them to calculations for different spherical harmonic degree, $l$. We find that the theoretical models for $f_1$ and $f_2$ are in a best agreement with the observations and lead to value of l = 1 modes, but the mode identification of $f_3$ is uncertain due to its low amplitude. Non-adiabatic pulsation models show that frequencies below 6\,d$^{-1}$ are stable, which means that the low frequency of $f_1$ cannot be reproduced. This is further confirmation that current models predict a narrower pulsation frequency range than actually observed.

We present X-ray, UV/optical, and radio observations of the stripped-envelope, core-collapse supernova (SN) 2011ei, one of the least luminous SNe IIb or Ib observed to date. Our observations begin with a discovery within 1 day of explosion and span several months afterward. Early optical spectra exhibit broad, Type II-like hydrogen Balmer profiles that subside rapidly and are replaced by Type Ib-like He-rich features on the timescale of one week. High-cadence monitoring of this transition suggests that absorption attributable to a high velocity (> 12,000 km/s) H-rich shell is not rare in Type Ib events. Radio observations imply a shock velocity of v = 0.13c and a progenitor star mass-loss rate of 1.4 x 10^{-5} Msun yr^{-1} (assuming wind velocity v_w=10^3 km/s). This is consistent with independent constraints from deep X-ray observations with Swift-XRT and Chandra. Overall, the multi-wavelength properties of SN 2011ei are consistent with the explosion of a lower-mass (3-4 Msun), compact (R* <= 1x10^{11} cm), He core star. The star retained a thin hydrogen envelope at the time of explosion, and was embedded in an inhomogeneous circumstellar wind suggestive of modest episodic mass-loss. We conclude that SN 2011ei's rapid spectral metamorphosis is indicative of time-dependent classifications that bias estimates of explosion rates for Type IIb and Ib objects, and that important information about a progenitor star's evolutionary state and mass-loss immediately prior to SN explosion can be inferred from timely multi-wavelength observations.

A considerable fraction of the central stars of planetary nebulae (CSPNe) are hydrogen-deficient. Almost all of these H-deficient central stars (CSs) display spectra with strong carbon and helium lines. Most of them exhibit emission line spectra resembling those of massive WC stars. Therefore these stars are classed as CSPNe of spectral type [WC]. Recently, quantitative spectral analysis of two emission-line CSs, PB 8 and IC 4663, revealed that these stars do not belong to the [WC] class. Instead PB 8 has been classified as [WN/WC] type and IC 4663 as [WN] type. In this work we report the spectroscopic identification of another rare [WN] star, the CS of Abell 48. We performed a spectral analysis of Abell 48 with the Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres. We find that the expanding atmosphere of Abell 48 is mainly composed of helium (85 per cent by mass), hydrogen (10 per cent), and nitrogen (5 per cent). The residual hydrogen and the enhanced nitrogen abundance make this object different from the other [WN] star IC 4663. We discuss the possible origin of this atmospheric composition.

We report on a multi-site photometric campaign on the high-amplitude $\delta$ Scuti star V2367 Cyg in order to determine the pulsation modes. We also used high-dispersion spectroscopy to estimate the stellar parameters and projected rotational velocity. Time series multicolour photometry was obtained during a 98-d interval from five different sites. These data were used together with model atmospheres and non-adiabatic pulsation models to identify the spherical harmonic degree of the three independent frequencies of highest amplitude as well as the first two harmonics of the dominant mode. This was accomplished by matching the observed relative light amplitudes and phases in different wavebands with those computed by the models. In general, our results support the assumed mode identifications in a previous analysis of Kepler data.

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) stars via detection of circular shells (typical of known confirmed and cLBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 um in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered cLBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R- and I-bands by 0.68\pm0.10 mag and 0.61\pm0.04 mag, respectively, during the last 13-18 years, while the second one (known as Hen 3-1383) varies its B,V,R,I and K_s brightnesses by \simeq 0.5-0.9 mag on time-scales from 10 days to decades. We also found significant changes in the spectrum of Hen 3-1383 on a timescale of \simeq 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered cLBVs might be directly related to their possible runaway status.

The outburst and subsequent brightness evolution of Comet 17P/Holmes has been observed using the MMT Observatory's All-Sky Camera (Pickering 2006) on Mt. Hopkins near Tucson, Arizona, USA. The comet was picked up at the limiting visual magnitude of 5.5 on October 24.38 and tracked by the camera continuously until sunrise four hours later. During this time the comet brightened to visual magnitude 3.5. Comet Holmes was next observed just after sunset on October 25.23 at visual magnitude 2.5 where it remained approximately constant over the next three days. The comet then began to dim slowly and was followed into the early months of 2008 with periods of dense time coverage.

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of the Sloan Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with deeper imaging and spectroscopy on the Multiple Mirror Telescope. From this, we derive a luminosity of M_v = -2.5, a half-light radius of 34 pc and a systemic velocity of -40$ km/s. Our MMT data also provides evidence for a stream around Segue 2 at a similar heliocentric velocity, and the SEGUE data show that it is also present in neighboring fields. We resolve the velocity dispersion of Segue 2 as 3.4 km/s and the possible stream as about 7 km/s. This object shows points of comparison with other recent discoveries, Segue 1, Boo II and Coma. We speculate that all four objects may be representatives of a population of satellites of satellites -- survivors of accretion events that destroyed their larger but less dense parents. They are likely to have formed at redshifts z > 10 and are good candidates for fossils of the reionization epoch.

We examine the infrared properties of five low surface brightness galaxies (LSBGs) and compare them with related but higher surface brightness galaxies, using Spitzer Space Telescope images and spectra. All the LSBGs are detected in the 3.6 and 4.5um bands, representing the stellar population. All but one are detected at 5.8 and 8.0um, revealing emission from hot dust and aromatic molecules, though many are faint or point-like at these wavelengths. Detections of LSBGs at the far-infrared wavelengths, 24, 70, and 160um, are varied in morphology and brightness, with only two detections at 160um, resulting in highly varied spectral energy distributions. Consistent with previous expectations for these galaxies, we find that detectable dust components exist for only some LSBGs, with the strength of dust emission dependent on the existence of bright star forming regions. However, the far-infrared emission may be relatively weak compared with normal star-forming galaxies.

We present optical and IR spectroscopy of the first two months of evolution of the Type II SN 1999em. We combine these data with high-quality optical/IR photometry beginning only three days after shock breakout, in order to study the performance of the ``Expanding Photosphere Method'' (EPM) in the determination of distances. With this purpose we develop a technique to measure accurate photospheric velocities by cross-correlating observed and model spectra. The application of this technique to SN 1999em shows that we can reach an average uncertainty of 11% in velocity from an individual spectrum. Our analysis shows that EPM is quite robust to the effects of dust. In particular, the distances derived from the VI filters change by only 7% when the adopted visual extinction in the host galaxy is varied by 0.45 mag. The superb time sampling of the BVIZJHK light-curves of SN 1999em permits us to study the internal consistency of EPM and test the dilution factors computed from atmosphere models for Type II plateau supernovae. We find that, in the first week since explosion, the EPM distances are up to 50% lower than the average, possibly due the presence of circumstellar material. Over the following 65 days, on the other hand, our tests lend strong credence to the atmosphere models, and confirm previous claims that EPM can produce consistent distances without having to craft specific models to each supernova. This is particularly true for the VI filters which yield distances with an internal consistency of 4%. From the whole set of BVIZJHK photometry, we obtain an average distance of 7.5+/-0.5 Mpc, where the quoted uncertainty (7%) is a conservative estimate of the internal precision of the method obtained from the analysis of the first 70 days of the supernova evolution.

Upper limits for the disk mass-to-light ratios for the low surface brightness galaxies UGC 6614 and F568-6 (Malin 2) are estimated by considering the minimum velocity perturbations in the HI velocity field that should result from the spiral structure observed in the R band images. The weak observed response in the $\phi$ velocity component limits the mass-to-light ratios of the disk within a scale length to M/L <~ 3 and 6 for UGC 6614 for F568-6 respectively (in solar units) based upon azimuthal variations observed in the R band images. These limits are sufficiently strong to require a significant dark matter component even in the central regions of these galaxies. Our limits furthermore imply that this dark matter component cannot be in the form of a cold disk since a cold disk would necessarily be involved in the spiral structure. However, a more massive disk could be consistent with the observations because of a non-linear gas response or if the gas is driven by bar-like distortions instead of spiral structure. To produce the large observed arm/interarm HI density variations it is likely that the spiral arm potential perturbation is sufficiently strong to produce shocks in the gas. For a forcing that is greater than 2% of the axisymmetric force, M/L >~ 1 is required in both galaxies in the outer regions. These lower limits imply that the stellar surface density is at least of the same order as the gas surface density. This is consistent with the large scale morphology of the spiral structure, and the stability of the gas disk, both which suggest that a moderate stellar component is required to produce the observed spiral structure. (Shortened abstract)

We present near-Infrared images of low surface brightness galaxies: H band images of UGC 6614 and F568-6 (Malin 2). The H band images show spiral structure so we are confident that we are seeing the galaxy disks as well as the central bulges. The optical-IR colors of these galaxies, $R-H = 2.2\pm 0.2$, $B-H = 3.5-4.2$, are extremely red and are similar to those of many S0 and elliptical galaxies. This represents strong evidence for a component of old stars. Recent studies find that low surface brightness galaxies span a wide range of morphologies and colors. In this context UGC 6614 and F568-6 could be examples of low surface brightness galaxies which have undergone a past epoch of more vigorous star formation.