author: Pushpa Khare
*Physics Department, Utkal University, Bhubaneswar, 751004, India*
and
S. Ikeuchi
*Department of Earth and Space Science, Osaka University, Toyonaka, Osaka 560, Japan*
title: The ionization and abundance of C and Si in QSO absorbers

C and Si in QSO absorbers

# Introduction

Quasar absorption lines have proved to be extremely useful probes of the high redshift Universe. In recent years important information has been obtained, among other things, about the chemical abundances in galaxies at high redshifts. Evidence for evolution in chemical abundances in these galaxies, the abundance increasing with decreasing redshift, has been obtained (i) directly, through abundance determination in damped Lyman alpha systems (DLAS) at different redshifts ( Pettini etal 1995; Lu etal 1996), and (ii) indirectly, from the variation of the number of C IV systems and Lyman limit systems (LLS) per unit redshift interval, per line of sight as a function of redshift ( Steidel, Sargent & Boksenberg, 1988; Khare & Rana, 1993). The absorption line studies have also given clues for understanding the history of stellar nucleosynthesis (Lauroesch etal 1996). Lu etal (1996), through an analysis of the observed abundances of 23 DLAS, have shown that the relative abundance patterns of several elements are consistent with their formation in Type II supernovae. In particular, they find evidence for \[Si/Fe\] $\simeq$ 0.4, which is very similar to the overabundance of Si found in the Galactic halo stars. Evidence for overabundance of Si with respect to C, by a factor of three, has also been found in Lyman alpha forest clouds (LAFCs) (Songaila & Cowie, 1996; hereafter SC96) as well as in some intervening and associated systems (Petitjean, Rauch & Carswell 1994, Savaglio etal 1996).

The shape of the UV background radiation is important in deciding the ionization balance of various elements in QSO absorbers. This radiation most likely originates in AGNs, however, a significant contribution from young star forming galaxies can not be ruled out (Bechtold etal 1987; Bajtlik, Duncan $\&$ Ostriker, 1988, Madau & Shull, 1996). The shape of the background radiation not only depends on the relative contribution from these two sources but is also decided by the absorption by the material giving rise to QSO absorption lines. Recent observation of a significant Gunn-Peterson optical depth ($\tau\ge$`<!-- -->`{=html}1.7 for redshift $>$ 3.0) at the wavelength of the Lyman alpha line of He II towards a QSO (Jackobsen etal 1994) indicates the presence of a large break by factors $>$ 25 (Madau & Meiksin, 1994) at the He II ionization edge. Similar break was also found necessary in order to explain the observations of column density of C and Si ions in the Lyman alpha forest lines (SC96) as well as in some heavy element systems (Savaglio etal 1996). The shape of the background UV field due to AGNs, at various redshifts, taking into account the absorption and reemission from the QSO absorbers has been recently determined by Haardt & Madau (1996).

A lower He II optical depth ($\tau\simeq$`<!-- -->`{=html}1.0) at somewhat lower redshift (z $\simeq$ 2.5) has been observed by Davidsen etal ( 1996) towards HS1700+64, indicating a higher degree of ionization of He below z=3.0. Evidence for the lowering of the optical depth of the Universe to the He II ionizing photons below z=3.1 has also been obtained by SC96, from the observations of the column density ratios of Si IV and C IV in the LAFCs. They find an abrupt change in the values of this ratio at z = 3.1, the ratio being higher at higher redshifts. As the absorbers producing heavy element systems in the QSO spectra are also most likely ionized by the intergalactic UV background (Srianand & Khare, 1995), it may be of interest to look for signatures of a change in opacity of the Universe in the absorption line data of heavy element systems.

A correlation with redshift of the ratio, of equivalent widths of lines of Si IV ($\lambda$`<!-- -->`{=html}1393) and C IV ($\lambda$`<!-- -->`{=html}1548) in the heavy element systems, observed at intermediate resolution, the ratio increasing with increasing redshift, was found by Bergeron & Ikeuchi (1990). High resolution data (FWHM $\le$ 23 km s$^{-1}$) for the column densities of C IV and Si IV are now available for a number of QSOs. In this paper we present the analysis of these data in order to understand the ionization of these elements at different redshifts as well as to study the overabundance of Si with respect to C in these systems. In section 2 we investigate the presence of a correlation between ratios of column densities of Si IV and C IV and redshift. In section 3 we present the results of photoionization models and compare these with the observations. Conclusions are presented in section 4.

# Correlation Analysis

We have collected from the literature (Cristiani etal 1995, Savaglio etal 1996, Petitjean, Rauch & Carswell 1994, Petitjean & Bergeron 1994, Giallongo etal 1993, Fan & Tytler 1994, Tripp, Lu & Savage 1996, Prochaska & Wolfe 1996, Wampler 1991, Wampler, Petitjean & Bergeron 1993) the column densities of C IV and Si IV observed in the heavy element systems towards several QSOs observed with FWHM between 8 to 23 km s$^{-1}$. The total sample consists of 30 non-DLA intervening systems including 10 for which only upper limits on the Si IV column density are available. In addition we have 23 components of intervening DLAS.

We note that the data are rather inhomogeneous in the sense that the observations have been made at somewhat different resolutions and with different S/N values. This in principle can introduce an incompleteness in the sample as the lower limit for detection of lines for different QSOs may be different. This could affect the correlation analysis if the incompleteness introduces a redshift dependence of the minimum detectable column density, as a result, introducing an artificial redshift dependence of the column densities in the sample. In order to check this we performed Spearman rank correlation tests for redshift dependence of column densities of C IV for the non-DLA and DLA intervening systems separately. No such correlation was found, the chance probabilities being 0.53 and 0.38 respectively. Similar exercise for column densities of Si IV also rules out any correlation with redshift, the chance probabilities being 0.84 and 0.74 respectively. We thus believe that our data, though inhomogeneous, are not biased and can be used for the correlation analysis.

In Fig.1, we have plotted the ratios of column densities of Si IV and C IV, R, and the upper limits on R as a function of redshift for both categories of intervening absorption lines mentioned above. The non-DLA intervening systems are believed to be produced by the gas in the halos of absorbing galaxies and are most likely irradiated by intergalactic UV background radiation (Srianand & Khare, 1995). Most DLAS have components with high as well as low ionization states. These presumably arise in the halo and disk components of the absorbing galaxy (Wolf etal 1996). Lu etal (1996) have pointed out that the Si IV and C IV absorption profiles of DLAS always resemble each other and have a different appearance from the low-ion absorption lines. They suggest that the bulk of the high-ions could arise from the halo clouds. The DLA components showing high-ions thus, very likely, belong to the same population as that of the non-DLA intervening systems. We can therefore combine the non-DLA intervening systems and the DLA components showing high-ions together, while looking for the redshift dependence of column density ratios of high-ions.

Spearman rank correlation tests for both categories of intervening systems taken together as well as taken separately rule out any correlation between R and z, the values of chance probability being 0.61, 0.229 and 0.213 for all intervening systems, non-DLA intervening systems and DLA intervening systems respectively. Generalized Kendall test (Isobe, Feigelson & Nelson, 1986) applied to the 53 values, including upper limits, for all intervening systems and separately to the 30 values, including upper limits, for the non-DLA intervening systems gives probability of 0.684 and 0.174 respectively, which again indicate the absence of any correlation. We also tried to study the correlation for systems restricted to smaller redshift ranges. No correlation was however seen. KS tests rule out any abrupt change in R values at any particular redshift. Note that Spearman rank correlation test applied to the values of R observed in LAFCs by SC96 gives chance probability equal to 4.27$\times10^{-3}$, showing a good correlation.

# Ionization State of the Absorbers and the Background Field

For 14 systems ( 7 DLA intervening, 5 non-DLA intervening and 2 associated systems) the column density of Si II is also known. We can thus try to investigate the shape of the UV radiation field as well as the Si overabundance with respect to C in these systems. The neutral hydrogen column density is not known for several of these systems. However, as will be seen below, the exact value of this column density is not very important. We have constructed a grid of photoionization models using the code 'cloudy 84' written by Prof. G. Ferland, for three values of neutral hydrogen column density, N$\rm_{H\;I}$, typical of DLAS (10$^{20}$ cm$^{-2}$), of LLS (3.0$\times$`<!-- -->`{=html}10$^{17}$ cm$^{-2}$) and of high ionization intervening systems ($\le 10 ^{16}$ cm$^{-2}$), for different shapes of UV radiation field. Heavy element abundance was assumed to be one thirtieth of the solar value, with solar values of relative abundances of different chemical elements. The results, which we present in the form of column density ratios, however, are not sensitive to the heavy element abundances as well as the particle density. The ratios are also independent of N$\rm_{H\; I}$ for N$\rm_{H\; I}\;\le\; 10^{16}\; cm\;s^{-1}$. The ratios for these three values of N$\rm_{H\; I}$ bracket the ion ratios for all the intervening and associated systems considered here.

Heavy element absorbers have complex structures and it has been argued that the lines of different ions may be produced in different regions of the absorbers and also it is possible that a hot collisionally ionized phase may exist in the absorbers (Giroux, Sutherland & Shull, 1994). However, here, we are restricting our analysis to high resolution observations. We can therefore assume that lines produced in individual clouds have been resolved and the the results of 'cloudy' models, which take into account the radiation transfer inside a cloud, can be applied to the column densities seen in individual clouds. Also as we are restricting to the ions of Si II, Si IV, C II and C IV only, the contributions from the collisionally ionized phase may not be very important.

In Fig 2 we have plotted R vs. the column density ratios of Si II to Si IV, S, for the three values of neutral hydrogen column densities mentioned above. Fig 2a, 2b and 2c are for different shapes of the background UV radiation field; these are (a) power law with a slope of -1.5, which corresponds to typical unprocessed AGN spectra, (b) AGN background filtered through the intervening galactic and intergalactic absorbers as given by Haardt & Madau (1996, hereafter HM96) for the redshift of 2.5 and (c) power law with slope of -1.5 and with a cutoff at 4 Ryd, which will be the case if the He II ionization fronts of different QSOs have not yet overlapped at the redshift of the absorbers as suggested by SC96. We have also plotted the observed ratios with their error bars. The observed values have large errors. These are essentially a result of the fact that the lines are often saturated and profile fitting analysis can yield acceptable fits to the observed profiles over a range of column density values. However, as will be seen below, definite conclusions regarding the overabundance and the background can be drawn from the comparison of the observed values with the results of the photoionization models. As seen from the figure the shape of the spectrum makes significant difference to the lower values of the ratios. Pure power law produces very low values of R for S$\;<$`<!-- -->`{=html}0.1. Values of R for S$\;<$`<!-- -->`{=html}0.1 increase with increase in the magnitude of the break at the He II ionization edge, however, the maximum value of R is around 0.1 for infinite break. For larger values of S, spectra of HM96 gives higher values of R compared to all other spectral shapes. This is due to the decrement in spectra of HM96 short wards of hydrogen ionization edge.

