---
title: 'Echelle Spectrograph: Design & Applications'
url: https://www.emergentmind.com/topics/echelle-spectrograph
type: topic
---

# Echelle Spectrograph: Design & Applications

An echelle spectrograph is a high-dispersion, cross-dispersed grating spectrograph designed for high-resolution, broad-wavelength astrophysical spectroscopy. It employs an echelle grating operated in high diffraction orders as the principal disperser and a cross-disperser (prism or grating) to spatially separate overlapping echelle orders onto a 2D detector. Echelle spectrographs are foundational for modern exoplanet, stellar, and ISM research, enabling resolving powers $R = \lambda/\Delta\lambda$ in excess of $10^4$–$10^5$ over broad optical and UV bands, and are engineered for stability at the sub-m/s radial-velocity precision level. Their architectures, performance metrics, reduction algorithms, and calibration strategies dictate the limits of ground- and space-based high-fidelity spectroscopy.

## 1. Optical Principles, Architecture, and Design Concepts

The defining element of an echelle spectrograph is the echelle grating—a coarse-ruling (typ. 20–80 lines/mm) grating with high blaze angle (typically $>60^\circ$) used at high diffraction orders ($m\gtrsim 20$), producing high angular dispersion and thus high $\Delta\lambda/\Delta x$. The grating equation is
\[
m\lambda = d (\sin\alpha + \sin\beta)
\]
where $m$ is the diffraction order, $d$ is the groove spacing, and $\alpha$, $\beta$ are the incidence and diffraction angles respectively. In Littrow or quasi-Littrow mounts ($\alpha\approx\beta$), peak efficiency is achieved near the blaze wavelength.

Because high orders overlap in wavelength, a cross-disperser (typically a prism or low-dispersion grating) introduces an orthogonal angular separation. Orders are arranged parallel on a 2D detector, with each order covering a free spectral range (FSR) $\Delta\lambda_\mathrm{FSR} = \lambda/m$.

White-pupil configurations, common in modern designs, refocus the aperture stop (pupil) after each main dispersing element, reducing stray light and vignetting compared to classical double-pass layouts. The entire beam path is engineered for stability, with monolithic optical tables, thermal enclosure, fiber-feeds (often octagonal for modal scrambling), and precisely mounted collimators and cameras.

Illustrative instrument architectures include:
- WES (Weihai Echelle Spectrograph): R = $40,600$–$57,000$, $371$–$1,100$ nm, white-pupil with pair of prisms for cross-dispersion, all-Fiber feed [1611.00116].
- VUES: Modular optomechanics for repeatable alignment, 100 $\mu$m octagonal fiber, selectable slit widths yielding $R=37,000$–$67,000$ over $400$–$880$ nm [1601.06024].
- SES-VIS (STELLA): White-pupil, R4 echelle ($41.6$ lines/mm), R $\simeq 55,000$ at $470$–$690$ nm, double-pass prism cross-disperser, optimized for sub-m/s stability [2012.08218].

Typical spatial/diameter beam sizes range from $50$–$120$ mm, constrained by grating and detector formats as well as desired étendue.

## 2. Performance Metrics: Resolving Power, Coverage, and Stability

Echelle spectrographs are engineered to maximize the following key parameters:
- **Resolving power**: $R = \lambda / \Delta\lambda$. In slit- or fiber-limited regimes, $R$ is set by the projected slit width or fiber image:
  \[
  R \approx mN = m g D
  \]
  where $g$ is groove density (lines/mm), $D$ is beam diameter, $m$ is the order.
- **Wavelength bandpass and order count**: Set by grating size, detector size, cross-disperser, and beam geometry.
- **Sampling**: Typically $\sim2-3$ pixels per FWHM to satisfy Nyquist criteria.
- **Signal-to-noise (SNR)** scaling and throughput: $\mathrm{SNR}(V, t) \propto 10^{-0.2V} \sqrt{t}$. System throughput depends on coatings, fiber losses, detector QE, and atmospheric transmission.
- **Radial-velocity stability and precision**: Influenced by environmental and mechanical stability, with limiting precision determined by:
  \[
  \sigma_v \sim \frac{c}{R \cdot \mathrm{SNR}}
  \]
  Best-case photon-limited precisions approach $1$ m/s (HARPS; R $\sim 115,000$, $S/N>150$ [1312.5471]); practical long-term stabilities depend critically on vacuum/thermal control (PARAS: $R=67,000$, $1.7$–$2.1$ m/s, $T=0.01$–$0.02$C rms at $0.1$ mbar [1312.5471]).

Representative performance comparison (select instruments):

| Instrument    | Resolving Power   | Bandpass (nm) | Typ. RV Stability   |
|---------------|------------------|---------------|--------------------|
| HARPS         | $115,000$        | $380$–$690$   | $<1$ m/s           |
| WES           | $40,600$–$57,000$| $371$–$1,100$ | $<15$ m/s (4 yr)   |
| VUES          | $37,000$–$67,000$| $400$–$880$   | $<10$ m/s (night)  |
| CAFE          | $63,000$–$70,000$| $3,960$–$9,500$| $20$ m/s (long-term)|
| Waltz         | $65,000$          | $450$–$800$    | $\leq5$ m/s (goal) |

