Indirect dark matter (DM) detection typically involves the observation of standard model (SM) particles emerging from DM annihilation/decay inside regions of high dark matter concentration. We consider an annihilation scenario in which this reaction has to be initiated by one of the DMs involved being boosted while the other is an ambient non-relativistic particle. This "trigger" DM must be created, for example, in a previous annihilation or decay of a heavier component of DM. Remarkably, boosted DM annihilating into gamma-rays at a specific point in a galaxy could actually have traveled from its source at another point in the same galaxy or even from another galaxy. Such a "non-local" behavior leads to a non-trivial dependence of the resulting photon signal on the galactic halo parameters, such as DM density and core size, encoded in the so-called "astrophysical" $J$-factor. These non-local $J$-factors are strikingly different than the usual scenario. A distinctive aspect of this model is that the signal from dwarf galaxies relative to the Milky Way tends to be suppressed from the typical value to various degrees depending on their characteristics. This feature can thus potentially alleviate the mild tension between the DM annihilation explanation of the observed excess of $\sim$ GeV photons from the Milky Way's galactic center vs. the apparent non-observation of the corresponding signal from dwarf galaxies.

It has been suggested that late-universe dark matter decays can alleviate the tension between measurements of $H_0$ in the local universe and its value inferred from cosmic microwave background fluctuations. Decaying dark matter can potentially account for this discrepancy as it reshuffles the energy density between matter and radiation and as a result allows dark energy to become dominant at earlier times. We show that the low multipoles amplitude of the cosmic microwave background anisotropy power spectrum severely constrains the feasibility of late-time decays as a solution to the $H_0$ tension.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Recent analyses of 21-cm neutral hydrogen (HI) emission have demonstrated that HI gas is organized into linear filamentary structures that are preferentially aligned with the local magnetic field, and that the coherence of these structures in velocity space traces line-of-sight magnetic field tangling. On this basis, we introduce a paradigm for modeling the properties of magnetized, dusty regions of the interstellar medium, using the orientation of HI structure at different velocities to map "magnetically coherent" regions of space. We construct three-dimensional (position-position-velocity) Stokes parameter maps using HI4PI full-sky spectroscopic HI data. We compare these maps, integrated over the velocity dimension, to Planck maps of the polarized dust emission at 353 GHz. Without any free parameters governing the relation between HI intensity and dust emission, we find that our $Q$ and $U$ maps are highly correlated ($r > 0.75$) with the 353 GHz $Q$ and $U$ maps of polarized dust emission observed by Planck and reproduce many of its large-scale features. The $E/B$ ratio of the dust emission maps agrees well with the HI-derived maps at large angular scales ($\ell \lesssim 120$), supporting the interpretation that this asymmetry arises from the coupling of linear density structures to the Galactic magnetic field. We demonstrate that our 3D Stokes parameter maps constrain the 3D structure of the Galactic interstellar medium and the orientation of the interstellar magnetic field.

The spatial distribution of neutral hydrogen (H I) emission is a powerful probe of interstellar medium physics. The small-scale structure in H I channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H I observations and 50,985 quasar spectra from SDSS. We measure the equivalent widths of interstellar Na I D$_1$ and Na I D$_2$ absorption, and robustly conclude that they depend more strongly on the column density of small-scale structure in H I than on either the large-scale H I structure or the total H I column. This is inconsistent with the hypothesis that small-scale channel map structure is driven by velocity crowding. Instead, the data favor the interpretation that this emission structure predominantly originates in cold, dense interstellar material, consistent with a clumpy cold neutral medium.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Sharpless 2-27 (Sh2-27) is a nearby ${\rm H\,{\small II}}$ region excited by $\zeta$Oph. We present observations of polarized radio emission from 300 to 480$\,$MHz towards Sh2-27, made with the Parkes 64$\,$m Radio Telescope as part of the Global Magneto-Ionic Medium Survey. These observations have an angular resolution of $1.35^{\circ}$, and the data are uniquely sensitive to magneto-ionic structure on large angular scales. We demonstrate that background polarized emission towards Sh2-27 is totally depolarized in our observations, allowing us to investigate the foreground. We analyse the results of Faraday tomography, mapping the magnetised interstellar medium along the 165$\,$pc path to Sh2-27. The Faraday dispersion function in this direction has peaks at three Faraday depths. We consider both Faraday thick and thin models for this observation, finding that the thin model is preferred. We further model this as Faraday rotation of diffuse synchrotron emission in the Local Bubble and in two foreground neutral clouds. The Local Bubble extends for 80$\,$pc in this direction, and we find a Faraday depth of $-0.8 \pm 0.4\,$rad$\,$m$^{-2}$. This indicates a field directed away from the Sun with a strength of $-2.5\pm1.2\,\mu$G. The near and far neutral clouds are each about 30$\,$pc thick, and we find Faraday depths of $-6.6\pm0.6\,$rad$\,$m$^{-2}$ and $+13.7\pm0.8\,$rad$\,$m$^{-2}$, respectively. We estimate that the line-of-sight magnetic strengths in the near and far cloud are $B_{\parallel, \text{near}} \approx -15\,\mu\text{G}$ and $B_{\parallel, \text{far}} \approx +30\,\mu\text{G}$. Our results demonstrate that Faraday tomography can be used to investigate the magneto-ionic properties of foreground features in front of nearby ${\rm H\,{\small II}}$ regions.

