Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of protostars by following the gravitational collapse of 1~M$_\odot$ gas clouds through the first hydrostatic core to stellar densities. The clouds are imposed with both rotational and turbulent velocities, and are threaded with a magnetic field that is parallel/anti-parallel or perpendicular to the rotation axis; we investigate two rotation rates and four Mach numbers. The initial radius and mass of the stellar core are only weakly dependent on the initial parameters. In the models that include ideal MHD, the magnetic field strength implanted in the protostar at birth is much higher than observed, independent of the initial level of turbulence; only non-ideal MHD can reduce this strength to near or below the observed levels. This suggests that not only is ideal MHD an incomplete picture of star formation, but that the magnetic fields in low mass stars are implanted later in life by a dynamo process. Non-ideal MHD suppresses magnetically launched stellar core outflows, but turbulence permits thermally launched outflows to form a few years after stellar core formation.

Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1~M$_\odot$ gas cloud that is threaded with a magnetic field and imposed with both rotational and turbulent velocities. We investigate magnetic fields that are parallel/anti-parallel and perpendicular to the rotation axis, two rotation rates and four Mach numbers. Disc formation occurs preferentially in the models that include non-ideal MHD where the magnetic field is anti-parallel or perpendicular to the rotation axis. This is independent of the initial rotation rate and level of turbulence, suggesting that subsonic turbulence plays a minimal role in influencing the formation of discs. Aside from first core outflows which are influenced by the initial level of turbulence, non-ideal MHD processes are more important than turbulent processes during the formation of discs around low-mass stars.

As part of a study funded by NASA Headquarters, we are developing a Probe-class mission concept called the Cosmic Evolution Through UV Spectroscopy (CETUS). CETUS includes a 1.5-m aperture diameter telescope with a large field-of-view (FOV). CETUS includes three scientific instruments: a Far Ultraviolet (FUV) and Near Ultraviolet (NUV) imaging camera (CAM); a NUV Multi-Object Spectrograph (MOS); and a dual-channel Point Source Spectrograph (PSS) in the Lyman Ultraviolet (LUV), FUV, and NUV spectral regions. The large FOV Three Mirror Anastigmatic (TMA) Optical Telescope Assembly (OTA) simultaneously feeds the three separate scientific instruments. That is, the instruments view separate portions of the TMA image plane, enabling parallel operation of the three instruments. The field viewed by the MOS, whose design is based on an Offner-type spectrographic configuration to provide wide FOV correction, is actively configured to select and isolate numerous field sources using a next-generation Micro-Shutter Array (MSA). The two-channel camera design is also based on an Offner-like configuration. The Point Source Spectrograph (PSS) performs high spectral resolution spectroscopy on unresolved objects over the NUV region with spectral resolving power, R~ 40,000, in an echelle mode. The PSS also performs long-slit imaging spectroscopy at R~ 20,000 in the LUV and FUV spectral regions with two aberration-corrected, blazed, holographic gratings used in a Rowland-like configuration. The optical system also includes two Fine Guidance Sensors (FGS), and Wavefront Sensors (WFS) that sample numerous locations over the full OTA FOV. In-flight wavelength calibration is performed by a Wavelength Calibration System (WCS), and flat-fielding is also performed, both using in-flight calibration sources. This paper will describe the current optical design and the major trade studies leading to the design.

Pluto's surface is geologically complex because of volatile ices that are mobile on seasonal and longer time scales. Here we analyzed New Horizons LEISA spectral data to globally map the nitrogen ice, including nitrogen with methane diluted in it. Our goal was to learn about the seasonal processes influencing ice redistribution, to calculate the globally averaged energy balance, and to place a lower limit on Pluto's N2 inventory. We present the average latitudinal distribution of nitrogen and investigate the relationship between its distribution and topography on Pluto by using maps that include the shifted bands of methane in solid solution with nitrogen to more completely map the distribution of the nitrogen ice. We find that the global average bolometric albedo is 0.83 +\- 0.11, similar to that inferred for Triton, and that a significant fraction of Pluto's N2 is stored in Sputnik Planitia. Under the assumption that Pluto's nitrogen-dominated 11.5 microbar atmosphere is in vapor pressure equilibrium with the nitrogen ice, the ice temperature is 36.93 +/- 0.10 K, as measured by New Horizons. Combined with our global energy balance calculation, this implies that the average bolometric emissivity of Pluto's nitrogen ice is probably in the range 0.47 - 0.72. This is consistent with the low emissivities estimated for Triton based on Voyager, and may have implications for Pluto's atmospheric seasonal variations, as discussed below. The global pattern of volatile transport at the time of the encounter was from north to south, and the transition between condensation and sublimation within Sputnik Planitia is correlated with changes in the grain size and CH4 concentration derived from the spectral maps. The low emissivity of Pluto's N2 ice suggests that Pluto's atmosphere may undergo an extended period of constant pressure even as Pluto recedes from the Sun in its orbit.

Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade.