As the heavy element absorption systems are associated with galaxies, it is possible that the radiation field incident on the absorbers may get a significant contribution from the stellar sources. As the stellar radiation field has very few photons beyond 3 Ryd, high values of R are possible if the UV field is dominated by galactic contribution. We have constructed photoionization models for clouds irradiated by different proportions of galactic and UV background (HM96) fields. Steidel (1995), from his study of a large sample of galaxies associated with QSO absorption lines, finds these galaxies to be normal in the sense of their star formation rates. We have therefore taken the shape of the galactic field to be that given by Bruzal (1983). Pure galactic background can not reproduce the observed ratios as it produces much higher values of R, compared to the observed values. A contribution to the radiation field from the intergalactic AGN background is necessary. In Figure 2d we have plotted the column density ratios calculated for the case when the ratio of the flux due to galactic radiation is 90% of the total flux at 1 Ryd, the rest 10% being contributed by the AGN background (HM96). The observed values of R for all the intervening systems are lower than the results of this model. We can thus conclude that a minimum of 10 % (and possibly a much larger fraction, as will be seen below) of the background radiation incident on the intervening absorbers is contributed by the AGNs.

## Non-DLA intervening systems

Assuming the shape of the radiation field incident on the absorbers producing these lines to be that given by HM96, the observed ion ratios for two of the systems, with redshifts 2.1 and 2.77, are consistent with the results of the photoionization models. For other three systems, one with redshifts 2.1 and two with redshifts 2.7, the observed values of R, even after allowing for the errors are considerably smaller than the model predictions. Adding galactic radiation field to that of HM96 can only increase R, thereby increasing the discrepancy. For other shapes of the background radiation the observed ratios for four of the five systems are either consistent with or are lower than the model values. For three of these systems we can definitely rule out overabundance of Si w.r.t. C for any shape and intensity of the background radiation. For these systems we can also rule out any contribution of the galactic field to the radiation field. This also holds for the remaining two systems if the shape of the background is that given by HM96. Overabundance of Si by factors $>$ 1.5 is necessary only for one system for other shapes of the background. Alternatively a small contribution from the stellar sources is required.

## DLA intervening systems

For two of the systems, with redshifts 3.08 and 3.39, we can rule out the overabundance of Si for any shape of the background. Between the three shapes of the background radiation, the observed ion ratios for other 3 sytems with redshifts between 1.7 and 3.38, are consistent with the results of the photoionization models, without requiring any overabundance of Si or any contribution from the galactic sources. For the remaining two systems, with redshifts 2.84 and 3.39, the observed values of R are higher than model results for all the three shapes of the background. An overabundance of Si by factors $>$ 1.5 or a significant contribution (Fig.2d) from the galactic sources is needed.

## Associated systems

For the two associated systems with redshifts around 3.0, S is smaller than 0.1 and R is larger than 0.3. Associated systems being close to the QSO are expected to be ionized by pure power law radiation field. Such a field, produces (Fig. 2a) smaller values of R. The errors in the observed ratios are, however, too large and within those error bars the ratios may be consistent with the model results. However, though no overabundance is warranted by the observations due to the large uncertainty in the observed values, one can not rule out the possibility of Si being overabundant w.r.t. C by a large ($>$`<!-- -->`{=html}10) factor in these systems.

## Lyman alpha forest clouds

In Fig 3 we have plotted the observed ratios, R, vs. the ratio of column densities of C II to C IV in Lyman alpha forest lines observed by SC96. The theoretical results, assuming the UV background near the Lyman alpha absorbers is purely from AGNs, for four spectral shapes (i) power law with a slope of -1.0 (ii) power law with a slope of -1.5 (iii) power law with a slope of -1.5 with a cutoff at 4 Ryd and (iv) field of HM96, are also plotted for N$\rm_{H\;I}\;=\;10^{16}$ km s$^{-1}$. These results are independent of N$\rm_{H\;I}$ for N$\rm_{H\;I}\;\le\;10^{16}$ km s$^{-1}$. No overabundance of Si is required for half of the LAFCs which are mostly at redshifts smaller than 3.1 (SC96). The column densities for these systems are consistent with radiation field without any break at 4 Ryd. The amount of overabundance of Si needed for the other systems (mostly at redshifts lager than 3.1) does depend on the shape of the background and is smaller than a factor of 2 if a complete cutoff beyond 4 Ryd is assumed. We thus agree with the conclusion of SC(96) that the data do indicate a break in the background at 4 Ryd at high redshifts indicating a higher He II opacity at these redshifts.

# Summary

We have analyzed high resolution observations of Si and C absorption lines in the QSO spectra in order to study the ionization state of absorbers and its change with redshift as well as to understand the overabundance of Si w.r.t. C.

The column density ratios of Si IV to C IV in LAFCs, show a correlation with redshift, as already noted by SC96. The observations of this ratio in 30 non-DLA intervening absorbers, as well as in 23 intervening DLAS, however, fail to show any correlation with redshift. The data do not show any abrupt change in the ionic ratios at any particular redshift, unlike the case of LAFCs noted by SC96. Thus there is no evidence for the change in the opacity of the Universe beyond 4 Ryd from the column density ratios of the heavy element line systems in the QSO spectra. It may be argued that the radiation field incident on the heavy element absorbers gets a significant contribution from local stellar sources, thereby diluting the effect which is seen in the LAFCs. Our analysis of ion ratios, presented in the previous section, argues against such a possibility. Srianand and Khare (1995) have also presented several arguments against such a possibility.

Observed column density ratios of Si IV to C IV and of Si II to Si IV in several intervening and associated systems have been compared with the results of photoionization models with different shapes of the incident radiation. In spite of a large uncertainty in the observed values, definite conclusions can be drawn about the overabundance. We find that for all the non-DLA intervening systems the overabundance of Si can be ruled out if the shape of the background is as given by HM96. For other shapes of the background also, the overabundance can be ruled out for three of the five systems. The other two systems may allow an overabundance by factors $>$ 1.5 if no contribution from galactic sources to the background is assumed to be present. For two of the DLA systems also, overabundance can be ruled out for any shape of the background. Three other DLA systems are consistent with the observations for shape of the background given by HM96. The remaining two DLA systems, however, either require an overabundance by factors $>$ 1.5 or a significant contribution from the galactic radiation to the background. The possibility of overabundance by factors $>$ 10 can not be ruled out for the associated systems. Lyman alpha forest clouds at high ($>$`<!-- -->`{=html}3) redshifts do indicate an overabundance of Si over C as well as higher opacity of the Universe to radiation beyond He II ionization edge at these redshifts.

The authors are greatful to the Japan Society for Promotion of Science and the Department of Science and Technology (Government of India) for sponsoring the collaboration. PK thanks the members of the Theoretical Astrophysics group of the University of Osaka for warm hospitality and R. Srianand for discussion. This work was partially supported by a grant (No. SP/S2/013/93) by the Department of Science and Technology, Government of India.

# References

Bahcall, J. N. et al 1993, ApJ, 87, 1 Bajtlik, S., Duncan R.C., Ostriker, J. P. 1988, ApJ., 327, 57 Bechtold, J., Weyman, R. J., Lin, Z. , Malkan, M. A. 1987, ApJ, 315, 180 Bergeron, J. & Ikeuchi, S. 1990 A&A, 235, 8 Bruzal, G. 1983, ApJS, 53, 497 Cristiani, S., D'Odorico, S., Fontana, A., Giallongo, E. , Savaglio, S. 1995, MNRAS, 273, 1016 Davidsen, A. F., Kriss, G. A. , Zheng, W. 1996, Nature, 380, 47 Fan, X. , Tytler, D., 1994, ApJS, 94, 17 Giallongo, E., Cristiani, S., Fontana, A. , Trevese, D. 1993, ApJ, 416, 137 Giroux, M. L., Sutherland, R. S. , Shull, J. M. 1994, ApJ, 435, L101 Haardt, F. , Madau, P. 1996, ApJ, 461, 20 Isobe, Feigelson , Nelson, 1986 ApJ, 306, 490 Jackobsen, P., Boksenberg, A., Deharveng, J. M., Greenfield, P., Jedrzejewski, R. , Paresce, F. 1994, Nature, 370, 35 Khare, P. , Rana, N. C. 1993, Journal of Astron & Astrophys, 14, 83 Lauroesch, J.T., Truran, J.W., Welty, D.E. , York, D.G. 1996, PASP, 108, 641 Lu, L., Sargent, W.L.W., Barlow, T.A., Churchill C. W. , Vogt, S. S. 1996, ApJS, 107, 475 Madau, P. , Meiksin, A. 1994, ApJ, 433, L53 Madau, P. , Shull, J. M. 1996, ApJ, 457, 551 Petitjean, P. , Bergeron, J. 1994, A&A, 283, 759 Petitjean, P., Rauch, M. , Carswell, R. F. 1994, A&A, 291, 29 Pettini, M., King, D. L., Smith, L. J. , Hunstead, R. W. 1995, 'QSO Absorption Lines' Ed: G. Meylan, (Springer), P 71 Prochaska, J. X. , Wolfe, A. M. 1996, preprint Sargent, W. L. W., Boksenberg, A. , Steidel, C. C. 1988, ApJS, 68, 539 Savaglio, S., Cristiani, S., D'Odorico, S., Fontana, A., Giallongo, E. , Molaro, P. 1996, preprint Songaila, A. , Cowie, L. L. 1996, AJ, 112, 335 Srianand, R. Khare, P. 1995 ApJ, 444, 643 Steidel, C. C., 1995, in 'QSO Absorption Lines' Ed: G. Meylan, (Springer), 139 Steidel, C. C., Sargent, W. L. W. , Boksenberg, A. 1988, ApJ, 333, L5 Storrie-Lombardi, L. J., McMahon, M. J., Irwin, M. J. , Hazard, C. 1996, ApJ, 468, 121 Tripp, T. M., Lu, L. , Savage, B. 1996, ApJS, 102, 239 Wampler, E. J., 1991, ApJ, 368, 40 Wampler, E. J., Petitjean, P. , Bergeron, J. 1993, AA, 273, 15 Wolfe, A., Fan, X., Tytler, D., Vogt, S., Keane, M. J. , Lanzetta, K. M. 1994, ApJ, 435, 101

