Wind is the process that connects Mars' climate system. Measurements of Mars atmospheric winds from orbit would dramatically advance our understanding of Mars and help prepare for human exploration of the Red Planet. Multiple instrument candidates are in development and will be ready for flight in the next decade. We urge the Decadal Survey to make these measurements a priority for 2023-2032.

We present a model for a Noise Matched Phased Array Feed (PAF) system and compare model predictions with the measurement results. The PAF system consists of an array feed, a receiver, a beamformer and a parabolic reflector. The novel aspect of our model is the characterization of the {\em PAF system} by a single matrix. This characteristic matrix is constructed from the open-circuit voltage covariance at the output of the PAF due to signal from the observing source, ground spillover noise, sky background noise and (low-noise) amplifier (LNA) noise. The best signal-to-noise ratio on the source achievable with the PAF system will be the maximum eigenvalue of the characteristic matrix. The voltage covariance due to signal and spillover noise are derived by applying the Lorentz reciprocity theorem. The receiver noise covariance and noise temperature are obtained in terms of Lange invariants such that they are suitable for noise matching the array feed. The model predictions are compared with the measured performance of a 1.4 GHz, 19-element, dual-polarized PAF on the Robert C. Byrd Green Bank Telescope. We show that the model predictions, obtained with an additional noise contribution due to the measured losses ahead of the low-noise amplifier, compare well with the measured ratio of system temperature to aperture efficiency as a function of frequency and as a function of offset from the boresight. Further, our modeling indicates that the bandwidth over which this ratio is optimum can be improved by a factor of at least two by noise matching the PAF with the LNA.

In this report we present a model for phased array feed (PAF) and compare the model predictions with measurements. A theory for loss-less PAF is presented first. To develop the theory we ask the question -- what is the best $T_{sys}/\eta_{ap}$ that can be achieved when a PAF is used on a telescope to observe a source at an angle $\theta_s, \phi_s$ from the boresight direction ? We show that a characteristic matrix for the {\em system} (i.e. PAF+telescope+receiver) can be constructed starting from the signal-to-noise ratio of the observations and the best $T_{sys}/\eta_{ap}$ can be obtained from the maximum eigenvalue of the characteristic matrix. For constructing the characteristic matrix, we derive the open-circuit voltage at the output of the antenna elements in the PAF due to (a) radiation from source, (b) radiation from ground (spillover), (c) radiation from sky background and (d) noise due to the receiver. The characteristic matrix is then obtained from the correlation matrices of these voltages. We then describe a modeling program developed to implement the theory presented here. Finally the model predictions are compared with results from test observations made toward Virgo A with a prototype PAF (Kite array) on the GBT (Roshi et al. 2015).

A new 1.4 GHz 19-element, dual-polarization, cryogenic phased array feed (PAF) radio astronomy receiver has been developed for the Robert C. Byrd Green Bank Telescope (GBT) as part of FLAG (Focal L-band Array for the GBT) project. Commissioning observations of calibrator radio sources show that this receiver has the lowest reported beamformed system temperature ($T_{\rm sys}$) normalized by aperture efficiency ($\eta$) of any phased array receiver to date. The measured $T_{\rm sys}/\eta$ is $25.4 \pm 2.5$ K near 1350 MHz for the boresight beam, which is comparable to the performance of the current 1.4 GHz cryogenic single feed receiver on the GBT. The degradation in $T_{\rm sys}/\eta$ at $\sim$ 4 arcmin (required for Nyquist sampling) and $\sim$ 8 arcmin offsets from the boresight is, respectively, $\sim$ 1\% and $\sim$ 20\% of the boresight value. The survey speed of the PAF with seven formed beams is larger by a factor between 2.1 and 7 compared to a single beam system depending on the observing application. The measured performance, both in frequency and offset from boresight, qualitatively agree with predictions from a rigorous electromagnetic model of the PAF. The astronomical utility of the receiver is demonstrated by observations of the pulsar B0329+54 and an extended HII region, the Rosette Nebula. The enhanced survey speed with the new PAF receiver will enable the GBT to carry out exciting new science, such as more efficient observations of diffuse, extended neutral hydrogen emission from galactic in-flows and searches for Fast Radio Bursts.

