The next-generation Very Large Array (ngVLA) is an astronomical observatory planned to operate at centimeter wavelengths (25 to 0.26 centimeters, corresponding to a frequency range extending from 1.2 to 116 GHz). The observatory will be a synthesis radio telescope constituted of approximately 244 reflector antennas each of 18 meters diameter, and 19 reflector antennas each of 6 meters diameter, operating in a phased or interferometric mode. We provide a technical overview of the Reference Design of the ngVLA. This Reference Design forms a baseline for a technical readiness assessment and the construction and operations cost estimate of the ngVLA. The concepts for major system elements such as the antenna, receiving electronics, and central signal processing are presented.

The next-generation Very Large Array (ngVLA) is an astronomical observatory planned to operate at centimeter wavelengths (25 to 0.26 centimeters, corresponding to a frequency range extending from 1.2 GHz to 116 GHz). The observatory will be a synthesis radio telescope constituted of approximately 214 reflector antennas each of 18 meters diameter, operating in a phased or interferometric mode. We provide an overview of the current system design of the ngVLA. The concepts for major system elements such as the antenna, receiving electronics, and central signal processing are presented. We also describe the major development activities that are presently underway to advance the design.

We describe an image-based method that uses two radio criteria, compactness and spectral index, to identify promising pulsar candidates among Fermi Large Area Telescope (LAT) unassociated sources. These criteria are applied to those radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150 MHz (TGSS ADR1) found within the error ellipses of unassociated sources from the 3FGL catalog and a preliminary source list based on 7 years of LAT data. After follow-up interferometric observations to identify extended or variable sources, a list of 16 compact, steep-spectrum candidates is generated. An ongoing search for pulsations in these candidates, in gamma rays and radio, has found six millisecond pulsars and one normal pulsar. A comparison of this method with existing selection criteria based on gamma-ray spectral and variability properties suggests that the pulsar discovery space using Fermi may be larger than previously thought. Radio imaging is a hitherto underutilized source selection method that can be used, as with other multi-wavelength techniques, in the search for Fermi pulsars.

We summarize the design, capabilities, and some of the priority science goals of a next generation Very Large Array (ngVLA). The ngVLA is an interferometric array with 10x larger effective collecting area and 10x higher spatial resolution than the current VLA and the Atacama Large Millimeter Array (ALMA), optimized for operation in the wavelength range 0.3cm to 3cm. The ngVLA opens a new window on the Universe through ultra-sensitive imaging of thermal line and continuum emission down to milliarcecond resolution, as well as unprecedented broad band continuum polarimetric imaging of non-thermal processes. The continuum resolution will reach 9mas at 1cm, with a brightness temperature sensitivity of 6K in 1 hour. For spectral lines, the array at 1" resolution will reach 0.3K surface brightness sensitivity at 1cm and 10 km/s spectral resolution in 1 hour. These capabilities are the only means with which to answer a broad range of critical scientific questions in modern astronomy, including direct imaging of planet formation in the terrestrial-zone, studies of dust-obscured star formation and the cosmic baryon cycle down to pc-scales out to the Virgo cluster, making a cosmic census of the molecular gas which fuels star formation back to first light and cosmic reionization, and novel techniques for exploring temporal phenomena from milliseconds to years. The ngVLA is optimized for observations at wavelengths between the superb performance of ALMA at submm wavelengths, and the future SKA1 at few centimeter and longer wavelengths. This memo introduces the project. The science capabilities are outlined in a parallel series of white papers. We emphasize that this initial set of science goals are simply a starting point for the project. We invite comment on these programs, as well as new ideas, through our public forum link on the ngVLA web page https://science.nrao.edu/futures/ngvla

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

The radio pulse from a pulsar can be temporally broadened by multipath scattering in the interstellar medium and by instrumental effects within the radio telescope. The observed pulse shape is a convolution of the intrinsic one with the impulse responses of the scattering medium and instrumentation. Until recently, common methods used to model the observed shape make assumptions regarding the intrinsic pulse shape and impulse responses, compute the convolution numerically, and solve for the pulse width and scattering timescale iteratively. An analytical solution is shown to exist for the specific case of the temporal broadening of a Gaussian-shaped pulse by a thin scattering screen. The solution is applied to multi-frequency observations of PSR B1834-10 to characterize the frequency dependence of its intrinsic pulse width and scattering timescale.

