High passive-stability diode-laser design for use in atomic-physics experiments

Abstract: We present the design and performance characterization of an external cavity diode-laser system optimized for high stability, low passive spectral linewidth, low cost, and ease of in-house assembly. The main cavity body is machined from a single aluminum block for robustness to temperature changes and mechanical vibrations, and features a stiff and light diffraction-grating arm to suppress low-frequency mechanical resonances. The cavity is vacuum-sealed, and a custom-molded silicone external housing further isolates the system from acoustic noise and temperature fluctuations. Beam shaping, optical isolation, and fiber coupling are integrated, and the design is easily adapted to many commonly used wavelengths. Resonance data, passive-linewidth data, and passive stability characterization of the new design demonstrate that its performance exceeds published specifications for commercial precision diode-laser systems. The design is fully documented and freely available.

\section{Introduction\label{sec:introduction}}

Atomic physics experiments require narrow-linewidth laser sources tuned to specific atomic resonances, and tunable external-cavity diode lasers (ECDLs) are often the light source of choice as a compact, reliable, low-cost, energy-efficient option.\cite{Wieman91}  Typical spectral linewidths on the order of 1~MHz are sufficient for many laser cooling and trapping experiments, but precision spectroscopy of narrow-linewidth transitions requires narrower-linewidth
sources. Linewidths can be narrowed by servo-locking to high-finesse external reference cavities,\cite{Young99} but the requirements are less stringent and the servo is more robust when starting with an inherently narrow and stable source.

Commercial precision diode-laser systems can cost \$20,000--\$30,000; as a result, many designs for Littrow configuration ECDLs that can be built in-house already exist.\cite{MacAdam92, Arnold98, Ricci95, Vassiliev06}  In the standard Littrow configuration, the lasing cavity is formed by the rear facet of the laser diode and a diffraction grating reflecting in first order, whose angle can be changed to tune the emission wavelength.  Temperature fluctuations and acoustic vibrations that couple to the diffraction-grating mounting arm affect the cavity length and thereby change the emitted frequency, so several groups have pursued novel designs to avoid these issues.\cite{Zorabedian88, Baillard06, Papp07thesis}  A promising example is the ``cat-eye'' laser,\cite{Zorabedian88, Baillard06} which abandons the diffraction grating in favor of a tiltable interference filter as the frequency discriminating element, and the external cavity is composed of two lens-mirror ``cat-eye'' reflectors.  One model with a long (10~cm) cavity has demonstrated a passive linewidth of 155~kHz.\cite{Baillard06}  Another design machines a Littrow-configuration external cavity out of a single block of aluminum to reduce vibrational coupling to environmental noise, and the diffraction grating is mounted on an adjustable kinematic mirror mount integrated into the main cavity block.\cite{Papp07thesis}

\begin{figure}[ht]
\centering
\includegraphics[width=9 cm]{fig1.eps}
\caption{Diagram of the (a) main laser body, (b) diode
collimation-tube mount, and (c) optical isolator mount (``aspirin tablet'' size). Not pictured are the lids and baseplate. All pieces are machined from 6061-T6 aluminum.
\label{fig:cavity}}
\end{figure}

\begin{figure*}[ht]
\centering
\includegraphics[width=14 cm]{fig2.eps}
\caption{Laser assembly:
cutaway top view (with all peripherals but no electrical connections) and front view (with lids and baseplate).
\label{fig:diagram}}
\end{figure*}

Our design improves upon
the basic single-block design in several ways.  In particular, the diffraction-grating arm itself is also machined as part of the main cavity block. The arm is made as stiff and light as possible to push its mechanical resonance frequency above most background noise sources.
The main cavity is hermetically/vacuum sealed to allow evacuation or backfilling with an inert gas to prevent condensation and reduce sensitivity to humidity and pressure fluctuations.  Beam shaping, optical isolation, and fiber coupling are all integrated into the assembly.  This design is readily adaptable to many different wavelengths, simply by switching the laser diode and diffraction grating (with an added shim for some wavelengths), and adjusting the beam-shaping prisms as necessary.  Gratings that would work in the geometry of our design are commercially available for many relevant atomic-physics wavelengths.

\section{Construction\label{sec:construction}}

The main cavity (Fig.~\ref{fig:cavity}a) is machined on a computer numerical controlled (CNC) milling
machine from a $6"\times 5"\times 1"$ block of 6061-T6 aluminum.  Use of CNC machinery
allows easy replication of the design, and complicated geometries like the diffraction-grating arm and a
ramped surface to mount a Brewster-angle window are quickly and easily machined into the main laser body.  A pump-out port mimicking a commercial seal-off valve (DLH Industries V1021-1, to be used in conjunction with valve operator V1025-3-25) is also machined into the front face for evacuation, and all potential air traps are vented.
Additional modules to hold the diode collimating tube (Fig.~\ref{fig:cavity}b) and the optical isolator (Fig.~\ref{fig:cavity}c) are machined separately.  The top and bottom lids are machined $0.75"$ thick, and each contains a groove for a Viton O-ring (McMaster 1201T827) to allow vacuum sealing of the
laser cavity.  A baseplate provides mechanical support and a thermal reservoir for the thermoelectric coolers.

The main cavity body is designed to be machined from only three orientations to reduce setup labor and overall machining costs.  The diffraction-grating arm contains four 2-56 counterbored clearance holes, irregularly spaced along the length of the arm. These holes serve both to reduce the mass of the diffraction-grating arm and to provide a means to stabilize the arm during machining of main cavity body.  Additional details of the machining process are in the Appendix.

Figure~\ref{fig:diagram} shows the laser assembly with all peripherals, and Fig.~\ref{fig:unilaser_pic} shows a photograph of the laser interior, illustrating the electrical connections.  Assembly begins with a thorough (vacuum-quality) cleaning of all machined parts. Throughout construction, cleanliness is emphasized, not only to produce a sterile, vacuum-friendly environment, but also to ensure an absence of dust and foreign material. This is especially important if the laser diode has no protective housing (e.g., for antireflection-coated diodes).  Vacuum-compatible epoxy (Loctite Hysol 1C) is used to strengthen delicate solder joints and affix the diffraction grating to the arm.

The diode is held in a collimating lens tube (Thorlabs LT230P-B) which is in turn clamped in the
collimation-tube mount (Fig.~\ref{fig:cavity}b).  This module is firmly bolted to the wall of the main cavity body, with access to the mounting bolts closed and sealed using 3/8-24 screws and Viton O-rings (McMaster 1201T23).  The face that contacts the main cavity body is lightly coated with a vacuum grease (Apiezon M) to aid in heat transfer and to help prevent long-term creep, which could otherwise occur from mechanical stresses created when the module is fixed to the main cavity body.

Prior to firmly securing the diode collimation tube module, the vertical tilt of this module can be adjusted with the help of two temporary screws. These screws have rounded ends, and are threaded vertically through the front and back of the mount, supporting the module against the bottom lid.
By careful differential adjustment of the screws, small changes in vertical tilt can be made. As the alignment is optimized, the lasing threshold is observed until a minimum is found, and this minimum is monitored as the module is secured and the tilt-adjust screws are removed.  The vertical degree of freedom is then fixed, and the collimation-tube mount is effectively a part of the rigid cavity body: tapping and pushing on the module has no effect on cavity alignment.

\begin{figure}[ht]
\centering
\includegraphics[width=8.5cm]{fig3.eps}
\caption{
Assembled laser: close-in view of the main cavity of the laser with the top lid removed.
}
\label{fig:unilaser_pic}
\end{figure}

The diffraction-grating mounting arm is a critical element of the laser, as its mechanical stability is a typical weak point of most diode-laser designs.  It is carved directly into the main cavity block, and designed to be as stiff and light as possible. The grating itself is custom-manufactured to a small size (6~mm high, 12~mm wide, 3~mm thick) in order to decrease the total mass of the grating arm, keeping mechanical resonance frequencies above most of the ambient acoustic noise.  Gross adjustments of the diffraction grating angle can be made using a precision 1/4-100 screw (New Focus 9376-K), with an access port sealed using another 3/8-24 screw and O-ring.  Fine wavelength adjustments are possible using a low-voltage piezo (Noliac SCMAP02/S2/A/5/5/10/60/10.6/1000) epoxied between two sapphire discs (Swiss Jewel Company W7.87) and inserted between the precision screw and the back of the grating mounting arm.  The piezo is driven by an ordinary, 15~V op-amp circuit (buffered by a high-current, unity-gain amplifier, BUF634), which typically tunes the laser $\sim\!2$~GHz/V.   A small O-ring sits in a groove in the bottom lid, compressed against the bottom of the tip of the diffraction-grating arm to provide shear damping of the arm motion.

Electrical connections enter through a D-sub 15 female connector on the front face of the laser.
Polyimide-coated wires (MWS 23 AWG, HML insulation) pass through channels potted with a high-temperature epoxy (Epotek 353ND) before proceeding to two temperature sensors, the piezo, and the diode (through a simple protection circuit\cite{MacAdam92, Meyrath03notes}). We use a 50-k$\Omega$ thermistor for feedback control of temperature and an AD590 sensor
to drive an LCD temperature display.  Two thermoelectric coolers (Laird 56460-501) are positioned between the bottom lid and the baseplate, with a thin layer of thermal conductive paste (Arctic Silver 5) to aid heat transfer.

\begin{figure}[ht]
\centering
\includegraphics[width=8.5cm]{fig4.eps}
\caption{Assembled laser mounted on table: A composite block and viscoelastic cushioning material isolate the laser from table vibrations, and a flexible, molded silicone cover offers further acoustic isolation and temperature regulation.}
\label{fig:mounting}
\end{figure}

The main cavity is sealed at the beam output by a window cut from a microscope slide. The window is glued with the thermal epoxy to a ramp cut at Brewster's angle in the main cavity body. The exiting beam then passes through an anamorphic prism pair (Thorlabs PS871-B), with the prism positions calculated to optimize the symmetry of the exiting beam and to center it on the correct output location.  Beam paths and prism locations are laser-etched onto a 0.064"-thick acrylic alignment template and fixed to the bottom of the prism cavity with five-minute epoxy (the guide can also be printed on sturdy paper).  The prisms are then attached in place on the template (after confirming proper alignment of the output beam) with small drops of five-minute epoxy.
The beam then passes through an optical isolator mounted in a separate aluminum block bolted to the main cavity body (Fig.~\ref{fig:cavity}c).  We have two versions of the isolator mount---one for small ``aspirin-tablet'' isolators (OFR  body type D) and another for larger versions (OFR body type II). Finally,  the light is output to a fiber coupler (Oz Optics HPUC-23A series) which is attached directly to the isolator mount. This minimizes the need for future adjustments of the fiber coupling and ensures that failed laser diodes can be replaced without the hassle of realigning an entire beam line.