**Figure Captions**

Fig 1: Plot of the ratio of column densities of Si IV and C IV with redshift. Triangles and circles represent non-DLA and DLA intervening systems respectively. Upper limits are indicated by T.
Fig. 2: Theoretical and observed column density ratios of Si IV and C IV vs that of Si II and Si IV. The dotted line is for N$_{\rm
H\;I}\le10^{16}$ cm$^{-2}$, solid line is for N$_{\rm H\;I}=
3.0\times10^{17}$ cm$^{-2}$ and the dash-dotted line is for N$_{\rm
H\;I}= 10^{20}$ cm$^{-2}$. Circles correspond to the DLAS, squares correspond to the associated systems and triangles correspond to the non-DLA intervening systems. The shape of the UV continuum for Fig 2a, 2b and 2c is power law with slope =-1.5, Haardt & Madau (96) spectra for redshift of 2.5 and power law with slope =-1.5, with a cutoff at 4 Ryd, respectively. Fig 2d shows the results for a combined background field due to galaxies (90$\%$ at 1 Ryd) and that given by Haardt & Madau (96) at the redshift of 2.5 (10$\%$ at 1 Ryd). Horizontal dashed lines indicate the range of observed values.
Fig. 3: Theoretical and observed (SC96) column density ratios of Si IV and C IV vs that of C II and C IV for Lyman alpha forest lines for different shapes of UV background for N$_{\rm H\;I}\le10^{16}$ cm$^{-2}$. Solid line is for Haardt & Madau (96) spectra for redshift of 2.5; dashed line is for power law with slope of -1.0; dotted line is for power law with slope of -1.5 and dash-dotted line is for power law with slope of -1.5 with a cutoff at 4 Ryd.


author: F. Lignières ; P. Petit ; T. Böhm ; M. Aurière
date: Received March 6, 2009/ Accepted April 29, 2009
subtitle: Towards a new class of magnetic A-type stars[^1]
title: First evidence of a magnetic field on Vega

We report the detection of a magnetic field on Vega through spectropolarimetric observations. We acquired 257 Stokes V, high signal-to-noise and high-resolution echelle spectra during four consecutive nights with the NARVAL spectropolarimeter at the 2-m Telescope Bernard Lyot of Observatoire du Pic du Midi (France). A circularly polarized signal in line profiles is unambiguously detected after combining the contribution of about $1200$ spectral lines for each spectrum and summing the signal over the 257 spectra. Due to the low amplitude of the polarized signal, various tests have been performed to discard the possibility of a spurious polarized signal. They all point towards a stellar origin of the polarized signal. Interpreting this polarization as a Zeeman signature leads to a value of $-0.6 \pm 0.3$ G for the disk-averaged line-of-sight component of the surface magnetic field. This is the first strong evidence of a magnetic field in an A-type star which is not an Ap chemically peculiar star. Moreover, this longitudinal magnetic field is smaller by about two orders of magnitude than the longitudinal magnetic field (taken at its maximum phase) of the most weakly magnetic Ap stars. Magnetic fields similar to the Vega magnetic field could be present but still undetected in many other A-type stars.

# Introduction

Despite recent progress in stellar magnetic field measurements, spectropolarimetric surveys of early-type stars indicate that photospheric magnetic fields can only be detected in a small fraction of these stars. Without direct constraints on the magnetic field of the vast majority of early-type stars, our understanding of the role of magnetic fields on the structure and evolution of intermediate mass and massive stars is necessarily limited. In this Letter, we report the detection of a magnetic field on Vega and argue that Vega is probably the first member of a new class of yet undetected magnetic A-type stars.

The proportion of stars hosting a detectable magnetic field is more firmly established for main sequence stars of intermediate mass (late-B and A-type stars) than for massive stars (early B and O-type stars) or intermediate mass pre-main-sequence stars (Herbig Ae/Be stars). Magnetic A-type stars are indeed identified with the group of Ap-Bp chemically peculiar stars (excluding the subgroup of HgMn stars) since all known magnetic A-type stars belong to this group and, when observed with sufficient precision, Ap/Bp stars always show photospheric magnetic fields . The incidence of the Ap/Bp chemical peculiarity among A-type stars then leads to a $5-10\%$ estimate of magnetic stars . Note that magnetic field detections have been reported for a few Am and HgMn stars but remain debated because they could not be confirmed by further investigations . Thanks to new high-resolution spectropolarimeters, magnetic fields are now also detected in pre-main-sequence stars and in massive stars. According to recent surveys, the fraction of magnetic stars among Herbig Ae/Be stars is $7 \%$ , while the rate of detection for early B and O-type stars is also small .

The magnetic fields of Ap/Bp stars are characterized by a strong dipolar component, a long-term stability and dipolar strengths ranging from a lower limit of about 300 Gauss to tens of kilo-Gauss . Thus, if a population of weak dipolar-like fields corresponding to a weak field continuation of Ap/Bp stars exists, a longitudinal component of the magnetic field in the range of $10$ to $100$ Gauss should have been detected by recent spectropolarimetric surveys of non Ap/Bp stars . Instead, these surveys suggest there is a dichotomy between the population of strong, stable and dipolar-like magnetic fields corresponding to the Ap/Bp stars and the rest of A-type stars, whose magnetic properties remain unknown, except that their surface longitudinal magnetic field should be very small.

Vega is well suited for the search of magnetic fields among A-type non Ap/Bp stars. Its brightness and its low equatorial projected velocity ensure high signal-to-noise V spectra, while the number of spectral lines of an A0-type star is important enough to allow a very large multiplex gain by gathering the polarimetric signal of all the lines using a cross-correlation technique (Least-Squares Deconvolution, , LSD hereafter). Another advantage of Vega's brightness is that its fundamental parameters are well known relative to other more anonymous stars . In particular, spectral analysis and interferometric observations have shown that Vega is a rapidly rotating star seen nearly pole-on .

Vega was already included in a previous spectropolarimetric survey of A-type non Ap/Bp stars using NARVAL at the Telescope Bernard Lyot of Pic du Midi, but the analysis of its 11 Stokes V spectra was not conclusive. Here we present the results of a four night observing run fully dedicated to Vega, during which more than 300 Stokes V spectra were obtained. Summing the information over a large number of these spectra leads to an unambiguous detection of a polarized signal.

The observations are described and interpreted in the next section. In section 3, the origin of a $\sim \!1$ G longitudinal magnetic field in an A-type non Ap/Bp star is discussed and some of the perspectives opened by this field detection are considered. Our conclusions are given in section 4.

# Instrumental setup, data reduction and multi-line extraction of Zeeman signatures

The observing material was gathered at the Telescope Bernard Lyot (Observatoire du Pic du Midi, France) using the NARVAL spectropolarimeter. As a strict copy of ESPaDOnS , NARVAL spectra provide a simultaneous coverage of the whole optical domain (from 370 nm to 1,000 nm) at high spectral resolution ($R=$`<!-- -->`{=html}65,000). The instrument consists of a bench-mounted spectrograph and a Cassegrain-mounted polarimeter, with an optical fiber carrying the light between the two units. A series of 3 Fresnel rhombs (two half-wave rhombs that can rotate about the optical axis and one fixed quarter-wave rhomb) are used in the polarimeter, followed by a Wollaston prism which splits the incident light into two beams, respectively containing light linearly polarized perpendicular/parallel to the axis of the prism. Each Stokes V spectrum is obtained from a combination of four sub-exposures taken with the half-wave rhombs oriented at different azimuths . The data reduction is performed by Libre-Esprit, a dedicated, fully automated software described by .

The data were collected during 4 consecutive nights in the summer of 2008, from July 25 to July 28, using 6 sec integration times for each sub-exposure of the Stokes V sequences (except the first two sequences of the run, for which exposure times of 15 and 10 sec were adopted). We retained the 257 Stokes V spectra with a typical peak signal-to-noise ratio (S/N hereafter) of 1,500 per 1.8  km s$^{-1}$, around $\lambda = 600$ nm.

For each spectrum, both Stokes I & V parameters were processed using the LSD cross-correlation method . Using a line mask computed from a stellar atmospheric model with T$_{\rm eff}$=10,000 K and $\log g$ = 4.0 , we calculated LSD line profiles from a total of 1,200 photospheric lines. The multiplex gain in the S/N from the raw spectra to the LSD mean profiles is about 30, reducing the noise level of the cross-correlation profiles to between $\sigma = 3$ and $7\times 10^{-5}I_{c}$, where $\sigma$ is the standard deviation of the noise and $I_{c}$ stands for the intensity of continuum. Since no signature was observed above noise level in individual Stokes V LSD profiles, we then calculated an average of the 257 profiles, where each profile is weighted by the square of its S/N. In this global profile, the noise in the Stokes V parameter is further decreased to $\sigma = 2\times 10^{-6}I_{c}$ (Fig. ) and a signature is now observed in circular polarization with an amplitude of $10^{-5}I_{c}$ (that is $5$ times the noise level). Running a $\chi^2$ test on the signature , we found a reduced $\chi^2$ of 3.5, which corresponds to a false-alarm probability of $3\times 10^{-11}$.