Radio polarimetry is a three-dimensional statistical problem. The three-dimensional aspect of the problem arises from the Stokes parameters Q, U, and V, which completely describe the polarization of electromagnetic radiation and conceptually define the orientation of a polarization vector in the Poincar'e sphere. The statistical aspect of the problem arises from the random fluctuations in the source-intrinsic polarization and the instrumental noise. A simple model for the polarization of pulsar radio emission has been used to derive the three-dimensional statistics of radio polarimetry. The model is based upon the proposition that the observed polarization is due to the incoherent superposition of two, highly polarized, orthogonal modes. The directional statistics derived from the model follow the Bingham-Mardia and Fisher family of distributions. The model assumptions are supported by the qualitative agreement between the statistics derived from it and those measured with polarization observations of the individual pulses from pulsars. The orthogonal modes are thought to be the natural modes of radio wave propagation in the pulsar magnetosphere. The intensities of the modes become statistically independent when generalized Faraday rotation (GFR) in the magnetosphere causes the difference in their phases to be large. A stochastic version of GFR occurs when fluctuations in the phase difference are also large, and may be responsible for the more complicated polarization patterns observed in pulsar radio emission.

The location of objects on the celestial sphere is a fundamental measurement in astronomy, and the distribution of these objects within the Milky Way is important for understanding their evolution as well as the large scale structure of the Galaxy. Here, physical concepts in Galactic astronomy are illustrated using straightforward mathematics and simplifying assumptions regarding the geometry of the Galaxy. Specifically, an analytical model for a smooth distribution of particles in an oblate ellipsoid is used to replicate the observed distributions of the Galactic coordinates for pulsars and supernova remnants. The distributions and the Lambert equal area projections (LEAPs) of the coordinates suggest that the dominant factors determining the general shape of the distributions are the heavy concentration of objects in the Galactic plane and the offset of the Galactic center from the coordinate system origin. The LEAPs and the distributions also show that the dispersion of pulsars about and along the plane are much larger than that for their progenitor supernovae. Additionally, the model can be used to derive an analytical expression for the dispersion measure along any line of sight within the Galaxy. The expression is used to create a hypothetical dispersion measure-distance map for pulsars in the Galaxy.

In almost 30 years of operation, the Very Large Array (VLA) has proved to be a remarkably flexible and productive radio telescope. However, the basic capabilities of the VLA have changed little since it was designed. A major expansion utilizing modern technology is currently underway to improve the capabilities of the VLA by at least an order of magnitude in both sensitivity and in frequency coverage. The primary elements of the Expanded Very Large Array (EVLA) project include new or upgraded receivers for continuous frequency coverage from 1 to 50 GHz, new local oscillator, intermediate frequency, and wide bandwidth data transmission systems to carry signals with 16 GHz total bandwidth from each antenna, and a new digital correlator with the capability to process this bandwidth with an unprecedented number of frequency channels for an imaging array. Also included are a new monitor and control system and new software that will provide telescope ease of use. Scheduled for completion in 2012, the EVLA will provide the world research community with a flexible, powerful, general-purpose telescope to address current and future astronomical issues.

A variety of intriguing polarization patterns are created when polarization observations of the single pulses from radio pulsars are displayed in a two-dimensional projection of the Poincare sphere. In many pulsars, the projections produce two clusters of data points that reside at antipodal points on the sphere. The clusters are formed by fluctuations in polarization amplitude that are parallel to the unit vectors representing the polarization states of the wave propagation modes in the pulsar magnetosphere. In other pulsars, however, the patterns are more complex, resembling annuli and bow ties or bars. The formation of these complex patterns is not understood and largely unexplored. An empirical model of pulsar polarization is used to show that these patterns arise from polarization fluctuations that are perpendicular to the mode vectors. The model also shows that the modulation index of the polarization amplitude is an indicator of polarization pattern complexity. A stochastic version of generalized Faraday rotation can cause the orientation of the polarization vectors to fluctuate and is a possible candidate for the perpendicular fluctuations incorporated in the model. Alternative models indicate that one mode experiences perpendicular fluctuations and the other does not, suggesting that the fluctuations could also be due to a mode-selective random process, such as scattering in the magnetosphere. A polarization stability analysis of the patterns implies that processes intrinsic to the emission are more effective in depolarizing the emission than fluctuations in the orientation of its polarization vector.

A statistical model of the polarization of pulsar radio emission is used to derive the general statistics of a polarization vector's orientation angles. The theoretical distributions are compared with orientation angle histograms computed from single-pulse, polarization observations of PSR B2020+28. The favorable agreement between the theoretical and measured distributions lends support to the underlying assumptions of the statistical model, and demonstrates, like recent work on other pulsars, that the handedness of circular polarization is associated with the radiation's orthogonally polarized modes. Comprehensive directional statistics of the vector's orientation angles are also derived, and are shown to follow the Watson bipolar and Fisher distributions in its limiting forms.

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical questions. We have recently finished a VLBI astrometry project to determine the proper motions and parallaxes of 27 pulsars, thereby doubling the total number of pulsar parallaxes. Here we summarise our astrometric technique and present the discovery of a pulsar moving in excess of 1000 km/s. As an example of the application of high precision pulsar astrometry we also infer the identification of 2 pulsars originating from a disrupted binary in the Cygnus Superbubble.

