
Inflation is the leading theory of the first instant of the universe.
Inflation, which postulates that the universe underwent a period of rapid
expansion an instant after its birth, provides convincing explanation for
cosmological observations. Recent advancements in detector technology have
opened opportunities to explore primordial gravitational waves generated by the
inflation through Bmode (divergentfree) polarization pattern embedded in the
Cosmic Microwave Background anisotropies. If detected, these signals would
provide strong evidence for inflation, point to the correct model for
inflation, and open a window to physics at ultrahigh energies.
LiteBIRD is a satellite mission with a goal of detecting
degreeandlargerangularscale Bmode polarization. LiteBIRD will observe at
the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors.
It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to
measure and subtract foregrounds.
The U.S. LiteBIRD team is proposing to deliver subKelvin instruments that
include detectors and readout electronics. A lensletcoupled sinuous antenna
array will cover lowfrequency bands (40 GHz to 235 GHz) with four frequency
arrangements of trichroic pixels. An orthomodetransducercoupled corrugated
horn array will cover highfrequency bands (280 GHz to 402 GHz) with three
types of single frequency detectors. The detectors will be made with Transition
Edge Sensor (TES) bolometers cooled to a 100 milliKelvin base temperature by
an adiabatic demagnetization refrigerator.The TES bolometers will be read out
using digital frequency multiplexing with Superconducting QUantum Interference
Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a
single SQUID amplidier.
We report on the subKelvin instrument design and ongoing developments for
the LiteBIRD mission.

Fluctuations in the intensity and polarization of the cosmic microwave
background (CMB) and the largescale distribution of matter in the universe
each contain clues about the nature of the earliest moments of time. The next
generation of CMB and largescale structure (LSS) experiments are poised to
test the leading paradigm for these earliest momentsthe theory of cosmic
inflationand to detect the imprints of the inflationary epoch, thereby
dramatically increasing our understanding of fundamental physics and the early
universe.
A future CMB experiment with sufficient angular resolution and frequency
coverage that surveys at least 1% of the sky to a depth of 1 uKarcmin can
deliver a constraint on the tensortoscalar ratio that will either result in a
5sigma measurement of the energy scale of inflation or rule out all
largefield inflation models, even in the presence of foregrounds and the
gravitational lensing Bmode signal. LSS experiments, particularly
spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will
complement the CMB effort by improving current constraints on running of the
spectral index by up to a factor of four, improving constraints on curvature by
a factor of ten, and providing nonGaussianity constraints that are competitive
with the current CMB bounds.

This is a report on the status and prospects of the quantification of
neutrino properties through the cosmological neutrino background for the Cosmic
Frontier of the Division of Particles and Fields Community Summer Study
longterm planning exercise. Experiments planned and underway are prepared to
study the cosmological neutrino background in detail via its influence on
distanceredshift relations and the growth of structure. The program for the
next decade described in this document, including upcoming spectroscopic galaxy
surveys eBOSS and DESI and a new StageIV CMB polarization experiment CMBS4,
will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass
measurement will produce a high significance detection of nonzero sum m_nu,
whose lower bound derived from atmospheric and solar neutrino oscillation data
is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this
measurement will definitively rule out the inverted neutrino mass hierarchy,
shedding light on one of the most puzzling aspects of the Standard Model of
particle physics  the origin of mass. This precise a measurement of N_eff
will allow for high sensitivity to any light and dark degrees of freedom
produced in the big bang and a precision test of the standard cosmological
model prediction that N_eff = 3.046.