Various tests have been performed to make sure that the polarization is of stellar origin and not due to an artefact of the instrument or the reduction process. This is particularly important in the present case, as the amplitude of the polarized signal is the lowest detected by NARVAL to date. A strong test to discard the possibility of a spurious signal is the "null" profile calculated from a different combination of the four sub-exposures constituting the polarimetric sequence . As shown in Fig. , no detectable counterpart of the Stokes V signal is seen in the "null" profile (note that a similar conclusion is reached by calculating another null profile (not shown here) from another possible combination of the sub-exposures). We then checked that the signal possesses the expected properties of a stellar polarized signal. First, we split the whole time-series into two independent subsets, containing respectively the first and second half of the observing run, both subsets having equivalent signal-to-noise ratios. As can be seen in Fig. a, the polarized signal is present in both sets, the false-alarm probabilities based on the $\chi^2$ test being $5\times 10^{-6}$ and $6\times 10^{-3}$, respectively. Second, we built two line-lists from the atmospheric model, containing all spectral lines with Landé factors respectively higher and lower than $g_c = 1.2$. The Stokes V profiles computed from the two line-lists are plotted in Fig. b. As expected, the amplitude of the polarized signal appears higher for the high Landé factor lines than for the low Landé factor lines. The peak-to-peak amplitudes of the polarized signals taken inside the line profile are respectively $2.8
\times 10^{-5}I_c$ and $2.0\times 10^{-5}I_c$. Their difference slightly exceeds the noise level ($\sigma = 6\times 10^{-6}I_c$ and $5\times 10^{-6}I_c$, respectively) and their ratio is roughly consistent with the ratio of the average Landé factors of the line-lists ($g_m=1.51$ and $g_m=0.94$ respectively), taking into account that the corresponding Stokes I LSD profiles have a similar depth, within 10%. Third, we checked that the signal was still consistently recovered when other ways of splitting our line-list were considered (low versus high excitation potential or low versus high wavelengths). Finally, we tested the effect of changing the line mask. Indeed, spectroscopic and interferometric studies of Vega have shown that its surface temperature is inhomogeneous, due to the gravity darkening effect induced by its rapid rotation. According to the latest model based on interferometric results , the equatorial velocity of Vega is 274 km s$^{-1}$ and the effective temperature and the gravity decrease from 9988 K and $\log g = 4.07$ at the pole to 7600 K and $\log g = 3.5$ at the equator. There is a significant discrepancy with the spectroscopic analysis of , who find that the polar to equator temperature difference is only $\sim\!900$ K while the equatorial velocity is reduced to 175 km s$^{-1}$. Taking the extreme case of a T$_{\rm eff}$=7500 K and $\log g$ = 3.5 stellar atmosphere, we computed the LSD line profiles (not shown here) for the associated line mask and this time obtained a false-alarm probability of $10^{-1}$. The polarized signal is consistent with the one derived from our previous (hotter) atmospheric model but the low significance of the detection is due to the inadequacy of the line mask that results in a much higher noise level.

These tests strongly support that the signal is of stellar origin and therefore that Vega possesses a magnetic field.

The circularly polarized signal has the typical anti-symmetric shape of a Zeeman signature (Fig. ). However, as compared to the width of the Stokes I line profile, it only shows up within a limited range of radial velocities about the line-center. This suggests that the magnetic field distribution is axisymmetric and confined in the polar region. However, a more detailed analysis will be needed to specify the surface field distribution of Vega. First, as the 257 spectra at our disposal cover a range of rotation phases, Zeeman signatures from non-axisymmetric magnetic features, if any, are mostly averaged out from the time-averaged line profile. Second, due to Vega's temperature inhomogeneities, the weak line profiles range from flat-bottomed to "V" shapes . As the LSD profile is obtained by assuming that all lines have a common profile, its interpretation in terms of the surface field distribution is not straightforward in the present context.

We use the center-of-gravity method to estimate the longitudinal magnetic field $B_l$ : $$B_l = -2.14\times10^{11}\frac{\int vV(v)dv}{\lambda_0 g_m c\int(I_c - I(v))dv}$$ where $v$ (km s$^{-1}$) is the radial velocity, $\lambda_0$ (nm) the mean wavelength of the line-list used to compute the LSD profiles, $g_m$ the mean Landé factor and $c$ (km s$^{-1}$) the light velocity. The integration limits cover a $\pm30$ km s$^{-1}$ velocity range around the line centroid. Using this equation, we obtain $B_l = -0.6 \pm 0.3$ G.

# Discussion

Three basic features distinguish the present detection from previous measurements of magnetic fields in main-sequence stars of intermediate mass : (i) It is the first time a magnetic field is detected in an A-type star which is not an Ap/Bp chemically peculiar star (if we exclude the debated field detections in a few Am and HgMn stars discussed in ). (ii) The longitudinal magnetic field of Vega is smaller by about two orders of magnitude than the field of the most weakly magnetic Ap/Bp stars. Indeed, the longitudinal field of a 300 G dipolar field aligned with the stellar rotation axis and viewed pole-on is close to 100 G, that is about two orders of magnitude larger than the $0.6$ G field of Vega. The longitudinal component of a dipolar field actually depends on its angle with respect to the rotation axis. But, whatever this angle, the amplitude of the circular polarization in the LSD Stokes V profile of a 300 G dipolar field will be more than one order of magnitude larger than that of Vega. (iii) The LSD Stokes V profile of Vega is also qualitatively distinct from LSD Stokes V profiles of Ap/Bp stars since the polarized signal of Vega is concentrated in the weakly Doppler shifted regions of the projected stellar disk.

These marked observational differences between Vega and the Ap/Bp magnetic stars suggest that we should consider Vega as a new type of magnetic A-type star. As there is no reason to believe Vega is unique among A-type stars, Vega should be considered as the first member of a new class of magnetic A-type stars.

The existence of such a new class might help understand some otherwise puzzling observations of the pre-main-sequence and post-main-sequence intermediate mass stars. The Herbig Ae/Be stars show a strong activity which has led investigators to suspect a widespread presence of magnetic fields in these stars . Nevertheless, these magnetic fields have not been found, since only a small fraction of Herbig Ae/Be stars appears to host one . Note that a similar discrepancy between widespread activity and a small fraction of detected fields exists in OB stars . A new class of magnetic A-type stars would shed new light on this issue. Indeed, the progenitors of these magnetic A-type stars could be the Herbig Ae/Be stars where magnetic fields have not been detected yet, this non-detections being compatible with the fact that magnetic fields of the same intensity are much more difficult to detect in the faint Herbig Ae/Be stars than in a bright A-type star like Vega. On the post-main-sequence side, the study of the white dwarf magnetic fields suggests that Ap/Bp stars are not sufficient to be the progenitors of magnetic white dwarfs . A new class of magnetic A-type stars might also help to resolve this issue.

The consequences for Vega itself should also be considered. Its magnetic field could indeed trigger active phenomena in its atmosphere. Signs of spectroscopic variability have been reported, but have not been confirmed since . On the other hand, despite its status as a photometric standard, a photometric variability of $1-2 \%$ with occasional excursions to $4\%$ has been reported . This might be produced by photospheric temperature inhomogeneities induced by its magnetic field. However, because of the near pole-on configuration of Vega, the variability would rather be due to intrinsic changes of the magnetic field than to rotational modulation.

The origin of Vega's magnetic field could be attributed to one of the three mechanisms generally invoked for early-type stars, namely (i) the fossil field hypothesis, (ii) the envelope dynamo, (iii) the convective core dynamo. Let us first consider the fossil field hypothesis, whereby the ISM magnetic field is confined and amplified during stellar formation. It is regarded as the most consistent explanation of the magnetic fields observed in Ap/Bp stars , but, as proposed by , it could also account for another population of stars hosting weak longitudinal magnetic fields. Their argument is based on the fact that large-scale, organized magnetic field configurations are subjected to a pinch-type instability driven by differential rotation when the magnetic field drops below a critical value. Consequently, for a distribution of large scale organized fields of different strengths issued from the star formation process, the instability would produce a magnetic dichotomy between a population of strong and stable large scale fields like in Ap/Bp stars and another population where the destabilized configuration is now structured at small length scales, thus resulting in a weak longitudinal field. A simple estimate of the critical field has been found to be consistent with the reported lower limit of Ap/Bp stars. Here, both the detection of a very small longitudinal field in Vega and the gap between this field and the lowest magnetic fields of Ap/Bp stars reinforce this scenario. Nevertheless, this scenario is not complete as it does not say what happens to the destabilized field configuration, which could either decay or be regenerated by a dynamo.

The magnetic field of Vega could indeed be generated by an envelope dynamo where the energy source is the rotation of the star. Following , the dynamo loop initiated by the differential rotation could be closed by the pinch-type instability mentioned above. This interesting possibility has been investigated by numerical simulations in a simplified cylindrical configuration and in a solar context , leading to opposite outcomes. Simulations in more realistic conditions for A-type stars are clearly needed to test this envelope dynamo. An important issue concerns the origin of the envelope differential rotation, which is a basic ingredient of this dynamo but which is not forced by a strong stellar wind in A-type stars, contrary to what is expected to occur in OB and Herbig Ae/Be stars . The third possibility is a dynamo in the convection core. While magnetic fields are likely to be generated there, an efficient mechanism to transport it throughout the radiative envelope to the star surface has not yet clearly been identified . We note that in the three cases considered, the magnetic field is expected to be structured at small scales and also probably variable in time. This calls for a spectropolarimetric monitoring of Vega that will investigate the surface distribution and the temporal variation of its magnetic field.

# Conclusion

A circularly polarized signal has been detected by accumulating a large number of high-quality echelle spectra of Vega with the NARVAL spectropolarimeter. The data analysis strongly supports a stellar origin of the polarization and thus the presence of a magnetic field on Vega. Due to the unprecedented low level of the detected polarization, new independent measurements will still be important to confirm this result. A magnetic field on Vega suggests that other A-type stars which are not Ap/Bp stars host weak magnetic fields and that their study can shed a new light on early-type star magnetism. While a spectropolarimetric survey of bright A-type stars will be necessary to find these stars, a detailed investigation of Vega's magnetic field should also provide clues to the origin of this magnetism.

[^1]: Based on observations at Telescope Bernard Lyot of Observatoire du Pic du Midi, CNRS/INSU and Université de Toulouse, France


author: Donald Arseneau [^1]
date: 2013-09-16
title: `url.sty` version 3.4

The package defines a form of command that allows linebreaks at certain characters or combinations of characters, accepts reconfiguration, and can usually be used in the argument to another command. It is intended for formatting email addresses, hypertext links, directories/paths, etc., which normally have no spaces. The font used may be selected using the command, and new url-like commands may be defined using . This package does not make hyper-links! For that purpose, see the hyperref package (or some other deprecated ones).