The polarization position angles (PA) of pulsar radio emission occupy a distribution that can be much wider than what is expected from the average linear polarization and the off-pulse instrumental noise. Contrary to our limited understanding of the emission mechanism, the excess dispersion in PA implies that pulsar PAs vary in a random fashion. An eigenvalue analysis of the measured Stokes parameters is developed to determine the origin of the excess PA dispersion. The analysis is applied to sensitive, well-calibrated polarization observations of PSR B1929+10 and PSR B2020+28. The analysis clarifies the origin of polarization fluctuations in the emission and reveals that the excess PA dispersion is caused by the isotropic inflation of the data point cluster formed by the measured Stokes parameters. The inflation of the cluster is not consistent with random fluctuations in PA, as might be expected from random changes in the orientation of the magnetic field lines in the emission region or from stochastic Faraday rotation in either the pulsar magnetosphere or the interstellar medium. The inflation of the cluster, and thus the excess PA dispersion, is attributed to randomly polarized radiation in the received pulsar signal. The analysis also indicates that orthogonal polarization modes (OPM) occur where the radio emission is heavily modulated. In fact, OPM may only occur where the modulation index exceeds a critical value of about 0.3.

Polarization observations of the radio emission from PSR B2016+28 at 1404 MHz reveal properties that are consistent with two, very different, interpretations of the pulsar's viewing geometry. The pulsar's average polarization properties show a rapid change in position angle (PA) near the pulse center, suggesting that the observer's sightline nearly intersects the star's magnetic pole. But single pulse, polarization observations of the pulsar show nearly orthogonal modes of polarization following relatively flat and parallel PA trajectories across the pulse, suggesting that the sightline is far from the pole. Additionally, PA histograms reveal a "modal connecting bridge", of unknown origin, joining the modal PA trajectories over much of the pulse and following the rapid PA change shown in the average data. The nonorthogonality of polarization modes is incorporated in a statistical model of radio polarization to account for the deviations from mode orthogonality that are observed in the pulsar. The model is used to interpret the rapid PA change and modal connecting bridge as a longitudinally-resolved transition between modes of nonorthogonal polarization. Thus, the modal PA trajectories are argued to reflect the pulsar's true viewing geometry. This interpretation is consistent with the pulsar's morphological classification, preserves the Radhakrishnan & Cooke model of pulsar radio emission, and avoids the complication that the modal connecting bridge might be produced by some other emission mechanism. The statistical model's ability to simulate the rich variety of polarization properties observed in the emission lends additional support to the model's applicability and its underlying assumption that the polarization modes occur simultaneously.

The measurement of radio polarization may be regarded as a three-dimensional statistical problem because the large photon densities at radio wavelengths allow the simultaneous detection of the three Stokes parameters which completely describe the radiation's polarization. The statistical nature of the problem arises from the fluctuating instrumental noise, and possibly from fluctuations in the radiation's polarization. A statistical model is used to derive the general, three-dimensional statistics that govern radio polarization measurements. The statistics are derived for specific cases of source-intrinsic polarization, with an emphasis on the orthogonal polarization modes in pulsar radio emission. The statistics are similar to those commonly found in other fields of the physical, biological, and Earth sciences. Given the highly variable linear and circular polarization of pulsar radio emission, an understanding of the three-dimensional statistics may be an essential prequisite to a thorough interpretation of pulsar polarization data.

The measurement or constraint of the masses of neutron stars and their binary companions tests theories of neutron star structure and of pulsar formation and evolution. We have measured the rate of the general relativistic advance of the longitude of periastron for the pulsar PSR B1802$-$07: $\dot\omega=0\fdg060\pm0\fdg009\,\mbox{yr}^{-1}$, which implies a total system mass, pulsar plus companion star, of $M=1.7\pm0.4\,\Msun$. We also present a much improved measurement of the rate of periastron advance for PSR B2303+46: $\dot\omega = 0\fdg0099\pm0\fdg0002\,\mbox{yr}^{-1}$, implying $M=2.53\pm0.08\,\Msun$ for this system. We discuss the available constraints on distribution of mass between the pulsars and their companions, and we compare the pulsar masses with other determinations of neutron star masses.