Weak gravitational lensing is a valuable probe of galaxy formation and
cosmology. Here we quantify the effects of using photometric redshifts
(photoz) in galaxygalaxy lensing, for both sources and lenses, both for the
immediate goal of using galaxies with photoz as lenses in the Sloan Digital
Sky Survey (SDSS) and as a demonstration of methodology for large, upcoming
weak lensing surveys that will by necessity be dominated by lens samples with
photoz. We calculate the bias in the lensing mass calibration as well as
consequences for absolute magnitude (i.e., kcorrections) and stellar mass
estimates, for a large sample of SDSS Data Release 8 (DR8) galaxies. The
redshifts are obtained with the template based photoz code ZEBRA on the SDSS
DR8 ugriz photometry. We assemble and characterise the calibration samples (~9k
spectroscopic redshifts from four surveys) to obtain photometric redshift
errors and lensing biases corresponding to our full SDSS DR8 lens and source
catalogues. Our tests of the calibration sample also highlight the impact of
observing conditions in the imaging survey when the spectroscopic calibration
covers a small fraction of its footprint; atypical imaging conditions in
calibration fields can lead to incorrect conclusions regarding the photoz of
the full survey. For the SDSS DR8 catalogue, we find sigma_z/(1+z)=0.096 and
0.113 for the lens and source catalogues, with flux limits of r=21 and r=21.8,
respectively. We also explore the systematic uncertainty in the lensing signal
calibration when using source photoz, and both lens and source photoz; given
the size of existing training samples, we can constrain the lensing signal
calibration (and therefore the normalization of the surface mass density) to
within 2 and 4 per cent, respectively. [ABRIDGED]

BigBOSS is a Stage IV groundbased dark energy experiment to study baryon
acoustic oscillations (BAO) and the growth of structure with a widearea galaxy
and quasar redshift survey over 14,000 square degrees. It has been
conditionally accepted by NOAO in response to a call for major new
instrumentation and a highimpact science program for the 4m Mayall telescope
at Kitt Peak. The BigBOSS instrument is a roboticallyactuated, fiberfed
spectrograph capable of taking 5000 simultaneous spectra over a wavelength
range from 340 nm to 1060 nm, with a resolution R = 30004800.
Using data from imaging surveys that are already underway, spectroscopic
targets are selected that trace the underlying dark matter distribution. In
particular, targets include luminous red galaxies (LRGs) up to z = 1.0,
extending the BOSS LRG survey in both redshift and survey area. To probe the
universe out to even higher redshift, BigBOSS will target bright [OII] emission
line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are
obtained to measure the BAO feature, trace the matter power spectrum at smaller
scales, and detect redshift space distortions. BigBOSS will provide additional
constraints on early dark energy and on the curvature of the universe by
measuring the Lyalpha forest in the spectra of over 600,000 2.2 < z < 3.5
quasars.
BigBOSS galaxy BAO measurements combined with an analysis of the broadband
power, including the Lyalpha forest in BigBOSS quasar spectra, achieves a FOM
of 395 with Planck plus Stage III priors. This FOM is based on conservative
assumptions for the analysis of broad band power (kmax = 0.15), and could grow
to over 600 if current work allows us to push the analysis to higher wave
numbers (kmax = 0.3). BigBOSS will also place constraints on theories of
modified gravity and inflation, and will measure the sum of neutrino masses to
0.024 eV accuracy.

Galaxy surveys are one of the most powerful means to extract the cosmological
information and for a given volume the attainable precision is determined by
the galaxy shot noise sigma_n^2 relative to the power spectrum P. It is
generally assumed that shot noise is white and given by the inverse of the
number density n. In this paper we argue one may be able to considerably
improve upon this: in the halo picture of cosmological structure all of the
dark matter is in halos of varying mass and galaxies are formed inside these
halos, but for the dark matter mass and momentum conservation guarantee that
nonlinear effects cannot develop a white noise in the dark matter power
spectrum on large scales. This suggests that with a suitable weighting a
similar effect may be achieved for galaxies, suppressing their shot noise. We
explore this idea with Nbody simulations by weighting central halo galaxies by
halo mass and find that the resulting shot noise can be reduced dramatically
relative to expectations, with a 1030 suppression at the highest number
density of n=4*10^3 (Mpc/h)^3 resolved in our simulations. For specific
applications other weighting schemes may achieve even better results and for
n=3*10^4(Mpc/h)^3 we can reduce sigma_n^2/P by up to a factor of 10 relative
to uniform weighting. These results open up new opportunities to extract
cosmological information in galaxy surveys, such as the recently proposed
multitracer approach to cancel sampling variance, and may have important
consequences for the planning of future redshift surveys. Taking full advantage
of these findings may require better understanding of galaxy formation process
to develop accurate tracers of the halo mass.