## 3. Calibration Frameworks and Wavelength Solutions

Wavelength calibration in echelle spectrographs entails establishing the mapping $\lambda(x, i)$ from pixel $x$ and order index $i$ to physical wavelength. The standard approach is to fit 2D polynomial models or physical-optics models to arc-lamp exposures (typically Th–Ar).

A robust, general model—critical for high-precision work—is:
\[
\lambda(x, i) = (m_0 + i)^{-1} \sum_{q=0}^{N_i} \sum_{c=0}^{N_x} a_{qc}\,P_q(i)\,P_c(x)
\]
where $a_{qc}$ are free coefficients, $P_q(i)$ and $P_c(x)$ are polynomial bases, and $m_0$ is a base order index [1910.08079].

Advanced packages (e.g., xwavecal [1910.08079]) automate this process:
- Exploit duplicated lines in overlapping adjacent orders to constrain the wavelength solution scale up to a multiplicative factor.
- Then solve for the scale and higher-order terms by matching observed line centers to a laboratory reference.
- Employ iteratively reweighted least-squares (IRWLS) to robustly fit the 2D model, excluding outliers at each iteration using MAD-clipping.
- Achieve RMS residuals as low as $5\times10^{-4}$ Å (HARPS), corresponding to $<1$ m/s in velocity units.

Physical forward models, such as paraxial-plus-corrections propagation plus Buchdahl/polynomial aberration corrections [1505.00397], enable sub-pixel ($\leq0.08$ px, HESP) predictions of spectral feature locations, facilitating model-based alignment and drift correction [1906.05064].

## 4. Extraction and Data Reduction Methodologies

Optimal extraction algorithms, essential for maximizing SNR and preserving spectral resolution, have superseded simple boxcar or column summation. The generalized model for observed echelle data is:
\[
D_{ij} = S_j P_{ij} + B_{ij} + N_{ij}
\]
where $D_{ij}$ is the detector value at spatial pixel $i$, dispersion pixel $j$; $S_j$ is the integrated source flux at $j$; $P_{ij}$ is the normalized spatial profile; $B_{ij}$ is the background; $N_{ij}$ is the noise [2008.05827].

Maximum-likelihood solutions for $S_j$ are obtained by solving:
\[
S_j^* = \frac{\sum_i P_{ij} D_{ij} / \sigma_{ij}^2}{\sum_i P_{ij}^2 / \sigma_{ij}^2}
\]
incorporating flat-fielded background subtraction and variance modeling.

Modern algorithms (see PyReduce/REDUCE) address tilted and curved slit images *without prior assumptions* on slit illumination, decomposing each swath into cross-dispersion profiles and propagating oversampled slit traces, with sparse lookup tables to track subpixel mapping [2008.05827]. Bad/cosmic pixels are iteratively masked based on local residuals.

Order splicing and continuum normalization utilize overlap-aware co-addition:
\[
\bar s(x) = \frac{
\frac{x_r - x}{\Delta x}\frac{s_l(x)}{\sigma_l^2(x)} +
\frac{x-x_l}{\Delta x}\frac{\tilde s_r(x)}{\sigma_r^2(x)}
}{
\frac{1}{\sigma_l^2(x)} + \frac{1}{\sigma_r^2(x)}
}
\]
and continuum fitting is regularized via minimization of:
\[
\sum_k w_k[s_k - f_k]^2 + \Lambda_1\sum_k \left(\frac{df}{d\lambda}\right)_k^2 + \Lambda_2\sum_k \left(\frac{d^2 f}{d\lambda^2}\right)_k^2
\]
where $\Lambda_1, \Lambda_2$ tune the smoothness; outliers are iteratively rejected.

## 5. Instrumental Stability, Environmental Control, and Drift Mitigation

Radial-velocity precision at the sub-m/s level is fundamentally limited by environmental drift, mechanical flexure, and optomechanical instability.
- **Temperature and pressure effects**: Refractive index variations of air ($\delta n$) introduce wavelength shifts, $\delta v = c\,\delta n$. Achieved by vacuum enclosures ($<0.1$ mbar) and tight thermal control ($\Delta T < 0.01$ K).
- **Mechanical/thermal design**: Use of low-CTE materials (Invar, Zerodur), vibration isolation, and bench mounting mitigate mechanical drifts and flexure.
- **Instrumental profile (IP) variability**: Fiber scrambling (octagonal fibers, agitation), image slicers, and active field stabilization reduce near- and far-field PSF changes.
- **Simultaneous reference**: Use of a secondary fiber feeding Th–Ar, laser-frequency combs, or Fabry–Pérot etalons to trace drift during observations [1312.5471, 2012.08218].

