Magnetic fields and gas flows around circumnuclear starbursts

Abstract: Radio continuum observations of barred galaxies revealed strong magnetic fields of >= 50-100 muG in the circumnuclear starbursts. Such fields are dynamically important and give rise to magnetic stress that causes inflow of gas towards the center at a rate of several solar masses per year, possibly along the spiral field seen in radio polarization and as optical dust lanes. This may solve the long-standing question of how to feed active nuclei, and explain the relation between the bolometric luminosity of AGN nuclei and the star-formation rate of their hosts. The strong magnetic fields generated in young galaxies may serve as the link between star formation and accretion onto supermassive black holes. -- Magnetic fields of >= 160 muG strength were measured in the central region of the almost edge-on starburst galaxy NGC 253. Four filaments emerging from the inner disk delineate the boundaries of the central outflow cone of hot gas. Strong Faraday rotation of the polarized emission from the background disk indicates a large-scale helical field in the outflow walls.

\section{Introduction}

Magnetic fields are a major agent in the ISM and also control the
density and propagation of cosmic rays. The radio--infrared
correlation indicates that turbulent fields are strongest in
star-forming regions. Magnetic fields and cosmic rays can provide
the pressure to drive galactic outflows. Outflows from starburst
galaxies in the early Universe may have magnetized the intergalactic
medium. In spite of our increasing knowledge of cosmic magnetic
fields, many important questions, especially their origin, strength
in intergalactic space, first occurrence in young galaxies and their
dynamical importance for galaxy evolution remain unanswered.

\section{Origin of galactic magnetism}

The most promising mechanism to sustain magnetic fields in the
interstellar medium of galaxies is the {\em dynamo}. In young
galaxies a small-scale dynamo probably amplified seed fields from
the protogalactic phase to the energy density level of turbulence
within less than $10^9$~yr (Schleicher et al. \cite{schleicher10}).
To explain the generation of large-scale fields in galaxies, the
mean-field dynamo has been developed (Beck et al. \cite{beck96}). It
is based on turbulence, differential rotation and helical gas flows,
driven by supernova explosions (Gressel et al. \cite{gressel08}) and
cosmic rays (Hanasz et al. \cite{hanasz09}). The mean-field dynamo
in galaxy disks predicts that within a few $10^9$~yr large-scale
regular fields are excited from the seed fields (Arshakian et al.
\cite{arshakian09}), forming patterns (``modes'') with different
azimuthal symmetries in the disk and vertical symmetries in the
halo.

\section{Observation of galactic magnetism}

Most of what we know about galactic magnetic fields comes through
the detection of radio waves. The intensity of {\em synchrotron
emission}\ is a measure of the number density of cosmic-ray
electrons and of the strength of the total magnetic field component
in the sky plane. The assumption of energy equipartition between the
total cosmic rays and total magnetic fields allows us to calculate
the total magnetic field strength from the synchrotron intensity
(Beck \& Krause \cite{beck+krause05}).

Polarized emission originates from ordered fields. As polarization
``vectors'' are ambiguous by $180^\circ$, they cannot distinguish
{\em regular (coherent) fields}, defined to have a constant
direction within the telescope beam, from {\em anisotropic fields},
which are generated from turbulent fields by compressing or shearing
gas flows, so that their direction frequently reverses perpendicular
to the flow direction. Unpolarized synchrotron emission indicates
{\em turbulent fields}\, which have random directions.

The intrinsic degree of linear polarization of synchrotron emission
is about 75\%. The observed degree of polarization is smaller due to
the contribution of unpolarized thermal emission, by {\em Faraday
depolarization}\ along the line of sight and across the beam
(Sokoloff et al. \cite{sokoloff98}) and by geometrical
depolarization due to variations of the field orientation within the
beam.

At radio wavelengths of a few centimeters and below, the orientation
of the observed B--vector is parallel to the field orientation, so
that the magnetic patterns of many galaxies can be mapped directly
(e.g. Beck \cite{beck05}). At longer wavelengths, the observed
polarization vector is rotated in a magnetized thermal plasma by
{\em Faraday rotation}. The rotation angle increases with the square
of the wavelength $\lambda^2$ and with the {\em Rotation Measure
(RM)}, which is the integral of the plasma density and the strength
of the component of the field along the line of sight. As the
rotation angle is sensitive to the sign of the field direction, only
regular fields give rise to Faraday rotation, while anisotropic and
random fields do not.

\section{Magnetic fields in barred galaxies}

\begin{figure*}[t]
\vspace*{2mm}
\begin{minipage}[t]{8cm}
\begin{center}
\includegraphics[bb = 48 49 522 552,width=8cm,clip=]{n1097cm6.eps}
\caption{Total radio emission (contours) and $B$--vectors of the
barred galaxy NGC~1097, observed at 6~cm wavelength with the {\it
VLA}\ and smoothed to 10'' resolution (Beck et al. \cite{beck05a}).
The background optical image is from Halton Arp ({\it Cerro Tololo
Observatory}).} \label{fig:n1097}
\end{center}
\end{minipage}\hfill
\begin{minipage}[t]{7cm}
\begin{center}
\includegraphics[bb = 50 51 519 503,width=7cm,clip=]{n1097cm3_center.eps}
\caption{Total radio emission (colors) and $B$--vectors of the
circumnuclear ring of the barred galaxy NGC~1097, observed at 3~cm
wavelength with the {\it VLA}\ and smoothed to 3'' resolution (Beck
et al. \cite{beck05a}).} \label{fig:n1097_center}
\end{center}
\end{minipage}
\end{figure*}

The gravitational potential of strongly barred galaxies drives
noncircular orbits and gas inflow. In many barred galaxies a
circumnuclear ring is formed. Numerical models show that gas
streamlines are deflected in the bar region along shock fronts,
behind which the cold gas is compressed in a fast shearing flow
(Athanassoula \cite{athanassoula}). The compression regions traced
by massive dust lanes develop along the edge of the bar that is
leading with respect to the galaxy's rotation because the gas
rotates faster than the bar pattern.

20 galaxies with large bars were observed with the Very Large Array
(VLA) and with the Australia Telescope Compact Array (ATCA) (Beck et
al. \cite{beck02,beck05a}). The total radio luminosity (a measure of
the total magnetic field strength) is strongest in galaxies with
high far-infrared luminosity (indicating high star-formation
activity), a result similar to that in non-barred galaxies. The
average radio intensity, radio luminosity and star-formation
activity all correlate with the relative bar length. Polarized
emission was detected in 17 of the 20 barred galaxies.

NGC~1097 (Fig.~\ref{fig:n1097}) is one of the nearest barred
galaxies and hosts a huge bar of about 16\,kpc length. The total and
polarized radio intensities are strongest in the downstream region
of the dust lanes (southeast of the center) by compression of
turbulent fields in the bar's shock.
The pattern of field lines in NGC~1097 is similar to that of the gas
streamlines as obtained in numerical simulations (Athanassoula
\cite{athanassoula}). This suggests that the ordered magnetic field
is aligned with the flow and amplified by shear. Remarkably, the
optical image of NGC~1097 shows dust filaments in the upstream
region which are aligned with the ordered field
(Fig.~\ref{fig:n1097}). Between the region upstream of the southern
bar and the downstream region the field lines smoothly change their
orientation by almost $90^\circ$. The ordered field is also hardly
compressed, probably coupled to the diffuse gas and strong enough to
affect its flow (Beck et al. \cite{beck05a}). The polarization
pattern in barred galaxies can be used as a tracer of shearing gas
flows in the sky plane and complements spectroscopic measurements of
radial velocities.

\section{Magnetic fields in circumnuclear rings and gas inflows}

Barred galaxies can create a gas reservoir and starburst in the
central kiloparsec, fuelling an AGN needs other processes (e.g.
Davies et al. \cite{davies10}). The central regions of barred
galaxies are often sites of ongoing intense star formation and
strong magnetic fields. Bright radio emission has been observed in
several galaxies with the following equipartition strengths of the
total magnetic fields (Beck et al. \cite{beck05b}):

\begin{itemize}
\item NGC~1097: $\simeq50$\,$\mu$G in the circumnuclear ring,
50\,$\mu$G in the ring knots
\item NGC~1365: $\simeq60$\,$\mu$G in the dust lanes (no ring)
\item NGC~1672: $\simeq70$\,$\mu$G in the ring knots
\item NGC~7552: $\simeq105$\,$\mu$G in the ring
\end{itemize}

Note that these values are lower limits because the energy densities
of cosmic-ray protons and electrons in starburst regions are reduced
by various losses (Thompson et al. \cite{thompson06}). This is
supported by an investigation of the $\gamma$-ray and radio emission
from the starburst galaxies M~82 and NGC~253 which suggests that
most of the radio emission is from secondaries undergoing strong
bremsstrahlung and ionization losses (Lacki et al. \cite{lacki11}).

One of the key problems of AGN physics is how to fuel the nucleus
with gas from the central region, which needs outwards transport of
angular momentum (Davies et al. \cite{davies09}). Magnetic fields
can help here. The circumnuclear ring of NGC~1097 is bright in the
optical, IR and radio spectral ranges. It has a diameter of about
1.5\,kpc and an active nucleus in its center
(Fig.~\ref{fig:n1097_center}). The nonthermal and weakly polarized
radio emission is a signature of strong turbulent magnetic fields
(Beck et al. \cite{beck05a}). Magnetic stress in the differentially
rotating circumnuclear ring can drive mass inflow (Balbus \& Hawley
\cite{balbus98}) at a rate of:

\begin{eqnarray}
dM/dt = 2\pi \, \sigma \, T_{r\phi} \, / \, \Omega
\end{eqnarray}
where $\sigma$ is the surface mass density, $\Omega$ its angular
rotation velocity and $T_{r \, \phi}$ is the magnetic stress tensor.
$r$ and $\phi$ denote the radial and azimuthal field components. Its
dominant component can be written in terms of the Alfv{\'e}n
velocity $v_A$ as $T_{r\phi} \simeq - < v_{A,r} \, v_{A,\phi} >$.
This yields:

\begin{eqnarray}
dM/dt = - (<b_r \, b_{\phi}> + B_r \, B_{\phi}) \, h \, / \, \Omega
\end{eqnarray}
where $h$ is the scale height of the gas, $b$ the strength of the
turbulent field and $B$ that of the ordered field. The correlation
between $b_r$ and $b_{\phi}$ is generated by shear from differential
rotation. For NGC~1097, $h\simeq100$\,pc, $v\simeq450$\,km/s at
1\,kpc radius, $b_r \simeq b_\phi \simeq$\,50\,$\mu$G and $B_r
\simeq B_\phi \simeq$\,10\,$\mu$G give an inflow rate of several
$M_\odot$/yr, which is sufficient to fuel the activity of the
nucleus of this galaxy (Beck et al. \cite{beck05a}). This mechanism
can be expected to operate in many galaxies.

Equipartition between the energy densities of cosmic rays and
magnetic fields yields:

\begin{eqnarray}
dM/dt \simeq - B_{tot}^2 \, \, h / \, \Omega \, \, \propto -
I_{sync}^{0.5} \, h \, / \, \Omega \, \, \propto - I_{IR}^{0.5} \, h
\, / \, \Omega \, \, \propto - SFR^{0.5} \, h \, / \, \Omega
\end{eqnarray}
where $B_{tot}$ is the total field strength, $I_{sync}$ is the radio
synchrotron intensity which is proportional to the infrared
intensity $I_{IR}$ (the radio-infrared correlation) and hence to the
star-formation rate $SFR$. The relation between the bolometric
luminosity of AGNs $L_{bol}$ and the $SFR$ of their hosts (e.g.
Trakhtenbrot \& Netzer \cite{trakhtenbrot10}) can be explained by
this scenario. The relation between the accretion rate and SFR
(Eq.~3) should be tested with a sample of galaxies.

The ordered field in the ring of NGC~1097 has a spiral pattern and
extends towards the nucleus. Its pitch angle agrees with that of the
spiral filaments seen in dust and H$_2$ (Prieto et al.
\cite{prieto05}, Davies et al. \cite{davies09}), suggesting gas
inflow along the magnetic field lines.

Radio emission is a tracer of star formation also in distant
galaxies. The radio--infrared correlation holds to redshifts of at
least $z\simeq3$ (Murphy \cite{murphy09}), indicating strong
magnetic fields in young galaxies, probably generated by the
turbulent dynamo mechanism (Schleicher et al. \cite{schleicher10}).
Infrared observations indicate that star formation in the hosts of
the most luminous AGNs peaked at $z\simeq3$ and decreased at later
epochs (Serjeant et al. \cite{serjeant10}), which is consistent with
the evolution of the accretion rate onto supermassive black holes
(Trakhtenbrot et al. \cite{trakhtenbrot11}). The strong magnetic
fields in young galaxies may serve as the link between star
formation and accretion.

\section{The nuclear outflow of the starburst galaxy NGC~253}

\begin{figure*}[t]
\vspace*{2mm}
\begin{minipage}[t]{7cm}
\begin{center}
\includegraphics[bb = 19 19 593 570,width=7cm,clip=]{n253cm6.eps}
\caption{Total radio emission (contours and colors) and polarization
$B$--vectors of the spiral galaxy NGC~253, combined from
observations at 6~cm wavelength with the {\it VLA}\ and the {\it
Effelsberg}\ telescope and smoothed to 30'' resolution (Heesen et
al. \cite{heesen09}).} \label{fig:n253}
\end{center}
\end{minipage}\hfill
\begin{minipage}[t]{8cm}
\begin{center}
\includegraphics[bb = 30 30 198 173,width=8cm,clip=]{n253cm20.eps}
\caption{Total radio emission (greyscale) of the inner starburst
region of NGC~253, observed at 20~cm with the {\it VLA}\ at a
resolution of 1.3'' x 2.2'' (Heesen et al. \cite{heesen11}). The
radio filaments (F1--4) and the minimum (M1) are indicated. The
contours show the 12CO(2--1) emission, observed with the {\it SMA}
(Sakamoto et al. \cite{sakamoto06}).} \label{fig:n253_center}
\end{center}
\end{minipage}
\end{figure*}

NGC~253 is a prototypical starburst galaxy, hosting a starburst
nucleus (Brunthaler et al. \cite{brunthaler09}). Its high
inclination of $78.5^{\circ}$ allows observing the extraplanar
emission. The bright radio halo (Fig.~\ref{fig:n253}) is probably a
result of a galactic wind with a speed of about 300\,km/s (Heesen et
al. \cite{heesen09}).

Radio observations of the inner starburst region with high
resolution revealed four filaments with widths of less than 40\,pc
and lengths of at least 500\,pc (Fig.~\ref{fig:n253_center}). These
are interpreted as the boundaries of the outflow cone of hot gas
interacting with the halo gas (Heesen et al. \cite{heesen11}). The
equipartition strength of the total field in the filaments is at
least 40\,$\mu$G, that in the central region $\ge160\,\mu$G.

\begin{figure*}[t]
\vspace*{2mm}
\begin{minipage}[t]{8.5cm}
\begin{center}
\includegraphics[bb = 14 14 566 428,width=8cm,clip=]{n253cm20+chandra.eps}
\caption{Total radio emission (contours) of the inner starburst
region of the galaxy NGC~253, observed at 20~cm wavelength with the
{\it VLA}\ at a resolution of 1.3'' x 2.2'' (Heesen et al.
\cite{heesen11}). The colour image shows the X-ray emission
(0.5-5\,keV) measured with the {\it CHANDRA}\ satellite.}
\label{fig:n253_X}
\end{center}
\end{minipage}\hfill
\begin{minipage}[t]{6.5cm}
\begin{center}
\includegraphics[bb = 34 35 350 402,width=5.4cm,clip=]{n253_model.eps}
\caption{Model of the helical magnetic field in the walls of the
outflow cone (Heesen et al. \cite{heesen11}).}
\label{fig:n253_model}
\end{center}
\end{minipage}
\end{figure*}

Faraday rotation measures (RM) between 3\,cm and 6\,cm wavelengths
are consistent with an axisymmetric, even-parity, spiral magnetic
field in the disk at radii of more than 2\,kpc, while the field in
the inner disk is anisotropic (Heesen et al. \cite{heesen11}). Near
the SE minor axis, in the region of the X-ray outflow cone
(Fig.~\ref{fig:n253_X}), RM jumps between +300\,rad\,m$^{-2}$ and
-300\,rad\,m$^{-2}$. This is interpreted as the Faraday rotation of
the polarized emission from the background disk by the helical field
in the walls of the outflow extending to at least 1.2\,kpc from the
center (Fig.~\ref{fig:n253_model}). Only the inner part of this
helical field structure is seen in emission
(Fig.~\ref{fig:n253_center}). The northern outflow cone is located
behind the disk, so that its X-ray emission is mostly absorbed and
cannot Faraday-rotate the radio emission.

The galaxy NGC~3079 hosts an energetic wind-driven outflow which
extends 3\,kpc above and below the plane (Duric \& Seaquist
\cite{duric88}) and emits strong and highly polarized synchrotron
emission, in contrast to NGC~253. Radio-bright nuclear outflows seem
to be rare.

\section{Summary and outlook}

Magnetic fields in the central regions of galaxies are dynamically
important. Magnetic stress drives gas inflow from the circumnuclear
regions towards the nuclei. This process connects the star formation
activity to the accretion rate and could be important for the
formation of supermassive black holes. Ordered fields as observed by
polarized radio emission may trace the direction of gas in- and
outflows.

Results have been obtained so far only for radio-bright and nearby
galaxies. A better insight into the interaction between gas and
magnetic fields, especially near galactic nuclei, needs higher
resolution and higher sensitivity. The {\it EVLA}\ has largely
improved the sensitivity for radio continuum observations. The {\it
EVLA}\ project {\it CHANGES}\ (PI: Judith Irwin) will search for
outflows in a large sample of edge-on galaxies. The goal for the
next decade is the {\it Square Kilometre Array (SKA)}\ with a
collecting area of about 100x that of present-day telescopes. The
SKA will open a new era in the observation of cosmic magnetic fields
(Beck \cite{beck10}). Starburst galaxies will be observable out to
redshifts of $z\simeq3$ (Murphy \cite{murphy09}).

\section{References}


Robust narrow-gap semiconducting behavior in square-net La$_{3}$Cd$_{2}$As$_{6}$

Abstract: ABSTRACT: Narrow-gap semiconductors are sought-after materials due to their potential for long-wavelength detectors, thermoelectrics, and more recently non-trivial topology. Here we report the synthesis and characterization of a new family of narrow-gap semiconductors, $R$$_{3}$Cd$_{2}$As$_{6}$ ($R=$ La, Ce). Single crystal x-ray diffraction at room temperature reveals that the As square nets distort and Cd vacancies order in a monoclinic superstructure. A putative charge-density ordered state sets in at 279~K in La$_{3}$Cd$_{2}$As$_{6}$ and at 136~K in Ce$_{3}$Cd$_{2}$As$_{6}$ and is accompanied by a substantial increase in the electrical resistivity in both compounds. The resistivity of the La member increases by thirteen orders of magnitude on cooling, which points to a remarkably clean semiconducting ground state. Our results suggest that light square net materials within a $I4/mmm$ parent structure are promising clean narrow-gap semiconductors.

\section{INTRODUCTION}

Electronic instabilities are ubiquitous in quantum materials. Notable examples
include spin-density waves, superconductivity, and charge-density waves (CDWs).
The latter is a modulation of conduction electron density
accompanied by a periodic distortion of the crystal lattice \cite{ReviewCDW,ReviewCDW2,CDWlowD2}.
Interestingly, CDWs have been recently observed in close proximity
to exotic correlated phenomena, such as high-temperature superconductivity and axionic topological states \cite{CDW,cuprates2019,axions2019,axions2021}.

In addition, the formation of a CDW phase minimizes
the energy of the system by opening a gap
at the Fermi level, which may drive the material towards a
narrow-gap semiconducting state \cite{ReviewCDW,ReviewCDW2,CDWlowD2}. Disorder-free narrow-gap semiconducting materials are of interest due to their potential functionalities,
which include thermoelectricity, dark-matter detection, and nontrivial topology \cite{Thermo,DarkM,PbTeThermo}. Their narrow gap, however, often suffers from conducting in-gap states created by defects or impurities \cite{Book,Book2}.

One way to design narrow-gap semiconductors is
to search for low-dimensional crystal structures, which are fundamentally unstable to the formation
of energy gaps from small distortions \cite{CDWlowD2,ReviewCDW,ReviewCDW2}.
In fact, CDW phases were initially identified in one-dimensional (1D) transition metal
trichalcogenides $MX_{3}$ ($M =$ Nb and Ta, $X =$ S, Se, or Te) containing chains of NbSe$_{6}$ prisms
and in quasi-2D layered transition metal
dichalcogenides \cite{ReviewCDW,ReviewCDW2,CDW}.

Here we consider layered tetragonal materials with the ``112" general formula
$RMX_{2}$ ($R =$ lanthanide, $M =$ Au, Ag, Cu, Cd, and Zn, $X =$ Bi, Sb, and As), which contain pnictide square nets.
Previous band structure analysis and molecular orbital models have shown that square nets built from
more electronegative elements are prone to distortions that lift the tetragonal symmetry \cite{Sqrnets_0}.
This tendency occurs for two reasons. First, the $R$-square net $p$-bands close to the Fermi level ($E_{F}$) and the $M$ $d$ bands become well separated, which allows electron transfer and the formation of anionic units, that form covalent bonds leading to Zintl phases \cite{Zintl_1,Zintl_2,Zintl_Review,Zintl_Rosa}. Second, there is greater mixing between $s$ and $p$ orbitals in pnictide elements with small atomic number ($Z$), which (again) favors the classic
octet rule in a distorted structure over hypervalent bonds in an undistorted square net \cite{Sqrnets_1,Sqrnets_2}.

Important insights into distorted structures come from chalcogenide-based square nets, and numerous polytellurides host CDW phases with distorted Te square nets. $RE$Te$_{3}$ ($RE$ = rare-earth) \cite{ReTe3,SmTe3,CeTe3}, Pb$_{3-x}$Sb$_{1+x}$S$_{4}$Te$_{2 - \delta}$ \cite{PbSbSTe}, Sm$_{2}$Te$_{5}$ \cite{Sm2Te5}, Cu$_{0.63}$EuTe$_{2}$ \cite{CuEuTe}, $AMR$Te$_{4}$ ($A$ = K, Na and $M$ = Cu, Ag) \cite{AMRETe4}, and K$_{1/3}$Ba$_{2/3}$AgTe$_{2}$ \cite{KBaAgTe2} are a few examples. In the latter, cation sites are partially occupied in order to satisfy the Zintl rule in the distorted structure. Importantly, the primary cause of the observed vacancy superstructure is argued to be a distortion in the tellurium square net, rather than the partial or complete ordering of the cations \cite{KBaAgTe2}. Yet many distorted polytellurides exhibit metallic behavior \cite{RTe3metal,SmTemetal,CuEuTe,AMRETe4}. Recently, GdTe$_{3}$ has been shown to display a remarkably high electron mobility of 61,200~cm$^{2}$V$^{-1}$s$^{-1}$ at low temperatures \cite{GdTe3}. Ultrahigh mobilities ($>$ 10,000~cm$^{2}$V$^{-1}$s$^{-1}$) and an anomalous Hall effect are also observed in pnictide square net materials such as EuMnBi$_{2}$ \cite{EuMnBi2}.

\begin{figure}[!t]
\includegraphics[width=0.48\textwidth]{Fig01.pdf}
\caption{(a) Structural relationship of the tetragonal CeCdAs$_2$ parent ($I4/mmm$, red), the conceptual orthorhombic 3$\times$1$\times$3 vacancy superstructure (green) and the true, monoclinic, Ce$_3$Cd$_2$As$_6$ unit cell ($C2/m$, blue). (b) Detailed view of the As square-net layers. To emphasize the stripe distortion, only the short As6--As7 bonds are drawn. (c) Detailed view of the CdAs layers. For reference, the nominal undistorted structure is drawn as a shaded background. Cd ions are significantly displaced towards the vacancy stripes. As a consequence, the apical As ions As1 and As3 are slightly drawn towards the Cd layer, while the As2 ions dangling above/below the Cd-vacancy stripes are slightly pushed away from the Cd layer. }
\label{Fig1}
\end{figure}

Here we highlight a route for
clean narrow-gap semiconducting behavior in distorted arsenic square nets. We investigate a rather underexplored phase space in which $M=$ post-transition metal \cite{layeredAs1,layeredAs2}. The
substantial electronegativity of low-$Z$ arsenic combined with the decrease in metallic character of post-transition metals
provides a promising route to the realization of distortion-driven charge-density waves.

In this paper, we report the synthesis and characterization of distorted LaCd$_{2/3}$As$_{2}$ and CeCd$_{2/3}$As$_{2}$.
At room temperature, both compounds crystallize in a vacancy superstructure of the $I4/mmm$ ``112" parent
structure and show semiconducting behavior. The ordered superstructure thus corresponds to the ``326" stoichiometry. As temperature is decreased, both compounds display
a phase transition at $T_{\mathrm{CDW}} = $ 136~K and 278~K for Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$, respectively. Given the many structural degrees of freedom due to the broken tetragonal symmetry, we attribute these transitions to CDW-driven structural distortions. This picture also explains the opening of a gap in the electronic density of states. Remarkably, the electrical
resistivity of La$_{3}$Cd$_{2}$As$_{6}$ increases by thirteen orders of magnitude on cooling below $T_{\mathrm{CDW}}$, in agreement
with the vanishingly small Sommerfeld coefficient from
specific heat measurements. The estimated activation energies from Arrhenius plots are
105(1)~meV and 74(1)~meV for La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$, respectively.

 \begin{figure*}[!t]
	\includegraphics[width=\textwidth]{Fig2_lowercase.pdf}%
	\caption{\label{Fig2} Comparison of experimental and calculated x-ray diffraction intensity maps of reciprocal space. The left and right halves of each panel show the measured and calculated intensity, respectively. For convenience, the data is indexed in the $I4/mmm$ parent cell. (a,b) $(hk\frac23)$ slices reveal the presence of (a) two stripe domains (i.e., twins) in LaCd$_{0.67}$As$_2$ and (b) four in CeCd$_{0.67}$As$_2$. (c,d)  ($h0l$) slices show a large number of $q_\text{Cd}=\langle\frac23 0 \frac23\rangle$ superstructure peaks. The modulations of the intensities of these peaks is highly sensitive of small distortions of the ionic arrangement away from the nominal positions in the tetragonal parent.}
\end{figure*}

\section{EXPERIMENTAL DETAILS}

Black plate-like single crystals of Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$ were grown by the vapor transport technique.
First, a polycrystalline seed of 1La:0.7Cd:2As was prepared via solid state reaction at 800~$^{\circ}$C.
Then the polycrystalline powder was loaded along with iodine in a quartz tube, which was sealed in vacuum.
The tube was kept in a temperature gradient from 830~$^{\circ}$C to 720~$^{\circ}$C for a week. The initial polycrystalline material was kept in the hot zone, and single crystals precipitated in the cold zone. The synthesized phase was investigated by Mo $K_\alpha$ single crystal x-ray diffraction and energy-dispersive x-ray spectroscopy (EDX) both at room temperature. EDX measurements confirmed the 3:2:6 stoichiometry within experimental error. The specific heat of a collection of approximately 10 single crystals (total mass $\approx 1$~mg) was measured using a Quantum Design PPMS that employs a quasi-adiabatic thermal relaxation technique. The electrical resistivity ($\rho$) was characterized with the same instrument in the standard four-probe configuration.
The current was applied in the basal plane of the crystal. For Ce$_{3}$Cd$_{2}$As$_{6}$, an AC bridge was used to measure $\rho$ in all temperature range. While, for La$_{3}$Cd$_{2}$As$_{6}$, $\rho$ was measured with an AC bridge, at high temperatures, whereas at low temperatures a DC method was required due to the large resistance of the sample.

\section{RESULTS}

\begin{table}[b]
	\begin{center}
		\begin{tabular}{ l c c c c c }
			\multicolumn{6}{l}{\textbf{$\boldsymbol{R}$Cd$_{2/3}$As$_{2}$}, $I4/mmm$ (\#139), $Z=4$ }\\
			\multicolumn{6}{l}{$a_t=b_t\sim4.1\,$\AA, $c_t\sim21.3\,$\AA}\\ [5pt]
			\multicolumn{2}{l}{ion, Wyck.} & $x$ & $y$ & $z$& occ. (\%)\\
			\hline
			$R$  & $4e$ &    $0$     &    $0$     &  0.11         & 100 \\
			Cd  & $4d$ &    $0$     & $\sfrac12$ & $\sfrac14$    & 66  \\
			As1 & $4c$ &    $0$     & $\sfrac12$ &    $0$        & 100 \\
			As2 & $4e$ &    $0$     &    $0$     &  0.34         & 100 \\
			\hline
		\end{tabular}
		\caption{ \label{Tab01} Structural parameters of $R_{3}$Cd$_{2}$As$_{6}$, inferred from refinements in a tetragonal cell (neglecting superstructure reflections). A detailed summary of refined parameters and uncertainties for either compound ($R=$ La, Ce) is provided in the Supporting Information \cite{sm}.}
	\end{center}
\end{table}

We start with the structural characterization of Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$. For both materials, Bragg reflections can be indexed in a tetragonal unit cell corresponding to the respective parent compounds with ``112'' stoichiometry. The approximate parameters of this $I4/mmm$ structure are given in Table~\ref{Tab01}. The corresponding tetragonal cell is shaded red in Fig.~\ref{Fig1}(a). A detailed summary of the refined parameters of either compound is provided in the Supporting Information \cite{sm}. The occupancy of the Cd site ($4d$ in $I4/mmm$) converges close to 66\% in both materials, confirming the stoichiometry measured in EDX.

Interestingly, the Cd vacancies are not randomly distributed, but rather order in a stripe pattern. Fig.~\ref{Fig2} shows x-ray intensity maps interpolated to the $(h,k,\frac23)$ and $(h0l)$ planes of reciprocal space (with reference to the $I4/mmm$ cell). Aside from the integer-index Bragg peaks, a pattern of weak satellite reflections is observed, which is described by the propagation vector $\mathbf{q}_\mathrm{Cd}=\langle\frac23,0,\frac23\rangle$. This indicates an ordering of the vacancies in a 3$\times$1$\times$3 superstructure [green cell in \ref{Fig1}(a)], which can be reduced to a base-centered monoclinic structure of space group $C2/m$ [blue cell in \ref{Fig1}(a)].

\begin{table}[b]
 \begin{center}
	\begin{tabular}{llllllll}
		\multicolumn{8}{l}{\textbf{\textit{R}$_3$Cd$_2$As$_6$}, $C2/m$ (\#12, unique axis $b$, cell choice 1)}\\
		\hline
		\multicolumn{2}{l}{$\vec{a}=\vec{a}_t-\vec{c}_t$ } & ,~~ i.e.~~$a=\sqrt(a_t^2+c_t^2)$ & \multicolumn{5}{l}{ $\sim21.6$\,\AA  }    \\
		\multicolumn{2}{l}{$\vec{b}=-\vec{b}_t$ }          & ,~~ i.e.~~$b= a_t$              & \multicolumn{5}{l}{ $\sim4.1$\,\AA }       \\
		\multicolumn{2}{l}{$\vec{c}=-3\,\vec{a}_t$}        & ,~~ i.e.~~$c= 3\,a_t$            & \multicolumn{5}{l}{ $\sim12.2$\,\AA }   \\
		\multicolumn{8}{l}{$\alpha=\gamma=90^\circ$}\\
		\multicolumn{8}{l}{$\beta= 90^\circ+\arctan(a_t/c_t) \sim 100.8^\circ$ } \\[5pt]
	\end{tabular}
	\begin{tabular}{>{\centering\arraybackslash}lc>{\centering\arraybackslash}p{13mm}>{\centering\arraybackslash}p{7mm}>{\centering\arraybackslash}p{19mm}>{\centering\arraybackslash}p{12mm}>{\centering\arraybackslash}p{12mm}}
		&       & \multicolumn{3}{>{\centering\arraybackslash}p{39mm}}{nominal position} & \multicolumn{2}{>{\centering\arraybackslash}p{24mm}}{~~deviation (\AA)}\\
		\multicolumn{2}{l}{ ion Wyck.}  & $x$ & $y$ & $z$ & d$x$  & d$z$ \\
		\hline
		$R$1 & $4i$  & $z_{\mathrm{t,}R}$  &     0       & $z_{\mathrm{t,}R}/3$         &  ~~~      & ~~~        \\
		$R$2 & $4i$  & $z_{\mathrm{t,}R}$  &     0       & $(z_{\mathrm{t,}R}+2)/3$     &  ~~~      & ~~~        \\
		$R$3 & $4i$  & $z_{\mathrm{t,}R}$  &     0       & $(z_{\mathrm{t,}R}+1)/3$     &  ~~~      & ~~~        \\
		Cd1 & $4i$  &   $\sfrac14$     &     0       & $\sfrac{11}{12}$                 &  ~~~      & $-0.27$    \\
		Cd2 & $4i$  &   $\sfrac14$     &     0       & $\sfrac{7}{12}$                  &  ~~~      & $+0.23$     \\
		As1 & $4i$  & $z_\mathrm{t,As}$  &     0       & $z_\mathrm{t,As}/3$            &  $-0.05$  & ~~~        \\
		As2 & $4i$  & $z_\mathrm{t,As}$  &     0       & $(z_\mathrm{t,As}+2)/3$        &  $+0.11$  & ~~~        \\
		As3 & $4i$  & $z_\mathrm{t,As}$  &     0       & $(z_\mathrm{t,As}+1)/3$        &  $-0.05$  & ~~~        \\
		As4 & $2b$  &        0         & $\sfrac12$  &     0                            &  ~~~     & ~~~        \\
		As5 & $4i$  &        0         &     0       & $\sfrac16$                       &  ~~~      & $+0.13$    \\
		As6 & $4i$  &    $\sfrac12$    &     0       & $\sfrac13$                       &  ~~~      & $+0.19$    \\
		As7 & $2c$  &        0         &    0        & $\sfrac12$                       &  ~~~      & ~~~        \\
		\hline
	\end{tabular}
		\caption{ \label{Tab02} Description of $R_3$Cd$_2$As$_6$ in the monoclinic setting. The parameters $a_\mathrm{t}$, $c_\mathrm{t}$, $z_{\mathrm{t, }R}\sim0.11$ and $z_\mathrm{t, As}\sim0.34$ refer to the $I4/mmm$ parent cells, as stated in Table~\ref{Tab01}. All sites are fully occupied. The refined $x$ and $z$ coordinates of the Wyckoff site $4i$ for either compound are given in the Supporting Information \cite{sm}. The right column states the resulting approximate deviation (in \AA) with respect to the nominal positions in the tetragonal parent cell.}
	\end{center}
\end{table}

 \begin{figure*}[!t]
	\includegraphics[width=\textwidth]{RxT_05_fixed_scale.pdf}
	\caption{(a) Electrical resistivity as a function of temperature for La$_{3}$Cd$_{2}$As$_{6}$. The inset shows the electrical resistivity of La$_{3}$Cd$_{2}$As$_{6}$ in the  low temperature range. (b) Electrical resistivity as a function of temperature for Ce$_{3}$Cd$_{2}$As$_{6}$.   Natural logarithm of the electrical resistivity as a function of 1000 over the temperature for La$_{3}$Cd$_{2}$As$_{6}$ (c) and for Ce$_{3}$Cd$_{2}$As$_{6}$ (d). The solid red lines are linear fits.}
	\label{RxT}
\end{figure*}

Similar vacancy order has previously been observed in PrZn$_{0.67}$As$_2$ ($P4/nmm$). In that case, the vacancy pattern also results in a stoichiometric ``326'' compound, but with orthorhombic (\textit{Pmmn}) symmetry~\cite{Pr326}. By contrast, the vacancy stripes in Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$ are staggered from one layer to the next in a sequence of six Cd layers. The transformation between the $I4/mmm$ and $C2/m$ cells, as illustrated in Fig.~\ref{Fig2}, is described in Table~\ref{Tab02}. The resulting staircase of Cd-vacancy stripes breaks the parent compounds' fourfold symmetry, which is accompanied by several structural modifications. While the $I4/mmm$ parent structure is determined by two ionic position parameters, twenty are required to define the $C2/m$ counterpart (cf. Table~\ref{Tab02}).

The satisfactory refinement within $I4/mmm$ when superstructure reflections are ignored shows that these distortions are weak. Nonetheless, small deviations of these parameters away from their value in the $I4/mmm$ setting cause significant variations of the structure factors of the superstructure peaks. This is illustrated in the reciprocal x-ray scattering intensity maps in Fig.~\ref{Fig2}. For comparison, each panel shows experimental intensities (left halves) along with our models (right halves). Details of this analysis are provided in the Supporting Information \cite{sm}. Major ionic displacements are summarized in Table~\ref{Tab02} in terms of the deviations (d$x$ and d$z$) relative to the ionic positions in the tetragonal parent compound. The full summary of refined parameters and uncertainties is given in the Supporting Information \cite{sm}.

The refinement adequately captures that the Bragg intensities at $(\frac230l)$ are significantly weaker than those along $(\frac430l)$. This is due to a ``relaxation'' of Cd ions by about $0.25$\,\AA~towards the vacancy stripes, as illustrated in Fig.~\ref{Fig1}(c). As expected, the displacement of Cd ions also leads to a measurable displacement of their ligands (As1, As2 and As3), which causes a characteristic modulation of superstructure intensities along the $l$ direction in reciprocal space.

Importantly, the arsenic square nets in the new ``326'' structure are also distorted. Here, the high-symmetry As1 site of $I4/mmm$ splits into As4 -- As7 (in $C2/m$). The As-As bond-lengths are then modulated by about $0.2$~\AA, as illustrated in Fig.~\ref{Fig1}(b), and become shorter than hypervalent Sb-Sb distances. In the distorted case, the tolerance factor defined by $t = d_{\textmd{As-As}}/d_{\textmd{Ce/La-As}}$ is no longer valid \cite{Sqrnets_1,Sqrnets_2}. Whether the square-net distortion leads to Cd vacancies, similar to K$_{1/3}$Ba$_{2/3}$AgTe$_{2}$ in Ref. [22], or vice-versa remains an open question.

Finally, we note that no major structural differences have been observed between La$_{3}$Cd$_{2}$As$_{6}$
and Ce$_{3}$Cd$_{2}$As$_{6}$, aside from the expected lattice contraction from La to Ce, namely the volume
of the Ce unit cell being smaller by 2.1\,\%. Nonetheless, the La$_3$Cd$_2$As$_6$ crystal investigated
featured only two of the four expected crystallographic twins corresponding to four propagation directions
of the vacancy staircase, as observed in Ce$_{3}$Cd$_{2}$As$_{6}$ [Figs.~\ref{Fig2}(a,b)]. This observation could stem from extrinsic reasons, such as finite strain during crystal growth.

We turn to the physical properties of La$_{3}$Cd$_{2}$As$_{6}$ single crystals. Figures~\ref{RxT}(a) and (b) display the in-plane electrical resistivity ($\rho$) of La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$, respectively. Upon cooling, a sudden increase in the resistivity is observed at $T_{\mathrm{CDW}} = $ 279 K for La$3$Cd$2$As$6$ and at $T_{\mathrm{CDW}} = $ 136 K for Ce$_{3}$Cd$_{2}$As$_{6}$. Below $T_{\mathrm{CDW}}$, the observed semiconducting behavior indicates the opening of an energy gap due to the CDW phase transition, similar to Sr$_{3}$Ir$_{4}$Sn$_{13}$, Ca$_{3}$Ir$_{4}$Sn$_{13}$ \cite{CDWSC}, Ce$_{3}$Co$_{4}$Sn$_{13}$ and La$_{3}$Co$_{4}$Sn$_{13}$ \cite{3413synth,Ce3413-CDW,Ce3413-pressure,La3413-pressure}. Notably,
the resistivity of La$_{3}$Cd$_{2}$As$_{6}$ increases up to
10$^{10}$~$\Omega$\ cm at 12~K, as presented in the inset of Fig.~\ref{RxT}(a).
This large increase of
thirteen orders in magnitude is a rare property in narrow-gap semiconductors, which are typically prone to  impurity bands that may dominate the electronic properties at low temperatures \cite{Book,Book2}. Remarkably, this increase was achieved in an as-grown single crystal, which suggests that sample
quality could be improved through further purification. The broad hump around 50~K for Ce$_{3}$Cd$_{2}$As$_{6}$ may be
associated with the depopulation of the first excited crystal-field state. The magnetic properties of Ce$_{3}$Cd$_{2}$As$_{6}$ will be the focus of a separate study.

Figures~\ref{RxT}(c) and (d) present the resistivity in $\Omega$~m on a logarithmic scale, as a function of $(1000/T)$ for La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$. An activated gap appears to open for temperatures below the CDW transition temperatures in both compounds.  The activation energies can be estimated from linear fits, considering the Arrhenius equation
ln$(\rho_{xx}) =\ $ln$(\rho_{0}) + \Delta/T$, yielding 105(1)~meV and 74(1)~meV for La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$, respectively. We note that the gap values are an estimation and should be taken with caution. Further, it is clear from Fig.~\ref{RxT}(c) that the gap value changes as a function of temperature. A temperature dependent activation energy may be caused by changes in the lattice parameters, due to thermal contraction, as observed in  Ce$_{3}$Bi$_{4}$Pt$_{3}$ \cite{Thermal}. Moreover, the presence of even small amounts of defects create extrinsic channels of conduction, such as variable range hopping and in-gap states \cite{Book,Book2}, which can also be a sensible scenario for a temperature dependent gap. Optical measurements are needed to directly probe the semiconducting gap of La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$.

\begin{figure}[!t]
	\includegraphics[width=0.48\textwidth]{cp_02.pdf}
	\caption{(a) Specific heat over the temperature ($C_{p}/T$) as a function of temperature for Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$. The insets show a zoom in view of the CDW phase transitions. (b) $C_{p}/T$ as a function of $T^{2}$ for La$_{3}$Cd$_{2}$As$_{6}$ The solid red line is a linear fit.}
	\label{cp}
\end{figure}

Figure~\ref{cp}(a) displays the specific heat per mole of lanthanide divided by temperature ($C_{p}/T$) for both Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$
as a function of temperature.  The insets of Fig.~\ref{cp}(a) present a detailed view of the transition temperatures,
which are characterized by peaks at 279~K and 136~K for La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$, respectively, in
agreement with the electrical resistivity measurements. At this temperature scale, the relevance of magnetic Ce correlations to this transition is unlikely. Therefore, given the structural degree of freedom of As ions in the monoclinic structure, we attribute these transitions to a CDW phase caused by a further distortion in the arsenic square net. Figure~\ref{cp}(b) displays $C_{p}/T$ per mole of lanthanide as a function of $T^{2}$ for La$_{3}$Cd$_{2}$As$_{6}$.
The solid red line is the linear fit used to estimate the Sommerfeld coefficient
($\gamma$) of 0.03(3)~mJ/molK$^{2}$ and a Debye temperature ($\theta_{\mathrm{D}}$) of 291(2)~K for La$_{3}$Cd$_{2}$As$_{6}$. A negligible value of $\gamma$ is expected in clean semiconducting compounds
and confirms the high quality of La$_{3}$Cd$_{2}$As$_{6}$ single crystals.

The strong suppression of the CDW transition temperature and activation energy going from La$_{3}$Cd$_{2}$As$_{6}$ to Ce$_{3}$Cd$_{2}$As$_{6}$ is an indication that the application of external pressure may be an effective tuning parameter, given the volume unit cell contraction by about 2~\% in Ce$_{3}$Cd$_{2}$As$_{6}$ as compared to La$_{3}$Cd$_{2}$As$_{6}$. High pressure experiments are needed to investigate whether the CDW phases in La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$ can be suppressed, which may give rise to a superconducting state as in Sr$_{3}$Ir$_{4}$Sn$_{13}$, Ca$_{3}$Ir$_{4}$Sn$_{13}$ \cite{CDWSC} and La$_{3}$Co$_{4}$Sn$_{13}$ \cite{La3413-pressure}. Moreover, chemical pressure may also be an effective tuning parameter, and hence single crystals with smaller lanthanide elements, such as Pr or Nd, are desirable. Finally, spectroscopic experiments are necessary to search for electronic instabilities, which could be driving the CDW phase in La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$.

\section{CONCLUSIONS}

In summary, we report the structural and electrical transport properties of narrow-gap
semiconductors Ce$_{3}$Cd$_{2}$As$_{6}$ and La$_{3}$Cd$_{2}$As$_{6}$, which crystallize
in a distorted variant of the $I4/mmm$ tetragonal structure $R$Cd$_{2/3}$As$_{2}$. The resulting base-centered monoclinic superstructure with $C2/m$ symmetry displays ordered Cd vacancies in a stripe pattern and a distorted As square net.
Notably, both compounds feature charge density wave
phases below 136~K (Ce$_{3}$Cd$_{2}$As$_{6}$) and 279~K (La$_{3}$Cd$_{2}$As$_{6}$), which suggests further distortion in the As square net.
The CDW phases in these materials
create gaps with estimated values of 105(1)~meV and 74(1)~meV for La$_{3}$Cd$_{2}$As$_{6}$ and
Ce$_{3}$Cd$_{2}$As$_{6}$. A remarkable increase of thirteen orders of
magnitude in the electrical resistivity was found for La$_{3}$Cd$_{2}$As$_{6}$ upon cooling from
room temperature to 12~K. Coupled to the vanishingly small specific heat Sommerfeld coefficient, this result not only points to a robust semiconducting ground state in this class of compounds but also provides a route to realize clean narrow-gap semiconductors in distorted arsenic square net materials.

\section{Supporting Information}

The Supporting Information consists of tables summarizing the crystallographic structure data for La$_{3}$Cd$_{2}$As$_{6}$ and Ce$_{3}$Cd$_{2}$As$_{6}$ and the complete crystallographic information files for both compounds at room temperature.

\begin{acknowledgments}
We acknowledge useful discussions with Peter Abbamonte. Work at Los Alamos was performed under the auspices of the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Science and Engineering under project ``Quantum Fluctuations in Narrow-Band Systems". This work was also supported by the S\~ao Paulo Research Foundation (FAPESP) grants 2015/15665-3, 2017/10581-1, 2017/25269-3, CNPq grant $\#$ 304496/2017-0 and CAPES, Brazil. Scanning electron microscope and energy dispersive X-ray measurements were performed at the Center for Integrated Nanotechnologies, an Office of Science User Facility operated for the U.S. Department of Energy (DOE) Office of Science. This research has been supported by the Deutsche Forschungsgemeinschaft through the SFB 1143 and the Würzburg-Dresden Cluster of Excellence EXC 2147 (ct.qmat). Work at Princeton was supported by the Arnold and Mabel Beckman Foundation through a Beckman Young Investigator grant awarded to L.M.S.
\end{acknowledgments}

Materials Science (cond-mat.mtrl-sci); Strongly Correlated Electrons (cond-mat.str-el)


Doping Dependence of the Pseudogap in La(2-x)Sr(x)CuO(4)

Abstract: We report the results of Raman scattering experiments on single crystals of La(2-x)Sr(x)CuO(4) that span the range from underdoped (x = 0.10) to overdoped (x =0.22). The spectra are consistent with the existence of a strong anisotropic quasiparticle interaction that results in a normal state depletion of spectral weight from regions of the Fermi surface located near the zone axes. The strength of the interaction decreases rapidly with increasing hole concentration and the spectral evidence for the pseudogap vanishes when the optimum doping level is reached. The results suggest that the pseudogap and superconducting gap arise from different mechanisms.

\section{Introduction}
It has been recognized for some time \cite{allen,bozovic}
that the normal state electronic properties of the high temperature
superconductors (HTS)
are very different from those of a conventional Fermi liquid. It is also
generally agreed that an understanding of the normal state properties is
important in identifying the physical mechanism
responsible for superconductivity in HTS.
As a result an increasing amount of attention has focused on
the temperature region T$\,>\,$T$_c$  and recent studies
\cite{randeria,puchkov,tatiana,slichter,ino,marshall,harris,ding,norman2,loram1,chen1,hackl,blumberg,quilty,tallon1}
of underdoped compounds have revealed remarkable deviations from Fermi
liquid behavior. These results
\cite{puchkov,tatiana,slichter,ino,marshall,harris,ding,norman2,loram1,chen1,hackl,blumberg,quilty,tallon1}
have been interpreted in terms of the opening of a normal state
pseudogap (PG), a term that is generally used to mean a large suppression
of low energy spectral weight.

Despite an extensive experimental effort the physical
origin of the PG remains undetermined. Some of the proposed mechanisms
involve precursor pairing \cite{randeria,emery,anderson} at T$\,>\,$T$_c$ .
For example Emery and Kivelson \cite{emery} have suggested that
incoherent pairs form above T$_c$, with phase coherence and hence
superconductivity occurring at T$_c$. Or if one has separation of charge
and spin \cite{anderson} d-wave pairing of
spinons  can take place at T$\,>\,$T$_c$, with charge condensation at T$_c$.
On the other hand in the nearly antiferromagnetic Fermi liquid (NAFL) model
proposed by Pines and coworkers \cite{pines1,pines2,pines3},
the PG and superconducting mechanisms are competing. In this model the PG or
spectral weight depletion (SWD) is caused by strong anisotropic
antiferromagnetic (AFM) interactions that are peaked near the AFM ordering
vector {\bf Q}$\,=\,(\pi,\pi)$.

To obtain additional, complementary information on the nature and extent of
the PG regime we have carried out electronic Raman scattering
experiments on La$_{\rm {2-x}}$Sr$_{{\rm x}}$CuO$_{4}$ [La214(x)].  This
compound is an excellent material for these studies for several reasons.
It has a relatively simple structure with a single CuO$_2$ plane in the
primitive cell. Thus the results are not influenced by
structural complications such as those introduced by the  chains
in YBa$_2$Cu$_3$O$_{\rm y}$ (Y123) or the structural modulations
in Bi$_2$Sr$_2$CaCu$_2$O$_{\rm z}$ (Bi2212).
Finally the hole concentration is determined simply by the Sr
concentration if oxygen stoichiometry is maintained.
As a result, one can obtain \cite{kishio} high quality and well characterized
single crystals of La214 that enable one to study the systematic
evolution of the electronic properties throughout the
complete doping range.  We have measured the
relative strengths and frequency distributions of the B$_{1g}$ and
B$_{2g}$ Raman continua as a function of temperature and doping.
The results of these experiments are consistent with a PG mechanism
that is intrinsic to a single CuO$_2$  layer and is present only in
the underdoped regime of La214. Our results also imply that the energy
scale E$_g$ associated with the PG is quite different from the
superconducting gap energy $\Delta$.

\section{Experiment}

The samples studied in this work were grown by a travelling
floating-zone method \cite{kishio}.  The crystals used here were
cut from larger single crystals and typically had dimensions of about
$2 \times 2 \times 0.5$ mm$^3$.  Characterization of the samples included
susceptibility and transport measurements of their critical temperatures
\cite{kishio}.  The sample surfaces were prepared for light scattering
experiments by polishing with diamond paste and etching \cite{werder} with a
bromine-ethanol solution.  Finally the samples were oriented in the basal
plane using Laue x-ray diffraction patterns.  The physical parameters
characterizing the La214 samples are summarized in Table I.

\begin{table}[htb]
\centering
\begin{tabular}{|c|c|lc|}
La$_{\rm {2-x}}$Sr$_{{\rm x}}$CuO$_{4}$&Nominal Sr content (x)&T$_c$
($^{\circ}$K)&\\ \hline
Underdoped&0.10&27&\\ \hline
Underdoped&0.13&35&\\ \hline
Optimally doped&0.17&37&\\ \hline
Overdoped&0.19&32&\\ \hline
Overdoped&0.22&30&\\
\end{tabular}
\vspace{0.10in}
\caption{The physical parameters characterizing the samples studied in
this paper. T$_c$ was determined magnetically and
x is the nominal hole concentration$^{23}$.}
\label{tab:xT$_c$}
\end{table}

The Raman spectra were obtained in a quasibackscattering geometry using
the 488.0 or 514.5 nm lines of an Ar$^+$ laser.  In most cases the
light was focussed onto the sample using a cylindrical lens, but on occasion
a combination of a spherical and cylindrical lenses was used to reduce the
dimensions of the laser beam at the sample.  In all cases the incident
power level was controlled to be less
than 10 W/cm$^2$.  With these power levels an estimate of about 11$\,$K has
been obtained for the amount of local heating by the laser.
This estimation was obtained from measurements of the power dependence of
the spectra at temperatures near T$_c$, and the Stokes/anti-Stokes ratio
at somewhat higher temperatures.  The temperatures reported
in this paper are the ambient value plus 11$\,$K,
or the estimated temperature in the excited region of the sample.

In a Raman spectrum the intensity of the scattered light is proportional
to the square of the Raman tensor $\gamma$, whose components can be
selected by appropriate choices of the incident and scattered polarizations
\cite{chen2}. In this paper spectra were obtained in the
$z(xy)\bar{z}$ and $z(x^{\prime}y^{\prime})\bar{z}$
scattering geometries where the letters inside the brackets denote
the polarizations of the incident and scattered light.
Here the $x$ and $y$ axes are parallel to the Cu-O bonds
and $x^{\prime}$ and $y^{\prime}$ axes are rotated by 45$^{\circ}$
with respect to $(x,y)$.
The xy component of the Raman tensor $\gamma_{xy}$ transforms as k$_x$k$_y$
or the B$_{2g}$ irreducible representation of the D$_{4h}$ point group,
and $\gamma_{x^{\prime}y^{\prime}}$ as (k$^2_x-\,$k$^2_y$) or
B$_{1g}$.  Thus $\gamma$(B$_{1g}$) must vanish along the diagonal directions
($\pm 1,\pm 1$) in k-space and $\gamma$(B$_{2g}$) must vanish along the
k$_x$ and k$_y$ axes.  These two scattering geometries allow us to probe
complementary regions of the FS \cite{chen2} since the B$_{1g}$ channel will
sample regions of the FS located near the
k$_x$ and k$_y$ axes while the B$_{2g}$ channel will probe regions located
near the diagonal directions.

Throughout this paper we deal with the frequency distribution
of the Raman response function $\chi^{\prime\prime}(\omega,$T)
which is obtained by dividing the measured intensity by the thermal
factor [1--exp(-$\hbar \omega$/$k_B$T)]$^{-1}$.

\section{Results}

\begin{figure}[htb]
\centerline{\epsfxsize=2.5 in \epsffile{fig1.eps}}
\vspace{0.10in}
\caption{ The B$_{1g}$ (a) and B$_{2g}$ (b) response functions
of the underdoped La214(0.13) crystal for temperatures between 50K and 300K.
The spectra are offset vertically for clarity.}
\label{1}
\end{figure}

A detailed description of Raman spectra of La214(0.17) sample,
which we will consider to be representative of the optimally doped state,
has been published recently \cite{chen2}. In this paper
we will focus on the changes that occur in the spectra of the underdoped
and overdoped compounds. We will first discuss spectra
obtained from the underdoped samples (x = 0.10 and 0.13). The normal state
response functions are shown in Fig. 1 for the 0.13 sample and the low
temperature response functions for both samples are presented in Fig. 2.
As is shown in Fig. 1a, the magnitude of $\chi^{\prime\prime}$(B$_{1g}$)
is independent  of temperature, to within the accuracy of our measurements
(15\%). From figures 2a and 2b, it is evident that the B$_{1g}$ spectra
obtained for each sample, above and below T$_c$, are essentially identical
and are thus unaffected by the superconducting (SC) transition.
In direct contrast, the low energy B$_{2g}$ spectra of both
underdoped compounds (Figs. 2c and 2d) undergo a superconductivity
induced renormalization (SCIR) as the temperatures of the samples are lowered
through  T$_c$.
The observed immunity of the B$_{1g}$ spectra to passage into the
SC state, and contrasting SCIR in the B$_{2g}$ channel,
is analogous to the behavior observed  in underdoped Y123 \cite{chen1}.

\begin{figure}[htb]
\centerline{\epsfxsize=2.5in \epsffile{fig2.eps}}
\vspace{0.10in}
\caption{The B$_{1g}$ response above and below T$_c$  for the La214(0.13)
crystal (a) and the La214(0.10) crystal (b).  The B$_{2g}$ response for the
same crystals are shown in (c) and (d).}
\label{2}
\end{figure}
\begin{figure}[htb]
\centerline{\epsfxsize=2.7 in \epsffile{fig3.eps}}
\vspace{0.12in}
\caption{Direct Comparison of the low energy B$_{1g}$ and
B$_{2g}$ Raman response functions measured at temperatures
below T$_c$ in La214(x) for (a) x = 0.10; (b) x = 0.13; (c) x = 0.17;
(d) x = 0.19; and (e) x = 0.22.  The scale of $\chi^{\prime\prime}$ is
the same for all five frames.}
\label{3}
\end{figure}

We have also obtained spectra from two overdoped crystals (x = 0.19
and 0.22) and we have found that the magnitude of the response functions
$\chi^{\prime\prime}$(B$_{1g}$) in both samples
(for $\omega \geq \,300\,cm^{-1}$)
are again approximately independent \cite{jafar1,jafar2} of temperature.
In contrast to the situation in the underdoped samples, there is a pronounced
SCIR that occurs at T$_c$ in {\em both} the B$_{1g}$ and B$_{2g}$ channels
for both the optimally doped \cite{chen2} and the overdoped samples
\cite{jafar1,jafar2}.  Finally the low temperature spectra of the overdoped
crystals, in both  B$_{1g}$ and B$_{2g}$ channels, are compared to the
underdoped and optimally doped spectra in Fig. 3.

\begin{figure}[htb]
\centerline{\epsfxsize=2.4 in \epsffile{fig4.eps}}
\vspace{0.12in}
\caption{The ratio R of the integrated Raman response functions plotted as
a function of hole concentration for La214 ($\bullet$) and Y123 ($\circ$).
The ratio was determined for $0 \leq \omega \leq 600\,cm^{-1}$
in all cases. The dashed line represents the calculated R using equation (1)
for different values of $\phi$ (Fig. 5) superimposed using the top scale.
The dotted line serves only as a guide to the eye.}
\label{4}
\end{figure}

As a measure of the strength of the response functions we will use the
integral of $\chi^{\prime\prime}(\omega)$ over the low frequency
spectral region.
In this context the strength of $\chi^{\prime\prime}$(B$_{2g}$) is
approximately the same (Fig. 3) for all the crystals studied here,
while the strength of $\chi^{\prime\prime}$(B$_{1g}$) changes dramatically
as we proceed  from the underdoped to the overdoped samples.
The ratio R = $\int{\rm d}\omega\, \chi^{\prime\prime}$(B$_{1g}$) /
$\int{\rm d}\omega \,\chi^{\prime\prime}$(B$_{2g}$),
and therefore the strength of the B$_{1g}$ Raman response, increases
by a factor of about 16 as we go from the underdoped x = 0.10 sample to
the optimally doped sample with x = 0.17.
The ratio R continues to increase, but at a slower rate as one proceeds into
the overdoped regime. This result is summarized quantitatively
in Fig. 4, where we have plotted the ratio R as a function of hole
concentration. For comparison purposes, data obtained previously \cite{chen1}
on Y123  are also included \cite{tallon2} in Fig 4. It is evident that the
ratio R for the two compounds exhibits a similar trend with doping.

\section{Discussion}

\subsection{Spectral Weight Depletion}

Relatively few Raman experiments have been carried out to
investigate the nature of the PG. Nemetschek  {\em et al}. \cite{hackl}
have attributed a small loss of spectral weight at $\omega \leq 700\,cm^{-1}$
in the B$_{2g}$ channel to the presence of a PG.  Blumberg {\em et al}.
\cite{blumberg}, in spectra obtained from underdoped Bi2212 have observed
a sharp feature at about 600$\,cm^{-1}$ which they have attributed to a bound
state associated with precursor pair formation above T$_c$.
In our studies of underdoped Y123 and La214 we have focused on
the striking loss of spectral weight in the B$_{1g}$ channel,
and the contrasting constancy of the B$_{2g}$ spectra, features which appear
to be common \cite{chen1,katsufuji,opel} to the hole doped cuprates.

\begin{figure}[htb]
\centerline{\epsfxsize=2.6 in \epsffile{fig5.eps}}
\vspace{0.14in}
\caption{Calculated B$_{1g}$ (thin lines) and B$_{2g}$ (thick lines)
Raman response functions in La214 (a) by integrating over the full Fermi
surface (dashed lines in the inset) and (b-d) by integrating only over region
lying between $\phi$ and $90^{\circ} - \phi$ in each quadrant. The inset
shows polar plots of the angular dependence of the
B$_{1g}$ (thin lines) and B$_{2g}$ (thick lines) components of the
Raman tensor $\gamma({\bf k})$.}
\label{5}
\end{figure}

To illustrate the dramatic change in the intensity in the B$_{1g}$ channel
and the relative constancy of the B$_{2g}$ channel, we can calculate the
unscreened Raman response functions in the superconducting state
using the conventional model of light scattering \cite{klein}
\begin{equation}
\chi^{\prime\prime}_{\gamma}(\omega)
\propto \, \langle\, \frac{\gamma^{2}({\bf k}) \Delta^{2}({\bf k})}
{\omega \sqrt{\omega^2-4\Delta^{2}({\bf k})}}\,\rangle,
\label{t0}
\end{equation}
where $\langle \cdots\rangle$ imply an average over the FS for k-vectors such
that $\omega > 2|\Delta({\bf k})|$. The results shown in Fig. 5 were obtained
with $\gamma({\bf k})$ calculated using a second neighbor tight binding
model with the parameters described previously \cite{chen2} and
$\Delta({\bf k})=\Delta_{0}[\cos(k_{x}a)-\cos(k_{y}a)]$. The response
functions obtained by integrating (1) over the full FS (dashed lines in the
inset of Fig. 5) are shown in Fig. 5a. To obtain a representation of
the Raman response in the underdoped case we integrate (1) over portions of
the FS lying between $\phi$ and $90^{\circ} - \phi$ in each quadrant.
That is, regions of the FS near the axis (defined by angle $\phi$) were
excluded from the integration, and the resulting response functions are
shown in Figs. 5b-5d for $\phi = \,$10$^{\circ}$, 20$^{\circ}$, and
30$^{\circ}$. As is evident an increasing truncation of the FS leads to
a rapid depletion of the B$_{1g}$ spectrum, but has much smaller
effect on the B$_{2g}$ spectrum.

We have also calculated the ratio R of the integrated response functions
as a function of $\phi$ and the results are superimposed (dashed line)
on the data points in Fig. 4. For the purpose of this comparison
both the measured and calculated response functions were integrated
from $\omega = 0$ to 600$\,cm^{-1}$.
It should be noted that the apparent agreement between calculation and
experiment is quite fortuitous given the simplistic nature of the model.
The results do, however, provide plausibility to the existence of
localized depletion of spectral weight from regions of the FS located
near ($\pi,0$).
Finally we should note that the effective gapping of quasiparticles (QP)
from these same region, at T $>$ T$_c$, is  consistent with the absence of
a superconductivity induced renormalization in the B$_{1g}$ spectra of the
underdoped compounds (Figs. 2a and 2b).

The observed dependence on scattering geometry thus implies that the spectral
weight loss is confined to regions of the FS located near ($\pm \pi,0$) and
($0,\pm \pi$) and that QP located near the diagonals are essentially
unaffected by a reduction in doping level.
Given that ARPES experiments \cite{marshall,harris,ding} also find a SWD that
reaches a maximum in regions near ($\pi,0$), we will assume that this feature
is characteristic of the PG, and that the PG thus manifests itself in the
Raman spectra {\em primarily} as a SWD in the B$_{1g}$ channel.
We should note however, that the B$_{2g}$ spectra is also
approximately independent of temperature (for T$\,>\,{\rm T}_c$), which
is difficult to reconcile with the FS contraction envisaged by Norman
{\em et al.} \cite{norman2}.

According to Fig. 4 the SWD or pseudogap has approximately
the same strength in La214 and Y123. This is not consistent with the
suggestion \cite{mm} that the PG arises from strong AFM interactions between
the closely spaced CuO$_2$ layers of bilayer compounds such as Y123.
The present results thus suggest that the PG arises from intralayer
interactions and in this regard are in agreement with IR
\cite{puchkov,tatiana}, ARPES
\cite{ino}, specific heat \cite{loram3}, and NMR \cite{yasuoka} measurements.

\subsection{Magnitude (E$_g$) of the Pseudogap}

In the spectra we have obtained from La214 there is no clear indication of
an onset temperature T$^*$ for the PG. Estimates \cite{pines2,yasuoka}
of T$^* \leq 200\,$K have been obtained from the analysis of NMR experiments
carried out on slightly underdoped La214 samples (x = 0.14 and 0.15).
On the other hand transport \cite{batlogg}, IR reflectivity \cite{tatiana},
and specific heat \cite{loram3} measurements on underdoped La214
have all found T$^*  > 300\,$K.
Our spectra (Fig. 1) indicate some minor changes at about
200$\,$K, that include the appearance of a broad and weak feature at about
600$\,cm^{-1}$ in both channels. In fact the variations that occur in
the B$_{2g}$ spectra below about 200$\,$K reflect a decrease in intensity
at low energies similar to that observed by Nemetschek  {\em et al}.
\cite{hackl} in Bi2212. However, no significant changes in spectral weight,
which one might expect to mark the onset of the PG, occur at any
temperature below $300\,$K.

The lack of an experimental consensus for T$^*$ in La214 is perhaps
supportive of the suggestion \cite{tallon3} that one cannot associate a well
defined T$^*$ with the PG, but only an energy scale E$_g$. In this regard,
however, there is also no clear evidence of a gap edge in the B$_{1g}$
spectra we have obtained (Fig. 1) from underdoped La214. This is perhaps
not surprising if E$_g$ is of order of J in underdoped samples and decreases
with increasing hole concentration  as has been suggested by photoemission
\cite{ino} and specific heat \cite{loram3} measurements.
This magnitude for the PG might be expected if short range AFM
correlation were responsible for the spectral weight depletion \cite{pines1}.
One possible scenario involves the existence of
AFM ordered regions and separate hole rich regions in the underdoped
samples. Evidence for this type of phase separation was obtained recently
in the muon spin resonance experiments on underdoped La214 and Ca doped Y123
\cite{ndmyr}.

\begin{figure}[htb]
\centerline{\epsfxsize=2.8in \epsffile{fig6.eps}}
\vspace{0.12in}
\caption{The high energy B$_{1g}$ Raman response functions measured at
50$\,$K and 100$\,$K in La214(x) for (a) x = 0.22, (b) x = 0.13, and (c)
x = 0.10
}
\label{6}
\end{figure}

In Raman spectra one-magnon excitations are not observed and hence the
PG energy would be reflected in the two-magnon peak that occurs at about
3J $\approx$ 3000$\,$cm$^{-1}$ in undoped La214 \cite{lyons,sugai} and
insulating Y123 \cite{lyons,reznik,blum}. The amplitude and energy of the
two-magnon peak decreases \cite{sugai,rubhausen} as the compounds are doped
and such peaks have not been observed for doping levels greater than or equal
to optimum. Although our spectrometer range is limited,
we have carried out preliminary measurements of the high energy B$_{1g}$
spectra for three samples (x = 0.10, 0.13, and 0.22).
The high energy ($\omega \,<\, 3000\,cm^{-1}$) response functions measured
at 50$\,$K and 100$\,$K for the three samples are shown in Fig. 6.
From the figure one can see a broad, weak peak at about 2000$\,$cm$^{-1}$
(J$\,\sim 700\,$cm$^{-1}$) in the x = 0.10 spectrum. This peak
essentially vanishes as we proceed to higher doping levels (Figs. 6a and 6b).
We must emphasize that these spectra have not been corrected for
variations in the optical constants or for
spectrometer response and are shown only to illustrate the effect of doping.
Our spectrometer response decreases in the red and correction might
push the peak in Fig. 6c to higher energies and also decrease its relative
amplitude. Such corrections would not however alter the observed trend that is
induced by changes in doping.

It appears plausible that the pseudogap or SWD in the B$_{1g}$
channel is due to short range AFM correlations and the formation of
``hot spots'' \cite{pines1} on regions of the FS located near ($\pm\pi,0$)
and ($0,\pm\pi$).
As the doping level is increased the AFM interactions become weaker,
as is evidenced by a weakening (Fig. 6) and shift of the two-magnon peak
to lower energies \cite{sugai,rubhausen}. Thus as the doping level is
increased, E$_g$ decreases and some spectral weight is shifted from higher
to lower energies.  At optimum doping, where a SCIR is observed in the
B$_{1g}$ channel (Fig. 3c), the PG is assumed to be filled.
Furthermore on the basis of their specific heat measurements in several
cuprates Loram {\em et al.} \cite{loram3} have proposed
that E$_g \approx\,$J$\,(1-\,p/p_{cr})$ where $p$ is the hole concentration
and $p_{cr}$ the value at which E$_g$ ``closes''. Given that
J$\,\approx 1000\,$cm$^{-1}$ in La$_2$CuO$_4$ and E$_g\,\approx\,
700\,$cm$^{-1}$ (Fig. 6c) suggests $p_{cr}\, > \,p_{opt}$. We would
thus conclude that the pseudogap ``fills'' before it ``closes''.

It is interesting to compare the PG energy scale E$_g$
to that of the superconducting gap $\Delta$. In Raman spectra the energy of
the pair-breaking peak in the B$_{1g}$ channel can be used \cite{tom1}
as a measure of $\Delta$. In optimally doped La214 the maximum value of
the superconducting gap \cite{chen2} is given by
$\omega$(B$_{1g}) \approx 2\Delta_{max} \approx 200\,$cm$^{-1}$,
as is also evident from Fig. 3.
In the overdoped region $\Delta_{max}$ decreases (Fig. 3), consistent with
the behavior in other cuprates \cite{chen3,white,kendziora}.
In the underdoped regime the absence of SCIR in the B$_{1g}$ channel
means that we do not have any direct measure of the gap maxima.
We can note however that the frequency of the pair-breaking peak in the
B$_{2g}$ channel remains approximately constant
[$\omega$(B$_{2g}) \approx 130 \pm 20\,$cm$^{-1}$] as the doping level is
decreased from optimum.  We thus infer that $\Delta$ remains approximately
constant throughout the doping range $0.10 \leq {\rm x} \leq 0.17$.
This behavior is also consistent with that found in ARPES and tunneling
measurements \cite{harris,renner} on Bi2212. Therefore we have an amplitude
$\Delta_{max} \approx 100\,$cm$^{-1}$ for the SC gap and E$_g \sim\,$J$\,\sim
700\,$cm$^{-1}$ for the PG. These very different energy scales suggest
that different physical mechanisms are involved.
\subsection{Summary}

The most noticeable aspect of the spectra shown in Fig. 3 is the
SWD that occurs in the B$_{1g}$ channel when the doping level is reduced
below optimum, while the strength of the B$_{2g}$ spectrum is essentially
unaffected by changes in doping.
These results  suggest that the QP scattering rate is highly anisotropic
\cite{tom2} and leads to  the existence of ``hot spots''
on regions of the FS located near ($\pm\pi,0$) and ($0,\pm\pi$).
This observation is consistent with the predictions of Pines and
coworkers \cite{pines1} who have pointed out that QP located near these
regions can be coupled by
the AFM ordering vector {\bf Q} = ($\pi,\pi$) and hence interact strongly
in underdoped compounds. Conversely, the QP on regions of the FS located
near  $|{\bf k}_x|=|{\bf k}_y|$ are weakly coupled
\cite{pines1}, and have been designated as ``cold''  \cite{pines1,ioffe}.
As the doping level is increased toward optimum the SWD decreases rapidly,
consistent with a weakening of the AFM interactions between QP \cite{pines1}.
In this language our results then suggest that the hot QP are gapped above
T$_c$ and hence are not affected by the SC transition. Only the cold QP,
as evidenced by the SCIR in the B$_{2g}$ channel at T$_c$,
participate in the transition to the SC state.

Thus far our attention has been focused on the most dominant feature of the
spectra obtained from underdoped compounds, namely the SWD that is confined
to the B$_{1g}$ channel. As mentioned previously Nemetschek {\em et al.}
\cite{hackl} found that a relatively small depletion of spectral weight
also occurs in the B$_{2g}$ channel of Y123 and Bi2212 below
$\omega \approx 700\,$cm$^{-1}$ and T$\approx 250\,$K. This energy and onset
temperature are in good agreement with IR measurements on the same compounds
\cite{puchkov}. This appears reasonable since, given the SWD in
the B$_{1g}$ channel one would expect the results of IR and transport
measurements, for example, to be determined mainly by the properties of the
cold QP \cite{hackl,pines1,ioffe}, or those on regions of the FS located near
the diagonal directions. Furthermore, one can see that a similar small
depression of spectral weight appears below 700$\,$cm$^{-1}$ and 200$\,$K
in the B$_{2g}$ spectra of the x = 0.13 crystal (Fig. 1b), in agreement with
IR measurements carried out \cite{tatiana} on underdoped La214.
Thus the loss of spectral weight in the B$_{2g}$ channel \cite{hackl},
the depression of the scattering rate observed in IR experiments, and
the spectral weight depletion observed in the B$_{1g}$ channel (Figs. 3 and 4)
are all characterized  by the same energy scale. This suggests that there
is a common physical mechanism underlying the apparently different phenomena
that are observed in the two channels.
However, the identification of a possible interrelation between the
two channels will require a very careful analysis of
temperature dependence of the spectra for 50$\,$K$\leq\,$T$\,\leq 300\,$K.
Such an analysis should also provide important information about the existence
and magnitude of crossover temperatures T$^*$ that have
been observed in other experiments \cite{tatiana,pines2,yasuoka,batlogg}
and are important features of some models \cite{emery,pines1}.
It is also possible that the SWD in the two channels arises from different
mechanisms and that there are different pseudogaps as have been suggested
 \cite{ino} by  Ino {\em et al}.

\section{Conclusions}
In the Raman spectra obtained from La214 there is a loss of low energy
spectral weight in the B$_{1g}$ channel as the doping level is decreased
below optimum, while the intensity in the B$_{2g}$ channel remains unchanged.
These results imply the existence of a strong anisotropic
scattering mechanism and the formation of hot spots \cite{pines1} on regions
of the Fermi surface located near k$_x$ and k$_y$ axes.
Furthermore, in underdoped compounds the B$_{1g}$ channel is unaffected by
the superconducting transition.
The pseudogap is thus manifested in Raman spectra both by a depletion
of low energy spectral weight or reduced intensity and by the absence of a
superconductivity induced renormalization in the B$_{1g}$ channel.
In this context the pseudogap in La214 and Y123 \cite{chen1} is present
only in the underdoped regime. Since the spectral weight depletion
is similar in the two compounds, we conclude that the
pseudogap is intrinsic to a single CuO$_2$ plane and does not arise
from interlayer interactions.
Furthermore the results suggest that the loss of low energy
intensity in the B$_{1g}$ channel is due to strong antiferromagnetic
interactions which in turn implies E$_g \approx\,$J$\,\approx 700\,$cm$^{-1}$
for the pseudogap, an energy scale that is much larger than the energy
of the superconducting gap ($\Delta_{max} \approx 100\,$cm$^{-1}$).
These energy scales imply that the pseudogap and superconducting gap
arise from different physical mechanisms. Since superconductivity grows
while the spectral weight loss or  pseudogap diminishes, it appears that
the two mechanisms compete with each other for the available quasiparticles.

\section{Acknowledgements}
The authors gratefully acknowledge useful conversations with T. P. Devereaux,
T. Startseva, T. Timusk, and J. Schmalian. One of us (JCI) has also benefited greatly
from discussions with R. G. Buckley, J. L. Tallon, A. J. Trodhal, and
G. V. M. Williams during a visit to the Industrial Research Laboratory in
Lower Hutt, New Zealand. The financial support of the Natural Sciences and
Engineering Research Council of Canada is gratefully acknowledged.