The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background.

Several post-detection approaches to the mitigation of radio-frequency interference (RFI) are compared by applying them to the strong RFI from the Iridium satellites. These provide estimates for the desired signal in the presence of RFI, by exploiting distinguishing characteristics of the RFI, such as its polarization, statistics, and periodicity. Our data are dynamic spectra with full Stokes parameters and 1 ms time resolution. Moreover, since most man-made RFI is strongly polarized, we use the data to compare its unpolarized component with its Stokes I. This approach on its own reduces the RFI intensity by many tens of dBs. A comprehensive approach that also recognizes non-Gaussian statistics, and the time and frequency structure inherent in the RFI permits exceedingly effective post-detection excision provided full Stokes intensity data are available.

Many spiral galaxies have extended outer H~I disks and display low levels of star formation, inferred from the far-ultraviolet emission detected by {\it GALEX}, well beyond the optical radius. Here, we investigate the supernova (SN) rate in the outskirts of galaxies, using the largest and most homogeneous set of nearby supernovae (SNe) from the Lick Observatory Supernova Search (LOSS). While SN rates have been measured with respect to various galaxy properties, such as stellar mass and metallicity, their relative frequency in the outskirts versus the inner regions has not yet been studied. Understanding the SN rate as a function of intragalactic environment has many ramifications, including the interpretation of LIGO observations, the formation of massive stars, and the puzzlingly high velocity dispersion of the outer H~I disk. Using data from the LOSS survey, we find that the rate beyond the optical radius of spiral galaxies is $2.5 \pm 0.5$ SNe per millennium, while dwarf galaxies host $ 4.0 \pm 2.2$ SNe per millennium. The rates of core-collapse SNe (that may collapse to form the massive black holes detected by LIGO/Virgo) in the outer disks of spirals is $1.5 \pm 0.15$ SNe per millennium and in dwarf galaxies is $2.6 \pm 1.5$ SNe per millennium. Core-collapse SNe in spiral outskirts contribute $7600 \pm 1700$\,SNe\,~Gpc$^{-3}$\,yr$^{-1}$ to the volumetric rate, and dwarf galaxies have a rate of $31,000 \pm 18,000$\,SNe\,Gpc$^{-3}$\,yr$^{-1}$. The relative ratio of core-collapse to Type Ia SNe is comparable in the inner and outer parts of spirals, and in dwarf galaxies.

We present the results of eighteen magnetohydrodynamical calculations of the collapse of a molecular cloud core to form a protostar. Some calculations include radiative transfer in the flux-limited diffusion approximation while others employ a barotropic equation of state. We cover a wide parameter space, with mass-to-flux ratios ranging from $\mu = 5$ to $20$; initial turbulent amplitudes ranging from a laminar calculation (i.e. where the Mach number, $\mathscr{M} = 0$) to transonic $\mathscr{M} = 1$; and initial rotation rates from $\beta_\mathrm{rot} = 0.005$ to $0.02$. We first show that using a radiative transfer scheme produces warmer pseudo-discs than the barotropic equation of state, making them more stable. We then `shake' the core by increasing the initial turbulent velocity field, and find that at all three mass-to-flux ratios transonic cores are weakly bound and do not produce pseudo-discs; $\mathscr{M} = 0.3$ cores produce very disrupted discs; and $\mathscr{M} = 0.1$ cores produce discs broadly comparable to a laminar core. In our previous paper (arXiv:1701.08741), we showed that a pseudo-disc coupled with sufficient magnetic field is necessary to form a bipolar outflow. Here we show that only weakly turbulent cores exhibit collimated jets. We finally take the $\mathscr{M} = 1.0$, $\mu = 5$ core and `stir' it by increasing the initial angular momentum, finding that once the degree of rotational energy exceeds the turbulent energy in the core the disc returns, with a corresponding (though slower), outflow. These conclusions place constraints on the initial mixtures of rotation and turbulence in molecular cloud cores which are conducive to the formation of bipolar outflows early in the star formation process.

