Small-scale magnetic flux ropes (SFRs) are a type of structures in the solar wind that possess helical magnetic field lines. In a recent report (Chen & Hu 2020), we presented the radial variations of the properties of SFR from 0.29 to 8 au using in situ measurements from the Helios, ACE/Wind, Ulysses, and Voyager spacecraft. With the launch of the Parker Solar Probe (PSP), we extend our previous investigation further into the inner heliosphere. We apply a Grad-Shafranov-based algorithm to identify SFRs during the first two PSP encounters. We find that the number of SFRs detected near the Sun is much less than that at larger radial distances, where magnetohydrodynamic (MHD) turbulence may act as the local source to produce these structures. The prevalence of Alfvenic structures significantly suppresses the detection of SFRs at closer distances. We compare the SFR event list with other event identification methods, yielding a dozen well-matched events. The cross-section maps of two selected events confirm the cylindrical magnetic flux rope configuration. The power-law relation between the SFR magnetic field and heliocentric distances seems to hold down to 0.16 au.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins.

First results from the Parker Solar Probe (PSP) mission have revealed ubiquitous coherent ion-scale waves in the inner heliosphere, which are signatures of kinetic wave-particle interactions and fluid-scale instabilities. However, initial studies of the circularly polarized ion-scale waves observed by PSP have only thoroughly analyzed magnetic field signatures, precluding a determination of solar-wind frame propagation direction and intrinsic wave-polarization. A comprehensive determination of wave-properties requires measurements of both electric and magnetic fields. Here, we use full capabilities of the PSP/FIELDS instrument suite to measure both the electric and magnetic components of circularly polarized waves. Comparing spacecraft frame magnetic field measurements with the Doppler-shifted cold-plasma dispersion relation for parallel transverse waves constrains allowable plasma frame polarizations and wave-vectors. We demonstrate that the Doppler-shifted cold-plasma dispersion has a maximum spacecraft frequency $f_{sc}^{*}$ for which intrinsically right-handed fast-magnetosonic waves (FMWs) propagating sunwards can appear left-handed in the spacecraft frame. Observations of left-handed waves with $|f|>f_{sc}^{*}$ are uniquely explained by intrinsically left-handed, ion-cyclotron, waves (ICWs). We demonstrate that electric field measurements for waves with $|f|>f_{sc}^{*}$ are consistent with ICWs propagating away from the sun, verifying the measured electric field. Applying the verified electric field measurements to the full distribution of waves suggests that, in the solar wind frame, the vast majority of waves propagate away from the sun, indicating that the observed population of coherent ion-scale waves contains both intrinsically left and right hand polarized modes.

Parker Solar Probe (PSP), NASA's latest and closest mission to the Sun, is on a journey to investigate fundamental enigmas of the inner heliosphere. This paper reports initial observations made by the Solar Probe Analyzer for Ions (SPAN-I), one of the instruments in the Solar Wind Electrons Alphas and Protons (SWEAP) instrument suite. We address the presence of secondary proton beams in concert with ion-scale waves observed by FIELDS, the electromagnetic fields instrument suite. We show two events from PSP's 2nd orbit that demonstrate signatures consistent with wave-particle interactions. We showcase 3D velocity distribution functions (VDFs) measured by SPAN-I during times of strong wave power at ion-scales. From an initial instability analysis, we infer that the VDFs departed far enough away from local thermodynamic equilibrium (LTE) to provide sufficient free energy to locally generate waves. These events exemplify the types of instabilities that may be present and, as such, may guide future data analysis characterizing and distinguishing between different wave-particle interactions.

The shape of the electron velocity distribution function plays an important role in the dynamics of the solar wind acceleration. Electrons are normally modelled with three components, the core, the halo, and the strahl. We investigate how well the fast strahl electrons in the inner heliosphere preserve the information about the coronal electron temperature at their origin. We analysed the data obtained by two missions, Helios spanning the distances between 65 and 215 R$_S$, and Parker Solar Probe (PSP) reaching down to 35 R$_S$ during its first two orbits around the Sun. The electron strahl was characterised with two parameters, pitch-angle width (PAW), and the strahl parallel temperature (T$_{s\parallel}$). PSP observations confirm the already reported dependence of strahl PAW on core parallel plasma beta ($\beta_{ec\parallel}$)\citep{Bercic2019}. Most of the strahl measured by PSP appear narrow with PAW reaching down to 30$^o$. The portion of the strahl velocity distribution function aligned with the magnetic field is for the measured energy range well described by a Maxwellian distribution function. T$_{s\parallel}$ was found to be anti-correlated with the solar wind velocity, and independent of radial distance. These observations imply that T$_{s\parallel}$ carries the information about the coronal electron temperature. The obtained values are in agreement with coronal temperatures measured using spectroscopy (David et al. 2998), and the inferred solar wind source regions during the first orbit of PSP agree with the predictions using a PFSS model (Bale et al. 2019, Badman et al. 2019).

Observations by the Parker Solar Probe mission of the solar wind at about 35.7 solar radii reveal the existence of whistler wave packets with frequencies below 0.1 f/fce (20-80 Hz in the spacecraft frame). These waves often coincide with local minima of the magnetic field magnitude or with sudden deflections of the magnetic field that are called switchbacks. Their sunward propagation leads to a significant Doppler frequency downshift from 200-300 Hz to 20-80 Hz (from 0.2 f/fce to 0.5 f/fce). The polarization of these waves varies from quasi-parallel to significantly oblique with wave normal angles that are close to the resonance cone. Their peak amplitude can be as large as 2 to 4 nT. Such values represent approximately 10% of the background magnetic field, which is considerably more than what is observed at 1 a.u. Recent numerical studies show that such waves may potentially play a key role in breaking the heat flux and scattering the Strahl population of suprathermal electrons into a halo population.

