abstract: A new method for measuring the large-scale structure of the Galactic magnetic field is presented. The Galactic magnetic field has been probed through the Galactic disk with near-infrared starlight polarimetry, however the distance to each background star is unknown. Using red clump stars as near-infrared standard candles, this work presents the first attempt to decompose the line of sight structure of the sky-projected Galactic magnetic field. Two example lines-of-sight are decomposed: toward a field with many red clump stars and toward a field with few red clump stars. A continuous estimate of magnetic field orientation over several kiloparsecs of distance is possible in the field with many red clump stars, while only discrete estimates are possible in the sparse example. toward the Outer Galaxy, there is a continuous field orientation with distance that shows evidence of perturbation by the Galactic warp. toward the Inner Galaxy, evidence for a large-scale change in the magnetic field geometry is consistent with models of magnetic field reversals, independently derived from Faraday rotation studies. A photo-polarimetric method for identifying candidate intrinsically polarized stars is also presented. The future application of this method to large regions of the sky will begin the process of mapping the Galactic magnetic field in a way never before possible.
author: Michael D. Pavel
title: Using Red Clump Stars to Decompose the Galactic Magnetic Field with Distance

# Introduction

Observationally, our knowledge of the large-scale structure of the Galactic magnetic field has primarily been derived from studies of background starlight polarimetry , Faraday rotation , and synchrotron emission . Yet only Faraday rotation has successfully resolved large-scale magnetic field structures with distance. This work presents a new method for resolving the large-scale magnetic field geometry with distance by using red clump stars as standard candles in the Galaxy. Combined with near-infrared (NIR) starlight polarimetry, the magnetic field properties of a line of sight can be decomposed with distance to map variations in the plane-of-the-sky magnetic field.

Previous Faraday rotation studies have applied similar methods to studies of the large-scale structure of the Galactic magnetic field. By combining observations of Galactic pulsars with distance estimates (that sample discrete portions of the Galactic magnetic field) with observations of polarized extragalactic sources (that sample the entire Galactic line of sight towardh that object) detailed models of the Galactic magnetic field have been generated . However, these studies suffer from the limitations of Faraday rotation, specifically that they only sense the line of sight component of magnetic fields. Also, pulsar distances can be uncertain. Few parallax measurements are available and only for nearby pulsars, HI kinematic distances typically only give upper or lower distance limits, and those distances have typical uncertainties of 0.5 kpc or greater .

Starlight polarization arises as an integral quantity of the line of sight between the background star and the observer. Regions of magnetically-aligned dust grains impose weak polarizations parallel to the local magnetic field direction on background starlight. Via this mechanism the orientation of magnetic fields are directly measurable , the plane-of-sky magnetic field strength can be indirectly inferred , though the field direction is not obtainable. Statistical studies of observed starlight polarization have shed light on the general symmetry and spiral-type pitch angle of the Galactic magnetic field , but even kpc-scale details of the morphology of the magnetic field remain hidden. To measure changes in the sky-projected magnetic field with distance, as is done with Faraday rotation for the line of sight field, reliable stellar distance markers are needed.

The first steps in using starlight polarimetry to disentangle magnetic fields with distance have already been taken. attempted to measure magnetic fields in the Galactic Center with NIR polarimetry by subtracting the polarization associated with foreground "bluer" field stars from the polarization of "redder" stars located near the Galactic Center and showed that both populations were polarimetrically indistinguishable. This work was followed by who used the same method to observe a transition from toroidal to poloidal magnetic field above and below the Galactic center. use a similar foreground subtraction method to remove a foreground Galactic contribution and resolve the magnetic field of a galaxy seen through the Taurus molecular cloud. This work extends magnetic field decomposition by moving away from statistical treatments of foreground stars, and instead uses red clump stars as standard candles to identify and characterize changes in magnetic field orientation with distance.

In color-magnitude diagrams (CMDs) of star clusters, the red clump is an overdensity of core helium burning (horizontal branch) stars seen in intermediate-age ($\sim 1-10$ Gyr) clusters and old ($\geq 10$ Gyr), metal-rich clusters . The relative invariance of red clump star K-band ($2.2\mu$m) absolute magnitude with age or metallicity allows it to serve as a K-band standard candle . In a NIR CMD of field stars, red clump stars located at different distances form an easily distinguishable feature. The photometric properties and methods for identifying this feature have been studied in the Galactic disk , the Large Magellanic Cloud , and in star clusters . These studies have further refined our knowledge of the intrinsic properties of these objects, and limitations to their use.

This work implements the NIR red clump selection algorithm described in to identify red clump stars at different distances. For a population of stars with a fixed absolute magnitude, the apparent magnitude is an excellent proxy for distance order, since extinction monotonically increases with distance along a given line of sight. Naively assuming a constant dust density along a sightline, the physical distances to these stars are estimated below. These stars are then used as distance markers along a line of sight to probe magnetic field properties between them via NIR starlight polarimetry.

The details of a new method for measuring the Galactic magnetic field structure with distance is presented in Section 2. The method is applied to example lines-of-sight toward the inner Galaxy in Section 3 and the outer Galaxy in Section 4. Section 5 presents new results on the morphology of the Galactic magnetic field toward these two lines-of-sight and also discusses evidence for intrinsically polarized red clump stars. A summary and conclusions are presented in Section 7.

# Decomposing Sky-Projected Magnetic Field Orientation Along the line of sight

