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

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 {\mu}K degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12\mu m). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are detected at 353 GHz and 143 GHz, with amplitudes around -7 % and +13 %, respectively, indicating that the measured spectra are likely more complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 100 GHz bands, mainly caused by to the contribution of the J=2-1 and J=1-0 12CO and 13CO emission lines. We derive maps of the dust temperature T, the dust spectral emissivity index beta, and the dust optical depth at 250 microns tau. The temperature map illustrates the cooling of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant T-beta anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest molecular regions to the faint diffuse regions. We use near-infrared extinction and HI data at 21-cm to perform a quantitative analysis of the spatial variations of the measured dust optical depth at 250 microns per hydrogen atom tau/NH. We report an increase of tau/NH by a factor of about 2 between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles.

Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively.

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

The data reported in Planck's Early Release Compact Source Catalogue (ERCSC) are exploited to measure the number counts (dN/dS) of extragalactic radio sources at 30, 44, 70, 100, 143 and 217 GHz. Due to the full-sky nature of the catalogue, this measurement extends to the rarest and brightest sources in the sky. At lower frequencies (30, 44, and 70 GHz) our counts are in very good agreement with estimates based on WMAP data, being somewhat deeper at 30 and 70 GHz, and somewhat shallower at 44 GHz. Planck's source counts at 143 and 217 GHz join smoothly with the fainter ones provided by the SPT and ACT surveys over small fractions of the sky. An analysis of source spectra, exploiting Planck's uniquely broad spectral coverage, finds clear evidence of a steepening of the mean spectral index above about 70 GHz. This implies that, at these frequencies, the contamination of the CMB power spectrum by radio sources below the detection limit is significantly lower than previously estimated.

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

We construct an all-sky map of the apparent temperature and optical depth of thermal dust emission using the Planck-HFI and IRAS data. The optical depth maps are correlated to tracers of the atomic and molecular gas. The correlation is linear in the lowest column density regions at high galactic latitudes. At high NH, the correlation is consistent with that of the lowest NH. In the intermediate NH range, we observe departure from linearity, with the dust optical depth in excess to the correlation. We attribute this excess emission to thermal emission by dust associated with a Dark-Gas phase, undetected in the available HI and CO measurements. We show the 2D spatial distribution of the Dark-Gas in the solar neighborhood and show that it extends around known molecular regions traced by CO. The average dust emissivity in the HI phase in the solar neighborhood follows roughly a power law distribution with beta = 1.8 all the way down to 3 mm, although the SED flattens slightly in the millimetre. The threshold for the existence of the Dark-Gas is found at NH = (8.0\pm 0.58) 10^{20} Hcm-2. Assuming the same dust emissivity at high frequencies for the dust in the atomic and molecular phases leads to an average XCO = (2.54\pm0.13) 10^{20} H2cm-2/(K km s-1). The mass of Dark-Gas is found to be 28% of the atomic gas and 118% of the CO emitting gas in the solar neighborhood. A possible origin for the Dark-Gas is the existence of a dark molecular phase, where H2 survives photodissociation but CO does not. The observed transition for the onset of this phase in the solar neighborhood (AV = 0.4 mag) appears consistent with recent theoretical predictions. We also discuss the possibility that up to half of the Dark-Gas could be in atomic form, due to optical depth effects in the HI measurements.

The all-sky coverage of the Planck Early Release Compact Source Catalogue (ERCSC) provides an unsurpassed survey of galaxies at submillimetre (submm) wavelengths, representing a major improvement in the numbers of galaxies detected, as well as the range of far-IR/submm wavelengths over which they have been observed. We here present the first results on the properties of nearby galaxies using these data. We match the ERCSC catalogue to IRAS-detected galaxies in the Imperial IRAS Faint Source Redshift Catalogue (IIFSCz), so that we can measure the spectral energy distributions (SEDs) of these objects from 60 to 850 microns. This produces a list of 1717 galaxies with reliable associations between Planck and IRAS, from which we select a subset of 468 for SED studies, namely those with strong detections in the three highest frequency Planck bands and no evidence of cirrus contamination. The SEDs are fitted using parametric dust models to determine the range of dust temperatures and emissivities. We find evidence for colder dust than has previously been found in external galaxies, with T<20K. Such cold temperatures are found using both the standard single temperature dust model with variable emissivity beta, or a two dust temperature model with beta fixed at 2. We also compare our results to studies of distant submm galaxies (SMGs) which have been claimed to contain cooler dust than their local counterparts. We find that including our sample of 468 galaxies significantly reduces the distinction between the two populations. Fits to SEDs of selected objects using more sophisticated templates derived from radiative transfer models confirm the presence of the colder dust found through parameteric fitting. We thus conclude that cold (T<20K) dust is a significant and largely unexplored component of many nearby galaxies.

We discuss methods to compute maps of the CMB in models featuring active causal sources and in non-Gaussian models ofinflation. We show our large angle results as well as some preliminary results on small angles. We conclude by discussing on-going work.

Much of the foundational work on quantum cosmology employs a simple minisuperspace model describing a Friedmann-Robertson-Walker universe containing a massive scalar field. We show that the classical limit of this model exhibits deterministic chaos and explore some of the consequences for the quantum theory. In particular, the breakdown of the WKB approximation calls into question many of the standard results in quantum cosmology.