Electrostatic analyzers of different designs have been used since the earliest days of the space age, beginning with the very earliest solar wind measurements made by Mariner 2 en route to Venus in 1962. The Parker Solar Probe (PSP) mission, NASA's first dedicated mission to study the innermost reaches of the heliosphere, makes its thermal plasma measurements using a suite of instruments called the Solar Wind Electrons, Alphas, and Protons (SWEAP) investigation. SWEAP's electron Parker Solar Probe Analyzer (SPAN-E) instruments are a pair of top-hat electrostatic analyzers on PSP that are capable of measuring the electron distribution function in the solar wind from 2 eV to 30 keV. For the first time, in-situ measurements of thermal electrons provided by SPAN-E will help reveal the heating and acceleration mechanisms driving the evolution of the solar wind at the points of acceleration and heating, closer than ever before to the Sun. This paper details the design of the SPAN-E sensors and their operation, data formats, and measurement caveats from Parker Solar Probe's first two close encounters with the Sun.

The Parker Solar Probe was launched on 2018 August 12 and completed its second orbit on 2019 June 19 with perihelion of 35.7 solar radii. During this time, the Energetic particle Instrument-Hi (EPI-Hi, one of the two energetic particle instruments comprising the Integrated Science Investigation of the Sun, ISOIS) measured seven proton intensity increases associated with stream interaction regions (SIRs), two of which appear to be occurring in the same region corotating with the Sun. The events are relatively weak, with observed proton spectra extending to only a few MeV and lasting for a few days. The proton spectra are best characterized by power laws with indices ranging from -4.3 to -6.5, generally softer than events associated with SIRs observed at 1 au and beyond. Helium spectra were also obtained with similar indices, allowing He/H abundance ratios to be calculated for each event. We find values of 0.016-0.031, which are consistent with ratios obtained previously for corotating interaction region events with fast solar wind < 600 km s-1. Using the observed solar wind data combined with solar wind simulations, we study the solar wind structures associated with these events and identify additional spacecraft near 1 au appropriately positioned to observe the same structures after some corotation. Examination of the energetic particle observations from these spacecraft yields two events that may correspond to the energetic particle increases seen by EPI-Hi earlier.

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfv\'enic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfv\'en velocity shear, whereas the median shear across 43 Alfv\'enic current sheets examined was 71% of the Alfv\'en velocity shear. This finding could suggest that large, albeit sub-Alfv\'enic, velocity shears suppress reconnection. An alternative interpretation is that the Alfv\'enic current sheets are isolated rotational discontinuities which do not undergo local reconnection.

We perform a statistical study of the turbulent power spectrum at inertial and kinetic scales observed during the first perihelion encounter of Parker Solar Probe. We find that often there is an extremely steep scaling range of the power spectrum just above the ion-kinetic scales, similar to prior observations at 1 AU, with a power-law index of around $-4$. Based on our measurements, we demonstrate that either a significant ($>50\%$) fraction of the total turbulent energy flux is dissipated in this range of scales, or the characteristic nonlinear interaction time of the turbulence decreases dramatically from the expectation based solely on the dispersive nature of nonlinearly interacting kinetic Alfv\'en waves.

Radio waves are strongly scattered in the solar wind, so that their apparent sources seem to be considerably larger and shifted than the actual ones. Since the scattering depends on the spectrum of density turbulence, better understanding of the radio wave propagation provides indirect information on the relative density fluctuations $\epsilon=\langle\delta n\rangle/\langle n\rangle$ at the effective turbulence scale length. Here, we have analyzed 30 type III bursts detected by Parker Solar Probe (PSP). For the first time, we have retrieved type III burst decay times $\tau_{\rm{d}}$ between 1 MHz and 10 MHz thanks to an unparalleled temporal resolution of PSP. We observed a significant deviation in a power-law slope for frequencies above 1 MHz when compared to previous measurements below 1 MHz by the twin-spacecraft Solar TErrestrial RElations Observatory (STEREO) mission. We note that altitudes of radio bursts generated at 1 MHz roughly coincide with an expected location of the Alfv\'{e}n point, where the solar wind becomes super-Alfv\'{e}nic. By comparing PSP observations and Monte Carlo simulations, we predict relative density fluctuations $\epsilon$ at the effective turbulence scale length at radial distances between 2.5$R_\odot$ and 14$R_\odot$ to range from $0.22$ and $0.09$. Finally, we calculated relative density fluctuations $\epsilon$ measured in situ by PSP at a radial distance from the Sun of $35.7$~$R_\odot$ during the perihelion \#1, and the perihelion \#2 to be $0.07$ and $0.06$, respectively. It is in a very good agreement with previous STEREO predictions ($\epsilon=0.06-0.07$) obtained by remote measurements of radio sources generated at this radial distance.

The physical mechanisms that produce the slow solar wind are still highly debated. Parker Solar Probe's (PSP's) second solar encounter provided a new opportunity to relate in situ measurements of the nascent slow solar wind with white-light images of streamer flows. We exploit data taken by the Solar and Heliospheric Observatory (SOHO), the Solar TErrestrial RElations Observatory (STEREO) and the Wide Imager on Solar Probe to reveal for the first time a close link between imaged streamer flows and the high-density plasma measured by the Solar Wind Electrons Alphas and Protons (SWEAP) experiment. We identify different types of slow winds measured by PSP that we relate to the spacecraft's magnetic connectivity (or not) to streamer flows. SWEAP measured high-density and highly variable plasma when PSP was well connected to streamers but more tenuous wind with much weaker density variations when it exited streamer flows. STEREO imaging of the release and propagation of small transients from the Sun to PSP reveals that the spacecraft was continually impacted by the southern edge of streamer transients. The impact of specific density structures is marked by a higher occurrence of magnetic field reversals measured by the FIELDS magnetometers. Magnetic reversals originating from the streamers are associated with larger density variations compared with reversals originating outside streamers. We tentatively interpret these findings in terms of magnetic reconnection between open magnetic fields and coronal loops with different properties, providing support for the formation of a subset of the slow wind by magnetic reconnection.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