Vibrational isolation of the laser is enhanced via additional damping components (Fig.~\ref{fig:mounting}).  For decoupling from the optical table, the assembled laser is mounted on an aggregate block (custom-molded from Castinite, with 1/4-20 threaded inserts to match the laser
baseplate), which in turn is placed on a viscoelastic cushioning material (E-A-R Specialty Composites, Isodamp C-1002-06). A flexible silicone housing (custom-molded in-house from Smooth-On Mold Max 15T) encases the laser, reducing sensitivity to acoustic and thermal environmental perturbations.

The design described above and shown in Figs.~\ref{fig:cavity} and \ref{fig:diagram} is the
third  generation of this general design.  To minimize resonances and push them to high frequencies, we decreased the mass of the grating arm by removing more material and using a smaller diffraction grating.  We also added a Viton O-ring under the arm for shear damping (this replaced a magnet introduced for eddy-current damping, which was ineffective), and we thickened the lids from 1/2" to 3/4" to suppress lid-dependent resonances.  We also lowered the position of the precision screw to push in the vertical center of the grating arm, after finding that some coupling between vertical tilt and the desired horizontal adjustment was leading to mode-hopping and alignment problems. The valve was not machined into the front face of the laser in the previous versions: we attempted to shrink/press-fit the commercial aluminum seal-off valve
into a socket, but this did not seal consistently.  Finally, the original version lacked a Brewster window: the back facet of the first beam-shaping prism was glued to the exit window to create the seal, but this geometry would have made it difficult to use the same cavity body for a variety of laser wavelengths.

\section{Characterization}
Our diode-laser system is designed primarily with the goal of long-term stability in mind, but the reduction in sensitivity to acoustic perturbations ensures that the passive short-term stability and linewidth performance exceed published specifications for commercial, high-precision diode-laser sources.

The short-term stability was characterized in three ways: first, we measured the cavity resonances by performing resonance spectroscopy on the diffraction grating arm; second, we examined the power spectral density of the laser spectrum; third, we measured the delayed self-heterodyne beat note.

For comparison, the same measurements were performed on an older 780-nm ECDL constructed in our lab and used to cool $^{87}$Rb (Fig.~\ref{fig:old780}). This laser is a
descendant\cite{Steck01} of the well-known design of Ricci \textit{et al}.\cite{Ricci95}
Compared to the original design, the cavity length is increased, the diffraction angle is decreased, and the pivot of the diffraction grating arm is moved to the plane of the diode emission to optimize the mode-hop free tuning range.\cite{McNicholl85} Machined from a fatigue-resistant bronze alloy with high thermal conductivity, the laser is enclosed in an acrylic housing.  An auxiliary aluminum plate and a layer of Sorbothane are clamped directly to the aluminum baseplate for vibration damping.   Both the ``bronze'' 780-nm laser and our new test version at 780~nm use the same diode (Sharp GH0781JA2C) and collimating lens tube.

\begin{figure}[ht]
\centering
\includegraphics[width=8.5cm]{fig5.eps}
\caption{Older ``bronze'' diode laser in use in our lab: We compare this model to the new design in all stability measurements.}
\label{fig:old780}
\end{figure}

\subsection{Cavity resonance spectroscopy}
We investigated the susceptibility of the cavity to acoustic noise by the measuring the resonances of the most vibration-susceptible and critical element
for frequency stability---the diffraction grating arm.  An external Fabry--Perot cavity was slowly servo-locked (100-ms integration time constant) to keep the laser near the half-maximum point of the cavity resonance.  A small-amplitude (e.g., 0.03~V$_\mathrm{p-p}$ through a 20.2-dB attenuator) chirped sine wave was swept from 0 to 20~kHz over 20~s and applied to the laser tuning piezo. We monitored the error signal from the Fabry--Perot servo and plot the root-mean-square (rms) frequency excursions of the laser averaged over 200 sweeps in Fig.~\ref{fig:resonances}.

For comparison, we have included the same resonance measurements performed on the
bronze design described above.  At low frequencies, the modulation simply displaces the grating, tuning the laser directly; since this tuning response varies with design (the diffraction-grating arm mass, stiffness, material, and incident angle are all different), we scale the rms
excursions so that the high-frequency noise floors (which fall off at
6~dB/octave, as expected for a $1/f^2$ driven response) are consistent.
\begin{figure}[ht]
	\begin{center}
		\includegraphics[width=8.5cm]{fig6.eps}
	\end{center}
	\vspace{-3mm}
	\caption{%
Cavity resonance spectroscopy: rms error signal as a function of the driving frequency of the laser's piezo. Included are plots of the new design both with (solid red line) and without (dashed green line) shear damping as well as our lab's original bronze Littrow laser model (dotted blue line). In the inset, the same data are plotted on a logarithmic vertical scale.
	\label{fig:resonances}}
\end{figure}

Due to the stiff, lightweight grating arm,
we do not observe any significant resonances below around 12 kHz in our new design, and the addition of the damping
O-ring beneath the end of the grating arm effectively reduces the amplitude of the fundamental resonance without hampering the ability to tune the laser at low frequencies.
We confirmed that the visible peaks are associated with the grating arm by attaching various weights to the arm and observing the shifts of the resonances.  In some cases we have also observed peaks at 3.5 and 7~kHz whose amplitudes are nonlinear as a function of drive amplitude. These are evidence of an improperly secured grating and can be eliminated by using a sufficiently strong epoxy bond.

The previous version of our new design had a larger diffraction grating and
more mass in the grating arm, as described above; we found that decreasing the grating arm mass and inserting the O-ring under the arm for shear damping helped to minimize resonances and push them higher in frequency.  By reducing the mass of the arm and grating, we were able to increase the frequency of the lowest resonance by nearly 30\%.  We also observed an apparent coupling between the arm resonances and the lids in the previous version, which was eliminated by increasing the lid thickness.

\subsection{Frequency-noise spectrum}