| Usage | Conditions |
|:---|:---|
| `\url{ }` | The argument must not contain unbalanced braces. If used in the argument to another command, the  argument cannot contain any "`%`",+", "`#`", or "`^^`+", or end with "`\`+". |
| `\url|  |` | where "`|`" is any character not used in the argument and not "`{`" or a space. The same restrictions apply as above except that the argument may contain unbalanced braces. |
|  | for the defined-url ""; such a command can be used anywhere, no matter what characters it contains. |

The "" command is fragile, and its argument is likely to be very fragile, but a defined-url is robust.

# Package options

Package Option: `obeyspaces`

Ordinarily, all spaces are ignored in the url-text. The "`[obeyspaces]`" option allows spaces, but may introduce spurious spaces when a url containing "" characters is given in the argument to another command. So if you need to obey spaces you can say "`[obeyspaces]``'{url``'}`", and if you need both spaces and backslashes, use a defined-url.

Package Option: `hyphens`

Ordinarily, breaks are not allowed after "`-`" characters because this leads to confusion. (Is the "`-`" part of the address or just a hyphen?) The package option "`[hyphens]`" allows breaks after explicit hyphen characters. The command will **never ever** hyphenate words.

Package Option: `spaces`

Likewise, given the "`[obeyspaces]`" option, breaks are not usually allowed after the spaces, but if you give the options "`[obeyspaces,spaces]`", will allow breaks at those spaces.

> Note that it seems logical to allow the sole option "`[spaces]`" to let input spaces indicate break points, but not to display them in the output. This would be easy to implement, but is left out to avoid(?) confusion.

Package Option: `lowtilde`

Normal treatment of the `~` character is to use the font's "" character, if it has one (or claims to). Otherwise, the character is faked using a mathematical "". The "`[lowtilde]`" option causes a faked character to be used always (and a bit lower than usual).

Package Option: `allowmove`

This option suppresses the test for being used in a so-called moving argument (check "fragile command"). Using it will enable to function in more contexts, but when it does fail, the error message may be incomprehensible.

# Defining a defined-url

Take for example the email address "firstname.lastname@example.com" which could not be given (using "" or "") in a caption or parbox due to the percent sign. This address can be predefined with

> `\urldef{\myself}\email@example.com}` or
> `\urldef{\myself}\firstname.lastname@example.com|`

and then you may use "" instead of "`\firstname.lastname@example.com}`" in an argument, and even in a moving argument like a caption because a defined-url is robust.

# Style

You can switch the style of printing using "`\urlstyle{`$xx$`}`", where "$xx$" can be any defined style. The pre-defined styles are "`tt`", "`rm`", "`sf`" and "`same`" which all allow the same linebreaks but use different fonts — the first three select a specific font and the "`same`" style uses the current text font. You can define your own styles with different fonts and/or line-breaking by following the explanations below. The "" command follows whatever the currently-set style dictates.

# Alternate commands

It may be desireable to have different things treated differently, each in a predefined style; e.g., if you want directory paths to always be in typewriter and email addresses to be roman, then you would define new url-like commands as follows:

> `\DeclareUrlCommand``{``}`
> `\DeclareUrlCommand\email{\urlstyle{rm}}`
> `\DeclareUrlCommand\directory{\urlstyle{tt}}`.

In fact, this example is exactly the definition which might be pre-defined by the package. Furthermore, basic is defined with

> `\DeclareUrlCommand\url{}`,

without any *settings*, so it uses whatever and other settings are already in effect.

You can make a defined-url for these other styles, using the usual command as in this example:

> `\urldef{\myself}{\firstname.lastname@example.com}`

which makes act like `\firstname.lastname@example.com}`, if the command is defined as above. The command would then be robust.

# Defining styles

Before describing how to customize the printing style, it is best to mention something about the unusual implementation of . Although the material is textual in nature, and the font specification required is a text-font command, the text is actually typeset in *math* mode. This allows the context-sensitive linebreaking, but also accounts for the default behavior of ignoring spaces. (Maybe that underlying design will eventually change.) Now on to defining styles.

To change the font or the list of characters that allow linebreaks, you could redefine the commands , , , etc., directly in the document, but it is better to define a new 'url-style' (following the example of and ) which defines all of , , , , and .

## Changing font

The command selects the font. The definition of done by the pre-defined styles varies to cope with a variety of LaTeX font selection schemes, but it could be as simple as `\def\UrlFont{\tt}`. Depending on the font selected, some characters may need to be defined in the list because many fonts don't contain all the standard input characters.

## Changing linebreaks

The list of characters after which line-breaks are permitted is given by the two commands (list macros) and . They consist of repeating for each relevant character `c`.

The differences are that 'BigBreaks' typically have a lower penalty (more easily chosen) and do not break within a repeating sequence (e.g., "`DEC::NODE`"). (For gurus: 'BigBreaks' are treated as mathrels while 'Breaks' are mathbins; see *The TeXbook*, p. 170.) The result is that a series of consecutive 'BigBreak' characters will break at the end and only at the end; a series of 'Break' characters will break after the first and after every following *pair*; there will be no break between a 'Break' character and a following 'BigBreak' char; breaks are permitted when a 'BigBreak' character is followed by 'Break' or any other char. In the case of `http://` it doesn't matter whether `:` is a 'Break' or 'BigBreak' — the breaks are the same in either case; but for (now ancient) *DECnet* addresses using `::` it was important to prevent breaks *between* the colons, and that is why colons are 'BigBreaks'. (The only other 'BigBreak' character is, optionally, the hyphen; slashes are regular 'Break's.)

It is possible for characters to prevent breaks after the next following character (this is used for parentheses). Specify these in .

You can allow some spacing around the breakable characters by assigning

> `\Urlmuskip = 0mu plus 1mu`

(with `mu` units because of math mode). You can change the penalties used for BigBreaks and Breaks by assigning

> `\mathchardef\UrlBreakPenalty=100`
> `\mathchardef\UrlBigBreakPenalty=100`

The default penalties are and . These have such odd non-LaTeX syntax because I don't expect people to need to change them often. (The `\mathchardef` does not relate to math mode; it is only a way to store a number without consuming registers.)

### Arbitrary character actions

You can do arbitrarily complex things with characters by specifying their definition(s) in . This makes them 'active' in math mode (mathcode `"8000`). The format for setting each special character `c` is: `\do\c{``}`, but other definitions not following this style can also be included.