We consider 2D joint distributions of normalised residual energy $\sigma_r(s,t)$ and cross helicity $\sigma_c(s,t)$ during one day of Parker Solar Probe's (PSP's) first encounter as a function of wavelet scale $s$. The broad features of the distributions are similar to previous observations made by HELIOS in slow solar wind, namely well correlated and fairly Alfv\'enic, except for a population with negative cross helicity which is seen at shorter wavelet scales. We show that this population is due to the presence of magnetic switchbacks, brief periods where the magnetic field polarity reverses. Such switchbacks have been observed before, both in HELIOS data and in Ulysses data in the polar solar wind. Their abundance and short timescales as seen by PSP in its first encounter is a new observation, and their precise origin is still unknown. By analysing these MHD invariants as a function of wavelet scale we show that MHD waves do indeed follow the local mean magnetic field through switchbacks, with net Elsasser flux propagating inward during the field reversal, and that they therefore must be local kinks in the magnetic field and not due to small regions of opposite polarity on the surface of the Sun. Such observations are important to keep in mind as computing cross helicity without taking into account the effect of switchbacks may result in spurious underestimation of $\sigma_c$ as PSP gets closer to the Sun in later orbits.

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity.

Data from the first two orbits of the Sun by Parker Solar Probe reveal that the solar wind sunward of 50 solar radii is replete with plasma waves and instabilities. One of the most prominent plasma wave power enhancements in this region appears near the electron cyclotron frequency (f_ce). Most of this wave power is concentrated in electric field fluctuations near 0.7 f_ce and f_ce, with strong harmonics of both frequencies extending above f_ce. At least two distinct, often concurrent, wave modes are observed, preliminarily identified as electrostatic whistler-mode waves and electron Bernstein waves. Wave intervals range in duration from a few seconds to hours. Both the amplitudes and number of detections of these near-f_ce waves increase significantly with decreasing distance to the Sun, suggesting that they play an important role in the evolution of electron populations in the near-Sun solar wind. Correlations are found between the detection of these waves and properties of solar wind electron populations, including electron core drift, implying that these waves play a role in regulating the heat flux carried by solar wind electrons. Observation of these near-f_ce waves is found to be strongly correlated with near-radial solar wind magnetic field configurations with low levels of magnetic turbulence. A scenario for the growth of these waves is presented which implies that regions of low-turbulence near-radial magnetic field are a prominent feature of solar wind structure near the Sun.

During Parker Solar Probe's first orbit, the solar wind plasma has been observed in situ closer than ever before, the perihelion on November 6th 2018 revealing a flow that is constantly permeated by large amplitude Alfv\'enic fluctuations. These include radial magnetic field reversals, or switchbacks, that seem to be a persistent feature of the young solar wind. The measurements also reveal a very strong, unexpected, azimuthal velocity component. In this work, we numerically model the solar corona during this first encounter, solving the MHD equations and accounting for Alfv\'en wave transport and dissipation. We find that the large scale plasma parameters are well reproduced, allowing the computation of the solar wind sources at Probe with confidence. We try to understand the dynamical nature of the solar wind to explain both the amplitude of the observed radial magnetic field and of the azimuthal velocities.

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

We present initial results from the Radio Frequency Spectrometer (RFS), the high frequency component of the FIELDS experiment on the Parker Solar Probe (PSP). During the first PSP solar encounter (2018 November), only a few small radio bursts were observed. During the second encounter (2019 April), copious Type III radio bursts occurred, including intervals of radio storms where bursts occurred continuously. In this paper, we present initial observations of the characteristics of Type III radio bursts in the inner heliosphere, calculating occurrence rates, amplitude distributions, and spectral properties of the observed bursts. We also report observations of several bursts during the second encounter which display circular polarization in the right hand polarized sense, with a degree of polarization of 0.15-0.38 in the range from 8-12 MHz. The degree of polarization can be explained either by depolarization of initially 100% polarized $o$-mode emission, or by direct generation of emission in the $o$ and $x$-mode simultaneously. Direct in situ observations in future PSP encounters could provide data which can distinguish these mechanisms.

Parker Solar Probe (PSP) observations during its first encounter at 35.7 $R_\odot$ have shown the presence of magnetic field lines which are strongly perturbed to the point that they produce local inversions of the radial magnetic field, known as switchbacks. Their counterparts in the solar wind velocity field are local enhancements in the radial speed, or jets, displaying (in all components) the velocity-magnetic field correlation typical of large amplitude Alfv\'en waves propagating away from the Sun. Switchbacks and radial jets have previously been observed over a wide range of heliocentric distances by Helios, WIND and Ulysses, although they were prevalent in significantly faster streams than seen at PSP. Here we study via numerical MHD simulations the evolution of such large amplitude Alfv\'enic fluctuations by including, in agreement with observations, both a radial magnetic field inversion and an initially constant total magnetic pressure. Despite the extremely large excursion of magnetic and velocity fields, switchbacks are seen to persist for up to hundreds of Alfv\'en crossing times before eventually decaying due to the parametric decay instability. Our results suggest that such switchback/jet configurations might indeed originate in the lower corona and survive out to PSP distances, provided the background solar wind is sufficiently calm, in the sense of not being pervaded by strong density fluctuations or other gradients, such as stream or magnetic field shears, that might destabilize or destroy them over shorter timescales.