We investigate the physical properties of structures seen in channel map observations of 21-cm neutral hydrogen (HI) emission. HI intensity maps display prominent linear structures that are well aligned with the ambient magnetic field in the diffuse interstellar medium (ISM). Some literature hold that these structures are "velocity caustics", fluctuations imprinted by the turbulent velocity field, and are not three-dimensional density structures in the ISM. We test this hypothesis by stacking probes of the density field -- broadband far infrared (FIR) emission and the integrated HI column density (NHI) -- at the locations of linear HI intensity structures. We find that the HI intensity features are real density structures and not velocity caustics. We generalize the investigation to all small-scale structure in HI channel maps, and analyze this correlation as a function of velocity channel width, finding no measurable contribution from velocity caustics to the HI channel map emission. Further, we find that small-scale HI channel maps structures have elevated FIR/NHI, implying that this emission originates from a colder, denser phase of the ISM than the surrounding material. The data are consistent with a multi-phase diffuse ISM in which small-scale structures in narrow HI channel maps are preferentially cold neutral medium (CNM) that is anisotropically distributed and aligned with the local magnetic field. The shallow spatial power spectrum (SPS) of narrow HI channels is often attributed to velocity caustics. We conjecture instead that the small-scale structure and narrow linewidths typical of CNM explain the observed relationship between the SPS and channel width.

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.

After a decade of great progress in understanding gas flow into, out of, and through the Milky Way, we are poised to merge observations with simulations to build a comprehensive picture of the multi-scale magnetized interstellar medium (ISM). These insights will also be crucial to four bold initiatives in the 2020s: detecting nanohertz gravitational waves with pulsar timing arrays (PTAs), decoding fast radio bursts (FRBs), cosmic B-mode detection, and imaging the Milky Way's black hole with the Event Horizon Telescope (EHT).

Magnetism is one of the most important forces on the interstellar medium (ISM), anisotropically regulating the structure and star formation that drive galactic evolution. Recent high dynamic range observations of diffuse gas and molecular clouds have revealed new links between interstellar structures and the ambient magnetic field. ISM morphology encodes rich physical information, but deciphering it requires high-resolution measurements of the magnetic field: linear polarization of starlight and dust emission, and Zeeman splitting. These measure different components of the magnetic field, and crucially, Zeeman splitting is the only way to directly measure the field strength in the ISM. We advocate a statistically meaningful survey of magnetic field strengths using the 21-cm line in absorption, as well as an observational test of the link between structure formation and field strength using the 21-cm line in emission. Finally, we report on the serendipitous discovery of linear polarization of the 21-cm line, which demands both theoretical and observational follow-up.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

For many dark matter models, the annihilation cross section to two-body final states is difficult to probe with current experiments because the dominant annihilation channel is velocity or helicity suppressed. The inclusion of gauge boson radiation for three-body final states can lift the helicity suppression, allowing a velocity-independent cross section to dominate the annihilation process, and providing an avenue to constrain these models. Here we examine experimental constraints on dark matter that annihilates to two leptons plus a bremsstrahlung boson, $\bar{\ell}+\ell+\gamma/W/Z$. We consider experimental constraints on photon final states from Fermi-LAT using both diffuse photon data and data from dwarf spheroidal galaxies, and compare to the implied constraints from 21 cm measurements. Diffuse photon line searches are generally the strongest over the entire mass regime. We in particular highlight the model in which dark matter annihilates to two neutrinos and a photon, and show that these models are more strongly constrained through photon measurements than through existing neutrino bounds.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