In this paper, we describe the development of a new L-band (1.4 GHz) Cryogenic Phased Array Feed (PAF) system, referred to as the GBT2 array. Results from initial measurements made with the GBT2 array are also presented. The PAF was developed for the Green Bank Telescope (GBT) as part of the Focal L-band Array for the GBT (FLAG) project. During the first stage of the development work (Phase I), a prototype cryogenic 19 element dual-polarized array with "Kite" dipole elements was developed and tested on the GBT. The measured system temperature over efficiency ($T_{sys}/\eta$) ratio for the bore sight beam of the Kite array was 45.5 K at 1.55 GHz. The off-boresight $T_{sys}/\eta$ shows an increase by 13 K at an offset equal to the half power beam width (7$^{'}$.2 at 1.7 GHz). Our measurements indicate that the off-boresight degradation and field-of-view (FoV) limitation of the Kite array is simply due to the fixed array size. To increase the FoV, a new 19-element GBT2 array with larger array spacing was developed during FLAG Phase II. The frequency response of the array was optimized from 1.2 to 1.6 GHz. A system with larger cryostat, new low noise amplifiers (LNAs), down-conversion and digitization close to the front end, unformatted digital transmission over fiber, ROACH II based polyphase filter banks (PFBs) with bandwidth 150 MHz and a data acquisition system that records voltage samples from one of the PFB channels were all developed. The data presented here is processed off-line. The receiver temperature measured with the new system is 17 K at 1.4 GHz, an improvement $>$ 8 K over the previous Kite array. Measurements with the GBT2 array on the telescope are in progress. A real time 150 MHz beamformer is also being developed as part of an NSF-funded collaboration between NRAO/GBO/BYU \& West Virginia University (Beamformer Project) to support science observations.

We use the Millennium simulation to show that halo clustering varies significantly with cosmic web type. Halos are classified as node, filament, sheet and void halos based on the eigenvalue decomposition of the velocity shear tensor. The velocity field is sampled by the peculiar velocities of a fixed number of neighbouring halos and spatial derivatives are computed using a kernel borrowed from smoothed particle hydrodynamics. The classification scheme is used to examine the clustering of halos as a function of web type for halos with masses larger than $10^{11}$. We find that node halos show positive bias, filament halos show negligible bias, and void and sheet halos are anti-biased independent of halo mass. Our findings suggest that the mass dependence of halo clustering is rooted in the composition of web types as a function of halo mass. The substantial fraction of node type halos for halo masses $\gtrsim 2\times10^{13}\,h^{-1}\rm M_\odot$ leads to positive bias. Filament type halos prevail at intermediate masses, $10^{12} - 10^{13}\,h^{-1}\rm M_\odot$, resulting in unbiased clustering. The large contribution of sheet type halos at low halo masses $\lesssim 10^{12}\,h^{-1}\rm M_\odot$ generates anti-biasing.

The cosmic web defines the large scale distribution of matter we see in the Universe today. Classifying the cosmic web into voids, sheets, filaments and nodes allows one to explore structure formation and the role environmental factors have on halo and galaxy properties. While existing studies of cosmic web classification concentrate on grid based methods, this work explores a Lagrangian approach where the V-web algorithm proposed by Hoffman et al. (2012) is implemented with techniques borrowed from smoothed particle hydrodynamics. The Lagrangian approach allows one to classify individual objects (e.g. particles or halos) based on properties of their nearest neighbours in an adaptive manner. It can be applied directly to a halo sample which dramatically reduces computational cost and potentially allows an application of this classification scheme to observed galaxy samples. Finally, the Lagrangian nature admits a straight forward inclusion of the Hubble flow negating the necessity of a visually defined threshold value which is commonly employed by grid based classification methods.

We have used 2D back-projection Doppler tomography as a tool to examine the influence of gravitational and magnetic phenomena in interacting binaries which undergo mass transfer from a magnetically-active star onto a non-magnetic main sequence star. This multi-tiered study of over 1300 time-resolved spectra of 13 Algol binaries involved calculations of the predicted dynamical behavior of the gravitational flow and the dynamics at the impact site, analysis of the velocity images constructed from tomography, and the influence on the tomograms of orbital inclination, systemic velocity, orbital coverage, and shadowing. The H$\alpha$ tomograms revealed eight sources: chromospheric emission, a gas stream along the gravitational trajectory, a star-stream impact region, a bulge of absorption or emission around the mass-gaining star, a Keplerian accretion disk, an absorption zone associated with hotter gas, a disk-stream impact region, and a hot spot where the stream strikes the edge of a disk. We described several methods used to extract the physical properties of the emission sources directly from the velocity images, including S-wave analysis, the creation of simulated velocity tomograms from hydrodynamic simulations, and the use of synthetic spectra with tomography to sequentially extract the separate sources of emission from the velocity image. In summary, the tomography images have revealed results that cannot be explained solely by gravitational effects: chromospheric emission moving with the mass-losing star, a gas stream deflected from the gravitational trajectory, and alternating behavior between stream state and disk state. Our results demonstrate that magnetic effects cannot be ignored in these interacting binaries.

Cosmicflows-2 is a compilation of distances and peculiar velocities for over 8000 galaxies. Numerically the largest contributions come from the luminosity-linewidth correlation for spirals, the TFR, and the related Fundamental Plane relation for E/S0 systems, but over 1000 distances are contributed by methods that provide more accurate individual distances: Cepheid, Tip of the Red Giant Branch, Surface Brightness Fluctuation, SNIa, and several miscellaneous but accurate procedures. Our collaboration is making important contributions to two of these inputs: Tip of the Red Giant Branch and TFR. A large body of new distance material is presented. In addition, an effort is made to assure that all the contributions, our own and those from the literature, are on the same scale. Overall, the distances are found to be compatible with a Hubble Constant H_0 = 74.4 +-3.0 km/s/Mpc. The great interest going forward with this data set will be with velocity field studies. Cosmicflows-2 is characterized by a great density and high accuracy of distance measures locally, falling to sparse and coarse sampling extending to z=0.1.

We report on the development of a prototype integrated receiver front-end that combines all conversions from RF to baseband, from analog to digital, and from copper to fiber into one compact assembly, with the necessary gain and stability suitable for radio astronomy applications. The emphasis in this article is on a novel digital data link over optical fiber which requires no formatting in the front-end, greatly reducing the complexity, bulk, and power consumption of digital electronics inside the antenna, facilitating its integration with the analog components, and minimizing the self-generated radio-frequency interference (RFI) which could leak into the signal path. Management of the serial data link is performed entirely in the back-end based on the statistical properties of signals with a strong random noise component. In this way, the full benefits of precision and stability afforded by conventional digital data transmission are realized with far less overhead at the focal plane of a radio telescope.

The neutral hydrogen properties of 1,822 galaxies are being studied with the Green Bank 100m and the Parkes 64m telescopes as part of the 'Cosmic Flows' program. Observed parameters include systemic velocities, profile line widths, and fluxes. The line width information can be combined with optical and infrared photometry to obtain distances. The 1,822 HI observations complement an inventory of archives. All told, HI line width information is available for almost all of five samples: (i) luminosity-line width correlation calibrators, (ii) zero-point calibrators for the supernova Ia scale, (iii) a dense local sample of spiral galaxies with M_{Ks} < -21 within 3,000 km/s, (iv) a sparser sample of 60 $\mu$m selected galaxies within 6,000 km/s that provides all-sky coverage of our extended supercluster complex, and (v) an even sparser sample of flat galaxies, extreme edge-on spirals, extending in a volume out to 12,000 km/s. The HI information for 13,941 galaxies, whether from the archives or acquired as part of the Cosmic Flows observational program, is uniformly re-measured and made available through the Extragalactic Distance Database web site.

In this paper we present a novel class of compact orthomode transducers which use digital calibration to synthesize the desired polarization vectors while maintaining high isolation and minimizing mass and volume. These digital orthomode transducers consist of an arbitrary number of planar probes in a circular waveguide, each of which is connected to an independent receiver chain designed for stability of complex gain. The outputs of each receiver chain are then digitized and combined numerically with calibrated, complex coefficients. Measurements on two prototype digital orthomode transducers, one with three probes and one with four, show better than 50 dB polarization isolation over a 10 C temperature range with a single calibration.

This article reports on the first step in a focused program to re-optimize radio astronomy receiver architecture to better take advantage of the latest advancements in commercial digital technology. Specifically, an L-Band sideband-separating downconverter has been built using a combination of careful (but ultimately very simple) analog design and digital signal processing to achieve wideband downconversion of an RFI-rich frequency spectrum to baseband in a single mixing step, with a fixed-frequency Local Oscillator and stable sideband isolation exceeding 50 dB over a 12 degree C temperature range.

Phased array feeds (PAFs) for reflector antennas offer the potential for increased reflector field of view and faster survey speeds. To address some of the development challenges that remain for scientifically useful PAFs, including calibration and beamforming algorithms, sensitivity optimization, and demonstration of wide field of view imaging, we report experimental results from a 19 element room temperature L-band PAF mounted on the Green Bank 20-Meter Telescope. Formed beams achieved an aperture efficiency of 69% and system noise temperature of 66 K. Radio camera images of several sky regions are presented. We investigate the noise performance and sensitivity of the system as a function of elevation angle with statistically optimal beamforming and demonstrate cancelation of radio frequency interference sources with adaptive spatial filtering.

An important component of the Extragalactic Distance Database (EDD) at http://edd.ifa.hawaii.edu is a group of catalogs related to the measurement of HI line profile parameters. One of these is the All Digital HI catalog which contains an amalgam of information from new data and old. The new data results from observations with Arecibo and Parkes telescopes and with the Green Bank Telescope (GBT), including continuing input since the award of the NRAO Cosmic Flows Large Program. The old data has been collected from archives, wherever available, particularly the Cornell University Digital HI Archive, the Nancay Telescope extragalactic HI archive, and the Australia Telescope archive. The catalog currently contains information on ~15,000 profiles relating to ~13,000 galaxies. The channel - flux per channel files, from whatever source, are carried through a common pipeline. The derived parameter of greatest interest is W_m50, the profile width at 50% of the mean flux. After appropriate adjustment, the parameter W_mx is derived, the linewidth which statistically approximates the peak to peak maximum rotation velocity before correction for inclination, 2 V_max sin(i).

Radio astronomy observations in the coming decade will require new levels of sensitivity while mapping large regions of space with much greater efficiency than is achieved with current telescopes. This requires new instrumentation with the greatest achievable sensitivity, dynamic range, and field of view. Receiver noise is quickly approaching fundamental limits at most radio wavelengths, so significant gains in sensitivity can only be made by increasing collecting area. Jointly, these requirements suggest using large arrays of smaller antennas, or many moderate-size antennas equipped with multi-beam arrays. The challenge is to develop receivers and wide bandwidth data transport systems which are lower cost, more compact, more reliable, lower weight, and more reproducible than the best current systems, with no compromise to performance. This can be achieved with a greater degree of component integration, extensive use of digital signal processing and transport, and replacement of functions currently performed in bulky waveguide and coaxial cable components with digital arithmetic and thin optical fibers. In this white paper, we outline the complete redesign and re-optimization of receiver architecture to take advantage of the latest advancements in commercial technology. This involves the seamless integration of the conversions from RF to baseband, from analog to digital, and from copper to fiber within a single receiver module.

This white paper offers cautionary observations about the planning and development of new, large radio astronomy instruments. Complexity is a strong cost driver so every effort should be made to assign differing science requirements to different instruments and probably different sites. The appeal of shared resources is generally not realized in practice and can often be counterproductive. Instrument optimization is much more difficult with longer lists of requirements, and the development process is longer and less efficient. More complex instruments are necessarily further behind the technology state of the art because of longer development times. Including technology R&D in the construction phase of projects is a growing trend that leads to higher risks, cost overruns, schedule delays, and project de-scoping. There are no technology breakthroughs just over the horizon that will suddenly bring down the cost of collecting area. Advances come largely through careful attention to detail in the adoption of new technology provided by industry and the commercial market. Radio astronomy instrumentation has a very bright future, but a vigorous long-term R&D program not tied directly to specific projects needs to be restored, fostered, and preserved.

Accreting X-ray pulsars provide us with laboratories for the study of extreme gravitational and magnetic fields, hence accurate descriptions of their observational properties contribute to our understanding of this group of objects. We aim to detect a cyclotron resonance scattering feature in the Be/X-ray binary Cep X-4 and to investigate pulse profile and spectral changes through the outburst. Spectral fitting and timing analysis are employed to probe the properties of Cep X-4 during an outburst in 2002 June. A previously announced cyclotron feature at 30.7 keV is confirmed, while the source shows spectral behaviour and luminosity related changes similar to those observed in previous outbursts. The long-term X-ray lightcurve shows a periodicity at 20.85 d, which could be attributed to the orbit in this Be system.

We study a sample of spectroscopic binaries (SBs) in the local solar neighbourhood to 100 pc and to an absolute magnitude of 4 in an attempt to find the distributions of the period, primary mass and the mass ratio as well as the IMF of the local population of field binaries. The sample was collated using available SB data and the Hipparcos catalogue, the latter being used for distances and to refer numbers of objects to fractions of the local stellar population as a whole. We use the better-determined double-lined SBs (SB2s) to calibrate a Monte-Carlo approach to modelling the mass ratio distribution of the single-lined SBs (SB1s) from their mass functions and primary masses. While a complete sample is not possible, given the data available, we are able to address important questions of incompleteness and parameter-specific biases by comparing subsamples of SBs with different ranges in parameter space. Our results show a clear peak in the mass ratio distribution of field binaries near unity. This is dominated by the SB2s, but the flat distribution of the SB1s is inconsistent with their components being chosen independently at random from a steep IMF.

The neutral hydrogen 21-cm spectral line (1420.4 MHz) and the four 18-cm lines of the hydroxyl molecule (1612-1720 MHz) are observable at redshifts which put their measured line frequencies well below their protected frequency bands. Part of the redshift ranges (z = 0.171-0.477 for HI and z = 0.37-0.73 for OH) fall in the 960 to 1215 MHz band that is allocated to aircraft navigation. Most of the signals in this band are pulsed emissions of low duty cycle so much of the time between pulses is interference free. This paper outlines the structure and measured properties of signals in this band and demonstrates a signal processing strategy that is effective at removing the pulsed signals from spectra at sensitivities produced by several hours of integration.