A series of solar energetic particle (SEP) events were observed at Parker Solar Probe (PSP) by the Integrated Science Investigation of the Sun (\ISOIS) during the period from April 18, 2019 through April 24, 2019. The PSP spacecraft was located near 0.48 au from the Sun on Parker spiral field lines that projected out to 1 au within $\sim 25^\circ$ of near Earth spacecraft. These SEP events, though small compared to historically large SEP events, were amongst the largest observed thus far in the PSP mission and provide critical information about the space environment inside 1 au during SEP events. During this period the Sun released multiple coronal mass ejections (CMEs). One of these CMEs observed was initiated on April 20, 2019 at 01:25 UTC, and the interplanetary CME (ICME) propagated out and passed over the PSP spacecraft. Observations by the Electromagnetic Fields Investigation (FIELDS) show that the magnetic field structure was mostly radial throughout the passage of the compression region and the plasma that followed, indicating that PSP did not directly observe a flux rope internal to the ICME, consistent with the location of PSP on the ICME flank. Analysis using relativistic electrons observed near Earth by the Electron, Proton and Alpha Monitor (EPAM) on the Advanced Composition Explorer (ACE) demonstrates the presence of electron seed populations (40--300 keV) during the events observed. The energy spectrum of the \ISOIS~ observed proton seed population below 1 MeV is close to the limit of possible stationary state plasma distributions out of equilibrium. \ISOIS~ observations reveal the \revise{enhancement} of seed populations during the passage of the ICME, which \revise{likely indicates a key part} of the pre-acceleration process that occurs close to the Sun.

The Solar Probe Cup (SPC) is a Faraday Cup instrument onboard NASA's Parker Solar Probe (PSP) spacecraft designed to make rapid measurements of thermal coronal and solar wind plasma. The spacecraft is in a heliocentric orbit that takes it closer to the Sun than any previous spacecraft, allowing measurements to be made where the coronal and solar wind plasma is being heated and accelerated. The SPC instrument was designed to be pointed directly at the Sun at all times, allowing the solar wind (which is flowing primarily radially away from the Sun) to be measured throughout the orbit. The instrument is capable of measuring solar wind ions with an energy/charge between 100 V and 6000 V (protons with speeds from $139-1072~km~s^{-1})$. It also measures electrons with an energy between 100 V and 1500 V. SPC has been designed to have a wide dynamic range that is capable of measuring protons and alpha particles at the closest perihelion (9.86 solar radii from the center of the Sun) and out to 0.25 AU. Initial observations from the first orbit of PSP indicate that the instrument is functioning well.

Since the launch on 2018/08/12, Parker Solar Probe (PSP) has completed its first and second orbits around the Sun, having reached down to 35.7 solar radii at each perihelion. In anticipation of the exciting new data at such unprecedented distances, we have simulated the global 3D heliosphere using an MHD model coupled with a semi-empirical coronal model using the best available photospheric magnetograms as input. We compare our heliospheric MHD simulation results with in situ measurements along the PSP trajectory from its launch to the completion of the second orbit, with particular emphasis on the solar wind structure around the first two solar encounters. Furthermore, we show our model prediction for the third perihelion, which occurred on 2019/09/01. Comparison of the MHD results with PSP observations provides a new insight on the solar wind acceleration. Moreover, PSP observations reveal how accurately the ADAPT-WSA predictions work throughout the inner heliosphere.

Stochastic heating is a non-linear heating mechanism driven by the violation of magnetic moment invariance due to large-amplitude turbulent fluctuations producing diffusion of ions towards higher kinetic energies in the direction perpendicular to the magnetic field. It is frequently invoked as a mechanism responsible for the heating of ions in the solar wind. Here, we quantify for the first time the proton stochastic heating rate $Q_\perp$ at radial distances from the Sun as close as $0.16$ au, using measurements from the first two Parker Solar Probe encounters. Our results for both the amplitude and radial trend of the heating rate, $Q_\perp \propto r^{-2.5}$, agree with previous results based on the Helios data set at heliocentric distances from 0.3 to 0.9 au. Also in agreement with previous results, $Q_\perp$ is significantly larger in the fast solar wind than in the slow solar wind. We identify the tendency in fast solar wind for cuts of the core proton velocity distribution transverse to the magnetic field to exhibit a flat-top shape. The observed distribution agrees with previous theoretical predictions for fast solar wind where stochastic heating is the dominant heating mechanism.

The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfv\'enic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfv\'en speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfv\'en point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

We compare magnetic field measurements taken by the FIELDS instrument on Parker Solar Probe (PSP) during its first solar encounter to predictions obtained by Potential Field Source Surface (PFSS) modeling. Ballistic propagation is used to connect the spacecraft to the source surface. Despite the simplicity of the model, our results show striking agreement with PSPs first observations of the heliospheric magnetic field from 0.5 AU (107.5 Rs) down to 0.16 AU (35.7 Rs). Further, we show the robustness of the agreement is improved both by allowing the photospheric input to the model to vary in time, and by advecting the field from PSP down to the PFSS model domain using in situ PSP/SWEAP measurements of the solar wind speed instead of assuming it to be constant with longitude and latitude. We also explore the source surface height parameter (RSS) to the PFSS model finding that an extraordinarily low source surface height (1.3-1.5Rs) predicts observed small scale polarity inversions which are otherwise washed out with regular modeling parameters. Finally, we extract field line traces from these models. By overlaying these on EUV images we observe magnetic connectivity to various equatorial and mid-latitude coronal holes indicating plausible magnetic footpoints and offering context for future discussions of sources of the solar wind measured by PSP.