Bifurcated fiber technology and real-time closed-loop stabilization (e.g., piezo-actuator controlled pickoff mirrors, as in EXOhSPEC [2011.10526]) deliver 4 m/s class stability for compact laboratory spectrographs. Iodine cell reference (stabilized to $\pm0.01$C) enables forward modeling of the instrumental profile and sub-m/s RV calibration [1608.06090].

## 6. Applications, Challenges, and Current Research Directions

Echelle spectrographs enable high-precision exoplanet RV searches, stellar and binary property characterization, ISM studies, and time-domain astrophysics. Their main challenges are:
- **Balancing $R$, throughput, and coverage**: Higher resolving power trades off with lower throughput and higher detector/optical costs.
- **Achieving ultimate environmental and mechanical stability**: mK-level thermal stabilization and high-vacuum operation are expensive and complex.
- **Modeling and control**: Advanced instrument modeling (Zemax, physical raytrace plus polynomials, or machine-learned surrogates) is critical for extracting maximum precision [1505.00397, 1906.05064].
- **Automation in extraction and calibration**: Modern pipelines are expected to provide robust, hands-off optimization (xwavecal, PyReduce/REDUCE), reducing human error and enabling robotic observing networks [1910.08079, 2008.05827].

Pathfinder instruments employing catalog optics with active control loops (EXOhSPEC [2011.10526]), compact educational designs (BACHES [1406.2492]), and FUV space-based instruments (CHESS [1606.04979]) demonstrate the flexibility and reach of the echelle spectrograph concept.

## 7. Comparative Table of Select Echelle Spectrographs

| Instrument         | Resolving Power ($R$) | Wavelength Range (nm) | Throughput      | RV Stability        | Configuration/Notable Features                                 |
|--------------------|----------------------|-----------------------|-----------------|---------------------|---------------------------------------------------------------|
| HARPS              | 115,000              | 380–690               | $\sim$20%       | $<1$ m/s            | Vacuum, double-fiber, R4 grating                              |
| PARAS              | 67,000               | 380–950               | 7–10% (total)   | 1–2 m/s             | Vacuum, dual thermal chamber, fiber-fiber simultaneous ThAr    |
| WES                | 40,600–57,000        | 371–1100              | --              | 10–15 m/s           | White-pupil, dual-prism XF, adjustable slit at fiber exit      |
| CAFE               | 63,000–70,000        | 3,960–9,500           | 20–30% (peak)   | $\sim$20 m/s (long) | Single fiber, white-pupil, cemented stability, fixed central λ |
| BACHES             | 20,000–21,000        | 440–758               | --              | 1–2 km/s            | Compact, slit-fed, remote/automated, 0.5 m telescope           |
| Waltz              | 65,000               | 450–800               | 46.5%           | $<5$ m/s            | White-pupil, image slicer, stabilized iodine cell              |
| VUES               | 37,000–67,000        | 400–880               | 25% (end-to-end)| 10 m/s (night)      | Modular, octagonal fiber, white-pupil                          |
| SES-VIS            | 55,000               | 468–691               | --              | $1$ m/s (goal)      | Vacuum, Fabry-Perot, thermal control to 0.01K                  |

## References

- “Optimal extraction of echelle spectra: getting the most from observations” [2008.05827]
- “WES - Weihai Echelle Spectrograph” [1611.00116]
- “Automatic Echelle Spectrograph Wavelength Calibration” [1910.08079]
- “A small actively-controlled high-resolution spectrograph based on off-the-shelf components” [2011.10526]
- “Spectroscopic Survey of Eclipsing Binaries with a Low Cost Échelle Spectrograph -- Scientific Commissioning” [1602.01238]
- “Stability analysis of VBT Echelle spectrograph for precise radial velocity measurements” [1706.09202]
- “The re-flight of the Colorado high-resolution Echelle stellar spectrograph (CHESS): improvements, calibrations, and post-flight results” [1606.04979]
- “Modelling high resolution Echelle spectrographs for calibrations: Hanle Echelle spectrograph, a case study” [1505.00397]
- “First Light results from PARAS: The PRL Echelle Spectrograph” [1007.4280]
- “CAFE: Calar Alto Fiber-fed Echelle spectrograph” [1301.2066]
- “Design and Construction of VUES: the Vilnius University Echelle Spectrograph” [1601.06024]
- “A High Resolution Spectrograph for the 72 cm Waltz Telescope at Landessternwarte, Heidelberg” [1608.06090]
- “The PRL Stabilized High Resolution Echelle Fiber-fed Spectrograph: Instrument Description & First Radial Velocity Results” [1312.5471]
- “Using raytracing to derive the expected performance of STELLA's SES-VIS spectrograph” [2012.08218]
- “BACHES - a compact échelle spectrograph for radial velocity surveys with small telescopes” [1406.2492]
- “Second generation spectroscopic instrumentation for the STELLA robotic observatory” [2012.08210]
- “Estimation and correction of the instrumental perturbations of Vainu Bappu Telescope Echelle spectrograph using a model-based approach” [1906.05064]

Source: https://www.emergentmind.com/topics/echelle-spectrograph