We report results from twelve simulations of the collapse of a molecular cloud core to form one or more protostars, comprising three field strengths (mass-to-flux ratios, {\mu}, of 5, 10, and 20) and four field geometries (with values of the angle between the field and rotation axes, {\theta}, of 0{\deg}, 20{\deg}, 45{\deg}, and 90{\deg}), using a smoothed particle magnetohydrodynamics method. We find that the values of both parameters have a strong effect on the resultant protostellar system and outflows. This ranges from the formation of binary systems when {\mu} = 20 to strikingly differing outflow structures for differing values of {\theta}, in particular highly suppressed outflows when {\theta} = 90{\deg}. Misaligned magnetic fields can also produce warped pseudo-discs where the outer regions align perpendicular to the magnetic field but the innermost region re-orientates to be perpendicular to the rotation axis. We follow the collapse to sizes comparable to those of first cores and find that none of the outflow speeds exceed 8 km s$^{-1}$. These results may place constraints on both observed protostellar outflows, and also on which molecular cloud cores may eventually form either single stars and binaries: a sufficiently weak magnetic field may allow for disc fragmentation, whilst conversely the greater angular momentum transport of a strong field may inhibit disc fragmentation.

Obtaining a stable magnetohydrodynamical (MHD) formalism in SPH - i.e. smoothed particle magnetohydrodynamics (SPMHD) - has proven remarkably difficult. To implement MHD requires two steps: a modification to the momentum equation and an induction equation, and both present challenges. We first provide an overview of how SPMHD is implemented, and then discuss how this implementation fails and the limitation of various corrective methods - with particular reference to the effects of particle disorder. Although there are many problems for which, with careful choice of corrective measures, good results can be obtained, we then show that, at the very limits of the state of the art, the ability to perform stable MHD calculations in SPH is curtailed by numerical issues.