The Solar Wind Electrons Alphas and Protons experiment on the Parker Solar Probe (PSP) mission measures the three-dimensional electron velocity distribution function. We derive the parameters of the core, halo, and strahl populations utilizing a combination of fitting to model distributions and numerical integration for $\sim 100,000$ electron distributions measured near the Sun on the first two PSP orbits, which reached heliocentric distances as small as $\sim 0.17$ AU. As expected, the electron core density and temperature increase with decreasing heliocentric distance, while the ratio of electron thermal pressure to magnetic pressure ($\beta_e$) decreases. These quantities have radial scaling consistent with previous observations farther from the Sun, with superposed variations associated with different solar wind streams. The density in the strahl also increases; however, the density of the halo plateaus and even decreases at perihelion, leading to a large strahl/halo ratio near the Sun. As at greater heliocentric distances, the core has a sunward drift relative to the proton frame, which balances the current carried by the strahl, satisfying the zero-current condition necessary to maintain quasi-neutrality. Many characteristics of the electron distributions near perihelion have trends with solar wind flow speed, $\beta_e$, and/or collisional age. Near the Sun, some trends not clearly seen at 1 AU become apparent, including anti-correlations between wind speed and both electron temperature and heat flux. These trends help us understand the mechanisms that shape the solar wind electron distributions at an early stage of their evolution.

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes.

Understanding the physical processes in the solar wind and corona which actively contribute to heating, acceleration, and dissipation is a primary objective of NASA's Parker Solar Probe (PSP) mission. Observations of coherent electromagnetic waves at ion scales suggests that linear cyclotron resonance and non-linear processes are dynamically relevant in the inner heliosphere. A wavelet-based statistical study of coherent waves in the first perihelion encounter of PSP demonstrates the presence of transverse electromagnetic waves at ion resonant scales which are observed in 30-50\% of radial field intervals. Average wave amplitudes of approximately 4 nT are measured, while the mean duration of wave events is of order 20 seconds; however long duration wave events can exist without interruption on hour-long timescales. Though ion scale waves are preferentially observed during intervals with a radial mean magnetic field, we show that measurement constraints, associated with single spacecraft sampling of quasi-parallel waves superposed with anisotropic turbulence, render the measured quasi-parallel ion-wave spectrum unobservable when the mean magnetic field is oblique to the solar wind flow; these results imply that the occurrence of coherent ion-scale waves is not limited to a radial field configuration. The lack of strong radial scaling of characteristic wave amplitudes and duration suggests that the waves are generated {\em{in-situ}} through plasma instabilities. Additionally, observations of proton distribution functions indicate that temperature anisotropy may drive the observed ion-scale waves.

The Radiation Assessment Detector (RAD), on board Mars Science Laboratory's (MSL) Curiosity rover, has been measuring ground level particle fluxes along with the radiation dose rate at the surface of Mars since August 2012. Similar to neutron monitors at Earth, RAD sees many Forbush decreases (FDs) in the galactic cosmic ray (GCR) induced surface fluxes and dose rates. These FDs are associated with coronal mass ejections (CMEs) and/or stream/corotating interaction regions (SIRs/CIRs). Orbiting above the Martian atmosphere, the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft has also been monitoring space weather conditions at Mars since September 2014. The penetrating particle flux channels in the Solar Energetic Particle (SEP) instrument onboard MAVEN can also be employed to detect FDs. For the first time, we study the statistics and properties of a list of FDs observed in-situ at Mars, seen both on the surface by MSL/RAD and in orbit detected by the MAVEN/SEP instrument. Such a list of FDs can be used for studying interplanetary CME (ICME) propagation and SIR evolution through the inner heliosphere. The magnitudes of different FDs can be well-fitted by a power-law distribution. The systematic difference between the magnitudes of the FDs within and outside the Martian atmosphere may be mostly attributed to the energy-dependent modulation of the GCR particles by both the pass-by ICMEs/SIRs and the Martian atmosphere.

We consider the effectiveness of foreground cleaning in the recovery of Cosmic Microwave Background (CMB) polarization sourced by gravitational waves for tensor-to-scalar ratios in the range $0<r<0.1$. Using the planned survey area, frequency bands, and sensitivity of the Cosmology Large Angular Scale Surveyor (CLASS), we simulate maps of Stokes $Q$ and $U$ parameters at 40, 90, 150, and 220 GHz, including realistic models of the CMB, diffuse Galactic thermal dust and synchrotron foregrounds, and Gaussian white noise. We use linear combinations (LCs) of the simulated multifrequency data to obtain maximum likelihood estimates of $r$, the relative scalar amplitude $s$, and LC coefficients. We find that for 10,000 simulations of a CLASS-like experiment using only measurements of the reionization peak ($\ell\leq23$), there is a 95% C.L. upper limit of $r<0.017$ in the case of no primordial gravitational waves. For simulations with $r=0.01$, we recover at 68% C.L. $r=0.012^{+0.011}_{-0.006}$. The reionization peak corresponds to a fraction of the multipole moments probed by CLASS, and simulations including $30\leq\ell\leq100$ further improve our upper limits to $r<0.008$ at 95% C.L. ($r=0.01^{+0.004}_{-0.004}$ for primordial gravitational waves with $r=0.01$). In addition to decreasing the current upper bound on $r$ by an order of magnitude, these foreground-cleaned low multipole data will achieve a cosmic variance limited measurement of the E-mode polarization's reionization peak.

We examine the consistency of WMAP9 and Planck data. We compare sky maps, power spectra, and inferred LCDM cosmological parameters. Residual dipoles are seen in the WMAP and Planck sky map differences, but are consistent within the uncertainties and are not large enough to explain the widely-noted differences in angular power spectra at higher l. After removing residual dipoles and galactic foregrounds, the residual difference maps exhibit a quadrupole and other large-scale systematic structure. We identify this structure as possibly originating from Planck's beam sidelobe pick-up, but note that it appears to have insignificant cosmological impact. We develop an extension of the internal linear combination technique and find features that plausibly originate in the Planck data. We examine LCDM model fits to the angular power spectra and conclude that the ~2.5% difference in the spectra at multipoles greater than l~100 are significant at the 3-5 sigma level. We revisit the analysis of WMAP's beam data and conclude that previously-derived uncertainties are robust and cannot explain the power spectrum differences. Finally, we examine the consistency of the LCDM parameters inferred from each data set taking into account the fact that both experiments observe the same sky, but cover different multipole ranges, apply different sky masks, and have different noise. While individual parameter values agree within the uncertainties, the 6 parameters taken together are discrepant at the ~6 sigma level, with chi2=56 for 6 dof (PTE = 3e-10). Of the 6 parameters, chi2 is best improved by marginalizing over Omega_c h^2, giving chi2=5.2 for 5 degrees of freedom. We find that perturbing the WMAP window function by its dominant beam error profile has little effect on Omega_c h^2, while perturbing the Planck window function by its corresponding error profile has a much greater effect on Omega_c h^2.