Using a digital spectrum analyzer, we monitored the transmission of a 780-nm laser through a Rb vapor cell to measure the frequency-noise spectrum.  A slow, integrating servo (3~s time constant) maintained the laser frequency near the half-maximum point of a saturated-absorption feature ($F~{=}~3\longrightarrow F'~{=}~3,4$ hyperfine crossover) of
$^{85}$Rb.  The intensity fluctuations of the detected light are converted to frequency noise using the
well-known hyperfine structure.

We observed that the current controller used to power the laser had the largest effect on the overall baseline of the frequency-noise spectrum. For most lasers in our lab we use a home-built supply\cite{Meyrath03notes} with $\sim\!\!10$-$\mu$A rms noise,
but initial measurements demonstrated that one feature observed near 700~kHz could be identified as arising from noise in the current supply.  Low-pass filtering the current supply output decreased the frequency noise significantly, and a modification of the circuit design to increase the
sense-resistor value and reduce gain in the feedback loop further improved the characteristics.  A commercial (Vescent Photonics D2-105) version of the Libbrecht-Hall current controller\cite{Libbrecht93} offers superior noise performance (100-nA rms noise), and was used to power the lasers during all measurements presented here (except for the heterodyne measurement in Fig.~\ref{fig:780het}).

Figure~\ref{fig:PSD} shows the frequency-noise spectra for the two lasers, as well as a background noise floor (measured with the detector blocked), with each curve averaged for several minutes. Many features can be identified as electrical artifacts that appear in the background as well, but for the older bronze design, there are some peaks that correspond to the grating arm resonances in Fig.~\ref{fig:resonances} (e.g., at 1.8, 2.9, 5.5, and 6.6~kHz), while the $\sim\!\!14$~kHz grating arm resonance of the new design is barely visible.  No features were observed in the 20--100~kHz band in the spectrum of either laser.  The high-frequency baselines were similar for both the bronze laser and the new design, but the new design was much quieter in the acoustic (0--20~kHz) band.

Integration of the frequency-noise spectral density gives a measure of the laser's
rms frequency noise $\Delta\nu_{\rm{rms}}$ within a particular frequency band.  We compute
\begin{equation}
\Delta\nu_{\rm{rms}}~{=}~\eta\sqrt{R\int_{f_1}^{f_2}\overline{P}(f)\,df}
\end{equation}
where $\overline{P}(f)$ is the spectrum analyzer output in
watts (normalized by subtracting the background), $R$ is the 50~$\Omega$ internal impedance of the spectrum analyzer, and $\eta$ is the voltage-to-frequency conversion factor derived from the hyperfine spectrum.
Over the 20~Hz--100~kHz band, we find $\Delta\nu_{\rm{rms}}~{=}~296(26)$~kHz for the bronze laser and 146(8)~kHz for the new design (the lower 20-Hz cutoff was chosen to remove the contribution of the dc bias in the measurement, which produces a sharp dc peak broadened by the resolution bandwidth  of the analyzer).  The measured contribution from noise in the 20--100~kHz band was negligible for both lasers.   The contribution from noise above 100~kHz, computed by extrapolating the high-frequency baseline of the spectrum, also has a negligible effect (though this ignores the contribution of relaxation oscillations in the GHz range).

\begin{figure}[ht]
	\begin{center}
		\centering
		\includegraphics[width=8.5cm]{fig7.eps}
	\end{center}
	\vspace{-3mm}
	\caption{%
		Frequency-noise spectral density for the old bronze design (green), new stable design (red), and detector background (blue). Inset are the same data on a logarithmic scale from 100~Hz to 10~kHz.
	\label{fig:PSD}}
\end{figure}

As an extreme test of the lasers' response to environmental perturbations, we also subjected
them to a harsh environment by sounding an air horn from across the optical table (a distance of
4~m) whose intensity in each laser's vicinity was 105 dBA. The new laser responded very little, and we were able to use the same $F~{=}~3\longrightarrow F'~{=}~3,4$ saturated-absorption feature as a frequency reference.  However, the bronze laser was quite sensitive to the air horn, exhibiting frequency excursions that far exceeded the width of this feature.  Instead we used the side of the Doppler-broadened $F=3\longrightarrow F'$ collective resonance of $^{85}$Rb; this gave a reduced sensitivity, allowing us to measure the bronze laser's response. The reduced sensitivity caused the unperturbed laser to be indistinguishable from the noise floor, and leads to a high-frequency floor which (after conversion to a frequency-noise density) is offset from the floor of the crossover reference. We
attempted to use the same, Doppler-broadened feature with the new laser, but did not observe a significant response to the perturbation.

Figure~\ref{fig:airhorn} shows the laser frequency-noise spectrum as the air horn is sounded. Sensitivity to the perturbation is reduced in our new design by two orders of magnitude. The response of the $\sim\!\!14$~kHz grating-arm resonance can be seen, along with lower-frequency peaks, many of which are common to both lasers, and therefore reflect either details of the air-horn spectrum or acoustic room resonances.  To quantify the effect of the perturbation, we measure $\Delta\nu_{\rm{rms}}~\sim~1$~MHz for the new design, versus 40~MHz for the old design (an increase by factors of 7 and 135 relative to the non-perturbed values, respectively).

\begin{figure}[ht]
	\begin{center}
		\includegraphics[width=8.5cm]{fig8.eps}
	\end{center}
	\vspace{-3mm}
	\caption{%
Spectral response exhibited by the old bronze design (green) and the new stable design (red)
for an air-horn perturbation (105~dBA local sound intensity).  In the inset, the same data are plotted on a logarithmic vertical scale.
	\label{fig:airhorn}}
\end{figure}

During these measurements, we also used a more robust servo (with both proportional and integral feedback) to lock each laser explicitly to the crossover transition. We observed that common laboratory perturbations (e.g., bumping the optical table or laser and dropping ball drivers or Allen keys on the table) were occasionally catastrophic to the bronze laser's stability, but had little effect on our new design.

\subsection{Delayed self-heterodyne linewidth characterization}
 We characterized the lasers' short-term spectral width via the standard self-heterodyne
technique,\cite{Gallion82, Gallion84, Mercer91, Ludvigsen98} where
the laser output
splits into two components---one with a frequency shift, and one with a time delay---which
beat on a detector. The
resulting signal is recorded through a spectrum analyzer.
Our setup used delays in the range of
$\sim\!\!11$--44~$\mu$s, corresponding to 2--9~km of SMF-28 fiber,
and a nominal frequency shift of $\omega_\mathrm{m}/2\pi~{=}~80$~MHz, using first-order diffraction by an acousto-optic modulator.

To model the spectrum, we assume
a frequency-noise spectrum consisting of white-noise and $1/f$-noise components,\cite{Mercer91}
parameterized by $\gamma$ and $k$, respectively:
\begin{equation}
  S_\omega(\omega)~{=}~\gamma + \frac{k}{\omega}.
\end{equation}
Here, $S_\omega(\omega)$ is the one-sided spectral density of frequency fluctuations, which we may
define in terms of the laser phase $\phi(t)$ by
\begin{equation}
  S_\omega(\omega):=\int_{-\infty}^\infty\!\! \left\langle{\dot{\phi}(t)\,\dot{\phi}(t+\tau)} \right\rangle\cos\omega\tau\,d\tau .
\end{equation}

With this noise model, the normalized spectrum from the spectrum analyzer is given
by\cite{Mercer91}
\begin{align}
    s_\mathrm{\scriptscriptstyle{SH}}(\omega)
    =  \frac{1}{2\pi}\int_{-\infty}^\infty\!\! d\tau \, & \Big\{ \,\,\,e^{i\Delta\tau-\gamma\min\{|\tau|,|\tau_\mathrm{d}|\}}
    \label{bigfitfunc}\\
    &\times |\tau|^{k\tau^2/\pi}
    |\tau+\tau_\mathrm{d}|^{-k(\tau+\tau_\mathrm{d})^2/2\pi}\nonumber\\
    &\times|\tau_\mathrm{d}|^{k\tau_\mathrm{d}^2/\pi}
    |\tau-\tau_\mathrm{d}|^{-k(\tau-\tau_\mathrm{d})^2/2\pi}\nonumber\\
    &\times e^{-(\delta\omega_\mathrm{\scriptscriptstyle RB}\tau)^2/16\log 2
    -\delta\omega_\mathrm{drift} \tau/2} \,\,\Big\}, \nonumber
\end{align}
where $\Delta:=\omega-\omega_\mathrm{m}$, and $\tau_\mathrm{d}$ is the optical-fiber time delay.
We have
included two cutoffs in the integrand corresponding to the Gaussian resolution bandwidth
$\delta\omega_\mathrm{\scriptscriptstyle RB}$ (full width at half maximum) of the spectrum analyzer,
and also a phenomenological exponential cutoff $\delta\omega_\mathrm{drift}$ to model
drifts of $\omega_\mathrm{m}$ throughout the acquisition and averaging of the spectra over many minutes.
These cutoffs convolve the expected self-heterodyne spectrum with Gaussian and Lorentzian
functions, respectively. The resolution bandwidth
$\delta\omega_\mathrm{\scriptscriptstyle RB}/2\pi$ is nominally 1~kHz, and
$\delta\omega_\mathrm{drift}$ is of the same order or much smaller.

We measure the self-heterodyne spectra with multiple delays, and fit them to the the above spectral model,
taking $\gamma$, $k$, $\delta\omega_\mathrm{\scriptscriptstyle RB}$,
and
$\delta\omega_\mathrm{drift}$ to be fitting parameters.  We also take the overall scale and an overall
offset (corresponding to a detector noise floor of typically $-92$~dBm) as fitting parameters,
and we
allow for an extra fitting parameter, the amplitude of a
coherent ``spike'' (delta function) at $\Delta=0$ (in addition to any expected in the spectrum),
which is of course also broadened by
 $\delta\omega_\mathrm{\scriptscriptstyle RB}$ and
$\delta\omega_\mathrm{drift}$.
Thus,  $\delta\omega_\mathrm{\scriptscriptstyle RB}$,
$\delta\omega_\mathrm{drift}$, and the last amplitude serve to fit the center spike, while the other
parameters characterize the spectral tails.
The fiber delay $\tau_\mathrm{d}$ is also a fitting parameter, but only for the shortest fiber delays
where delay-dependent fringes are clearly visible in the spectral tails. This allows us to properly account for the wavelength-dependent fiber delay.

Once the $\gamma$ and $k$ parameters of the frequency-noise spectrum are determined, inferring the width of the
laser spectral line is somewhat subtle.  While the $\gamma$-dependent white noise leads to
a Lorentzian spectrum, the $k$-dependent, $1/f$ noise leads to a divergent spectral width.
Thus, the spectral width must be computed with respect to an observation time\cite{Cutler66}
$T\obs$,
since the spectral width diverges with $T\obs$.  (The spectral width also diverges
as $T\obs\longrightarrow 0$ owing to time-frequency uncertainty.)
To compute the spectral width, we thus numerically evaluate an
integral expression\cite{Steck06notes} analogous to Eq.~(\ref{bigfitfunc}) for the laser spectrum, incorporating the observation time, and
compute its full width at half maximum.
This can be computed for uniform observation of the signal over a time $T\obs$, or a Gaussian-weighted
observation in time, corresponding to a Gaussian resolution bandwidth of a
(hypothetical) spectrum analyzer for the optical frequencies.

While the frequency-noise spectrum characterization could only be performed for 780-nm model diode lasers due to the availability of the Rb vapor cells for calibration, we measured the delayed self-heterodyne signal for both 780-nm and 922-nm diode lasers, as well as for two diode lasers at 689~nm, one with the usual
2~cm cavity length and one with a longer, 10~cm cavity (otherwise identical to the 2~cm models).  Table~\ref{lasers_built} summarizes the laser parameters and measurement results.

\begin{figure*}[tb]
	\begin{center}
		\includegraphics[scale=1]{fig9.eps}
	\end{center}
	\vspace{-3mm}
	\caption{%
		Measured self-heterodyne spectra (solid green lines) and fits (dashed blue lines)
		to the model described in the text, for four different fiber-optic delay lengths, for the
		922~nm laser.
		Averaged over the four fits, the fitted spectral parameters are
		$\gamma/2\pi~{=}~1.86(29)$~kHz and $\smash{\sqrt{k}/2\pi}~{=}~27.0(15)$ kHz.
	\label{fig:922fit}}
\end{figure*}

The measured spectra for the 922-nm laser are plotted in
Fig.~\ref{fig:922fit} with the fitted spectra for four different delays.
The shortest (2.260~km) delay line gives a fitted delay of 11.116~$\mu$s;
the longest delay (8.872~km) thus corresponds to a 43.637-$\mu$s delay.
Averaging the fitted parameters from the four spectra, we find a white-noise width of
$\gamma/2\pi~{=}~1.86(29)$~kHz, and a $1/f$-noise parameter of $\smash{\sqrt{k}/2\pi}~{=}~ 27.0(15)$~kHz.
These parameters give spectral widths of
19, 32, and 75 kHz for $T\obs~{=}~50$, 100, and
1,000~$\mu$s, respectively.
Alternately, the spectral width is 40 or 71 kHz for a Gaussian resolution bandwidth of
$\delta\omega/2\pi~{=}~10$
or 1~kHz, respectively.

The measured spectra for the 780-nm models indicate a white-noise width of $\gamma/2\pi~{=}~1.8(3)$~kHz for the new laser and 15(1)~kHz for the older bronze laser, and $1/f$-noise parameters of
$\smash{\sqrt{k}}/2\pi~{=}~$
93.3(5) and 85(2)~kHz, respectively.  While the white-noise width of the new 780-nm laser was substantially narrower than the older bronze laser, and in fact was similar to that found for the new 922-nm laser, the $1/f$ parameter dominates and caused the two lasers to have
similar overall spectral widths of $\sim\!200$~kHz for $T\obs~{=}~100~\mu$s.  The narrower white-noise width of the new design is due partially to improved mechanical stability, but increased diffraction-grating resolution probably dominates.
The new design uses a grating with 1800 lines/mm, compared to 1200 lines/mm in the old bronze laser, and the incident angle is greater (45$^\circ$ versus 28$^\circ$), resulting in a larger illuminated grating area.

For comparison, we also examined the heterodyne spectrum (acquired with a 20-kHz Gaussian resolution bandwidth) between two identical new-model 780-nm lasers (Fig.~\ref{fig:780het}), with one laser powered by the Vescent current controller and the other by the improved, filtered home-built supply.  We fit the data with the heterodyne analog of Eq.~(\ref{bigfitfunc}),
and infer linewidth parameters of $\gamma/2\pi~{=}~20(10)$~kHz and $\smash{\sqrt{k}}/2\pi~{=}~$90(40)~kHz.  The uncertainty in these values is large due to the fact that averaging the heterodyne spectrum of free-running lasers is problematic, but the overall spectral width of $\sim\!\!200$~kHz for $T\obs~{=}~100~\mu$s is consistent within error with the self-heterodyne measurement.

\begin{figure}[ht]
	\begin{center}
		\includegraphics[width=8.5cm]{fig10.eps}
	\end{center}
	\vspace{-3mm}
	\caption{Heterodyne signal (solid green line) and fit (dashed blue line) between two identical 780~nm lasers, acquired using a 20-kHz Gaussian resolution bandwidth. The fitted spectral parameters are $\gamma/2\pi~{=}~20(10)$~kHz and $\smash{\sqrt{k}/2\pi}~{=}~90(40)$ kHz.
	\label{fig:780het}}
\end{figure}

The measured spectra for the 689-nm model with the same short 2~cm cavity length
indicates a linewidth  wider than either the 780 or 922-nm models: we find a white-noise width of $\gamma/2\pi~{=}~$93(36)~kHz and a $1/f$-noise parameter of $\smash{\sqrt{k}/2\pi}~{=}~$87.0(7) kHz, leading to an overall spectral width of 254~kHz for $T\obs~{=}~100~\mu$s.  Lengthening the cavity did, however, decrease the linewidth significantly: for the 10~cm cavity length, we find a white-noise width of $\gamma/2\pi~{=}~$4(1)~kHz and a $1/f$-noise parameter of $\smash{\sqrt{k}/2\pi}~{=}~$22(1) kHz, leading to an overall spectral width of 11.7~kHz for $T\obs~{=}~100~\mu$s.

The observed linewidth reduction due to the cavity-length extension is consistent with theory. The Schawlow-Townes linewidth is proportional to the cavity bandwidth squared,\cite{Schawlow58} so we expect the white-noise linewidths of the two lasers (which use identical diodes and have negligible front-facet reflectivity) to be related:
\begin{equation}
\frac{\gamma_{\rm{short}}}{\gamma_{\rm{long}}} = \left(\frac{nd+L_{\rm{long}}}{nd+L_{\rm{short}}}\right)^2,
\end{equation}
where $n$ is the semiconductor index of refraction, $d$ is the length of the solitary diode, and the $L_i$ are the lengths of the passive, extended cavities.
For typical diode parameters ({$nd\sim600~\mu$m}) and our external cavity lengths, we find that the measured linewidths are consistent with this relation.  Furthermore, the extended length reduces the sensitivity of the laser frequency to changes in diode index of refraction, and the above linewidth reduction is also equal to the square of the reduction in the frequency chirp $\partial\nu/\partial I$.\cite{Kazarinov87} We saw evidence of this effect, measuring current tuning values of 57(15) and 15(3)~MHz/mA for the short and long lasers, respectively, where the uncertainty
corresponds to systematic variations that depend on the overall current value and proximity to mode hops.

\begin{table*}
\begin{tabular}{|c|c||c|c||c|c||c|c|c|}
\hline
\multicolumn{2}{|c||}{Diode} & \multicolumn{2}{|c||}{Grating} & \multicolumn{2}{|c||}{Cavity} &  \multicolumn{3} {|c|}{Linewidth Parameters (kHz)}\\
\hline
\hline
\ \ wavelength \ \ & $R_{\rm{d}}$ & \ \ \ lines/mm\ \ \  & $R_{\rm{g}}$ & \ \ $L_{\rm{c}}$ (cm) \ \ & \ \ Design \ \ & \ \ \ \ \ $\gamma/2\pi$ \ \ \ \ \ & $\ \ \ \ \sqrt{k}/2\pi$\ \ \ \  &\ \ \ FWHM\ \ \  \\
 \hline
\ \ \ 780 nm \ \ \    &\ \ \  N/A \ \ \ & 1200 & \ \ 0.45 \ \ & 2.5 & old & 15(1) & 85(2) & 200\\
\ \ \ 780 nm \ \ \    &\ \ \  N/A \ \ \ & 1800 & \ \ 0.25 \ \ & 2.2 & new & 1.8(3) & 93.3(5) & 200\\

922 nm  &  10$^{-6}$  &  1500    &  0.19 & 2.2 & new & 1.86(29) & 27.0(15) & 32 \\
689 nm & 10$^{-5}$    & 2000   & 0.16  & 2.2 & new & 93(36) & 87.0(7) & 254 \\
689 nm & 10$^{-5}$    & 2000   & 0.16  & 10.0 & new & 4(1) & 22(1) & 11.7 \\
\hline
\end{tabular}
\caption{Summary of the lasers characterized by the self-heterodyne technique. $R_{\rm{d}}$ is the diode front-facet reflectivity, $R_{\rm{g}}$ is the measured diffraction grating efficiency, and $L_{\rm{c}}$ is the external cavity length. Also included are measured linewidth parameters (white-noise and $1/f$ components) and calculated FWHM for $T\obs$~{=}~100~$\mu$s. The diode and diffraction grating combination used for each wavelength are, respectively: old 780~nm---Sharp Microelectronics GH0781JA2C and Edmund NT43-751; new 780~nm---Sharp Microelectronics GH0781JA2C and Richardson 53-*-330H; 922~nm---Sacher SAL-920-60 and Richardson 53-*-239H; and 689~nm---Sacher SAL-690-25 and Richardson 53-*-059H.
\label{lasers_built}}
\end{table*}

There are several reasons why different diodes used in the same short laser cavity might display very different linewidths.  As a result of coupling between phase and intensity noise, the fundamental quantum Schawlow-Townes linewidth is known to be broadened in semiconductor lasers by a factor\cite{Henry82} that typically ranges from \mbox{3--7}.  While this factor is not often quoted by laser-diode manufacturers, since the technical $1/f$ noise component typically dominates ultimate linewidths, it could account for the increased white-noise parameter $\gamma$ observed in the 689-nm diode lasers.

Another factor to consider is that the diode-front-facet and diffraction-grating reflectivities for each model were very different:  The 780-nm diode lasers used diodes without a special front-facet anti-reflection (AR) coating, while the 922-nm and 689-nm diodes were AR-coated (with reflectivities $\le\!\!10^{-6}$ and $\le\!\!10^{-5}$, respectively).  The grating reflectivities were measured to be 0.45, 0.25, 0.19, and 0.16, respectively, for the gratings used in the old 780-nm model, the new 780-nm model, and the 922-nm and 689-nm models.  Although previous groups have found that grating resolution is a much more important factor than grating or front-facet reflectivities,\cite{Loh06, Saliba09} the relative resolution (ratio of wavelength to groove separation) of all of our gratings are similar in the new design, although the resolution of the old 780-nm model is lower.  However, the different reflectivity parameters suggest that competition between the external cavity modes and those of the diode itself may play a role in increasing the $1/f$-noise parameters of the new 780-nm and 689-nm models.

\section{Conclusion}\label{sec:conclusion}

We have constructed and characterised five
lasers using this new single-block design: two prototypes at 780~nm for comparison to existing homebuilt diode lasers, two more with the same cavity length at 922 and 689 nm, and one extra-long cavity (10-cm cavity length, as opposed to 2~cm for the shorter design) for a narrower passive linewidth (see Table~\ref{lasers_built} for a list of laser parameters and linewidth-measurement results).  In addition, we have also constructed two other lasers at 994 and 914 nm to complete the set of lasers needed for precision measurements using magneto-optically trapped, neutral strontium.
The need for so many different wavelengths, and for one ultra-narrow laser to probe strontium's ``forbidden''  $5\mathrm{s}^2\,{^1}\mathrm{S}_0\longrightarrow 5\mathrm{p}\,{^3}\mathrm{P}_1$
``clock'' transition at 689~nm, motivated our interest in designing a high-passive-stability diode-laser design adaptable to a wide range of diodes and cavity lengths.

Our diode laser design offers passive stability and linewidths superior to published specifications of commercial precision diode-laser systems.
For an observation time of $T\obs~{=}~100~\mu$s, we measured linewidths of 32~kHz for an AR-coated diode at 922~nm, and $\sim$200~kHz for a non-AR coated diode at 780 nm.
For AR-coated diodes at 689~nm in the regular and the extra-long cavities, we measured linewidths of 254~kHz and 11.7~kHz, respectively.  We have also demonstrated that our laser exhibits no low-frequency mechanical resonances and is extremely insensitive to acoustic noise.

For comparison, the New Focus flagship ECDL, the Vortex II (TLB-6900 series), displays several resonances under 10~kHz (with the largest feature located at $\sim $500 Hz), and specifies
a linewidth of $\le\!\!300$~kHz from a simple full-width-half-max (FWHM) measurement of the heterodyne spectrum.\cite{NewFocus11brochure}  The Toptica DL Pro does not provide mechanical-vibration measurements, and specifies a 100~kHz typical linewidth, associated with a 5~$\mu$s time figure whose meaning is not stated.\cite{Toptica11brochure} Toptica also presents a Lorentzian fit to a heterodyne spectrum for two of these lasers, with a FWHM of 300~kHz, deconvolved to $150$~kHz;\cite{Toptica11cat} however, the spectrum itself is not Lorentzian, and this figure underestimates the deconvolved linewidth by a significant margin.  Vescent Photonics also uses a Lorentzian fit to characterize a non-Lorentzian spectrum, similarly underestimating to specify a fast linewidth on the order of 50 kHz for their Chip External-Cavity Laser (CECL).\cite{VescentCECL, VescentCECL2} Finally, Sacher's Lynx ECDL has a typical linewidth of 1~MHz for $T\obs~{=}~1~$ms.\cite{SacherLynx}

A set of components for a diode-laser system using our design costs under
\$6,000. Machining, materials, and a set of peripherals (including diffraction grating and fiber coupler)
 costs under \$3,000.  In addition, an appropriate diode and optical isolator must be purchased; costs for these components range from \$1,000--\$3,000, depending on diode wavelength, AR-coating, isolator model, etc.  In our lab this past summer, a talented undergraduate student---who had no prior laser-assembly experience---assembled and aligned four laser systems in under a month.
All machine drawings and CAD files, detailed assembly instructions, and all circuit schematics and board layouts for the relevant electronics are available on our project web site.\cite{Unilaser_web}

Optics (physics.optics); Atomic Physics (physics.atom-ph); Instrumentation and Detectors (physics.ins-det)


Correcting for the Effects of Interstellar Extinction

Abstract: This paper addresses the issue of how best to correct astronomical data for the wavelength-dependent effects of Galactic interstellar extinction. The main general features of extinction from the IR through the UV are reviewed, along with the nature of observed spatial variations. The enormous range of extinction properties found in the Galaxy, particularly in the UV spectral region, is illustrated. Fortunately, there are some tight constraints on the wavelength dependence of extinction and some general correlations between extinction curve shape and interstellar environment. These relationships provide some guidance for correcting data for the effects of extinction. Several strategies for dereddening are discussed along with estimates of the uncertainties inherent in each method. In the Appendix, a new derivation of the wavelength dependence of an average Galactic extinction curve from the IR through the UV is presented, along with a new estimate of how this extinction law varies with the parameter R = A(V)/E(B-V). These curves represent the true monochromatic wavelength dependence of extinction and, as such, are suitable for dereddening IR--UV spectrophotometric data of any resolution, and can be used to derive extinction relations for any photometry system.

\section{Introduction}

A precise knowledge of the wavelength dependence of interstellar
extinction (i.e., the absorption and scattering of light by
interstellar dust grains) and of any spatial variability in this
dependence is important for two distinct reasons.  First, the
extinction depends on the optical properties of the dust grains along a
line-of-sight and potentially can reveal information about the
composition and size distribution of the grains.  Further, changes in
the extinction from place to place may reveal the degree and nature of
dust grain processing occurring in the ISM.  Second, the wavelength
dependence of extinction is required to remove the effects of dust
obscuration from observed energy distributions, since most astronomical
objects are viewed through at least some small amount of interstellar
dust.  Spatial variations in the extinction potentially limit the
accuracy to which energy distributions can be ``dereddened.'' Such
uncertainties might be acceptably small for very lightly reddened
objects, but quickly can become debilitating along modestly reddened
sightlines.

This paper follows from a talk given at the meeting ``Ultraviolet
Astrophysics Beyond the IUE Final Archive'' (Fitzpatrick 1998) and
addresses the second point raised above, i.e., the correction of energy
distributions for the effects of interstellar extinction.  The goals
are to provide a summary of what is known about the wavelength
dependence and spatial variability of interstellar extinction and to
present strategies for the best removal of the effects of extinction
from astronomical data.  In \S 2, I briefly review the chief features
of IR-through-UV extinction and the nature of the known spatial
variations.  These variations are illustrated with a set of optical/UV
extinction curves which span the (currently) known extreme limits of
extinction variations.  Constraints on the wavelength dependence of
extinction and some general correlations between extinction curve shape
and interstellar environment are also noted in \S 2.  Strategies for
dereddening are discussed in \S 3 along with estimates of the
appropriate uncertainties which ought to be incorporated into an error
analysis.  Some final comments regarding extinction along extragalactic
sightlines and at wavelengths shortward of 1150 \AA\/ are given in \S
4.

The Appendix of this paper describes the construction of a new estimate
of the shape of a mean and an environment-dependent IR-through-UV
extinction curve.  These curves, which may be obtained from the author,
are suitable for dereddening multiwavelength spectrophotometric data,
such as has become available with the {\it Hubble Telescope's} Faint
Object Spectrograph, and for deriving extinction relationships for any
photometric system.

\section{Spatial Variations in Galactic Extinction}
\subsection{Variations Abound ...}

Figure 1 shows several estimates of the shape of a ``mean'' Galactic
extinction curve from the far-IR through the UV region.  The curves are
presented in a commonly used normalization scheme,
$E(\lambda-V)/E(B-V)$, and plotted against inverse wavelength.  The
significance of the different curves will be discussed later in this
paper.  For the present, these curves serve to illustrate the overall
characteristics of interstellar dust extinction.  Typically, the
extinction rises through the IR with a power law-like  dependence (see
Appendix), rolls over slightly in the optical region (``knee''), shows
a prominent feature at 2175 \AA\/ in the near-UV (``bump''), and has a
sometimes steep rise in the far-UV (``fuv rise'').  Figure 1 does not
attempt to show the relatively narrow optical/IR features which may be
associated with interstellar dust, such as the diffuse interstellar
bands in the optical or the solid state absorption (or emission)
features in the IR (e.g., the 9.8 \micron\/ silicate feature).

IR and optical extinction traditionally have been studied using ground
based broad- or intermediate- band photometry and
sightline-to-sightline variations have been recognized for a long time
(e.g., Johnson 1965).  These are often characterized using the
parameter $R \equiv A(V)/E(B-V)$, i.e., the ratio of total to selective
extinction at $V$.  The value of $R$ ranges between about 2.2 and 5.8
for sightlines along which UV extinction has also been measured, with a
mean value of $\sim$3.1 for the diffuse interstellar medium.  It is
believed that the shape of the IR extinction law longward of $\sim$8000
\AA\/ (the ``power law'' region) may be invariant and that the observed
range in $R$ is due to spatial variations in the steepness of the
extinction in the optical region around the ``knee''  (e.g., Cardelli,
Clayton, \& Mathis 1989;  Martin \& Whittet 1990).

Measurements of UV interstellar extinction have relied on satellite or
rocket spectrophotometric data, and date back to the original discovery
of the 2175 \AA\/ bump nearly 30 years ago (Stecher 1969; Bless \&
Savage 1970).  These studies immediately revealed that
sightline-to-sightline differences exist within the Galaxy in the
detailed properties of UV extinction curves.  Important results on the
degree of spatial variability of UV extinction were subsequently
derived from the large photometric databases accumulated by the {\it
TD-1} and {\it ANS} satellites (e.g.\, Koornneef 1978; Kester 1981;
Meyer \& Savage 1981).

The most graphic illustrations of spatial variations in UV extinction
were provided by the {\it IUE} satellite.  With its complete spectral
coverage between 1150 and 3200 \AA\/ and relatively high spectral
resolution  ($\sim$6 \AA\/ in the low-resolution mode) {\it IUE} was
well-suited for the study of UV extinction and many studies have
utilized {\it IUE} data to determine the wavelength dependence of
extinction along one or more Milky Way lines-of-sight (e.g.\, Bohlin \&
Savage 1980; Witt, Bohlin, \& Stecher 1980; Seab, Snow, \& Joseph 1981;
Morgan, McLachlan, \& Nandy 1982; Hecht et al. 1982; Massa, Savage, \&
Fitzpatrick 1983; Witt, Bohlin, \& Stecher 1984; Massa \& Savage 1984;
Franco, Magazz\`{u}, \& Stalio 1985; Aiello et al. 1988; and many
others).   Figure 2 shows analytical representations for $\sim$80 UV
extinction curves derived from {\it IUE\/} data (solid curves),
illustrating the wide range of properties observed for Galactic UV
extinction.  These curves are taken from the catalog of Fitzpatrick \&
Massa (1990), with a few additions, including the lines-of-sight toward
HD 210121 (Welty \& Fowler 1992) and HD 62542 (Cardelli \& Savage
1988).  The basis for the analytical representation will be discussed
further below.

The overall impression that has come from the last 30 years of
extinction studies is that there is a bewildering variety of
IR-through-UV extinction curves and that what has often been called
``peculiar'' extinction is actually quite common.  This bodes poorly
for the prospects of precisely dereddening astronomical energy
distributions.  In particular, the range of extinction properties seen
in the UV (Figure 2), where the effects of extinction are strongest,
may represent a bonanza of information for those who study the
properties of interstellar dust grains, but is a potential disaster for
those whose main interest in UV extinction is the removal of its
effects.

\subsection{... But There is Some Order Amid the Chaos}

The recognition that interstellar extinction is spatially highly
variable is critically important for both the study of dust grains and
for the purposes of dereddening.  Fortunately for the latter, it has
been realized that within the seemingly endless variety of observed
extinction curves there are actually some constraints on the wavelength
dependence and some links between extinction in the UV, optical, and
IR.

This began with Savage (1975) who showed, using {\it OAO-2} satellite
data, that the shape of 2175 \AA\/ bump can be reproduced well with a
Lorentzian profile.  This result --- which provides important clues
about the nature of the 2175 \AA\/ feature --- allows a functional
representation of the bump, and permitted Seaton (1979) to present an
estimate of the average Galactic UV extinction curve as a simple
analytical formula.  This curve is shown in Figure 1 by the dashed
curve, extending from 2.7 to 8.7 \micron$^{-1}$.

Following the lead of Savage and Seaton,  Fitzpatrick \& Massa (1986,
1988, 1990; collectively referred to hereafter as FM) examined {\it
IUE\/} extinction curves for many lines-of-sight and found that all
these curves could be fitted extremely well by a single analytical
expression with a small number of free parameters.  This expression
consists of (1) a Lorentzian-like bump term (requiring three
parameters, corresponding to bump width $\gamma$, position $x0$, and
strength $c3$), (2) a far-UV curvature term (one parameter $c4$), and
(3) a linear term underlying the bump and the far-UV (two parameters
$c1$ and $c2$).  This set of basis functions is shown in Figure 3.
With the proper choice of the 6 parameters, essentially all UV
extinction curves can be reproduced to within the uncertainties
inherent in the data.  For most sightlines, this number of free
parameters can be reduced to 4, since it has been seen that the
position of the 2175 \AA\/ bump is nearly invariant and that the
intercept and slope of the linear term are tightly correlated (FM; see
also Carnochan 1986 and Jenniskens \& Greenberg 1993), such that the
linear term appears to pivot about the point
[$1/\lambda,\frac{E(\lambda-V)}{E(B-V)}$] $\simeq$ [3,2].  The UV
portions of the curves shown in Figure 2 ($\lambda < 2700$ \AA) were
all constructed using these functions.  There is, thus far, no convincing
evidence for substructure in UV extinction curves.

The FM results show that some order underlies the wide spectrum of
wavelength dependences seen in the UV, but do not necessarily permit a
more accurate dereddening of an object along any given line of sight.
Fortunately it has been observed that extinction curves do not vary
randomly across the sky, but show distinct regional signatures (e.g.,
Meyer \& Savage 1981; Morgan, McLachlan, \& Nandy 1982; Panek 1983;
Massa \& Savage 1984; Clayton \& Fitzpatrick 1987) and it is clear that
extinction properties reflect the physical conditions and past
processing histories of the environments in which the dust grains
reside.

The work of Cardelli, Clayton, \& Mathis (1988, 1989; hereafter CCM)
provides a link between one measure of dust grain environment and the
wavelength dependence of UV-through-IR extinction.  CCM found that the
shape of UV extinction curves correlates with the parameter $R$
(defined in \S 2.1 above).  This suggests that --- although there is
considerable scatter --- extinction curves from the UV through the IR
can be characterized as approximately a one-parameter family dependent
on $R$.  Thus, if the value of $R$ can be determined (from optical and
IR photometry), then the properties of the the entire UV--IR extinction
curve can be predicted.  The essence of the CCM result is illustrated
in Figure 4.  Four representative extinction curves are shown, each
determined by the value of $R$ listed at the righthand side of the
curves.   This figure shows that extinction curves which are seen to be
very ``flat,'' or ``grey,'' in the UV roll over strongly in the optical
and are characterized by large values of $R$.  Steep UV curves remain
steep in the optical and are characterized by small values of $R$.  The
CCM curve for $R=3.1$ is shown in Figure 1 by the dotted line.

The results of CCM are important in 3 ways:  (1) They indicate that
some of the spatial variations seen in extinction curves behave
coherently and systematically over a wide wavelength range, thus
potentially allowing for the consistent dereddening of multiwavelength
energy distributions.  (2) The observed dependence on $R$ allows for
the definition of a meaningful average Galactic extinction curve.
Existing datasets of extinction curves, particularly for the UV, are
biased toward ``interesting'' regions where extinction properties are
extreme.  Computing a simple mean curve from these data does not
necessarily yield a reasonable estimate of a Galactic mean.  Since it
is well-established, however,  that the appropriate mean value of $R$
for the diffuse ISM is $\sim$3.1 (e.g., Schultz \& Wiemer 1975; Whittet
\& van Breda 1980; Rieke \& Lebofsky 1985), a reasonable definition of
a mean Galactic extinction law is that which corresponds to the case
where $R = 3.1$.  (3)  They demonstrate a general correlation between
dust grain environment and the wavelength dependence of extinction,
since large values of $R$ are generally found in dense environments
where dust grain growth is thought to occur.  The ``greyness'' of the
UV/optical extinction curves for large $R$ is consistent with a larger
than normal dust grain population.   These results have an important
bearing on how to correct for the effects of extinction, discussed in
the following section.

\section{Correcting for the Effects of Extinction}

Given the observed complexity of Galactic UV extinction, what is the
best way to deredden an observed energy distribution?   There are three
different possibilities, which involve the use of a global mean
extinction curve, an $R$-dependent curve, and a sightline-specific
curve.  These three alternatives and their inherent uncertainties are
discussed below, and a new derivation of the shape of the average
IR-through-UV extinction curve is presented.  A more thorough
discussion of the uncertainties involved in UV extinction corrections
is given by Massa (1987).

\subsection{Using a Global Mean Curve}

If there is no specific information available about the wavelength
dependence of the extinction along a sightline of interest then the
only alternative is to adopt some globally defined mean curve and
perform a realistic error analysis.  This is the least attractive of
the three cases discussed here, but is by far the most commonly used.
The average Galactic extinction curves from either Seaton (1979; see
Figure 1) or Savage \& Mathis (1979) are often adopted for dereddening
UV data.  These should now be superseded by the use of an $R$-dependent
curve computed for the case $R = 3.1$ (see \S 2.2).  Figure 1 shows the
$R = 3.1$ curve from CCM (dotted line) and a new determination of the
$R = 3.1$ case (solid line) which is described in the Appendix to this
paper.  This new curve aims to reproduce the detailed wavelength
dependence of the $R = 3.1$ extinction law and has been constructed to
account for bandpass effects properly in optical/IR extinction data and
to reproduce the observed broad-, intermediate-, and narrow- band
extinction measurements.  The curve is thus suitable for dereddening
multiwavelength spectrophotometric observations and can be used to
derive the average extinction relationships in any photometric system.

A complete evaluation of the likely error in a dereddened relative
energy distribution $m(\lambda-V)$ requires the propagation of the
uncertainty in $E(B-V)$ and the often-neglected but often-dominant
uncertainty in the adopted mean curve.  This is given by
\begin{equation}
\sigma^2_{m(\lambda-V)} =  \sigma^2_{E(B-V)} \times k(\lambda-V)^2 + E(B-V)^2 \times \sigma^2_{k(\lambda-V)}
\end{equation}
where $k(\lambda-V)$ represents the normalized extinction curve
$E(\lambda-V)/E(B-V)$.  If absolute fluxes are desired, then the
uncertainty in the assumed value of $R$ must also be incorporated (two
additional terms on the righthand side of eq. 1, similar to the current
terms but with $R$ substituted for $k(\lambda-V)$).

The large database of {\it ANS} satellite extinction measurements from
Savage et al. (1985) shows that the 1-$\sigma$ scatter at 1500 \AA\/ is
$\sigma_{k(15-V)} = 0.74$, based on $\sim$400 sightlines with $E(B-V)
\ge 0.5$.  (Only relatively large values of $E(B-V)$ were considered in
order to minimize the effects of spectral mismatch error and random
noise on the $\sigma$ measurement.) To extend this analysis to other
wavelengths, we computed the standard deviation at each wavelength
point for the 80 curves shown in Figure 2, and then scaled the result
to match the {\it ANS} value at 1500 \AA.  (The actual standard
deviation at 1500 \AA\/ for the 80 {\it IUE} curves is somewhat higher
than the {\it ANS} result because of the bias in Figure 2 toward
extreme extinction curves).  The resultant values of
$\sigma_{k(\lambda-V)}$ are shown by the thick dotted curve near the
bottom of Figure 2 (labeled ``$\sigma$'') and are listed at selected
wavelengths in the third column of Table 1.   This estimate of
$\sigma_{k(\lambda-V)}$ should be adopted whenever the average Galactic
extinction curve is used for dereddening an observed energy
distribution.  The uncertainties approach zero for $1/\lambda < 3$ $\mu
m^{-1}$ due to the curve normalization.  The quantity $E(B-V)$ is
usually derived directly from photometry and often has an easily
quantifiable uncertainty; thus the computation of
$\sigma^2_{m(\lambda-V)}$ from eq. 1 is straightforward.

Sometimes $E(B-V)$ is not known {\it a priori}\/ and is estimated  by
``ironing out'' the 2175 \AA\/ extinction bump using an assumed
extinction curve shape (e.g., see Massa, Savage, \& Fitzpatrick 1983).
In many cases, uncertainties in $E(B-V)$ as small as 0.01-0.02 mag for
moderately reddened objects have been quoted from this process.  This
is incorrect!  The normalized ``height'' of the 2175 \AA\/ extinction
bump --- which is the quantity that is important in the ironing-out
process --- has a 1-$\sigma$ scatter of about $\pm$20\% around its mean
value (from the data in Figure 2 and Savage et al. 1985).  Thus the
uncertainty in an $E(B-V)$ measurement derived from the bump must be
considered to have a similar relative uncertainty.  Since bump strength
does not correlate well with other aspects of UV extinction, such as
the slope of the linear component or the strength of the far-UV rise
(FM), the uncertainty in an energy distribution dereddened this way can
be estimated by using eq.  1 with $\sigma_{E(B-V)} \simeq 20\%$ and the
values of $\sigma_{k(\lambda-V)}$ from Table 1.  The result of this
calculation needs to be modified somewhat by removing the signature of
the 2175 \AA\/ bump.  The final estimate of $\sigma^2_{m(\lambda-V)}$
computed this way is listed in the last column of Table 1.  Multiplying
these values by $E(B-V)$ gives the uncertainty in a dereddened relative
energy distribution.  The values in Table 1 agree well with estimates
$\sigma^2_{m(\lambda-V)}$ at 1500, 1800, 2200, 2500, and 3300 \AA\/
derived from the {\it ANS} data of Savage et al (1985), and are
consistent with the analysis of Massa (1987).

\subsection{Using a $R$-Dependent Mean Curve}

If the value of $R$ has been measured along a sightline of interest, or
if some general information is available about the dust grain
environment (e.g., dense region or diffuse region?), then the
uncertainties in dereddening can be reduced somewhat by adopting an
$R$-dependent extinction curve computed at the appropriate value of
$R$.  For the original 29 sightlines used by CCM, the standard
deviation of the individual observed extinction curves minus the
computed $R$-dependent curves is about a factor of 2 smaller than when
the individual curves are compared to the average Galactic ($R = 3.1$)
curve.  Thus the uncertainties in a dereddened energy distribution can
be estimated using eq. 1 and adopting 0.5 times the
$\sigma_{k(\lambda-V)}$ values from Table 1.  If $R$ is poorly measured
or only estimated from environmental factors, then allowance must be
made for this in the error analysis.

Using an appropriate $R$-dependent curve can slightly improve the
accuracy of results from ironing out the 2175 \AA\/ bump.  The values
of $\sigma^2_{m(\lambda-V)}$ can be estimated as above in \S 3.1, using
the same value of $\sigma_{E(B-V)} \simeq 20\%$ and 0.5 $\times$
$\sigma_{k(\lambda-V)}$ from Table 1.  The resultant values of
$\sigma^2_{m(\lambda-V)}$ are about 85\% as large as those listed in
the last column of Table 1.  The gain from using the $R$-dependent
curves is only marginal because the normalized bump heights are highly
variable and not correlated with $R$.

It is important to note that, even when $R$ and $E(B-V)$ are
well-determined, the use of an $R$-dependent extinction curve only
reduces --- but does not eliminate --- the wavelength dependent
dereddening error.   The standard deviation of the observed vs.
$R$-dependent curves does not go to zero, even for the sample of 29
sightlines used by CCM to define the $R$-dependence, because the
UV/optical curves are not really a one-parameter family dependent only on
$R$.   Curves derived for different sightlines with the same values of
$R$ show a wide range of properties, including differences in the
strengths of the bump and the far-UV rise.  Thus, as noted by CCM, the
CCM formula reproduces a general trend, but does not provide
particularly good fits to individual extinction curves, even when the
value of $R$ is well-determined.

Apart from this intrinsic scatter around the mean $R$ relation, the
accuracy of the CCM results is limited in the IR/optical region due to
bandpass effects with the broadband Johnson filters used to measure the
extinction.   As a result, CCM tends to overestimate the level of
extinction in the near-IR and blue-visible.  In the Appendix, a new
derivation of the $R$-dependence of IR-through-UV extinction is
presented.  This corrects for the systematic bandpass effects in CCM,
but is still plagued by the same uncertainties caused by the ``cosmic''
scatter of extinction properties around the mean $R$ relation.

\subsection{Using a Sightline-Specific Curve}

Without question, the best of all possible dereddening choices is to
use an extinction curve derived explicitly for the sightline of
interest.  This advice is not so unhelpful as it may sound, because the
archives of the {\it IUE\/} satellite  may well contain normal stars in
the direction (and at the distance) of the object of interest, from
which extinction curves may be derived (e.g., halo stars in the case of
extragalactic sightlines).  The accuracy of the dereddened energy
distribution may thus be limited chiefly by the uncertainties in the
``pair method'' of constructing the extinction curve (see, e.g., Massa
et al. 1983) and by the determination of $E(B-V)$ for the object of
interest.  For a well-determined curve, ``ironing out'' the bump --- in
the cases when $E(B-V)$ is not known explicitly --- should yield the
same degree of accuracy as for the cases when $E(B-V)$ is known {\it a
priori}.  An example of such an approach is the determination of the
intrinsic energy distributions of rare or exotic objects located in
clusters with more mundane objects.  Massa \& Savage (1985) utilize
this technique, using extinction curves derived from main sequence B
stars in open clusters to determine the continuum shapes of cluster O
stars.

\section{Final Comments}

The results discussed above strictly pertain to Milky Way extinction
for wavelengths longward of about 1150 \AA\/ (the {\it IUE} satellite
cutoff).  Observers of extragalactic objects must contend with
foreground Milky Way extinction as well as extinction due to dust
grains local to the object of interest.  It is reasonable to assume
that the properties of interstellar extinction in external galaxies are
fully as complex as found for the Milky Way.  However, only for the
nearby Large Magellanic Cloud (LMC) do observed extinction curves for
individual sightlines approach the accuracy needed to study spatial
extinction variations (see Fitzpatrick 1998 and references within).
The properties of Milky Way extinction along the high galactic latitude
sightlines relevant to most extragalactic observations are also not
well-determined.  One of the few quantitative studies of high latitude
extinction, that by Kiszkurno-Koziej \& Lequeux (1987) utilizing {\it
ANS} satellite data, suggests that halo extinction is slightly steeper
in the far-UV and has a weaker 2175 \AA\/ bump than the typical
extinction found in the plane of the Galaxy.  Attempts to eliminate
galactic foreground reddening by ironing out the bump may lead to
systematic underestimates of the total amount of foreground
extinction.

The properties of extinction in the wavelength range 912-1150 \AA\/ are
poorly known, but the work that has been done (e.g., York et al. 1973;
Snow, Allen, \& Polidan 1990) suggests that the far-UV rise continues
at least down to wavelengths of 950 \AA\/ or so.  An intrinsic problem
with measuring extinction in this wavelength range is the high density
of H I and H$_2$ absorption lines which mimic continuous absorption in
low-resolution data.  CCM suggest that the wavelength dependence of
extinction in this region is reasonably well-represented by an
extrapolation of the fitting functions used in the {\it IUE} wavelength
range.  At the present, the best dereddening strategy for 912-1150
\AA\/ data would be to adopt such an extrapolation, based if possible
on {\it IUE} extinction curves derived for the sightline of interest.
The uncertainties listed in Table 1 for 1000 \AA\/ are based on this
simple extrapolation.

\acknowledgements

I thank my partner-in-interstellar-grime, Derck Massa, for many years
of enjoyable collaboration.  I also acknowledge the contributions of
the late Jason Cardelli to the study of interstellar extinction.
Future progress in the field will proceed at a slower pace without
Jason's enthusiasm and scientific insight.  This work was supported in
part by NASA grant no. NAG5-7113.

\begin{appendix}
\section{Deriving a $R$-Dependent IR-through-UV Extinction Curve}

In this Appendix a new estimate of the wavelength dependence of Milky
Way extinction from the IR through the UV is derived.  The motivation
for this is the need for a multiwavelength extinction curve suitable
for dereddening IR-through-UV spectrophotometric data now available
from, for example, the  {\it Hubble Space Telescope's} Faint Object
Spectrograph (see, e.g., Fitzpatrick \& Massa 1998 and Guinan et al.
1998).  The specific goals are 1) to produce a curve for the $R=3.1$
case which reproduces existing narrow band extinction measurements {\it
and} the principal Johnson and Str\"{o}mgren photometric extinction
relations and 2) to define the $R$-dependence of the extinction law,
relying principally on spectral regions where this dependence is
well-determined.   Existing estimates of the average wavelength
dependence of optical/near-IR extinction do not accurately portray the
true shape of the curve due to bandpass effects arising from the
broadband photometry used to derive the curves.  This deficiency will
be corrected by using synthetic photometry.

\subsection{The UV Region}

CCM noted a correlation between the value of $R^{-1}$ (as determined
from IR/optical photometry) and the values of $A(\lambda)/A(V)$ (i.e.,
total extinction at $\lambda$ normalized by total extinction at V) at
UV wavelengths using parametrized extinction curves for a subset of 29
stars from the FM sample.  They consequently derived a complex
polynomial expression to reproduce the wavelength- and $R$-
dependences.  We approach the problem slightly differently and note
that the essence of the CCM result is two correlations involving the
coefficients of the FM fitting function originally used to parametrize
the extinction curves studied by CCM.  We thus derive the functional
forms of these correlations and retain the FM fitting function to
compute the values of $E(\lambda-V)/E(B-V)$ for wavelengths $\lambda <
2700$ \AA.

Figure 5 illustrates the two correlations.  The top panel shows a
plot of the FM parameter $c2$ (representing the slope of the linear UV
extinction component) vs. $R^{-1}$ (filled circles) along with a least
squares estimate of the linear relationship between the two (dashed
line).  The data are for a subset of the FM sample for which IR
photometry could be used to deduce the value of $R$ (see
below),  with the addition of the sightline toward HD 210121 (Welty \&
Fowler 1992; Larson, Whittet, \& Hough 1996).  This represents the only
convincing correlation between $R$ and the properties of the UV
extinction curve (see, e.g., Jenniskens \& Greenberg 1993).  The
equation of the least-squares fit to the data is
\begin{equation}
c2 = -0.824 + 4.717\times R^{-1}             \\
\end{equation}
No physical interpretation or significance is placed on the form of the
adopted functional relationship between $R$ and $c2$; it is merely that
which best reproduces the observed correlation.  We note, however, that
the intrinsic slope of extinction in the optical region is proportional
to $R^{-1}$, and that eq. A1 thus simply implies that the slopes of the
extinction curves in the UV and optical vary together in a linear
manner --- as optical extinction steepens, UV extinction steepens.

The bottom panel of Figure 5 shows the well-known relationship between
slopes ($c2$) and intercepts ($c1$) of the linear background component
(FM; Carnochan 1986; Jenniskens \& Greenberg 1993).  The data are for
the full set of 80 curves from FM, plus HD 210121, and minus the Orion
Nebula stars (HD 36982, 37022, 37023, and 37061) which suffer from
scattered light contamination (which mainly affects the linear intercept).
The equation of the least squares fit to this relationship (dashed
line), is
\begin{equation}
c1 =  2.030 - 3.007 \times c2                       \\
\end{equation}

We use the following mean values for the other four parameters required
to specify UV extinction curves with the FM formula:  $x0$ (bump
position) = 4.596 \micron$^{-1}$; $\gamma$ (bump width) = 0.99
\micron$^{-1}$; $c3$ (bump strength) = 3.23; and $c4$ (FUV curvature) =
0.41.

\subsection{The Optical/IR Region}

A common way of representing the wavelength dependence of extinction in
the optical and IR regions is to draw a smooth curve through the values
of $E(\lambda-V)/E(B-V)$ (or some other curve normalization) derived
from Johnson {\it UBVRIJHKLM} photometry, assuming that these represent
the monochromatic values of the extinction at the filter effective
wavelengths $\lambda_{eff}$ (e.g., CCM, Martin \& Whittet 1990).  This
procedure {\it does not} accurately yield the monochromatic extinction
curve because it ignores the wavelength dependence of extinction across
the width of an individual filter.  Extinction decreases towards longer
wavelengths in the optical/IR region and therefore the value of
$\lambda_{eff}$ for a reddened star is shifted to a longer wavelength
than for an identical unreddened star.  An extinction measurement made
by comparing photometric indices for two such stars (i.e., the ``pair
method'') will always overestimate the monochromatic extinction in the
neighborhood of $\lambda_{eff}$ because stellar emergent fluxes in the
optical/IR decrease towards longer wavelengths (for early-type stars).
The magnitude of this effect depends on how much the true extinction
curve varies across the filter.  A broadband $E(\lambda-V)/E(B-V)$
measurement will also be influenced by the intrinsic energy
distributions of the stars used in the pair method and by the total
amount of extinction.

The approach taken here is to find the wavelength dependence of the
optical/IR extinction curve which reproduces the photometric extinction
measurements when synthetic photometry of an artificially reddened
stellar energy distribution is compared with that of identical but
unreddened star.  To represent the stellar energy distribution, an
ATLAS9 model atmosphere from R.L.  Kurucz (with $T_{eff} = 30000$ K and
$\log g = 4.0$) is used, along with an adopted value of $E(B-V) =
0.5$.  Synthetic photometry is performed to yield measurements in the
broadband Johnson {\it UBVRIJHKLM} system and the intermediate
band Str\"{o}mgren {\it uvby} system.  The first step is to determine the
curve shape for the mean $R=3.1$ case and then to define the nature of
the variation with $R$.  The first two columns of Table 2 list the
photometric extinction ratios (in both the Johnson and Str\"{o}mgren
systems) and their observed values, which the $R=3.1$ extinction curve
is required to reproduce.  References to the observations are given in
Column 3.  In addition, the curve is constrained to reproduce the mean
narrowband measurements in the 3400--7900 \AA\/ region published by
Bastiaansen (1992), which are assumed to represent the monochromatic
values of the extinction for the case $R = 3.1$ at the filter central
wavelengths.

Figure 6 shows the shape of the monochromatic extinction curve which
best satisfies these requirements (thick solid curve), plotted as total
extinction $A(\lambda)$ normalized by $E(B-V)$.  The arbitrarily scaled
profiles of the Johnson and Str\"{o}mgren filters are indicated, and
the Bastiaansen data shown by the small plus signs (``+'').  The
agreement between the $R=3.1$ curve and the Bastiaansen data is clear
from the figure.  The fourth column of Table 4 gives the values of the
photometric extinction ratios produced by the curve, also in very good
agreement with the observations.

For wavelengths longward of 2700 \AA ($1/\lambda < 3.7 \mu$m$^{-1}$),
the $R=3.1$ curve is defined by a cubic spline interpolation between a
set of optical/IR anchor points (filled circles) and a pair of UV
anchor points (filled squares).  The values of the UV anchors (at 2700
\AA\/ and 2600 \AA)  are determined by the FM fitting function for the
case $R=3.1$ (see above) and assure a smooth junction between the
optical and UV regions at 2700 \AA.  The anchor at 0 \micron$^{-1}$ is
fixed at 0 (i.e., no extinction at infinite wavelength) and the values
of the other six optical/IR points were adjusted iteratively to find
the curve shape which best reproduced the extinction observations.  The
wavelengths chosen for the spline anchor points are somewhat arbitrary,
although points in the IR, the optical normalization region, and the
near-UV are clearly required.  The wavelengths and values of the spline
anchors for the $R=3.1$ curve are given in Table 3.

Note the {\em slope} of the derived extinction curve does not approach
zero as $1/\lambda$ approaches zero.  We sacrificed this physically
expected requirement in order to achieve a better fit to the Johnson IR
photometry and to preserve the simplicity of the fitting procedure.  A
zero slope could have been guaranteed by, for example, adopting a power
law to represent the extinction; but no single power law can reproduce
the IR photometry to an acceptable level.  The new curve should be
treated as very approximate at wavelengths beyond the limit of the
M band (i.e., at $\lambda > 6 \micron$).

The overall $R$-dependence of the optical/IR curve is relatively easy
to incorporate by the following adjustments in the spline anchors:  (1)
the UV points are computed by the FM fitting function using the
coefficient values given above, including the $R$-dependent value of
$c2$;  (2) the IR points at $1/\lambda < 1 $\micron$^{-1}$ are simply
scaled by $R$/3.1, since the shape of the far-IR extinction is believed
to be invariant (see \S 2.1 and below); and (3) the optical points are
vertically offset by an amount $R-3.1$, with slight corrections made to
preserve the normalization.  Without these corrections (which are less
than 0.015 over the range $R$ = 2 to 6) the extinction curves would
drift away from the standard normalization, i.e., $E(B-V) = 1$, by a
few hundredths of a magnitude as $R$ departed from the value 3.1.
Table 4 gives formulae for computing the $R$-dependent values of the
optical spline anchors.  The corrections just noted are manifested
in the departure of the linear term from a value of 1.0 and in the
higher order term for the 4110 \AA\/ point.

It is not obvious {\em a priori} that the above adjustments in the
spline should produce the proper curve shapes in the relatively large
gaps between the 4 points in the optical normalization region and the
UV or IR regions. However, it will be shown below that the resultant
curves are in good agreement with observations.

\subsection{The Full {\it R}-Dependent Curve}

Figure 7 shows the full wavelength range of the IR-through-UV
$R$-dependent extinction curves derived here, for four representative
values of $R$ (thick solid and dashed lines).  For $\lambda \le 2700$
\AA\/ the curves are computed using the FM fitting function and for
$\lambda >  2700$ \AA\/ the curves are spline interpolations
between the $R$-dependent spline anchors.   The CCM results are shown
for comparison by the thin dotted lines for the same four $R$ values.
In the UV region, the new results and the CCM curves are similar for $R
\le$ $\sim$4.0, but diverge for larger $R$.  The main reason for this
discrepancy lies in the UV linear extinction component.  While CCM
adopted a linear relationship between the slope $c2$ and intercept
$c1$,  their relation --- shown by the dotted line in the bottom panel
of Figure 5 --- does not agree well with the observations.  The large
discrepancy at small values of $c2$ (i.e., large $R$) produces the
difference seen in Figure 7.

A closeup comparison between the new results and CCM in the optical/IR
region for the $R=3.1$ case can be seen in Figure 6, where the CCM
curve is indicated by the dotted line.  The disagreement in the region
near the Johnson $R$ filter illustrates the bandwidth effect discussed
above.  The CCM curve is fixed at $A(\lambda)/E(B-V) = 2.32$ at 7000
\AA, which is taken as $\lambda_{eff}$ for the Johnson $R$.  The new
curve is constructed to yield the same value of $A(\lambda)/E(B-V)$ but
for synthetic photometry with the Johnson $R$.  The discrepancy is
particularly large for this filter because of the steep change in
extinction across its relatively broad profile.  This effect also
accounts for the lesser discrepancies in the regions of the $U$, $B$,
$I$, and $J$ filters.  O'Donnell (1994) presented a revision to the CCM
formula in the region of the $U$, $B$, and $V$ filters based on the
Str\"{o}mgren photometric indices.  This revision uses the same stategy
as CCM of fixing the extinction ratios at the filter effective
wavelengths, but more closely resembles the newly derived monochromatic
extinction curve because the bandwidth effects in the narrower
Str\"{o}mgren filters are smaller.

Note that the technique used here to produce $R$-dependent curves can
be used to construct ``customized'' UV/optical extinction curves for
sightlines with well-defined UV extinction properties.  It is simply
necessary to substitute the measured values of $A(\lambda)/E(B-V)$ at
2700 \AA\/ and 2600 \AA\/ for the UV spline anchors listed in Table 3
and then perform the cubic spline interpolation to determine the
optical portion of the curve, which smoothly joins the measured UV
curve at 2700 \AA.  This method was used in Figure 2 to extend the
{\it IUE} curves into the optical region.

The assertion that the shape of the IR extinction law is invariant can
be tested by comparing results obtained from the $R$-dependent curves
in Figure 7 --- which were constructed based on this assumption ---
with observations.  This comparison is made in Table 5 for a number of
sightlines which span the observed range in $R$.  The ``model'' data
listed for each sightline are the results from synthetic photometry on
model stellar energy distributions artificially reddened using the
$R$-dependent curves, with $E(B-V)$ values from column 2 of the table
and $R$ values from column 7.  The listed values of $R$ are those that
best reproduce the observed color excesses and the uncertainties result
from assuming uncertainties of $\pm$0.02 mag in $E(B-V)$ and $\pm$0.05
mag in the IR $E(V-\lambda)$.  The agreement between observations and
model values is excellent and there are no strong systematic trends
evident.  The results are thus consistent with an invariant IR
extinction curve.  At wavelengths greater than $\sim$1 \micron, the
extinction curve roughly resembles a power law with an index of $\sim$1.5.
This is similar to that adopted by CCM (1.6), but much flatter than
that of Martin \& Whittet (1990; $\sim$1.8).

Values of $R$ are sometimes estimated from the relation $R \simeq 1.1
\times E(V-K)/E(B-V)$.  Exact relationships between R and the IR color
excesses can be derived for the new $R$-dependent curves.  These are
given by
\begin{equation}
R = 1.39 \times \frac{E(V-J)}{E(B-V)} - 0.02                       \\
\end{equation}
\begin{equation}
R = 1.19 \times \frac{E(V-H)}{E(B-V)} + 0.04                       \\
\end{equation}
\begin{equation}
R = 1.12 \times \frac{E(V-K)}{E(B-V)} + 0.02                       \\
\end{equation}
\begin{equation}
R = 1.07 \times \frac{E(V-L)}{E(B-V)} - 0.01                       \\
\end{equation}
The coefficients in each of these equations actually depend on
the value of E(B-V) itself but, over the range $E(B-V) = 0$ -- 2.0,
vary by only several hundredths.  The values in the equations are the
results for $E(B-V) = 0.5$.

At the blue end of the optical region, the interpolation between the
optical and UV spline anchor points can be tested by comparing
predictions of the extinction indices $E(U-B)/E(B-V)$ and
$E(c1)/E(b-y)$ with photometric measurements.  Figure 8 shows this
comparison, with the difference between the observed values and
predicted values --- derived from synthetic photometry of
``customized'' extinction curves from the FM sample --- plotted against
the slope of the linear component $c2$.  The curves were produced as
described above, by adopting the observed values of $A(\lambda)/E(B-V)$
at 2600 \AA\/ and 2700 \AA\/ as the UV anchor points. The random scatter
(measurement noise) in both indices is large,  but no
systematic trends are seen and the model curves appear to reproduce at
least the integrated properties of extinction well in the near-UV
region. Over the range of $c2$ values shown, the model values of
$E(U-B)/E(B-V)$ range from about 0.6 to 0.8 and the values of
$E(c1)/E(b-y)$ from about 0.0 to 0.3.

In summary, the $R$-dependent curves derived here reproduce existing
photometric and spectrophotometric measurements and provide a good
estimate of  the true monochromatic wavelength dependence of
interstellar extinction in the IR-through-UV regions.  These curves
should be preferred for dereddening UV, optical, and near-IR
spectrophotometry.  Since the new curves give the detailed wavelength
dependence of extinction, they can be used to predict the extinction
relationships for any photometric system by using synthetic photometry
of artificially reddened energy distributions.  An IDL procedure to
produce the new curves at any desired value of $R$ and over any
wavelength range can be obtained from the author or via anonymous ftp
at {\it astro1.vill.edu}.  After logging in, change directories to {\it
pub/fitz/Extinction}, and download the file ``FMRCURVE.pro.''
Alternatively, this directory contains a series of compressed files
named ``FMRCURVEn.n.txt'' (e.g., ``FMRCURVE3.1.txt'') which contain
ASCII versions of the curves for various values of $R$ (``n.n'').  Two
columns of data are contained in each file; the first contains
wavelengths in \AA\/ and the second contains the extinction curve in
$A(\lambda)/E(B-V)$. There are 1099 sets of points in each file.

\end{appendix}

\newpage