The Galactic neutral hydrogen (HI) sky at high Galactic latitudes is suffused with linear structure. Particularly prominent in narrow spectral intervals, these linear HI features are well aligned with the plane-of-sky magnetic field orientation as measured with optical starlight polarization and polarized thermal dust emission. We analyze the coherence of the orientation of these features with respect to line-of-sight velocity, and propose a new metric to quantify this HI coherence. We show that HI coherence is linearly correlated with the polarization fraction of 353 GHz dust emission. HI coherence constitutes a novel method for measuring the degree of magnetic field tangling along the line of sight in the diffuse interstellar medium. We propose applications of this property for HI-based models of the polarized dust emission in diffuse regions, and for studies of frequency decorrelation in the polarized dust foreground to the cosmic microwave background (CMB).

Recently the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) reported the detection of a 21cm absorption signal stronger than astrophysical expectations. In this paper we study the impact of radiation from dark matter (DM) decay and primordial black holes (PBH) on the 21cm radiation temperature in the reionization epoch, and impose a constraint on the decaying dark matter and PBH energy injection in the intergalactic medium, which can heat up neutral hydrogen gas and weaken the 21cm absorption signal. We consider decay channels DM$\rightarrow e^+e^-, \gamma\gamma$, $\mu^+\mu^-$, $b\bar{b}$ and the $10^{15-17}$g mass range for primordial black holes, and require the heating of the neutral hydrogen does not negate the 21cm absorption signal. For $e^+e^-$, $\gamma\gamma$ final states and PBH cases we find strong 21cm bounds that can be more stringent than the current extragalactic diffuse photon bounds. For the DM$\rightarrow e^+e^-$ channel, the lifetime bound is $\tau_{\rm DM}> 10^{27}$s for sub-GeV dark matter. The bound is $\tau_{\rm DM}\ge 10^{26}$s for sub-GeV DM$\rightarrow \gamma\gamma$ channel and reaches $10^{27}$s at MeV DM mass. For $b\bar{b}$ and $\mu^+\mu^-$ cases, the 21 cm constraint is better than all the existing constraints for $m_{\rm DM}<20$ GeV where the bound on $\tau_{\rm DM}\ge10^{26}$s. For both DM decay and primordial black hole cases, the 21cm bounds significantly improve over the CMB damping limits from Planck data.

We consider dark matter annihilation into a general set of final states of Standard Model particles, including two-body and four-body final states that result from the decay of intermediate states. For dark matter masses ~10-10^5 GeV, we use updated data from Planck and from high gamma-ray experiments such as Fermi-LAT, MAGIC, and VERITAS to constrain the annihilation cross section for each final state. The Planck constraints are the most stringent over the entire mass range for annihilation into light leptons, and the Fermi-LAT constraints are the most stringent for four-body final states up to masses ~10^4 GeV. We consider these constraints in light of the recent AMS antiproton results, and show that for light mediators it is possible to explain the AMS data with dark matter, and remain consistent with Fermi-LAT Inner Galaxy measurements, for m_\chi ~ 60-100 GeV mass dark matter and mediator masses m_\phi / m_\chi ~< 1.

We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor-Couette flow. This is a multiscale perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor-Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg-Landau equation (GLE), while the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.

The magnetorotational instability (MRI) is a fundamental process of accretion disk physics, but its saturation mechanism remains poorly understood despite considerable theoretical and computational effort. We present a multiple scales analysis of the non-ideal MRI in the weakly nonlinear regime -- that is, when the most unstable MRI mode has a growth rate asymptotically approaching zero from above. Here, we develop our theory in a local, Cartesian channel. Our results confirm the finding by Umurhan et al. (2007) that the perturbation amplitude follows a Ginzburg-Landau equation. We further find that the Ginzburg-Landau equation will arise for the local MRI system with shear-periodic boundary conditions when the effects of ambipolar diffusion are considered. A detailed force balance for the saturated azimuthal velocity and vertical magnetic field demonstrates that even when diffusive effects are important, the bulk flow saturates via the combined processes of reducing the background shear and rearranging and strengthening the background vertical magnetic field. We directly simulate the Ginzburg-Landau amplitude evolution for our system and demonstrate the pattern formation our model predicts on long length and time scales. We compare the weakly nonlinear theory results to a direct numerical simulation of the MRI in a thin-gap Taylor Couette flow.