The determination of the Hubble constant has been a central goal in observational astrophysics for nearly 100 years. Extraordinary progress has occurred in recent years on two fronts: the cosmic distance ladder measurements at low redshift and cosmic microwave background (CMB) measurements at high redshift. The CMB is used to predict the current expansion rate through a best-fit cosmological model. Complementary progress has been made with baryon acoustic oscillation (BAO) measurements at relatively low redshifts. While BAO data do not independently determine a Hubble constant, they are important for constraints on possible solutions and checks on cosmic consistency. A precise determination of the Hubble constant is of great value, but it is more important to compare the high and low redshift measurements to test our cosmological model. Significant tension would suggest either uncertainties not accounted for in the experimental estimates, or the discovery of new physics beyond the standard model of cosmology. In this paper we examine in detail the tension between the CMB, BAO, and cosmic distance ladder data sets. We find that these measurements are consistent within reasonable statistical expectations, and we combine them to determine a best-fit Hubble constant of 69.6+/-0.7 km/s/Mpc. This value is based upon WMAP9+SPT+ACT+6dFGS+BOSS/DR11+H_0/Riess; we explore alternate data combinations in the text. The combined data constrain the Hubble constant to 1%, with no compelling evidence for new physics.

The Cosmology Large Angular Scale Surveyor (CLASS) experiment aims to map the polarization of the Cosmic Microwave Background (CMB) at angular scales larger than a few degrees. Operating from Cerro Toco in the Atacama Desert of Chile, it will observe over 65% of the sky at 38, 93, 148, and 217 GHz. In this paper we discuss the design, construction, and characterization of the CLASS 38 GHz detector focal plane, the first ever Q-band bolometric polarimeter array.

The Cosmology Large Angular Scale Surveyor (CLASS) is an experiment to measure the signature of a gravita-tional-wave background from inflation in the polarization of the cosmic microwave background (CMB). CLASS is a multi-frequency array of four telescopes operating from a high-altitude site in the Atacama Desert in Chile. CLASS will survey 70\% of the sky in four frequency bands centered at 38, 93, 148, and 217 GHz, which are chosen to straddle the Galactic-foreground minimum while avoiding strong atmospheric emission lines. This broad frequency coverage ensures that CLASS can distinguish Galactic emission from the CMB. The sky fraction of the CLASS survey will allow the full shape of the primordial B-mode power spectrum to be characterized, including the signal from reionization at low $\ell$. Its unique combination of large sky coverage, control of systematic errors, and high sensitivity will allow CLASS to measure or place upper limits on the tensor-to-scalar ratio at a level of $r=0.01$ and make a cosmic-variance-limited measurement of the optical depth to the surface of last scattering, $\tau$.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

[Abridged] We present updated estimates of Galactic foreground emission using seven years of WMAP data. Using the power spectrum of differences between multi-frequency template-cleaned maps, we find no evidence for foreground contamination outside of the updated (KQ85y7) foreground mask. We place a 15 microKelvin upper bound on rms foreground contamination in the cleaned maps used for cosmological analysis. We find no indication in the polarization data of an extra "haze" of hard synchrotron emission from energetic electrons near the Galactic center. We provide an updated map of the cosmic microwave background (CMB) using the internal linear combination (ILC) method, updated foreground masks, and updates to point source catalogs with 62 newly detected sources. Also new are tests of the Markov chain Monte Carlo (MCMC) foreground fitting procedure against systematics in the time-stream data, and tests against the observed beam asymmetry. Within a few degrees of the Galactic plane, WMAP total intensity data show a rapidly steepening spectrum from 20-40 GHz, which may be due to emission from spinning dust grains, steepening synchrotron, or other effects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepening synchrotron models fit the combination of WMAP and 408 MHz data equally well. ARCADE data appear inconsistent with the steepening synchrotron model, and consistent with the spinning dust model, though some discrepancies remain regarding the relative strength of spinning dust emission. More high-resolution data in the 10-40 GHz range would shed much light on these issues.

(Abridged) The 7-year WMAP data and improved astrophysical data rigorously test the standard cosmological model and its extensions. By combining WMAP with the latest distance measurements from BAO and H0 measurement, we determine the parameters of the simplest LCDM model. The power-law index of the primordial power spectrum is n_s=0.968+-0.012, a measurement that excludes the scale-invariant spectrum by 99.5%CL. The other parameters are also improved from the 5-year results. Notable examples of improved parameters are the total mass of neutrinos, sum(m_nu)<0.58eV, and the effective number of neutrino species, N_eff=4.34+0.86-0.88. We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis. We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z=1090 and the dominance of adiabatic scalar fluctuations. With the 7-year TB power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved to Delta(alpha)=-1.1+-1.4(stat)+-1.5(syst) degrees. We report significant detections of the SZ effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data. However, it is a factor of 0.5 to 0.7 times the predictions from "universal profile" of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically-expected SZ power spectrum recently measured by the South Pole Telescope collaboration.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