Accurate photometric redshifts are among the key requirements for precision
weak lensing measurements. Both the large size of the Sloan Digital Sky Survey
(SDSS) and the existence of large spectroscopic redshift samples that are
fluxlimited beyond its depth have made it the optimal data source for
developing methods to properly calibrate photometric redshifts for lensing.
Here, we focus on galaxygalaxy lensing in a survey with spectroscopic lens
redshifts, as in the SDSS. We develop statistics that quantify the effect of
source redshift errors on the lensing calibration and on the weighting scheme,
and show how they can be used in the presence of redshift failure and sampling
variance. We then demonstrate their use with 2838 source galaxies with
spectroscopy from DEEP2 and zCOSMOS, evaluating several public photometric
redshift algorithms, in two cases including a full p(z) for each object, and
find lensing calibration biases as low as 1% (due to fortuitous cancellation of
two types of bias) or as high as 20% for methods in active use (despite the
small mean photoz bias of these algorithms). Our work demonstrates that
lensingspecific statistics must be used to reliably calibrate the lensing
signal, due to asymmetric effects of (frequently nonGaussian) photoz errors.
We also demonstrate that largescale structure (LSS) can strongly impact the
photoz calibration and its error estimation, due to a correlation between the
LSS and the photoz errors, and argue that at least two independent degreescale
spectroscopic samples are needed to suppress its effects. Given the size of our
spectroscopic sample, we can reduce the galaxygalaxy lensing calibration error
well below current SDSS statistical errors.

We present the 3D real space clustering power spectrum of a sample of
\~600,000 luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey
(SDSS), using photometric redshifts. This sample of galaxies ranges from
redshift z=0.2 to 0.6 over 3,528 deg^2 of the sky, probing a volume of 1.5
(Gpc/h)^3, making it the largest volume ever used for galaxy clustering
measurements. We measure the angular clustering power spectrum in eight
redshift slices and combine these into a high precision 3D real space power
spectrum from k=0.005 (h/Mpc) to k=1 (h/Mpc). We detect power on gigaparsec
scales, beyond the turnover in the matter power spectrum, on scales
significantly larger than those accessible to current spectroscopic redshift
surveys. We also find evidence for baryonic oscillations, both in the power
spectrum, as well as in fits to the baryon density, at a 2.5 sigma confidence
level. The statistical power of these data to constrain cosmology is ~1.7 times
better than previous clustering analyses. Varying the matter density and baryon
fraction, we find \Omega_M = 0.30 \pm 0.03, and \Omega_b/\Omega_M = 0.18 \pm
0.04, The detection of baryonic oscillations also allows us to measure the
comoving distance to z=0.5; we find a best fit distance of 1.73 \pm 0.12 Gpc,
corresponding to a 6.5% error on the distance. These results demonstrate the
ability to make precise clustering measurements with photometric surveys
(abridged).

We analyze the SDSS Lyalpha forest P_F(k,z) measurement to determine the
linear theory power spectrum. Our analysis is based on fully hydrodynamic
simulations, extended using hydroPM simulations. We account for the effect of
absorbers with damping wings, which leads to an increase in the slope of the
linear power spectrum. We break the degeneracy between the mean level of
absorption and the linear power spectrum without significant use of external
constraints. We infer linear theory power spectrum amplitude
Delta^2_L(k_p=0.009s/km,z_p=3.0)=0.452_{0.0570.116}^{+0.069+0.141} and slope
n_eff=2.321_{0.0470.102}^{+0.055+0.131} (possible systematic errors are
included through nuisance parameters in the fit  a factor >~5 smaller errors
would be obtained on both parameters if we ignored modeling uncertainties). The
errors are correlated and not perfectly Gaussian, so we provide a chi^2 table
to accurately describe the results. The result corresponds to sigma_8=0.85,
n=0.94, for a LCDM model with Omega_m=0.3, Omega_b=0.04, and h=0.7, but is most
useful in a combined fit with the CMB. The inferred curvature of the linear
power spectrum and the evolution of its amplitude and slope with redshift are
consistent with expectations for LCDM models, with the evolution of the slope,
in particular, being tightly constrained. We use this information to constrain
systematic contamination, e.g., fluctuations in the UV background. This paper
should serve as a starting point for more work to refine the analysis,
including technical improvements such as increasing the size and number of the
hydrodynamic simulations, and improvements in the treatment of the various
forms of feedback from galaxies and quasars.

