--- title: Differential Chromatic Refraction url: https://www.emergentmind.com/topics/differential-chromatic-refraction-dcr type: topic --- # Differential Chromatic Refraction Searching arXiv for recent and foundational papers on Differential Chromatic Refraction. Differential Chromatic Refraction (DCR) is the wavelength-dependent displacement of an astronomical source’s apparent position caused by atmospheric refraction. Because the refractive index of air decreases with wavelength, shorter-wavelength light is bent more strongly than longer-wavelength light, so the apparent position, centroid, and effective point-spread function (PSF) of a source depend on its spectral energy distribution (SED), the bandpass, and the observing geometry. In broadband imaging, DCR appears as a color-dependent centroid shift toward the zenith and as a wavelength-dependent elongation of the PSF along the parallactic direction; in spectroscopy and high-contrast imaging the same underlying phenomenon is often discussed as atmospheric dispersion or differential atmospheric refraction (DAR) [2304.01858]. ## 1. Physical basis and terminology DCR arises because Earth’s atmosphere is a dispersive medium: its refractive index $n(\lambda)$ varies with wavelength, so light at shorter wavelengths is bent more strongly than light at longer wavelengths. The first-order description used across the literature is $$ R(\lambda, z) \approx [n(\lambda)-1]\tan z, $$ with differential refraction between two wavelengths given by $$ \Delta R(\lambda_1,\lambda_2,z)\approx [n(\lambda_1)-n(\lambda_2)]\tan z. $$ This dependence on $\tan z$ makes the effect grow with zenith distance, or equivalently with airmass, and the chromatic term is largest in bluer passbands [2304.01858]. The observed displacement is directed along the zenith direction, usually parameterized on the sky by the parallactic angle. In imaging work this means that DCR is not an arbitrary astrometric perturbation: it is highly anisotropic, aligned with the altitude or parallactic direction, and its sign and amplitude are controlled by observing geometry and source color. For objects observed near the meridian, that projection can become predominantly declinational, which is why DCR in Pan-STARRS1 mean positions was treated primarily as a declination bias [2206.09988]. The same physics is described under slightly different names in different subfields. Weak-lensing and time-domain imaging papers generally use “Differential Chromatic Refraction,” while high-contrast imaging and instrumentation papers often use “Differential Atmospheric Refraction” or “atmospheric dispersion.” The underlying phenomenon is the same wavelength-dependent atmospheric refraction, but the emphasis changes with application: centroid shifts and PSF anisotropy in survey imaging, residual dispersion after an atmospheric dispersion corrector in coronagraphy, and wavelength-dependent flux losses in multi-object spectroscopy [2404.06378]. ## 2. Mathematical description, broadband centroids, and geometry For broadband observations, DCR is not set by a single monochromatic wavelength but by an SED-weighted refraction across the full passband. A physically transparent description integrates the source SED and total throughput over wavelength. In the DES supernova analysis, the color-dependent centroid shift relative to a reference spectrum is written as $$ \delta \approx \frac{\int S(\lambda)\,T(\lambda)\,R(\lambda,z)\,d\lambda}{\int S(\lambda)\,T(\lambda)\,d\lambda} - \mathrm{reference}, $$ where $S(\lambda)$ is the source SED and $T(\lambda)$ is the system throughput [2304.01858]. Closely related formulations appear in weak-lensing analyses through the photon-weighted first moment $\bar R$ and variance $V$ of the DCR kernel [1409.6273]. This formalism makes explicit why DCR is source-dependent. Two objects observed in the same exposure, at the same airmass and through the same filter, can have different effective refractions because their SEDs weight the bandpass differently. The distinction is especially important for objects with time-variable spectra. Type Ia supernovae evolve in color and SED with phase, so the DCR-induced centroid and PSF shape evolve over the light curve; stellar flares become bluer during impulsive heating, shifting their effective wavelength and causing measurable zenith-ward centroid motion in a single band [2304.01858]. The sky projection is determined by the parallactic angle. In the DES-SN5YR appendix, the parallactic-angle relations are $$ \sin(q)=\frac{-\cos(\phi)\sin(A)}{\cos(\delta)}, \qquad \cos(q)=\frac{\sin(\phi)-\sin(\delta)\sin(a)}{\cos(\delta)\cos(a)}, $$ with $\phi$ the observer latitude, $A$ azimuth, $\delta$ declination, and $a$ altitude [2304.01858]. In astrometric and flare applications, the displacement is commonly projected into equatorial coordinates through $$ \Delta\alpha\cos\delta=-\Delta R\sin q,\qquad \Delta\delta=\Delta R\cos q, $$ or into a scalar zenith-directed component such as $d_{\parallel}$, which isolates the atmospheric signal along the vertical circle [2008.03446]. Several papers also define an effective wavelength to compress the band-integrated SED dependence into a single chromatic coordinate. In survey photometry this can be written as $$ \lambda_{\mathrm{eff}}=\frac{\int \lambda\,S(\lambda)\,F(\lambda)\,d\lambda}{\int S(\lambda)\,F(\lambda)\,d\lambda}, $$ while flare studies adopt a log-weighted form for $\lambda_{\mathrm{eff}}$ inside a filter [2304.01858]. These constructions are not identical across papers, but they serve the same purpose: mapping a broadband SED to an effective chromatic refraction. ## 3. Imaging manifestations: astrometry, PSF morphology, and photometric bias In imaging data, DCR manifests as a centroid shift of the source along the direction to the zenith and as a wavelength-dependent elongation of the PSF. In weak-lensing language, single-exposure DCR adds variance only along altitude, so the DCR-induced second-moment error is an altitude-aligned excess $\Delta V$ in $I_{yy}$; multi-epoch centroid differences then produce additional second-moment errors in stacked images through color-dependent misregistration [1409.6273]. The principal consequence is that a stellar PSF is not the correct effective PSF for a galaxy or transient with a different SED. This mismatch propagates directly into measurement bias. In DES Scene Modeling Photometry (SMP), the supernova flux estimator is written schematically as $$ F_{\mathrm{SN}}=\frac{\sum_i F_iP_i}{\sum_i P_i^2}, $$ so even small centroid errors or width mismatches change the PSF weights $P_i$ and bias the inferred flux [2304.01858]. Simple tests with Gaussian and Moffat PSFs show that a centroid error of $10\%$ of the PSF FWHM induces a $12$–$15$ mmag magnitude bias, and a PSF-width mismatch of $3\%$ induces $\approx 30$ mmag biases [2304.01858]. The DES-SN5YR analysis separated these effects into astrometric “COORD” corrections and PSF-shape “SHAPE” corrections. Ignoring the astrometric offsets and PSF shape changes causes average biases of $+0.2$ mmag and $-0.3$ mmag respectively, with standard deviations of $0.7$ mmag and $2.7$ mmag across all DES observing bands (\textit{griz}) throughout all redshifts used in the analysis [2304.01858]. The largest corrections appear in the bluest bands: for COORD, $1.3\%$ of $g$-band corrections exceed $0.005$ mag, while for SHAPE about $7.1\%$ of $g$ and $r$ corrections exceed $0.005$ mag [2304.01858]. Weak-lensing analyses reach the same conclusion at the level of shear systematics. Uncorrected DCR and wavelength-dependent seeing induce PSF moment errors that map into multiplicative and additive shear biases, and the paper “Impact of Atmospheric Chromatic Effects on Weak Lensing Measurements” shows that uncorrected atmospheric chromatic effects can bias shear beyond DES and LSST statistical limits [1409.6273]. A central result of that work is that catalog-level corrections are insufficient because a “chromatic model bias” appears when a galaxy convolved with its galactic effective PSF is fitted with a stellar effective PSF. ## 4. Empirical and forward-model correction strategies A major methodological division in the literature is between empirical astrometric calibration and explicit chromatic forward modeling. In Pan-STARRS1, DCR in mean positions was corrected empirically by iteratively subtracting color- and declination-dependent PS1/Gaia EDR3 declination residuals. For about $440$ million point-like reference objects cross-matched to Gaia EDR3, Gaia EDR3 provided a $\sim 3\%$ improvement over Gaia DR2, and DCR corrections provided an additional $\sim 5\%$ improvement; for blue objects with $01.2$ given LSST’s minimum specified requirement on single-visit relative astrometric accuracy [2402.06002]. That proposal was realized observationally in “First Temperature Profile of a Stellar Flare using Differential Chromatic Refraction,” which reports the first derivation of a stellar flare temperature profile from single-band photometry [2507.19584]. For a giant DECam $g$-band flare on an M7 dwarf observed at $1.45\lesssim X\lesssim1.75$, the paper measured zenith-ward centroid drift, modeled the flare SED as a blackbody plus line contributions, and—under its most physically motivated scenario—found a maximum temperature of $19{,}510\,\mathrm{K}$ and a duration of about $2.4$ minutes above $14{,}000\,\mathrm{K}$ [2507.19584]. The same study emphasizes that ignoring line emission can overestimate temperature. Supernova cosmology has also begun to use DCR constructively. “Astrometric Redshifts of Supernovae” models the DCR time series of LSST Deep Drilling Field SNe Ia and uses the resulting astrometric offsets to constrain redshift [2405.04522]. For a conservative $5$-mas systematic uncertainty floor, the paper finds that the astrometric redshift estimation is accurate at $z<0.6$, and that combining DCR-based astrometric redshifts with host-galaxy photometric redshifts and supernova light-curve redshifts considerably improves the overall redshift estimates [2405.04522]. This suggests that DCR is becoming a measurable transient observable, not only a calibration term. ## 6. Instrumentation, residual dispersion, and future survey requirements Instrumental design determines whether DCR is exposed, tolerated, or actively corrected. Rubin/LSST and DECam do not use atmospheric dispersion correctors, so DCR remains present in survey imaging and must be modeled in astrometry, PSF construction, and scene-modeling photometry [2402.06002]. By contrast, high-contrast imagers and fibre spectrographs often attempt active correction, but the residuals remain scientifically important. The Gemini Planet Imager study provides a direct example of incomplete correction. Its atmospheric dispersion corrector was designed to reduce residual DAR below $5$ mas, but measured averages are $\sim 7$ mas in $H$ band and $\sim 11$ mas in $J$ band [2404.06378]. The analysis attributes the $H$-band residual largely to a nearly constant undercorrection of about $7$ mas and the $J$-band residual largely to perpendicular dispersion induced by the ADC, which accounts for $86\%$ of the residual DAR in that band [2404.06378]. The same paper also shows that the Roe-like refractivity model used by the instrument underestimates humidity’s impact on DAR and introduces a new humidity-inclusive approximation for $n(\lambda)$ from $0.7$ to $1.36\,\mu$m [2404.06378]. In fibre spectroscopy, the dominant observable is wavelength-dependent flux loss rather than centroid bias alone. “Atmosphyre: Modelling Atmospheric Chromatic Dispersion for Multi-Object Spectrographs” models the chromatic displacement $$ \Delta R = 206265\,[n(\lambda)-n(\lambda_{\rm ref})]\tan z $$ and its impact on fibre coupling [2407.14356]. The paper concludes that the guiding wavelength should typically be bluer than the observing-band mid-wavelength, around $25$–$45\%$ of the way through the band, and finds for MOSAIC that differential losses greater than $10\%$ are unavoidable for $1$ h observations after a local hour angle of $2.5$ h or at declinations below $-60^\circ$ and above $10^\circ$ [2407.14356]. It therefore argues that an atmospheric dispersion corrector would significantly reduce spectral distortions and increase survey speed. Future survey guidance is consistent across these domains. The DES-SN5YR paper recommends deriving per-exposure, and even per-object, color estimates; incorporating atmospheric dispersion explicitly in PSF models; using improved atmospheric refractivity models with contemporaneous pressure, temperature, and humidity; and, where possible, using full SED information rather than a single color [2304.01858]. Weak-lensing work similarly emphasizes exposure-specific $z$ and $q$, PSF-level chromatic corrections, and accurate throughput knowledge [1409.6273]. A plausible implication is that DCR will increasingly be treated as a first-class component of forward image formation in high-precision survey pipelines, especially for $u$-band observations and high-airmass data, where its amplitude grows rapidly [2304.01858]. ## 7. Limitations, misconceptions, and scope A common misconception is that DCR is only an astrometric shift. The literature consistently shows that this is incomplete. In broadband data, DCR also alters second moments, induces PSF anisotropy aligned with the parallactic direction, and can create photometric bias when stellar PSFs are applied to objects with different SEDs [1409.6273]. In time-domain photometry, the resulting bias can remain at the millimagnitude level in DES yet still become important for future higher-precision surveys [2304.01858]. Another misconception is that a single broadband color fully characterizes DCR. Several papers use color proxies—such as $g-i$ for supernovae or BP-RP for astrometric calibration—because they are practical and often effective, but the limitations are explicit. The DES supernova analysis notes that using $g-i$ as a proxy for the full supernova SED ignores phase-dependent spectral features within the bandpass; host-galaxy light is assumed static and non-dispersed; interpolation is performed on precomputed grids rather than exact exposure conditions; and focal-plane optical variations such as lateral color are averaged over rather than modeled explicitly [2304.01858]. These are approximations, not exact equivalences. A further misconception is that DCR can always be removed at the catalog level after image measurement. Weak-lensing analyses show that this is false in the presence of chromatic model bias: if the wrong PSF is used during model fitting, later catalog corrections cannot reconstruct the measurement that would have been obtained with the correct chromatic PSF [1409.6273]. This is why PSF-level or scene-level chromatic modeling is repeatedly favored in high-precision applications. The current literature therefore places DCR in a dual role. It is simultaneously a calibration complication and a measurable atmospheric encoding of source SED. In astrometry and photometry it must often be corrected to avoid bias; in flare physics, subband reconstruction, and supernova redshift inference it can be exploited as a signal [2402.06002]. This duality is central to the modern treatment of DCR: its scientific meaning depends less on the underlying atmospheric physics—which is fixed—than on whether the observing program regards chromatic refraction as contamination, calibration leverage, or an observable in its own right. Source: https://www.emergentmind.com/topics/differential-chromatic-refraction-dcr