We describe a sampling method to estimate the polarized CMB signal from observed maps of the sky. We use a Metropolis-within-Gibbs algorithm to estimate the polarized CMB map, containing Q and U Stokes parameters at each pixel, and its covariance matrix. These can be used as inputs for cosmological analyses. The polarized sky signal is parameterized as the sum of three components: CMB, synchrotron emission, and thermal dust emission. The polarized Galactic components are modeled with spatially varying power law spectral indices for the synchrotron, and a fixed power law for the dust, and their component maps are estimated as by-products. We apply the method to simulated low resolution maps with pixels of side 7.2 degrees, using diagonal and full noise realizations drawn from the WMAP noise matrices. The CMB maps are recovered with goodness of fit consistent with errors. Computing the likelihood of the E-mode power in the maps as a function of optical depth to reionization, tau, for fixed temperature anisotropy power, we recover tau=0.091+-0.019 for a simulation with input tau=0.1, and mean tau=0.098 averaged over 10 simulations. A `null' simulation with no polarized CMB signal has maximum likelihood consistent with tau=0. The method is applied to the five-year WMAP data, using the K, Ka, Q and V channels. We find tau=0.090+-0.019, compared to tau=0.086+-0.016 from the template-cleaned maps used in the primary WMAP analysis. The synchrotron spectral index, beta, averaged over high signal-to-noise pixels with standard deviation sigma(beta)<0.25, but excluding ~6% of the sky masked in the Galactic plane, is -3.03+-0.04. This estimate does not vary significantly with Galactic latitude, although includes an informative prior.

We present an alternate interpretation of recent STEREO/STE observations that were originally attributed to energetic neutral atoms (ENA) from the heliosheath. The signal attributed to the diffuse ENA source instead shows the characteristics of a point source. We point out that the peak intensity seen by STEREO/STE is centered at the ecliptic longitude of the bright X-ray source Sco X-1. The observed energy spectrum and intensity are also consistent with the X-rays from Sco X-1. The problem of energy dissipation at the solar wind termination shock remains unsolved while current understanding of the interaction between the solar wind and interstellar wind awaits future observations.

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced chi-squared of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

This paper focuses on cosmological constraints derived from analysis of WMAP data alone. A simple LCDM cosmological model fits the five-year WMAP temperature and polarization data. The basic parameters of the model are consistent with the three-year data and now better constrained: Omega_b h^2 = 0.02273+-0.00062, Omega_c h^2 = 0.1099+-0.0062, Omega_L = 0.742+-0.030, n_s = 0.963+0.014- 0.015, tau = 0.087+-0.017, sigma_8 = 0.796+-0.036. With five years of polarization data, we have measured the optical depth to reionization, tau>0, at 5 sigma significance. The redshift of an instantaneous reionization is constrained to be z_reion = 11.0+-1.4 with 68% confidence. This excludes a sudden reionization of the universe at z=6 at more than 3.5 sigma significance, suggesting that reionization was an extended process. Using two methods for polarized foreground cleaning we get consistent estimates for the optical depth, indicating an error due to foreground treatment of tau~0.01. This cosmological model also fits small-scale CMB data, and a range of astronomical data measuring the expansion rate and clustering of matter in the universe. We find evidence for the first time in the CMB power spectrum for a non-zero cosmic neutrino background, or a background of relativistic species, with the standard three light neutrino species preferred over the best-fit LCDM model with N_eff=0 at >99.5% confidence, and N_eff > 2.3 (95% CL) when varied. The five-year WMAP data improve the upper limit on the tensor-to-scalar ratio to r < 0.43 (95% CL), for power-law models. With longer integration we find no evidence for a running spectral index, with dn_s/dlnk = -0.037+-0.028.

We present the list of point sources found in the WMAP 5-year maps. The technique used in the first-year and three-year analysis now finds 390 point sources, and the five-year source catalog is complete for regions of the sky away from the galactic plane to a 2 Jy limit, with SNR > 4.7 in all bands in the least covered parts of the sky. The noise at high frequencies is still mainly radiometer noise, but at low frequencies the CMB anisotropy is the largest uncertainty. A separate search of CMB-free V-W maps finds 99 sources of which all but one can be identified with known radio sources. The sources seen by WMAP are not strongly polarized. Many of the WMAP sources show significant variability from year to year, with more than a 2:1 range between the minimum and maximum fluxes.

(Abridged) The WMAP 5-year data strongly limit deviations from the minimal LCDM model. We constrain the physics of inflation via Gaussianity, adiabaticity, the power spectrum shape, gravitational waves, and spatial curvature. We also constrain the properties of dark energy, parity-violation, and neutrinos. We detect no convincing deviations from the minimal model. The parameters of the LCDM model, derived from WMAP combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO), are: Omega_b=0.0456+-0.0015, Omega_c=0.228+-0.013, Omega_Lambda=0.726+-0.015, H_0=70.5+-1.3 km/s/Mpc, n_s=0.960+-0.013, tau=0.084+-0.016, and sigma_8=0.812+-0.026. With WMAP+BAO+SN, we find the tensor-to-scalar ratio r<0.22 (95% CL), and n_s>1 is disfavored regardless of r. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and curvature. We provide a set of "WMAP distance priors," to test a variety of dark energy models. We test a time-dependent w with a present value constrained as -0.33<1+w_0<0.21 (95% CL). Temperature and matter fluctuations obey the adiabatic relation to within 8.9% and 2.1% for the axion and curvaton-type dark matter, respectively. The TE and EB spectra constrain cosmic parity-violation. We find the limit on the total mass of neutrinos, sum(m_nu)<0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large-scale structure data. The effective number of neutrino species is constrained as N_{eff} = 4.4+-1.5 (68%), consistent with the standard value of 3.04. Finally, limits on primordial non-Gaussianity are -9<f_{NL}^{local}<111 and -151<f_{NL}^{equil}<253 (95% CL) for the local and equilateral models, respectively.

Cosmology and other scientific results from the WMAP mission require an accurate knowledge of the beam patterns in flight. While the degree of beam knowledge for the WMAP one-year and three-year results was unprecedented for a CMB experiment, we have significantly improved the beam determination as part of the five-year data release. Physical optics fits are done on both the A and the B sides for the first time. The cutoff scale of the fitted distortions on the primary mirror is reduced by a factor of ~2 from previous analyses. These changes enable an improvement in the hybridization of Jupiter data with beam models, which is optimized with respect to error in the main beam solid angle. An increase in main-beam solid angle of ~1% is found for the V2 and W1-W4 differencing assemblies. Although the five-year results are statistically consistent with previous ones, the errors in the five-year beam transfer functions are reduced by a factor of ~2 as compared to the three-year analysis. We present radiometry of the planet Jupiter as a test of the beam consistency and as a calibration standard; for an individual differencing assembly, errors in the measured disk temperature are ~0.5%.

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that Ka band data (33 GHz) is suitable for use in cosmological analysis, after foreground cleaning. This significantly reduces the overall polarization uncertainty. With the 5 year WMAP data, we detect no convincing deviations from the minimal 6-parameter LCDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations, we find (68% CL uncertainties): Omega_bh^2 = 0.02267 \pm 0.00059, Omega_ch^2 = 0.1131 \pm 0.0034, Omega_Lambda = 0.726 \pm 0.015, n_s = 0.960 \pm 0.013, tau = 0.084 \pm 0.016, and Delta_R^2 = (2.445 \pm 0.096) x 10^-9. From these we derive: sigma_8 = 0.812 \pm 0.026, H_0 = 70.5 \pm 1.3 km/s/Mpc, z_{reion} = 10.9 \pm 1.4, and t_0 = 13.72 \pm 0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.22 (95% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and spatial curvature: -0.14 < 1+w < 0.12 and -0.0179 < Omega_k < 0.0081 (both 95% CL). The number of relativistic degrees of freedom (e.g. neutrinos) is found to be N_{eff} = 4.4 \pm 1.5, consistent with the standard value of 3.04. Models with N_{eff} = 0 are disfavored at >99.5% confidence.

We present the temperature and polarization angular power spectra of the cosmic microwave background (CMB) derived from the first 5 years of WMAP data. The 5-year temperature (TT) spectrum is cosmic variance limited up to multipole l=530, and individual l-modes have S/N>1 for l<920. The best fitting six-parameter LambdaCDM model has a reduced chi^2 for l=33-1000 of chi^2/nu=1.06, with a probability to exceed of 9.3%. There is now significantly improved data near the third peak which leads to improved cosmological constraints. The temperature-polarization correlation (TE) is seen with high significance. After accounting for foreground emission, the low-l reionization feature in the EE power spectrum is preferred by \Delta\chi^2=19.6 for optical depth tau=0.089 by the EE data alone, and is now largely cosmic variance limited for l=2-6. There is no evidence for cosmic signal in the BB, TB, or EB spectra after accounting for foreground emission. We find that, when averaged over l=2-6, l(l+1)C^{BB}_l/2\pi < 0.15 uK^2 (95% CL).

Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method.

We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.

We present a careful frequentist analysis of one- and two-point statistics of the hot and cold spots in the cosmic microwave background (CMB) data obtained by the Wilkinson Microwave Anisotropy Probe (WMAP). Our main result is the detection of a new anomaly at the 3-sigma level using temperature-weighted extrema correlation functions. We obtain this result using a simple hypothesis test which reduces the maximum risk of a false detection to the same level as the claimed significance of the test. We further present a detailed study of the robustness of our earlier claim (Larson and Wandelt 2004) under variations in the noise model and in the resolution of the map. Free software which implements our test is available online.

We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.

We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant. Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions.

This paper presents a frequentist analysis of the hot and cold spots of the cosmic microwave background data collected by the Wilkinson Microwave Anisotropy Probe (WMAP). We compare the WMAP temperature statistics of extrema (number of extrema, mean excursion, variance, skewness and kurtosis of the excursion) to Monte-Carlo simulations. We find that, on average, the local maxima (high temperatures in the anisotropy) are too cold and the local minima are too warm. In order to quantify this claim we describe a two-sided statistical hypothesis test which we advocate for other investigations of the Gaussianity hypothesis. Using this test we reject the isotropic Gaussian hypothesis at more than 99% confidence in a well-defined way. Our claims are based only on regions that are outside the most conservative WMAP foreground mask. We perform our test separately on maxima and minima, and on the north and south ecliptic and Galactic hemispheres and reject Gaussianity at above 95% confidence for almost all tests of the mean excursions. The same test also shows the variance of the maxima and minima to be low in the ecliptic north (99% confidence), but consistent in the south; this effect is not as pronounced in the Galactic north and south hemispheres.

We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (N_side = 128) and high (N_side = 512) resolutions, to study effects such as incomplete sky coverage and white vs. correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio, and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current super-computer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter.

We describe an efficient and exact method that enables global Bayesian analysis of cosmic microwave background (CMB) data. The method reveals the joint posterior density (or likelihood for flat priors) of the power spectrum $C_\ell$ and the CMB signal. Foregrounds and instrumental parameters can be simultaneously inferred from the data. The method allows the specification of a wide range of foreground priors. We explicitly show how to propagate the non-Gaussian dependency structure of the $C_\ell$ posterior through to the posterior density of the parameters. If desired, the analysis can be coupled to theoretical (cosmological) priors and can yield the posterior density of cosmological parameter estimates directly from the time-ordered data. The method does not hinge on special assumptions about the survey geometry or noise properties, etc. It is based on a Monte Carlo approach and hence parallelizes trivially. No trace or determinant evaluations are necessary. The feasibility of this approach rests on the ability to solve the systems of linear equations which arise. These are of the same size and computational complexity as the map-making equations. We describe a pre-conditioned conjugate gradient technique that solves this problem and demonstrate in a numerical example that the computational time required for each Monte Carlo sample scales as $n_p^{3/2}$ with the number of pixels $n_p$. We test our method using the COBE-DMR data and explore the non-Gaussian joint posterior density of the COBE-DMR $C_\ell$ in several projections.