Engineering Papers⌕ Search

Engineering topics

Wollack, E.

Publications and source records attributed to Wollack, E..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Horn Coupled Multichroic Polarimeters for the Atacama Cosmology Telescope Polarization Experiment

Multichroic polarization sensitive detectors enable increased sensitivity and spectral coverage for observations of the cosmic microwave background. An array optimized for dual frequency detectors can provide 1.7 times gain in sensitivity compared to a single frequency array. We present the design and measurements of horn coupled multichroic polarimeters encompassing the 90 and 150 GHz frequency bandsand discuss our plans to field an array of these detectors as part of the ACTPol project.

cosmic microwave background↗
NASA NTRS2013

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also 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 cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM 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 (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 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 the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −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 Lambda-CDM 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 Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗
NASA NTRS2013

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results

We present cosmological parameter constraints based on the final nine-year Wilkinson Microwave Anisotropy Probe (WMAP) data, in conjunction with a number of additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter Lambda-CDM model. When WMAP data are combined with measurements of the high-l cosmic microwave background anisotropy, the baryon acoustic oscillation scale, and the Hubble constant, the matter and energy densities Omega(sub b)h(exp 2), Omega(sub c)h(exp 2)and Omega(sub Lambda), are each determined to a precision of approx. 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 5 sigma 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 Lambda-CDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their Lambda-CDM 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 Omega(sub kappa) = (0.0027 (sub +0.0039) (sup −0.0038;) the summed mass of neutrinos is limited to Sigma M(sub nu) < 0.44 eV (95% CL); and the number of relativistic species is found to lie within N(sub eff) = 3.84 +/- 0+/-40, when the full data are analyzed. The joint constraint on N(sub eff) and the primordial helium abundance, Y(sub He), 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.

cosmology↗
NASA NTRS2012

Scientific Verification of Faraday Rotation Modulators: Detection of Diffuse Polarized Galactic Emission

The design and performance of a wide bandwidth linear polarization-modulator based on the Faraday effect is described. Faraday Rotation Modulators (FRMs) are solid-state polarization switches that are capable of modulation up to approx 10 kHz. Six FRMs were utilized during the 2006 observing season in the Background Imaging of Cosmic Extragalactic Polarization (BICEP) experiment; three FRMs were used at each of BICEP fs 100 and 150 GHz frequency bands. The technology was verified through high signal-to-noise detection of Galactic polarization using two of the six FRMs during four observing runs in 2006. The features exhibit strong agreement with BICEP fs measurements of the Galaxy using non-FRM pixels and with the Galactic polarization models. This marks the first detection of high signal-to-noise mm-wave celestial polarization using fast, active optical modulation. The performance of the FRMs during periods when they were not modulated was also analyzed and compared to results from BICEP fs 43 pixels without FRMs.

Moyerman, S.↗
NASA NTRS2011

The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectrum

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg(exp 2) with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500 < l < 10,000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. At l = 3000, about half the power at 148 GHz comes from primary CMB after masking bright radio sources. The power from thermal and kinetic SZ is estimated to be Beta(sub 3000) is identical to 6.8 +/- 2.9 mu K (exp 2), where Beta (sub l) is identical to l(l + 1) C(sub l)/2pi. The IR Poisson power at 148 GHz is Bewta(sub 3000) 7.8 +/- 0.7 muK(exp 2) (C(sub l) = 5.5 +/- 0.5 nK(exp 2)), and a clustered IR component is required with Beta (sub 3000) = 4.6 +/- 0.9 muK(exp 2), assuming an analytic model for its power spectrum shape. At 218 GHz only about 15% of the power, approximately 27 mu K(exp 2), is CMB anisotropy at l = 3000. The remaining 85% is attributed to IR sources (approximately 50% Poisson and 35% clustered), with spectral index alpha = 3.69 +/- 0.14 for flux scaling as S(nu) varies as nu(sup alpha). We estimate primary cosmological parameters from the less contaminated 148 GHz spectrum, marginalizing over SZ and source power. The ACDM cosmological model is a good fit to the data (chi square/dof = 29/46), and ACDM parameters estimated from ACT+Wilkinson Microwave Anisotropy Probe (WMAP) are consistent with the seven-year WMAP limits, with scale invariant n(sub s) = 1 excluded at 99.7% confidence level (CL) (3 sigma). A model with no CMB lensing is disfavored at 2.8 sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6 sigma detection of primordial helium, with Y(sub p) = 0.313 +/- 0.044, and a 4 sigma detection of relativistic species, assumed to be neutrinos, with N(sub eff) = 5.3 +/- 1.3 (4.6 +/- 0.8 with BAO+H(sub 0) data). From the CMB alone the running of the spectral index is constrained to be d(sub s) / d ln k = -0,034 +/- 0,018, the limit on the tensor-to-scalar ratio is r < 0,25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension G(sub mu) < 1.6 x 10(exp -7) (95% CL),

Dunkley, J.↗
NASA NTRS2011

Mu-Spec - A High Performance Ultra-Compact Photon Counting spectrometer for Space Submillimeter Astronomy

We have designed and are testing elements of a fully integrated submillimeter spectrometer based on superconducting microstrip technology. The instrument can offer resolving power R approximately 1500, and its high frequency cutoff is set by the gap of available high performance superconductors. All functions of the spectrometer are integrated - light is coupled to the microstrip circuit with a planar antenna, the spectra discrimination is achieved using a synthetic grating, orders are separated using planar filter, and detected using photon counting MKID detector. This spectrometer promises to revolutionize submillimeter spectroscopy from space. It replaces instruments with the scale of 1m with a spectrometer on a 10 cm Si wafer. The reduction in mass and volume promises a much higher performance system within available resource in a space mission. We will describe the system and the performance of the components that have been fabricated and tested.

Moseley, H.↗
NASA NTRS2011

Interpretation of the Arcade 2 Absolute Sky Brightness Measurement

We use absolutely calibrated data between 3 and 90 GHz from the 2006 balloon flight of the ARCADE 2 instrument, along with previous measurements at other frequencies to constrain models of extragalactic emission. Such emission is a combination of the cosmic microwave background (CMB) monopole, Galactic foreground emission, the integrated contribution of radio emission from external galaxies, any spectral distortions present in the CMB, and any other extragalactic source. After removal of estimates of foreground emission from our own Galaxy, and an estimated contribution of external galaxies, we present fits to a combination of the flat-spectrum CMB and potential spectral distortions in the CMB. We find 217 upper limits to CMB spectral distortions of u < 6x10(exp -4) and [Y(sub ff)] < 1x10(exp -4). We also find a significant detection of a residual signal beyond that, which can be explained by the CMB plus the integrated radio emission from galaxies estimated from existing surveys. This residual signal may be due to an underestimated galactic foreground contribution, an unaccounted for contribution of a background of radio sources, or some combination of both. The residual signal is consistent with emission in the form of a power law with amplitUde 18.4 +/- 2.1 K at 0.31 GHz and a spectral index of -2.57 +/- 0.05.

Seiffert, M.↗
NASA NTRS2011

ARCADE 2 Measurement of the Absolute Sky Brightness at 3-90 GHz

The ARCADE 2 instrument has measured the absolute temperature of the sky at frequencies 3, 8, 10, 30, and 90 GHz, uSing an open-aperture cryogenic instrument observing al balloon altitudes with no emissive windows between the beam-forming optics and the sky. An external blackbody calibrator provides an in situ reference. Systematic errors were greatly reduced by using differential radiometers and cooling all critical components to physical temperatures approximating the cosmic microwave background (CMB) temperature. A linear model is used to compare the output of each radiometer to a set of thermometers on the instrument. Small correction. are made for the residual emission from the flight train, balloon, atmosphere, and foreground Galactic emission. The ARCADE 2 data alone show an excess radio rise of 54 +/- 6 mK at 3.3 GHz in addition to a CMB temperature of 2.731 +/- 0.004 K. Combining the ARCADE 2 data with data from the literature shows an excess power-law spectrum of T = 24.1 +/- 2.1 (K)(v/v(sub o)(exp -2.599+/-0.036 from 22 MHz to 10 GHz (v(sub 0) = 310 MHz) in addition to a CMB temperature of 2.725 +/- 0.001 K.

Fixsen, D. J.↗
NASA NTRS2011

The Extragalactic Radio Background

The existence of an isotropic component of the high-latitude radio sky has been recognized for nearly fifty years, but has typically been assumed to be Galactic in origin. We use recent radio observations to test whether the observed high-latitude component could originate within either an extended Galactic halo or a more local "bubble" structure. The lack of significant polarization from the isotropic component, combined with the lack of significant correlation with the Galactic far-infrared emission, rule out an origin within the Galaxy. We conclude that an extragalactic origin is the only viable alternative for the bulk of the isotropic high-latitude emission. The extragalactic component is 2-3 times brighter than local (Galactic) emission towards the Galactic poles and is consistent with a power law in frequency with amplitude T(sub r) = 24.1 plus or minus 2.1 K and spectral index beta = -2.599 plus or minus 0.036 evaluated at reference frequency 310 MHz.

Kogut, A.↗
NASA NTRS2011

Measuring the Optical Properties of Astrophysical Dust Analogues: Instrumentation and Methods

Dust is found throughout the universe and plays an important role for a wide range of astrophysical phenomena. In recent years, new infrared facilities have provided powerful new data for understanding these phenomena. However, interpretation of these data is often complicated by a lack of complementary information about the optical properties of astronomically relevant materials. The Optical Properties of Astronomical Silicates with Infrared Techniques (OPASI-T) program at NASA's Goddard Space Flight Center is designed to provide new high-quality laboratory data from which we can derive the optical properties of astrophysical dust analogues. This program makes use of multiple instruments, including new equipment designed and built specifically for this purpose. The suite of instruments allows us to derive optical properties over a wide wavelength range, from the near-infrared through the millimeter, also providing the capability for exploring how these properties depend upon the temperature of the sample. In this paper, we discuss the overall structure of the research program, describe the new instruments that have been developed to meet the science goals, and demonstrate the efficacy of these tools.

Rinehart, S. A.↗
NASA NTRS2010

ARCADE 2 Observations of Galactic Radio Emission

We use absolutely calibrated data from the Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission (ARCADE 2) flight in July 2006 to model Galactic emission at frequencies 3, 8, and 10 GHz. The spatial structure in the data is consistent with a superposition of free-free and synchrotron emission. Emission with spatial morphology traced by the Haslam 408 MHz survey has spectral index beta_synch = -2.5 +/- 0.1, with free-free emission contributing 0.10 +/- 0.01 of the total Galactic plane emission in the lowest ARCADE 2 band at 3.15 GHz. We estimate the total Galactic emission toward the polar caps using either a simple plane-parallel model with csc|b| dependence or a model of high-latitude radio emission traced by the COBE/FIRAS map of CII emission. Both methods are consistent with a single power-law over the frequency range 22 MHz to 10 GHz, with total Galactic emission towards the north polar cap T_Gal = 0.498 +/- 0.028 K and spectral index beta = -2.55 +/- 0.03 at reference frequency 0.31 GHz. The well calibrated ARCADE 2 maps provide a new test for spinning dust emission, based on the integrated intensity of emission from the Galactic plane instead of cross-correlations with the thermal dust spatial morphology. The Galactic plane intensity measured by ARCADE 2 is fainter than predicted by models without spinning dust, and is consistent with spinning dust contributing 0.4 +/- 0.1 of the Galactic plane emission at 23 GHz.

Kogut, A.↗
NASA NTRS2010

Five-Year Wilkinson Microwave Anisotropy Probe Observations: Data Processing, Sky Maps, and Basic Results

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 Wilkinson Microwave Anisotropy Probe (WMAP) sky survey. The new maps are consistent with previous maps and are more sensitive. 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 W-band polarization data is not yet suitable for cosmological studies, but we suggest directions for further study. We do find that Ka-band data is suitable for use; in conjunction with the additional years of data, the addition of Ka band to the previously used Q- and V-band channels significantly reduces the uncertainty in the optical depth parameter, tau. Further scientific results from the five-year data analysis are presented in six companion papers and are summarized in Section 7 of this paper. With the five-year WMAP data, we detect no convincing deviations from the minimal six-parameter ACDM 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 in the galaxy distribution, we find (68% CL uncertainties): OMEGA(sub b)h(sup 2) = 0.02267(sup +0.00058)(sub -0.00059), OMEGA(sub c)h(sup 2) = 0.1131 plus or minus 0.0034, OMEGA(sub logical and) = 0.726 plus or minus 0.015, ns = .960 plus or minus 0.013, tau = 0.84 plus or minus 0.016, and DELTA(sup 2)(sub R) = (22.445 plus or minus 0.096) x 10(exp -9) at k = 0.002 Mpc(exp -1). From these we derive sigma(sub 8) = 0.812 plus or minus 0.026, H(sub 0) = 70.5 plus or minus 1.3 kilometers per second Mpc(exp -1), OMEGA(sub b) = 0.0456 plus or minus 0.0015, OMEGA(sub c) = .228 plus or minus 0.013, OMEGA(sub m)h(sup 2) = 0.1358(sup +0.0037)(sub -0.0036), z reion = 10.9 plus or minus 1.4, and t(sub 0) = 13.72 plus or minus 0.12 Gyr. The new limit on the tensor-to-scalar ration is r less than 0.22 (95% CL), while the evidence for a running spectral index is insignificant, dn(sub s)/d ln k = -0.028 plus or minus 0.020 (68% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and the spatial curvature of the universe: -0.14 less than 1 + w less than 0.12 (95% CL) and -0.0179 less than OMEGA(sub k) less than 0.0081 (95% CL). The number of relativistic degrees of freedom, expressed in units of the effective number of neutrino species, is found to be N(sub eff) = 4.4 plus or minus 1.5 (69% CL), consistent with the standard value of 3.04. Models with N(sub eff) = 0 are disfavored at greater than 99% confidence. Finally, new limits on physically motivated primordial non-Gaussianity parameters are -9 less than f(sup local)(sub NL) less than 111 (95% CL) and -151 less than f(sup equal)(sub NL) less than 253 (95% CL) for the local and equilateral models, respectively.

Hinshaw, G.↗
NASA NTRS2008

A Detector for Cosmic Microwave Background Polarimetry

We present preliminary design and development work on polarized detectors intended to enable Cosmic Microwave Background polarization measurements that will probe the first moments of the universe. The ultimate measurement will be challenging, requiring background-limited detectors and good control of systematic errors. Toward this end, we are integrating the beam control of HE-11 feedhorns with the sensitivity of transition-edge sensors. The coupling between these two devices is achieved via waveguide probe antennas and superconducting microstrip lines. This implementation allows band-pass filters to be incorporated on the detector chip. We believe that a large collection of single-mode polarized detectors will eventually be required for the reliable detection of the weak polarized signature that is expected to result from gravitational waves produced by cosmic inflation. This focal plane prototype is an important step along the path to this detection, resulting in a capability that will enable various future high performance instrument concepts.

Wollack, E.↗
NASA NTRS2008

The Wilkinson Microwave Anisotropy Probe (WMAP) Source Catalog

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 greater than 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:l range between the minimum and maximum fluxes.

Wright, E.L.↗
NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP1) Observations: Galactic Foreground Emission

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 X(sup 2) (sub v) 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.

Gold, B.↗
NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP)Observations: Beam Maps and Window Functions

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 approximately 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 approximately 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 approximately 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 approximately 0.5%.

Hill, R.S.↗
NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP Data

The Wilkinson Microwave Anisotropy Probe (WMAP), launched in 2001, has mapped out the Cosmic Microwave Background with unprecedented accuracy over the whole sky. Its observations have led to the establishment of a simple concordance cosmological model for the contents and evolution of the universe, consistent with virtually all other astronomical measurements. The WMAP first-year and three-year data have allowed us to place strong constraints on the parameters describing the ACDM model. a flat universe filled with baryons, cold dark matter, neutrinos. and a cosmological constant. with initial fluctuations described by nearly scale-invariant power law fluctuations, as well as placing limits on extensions to this simple model (Spergel et al. 2003. 2007). With all-sky measurements of the polarization anisotropy (Kogut et al. 2003; Page et al. 2007), two orders of magnitude smaller than the intensity fluctuations. WMAP has not only given us an additional picture of the universe as it transitioned from ionized to neutral at redshift z approx.1100. but also an observation of the later reionization of the universe by the first stars. In this paper we present cosmological constraints from WMAP alone. for both the ACDM model and a set of possible extensions. We also consider tlle consistency of WMAP constraints with other recent astronomical observations. This is one of seven five-year WMAP papers. Hinshaw et al. (2008) describe the data processing and basic results. Hill et al. (2008) present new beam models arid window functions, Gold et al. (2008) describe the emission from Galactic foregrounds, and Wright et al. (2008) the emission from extra-Galactic point sources. The angular power spectra are described in Nolta et al. (2008), and Komatsu et al. (2008) present and interpret cosmological constraints based on combining WMAP with other data. WMAP observations are used to produce full-sky maps of the CMB in five frequency bands centered at 23, 33, 41, 61, and 94 GHz (Hinshaw et al. 2008). With five years of data, we are now able to place better limits on the ACDM model. as well as to move beyond it to test the composition of the universe. details of reionization. sub-dominant components, characteristics of inflation, and primordial fluctuations. We have more than doubled the amount of polarized data used for cosmological analysis. allowing a better measure of the large-scale E-mode signal (Nolta et al. 2008). To this end we describe an alternative way to remove Galactic foregrounds from low resolution polarization maps in which Galactic emission is marginalized over, providing a cross-check of our results. With longer integration we also better probe the second and third acoustic peaks in the temperature angular power spectrum, and have many more year-to-year difference maps available for cross-checking systematic effects (Hinshaw et al. 2008).

Dunkey, J.↗