Here is an example to make "`!`" inside force a line break instead of being treated verbatim (it uses LaTeX's ):

> `\makeatletter \g@addto@macro\UrlSpecials{\do\!{\newline}}`

Here is another overly-complicated example to put extra flexible muglue around each "`/`" character, except when followed by another "`/`", as in "`http://`", where extra spacing looks poor.

> % what we'll insert before and after each (lone) slash:
>     \newmuskip\Urlslashmuskip
>     \Urlslashmuskip=2mu plus2mu minus2mu
>
>     % change what / does:
>     \g@addto@macro\UrlSpecials{\do\/{\Urlspaceyslash}}
>
>     % need to look ahead:
>     \def\Urlspaceyslash{\futurelet\Urlssnext\finishUrlspaceyslash}
>
>     \def\finishUrlspaceyslash{%
>       \mskip\Urlslashmuskip  % extra space before
>       \mathchar8239  % "202f, i.e., binary op, \fam, / char
>       % if we see //, eliminate the extra space to taste:
>       \ifx\Urlssnext/\mskip-\Urlslashmuskip
>       \else\mskip\Urlslashmuskip \fi
>      }

If this sounds confusing … well, it is! But I hope you won't need to redefine breakpoints — the default assignments seem to work well for a wide variety of applications. If you do need to make changes, you can test for breakpoints using regular math mode and the characters "`+=(a`".

# Yet more flexibility

You can also customize the presentation of verbatim text by defining and/or . An example for ISO formatting of urls surrounded by `<  >` is

> \DeclareUrlCommand\url{\def\UrlLeft{<url:\ }\def\UrlRight{>}%
>         \urlstyle{tt}}

The meanings of and are *not* reproduced verbatim. This lets you use formatting commands there, but you must be careful not to use TeX's special characters (`\^_%~#$&{}` etc.) improperly. You can also define to reprocess the verbatim text, but the format of the definition is special:

> `\def\UrlLeft#1\UrlRight{` … do things with `#1` … `}`

Yes, that is `#1` followed by then the definition (a TeX macro with delimited arguments). For example, to produce a hyperTeX hypertext link:

> `\def\UrlLeft#1\UrlRight{%`
> `      \special{html:<a href="#1">}#1\special{html:</a>}}`

Using this technique, can provide a convenient interface for performing various operations on verbatim text. You don't even need to print out the argument! For greatest efficiency in such obscure applications, you can define a null url-style where all the lists like are empty.

Please note that this method is *not* how the hyperref package manages urls for its command, even though it makes use of . Instead, hyperref's reads its argument in a less-verbatim manner than described above, produces its hyperlink, and invokes to format the text. is the command descibed herein.

[^1]: Thanks to Robin Fairbairns for documentation conversion!


abstract: We introduce new algebro-topological invariants of directed networks, based on the topological construction of the directed clique complex. The shape of the underlying directed graph is encoded in a way that can be studied mathematically to obtain network invariants such as the Euler characteristic and the Betti numbers.
  .
  Two different cases illustrate the application of the Euler characteristic. We investigate how the evolution of a Boolean recurrent artificial neural network is influenced by its topology in a dynamics involving pruning and strengthening of the connections, and to show that the topological features of the directed clique complex influence the dynamical evolution of the network. The second application considers the directed clique complex in a broader framework, to define an invariant of directed networks, the network degree invariant, which is constructed by computing the topological invariant on a sequence of sub-networks filtered by the minimum in- or out-degree of the nodes.
  .
  The application of the Euler characteristic presented here can be extended to any directed network and provides a new method for the assessment of specific functional features associated with the network topology.
address: Neuroheuristic Research Group, Faculty of Business and Economics (HEC), University of Lausanne, Quartier Dorigny, 1015 Lausanne, Switzerland
author: Paolo Masulli and Alessandro E. P. Villa
title: The topology of the directed clique complex as a network invariant

# Background

The main interest of algebraic topology is to study and understand the functional properties of spatial structures. Algebro-topological constructions have been applied successfully in the field of data science with the application of the framework of persistent homology, which has proved to be a powerful tool to understand the inner structure of a data set by representing it as a sequence of topological spaces. A network is a set of points satisfying precise properties of connectedness, which can be used to define a class of topological spaces. Network theory aims to understand and describe the shape and the structure of networks, and the application of the tools developed within the framework of algebraic topology can provide new insights of network properties in several research fields.

The directed clique complex is a rigorous way to encode the topological features of a network in the mathematical framework of a simplicial complex, allowing the construction of a class of invariants which have only been recently applied for the first time in the context of network theory . Active nodes are those nodes whose state depend on a set of precise rules that depend on network topology and dynamics. In a highly interconnected network of such nodes, the activity of each node is necessarily related to the combined activity of the afferent nodes transmitted by the connecting edges. Due to the presence of reciprocal connections between certain nodes, re-entrant activity occurs within such network. Hence, selected pathways through the network may emerge because of dynamical processes that may produce activity-dependent connection pruning. The overall goal of these studies is to understand the properties of a network given the topology described by its link structure.

Neuronal networks are a complex system characterized by coupled nonlinear dynamics. This topic is a long-standing scientific program in mathematics and physics . In general, the synchronization of two systems means that their time evolution is periodic, with the same period and, perhaps, the same phase . This notion of synchronization is not sufficient in a context where the systems are excited by non-periodic signals, representing their complex environment. Synchronization of chaotic systems has been discovered and since then it has become an important research topic in mathematics , physics and engineering . In interconnected cell assemblies embedded in a recurrent neural network, some ordered sequences of intervals within spike trains of individual neurons, and across spike trains recorded from different neurons, will recur whenever an identical stimulus is presented. Such recurring, ordered, and precise interspike interval relationships are referred to as "preferred firing sequences". One such example can be represented by brain circuits shaped by developmental and learning processes . The application of tools from algebraic topology to the study of these systems and networks will be of great use for determining deterministic chaotic behavior in experimental data and develop biologically relevant neural network models that do not wipe out temporal information

In the current study we introduce a mathematical object, called directed clique complex, encoding the link structure of networks in which the edges (or links) have a given orientation. This object is a simplicial complex that can be studied with the techniques of algebraic topology to obtain invariants such as the Euler characteristic and the Betti numbers. We propose general constructions valid for any directed network, but we present an application to evolvable Boolean recurrent neural networks with convergent/divergent layered structure with an embedded dynamics of synaptic plasticity. The Euler characteristic, which is defined given the network connectivity, is computed during the network evolution. We show evidence that this topological invariant predicts how the network is going to evolve under the effect of the pruning dynamics. Despite being just a toy-example of the dynamics observed in biological neuronal networks, we suggest that algebraic topology can be used to investigate the properties of more refined biologically-inspired models and their temporal patterns. We show also that, for a directed network, the Euler characteristic computed on a sequence of networks generated by filtrating its nodes by in- and out-degrees can provide a general metric helpful for a network classification. Hence, the topological invariants computed for each network in the filtration give a sequence of numbers that may be interpreted as a fingerprint of the complete network.

## Acknowledgements

This work was partially supported by the Swiss National Science Foundation grant CR13I1-138032. We wish to thank Prof. Kathryn Hess, Prof. Ran Levi and Paweł Dłotko for suggesting the idea of using oriented cliques in order to define simplices in a directed graph, as it was shown in the talk given at the SIAM DS15 Conference in Snowbird (USA) in May 2015.

# Results and discussion

## Dynamics of artificial neural networks

We considered a directed graph representing a simplified model of feedforward neural network with convergent/divergent layered structure with few embedded recurrent connections. In this model, the nodes represent individual neurons and the connections between them are oriented edges with a weight given by the connection strength. We have computed the Euler characteristic and its variation during the evolution of such networks, both for the entirety of the nodes in the network and for the sub-network induced by the nodes that are active at each time step in order to detect how the structure changes as the network evolves. The Betti numbers and their variation during the network evolution were also computed but we do not discuss further this topological measurement. Notice that activation of the networks follows a very simple dynamics. The nodes of the input layer are activated at regular time intervals, which is not meant to be biologically realistic, but has been adopted to favor the simplicity of the model. It was shown elsewhere that a stable activity level in a network like this could be achieved only with an appropriate balance of excitatory and inhibitory connections. The networks studied here are oversimplified and formed only by excitatory nodes. We selected the ranges of the parameters such that the simulations maintained a level of activity for $100$ steps without neither saturation nor extinction of the activity, thus suggesting that connection pruning enabled topological changes. However, notice that even within selected areas of the parameter space of the simulations we observed that the activity level tended either to increase towards paroxysmal activation (i.e., saturation) or to decrease towards complete inactivation (i.e, extinction).

We observed that the Euler characteristic of the entire network could detect the pruning activity during the neural network evolution (Figure A). In particular, the step to step variation of the Euler characteristic matched the number of connections pruned over time. If we considered only the sub-network of the active nodes, we observed that the Euler characteristic decreased or increased if the number of active nodes increased or decreased, respectively (Figure B). Then, the Euler characteristic is a good estimator of the activity level within the network. These results confirmed that the Euler characteristic gives a precise measure of the topological changes in a network associated with connection pruning for the entire network and associated with activation patterns for the active sub-network.

The type of dynamics undergoing the neural network evolution and the structure of the directed clique complex of that network at the very beginning of the simulation (i.e. before the occurrence of connection pruning) were correlated. This was possibly the most unexpected and significant result regarding the dynamics of artificial neural networks. In the simulations leading to the activation of at least $5\%$ of the nodes, the average number of active units was correlated to the number of simplices, in the directed clique complex, of dimension two (Pearson correlation coefficient $r_{(370)} = .560$, $p < 0.001$) and dimension three ($r_{(370)} = 0.445$,$p < 0.001$). This may appear surprising because the topology of the directed clique complex of a network *a priori* ignores any dynamics of pruning, the evolution of the network topology and how this is going to influence the activation level. However, the rationale is that directed cliques are fully connected sub-networks, i.e. sub-networks with an initial and a final node that are connected in the highest possible number of ways. Then, a high number of directed cliques leads to a higher chance of propagation of the activation through the network. Notice the fact that it is essential to consider here only *directed* cliques, because the activation of a node occurs only if the connected upstream nodes are activated. Activation is indeed a phenomenon happening in a directional way prescribed by the connectivity pattern. The invariant presented should also be considered as a complementary measurement of complexity for the assessment of the computational power of Boolean recurrent neural networks .

## Network filtrations and invariants

The in- and out-degrees of nodes are important factors in shaping the network topology. We applied our topological construction to devise invariants for any directed networks. We compute the Euler characteristic on a sequence of sub-networks defined by *directed degeneracy* of their nodes, or in other words the in- and out-degrees of the vertices, as described in detail in the Methods section. Two separate sequences are defined because in- and out-degrees represent different aspects in the network connectivity. The sequences of sub-networks of a given network are a *filtration* in the sense that each network appearing in the sequence is contained in all those that follow. The values of the Euler characteristic for each network of the sequences gives rise to two separate sequences of integers that give a measure of the shape and the topology of the complete network. We propose this invariant to describe general directed networks.

The sequences of the Euler characteristic of the in-degree and the out-degree filtrations are plotted as a function of the normalized minimum degree of vertices for representative types of a scale-free network (SF) , a random network (RN) , and a small world network (SW) (Figure ). This normalization is necessary to compare networks of different sizes at each filtration level posing the maximum degree of the vertices in the network to 1, as described in the Methods section. Each network type was simulated 50 times using different random seeds. In the case of a SF network the values of the Euler characteristic of the in-degree filtration (dotted red line) is always larger than the curve of the out-degree filtration (solid blue line, Figure A). Moreover, for SF networks the curve of the in-degree increases sharply at near $0.4$ and reaches a maximum value of the Euler characteristic approximately at $0.6$ of the normalized maximum value of the vertex minimum degree. The out-degree curve increases monotonically after this level of vertex minimum degree but does not reach the in-degree curve. In the case of a RN network the curves of in- and out-degrees overlap at all levels of the filtration (Figure B). It is interesting to notice that the maximum value of the Euler characteristic is observed for the smallest values of the vertex minimum degree. Then, for RN networks, both curves decrease to a minimum at approximately $0.8$ of the normalized maximum value of the vertex minimum degree, followed by a monotonic increase. For a SW network both curves of in- and out-degrees start from the minimal value of the Euler characteristic with the least vertex minimum degree, followed by a non monotonic increase and a tendency of overlap between the two curves (Figure C). The monotonicity of the curves and the differences between in- and out-degree filtration differ greatly for the three types of networks, thus suggesting that this invariant is a good descriptor of network topology.

A distinct topological invariant defined for non-directed networks, referred to the Betti curves, was recently proposed by following the idea of filtering the network by the weight of connections. This invariant appears well suited for continuously distributed connection weights, for instance when the weights are related to the distances of points and represent a symmetric relation between nodes. In the case of directed networks with modifiable values of connection weights restricted to a limited set , the network dynamics evolves towards a bimodal distribution of the connection weights densely grouped near the minimum and maximum values of the range. This is a general behaviour in neuronal networks . In this kind of networks, filtering the network by the connection weights following is not suitable, because most connections would have the same weight. Our approach for directed networks is to filter the connections by the in- and out-degrees separately in order to measure how the nodes of each degree shape the topology of the network. It is important to point out that other methods are based on spectral properties of the adjacency matrix and therefore only make sense if all the transformations of the network data are linear .

The results presented here open the way to further applications of the topological invariants. The analytic study of the values of the Euler characteristic in the filtrations framework can provide a metric of similarity between networks which is only dependent on their internal topology, thus allowing the application of clustering algorithms for the detection of distinct functional classes of networks. The study of brain complex networks in clinical neuroscience offers as a particularly promising field of application of the new topological invariant, as suggested by other studies using different techniques to the same aim . Another promising application is the study of the temporal dynamics in neural activity. The finding of precise and repeating firing sequences in experimental and simulated spike train recordings has been discussed with respect to the existence of synfire chains or chaotic attractors . In both cases the underlying network structure is assumed to be a directed graph. This hypothesis together with the assumption of spike-timing modifiable connections provide a rational basis for the application of topological invariants towards understanding the association between topological structures and neural coding.

# Conclusions

We have developed new invariants for directed networks using techniques derived from algebraic topology, showing that this subject provides a very useful set of tools for understanding networks and their functional and dynamical properties. Simple invariants such as the Euler characteristic can already detect the changes in the network topology. The promising results shown here are a contribution to the application of algebraic topology to the study of more complex networks and their dynamics, including models of neuronal networks that are biologically inspired. We believe that the framework present here may open the way to many computational applications to unveil data structures in data and network sciences.

# Methods

## Graphs and clique complexes

An *abstract oriented simplicial complex* $K$ is the data of a set $K_0$ of vertices and sets $K_n$ of lists $\sigma = (x_0, \dots, x_n)$ of elements of $K_0$ (called *$n$-simplices*), for $n \geq 1$, with the property that, if $\sigma = (x_0, \dots, x_n)$ belongs to $K_n$, then any sublist $(x_{i_0}, \dots, x_{i_k})$ of $\sigma$ belongs to $K_k$. The sublists of $\sigma$ are called *faces*.

We consider a finite directed weighted graph $G = (V,E)$ with vertex set $V$ and edge set $E$ with no self-loops and no double edges, and denote with $N$ the cardinality of $V$. Associated to $G$, we can construct its *(directed) clique complex* $K(G)$, which is the directed simplicial complex given by $K(G)_0 = V$ and $$\label{eq:directed_clique_complex}
K(G)_n = \{(v_0, \dots, v_n) \colon (v_i, v_j) \in E \textrm{ for all } i < j \} \quad \textrm{ for } n \geq1.$$ In other words, an $n$-simplex contained in $K(G)_n$ is a directed $(n+1)$-clique or a completely connected directed subgraph with $n+1$ vertices. Notice that an $n$-simplex is though of as an object of dimension $n$ and consists of $n+1$ vertices.

By definition, a directed clique (or a simplex in our complex) is a fully-connected directed sub-network (Figure ) this means that the nodes are ordered and there is one source and one sink in the sub-network, and the presence of the directed clique in the network means that the former is connected to the latter in all the possible ways within the sub-network.

## The topological invariants

The directed clique complex is the basic topological object that allows us to introduce invariants of the graph: the *Euler characteristic* of the directed clique complex $K(G)$ of $G$ is the integer defined by $$\chi(K(G)) = \sum_{n=0}^N (-1)^n \ \vert K(G)_n \vert,$$ or in other words the alternating sum of the number of simplices that are present in each dimension.

Let us now consider, for each $n$, the vector space $\mathbf{Z}/2\langle K(G)_n\rangle$ given by the linear combinations of $n$-simplices with coefficients in the field of two elements $\mathbf{Z}/2$. We can define the *boundary maps* $\partial_n \colon \mathbf{Z}/2\langle K(G)_n\rangle \to \mathbf{Z}/2\langle K(G)_{n-1}\rangle$ which are given by mapping each simplex to the sum of its faces. Then we can define the quantities: $$\beta_n(K(G)) = \mathrm{dim}(\mathrm{ker}\ \partial_n) - \mathrm{dim}(\mathrm{Im}\ \partial_{n+1}),$$ given by the difference of the dimension of the space of the $n$-simplices whose boundary is zero and the dimension of the space of boundaries of $(n+1)$-simplices. It can be checked that, if we apply a boundary map twice on any linear combination of simplices, we get zero, and so the quantities $\beta_n(K(G))$ are always non-negative integers. These classically known numbers take the name of *Betti numbers* and, for each $n$, the $n$-th Betti number $\beta_n(K(G))$ corresponds to the dimension of the $n$-th homology space (with $\mathbf{Z}/2$-coefficients) of the clique complex $K(G)$ of $G$.

The intuitive sense of this construction is to count the "holes" that remain in the graph after we have filled all the directed cliques. In particular, the $n$-th Betti number is counting the $n$-dimensional holes. One can also see that $\beta_0$ counts the number of connected components of the graph. A classical result in topology shows a connection between the Euler characteristic and the Betti numbers, expressed by the identity: $\chi(K(G)) = \sum_{n=0}^N (-1)^n \beta_n(K(G))$, which gives another way of computing the Euler characteristic.

Notice that the construction of the directed clique complex of a given network $G$ does not involve any choice, and therefore, since the Betti numbers and the Euler characteristic of a simplicial complex are well-defined quantities for a simplicial complex , our constructions produce quantities that are well-defined for the network $G$, and we shall refer to them simply as the Euler characteristic and the Betti numbers of $G$.

## Boolean recurrent artificial neural networks

### Network structure and dynamics

The artificial recurrent neural networks consist of a finite number of Boolean neurons organized in layers with a convergent/divergent connection structure . The networks are composed by $50$ layers, each of them with $10$ Boolean neurons. The first layer is the input layer and all its $10$ neurons get activated at the same time at a fixed frequency of $0.1$, i.e. every $10$ time steps of the history. Each neuron in a layer is connected to a randomly uniformly distributed number of target neurons $f$ belonging to the next downstream layer. The networks include *recurrence* in their structure, meaning that a small fraction $g$ of the neurons appears in two different layers. This means that a neuron $k$ that is also identified as neuron $l$, is characterized by the union of the input connections of neurons $k$ and $l$, as well as by the union of their respective efferent projections.

The state $S_i(t)$ of a neuron $i$ take values $0$ (inactive) or $1$ (active) and all Boolean neurons are set inactive at the beginning of the simulation. The state $S_i(t)$ is a function of the its activation variable $V_i(t)$ and a threshold $\theta$, such that $S_i(t)  =  \mathcal{H}(V_i(t)-\theta)$. $\mathcal{H}$ is the Heaviside function, $\mathcal{H}(x)=0 : x<0$, $\mathcal{H}(x)=1 : x\ge0$. At each time step, the value $V_i(t)$ of the activation variable of the $i^{th}$ neuron is calculated such that $V_i(t+1) =  \sum_{j} S_j(t) w_{ji}(t)$, where $w_{ji}(t)$ are the weights of the directed connections from any $j^{th}$ neuron projecting to neuron $i$. The connection weights can only take four values, i.e. $w_1 = 0.1$, $w_2 = 0.2$, $w_3 = 0.4$, $w_4 = 0.8$. At the begin of the simulations all connection weights are randomly uniformly distributed among the four possible values. The weights of all the neurons are computed synchronously at each time step.

The network dynamics implements activity-dependent plasticity of the connection weights. Whenever the activation of a connection does not lead to the activation of its target neuron during an interval lasting $a$ time steps, its weight is weakened to the level immediately lower than the current one. Whenever the weight of a connection reaches the lowest level, the connection is removed from the network . Then, the pruning of the connections provokes the selection of the most significant ones and changes the topology of the network. Similarly, whenever a connection with a weight $w_m$ is activated at least $m+1$ consecutive time steps, the connection weight is strengthened to the level immediately higher than the current one. Hence, the parameter space of our simulations was defined by four parameters: the number $f$ of layer-to-layer downstream connections in the range $3$–$10$ by steps of 1, the small fraction $g$ of the neurons appearing in two different layers in the range $1$–$3$% by steps of 1%, the threshold of activation $\theta$ in the range $0.8$–$1.4$ by steps of 0.1, and the interval $a$ of the weakening dynamics of the connections in the range $7$–$9$ by steps of 1.

### Implementation of the simulations

The simulation software was implemented from scratch in Python. The network evolved with the dynamics explained above and the program computed the directed clique complex at each change of the network topology. For the entire network, the directed clique complex was computed each time the connectivity changed because of pruning. For the sub-network of the active nodes, the computation was carried out at each step of the simulation.

The computed directed clique complexes were used to compute the Euler characteristic both for the complexes representing the entire network and for the sub-complexes of the active nodes. To compute the directed clique complex of a network we used the implementation of the algorithm of in the `igraph` Python package , adapted to find directed cliques. The experiments were run in parallel on several CPUs using the tool GNU Parallel .

## Network filtrations

### Network structures

Many essential topological features of a network are determined by the distribution of edges over its graph. Different types of distributions result in different types of networks. For instance, pure random networks (RN) are formed assuming that edges in the network are independent of each other and they are equally likely to occur . For RN we have used the algorithm implemented in the Python package 'NetworkX' (https://networkx.github.io/) with the function 'erdos_renyi_graph' with parameters number of nodes $n=40$ and the probability for edge creation $p=0.2$.

These simple construction assumptions are generally not followed in networks obtained experimentally from ecological or gene systems, telecommunication networks or the Internet which are characterized by short average path lengths and high clustering, resulting in the so called small-world topology (SW) . For SW we used the same Python package 'NetworkX' with the function 'newman_watts_strogatz_graph' with parameters number of nodes $n=40$ and the number of connected neighbours in ring topology $k=20$ and the probability for adding a new edge $p=0.4$.

Other real-world networks such as brain, social networks, power grids and transportation networks exhibit topologies where more connected nodes, hubs, are more likely to receive new edges. The presence of these hubs and a power law distribution for the degree of the nodes defines scale-free networks (SF) . For SF we used the same Python package 'NetworkX' with the function 'barabasi_albert_graph' with parameters number of nodes $n=40$ and the number of edges to attach from a new node $m=10$.

### Network degree invariant

Given a directed network $G$, we define two filtrations by sub-networks (ordered sequences of networks in which each network is a sub-network of all the following ones) using the in- and out-degree of nodes. Let $ODF(G)$ be the out-degree filtration of $G$: the $i$-th network $ODF(G)_i$ in this filtration is the sub-network of $G$ induced by the vertices having out-degree at least $i$ and all the target nodes of their outgoing connections. In the same way we define the in-degree filtration $IDF(G)$: the $i$-th network $IDF(G)_i$ in this filtration is the sub-network of $G$ induced by the vertices having in-degree at least $i$ and all the source nodes of their incoming connections.

We computed the Euler characteristic for each network of the two filtrations, obtaining two sequences of integers, which are plotted to display a measure of the network topology, as a function of the degree levels of the filtration, normalized by the maximum degree present in the network. For example, let us consider the case illustrated in Figure B: one of the random networks with $n=40$ vertices that we have generated with a parameter $p=0.20$, as described above, had a maximum out-degree of its vertices equal to $19$. Therefore all the filtration levels have been divided by this value to normalize them (between 0 and 1).

For each network family (SF, RN, SW), we generated $N=50$ distinct networks with different seeds for the random numbers generator (the seeds were uniformly distributed integers in the interval $[1, 10000]$). We calculated the network degree filtration invariant sequences for in- and out-degree, which were then averaged for each network family and represented in Figure with the $95\%$ pointwise confidence bands.

# Competing interests

The authors declare that they have no competing interests.

# Author's contributions

Conceived and designed the experiments: PM, AEPV. Developed the mathematical construction, implemented the simulation, analyzed the results and drafted the manuscript: PM, AEPV. Both authors reviewed, read and approved the final manuscript.


abstract: Multifrequency media access control has been well understood in general wireless ad hoc networks, while in wireless sensor networks, researchers still focus on single frequency solutions. In wireless sensor networks, each device is typically equipped with a single radio transceiver and applications adopt much smaller packet sizes compared to those in general wireless ad hoc networks. Hence, the multifrequency MAC protocols proposed for general wireless ad hoc networks are not suitable for wireless sensor network applications, which we further demonstrate through our simulation experiments. In this article, we propose MMSN, which takes advantage of multifrequency availability while, at the same time, takes into consideration the restrictions of wireless sensor networks. Through extensive experiments, MMSN exhibits the prominent ability to utilize parallel transmissions among neighboring nodes. When multiple physical frequencies are available, it also achieves increased energy efficiency, demonstrating the ability to work against radio interference and the tolerance to a wide range of measured time synchronization errors.
author: GANG ZHOU YAFENG WU TING YAN TIAN HE CHENGDU HUANG JOHN A. STANKOVIC TAREK F. ABDELZAHER
bibliography: bibliography.bib
title: A Multifrequency MAC Specially Designed for Wireless Sensor Network Applications

# Introduction

As a new technology, Wireless Sensor Networks (WSNs) has a wide range of applications \[Culler 2001,Bahl 2002,Akyildiz 2001\], including environment monitoring, smart buildings, medical care, industrial and military applications. Among them, a recent trend is to develop commercial sensor networks that require pervasive sensing of both environment and human beings, for example, assisted living \[Akyildiz 2002,Harvard 2001,CROSSBOW\] and smart homes \[Harvard 2001,Adya 2001,CROSSBOW\].

> "For these applications, sensor devices are incorporated into human cloths \[Natarajan 2001,Zhou 2006,Bahl 2002,Adya 2001\] for monitoring health related information like EKG readings, fall detection, and voice recognition".

While collecting all these multimedia information \[Akyildiz 2002\] requires a high network throughput, off-the-shelf sensor devices only provide very limited bandwidth in a single channel: 19.2Kbps in MICA2 \[Bahl 2002\] and 250Kbps in MICAz.

In this article, we propose MMSN, abbreviation for Multifrequency Media access control for wireless Sensor Networks. The main contributions of this work can be summarized as follows.

-   To the best of our knowledge, the MMSN protocol is the first multifrequency MAC protocol especially designed for WSNs, in which each device is equipped with a single radio transceiver and the MAC layer packet size is very small.

-   Instead of using pairwise RTS/CTS frequency negotiation \[Adya 2001,Culler 2001; Tzamaloukas 2001; Zhou 2006\], we propose lightweight frequency assignments, which are good choices for many deployed comparatively static WSNs.

-   We develop new toggle transmission and snooping techniques to enable a single radio transceiver in a sensor device to achieve scalable performance, avoiding the nonscalable "one control channel + multiple data channels" design \[Natarajan 2001\].

# MMSN Protocol

## Frequency Assignment

We propose a suboptimal distribution to be used by each node, which is easy to compute and does not depend on the number of competing nodes. A natural candidate is an increasing geometric sequence, in which $$\label{eqn:01}
P(t)=\frac{b^{\frac{t+1}{T+1}}-b^{\frac{t}{T+1}}}{b-1},$$ where $t=0,{\ldots}\,,T$, and $b$ is a number greater than $1$.

In our algorithm, we use the suboptimal approach for simplicity and generality. We need to make the distribution of the selected back-off time slice at each node conform to what is shown in Equation (). It is implemented as follows: First, a random variable $\alpha$ with a uniform distribution within the interval $(0, 1)$ is generated on each node, then time slice $i$ is selected according to the following equation: $$i=\lfloor(T+1)\log_b[\alpha(b-1)+1]\rfloor.$$ It can be easily proven that the distribution of $i$ conforms to Equation ().

So protocols \[Bahl 2002,Culler 2001,Zhou 2006,Adya 2001,Culler 2001; Tzamaloukas-01; Akyildiz-01\] that use RTS/CTS controls[^1] for frequency negotiation and reservation are not suitable for WSN applications, even though they exhibit good performance in general wireless ad hoc networks.

### Exclusive Frequency Assignment

In exclusive frequency assignment, nodes first exchange their IDs among two communication hops so that each node knows its two-hop neighbors' IDs. In the second broadcast, each node beacons all neighbors' IDs it has collected during the first broadcast period.

#### Eavesdropping

Even though the even selection scheme leads to even sharing of available frequencies among any two-hop neighborhood, it involves a number of two-hop broadcasts. To reduce the communication cost, we propose a lightweight eavesdropping scheme.

## Basic Notations

As states, for each frequency number, each node calculates a random number (${\textit{Rnd}}_{\alpha}$) for itself and a random number (${\textit{Rnd}}_{\beta}$) for each of its two-hop neighbors with the same pseudorandom number generator.

Bus masters are divided into two disjoint sets, $\mathcal{M}_{RT}$ and $\mathcal{M}_{NRT}$.

RT Masters

:   $\mathcal{M}_{RT}=\{ \vec{m}_{1},\dots,\vec{m}_{n}\}$ denotes the $n$ RT masters issuing real-time constrained requests. To model the current request issued by an $\vec{m}_{i}$ in $\mathcal{M}_{RT}$, three parameters—the recurrence time $(r_i)$, the service cycle $(c_i)$, and the relative deadline $(d_i)$—are used, with their relationships.

NRT Masters

:   $\mathcal{M}_{NRT}=\{ \vec{m}_{n+1},\dots,\vec{m}_{n+m}\}$ is a set of $m$ masters issuing nonreal-time constrained requests. In our model, each $\vec{m}_{j}$ in $\mathcal{M}_{NRT}$ needs only one parameter, the service cycle, to model the current request it issues.

Here, a question may arise, since each node has a global ID. Why don't we just map nodes' IDs within two hops into a group of frequency numbers and assign those numbers to all nodes within two hops?

# Simulator

If the model checker requests successors of a state which are not created yet, the state space uses the simulator to create the successors on-the-fly. To create successor states the simulator conducts the following steps.

1.  Load state into microcontroller model.

2.  Determine assignments needed for resolving nondeterminism.

3.  For each assignment.

    1.  either call interrupt handler or simulate effect of next instruction, or

    2.  evaluate truth values of atomic propositions.

4.  Return resulting states.

shows a typical microcontroller C program that controls an automotive power window lift. The program is one of the programs used in the case study described in . At first sight, the programs looks like an ANSI C program. It contains function calls, assignments, if clauses, and while loops.

## Problem Formulation

The objective of variable coalescence-based offset assignment is to find both the coalescence scheme and the MWPC on the coalesced graph. We start with a few definitions and lemmas for variable coalescence.

*Proof.* C-MWPC can be easily reduced to the MWPC problem assuming a coalescence graph without any edge or a fully connected interference graph. Therefore, each C-node is an uncoalesced live range after value separation and C-PC is equivalent to PC. A fully connected interference graph is made possible when all live ranges interfere with each other. Thus, the C-MWPC problem is NP-complete. ◻

*Proof.* Simply, any solution to the MWPC is also a solution to the C-MWPC. But some solutions to C-MWPC may not apply to the MWPC (if any coalescing were made). ◻

# Performance Evaluation

During all the experiments, the Geographic Forwarding (GF) \[Akyildiz 2001\] routing protocol is used. GF exploits geographic information of nodes and conducts local data-forwarding to achieve end-to-end routing. Our simulation is configured according to the settings in . Each run lasts for 2 minutes and repeated 100 times. For each data value we present in the results, we also give its 90% confidence interval.

```latex
\tbl{Simulation Configuration\label{tab:one}}{%
\begin{tabular}{|l|l|}
\hline
TERRAIN{$^a$}   & (200m$\times$200m) Square\\\hline
Node Number     & 289\\\hline
Node Placement  & Uniform\\\hline
Application     & Many-to-Many/Gossip CBR Streams\\\hline
Payload Size    & 32 bytes\\\hline
Routing Layer   & GF\\\hline
MAC Layer       & CSMA/MMSN\\\hline
Radio Layer     & RADIO-ACCNOISE\\\hline
Radio Bandwidth & 250Kbps\\\hline
Radio Range     & 20m--45m\\\hline
\end{tabular}}
```

# Conclusions

In this article, we develop the first multifrequency MAC protocol for WSN applications in which each device adopts a single radio transceiver. The different MAC design requirements for WSNs and general wireless ad-hoc networks are compared, and a complete WSN multifrequency MAC design (MMSN) is put forth. During the MMSN design, we analyze and evaluate different choices for frequency assignments and also discuss the nonuniform back-off algorithms for the slotted media access design.

# Typical references in new ACM Reference Format

A paginated journal article , an enumerated journal article , a reference to an entire issue , a monograph (whole book) , a monograph/whole book in a series (see 2a in spec. document) , a divisible-book such as an anthology or compilation followed by the same example, however we only output the series if the volume number is given (so Editor00a's series should NOT be present since it has no vol. no.), a chapter in a divisible book , a chapter in a divisible book in a series , a multi-volume work as book , an article in a proceedings (of a conference, symposium, workshop for example) (paginated proceedings article) , a proceedings article with all possible elements , an example of an enumerated proceedings article , an informally published work , a doctoral dissertation , a master's thesis: , an online document / world wide web resource , , , a video game (Case 1) and (Case 2) and and (Case 3) a patent , work accepted for publication , 'YYYYb'-test for prolific author and . Other cites might contain 'duplicate' DOI and URLs (some SIAM articles) . Boris / Barbara Beeton: multi-volume works as books and .

# APPENDIX

In this appendix, we measure the channel switching time of Micaz \[CROSSBOW\] sensor devices. In our experiments, one mote alternatingly switches between Channels 11 and 12. Every time after the node switches to a channel, it sends out a packet immediately and then changes to a new channel as soon as the transmission is finished. We measure the number of packets the test mote can send in 10 seconds, denoted as $N_{1}$. In contrast, we also measure the same value of the test mote without switching channels, denoted as $N_{2}$. We calculate the channel-switching time $s$ as $$\begin{aligned}
%
s=\frac{10}{N_{1}}-\frac{10}{N_{2}}. \nonumber
\end{aligned}$$ By repeating the experiments 100 times, we get the average channel-switching time of Micaz motes: 24.3$\mu$s.

# This is an example of Appendix section head

Channel-switching time is measured as the time length it takes for motes to successfully switch from one channel to another. This parameter impacts the maximum network throughput, because motes cannot receive or send any packet during this period of time, and it also affects the efficiency of toggle snooping in MMSN, where motes need to sense through channels rapidly.

By repeating experiments 100 times, we get the average channel-switching time of Micaz motes: 24.3 $\mu$s. We then conduct the same experiments with different Micaz motes, as well as experiments with the transmitter switching from Channel 11 to other channels. In both scenarios, the channel-switching time does not have obvious changes. (In our experiments, all values are in the range of 23.6 $\mu$s to 24.9 $\mu$s.)

# Appendix section head

The primary consumer of energy in WSNs is idle listening. The key to reduce idle listening is executing low duty-cycle on nodes. Two primary approaches are considered in controlling duty-cycles in the MAC layer.

[^1]: RTS/CTS controls are required to be implemented by 802.11-compliant devices. They can be used as an optional mechanism to avoid Hidden Terminal Problems in the 802.11 standard and protocols based on those similar to \[Akyildiz 2001\] and \[Adya 2001\].