Starlight linear polarization arises from dichroic extinction via aligned dust grains along the line of sight between a star and an observer. Spinning, non-spherical dust grains preferentially align with their long axes perpendicular to the local magnetic field orientation . These aligned dust grains preferentially extinct photons polarized perpendicular to the local magnetic field direction (parallel to a grain's long axis), giving rise to a weak linear polarization along the magnetic field direction. describes the equations of polarized radiative transfer and discusses the effects of magnetic field geometry on the observed polarization properties.

Summarizing , linear polarization will trace the average, extinction-weighted, sky projected orientation of the magnetic field between an unpolarized star and observer. Magnetic field components parallel to the line of sight do not contribute to this polarization, though these are the fields traced by Faraday rotation. The Galactic magnetic field consists of regular and turbulent components in approximate equipartition. Along long (kpc) lines-of-sight in the diffuse interstellar medium, the turbulent magnetic field contributions will tend to cancel, allowing linear polarimetry to trace the large-scale structure of the Galactic magnetic field.

If the distances to stars with measured polarizations can be estimated, then the average sky-projected magnetic field orientation in the material between those stars can be measured via differences in their Stokes parameters. These differences arise in the dusty, magnetized medium between the stars. This method of decomposition requires that the different stars sample the same line of sight so that polarizing foregrounds can be reliably removed. If the angular separation on the sky between stars is large, then their light may not be sampling the same polarizing layers. The typical size scales for significant variation in the random component of the Galactic magnetic field are $10-100$ pc . Two stars 5 kpc distant and separated by 10 arcmin on the sky would only be physically separated by 14.5 pc, which is smaller than the typical random fluctuations in the Galactic magnetic field, and much smaller than the size scales for variations in the regular magnetic field. As the light from these two stars approach an observer and converge, they propagate through almost identical paths. If instead one star is located between the observer and the other star, the conditions are favorable for measurement of the magnetic properties between the stars via Stokes subtraction.

## Identifying Red Clump Stars

Red clump stars are identified using the CMD method described in , illustrated in Figure with Two Micron All Sky Survey (2MASS) data for a $100' \times 100'$ region centered on $(\ell , b)= (19.49\degr,+0.56\degr)$. In Figure , two significant features are seen. Main sequence stars of various reddenings are located in the region $0.5 < (J-K) < 1.0$ and $K > 8 + 4(J-K)$ (approximately traced by the solid white triangle in Figure ), and the roughly diagonal feature running from the top-center toward the lower-right is dominated by red clump stars.

To objectively quantify the stars in the red clump CMD sequence, first grouped all stars by $K$ in 0.5 mag wide groups for $8<K<13$. For each $K$ group, a histogram of $(J-K)$ was calculated with 0.05 mag wide $(J-K)$ bins. This histogram was fit with a Gaussian to identify the $(J-K)$ location of the peak of the red clump feature at that $K$. The locus of points tracing the red clump stars was fit with a polynomial. Stars within 0.2 mag of the $(J-K)$ center of the feature are considered red clump stars. In Figure , a dashed line outlines this region.

The details of this fitting are shown in Figure . The well-separated red clump sequence in each $K$-band brightness bin (labeled in each panel) is shown. A similar histogram for the bluer main sequence is not shown in each brightness bin. The fitted Gaussian to each histogram is also shown. These centers of these fits are marked with vertical solid lines. Stars within 0.2 mag of each $(J-K)$ center are classified as red clump stars, and these color limits are shown by the vertical dashed lines in each panel. For fainter stars, the red end of this feature may suffer incompleteness (e.g., the lower right panel in Figure ), which shifts the fitted center up to 0.15 mag bluer in $(J-K)$ than would be identified by eye. However, this remains a conservative method for identifying red clump stars. Some actual red clump stars are excluded so as to minimize false positive identifications. Sample contamination will be discussed in more detail in the next subsection.

The polarimetric observations described below cover single $10' \times 10'$ regions of sky, which typically contain too few stars with 2MASS photometry to reliably identify the red clump feature. To assist with identification of red clump stars within this small region, a much larger sample of 2MASS photometry is needed. As described above, a $100' \times 100'$ region was used for the red clump star identification. Then, only those red clump stars within the much smaller $10' \times 10'$ region were flagged for use. As an example, the four red clump stars identified in the $10'\times 10'$ region toward $(\ell , b)= (19.49\degr,+0.56\degr)$ are marked with diamonds in Figure

## Sample Contamination

A key concern about the red clump identification method described above is the contamination of other stars into this region of the CMD. In particular, reddened main sequence stars could make up a significant fraction of the selected stars. comment on dwarf contamination and quantify the fraction of dwarfs in the total star counts as a function of magnitude (see Figure 5 in that paper). However, that analysis is for all stars in the field and does not quantify the effectiveness of the NIR color selection method at removing reddened interlopers.

To estimate the total contamination along a line of sight, we use the stellar population synthesis model of to generate sample star fields toward an example direction in the inner Galaxy, $(\ell , b)= (20\degr,0\degr)$. In addition to observed photometric properties, the intrinsic properties of each star are known. The average extinction per unit length was varied from $A_V /L = 0.7-1.8$ mag kpc$^{-1}$ in the different Monte Carlo samples. Using $A_V /L=$`<!-- -->`{=html}1.8 mag kpc$^{-1}$ is a reasonable assumption for the Galactic disk, but we also included lower values as a weak proxy for sightlines at higher Galactic latitudes.

The first quantity calculated was the dwarf fraction of all stars as a function of $K$-band apparent magnitude, the same quantity calculated in for the outer Galaxy. This quantity is plotted by the dashed line in Figure . As shown in that figure, the dwarf contamination fraction rises for fainter stars. This contamination should be dominated by nearby, cool main sequence stars . The contamination fractions calculated here are lower than the Outer Galaxy fractions from . This difference can be attributed to two effects: (1) lower extinction in the outer Galaxy allows intrinsically fainter stars to remain above a given brightness limit (reducing the Malmquist bias over a fixed volume), and (2) the distant (several kiloparsec) background in the outer Galaxy has a lower stellar space density of very bright sources (e.g., red clump stars) compared to the inner Galaxy.

Next, the contamination of the fitted red clump sequence was calculated. The photometric method described in Section was applied to each of the simulated fields to identify red clump candidates. In this case dwarfs and supergiants, though rare, were considered as interlopers. The average fraction of interlopers is shown by the solid line in Figure and the shaded region represents the uncertainty of that fraction. This uncertainty contains two components: the intrinsic uncertainty in the contamination fraction and the variation caused by different assumptions about the average extinction per unit length. In general, lower extinction per unit length is associated with more red clump contamination (toward the upper envelope of the shaded region), and larger extinction per unit length is associated with less contamination (toward the lower envelope of the shaded region). Later, the measured contamination fraction in Figure will be used to estimate the probability of false detections.

## Distance Estimation

To estimate the distance of a red clump star, the K-band absolute magnitude was assumed to be $M_K=-1.57\pm0.05$ . This is slightly dimmer than the value used by . K-band extinction was estimated from the apparent NIR (H-K) colors, assuming (H-K)$_o = 0.09\pm0.03$ , $R_V = 3.1$, and the extinction law.

## Isolating Polarizing Layers in the ISM

Together, the position angle (PA) and degree of polarization ($P$) fully characterize the linear polarization properties of observed starlight, and correspond to a unique pair of linear Stokes parameters, $U$ and $Q$. Consider a series of $n$ optically thin polarizing layers characterized by Stokes parameters ($U_1$, $Q_1$), ($U_2$, $Q_2$),...,($U_n$, $Q_n$). The linear polarization ($U_f$, $Q_f$) of an initially unpolarized beam of light propagated through these layers becomes: $$U_f = \displaystyle\sum\limits_{i=0}^n U_i$$ $$Q_f = \displaystyle\sum\limits_{i=0}^n Q_i,$$ under the assumption that all layers are optically thin. With measurements of ($U_f$, $Q_f$) through different numbers of polarizing layers, the properties of those polarizing layers can be measured by differencing Stokes parameters for stars at different distances. Similar techniques for the isolation of the polarization properties of stars in the Galactic center and a background galaxy seen through Taurus have been used before. However, those cases only considered the measurement and removal of a single, foreground polarization component.

In the following, this method is applied to the entire line of sight to decompose the sky-projected magnetic field orientation with distance. The first example below (Section ) demonstrates a coarse implementation of this procedure for an example line of sight with few identified red clump stars. In Section , a more sophisticated implementation is shown later for a line of sight with many red clump identifications.

## The Effect of Discrete Clouds

An individual (dense) molecular cloud along the line of sight probed by this method will have a significant effect on the observed starlight polarization. Since starlight polarimetry uses the *dust-weighted* magnetic field orientation, the presence of these higher density regions could contaminate the polarization signal expected from the large-scale Galactic magnetic field.

Comparative studies of sub-millimeter thermal dust polarization and optical starlight polarimetry have shown that the magnetic orientation of these clouds are generally aligned with the local magnetic field . This has also been observed in giant molecular clouds in M33 by . The polarization signal superimposed by these clouds also depends on the degree of turbulence within the cloud: a magnetically turbulent cloud will produce a lower polarization signal per unit extinction than a less turbulent cloud .

However, these high-density regions can be identified via star counts , long wavelength continuum emission , or line emission from high-density tracers . Velocity-resolved line emission at radio wavelengths is perhaps the best method for identifying individual high-density regions along a single line of sight. Once identified, the polarimetric effect of these clouds can be individually assessed. In the best case scenario, red clump stars with high polarimetric signal-to-noise on both sides of the cloud are found and the polarization signal from the cloud can be independently measured. In the more likely scenario the polarization signal from the cloud will dominate any signals from more distant stars making analysis beyond the cloud impossible. In this case the clouds will create "shadow zones" where the magnetic field cannot be measured by this method.

# Discrete Magnetic Field Decomposition with Few Red Clump Stars

$H$-band ($1.6\mu$m) starlight polarimetry toward $(\ell , b)= (19.49\degr,+0.56\degr)$ was drawn from the first data release of the Galactic Plane Infrared Polarization Survey . These observations were made with the Mimir instrument on the 1.8m Perkins telescope of Lowell Observatory outside Flagstaff, AZ. Each observation covered a $10'\times 10'$ field-of-view and data were reduced with custom IDL-based software . For each observation, a catalog of sky coordinates, starlight polarization properties, and measured $H$-band photometry were produced. These catalogs were matched to 2MASS and GLIMPSE point source catalogs. GPIPS polarimetric observations are complete to $H\approx 12$, which provides a faintness (and ultimately distance) limit on red clump stars with measurable polarimetry. Data from the $^{13}$CO Galactic Ring Survey (Jackson et al. 2006) shows no significant emission at any velocity toward this direction indicating that it is free from cloud contamination.

Using the photometric method described in Section , four candidate red clump stars were identified in this $10'\times 10'$ field. Using the estimated contamination fractions presented in Figure , the K-band magnitude maps directly to a false-positive probability for that star. These $K$-band magnitudes are 8.82, 10.57, 10.92, and 11.67; which become false-positive probabilities of 0.141, 0.049, 0.047, and 0.045. The probability that all four stars are actually red clump stars (a reliability probability) is only 0.743 with the brightest star being the least likely. Nevertheless, the following analysis assumes that all four stars are actually red clump stars.

```latex
\begin{deluxetable*}{ccccccc}
\tablewidth{0pt}
\tablecolumns{7}
\tablecaption{Inner Galaxy $H$-band Polarimetry\label{pol_table}}
\tablehead{
    \colhead{Star} &
    \colhead{$H_{2MASS}$} &
    \colhead{Distance} &
    \colhead{P\tablenotemark{a}} &
    \colhead{GPA\tablenotemark{b}} &
    \colhead{$U$} &
    \colhead{$Q$} \\
    \colhead{} &
    \colhead{(mag)} &
    \colhead{(kpc)} &
    \colhead{(\%)} &
    \colhead{(deg)} &
    \colhead{(\%)} &
    \colhead{(\%)} \\
    \colhead{(1)} &
    \colhead{(2)} &
    \colhead{(3)} &
    \colhead{(4)} &
    \colhead{(5)} &
    \colhead{(6)} &
    \colhead{(7)} \\
    }
\startdata
1 & $9.04\pm0.02$ & $1.07\pm0.23$ & $0.31\pm0.15$ & $83.8\pm12.5$ & $0.23\pm0.15$ & $0.25\pm0.15$ \\
2 & $10.89\pm0.04$ & $2.19\pm0.61$ & $2.87\pm0.69$ & $91.4\pm  6.9$ & $2.52\pm0.68$ & $1.54\pm0.72$ \\
3 & $11.31\pm0.02$ & $2.45\pm0.56$ & $2.08\pm0.68$ & $75.3\pm  9.2$ & $0.97\pm0.71$ & $1.96\pm0.68$ \\
4 & $12.09\pm0.03$ & $3.41\pm0.87$ & $4.67\pm1.94$ & $59.0\pm11.0$ & $-0.55\pm1.94$ & $5.03\pm1.94$ \\
\enddata
\tablenotetext{a}{These P values and the derived position angle uncertainties have been corrected for positive statistical bias via the equation $P=\sqrt{U^2+Q^2-\sigma_P^2}$ \citep{1974ApJ...194..249W,2006PASP..118.1340V}.}
\tablenotetext{b}{GPA values were derived from the equatorial Stokes parameters, then rotated to the Galactic coordinate system using the equation in \citet{1968ApJ...151..907A}.}
\end{deluxetable*}
```

For the four red clump stars identified toward $(\ell , b)= (19.49\degr,+0.56\degr)$, the Galactic position angle (GPA) is plotted against the 2MASS $K$-band magnitude in the top panel of Figure and all of the H-band polarization properties are listed in Table . GPA is given in the Galactic coordinate system, where GPA=$90\degr$ is parallel to the Galactic disk, shown by the dashed line in Figure . In the bottom panel of that Figure, the $x$-axis shows the estimated distance to each red clump star.

Qualitative estimates of the magnetic properties with distance can be deduced from Figure . The first two stars have GPAs consistent with $90\degr$, suggesting that the magnetic field in this direction is parallel to the Galactic plane out to a distance of $\sim$`<!-- -->`{=html}2 kpc. The next two stars show a monotonic change in GPA with apparent magnitude (i.e., distance). Beyond 2 kpc, the magnetic field orientation changes. An alternative interpretation is that there is a constant change in GPA with $K$ or distance (representable by a straight line with non-zero slope) instead of a knee. A more quantitative analysis can be made using the procedure described below.

```latex
\begin{deluxetable}{ccccc}
\tablewidth{0pt}
\tablecolumns{5}
\tablecaption{Inner Galaxy Magnetic Field Properties, Binned by Distance\label{mag_table}}
\tablehead{
    \colhead{Interval\tablenotemark{a}} &
    \colhead{Distance} &
    \colhead{GPA$_i$} &
    \colhead{$\Delta U_i$} &
    \colhead{$\Delta Q_i$} \\
    \colhead{} &
    \colhead{Range} &
    \colhead{} &
    \colhead{} &
    \colhead{} \\
    \colhead{} &
    \colhead{[kpc]} &
    \colhead{[deg]} &
    \colhead{[\%]} &
    \colhead{[\%]} \\
    \colhead{(1)} &
    \colhead{(2)} &
    \colhead{(3)} &
    \colhead{(4)} &
    \colhead{(5)} \\
    }
\startdata
0$\rightarrow$1 &  0$\rightarrow$1.07      & $83.8\pm12.5$ & $0.23\pm0.15$ & $0.25\pm0.15$  \\
1$\rightarrow$2 &  1.07$\rightarrow$2.19 & $92.4\pm7.9$   & $2.29\pm0.70$ & $1.29\pm0.74$  \\
2$\rightarrow$3 &  2.19$\rightarrow$2.45 & $24.7\pm17.6$ & $-1.55\pm0.98$ & $0.42\pm0.99$  \\
3$\rightarrow$4 &  2.45$\rightarrow$3.41 & $49.0\pm17.3$ & $-1.52\pm2.07$ & $3.07\pm2.06$  \\
\enddata
\tablenotetext{a}{These numbers refer to the stars listed in Table \ref{pol_table}, zero being the observer.}
\end{deluxetable}
```

The properties of the polarizing layers toward this direction were derived using Stokes subtraction. Table lists the derived properties for the polarizing intervals between each adjacent pairs of stars. Each interval is designated by the star identifiers in Table , zero referring to an observer on Earth. The distance span of each interval is bounded by the estimated distances of the stars, zero for an observer. The differences in Stokes parameters for the stars bounding each interval, $\Delta U_i$ and $\Delta Q_i$, and the corresponding polarization P.A. (rotated from equatorial to the Galactic coordinate system) are also listed. The corresponding degree of polarization is not listed because this quantity correlates more with the intervening dust column than with the magnetic field properties .

The derived dependence of the magnetic field orientation with distance along this direction is shown in Figure . There, the horizontal error bars represent the distance intervals listed in Table and the vertical error bars are the uncertainty in the average, dust-weighted magnetic field orientation in that interval. Out to a distance of 2 kpc, the magnetic field orientation is consistent with GPA$=90\degr$. Beyond 2 kpc, the magnetic field orientation is significantly different, being tilted approximately $45\degr$ to the direction of the Galactic plane.

From the $^{13}$CO data, it is unlikely that there is a cloud along this particular line of sight that could create a false signal. To further test this possibility the extinction per unit length (A$_v$/L) for the four stars toward this direction was calculated. If there was a dense cloud near 2 kpc there should be a large jump in A$_v$/L beyond that distance. A$_v$/L for the four stars are (in order of increasing distance): 1.88, 1.90, 1.67, and 1.46. The lack of a jump in A$_v$/L is further evidence against the presence of an intervening cloud. In fact, this result suggests that there is less extinction at further distances (probably because of increasing Galactic height).

This result is somewhat robust against the individual false-positive probabilities calculated above. The removal of any one star may increase the distance uncertainty of the change in magnetic orientation, but it would still remain detected. For example, the brightest star in Figure has the lowest probability of being a red clump star. Excluding this star from the previous analysis does not affect the presence of the change in magnetic field orientation beyond 2 kpc, but actually increases the red clump reliability from 0.743 to 0.865.

# Continuous Magnetic Field Decomposition with Many Red Clump Stars

An inherent drawback of this method, well demonstrated toward the inner Galaxy line of sight in the previous section, is the small number of identifiable red clump stars with measured NIR polarization in each $10'\times 10'$ GPIPS field. With so few stars, measurement errors could easily mimic a false signal and impede interpretation of the results, any intrinsic stellar polarization could also mask the signature of the Galactic magnetic field, and false detections would introduce noise into the results. Furthermore, unidentified dense regions along the line of sight can significantly affect interpretation of the measured polarizations. Larger numbers of red clump stars were available in a sample of starlight polarimetry drawn from and is here used to illustrate this method toward an example direction in the outer Galaxy, $(\ell , \it{b})$ = (150, -2.5).

Red clump stars were identified in this direction using the same procedure described in Section , and many more stars were found than in the previous inner Galaxy example. The larger number of red clump stars found in this outer Galaxy example is likely caused by lower extinction allowing larger distances to be probed, which can be seen in the bottom panels of Figures and . Plots of the measured polarization GPA with apparent magnitude and estimated distance are shown in Figure .

The larger number of stars available toward this outer Galaxy line of sight allows a more sophisticated analysis of changes in magnetic field orientation with distance. By fitting the change in Stokes parameters with estimated distance, the geometry of the sky-projected magnetic field can be measured without needing to invoke a specific model. Figure shows the variation of Stokes $U$ and $Q$ parameters with estimated distance. Excluding the two stars showing evidence for intrinsic polarization (discussed below), a least-squares fit to the data was calculated and is also plotted as solid lines. A statistical F-test showed that a linear fit was appropriate and that higher order fitting terms were insignificant in both cases. In Section , the entire line of sight was divided into four regions, here the sampled line of sight is divided into infinitesimally small regions and the instantaneous change in Stokes parameters is calculated from the fits in Figure . The change in Stokes parameters with distance maps directly to the orientation of the magnetic field via the equation: $${GPA_i}=\frac{1}{2} \, atan \left[ \frac{dU_i/ds}{dQ_i/ds} \right]$$ where $i$ refers to the $i$th polarizing region along the line of sight and $ds$ is an infinitesimal change in distance along this direction. For the fits in Figure , $dU/ds=-0.404 \pm 0.312$ kpc$^{-1}$ and $dQ/ds=0.497 \pm 0.582$ kpc$^{-1}$. Since these values are constants, this is interpreted as a constant magnetic field orientation in this direction between 2 and 5 kpc with equatorial PA$=160.5\degr \pm 19.6\degr$ and GPA$=120.9\degr \pm 19.6\degr$.

This result can be seen in Figure . The line of sight between the observer and the first star at 2 kpc imposes a net polarization with GPA$=99.0\degr \pm 2.1\degr$. Stars at larger distances show a weak monotonic increase in P.A. toward GPA$=118.6 \pm 3.2\degr$. The polarization of a star at a distance of 5 kpc is dominated by the polarizing layer between 2 and 5 kpc with a constant mean magnetic field orientation measured by the slopes of $dU/ds$ and $dQ/ds$.

# Identification of Intrinsically Polarized Stars

The application of this method assumes that the stars are intrinsically unpolarized. Any stars violating this assumption could be interpreted as changes in the magnetic field geometry in the polarizing medium between the stars. This could be particularly troublesome for the 'sparse' example in Section where the polarization of each star was used to measure the magnetic field properties between adjacent stars. The 'rich' example in Section uses a statistically more robust sample to fit the magnetic field properties with distance. The polarization properties of this more robust sample can be used to tentatively identify intrinsically polarized stars by identifying significant polarization deviants toward a given direction.

The 2MASS photometry and $H$-band degree of polarization of these Outer Galaxy red clump stars are shown in Figure . In the figure, the top panels show two CMDs: $K$ versus $J-K$ (upper left panel) and $K$ versus $H-K$ (upper right panel). If they are truly red clump stars then all should have the same absolute $K$-band magnitude and intrinsic colors. The overall trend in each CMD illustrates the red clump selection method in Section : more distant standard candles are fainter and redder. The lower panels plot the degree of $H$-band polarization against these same colors. In these lower panels, a linear fit to the data is shown by the solid line, which recovers the expected polarization properties for the magnetic field derived in Section , and the shaded regions represents the area enclosed by the standard error of the fit. In both bottom panels, the same two outliers are seen outside of the shaded region.

Both polarization and extinction (traced by NIR color) are integral quantities along the line of sight, meaning that the observed properties depend on the conditions between the object and the observer. While extinction monotonically increases with distance, the Stokes $U$ and $Q$ parameters are signed quantities, meaning that the degree of polarization is not necessarily a monotonic function of distance. In Figure , the photometric properties of the stars are 'well-behaved' in that there is a correlation between apparent $K$ magnitude and color, because both fundamentally depend on distance. For a line of sight with a constant magnetic field orientation, the degree of starlight polarization would also depend monotonically on distance, and therefore would also correlate with apparent brightness and reddening.

The degree of polarization has a more complicated dependence on distance if the line of sight suffers depolarization . Strong depolarization, for example occurring as a result of a change in the magnetic field direction at some distance, would cause the polarization trend to change and all stars beyond that distance would be affected. Similarly, the presence of a polarizing cloud along the line of sight would also affect all stars beyond its distance. This would not give rise to individual outliers, but to a discreet offset in polarization versus distance at that distance. A more likely explanation for individual outliers is that some intrinsic property of the star is contributing to its measured polarization.

The variation in Stokes $U$ and $Q$ parameters with $K$ are shown in Figure and can distinguish between intrinsically polarized stars and the presence of a change in field orientation or a polarizing cloud. In that figure the candidate intrinsically polarized red clump stars are marked with diamonds. A linear fit to the data, excluding the intrinsically polarized candidates, is shown by the solid line. The dashed lines represent the standard error of that fit. If there were a systematic change in the magnetic field orientation or a polarizing cloud along the line of sight there would be bend or offset in the trend of $U$ and $Q$ with distance.

The candidate stars appear as outliers in that figure as well, particularly in Stokes $U$ (lower panel). The intrinsic polarization of each star (assuming they are not interlopers) was estimated by the difference between the fit and the observed Stokes parameters. For the brighter candidate star at $K=11.01$ mag (2MASS J03534501+5026427), the intrinsic degree of polarization is $1.58\%\pm0.45\%$ with GPA$=114.0\pm14.9\degr$. The fainter candidate star (2MASS J03543283+5029199) at $K=11.75$ mag has an intrinsic polarization of $1.94\%\pm0.73\%$ with GPA$=116.3\pm21.4\degr$.

These stars were assumed to be red clump stars. However, they could be true interlopers in which case their distance estimates would be incorrect. If they were distant reddened supergiants or dwarfs then the measured polarizations could be consistent with the rest of that line of sight. No additional information is currently available in the literature for these two objects, so their true identity or the cause of this anomalous polarization remains unknown. Spectroscopic follow-up could easily determine whether these stars belong to the red clump or not. Some possible explanations for anomalous polarization of red clump stars are presented below.

# Discussion

Red clump stars were used as standard candles to probe the structure of the sky-projected magnetic field along two Galactic lines-of-sight. toward one field in the inner Galaxy, the magnetic field showed evidence for a large change in magnetic field orientation at a specific distance. toward the outer Galaxy, the magnetic field orientation was uniform in the region probed, however not parallel to the Galactic disk.

## Evidence for a Magnetic Reversal?

In Figure , the H-band polarimetry toward the inner Galaxy showed a change in the magnetic field orientation at a distance of $\sim$`<!-- -->`{=html}2 kpc. Apparently, the sky-projected magnetic field orientation has changed, but this could also be evidence for the presence of magnetic reversals in the Galactic magnetic field.

Considering the Galactic center located at the origin of a cylindrical coordinate system, the magnetic field at any point can be defined by the vector sum of radial, vertical, and azimuthal components. In the disk, theoretically and observationally , the azimuthal component dominates the total field strength in the disk of spiral galaxies, including the Milky Way.

In the Milky Way, the dominant azimuthal component of the magnetic field does not necessarily have the same sign everywhere. Evidence for discreet azimuthal magnetic field reversals have been seen in the Galactic disk in Faraday rotation measures . As described in , starlight polarimetry is insensitive to magnetic field direction, only its orientation, leading to a $180\degr$ ambiguity in the direction of the magnetic field aligning the dust. Formally, starlight polarimetry should be insensitive to full $180\degr$ magnetic reversals. However, if instead of simply reversing the azimuthal magnetic field direction, the azimuthal reversal is associated with a large-scale perturbation in the vertical and/or radial magnetic field components, then the projection of that field onto the sky would also change and the effects of magnetic reversals would be expected to be seen with starlight polarimetry.

Magnetic reversals are expected in two main situations: (1) there is a field discontinuity, separated by a current sheet , or (2) if the magnetic field is continuous, then the amplitude of the azimuthal magnetic field must decrease to zero before increasing in the opposite direction, though no limits are placed on the vertical and radial field components that may also change, possibly in such a way that the total field strength is constant. In the first situation the magnetic field direction is reversed on opposite sides of a current sheet. Magnetic field rotations of $180\degr$ have been observed by in situ satellites passing through the heliospheric current sheet in the solar system . The presence of a current sheet would have little effect on the observed polarization orientation since the magnetic fields on either side of an interstellar current sheet would be indistinguishable with this NIR background starlight observational tool, illustrated in Figure 1 of . In the second situation, as the dominant azimuthal magnetic field vanishes the other field components (radial and vertical) will be expressed. This would manifest as a change in the linear polarizing properties of the medium hosting this magnetic field.

A change in the magnetic field geometry toward the inner Galaxy has been modeled from observations of RMs of Galactic pulsars and polarized extragalactic objects. toward the direction shown in Figure ($\ell=20\degr$, $b=0\degr$), saw evidence for magnetic reversals at Galactocentric radii (R) of 7.5, 6, and 5 kpc (corresponding to line of sight distances of 1.1, 2.7, and 3.9 kpc, respectively) in their ASS+RING Galactic magnetic field model. In this same direction concluded that a magnetic reversal occurs at R=5.8 kpc (a line of sight distance of 3.0 kpc). It seems that the same magnetic feature seen in the RMs is being observed by NIR starlight polarimetry. If so, this would corroborate the presence of a large-scale magnetic field change toward that direction.

For now, the NIR starlight polarimetry results, while intriguing, are tentative. Spectroscopic follow-up of the stars along this line of sight is required to derive reliable spectrophotometric distances, but will require a substantial observational program to test the photometric results presented here.

## Magnetic Fields Following the Galactic Warp

Figure shows a line of sight toward ($\ell$, $b$) = ($150\degr$, $-2.5\degr$). Toward this direction in the outer Galaxy, the magnetic field orientation was shown to be uniform with distance, which agrees with the observationally-driven models of and . Observations and simulations conclude that in this direction the disk magnetic field should be roughly parallel the Galactic plane .

However, the deviation of the polarization P.A., as described in Section (GPA$=121\degr$), is inconsistent with a plane-parallel magnetic field in the Galactic disk (GPA$=90\degr$). Geometric effects cannot account for the discrepancy. This particular outer Galaxy line of sight samples the Galactic magnetic field below the disk ($b=-2.5$), and the combination of inclined viewing angle and the pitch angle of the Galactic disk magnetic field only contribute a 1.1$\degr$ twist to the projected magnetic field orientation for stars within $\sim 7$ kpc.

A possible explanation for the observed magnetic field orientation in the Outer Galaxy is the presence of the Galactic warp . Using the model of , the Galactic warp can be described by $z(\phi)\sim sin(\phi+5\degr)$ where $\phi = 0\degr$ is toward the Sun and $\ell=150\degr$ spans $\phi=0-12.3\degr$ from the Sun to a distance of 4 kpc. For these values of $\phi$, the Galactic plane is tilted upward $44\degr$ toward Galactic North. This leads to a predicted magnetic field orientation of GPA=$134\degr$, consistent with the measured value (GPA$=121\pm20\degr$). A similar result was seen in M33 where vertical magnetic fields were suggested to actually be disk fields perturbed by a Galactic warp .

## Intrinsically Polarized Red Clump Stars

One unexpected result of this study was the possible identification of intrinsically polarized red clump star candidates. These stars were photometrically identified as red clump stars, but their polarization properties (shown in Figures and ) were inconsistent with the trend of Stokes $U$ and $Q$ versus brightness (or distance) seen in the other red clump stars. The reason for intrinsic polarization is unclear from the data available, though there are several possibilities. Rapid rotation in old red clump stars is unlikely, but its presence could distort the star from spherical symmetry leading to a net polarization signature . This rotation would be easily seen in moderate resolution spectroscopy and also allow identification of the star's rotation axis. The presence of a disk or envelope could also give rise to a net polarization . If true, the observed NIR colors would suggest that this is not a red clump star but possibly a reddened main sequence star that is observed in the same color-magnitude space as true red clump stars. Future variability studies or spectroscopy could shed light on the physical mechanism creating these apparently anomalous polarizations.

The method described above can be used to measure the mean interstellar polarization signal with distance for some direction, then easily distinguish polarimetric deviants from stars following the nominal polarization trend. This allows for the possibility of surveying for intrinsically polarized stars over large areas of the sky. Intrinsic polarization is associated with a range of physical phenomena that may not otherwise be seen in large-scale photometric surveys. By looking for polarization deviants, a catalog of interesting targets can be developed for follow-up investigation.

# Summary and Conclusions

This work presents a new method for measuring the structure of the large-scale Galactic magnetic field. Photometrically identified red clump stars were used as standard candles along two example directions in the sky, one with few red clump stars and one with many red clump stars. The false-detection probability as a function of $K$-band brightness when using the method of was calculated and used to inform the reliability of the results. H-band starlight polarimetry of these stars from GPIPS provides information about the integrated magnetic field properties toward each star. By decomposing the Stokes $U$ and $Q$ parameters with distance, changes in the magnetic field orientation in the plane of the sky were measured with distance.

One line of sight, toward $(\ell , b)= (19.49\degr,+0.56\degr)$ in the inner Galaxy, shows evidence for a strong magnetic field perturbation at a distance of 2 kpc that cannot be attributed to dusty clouds along the line of sight. This is approximately the same distance as a magnetic field reversal seen in Faraday rotation surveys , which are sensitive to the electron-weighted magnetic field. Formally, starlight polarimetry should not be sensitive to magnetic field reversals , but changes in large-scale magnetic field geometry could manifest as both NIR polarization P.A. changes and as magnetic field reversals in Faraday rotation studies. While this remains circumstantial evidence for the presence of a magnetic field reversal toward this direction, the overall agreement between Faraday rotation studies of the Galactic magnetic field and this method is intriguing.

A different line of sight, toward $(\ell , b)= (150\degr,-2.5\degr)$ in the outer Galaxy, shows a constant magnetic field orientation between distances of 2-5 kpc. The larger number of identified red clump stars allowed a more detailed analysis of the relatively simple magnetic field in this direction. The measured magnetic field orientation (GPA=$121\degr$) does not fit the plane-parallel model typically assumed for the Galactic plane. However, this value can be explained by the Galactic warp toward this direction in the Galaxy. From this analysis, it appears that the large-scale Galactic magnetic field toward this direction follows the Galactic warp.

An unexpected outcome of the analysis toward the outer Galaxy was the possible detection of two intrinsically polarized red clump stars. Photometrically, these stars fall along the red clump sequence, but their polarization properties do not match the growth of percentage polarization with distance exhibited by the remaining stars. The deviations of the NIR polarization observations of these stars to the fits of the other stars' Stokes parameters allowed determination of these intrinsic polarizations. Moderate resolution spectra (to test whether these are red clump stars and to measure $v$ sin $i$) would provide additional insight into the physical natures of these objects.

This work represents a first step in decomposing with distance the magnetic fields probed by background starlight polarimetry. Additional work is required to fully develop this method. Photometric methods can statistically identify red clump stars; however, spectroscopic follow-up is needed. Also, while the relative distances, or at least distance order, of stars derived from photometry are reliable, physical distances require spectrophotometric follow-up. The future application of photometrically identified red clump stars to large, uniformly sampled background starlight polarimetry will allow the large-scale Galactic magnetic field to be unambiguously disentangled with distance for the first time.

The author would like to thank Dan P. Clemens, Brian W. Walsh, Kamen Kozarev, Daniel Jaffe, and Neal Evans for useful discussions. The author would also like to thank the anonymous referee for comments that greatly improved the manuscript. The author gratefully acknowledges financial support from McDonald Observatory. This research used data from the Boston University (BU) Galactic Plane Infrared Polarization Survey (GPIPS), funded in part by NSF grants AST 06-07500 and 09-07790. GPIPS used the Mimir instrument, jointly developed at BU and Lowell Observatory and supported by NASA, NSF, and the W.M. Keck Foundation.


abstract: We introduce a new scheme for generating short-wavelength coherent broadband radiation with well-controlled spectral characteristics. It is based on shaping long-wavelength femtosecond pulse to coherently control atomic resonance-mediated (2+1) three-photon excitation to a broad far-from-resonance continuum. Here, the spectrum (central frequency and bandwidth) of deep-ultraviolet coherent broadband radiation generated in Na vapor is experimentally controlled by tuning the linear chirp we apply to the driving phase-shaped near-infrared femtosecond pulse. This is a first step in implementing the full scheme for producing shaped femtosecond pulses at wavelengths down-to the vacuum-ultraviolet range.
author: Leonid Rybak, Lev Chuntonov, Andrey Gandman, Naser Shakour; Zohar Amitay
title: Multiphoton Femtosecond Control of Resonance-Mediated Generation of Short-Wavelength Coherent Broadband Radiation

Femtosecond coherent control is a field of both scientific and technological importance concerning a variety of physical and chemical photo-induced processes . It utilizes the broad coherent spectrum of femtosecond pulses, which provides a manifold of excitation pathways, to control the photo-induced dynamics and outcome in ways that other means cannot. The control handle is the pulse shape that is determined by the spectral characteristics of the broadband pulse, i.e., the amplitude, phase, and polarization of its different frequency components .

Indeed, over the past decade, coherent control using shaped femtosecond pulses has been successfully demonstrated in the laboratory with many systems and processes . However, although very successful, these experiments have been limited to atomic and molecular species that absorb visible (VIS) or near-infrared (NIR) light. The reason is the experimental difficulties in shaping femtosecond pulses outside these spectral regions, where the devices commonly-used for flexible pulse shaping (liquid-crystal modulators and TeO$_{2}$ acousto-optic modulators) are inappropriate. Extending the techniques to generate shaped pulses at short wavelengths down to the deep-ultraviolet (UV) and vacuum-ultraviolet (VUV) regions will significantly extend the variety of molecules accessible for coherent control, since most electronic transitions are in this spectral range.

Recently, significant progress has been made in producing shaped UV pulses based on generating the UV pulse from a source NIR femtosecond pulse via (nonresonant) nonlinear optical processes in crystals . For the UV pulse shaping the corresponding works have followed two approaches. One approach is of direct UV shaping , where the generated UV pulse undergoes direct shaping using new types of special devices that are appropriate to the UV spectral range. The other approach is of indirect UV shaping , where the NIR-to-UV nonlinear optical conversion is done with a (conventionally) shaped NIR pulse and, hence, the generated UV pulse is already a shaped one. The output shaped broadband UV spectrum is determined by the input shaped NIR spectrum following its NIR-to-UV nonlinear optical transformation. However, even though all these studies provide a feasible pathway to UV pulse shaping, their basic use of optical crystals actually limits the generated pulse to UV wavelengths above 200-250 nm.

In this Letter, we pursue a line of action that is generally applicable for generating short-wavelength shaped pulses down-to the VUV spectral range. It follows the above-mentioned indirect shaping approach, however, instead of using optical crystals, the shaped short-wavelength pulse is generated from a shaped long-wavelength pulse via an atomic resonance-mediated (2+1) three-photon excitation to a broad far-from-resonance continuum. The resulting transient polarization drives the generation of the shaped short-wavelength coherent broadband radiation. Resonance-mediated three-photon excitation in atomic vapors has actually been employed very successfully in the past for efficient and tunable generation of short-wavelength radiation, however it has been used only for generating narrowband (nanosecond or picosecond) UV and VUV coherent radiation by narrowband NIR and VIS excitation (for example, see ); The present broadband nature introduces a completely new dimension to the process. Our new scheme also utilizes previous works on femtosecond broadband coherent control of resonance-mediated (2+1) three-photon absorption , extending them from a case of a single (real) final state to a case of a continuous manifold of multiple (virtual) final states. Here, as a first step in implementing the full scheme, the spectrum (central frequency and bandwidth) of deep-ultraviolet coherent broadband radiation generated in Na vapor is experimentally controlled by tuning the linear chirp that we apply to the driving phase-shaped near-infrared femtosecond pulse.

Figure (a) shows the atomic broadband excitation scheme considered here, involving an initial ground state $\left|g\right\rangle$ and an excited state $\left|r\right\rangle$ of one symmetry and a manifold of excited states $\left|v\right\rangle$ of another symmetry. Irradiation with a shaped long-wavelength femtosecond pulse induces a resonance-mediated (2+1) three-photon excitation from $\left|g\right\rangle$ via $\left|r\right\rangle$ to a broad far-from-resonance continuum, resulting in a time-dependent atomic polarization that drives the generation of a corresponding shaped short-wavelength coherent broadband radiation. The last coherent emission step is associated with a corresponding de-excitation to $\left|g\right\rangle$. Along this whole photo-induced process the $\left|v\right\rangle$ states are accessed nonresonantly and provide the $\left|g\right\rangle$-$\left|r\right\rangle$ nonresonant two-photon couplings at the one- and three-photon excitation energies.

Expressing the temporal atomic polarization using third-order time-dependent perturbation theory and transforming it to the frequency domain, we have obtained the generated short-wavelength spectral field $E_{out}(\Omega)$ at frequency $\Omega$ to be given by $$\begin{aligned}
E_{out}(\Omega) & \propto & \mu^{2}_{g,r} \chi^{2}_{r,g}
\left[A^{(2+1)on-res}(\Omega) + A^{(2+1)near-res}(\Omega)\right] \; ,
\label{eq1} \\
%\end{eqnarray}
%where
%\begin{eqnarray}
A^{(2+1)on-res}(\Omega) & = & i \pi A^{(2)}(\omega_{r,g}) E(\Omega-\omega_{r,g}) \; ,
\label{eq2} \\
%\end{eqnarray}
%\begin{eqnarray}
A^{(2+1)near-res}(\Omega) & = & -\wp\int_{-\infty}^{\infty}\frac{1}{\delta}A^{(2)}(\omega_{r,g}-\delta)E(\Omega-\omega_{r,g}+\delta)d\delta \; ,
\label{eq3} \\
%\end{eqnarray}
%with %and
%\begin{eqnarray}
A^{(2)}(\xi) & = & \int_{-\infty}^{\infty} E(\omega)E(\xi-\omega)d\omega \; , %\hspace{2.75cm}
\label{eq4}
\end{aligned}$$ with $\omega_{r,g}$ being the $\left|g\right\rangle$-$\left|r\right\rangle$ transition frequency. The $E(\omega)$$\equiv$$\left|E(\omega)\right|\exp{[i\Phi(\omega)]}$ is the long-wavelength spectral field of the excitation pulse, with $\left|E(\omega)\right|$ and $\Phi(\omega)$ being, respectively, the spectral amplitude and phase of frequency $\omega$. For the (unshaped) transform-limited (TL) pulse, $\Phi(\omega)$=0. The $\mu^{2}_{r,g}$ and $\chi^{2}_{g,r}$ are the $\left|g\right\rangle$-$\left|r\right\rangle$ nonresonant two-photon couplings (due to the $\left|v\right\rangle$ states) at the $\omega_{c}$ and $\omega_{r,g}+\omega_{c}$ excitation regions, with $\omega_{c}$ being the carrier frequency of $E(\omega)$. The spectral intensity generated at $\Omega$ is $I_{out}(\Omega)\propto\left|E_{out}(\Omega)\right|^{2}$, and the total short-wavelength yield is $Y_{out}=\int_{-\infty}^{\infty}I_{out}(\Omega)d\Omega$.

As can be concluded from the above frequency-domain formulation and is illustrated in Fig. (a), the short-wavelength field $E_{out}(\Omega)$ generated at a frequency $\Omega$ results from the interferences among all the possible resonance-mediated three-photon excitation pathways from $\left|g\right\rangle$ to the excitation energy corresponding to $\Omega$. Each of these pathways is either on or near resonance with the intermediate state $\left|r\right\rangle$. A pathway of detuning $\delta$ involves a nonresonant absorption of two photons with a frequency sum of $\omega_{r,g}$$-$$\delta$ and the absorption of a third complementary photon of frequency $\Omega$$-$$(\omega_{r,g}$$-$$\delta)$. The on- ($\delta$=0) and near-resonant ($\delta$$\ne$`<!-- -->`{=html}0) pathways reaching the energy of $\Omega$ are interfered separately in $A^{(2+1)on-res}(\Omega)$ and $A^{(2+1)near-res}(\Omega)$ ($\wp$ is the Cauchy principal value). They are expressed using the parameterized amplitude $A^{(2)}(\xi)$ interfering all the nonresonant two-photon pathways of frequency sum $\xi$, each composed of photon pair $\omega$ and $\xi$$-$$\omega$. As the exciting long-wavelength field $E(\omega)$ controls the amplitudes and phases of the different interfering three-photon pathways, it also controls the amplitude and phase of the short-wavelength field $E_{out}(\Omega)$ generated at each $\Omega$.

The model system of the study is the sodium (Na) atom, with the $3s$ ground state as $\left|g\right\rangle$, the $4s$ state as $\left|r\right\rangle$, and the manifold of $p$ states as $\left|v\right\rangle$ \[see Fig. (a)\]. The $3s$-$4s$ transition frequency is $\omega_{r,g} \equiv \omega_{4s,3s} = 25740$ cm$^{-1}$ corresponding to two 777-nm photons. The sodium is irradiated with phase-shaped linearly-polarized near-infrared (NIR) femtosecond pulses having a spectrum peaked at $\omega_{0}$=12876 cm$^{-1}$ (776.6 nm) with 88-cm$^{-1}$ bandwidth (FWHM), as shown in Fig. (b). The pulse energy is 10 $\mu$J corresponding to a TL peak intensity of $I_{TL}$$\sim$`<!-- -->`{=html}3$\times$`<!-- -->`{=html}10$^{10}$ W/cm$^{2}$. The NIR spectrum results from blocking an original spectrum at its lower-frequency edge to prevent three-photon resonant access to the $7p$ state \[$\omega_{7p,4s} = 12801$ cm$^{-1}$ (781.2 nm)\], which is undesired for generating short-wavelength radiation of broadband nature. The interaction with the NIR excitation pulse leads to the generation of a coherent broadband UV radiation. Experimentally, a sodium vapor in a heated cell at 900 K with Ar buffer gas is irradiated with the NIR laser pulses after they undergo shaping in a setup incorporating a pixelated liquid-crystal spatial light phase modulator . The shaping resolution is 2.05 cm$^{-1}$. The coherent UV radiation emitted in the propagation direction of the NIR beam is measured using a UV-spectrometer coupled to a time-gated camera system. The corresponding overall UV spectral resolution is 35 cm$^{-1}$ (0.23 nm), with 5.8-cm$^{-1}$ spectral width of each camera pixel.

Figure (c) shows the measured coherent UV spectrum generated by the NIR TL pulse. As seen, it is of broadband nature with a bandwidth of 125 cm$^{-1}$ around its maximum at $\Omega$=38642 cm$^{-1}$ (258.78 nm). After deconvoluting the measurement resolution, the UV spectral bandwidth amounts to 120 cm$^{-1}$ ($\sim$`<!-- -->`{=html}125-fs TL duration). The peak UV frequency of 38642 cm$^{-1}$ corresponds to $\omega_{4s,3s}$ plus a frequency of 12902 cm$^{-1}$ that is around the peak NIR frequency $\omega_{0}$. As seen from Eqs. ()-(), the on-resonant amplitude $A^{(2+1)on-res}(\Omega)$ is by-definition maximal at $\Omega$=$\omega_{4s,3s}$+$\omega_{0}$, while the peak location of the near-resonant amplitude $A^{(2+1)near-res}(\Omega)$ can deviates from it according to the specific NIR excitation spectrum.

In terms of the absolute efficiency of the UV generation process, under our experimental conditions, at a temperature of 900 K corresponding to Na density of 3.7$\times$`<!-- -->`{=html}10$^{17}$ cm$^{3}$, the total energy of the emitted coherent broadband UV radiation reaches about 1% of the NIR pulse energy for the TL excitation of $I_{TL}$$\sim$`<!-- -->`{=html}3$\times$`<!-- -->`{=html}10$^{10}$ W/cm$^{2}$. As expected for a coherent generation process, the UV generation efficiency exhibits a quadratic dependence ($n_\textrm{Na}^{2}$) on the Na density $n_{\textrm{Na}}$ (controlled by the cell temperature). Up to the NIR TL peak intensity of $I_{TL}$$\sim$`<!-- -->`{=html}3$\times$`<!-- -->`{=html}10$^{10}$ W/cm$^{2}$ the UV generation efficiency exhibits a cubic dependence ($I_{TL}^{3}$) on $I_{TL}$, indicating the third-order perturbative regime associated with Eqs. ()-(). Higher values of $I_{TL}$ lead to a deviation from this regime, which we avoid in order to maintain the validity of the relatively simple theoretical description.

As a first control study on resonance-mediated generation of short-wavelength coherent broadband radiation, we demonstrate control over the UV generation with linearly-chirped NIR pulses shaped with quadratic spectral phase patterns $\Phi(\omega)=\frac{1}{2} \kappa_{chirp}(\omega-\omega_{0})^{2}$ of variable chirp parameter $\kappa_{chirp}$. The $\kappa_{chirp}$ value is tuned here from negative to positive values, with the zero value $\kappa_{chirp}$=0 corresponding to the NIR TL pulse.

Figure (a) presents the complete set of the experimentally generated chirp-controlled UV spectra. The results are presented as a color-coded map of the spectral intensity $I_{out}(\Omega; \kappa_{chirp})$ measured at each UV frequency $\Omega$ as a function of the $\kappa_{chirp}$ value applied to the generating linearly-chirped NIR pulse. Each UV spectrum at a particular value of $\kappa_{chirp}$ is normalized (independently) by its maximal intensity. The experimental chirp dependence of the peak frequency (i.e., the frequency of maximal intensity) and bandwidth (FWHM) of the generated UV spectrum, as obtained from the map data, are shown in Figs. (b) and (c) (circles). Also shown there are corresponding theoretical results (lines), which are the outcome of numerical calculations based on Eqs.()-() that have been followed by a convolution with the UV spectral measurement resolution. As seen, the theoretical results describe very well the experimental ones.

As $\kappa_{chirp}$ is experimentally tuned from negative to positive values, the peak UV frequency $\Omega_{0}$ changes from around 38580$-$`<!-- -->`{=html}38585 cm$^{-1}$ when $\kappa_{chirp}$$\le$$-$`<!-- -->`{=html}0.10 ps$^{2}$ to around 38635$-$`<!-- -->`{=html}38640 cm$^{-1}$ when $\kappa_{chirp}$$\ge$`<!-- -->`{=html}0.10 ps$^{2}$, with a value of 38642 cm$^{-1}$ when $\kappa_{chirp}$=0. The change of $\Omega_{0}$ takes place mostly at the region of small chirp values around zero chirp. A line of reference, which we refer to further below, is indicated in Figs. (a)-(b) at the UV frequency of 3$\times$$\omega_{4s,3s}/2$ (38610 cm$^{-1}$). The measured UV bandwidth $\Delta\Omega$ amounts to about 125 cm$^{-1}$ around zero chirp and decreases continuously down-to 55$-$`<!-- -->`{=html}75 cm$^{-1}$ at the larger negative and positive chirps of $|\kappa_{chirp}|$$\ge$`<!-- -->`{=html}0.1 ps$^{2}$. After deconvoluting the experimental resolution, the corresponding values are about 120 cm$^{-1}$ and 40$-$`<!-- -->`{=html}65 cm$^{-1}$, respectively. The chirp dependence of $\Delta\Omega$ is nearly symmetrical around $\kappa_{chirp}$=0. As explained above, all these specific chirp-dependent values of $\Omega_{0}$ and $\Delta\Omega$ and their tunability range result from the specific NIR excitation spectrum used here.

Highly-intuitive qualitative understanding of these chirp dependencies is actually obtained by combining the frequency-domain picture of Eqs. ()-() with the picture of the so-called joint time-frequency domain. Under this framework, zero chirp (TL pulse) corresponds to the simultaneous arrival of all the spectral pulse frequencies, while positive and negative chirps correspond to the cases where lower frequencies arrive, respectively, before or after higher frequencies. The temporal separation between two different frequencies increases as the chirp magnitude increases. With regard to the Na excitation, the resonance-mediated (2+1) three-photon excitation to a given total energy of $\Omega$ is expected to be of large amplitude, leading to high UV intensity at $\Omega$, once the third-photon absorption occurs shortly after the on- or near-resonant excitation of the intermediate $4s$ state. Based on Eq. () of $A^{(2)}$($\xi$$\approx$$\omega_{4s,3s}$), such $4s$ excitation occurs by absorbing all the NIR photon pairs having their frequency-sum close to $\omega_{4s,3s}$, i.e., their frequencies are symmetric around $\sim$$\omega_{4s,3s}$/2, and thus can be associated with a narrow time window around the arrival time of the pulse frequency $\omega_{4s,3s}$/2. Hence, overall, efficient resonance-mediated excitation is expected for the energies of $\Omega \approx \omega_{4s,3s}+\omega_{3}$ with $\omega_{3}$ either larger or smaller than $\omega_{4s,3s}/2$ according to whether the NIR chirp is positive or negative, respectively. Consequently, the positive and negative NIR chirps are expected to lead to efficient UV generation mostly at UV frequencies that are either larger or smaller than 3$\times$$\omega_{4s,3s}/2$, respectively. Also, since larger chirp magnitude leads to larger temporal spread of the pulse frequencies such that less frequencies arrive within the time window of efficient excitation, the generated UV bandwidth is expected to be the largest for NIR chirps around zero and to decrease toward a constant level as the chirp magnitude increases. Indeed, as seen in Figs. (a)-(c), the observed chirp-dependent results follow this qualitatively-predicated behavior.

Last, for completing the picture, Fig. (d) presents experimental (circles) and theoretical-numerical (lines) results for the chirp dependence of the total generated UV yield $Y_{out}$. The results are normalized by the yield generated by the NIR TL pulse. As $\kappa_{chirp}$ changes, the total UV yield experimentally changes continuously from about 10% to 105% of the TL-generated UV yield. From an applicative point of view, such a change in $Y_{out}$ is not optimal for using the generated shaped short-wavelength broadband radiation. Nevertheless, having in hand a corresponding calibration curve, such as the one shown in Fig. (d), removes this inconvenience. Also here, similar to the above spectral UV results, the experimental results of $Y_{out}$ are described very well by the theoretical ones.

In conclusion, we have introduced a new scheme for generating short-wavelength coherent broadband radiation with well-controlled spectral characteristics. It utilizes shaped long-wavelength femtosecond pulses to coherently control atomic resonance-mediated multiphoton excitation to a broad far-from-resonance continuum. As compared to a completely nonresonant atomic excitation, the resonance-mediated nature of the excitation provides a higher degree of control over the properties of the generated radiation together with a higher generation efficiency. In the present work, by tuning the linear chirp applied to the exciting phase-shaped NIR femtosecond pulse we have experimentally controlled the spectrum of the coherent broadband UV radiation generated via resonance-mediated (2+1) three-photon excitation of Na. The coherent UV bandwidth corresponds to a femtosecond UV pulse. This is a first step in implementing the full scheme for producing shaped femtosecond pulses at wavelengths down-to the VUV range, which currently is inaccessible by the present pulse shaping techniques. Inspired by the past demonstrations of narrowband coherent VUV generation (see above), the extension to the VUV range will be achieved by a proper change of the excited atomic species and the driving pulse wavelength and/or by utilizing resonance-mediated (N+1) excitations of higher order (i.e., N$>$`<!-- -->`{=html}3). Such shaped pulses will make many VUV-absorbing molecular species newly available for coherent control. This is important for nonlinear spectroscopy and microscopy as well as for chemical reaction photo-control.