The equations of smoothed particle magnetohydrodynamics (SPMHD), even with the various corrections to instabilities so far proposed, have been observed to be unstable when a very steep density gradient is necessarily combined with a variable smoothing length formalism. Here we consider in more detail the modifications made to the SPMHD equations in LBP2015 that resolve this instability by replacing the smoothing length in the induction and anisotropic force equations with an average smoothing length term. We then explore the choice of average used and compare the effects on a test `cylinder-in-a-box' problem and the collapse of a magnetised molecular cloud core. We find that, aside from some benign numerical effects at low resolutions for the quadratic mean, the formalism is robust as to the choice of average but that in complicated models it is essential to apply the average to both equations; in particular, all four averages considered exhibit similar conservation properties. This improved formalism allows for arbitrarily small sink particles and field geometries to be explored, vastly expanding the range of astronomical problems that can be modeled using SPMHD.

We have developed a modified form of the equations of smoothed particle magnetohydrodynamics which are stable in the presence of very steep density gradients. Using this formalism, we have performed simulations of the collapse of magnetised molecular cloud cores to form protostars and drive outflows. Our stable formalism allows for smaller sink particles (< 5 AU) than used previously and the investigation of the effect of varying the angle, {\theta}, between the initial field axis and the rotation axis. The nature of the outflows depends strongly on this angle: jet-like outflows are not produced at all when {\theta} > 30{\deg}, and a collimated outflow is not sustained when {\theta} > 10{\deg}. No substantial outflows of any kind are produced when {\theta} > 60{\deg}. This may place constraints on the geometry of the magnetic field in molecular clouds where bipolar outflows are seen.

Photodissociation of 14N2 and 14N15N occurs in interstellar clouds, circumstellar envelopes, protoplanetary discs, and other environments due to UV radiation from stellar sources and the presence of cosmic rays. This source of N atoms initiates the formation of complex N-bearing species and influences their isotopic composition. To study the photodissociation rates of 14N15N by UV continuum radiation and both isotopologues in a field of cosmic ray induced photons. To determine the effect of these on the isotopic composition of more complex molecules. High-resolution photodissociation cross sections of N2 are used from an accurate and comprehensive quantum- mechanical model of the molecule based on laboratory experiments. A similarly high-resolution spectrum of H2 emission following interactions with cosmic rays has been constructed. The spectroscopic data are used to calculate dissociation rates which are input into isotopically differentiated chemical models, describing an interstellar cloud and a protoplanetary disc. The dissociation rate of 14N15N in a Draine field assuming 30K excitation is 1.73x10-10s-1 and the rate due to cosmic rays assuming an H2 ionisation rate of 10-16s-1 is about 10-15s-1, with up to a factor of 10 difference between isotopologues. Shielding functions for 14N15N by 14N2, H2, and H are presented. Incorporating these into an interstellar cloud model, an enhancement of the atomic 15N/14N ratio is obtained due to the self-shielding of external radiation at an extinction of about 1.5 mag. This effect is larger where grain growth has reduced the opacity of dust to ultraviolet radiation. The transfer of photolytic isotopic fractionation N2 to other molecules is significant in a disc model, and is species dependent with 15N enhancement approaching a factor of 10 for HCN.

Molecular nitrogen is one of the key species in the chemistry of interstellar clouds and protoplanetary disks and the partitioning of nitrogen between N and N2 controls the formation of more complex prebiotic nitrogen-containing species. The aim of this work is to gain a better understanding of the interstellar N2 photodissociation processes based on recent detailed theoretical and experimental work and to provide accurate rates for use in chemical models. We simulated the full high-resolution line-by-line absorption + dissociation spectrum of N2 over the relevant 912-1000 \AA\ wavelength range, by using a quantum-mechanical model which solves the coupled-channels Schr\"odinger equation. The simulated N2 spectra were compared with the absorption spectra of H2, H, CO, and dust to compute photodissociation rates in various radiation fields and shielding functions. The effects of the new rates in interstellar cloud models were illustrated for diffuse and translucent clouds, a dense photon dominated region and a protoplanetary disk.

The first comprehensive high-resolution photoabsorption spectrum of 14N15N has been recorded using the Fourier-transform spectrometer attached to the Desirs beamline at the Soleil synchrotron. Observations are made in the extreme ultraviolet (XUV) and span 100,000-109,000 cm-1 (100-91.7 nm). The observed absorption lines have been assigned to 25 bands and reduced to a set of transition energies, f values, and linewidths. This analysis has verified the predictions of a theoretical model of N2 that simulates its photoabsorption and photodissociation cross section by solution of an isotopomer independent formulation of the coupled-channel Schroedinger equation. The mass dependence of predissociation linewidths and oscillator strengths is clearly evident and many local perturbations of transition energies, strengths, and widths within individual rotational series have been observed.

The IRAS source, 19312+1950, exhibits SiO maser emission, which is predominantly detected in evolved stars enshrouded by a cold molecular envelope. In fact, the mojority of the observational properties of IRAS 19312+1950 is consistent with the nature of an asymptotic giant branch (AGB) star or post-AGB star. Interestingly, however, some of the observational properties cannot be readily explained within the standard scheme of stellar evolution, and those are rather reminiscent of young stellar objects. In the present research we considered the evolutionary status of IRAS 19312+1950 as revealed by the VLBI and MERLIN observations in SiO, H2O and OH maser lines. The double-peaked profile of the 22 GHz H2O maser line is clearly detected, with the emission regions of its red and blue-shifted components separately located, leaving a space of about 10.9 mas between them. The kinematic properties of H2O maser emission region appear to be more consistent with a bipolar flow rather than other interpretations such as the Keplerian rotation of a disk. The red-shifted component of the SiO maser emission, which exhibits a double-peak profile in previous single-dish observations, is clearly detected in the present interferometry, while the 1612 MHz OH maser line exhibits a complicated line profile consisting of a single strong peak and many weak, high-velocity spikes. The structure of OH maser emission region is partially resolved, and the kinematic properties of the OH maser emission region are reminiscent observations of a spherically expanding shell, even though the evidence is scant. Collectively, the maser observations described here provide additional support for the evolved star hypothesis for IRAS 19312+1950.

We present a concept, developed at the National Astronomy and Ionosphere Center (NAIC) at Arecibo, Puerto Rico, for active suppression of Global Positioning System (GPS) signals in the 305 m dish radio receiver path prior to backend processing. The subsystem does not require an auxiliary antenna and is intended for easy integration with radio telescope systems with a goal of being a turnkey addition to virtually any facility. Working with actual sampled signal data, we have focused on the detection and cancellation of the GPS L3 signal at 1381.05 MHz which, during periodic test modes and particularly during system-wide tests, interfere with observations of objects in a range of redshifts that includes the Coma supercluster, for example. This signal can dynamically change modulation modes and our scheme is capable of detecting these changes and applying cancellation or sending a blanking signal, as appropriate. The subsystem can also be adapted to GPS L1 (1575.42 MHz), L2C (1227.6 MHz), and others. A follow-up is underway to develop a prototype to deploy and evaluate at NAIC.

We report the results of a search for 12.2-GHz methanol maser emission, targeted towards 113 known 6.7-GHz methanol masers associated with 1.2-mm dust continuum emission. Observations were carried out with the Australia Telescope National Facility (ATNF) Parkes 64-m radio telescope in the period 2008 June 20 - 25. We detect 68 12.2-GHz methanol masers with flux densities in excess of our 5-sigma detection limit of 0.55 Jy, 30 of which are new discoveries. This equates to a detection rate of 60 per cent, similar to previous searches of comparable sensitivity. We have made a statistical investigation of the properties of the 1.2-mm dust clumps with and without associated 6.7-GHz methanol maser and find that 6.7-GHz methanol masers are associated with 1.2-mm dust clumps with high flux densities, masses and radii. We additionally find that 6.7-GHz methanol masers with higher peak luminosities are associated with less dense 1.2-mm dust clumps than those 6.7-GHz methanol masers with lower luminosities. We suggest that this indicates that more luminous 6.7-GHz methanol masers are generally associated with a later evolutionary phase of massive star formation than less luminous 6.7-GHz methanol maser sources. Analysis of the 6.7-GHz associated 1.2-mm dust clumps with and without associated 12.2-GHz methanol maser emission shows that clumps associated with both class II methanol maser transitions are less dense than those with no associated 12.2-GHz methanol maser emission; suggesting that 12.2-GHz methanol masers are associated with a later evolutionary phase of massive star formation. We present an evolutionary sequence for masers in high-mass star formation regions, placing quantitative estimates on the relative lifetimes for the first time.

Observations show that 1612 MHz masers of OH/IR stars can fade on a timescale of a decade. This fading is probably associated with the switch from rapid mass loss, which is ultimately linked with an internal He-shell flash, to the much slower mass loss supported by more quiescent conditions. We study the observed maser decay with a composite computational model, comprising a time-dependent chemical model of the envelope, and a radiation transfer model which provides the maser pumping. Our combined model is able to reproduce the rapid decay of maser intensity, following a sudden drop in the stellar mass-loss rate. The explanation for the rapid fall in maser emission is not a fall in the OH number density, or the kinetic temperature in the inverted layers, but the loss of a radiative pump route which carries population from level 1 to level 4 via levels 16, 15 and 11. The loss of these pump routes is a result, in turn, of a greatly reduced energy density of 53 micron radiation.

An isolated, spherically-symmetric, self-gravitating, collisionless system is always a polytrope when it reaches equilibrium (Nakamura 2000). This strongly suggests as a corollary, however, that the same polytrope dominates its precursor states, since the dynamical equations for its constituents can be time-reversed. Moreover this assumption, which precludes a polytrope from ever accreting 100% of the mass from an infalling shell, as a subsequent state will eventually be a polytrope, is confirmed now by our finding that a collisionless polytrope cannot accrete 100% of an infalling shell while simultaneously guaranteeing that the entropy of the Universe as a whole increases. These strictures are only evaded by the shedding of some mass. A polytrope must lose mass to gain mass. We deduce from the time-reversible property of a collisionless polytrope that the scalar sum, P, over constituent momenta in its rest frame is an independent state variable that is conserved with respect to its surface radius through interactions between a polytrope and an infalling shell. This new constraint, together with conservation of energy, enables us (i) to show that an isolated polytrope is indeed stable against spherically-symmetric mass-loss, which is the essential content of our initial assumption; (ii) to calculate both the velocity and the fraction of infall-mass returned to infinity, provided the "accretion law" between the change in mass and surface radius is specified. Numerical results confirm a frequent empirical finding (Livio 2000) that the velocity of a mass outflow is of the same order of magnitude as the escape velocity from the system.

Gravitational contraction always generates a radially directed momentum flux. A particularly simple example occurs in the electron-degenerate cores of AGB stars, which contract steadily under the addition of helium ashes from shell hydrogen burning. The resulting momentum flux is quantified here. And since the cores of AGB stars lack efficient momentum cancellation mechanisms, they can maintain equilibrium by exporting their excess momentum flux to the stellar envelope, which disposes of much of it in a low velocity wind. Gravitational contraction easily accounts for the momentum flux in the solar wind, as well as the flux required to lift mass into the dust formation zone of every AGB star, whereon radiation pressure continues its ejection as a low velocity wind. This mechanism explains the dependence of the AGB mass-loss rate on core mass; its generalization to objects with angular momentum and/or strong magnetic fields suggests a novel explanation of why most planetary nebulae and proto planetary nebulae exhibit axial symmetry. Quasistatic contraction is inherently biased to the generation of the maximum possible momentum flux. Its formalism is therefore readily adapted to providing an upper limit to the momentum flux needed to sustain mass loss when this begins from a semi-continuous rather than impulsive process.

There is an approximately 60% detection rate for 1612 MHz masers in association with red, color-selected IRAS sources, though few are detected from the bluer circumstellar shells of Mira variables. On the other hand and complementarily, past, pre-IRAS 22 GHz surveys detected many water masers in association with Mira variables. This paper reports on a 22 GHz survey of blue, color-selected Miras at Haystack, wherein 18 new detections are found from 238 searched objects.