We investigate constraints on the abundance of primordial black holes (PBHs) in the mass range 10^{15}-10^{17} g using data from the Cosmic Microwave Background (CMB) and MeV extragalactic gamma-ray background (EGB). Hawking radiation from PBHs with lifetime greater than the age of the universe leaves an imprint on the CMB through modification of the ionization history and the damping of CMB anisotropies. Using a model for redshift dependent energy injection efficiencies, we show that a combination of temperature and polarization data from Planck provides the strongest constraint on the abundance of PBHs for masses \sim 10^{15}-10^{16} g, while the EGB dominates for masses \gtrsim 10^{16} g. Both the CMB and EGB now rule out PBHs as the dominant component of dark matter for masses \sim 10^{16}-10^{17} g. Planned MeV gamma-ray observatories are ideal for further improving constraints on PBHs in this mass range.

We present a proof-of-concept study of a method to estimate the inclination angle of compact high velocity clouds (CHVCs), i.e. the angle between a CHVC's trajectory and the line-of-sight. The inclination angle is derived from the CHVC's morphology and kinematics. We calibrate the method with numerical simulations, and we apply it to a sample of CHVCs drawn from HIPASS. Implications for CHVC distances are discussed.

The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {{\rm cm}}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud.

Using high-resolution data from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey, we show that linear structure in Galactic neutral hydrogen (HI) correlates with the magnetic field orientation implied by Planck 353 GHz polarized dust emission. The structure of the neutral interstellar medium is more tightly coupled to the magnetic field than previously known. At high Galactic latitudes, where the Planck data are noise-dominated, the HI data provide an independent constraint on the Galactic magnetic field orientation, and hence the local dust polarization angle. We detect strong cross-correlations between template maps constructed from estimates of dust intensity combined with either HI-derived angles, starlight polarization angles, or Planck 353 GHz angles. The HI data thus provide a new tool in the search for inflationary gravitational wave B-mode polarization in the cosmic microwave background, which is currently limited by dust foreground contamination.

As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $\lambda$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds.

We present observations of a new group of structures in the diffuse Galactic ISM: slender, linear HI features we dub "fibers" that extend for many degrees at high Galactic latitude. To characterize and measure the extent and strength of these fibers, we present the Rolling Hough Transform (RHT), a new machine vision method for parameterizing the coherent linearity of structures in the image plane. With this powerful new tool we show the fibers are oriented along the interstellar magnetic field as probed by starlight polarization. We find that these low column density (N(HI) ~ 5 x 10^18 cm^-2) fiber features are most likely a component of the local cavity wall, about 100 pc away. The HI data we use to demonstrate this alignment at high latitude are from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey and the Parkes Galactic All Sky Survey (GASS). We find better alignment in the higher resolution GALFA-HI data, where the fibers are more visually evident. This trend continues in our investigation of magnetically aligned linear features in the Riegel-Crutcher HI cold cloud, detected in the Southern Galactic Plane Survey (SGPS). We propose an application of the RHT for estimating the field strength in such a cloud, based on the Chandrasekhar-Fermi method. We conclude that data-driven, quantitative studies of ISM morphology can be very powerful predictors of underlying physical quantities.

We present 109-115 GHz (3 mm) wide-field spectral line observations of 12^CO, 13^CO and C^18O J=1-0 molecular emission and 5.5 and 8.8 GHz (6 and 3 cm) radio continuum emission towards the high-mass star forming complex known as G305. The morphology of G305 is dominated by a large evacuated cavity at the centre of the complex driven by clusters of O stars surrounded by molecular gas. Our goals are to determine the physical properties of the molecular environment and reveal the relationship between the molecular and ionised gas and star formation in G305. This is in an effort to characterise the star-forming environment and constrain the star formation history in an attempt to evaluate the impact of high-mass stars on the evolution of the G305 complex. Analysis of CO emission in G305 reveals 156 molecular clumps with the following physical characteristics. The 5.5 and 8.8GHz radio continuum emission reveals an extended low surface brightness ionised environment within which we identify 15 large-scale features with a further eight smaller sources projected within these features. By comparing to mid infrared emission and archival data, we identify nine HII regions, seven compact HII regions, one UC HII region, four extended regions. The total integrated flux of the radio continuum emission at 5.5 GHz is ~180 Jy corresponding to a Lyman continuum output of 2.4x10^50 photons s^-1. We compare the ionised and molecular environment with optically identified high-mass stars and ongoing star formation, identified from the literature. Analysis of this dataset reveals a star formation rate of 0.008--0.016 and efficiency of 7--12%, allows us to probe the star formation history of the region and discuss the impact of high-mass stars on the evolution of G305.

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