We combine the measurements of luminosity dependence of bias with the
luminosity dependent weak lensing analysis of dark matter around galaxies to
derive the galaxy bias and constrain nonlinear mass and cosmological
parameters. We take advantage of theoretical and simulation predictions that
predict that while halo bias is rapidly increasing with mass for high masses,
it is nearly constant in low mass halos. We use a new weak lensing analysis
around the same SDSS galaxies to determine their halo mass probability
distribution. These halo mass probability distributions are used to predict the
bias for each luminosity subsample and we find an excellent agreement with
observed values. The required nonlinear mass suggests slightly lower matter
density than usually assumed, Omegam=0.25+/ 0.03 for the simplest models. We
combine the bias constraints with those from the WMAP and the SDSS power
spectrum analysis to derive new constraints on bias and sigma_8. For the most
general parameter space we find sigma_8=0.88+/ 0.06 and b_*=0.99+/ 0.07. In
the context of spatially flat models we improve the limit on the neutrino mass
for the case of 3 degenerate families from m_nu<0.6eV without bias to
m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running
is allowed. The corresponding limit for 3 massless + 1 massive neutrino is
1.37eV.

We combine the constraints from the recent Lyalpha forest and bias analysis
of the SDSS with previous constraints from SDSS galaxy clustering, the latest
supernovae, and WMAP . Combining WMAP and the lya we find for the primordial
slope n_s=0.98\pm 0.02. We see no evidence of running, dn/d\ln k=0.003\pm
0.010, a factor of 3 improvement over previous constraints. We also find no
evidence of tensors, r<0.36 (95% c.l.). A positive correlation between tensors
and primordial slope disfavors chaotic inflation type models with steep slopes:
V \propto \phi^4 is outside the 3sigma contour. For the amplitude we find
sigma_8=0.90\pm 0.03 from the lyaf and WMAP alone. We find no evidence of
neutrino mass: for the case of 3 massive neutrino families with an inflationary
prior, \sum m_{\nu}<0.42eV and the mass of lightest neutrino is m_1<0.13eV at
95% c.l. For the 3 massless + 1 massive neutrino case we find m_{\nu}<0.79eV
for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions
compatible with the LSND results. We explore dark energy constraints in models
with a fairly general time dependence of dark energy equation of state, finding
Omega_lambda=0.72\pm 0.02, w(z=0.3)=0.98^{+0.10}_{0.12}, the latter changing
to w(z=0.3)=0.92^{+0.09}_{0.10} if tensors are allowed. We find no evidence
for variation of the equation of state with redshift,
w(z=1)=1.03^{+0.21}_{0.28}. These results rely on the current understanding
of the lyaf and other probes, which need to be explored further both
observationally and theoretically, but extensive tests reveal no evidence of
inconsistency among different data sets used here.

We explore several physical effects on the power spectrum of the Lymanalpha
forest transmitted flux. The effects we investigate here are usually not part
of hydrodynamic simulations and so need to be estimated separately. The most
important effect is that of high column density absorbers with damping wings,
which add power on large scales. We compute their effect using the
observational constraints on their abundance as a function of column density.
Ignoring their effect leads to an underestimation of the slope of the linear
theory power spectrum. The second effect we investigate is that of fluctuations
in the ionizing radiation field. For this purpose we use a very large high
resolution Nbody simulation, which allows us to simulate both the fluctuations
in the ionizing radiation and the small scale LyaF within the same simulation.
We find an enhancement of power on large scales for quasars and a suppression
for galaxies. The strength of the effect rapidly increases with increasing
redshift, allowing it to be uniquely identified in cases where it is
significant. We develop templates which can be used to search for this effect
as a function of quasar lifetime, quasar luminosity function, and attenuation
length. Finally, we explore the effects of galactic winds using hydrodynamic
simulations. We find the wind effects on the LyaF power spectrum to be be
degenerate with parameters related to the temperature of the gas that are
already marginalized over in cosmological fits. While more work is needed to
conclusively exclude all possible systematic errors, our results suggest that,
in the context of data analysis procedures where parameters of the LyaF model
are properly marginalized over, the flux power spectrum is a reliable tracer of
cosmological information.

We measure cosmological parameters using the threedimensional power spectrum
P(k) from over 200,000 galaxies in the Sloan Digital Sky Survey (SDSS) in
combination with WMAP and other data. Our results are consistent with a
``vanilla'' flat adiabatic LambdaCDM model without tilt (n=1), running tilt,
tensor modes or massive neutrinos. Adding SDSS information more than halves the
WMAPonly error bars on some parameters, tightening 1 sigma constraints on the
Hubble parameter from h~0.74+0.180.07 to h~0.70+0.040.03, on the matter
density from Omega_m~0.25+/0.10 to Omega_m~0.30+/0.04 (1 sigma) and on
neutrino masses from <11 eV to <0.6 eV (95%). SDSS helps even more when
dropping prior assumptions about curvature, neutrinos, tensor modes and the
equation of state. Our results are in substantial agreement with the joint
analysis of WMAP and the 2dF Galaxy Redshift Survey, which is an impressive
consistency check with independent redshift survey data and analysis
techniques. In this paper, we place particular emphasis on clarifying the
physical origin of the constraints, i.e., what we do and do not know when using
different data sets and prior assumptions. For instance, dropping the
assumption that space is perfectly flat, the WMAPonly constraint on the
measured age of the Universe tightens from t0~16.3+2.31.8 Gyr to
t0~14.1+1.00.9 Gyr by adding SDSS and SN Ia data. Including tensors, running
tilt, neutrino mass and equation of state in the list of free parameters, many
constraints are still quite weak, but future cosmological measurements from
SDSS and other sources should allow these to be substantially tightened.

We study the variation of the dark matter mass fraction of elliptical
galaxies as a function of their luminosity, stellar mass, and size using a
sample of 29,469 elliptical galaxies culled from the Sloan Digital Sky Survey.
We model ellipticals as a stellar Hernquist profile embedded in an
adiabatically compressed dark matter halo. This model allows us to estimate a
dynamical mass ($M_{dynm}$) at the halflight radius from the velocity
dispersion of the spectra, and to compare these to the stellar mass estimates
($M_{*}$) from Kauffmann et al (2003). We find that $M_{*}/L$ is independent of
luminosity, while $M_{dynm}/L$ increases with luminosity, implying that the
dark matter fraction increases with luminosity. We also observe that at a fixed
luminosity or stellar mass, the dark matter fraction increases with increasing
galaxy size or, equivalently, increases with decreasing surface brightness:
high surface brightness galaxies show almost no evidence for dark matter, while
in low surface brightness galaxies, the dark matter exceeds the stellar mass at
the half light radius. We relate this to the fundamental plane of elliptical
galaxies, suggesting that the tilt of this plane from simple virial predictions
is due to the dark matter in galaxies. We find that a simple model where
galaxies are embedded in dark matter halos and have a star formation efficiency
independent of their surface brightness explains these trends. We estimate the
virial mass of ellipticals as being approximately 730 times their stellar
mass, with the lower limit suggesting almost all of the gas within the virial
radius is converted into stars.

We investigate the promise of the Lyalpha forest for high precision
cosmology in the era of the Sloan Digital Sky Survey using low order Npoint
statistics. We show that with the existing data one can determine the
amplitude, slope and curvature of the slope of the matter power spectrum with a
few percent precision. Higher order statistics such as the bispectrum provide
independent information that can confirm and improve upon the statistical
precision from the power spectrum alone. The achievable precision is comparable
to that from the cosmic microwave background with upcoming satellites, and
complements it by measuring the power spectrum amplitude and shape at smaller
scales. Since the data cover the redshift range 2<z<4, one can also extract the
evolution of the growth factor and Hubble parameter over this range, and
provide useful constraints on the presence of dark energy at z>2.

We compare three independent, cosmological linear perturbation theory codes
to asses the level of agreement between them and to improve upon it by
investigating the sources of discrepancy. By eliminating the major sources of
numerical instability the final level of agreement between the codes was
improved by an order of magnitude. The relative error is now below 0.1% for the
dark matter power spectrum. For the cosmic microwave background anisotropies
the agreement is below the sampling variance up to l=3000, with close to 0.1%
accuracy reached over most of this range of scales. The same level of agreement
is also achieved for the polarization spectrum and the temperaturepolarization
crossspectrum. Linear perturbation theory codes are thus well prepared for the
present and upcoming high precision cosmological observations.