Observations from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey of the tail of Complex C are presented and the halo clouds associated with this complex cataloged. The properties of the Complex C clouds are compared to clouds cataloged at the tail of the Magellanic Stream to provide insight into the origin and destruction mechanism of Complex C. Magellanic Stream and Complex C clouds show similarities in their mass distributions (slope = -0.7 and -0.6, respectively) and have a common linewidth of 20 - 30 km/s (indicative of a warm component), which may indicate a common origin and/or physical process breaking down the clouds. The clouds cataloged at the tail of Complex C extend over a mass range of 10^1.1 to 10^4.8 solar masses, sizes of 10^1.2 to 10^2.6 pc, and have a median volume density of 0.065 cm^(-3) and median pressure of (P/k) = 580 K cm^{-3}. We do not see a prominent two-phase structure in Complex C, possibly due to its low metallicity and inefficient cooling compared to other halo clouds. From assuming the Complex C clouds are in pressure equilibrium with a hot halo medium, we find a median halo density of 5.8 x 10^(-4) cm^(-3), which given a constant distance of 10 kpc, is at a z-height of ~3 kpc. Using the same argument for the Stream results in a median halo density of 8.4 x 10^(-5) x (60kpc/d) cm^(-3). These densities are consistent with previous observational constraints and cosmological simulations. We also assess the derived cloud and halo properties with three dimensional grid simulations of halo HI clouds and find the temperature is generally consistent within a factor of 1.5 and the volume densities, pressures and halo densities are consistent within a factor of 3.

The compact association Cygnus OB2 is known to contain a large population of massive stars, but its total mass is currently a matter of debate. While recent surveys have uncovered large numbers of OB stars in the area around Cyg OB2, detailed study of the optically brightest among them suggests that most are not part of the association. We observed an additional sample of optically faint OB star candidates, with the aim of checking if more obscured candidates are correspondingly more likely to be members of Cyg OB2. Low resolution spectra of 9 objects allow the rejection of one foreground star and the selection of four O-type stars, which were later observed at higher resolution. In a subsequent run, we observed three more stars in the classification region and three other stars in the far red. We identify five (perhaps six) new evolved very massive stars and three main sequence O-type stars, all of which are likely to be members of Cyg OB2. The new findings allow a much better definition of the upper HR diagram, suggesting an age ~2.5Myr for the association and hinting that the O3-5 supergiants in the association are blue stragglers, either younger or following a different evolutionary path from other cluster members. Though the bulk of the early stars seems to belong to an (approximately) single-age population, there is ample evidence for the presence of somewhat older stars at the same distance. Our results suggest that, even though Cyg OB2 is unlikely to contain as many as 100 O-type stars, it is indeed substantially more massive than was thought prior to recent infrared surveys.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

A tensor description of perturbative Einsteinian gravity about an arbitrary background spacetime is developed. By analogy with the covariant laws of electromagnetism in spacetime, gravito-electromagnetic potentials and fields are defined to emulate electromagnetic gauge transformations under substitutions belonging to the gauge symmetry group of perturbative gravitation. These definitions have the advantage that on a flat background, with the aid of a covariantly constant timelike vector field, a subset of the linearised gravitational field equations can be written in a form that is fully analogous to Maxwell's equations (without awkward factors of 4 and extraneous tensor fields). It is shown how the remaining equations in the perturbed gravitational system restrict the time dependence of solutions to these equations and thereby prohibit the existence of propagating vector fields. The induced gravito-electromagnetic Lorentz force on a test particle is evaluated in terms of these fields together with the torque on a small gyroscope. It is concluded that the analogy of perturbative gravity to Maxwell's description of electromagnetism can be valuable for (quasi-)stationary gravitational phenomena but that the analogy has its limitations.

We report results from HST/NICMOS 2 micron imaging polarimetry of the central region of Centaurus A. In the vicinity of the nucleus we observe a complex polarization structure which we explain by a combination of scattering of nuclear light and dichroic polarization associated with the dust lane. The scattered nuclear radiation is found in an angular region which extends over ~ 70 degrees and thus it does not originate from a highly collimated beam, but is associated with more omni-directional nuclear illumination. These observations also show the presence of an unresolved, highly polarized (P = 11.1 %) nuclear source whose polarization angle PA = 148.2 degrees is perpendicular to the jet axis. We set an upper limit of 0.04'' (~0.8 pc) to its extent. The observed nuclear polarization is naturally accounted for if we are observing scattered light from an otherwise obscured nucleus provided that both the scattering region and the occulting torus are extremely compact, with an outer radius of less than ~ 1 pc. Alternatively, we might be directly seeing the infrared counterpart of the radio core, similar to those found in other low luminosity radio-galaxies observed with HST. We discuss these results in the framework of the FRI / BL Lac unifying model.