We present a survey of the cosmological applications of the next generation
of weak lensing surveys, paying special attention to the computational
challenges presented by the number of galaxies, $N_{gal} ~$ 10$^{5}$. We focus
on optimal methods with no pixelization and derive a multigrid $P^3M$ algorithm
that performs the relevant computations in $O(N_{gal} \log N_{gal})$ time. We
test the algorithm by studying three applications of weak lensing surveys 
convergence map reconstruction, cluster detection and $E$ and $B$ power
spectrum estimation using realistic 1 deg^{2} simulations derived from Nbody
simulations. The map reconstruction is able to reconstruct large scale features
without artifacts. Detecting clusters using only weak lensing is difficult
because of line of sight contamination and noise, with low completeness if one
desires low contamination of the sample. A power spectrum analysis of the
convergence field is more promising and we are able to reconstruct the
convergence spectrum with no loss of information down to the smallest scales.
The numerical methods used here can be applied to other data sets with same
$O(N\log N)$ scaling and can be generalised to a sphere.

Observations of galaxygalaxy lensing from Sloan Digital Sky Survey (SDSS)
are combined with the TullyFisher and fundamental plane relations to derive
constraints on galactic halo profiles.We show that both for early and late type
galaxies around L* the rotation velocity drops significantly from its peak
value at the optical radius to the virial radius, v_opt/v_200 \sim 1.8 with
about 20% uncertainty. Such a drop is expected in models in which the halo
profile is very concentrated, so that it declines steeper than isothermal at
large radii. This large drop can be explained as a result of both a
concentrated dark matter profile and a significant stellar contribution to the
rotation velocity at the optical radii. We model the stellar component with a
thin rotationally supported disk or a Hernquist profile and use adiabatic dark
matter response model to place limits on the halo concentration as a function
of the stellar mass to light ratio. For reasonable values of the latter we find
concentrations consistent with CDM predictions, suggesting there is no evidence
for low concentrations for the majority of halos in the universe. We also
discuss the origin of FaberJackson relation L \propto sigma^4 in light of L
\propto v_200^2.5 relation found for early type galaxies from galaxygalaxy
lensing. This leads to a decrease in v_opt/v_{200} with luminosity above L*, so
that at 7L_* the ratio is 1.4. This is expected from the fundamental plane
relation and implies that relations such as TullyFisher and FaberJackson are
not simply relations between the mass of dark matter halo and galaxy
luminosity, but are also significantly influenced by the baryonic effects on
the rotation velocity at optical radii.

We use a highaccuracy computational code to investigate the precision with
which cosmological parameters could be reconstructed by future cosmic microwave
background (CMB) experiments, in particular the two satellite missions MAP and
Planck Surveyor (COBRAS/SAMBA). We identify several parameter combinations that
could be determined with a few percent accuracy with MAP and the Planck
Surveyor, as well as some degeneracies among the parameters that cannot be
accurately resolved with the temperature data alone. These degeneracies can be
broken by other astronomical measurements. Polarization measurments can
significantly enhance the science return of both missions by allowing accurate
determination of cosmological parameters, by enabling the detection of gravity
waves and by probing the ionization history of the universe. We also address
the question of how gaussian the likelihood function is around the maximum and
whether gravitational lensing changes the constraints.

Galaxies modeled as singular isothermal ellipsoids with an axis ratio
distribution similar to the observed axis ratio distribution of E and S0
galaxies are statistically consistent with both the observed numbers of
twoimage and fourimage lenses and the inferred ellipticities of individual
lenses. However, no fourimage lens is well fit by the model (typical
$\chi^2/N_{dof} \sim 20$), the axis ratio of the model can be significantly
different from that of the observed lens galaxy, and the major axes of the
model and the galaxy may be slightly misaligned. We found that models with a
second, independent, external shear axis could fit the data well (typical
$\chi^2/N_{dof} \sim 1$), while adding the same number of extra parameters to
the radial mass distribution does not produce such a dramatic improvement in
the fit. An independent shear axis can be produced by misalignments between the
luminous galaxy and its dark matter halo, or by external shear perturbations
due to galaxies and clusters correlated with the primary lens or along the line
of sight. We estimate that the external shear perturbations have no significant
effect on the expected numbers of twoimage and fourimage lenses, but that
they can be important perturbations in individual lens models. However, the
amplitudes of the external shears required to produce the good fits are larger
than our estimates for typical external shear perturbations (1015% shear
instead of 13% shear) suggesting that the origin of the extra angular
structure must be intrinsic to the primary lens galaxy in most cases.