author: T. Tsujimoto; K. Bekki
date: Received 26 November 2010 / Accepted 14 February 2011
title: Implications of a non-universal IMF from C, N, and O abundances in very metal-poor Galactic stars and damped Ly$\alpha$ absorbers

# Introduction

Precise measurements of elemental abundances in the early Universe enable us to have direct access to the chemical enrichment process by studying nucleosynthesis in massive stars. Now we have two different kinds of available targets, i.e, extremely metal-poor (EMP) stars in the Galactic halo and high-z damped Ly$\alpha$ systems (DLAs) (e.g., Pettini et al. 2002, 2008). Their metallicity ranges are approximately -4$<$\[Fe/H\]$<$-2.5 and -3$<$\[Fe/H\]$<$-2, respectively. Several theoretical works claim that the abundances of EMP stars retain the nucleosynthesis yields of individual type II supernovae (SNe II) . In contrast, DLAs represent the abundances of the interstellar medium (ISM) integrated along the line of sight, which is equivalent to an average of all the mass of SNe II weighted with an initial mass function (IMF).

measured the N and O abundances of metal-poor DLAs . It is well-known that the N/O ratio of extragalactic H II regions exhibits a plateau at 12+$\log$ O/H .5ex8.0 \[Note that the solar O abundance is 8.69 (Asplund et al. 2009)\], which is considered to represent the primary production level of N. Their finding is, however, that some DLAs in their samples have a lower N/O ratio than the plateau level. This was firmly confirmed by with more additional samples (see left panel of Fig. 1). It is possible to infer from this that the observed low N/O ratio is attributable to the time delay of N ejection from asymptotic giant branch (AGB) stars with respect to O production by SNe II, provided that the primary N is mainly produced in AGB stars and star formation proceeds at a lower rate in DLAs than in extragalactic HII regions .

However, these predictions are incompatible with the observed N/O ratios of EMP stars. determined the N abundance of red giants, which is not influenced by any mixing inside evolved stars. They utilized the Li abundance measured in individual stars to diagnose the mixing level. Their result is shown by squares in the left panel of Fig. 1. We can see that an average N/O ratio of stars is broadly consistent with the plateau value. That implies that (i) massive stars synthesize a large amount of primary N, and (ii) the assembly of primary N from massive stars produces the plateau exhibited by extragalactic HII regions.

This high N abundance observed in EMP stars has inspired renewed interest in the nucleosynthesis of N in massive stars, which has led to the discovery of a new channel that enables high N production in rapidly rotating massive stars . Leaving more progress on this issue to future studies, the question is then how we interpret the low N/O ratio in DLAs if enough N is synthesized in massive stars to reach the level equivalent to the N/O ratios of EMP stars.

A clue to tackling this enigma can be found in the C abundance of DLAs, which is measured together with N and O abundances . As shown in the right panel of Fig. 1, the C/O ratios in DLAs (triangles) are somewhat higher than those of EMP stars (squares), which is the opposite of the trend found for the N/O ratio. The origin of C/O feature is indeed far easier to access than that of N/O because we have a much clearer understanding of the nucleosynthesis of C than N in massive stars .

In this paper, we propose a unified scheme that explains C, N, and O abundances in both DLAs and EMP stars, as well as in extragalactic HII regions. Our major claim is that if the IMF lacks very massive stars beyond $\sim$ 20-25 $M_{\odot}$ in metal-poor DLAs for some physical reasons that will be discussed in detail, both the observed high C/O and low N/O ratios are consistently explained. A theoretical scheme that does not invoke changes in the IMF based on the idea that DLAs have a different speed of chemical enrichment does not sufficiently explain their abundances, since this view predicts the same direction of change in both the C/O and N/O ratios against those of EMP stars. A model in which the major contribution of Population (Pop) III stars to chemical abundances is considered is not likely either, because the N /O ratios of EMP stars seem to reflect no clear difference in the production of N between Pop. II and Pop. III stars. In addition, a metallicity-dependent IMF, as well as metallicity-dependent yields among C, N, and O, may have to be considered as the origin of abundance trends in the Galactic halo but not as the mechanism of the contrast between DLAs and EMP stars.

The proposed IMF is supported by an observational result that suggests that the observed flux ratio of H$\alpha$ to the far-ultraviolet in low surface-brightness (LSB) galaxies favors such a truncated IMF , though there seems to be a wide range of other explanations for the observed trend (Lee et al., 2009, 2011a). In contrast to these truncated IMFs, there is observational evidence of a top-heavy IMF in massive galaxies at high redshifts . Accordingly, we insist on a non-universal IMF that varies among different galactic objects.

# How different are DLAs and stars?

The argument presented in this paper is based on the observationally implied contrast between the relative abundances of C, N, and O in Galactic halo stars and metal-poor DLAs. However, as shown in Fig. 1, mainly because of the large errors in each set of data, the extents of their differences seem insignificant and debatable. In this section, we do the statistical tests and also discuss the uncertainties involved in the abundance determination that may yield the systematic errors.

## C/O

We assess the extent of the differences in C/O between the two data groups, including the effect of the random error in the observed abundances of C/O and O/H of individual stars and DLAs. We compare all DLAs with the stars residing in $\log$ (O/H)+12$<$`<!-- -->`{=html}7.0 to ensure that both of the metallicity ranges broadly coincide. Since the C/O ratio of stars has a clear trend against O/H, we perform the least squares fit and see the C/O value at $\log$ (O/H)+12$=$`<!-- -->`{=html}6.5. We take this O/H value as the zero-point of a straight line. The same procedure is performed for DLAs and also for N/O in Sect. 2.2. As a result, we obtain $\log$ (C/O)=-0.68$\pm$`<!-- -->`{=html}0.04 for stars and -0.51$\pm$`<!-- -->`{=html}0.04 for DLAs.

In addition to random errors, there are systematic uncertainties because the stellar and DLA abundances are derived by different methods. Most notably, the efficiency of inelastic H collisions, whose cross-section remains essentially unknown, influences the values of stellar C/O and O/H. Here we adopt the case of a high collision efficiency that is adopted as a standard case in . If these collisions are unimportant, stellar and DLA abundances are in broad agreement . However, this results in an unusual trend of \[O/Fe\] against \[Fe/H\]: the \[O/Fe\] ratio of halo stars decreases towards lower \[Fe/H\] for \[Fe/H\]$<$-2.0 , which is inconsistent with the elemental ratios of other $\alpha$-elements to Fe. Moreover, a finding of a very high C abundance in a DLA with \[Fe/H\]=-3 is likely to agree with the contrast of C/O between DLAs and Galactic stars.

## N/O

We compare all EMP stars with the DLAs residing in $\log$ (O/H)+12$<$`<!-- -->`{=html}7.2 for the same reason as in the C/O case. At $\log$ (O/H)+12$=$`<!-- -->`{=html}6.5, $\log$ (N/O)=-1.64$\pm$`<!-- -->`{=html}0.1 for stars and -2.1$\pm$`<!-- -->`{=html}0.07 for DLAs are deduced. Here for two DLAs, we assume their mean values to be their upper limits. The difference between stars and DLAs becomes smaller when we exclude the four stars at $\log$ (O/H)+12$<$`<!-- -->`{=html}6.2, where there is no data for DLAs. This stellar ensemble then gives $\log$ (C/O)=-1.82$\pm$`<!-- -->`{=html}0.23.

For the systematic errors, we note that the stellar abundance of N is derived from the ultraviolet NH band, whose physical parameters are not yet well established . Thus, the possibility of the stellar N/O ratios being systematically somewhat smaller than those in Fig. 1 cannot be ruled out. The N/O ratio measured in unevolved dwarf stars includes three stars in our metallicity range, i.e., $\log$ (O/H)+12$<$`<!-- -->`{=html}7.2. Their average N/O value is $\log$ (N/O)=-1.49, which seems broadly consistent with the analysis of . Moreover, if we add three samples to Spite's data in the range of 6.2$<$$\log$ (O/H)+12$<$`<!-- -->`{=html}7.2, $\log$ (N/O) at $\log$ (O/H)+12$=$`<!-- -->`{=html}6.5 is estimated to be -1.59$\pm$`<!-- -->`{=html}0.18.

We perfome the Kolmogorov-Smirnov (KS) test, even though it is not likely to provide a meaningful statistical check because of the small number of samples. We determine the probabilities that the two data sets are drawn from the same population are less than 1% and 9% for each case of C/O and N/O, respectively. In summary, the statistical analysis including the random errors in the abundances, suggests a relatively higher C/O and lower N/O ratios for DLAs than for EMP stars. However, this systematic difference remains uncertain until the potential systematic errors inherent in the stellar abundances are erased. Therefore, we await future work to improve the abundance determination with much more data before validating the argument presented in the following sections. Bearing in mind that the contrast between the relative abundances of C, N, and O claimed here must be viewed with caution, we attempt a theoretical interpretation of the stellar and DLA abundances.

# Analysis of C/O ratio

As shown in the right panel of Fig. 1, the C/O ratio of stars clearly decreases as metallicity increases for $\log$ (O/H) + 12 $<$ 7.0. For theoretical interpretations of the origin of this feature, proposed a differential C production in accordance with the stellar rotation, while proposed a time-dependent IMF. However, the former scenario predicts a more strongly decreasing trend of N/O than that of C/O because the stellar rotation influences N production more efficiently, which is not the case for the observed N/O trend (left panel). As for the latter case, there is no supporting evidence for it in other elemental ratios. In addition, the claimed initially top-heavy IMF contradicts our present result, as we show in Sects. 5 and 6.

Here we present the theoretical interpretation of the C/O ratio observed in EMP stars, which gives a consistent view of the N/O trend and is developed into models of chemical evolution for the Galaxy and DLAs.

## Nucleosynthesis C/O

We first see the nucleosynthesis results for C and O in massive stars, and show that the nucleosynthesis C/O ratio has a clear trend with respect to the progenitor mass. The left panels of Fig. 2 show the C and O yields in massive stars as a function of the progenitor mass, as calculated by three groups (circles, Woosley & Weaver 1995; triangles, Nomoto et al. 1997; squares, Ekström et al. 2008). There exist some differences in their models. The two former groups present the results of SN models, while the calculation by does not include the stage of SN explosion. Since both C and O are not synthesized or broken by the SN shock, the three results can be equally compared. Another distinction is that the two SN models do not include the stellar rotation for the progenitor models, while the models by include very fast rotation, which leads to a high C yield compared with the SN models. These comparisons reveal (i) good agreement with the mass-dependent O yield, but (ii) a variation in the C yield, partly due to the difference in the stellar rotation.

Irrespective of the C variation in the different models, the resultant C/O ratios have a trend in common with the progenitor mass, i.e., a lower progenitor mass model results in a higher C/O ratio (right panel). Here, the C/O ratios denoted by \[C/O\] are normalized by a solar C/O ratio of -0.26. We emphasize that the range of the predicted \[C/O\] ratios with different progenitor masses covers an observed \[C/O\] range of -0.7 to 0 where there are EMP stars.

## Comparison with the observed C/O

According to the claim based on the analysis of elemental abundance patterns of EMP stars , we assume that the formation of EMP stars is triggered by a SN remnant (SNR), thus the stars formed retain the abundance pattern of this SN. This hypothesis leads to the comparison of the observed correlation of C/O and O/H for each EMP star with those of individual theoretical models. For this purpose, we calculate each model's \[O/H\] ratio obtained from the O yield in each SN model (lower left panel of Fig. 2) divided by the mass of hydrogen swept up by an SNR, as done by .

The mass $M_{\rm SW}$ swept up by an SNR is approximated by the formula

$$\begin{aligned}
M_{\rm SW}=5.1 &\times& 10^4 M_\odot \left (\frac{E_0}{10^{51} \ {\rm ergs}}\right )^{0.97}  \nonumber \\
&\times& n_1^{-0.062}\left (\frac{C_s}{10 \ {\rm km \ s}^{-1}}\right )^{-9/7} \ \ \ ,
\end{aligned}$$

where $E_0$, $n_1$, and $C_s$ denotes the explosion energy of SN, the number density of the ISM, and the speed of sound, respectively. In the models of , $E_0$ varies from 1.28$\times 10^{51}$ to 3.01$\times 10^{51}$ ergs, while a single value of $10^{51}$ ergs is assigned to all models in . Since $M_{\rm SW}$ is insensitive to $n_1$, we determine the values of $M_{\rm SW}$ by adjusting $C_s$ so as to fit with the observed correlation and our choice is $C_s$=17 km s$^{-1}$ for the two SN models. On the other hand, for the models of massive stars by , we adopt the values $E_0$= 10$^{51}$ ergs and $C_s$ = 6 km s$^{-1}$.

The resultant C/O-O/H correlations predicted by individual SN models are compared with the observation in Fig. 3. From the results given by each model, the points that reside in the range of -1$<$\[C/O\]$<+0.2$ and \[O/H\]$>$-3 are plotted. We see a good coincidence between the prediction and the observed data of EMP stars. Moreover, the average C/O ratio weighted by a Salpeter IMF utilizing the yields of gives a good agreement with the C/O ratios for stars with -1.5$<$\[O/H\]$<$-1. Accordingly, we conclude that the decrease in C/O with increasing O/H is a reflection of the elemental feature of an individual SNR. In the end, the overall evolution of C/O imprinted in stellar abundances (i.e., squares) in Fig. 1 can be interpreted as follows. We can roughly divide its evolution into three stages that result in an increase in O/H: (i) for $\log$ (O/H)+12$<$`<!-- -->`{=html}7.0, the decreasing trend reflects the nucleosynthesis of the C/O ratio for individual SNe II; (ii) for 7.0$<$ $\log$ (O/H)+12$<$`<!-- -->`{=html}8.0, there is a plateau as a result of a mixture of C/O synthesized in various SNe II with different masses; and (iii) for $\log$ (O/H)+12$>$`<!-- -->`{=html}8.0, an increasing trend appears owing to a time-delayed additional C ejection from AGB stars.

The claim that we highlight is that a C/O ratio expected from individual SNe II is, theoretically and observationally, likely to be higher for a less massive progenitor mass. In this case, emergence of a kind of IMF, such as a cut-off of the upper mass end of massive stars results in a higher C/O ratio in the ISM at an early evolutionary stage, when only SNe II contribute to chemical enrichment. This situation will be incorporated into the models of chemical evolution and shown to be a likely case for metal-poor DLAs.

# Chemical evolution model

We first construct the model for Galactic chemical evolution (GCE) as a standard case. Subsequently, we try to reproduce the chemical evolution of other systems such as DLAs and extragalactic H II regions, modifying the GCE model by changing the fundamental model parameters, including the IMF and/or the star formation rate (SFR).

## Galactic chemical evolution model

The basic picture of the model is that the disk was formed through a continuous infall of material from outside the disk region over 14 Gyr. As is common in chemical evolution studies, we assume that the gas is distributed uniformly and that the heavy elements are well-mixed within the zone. In this framework, the SFR is assumed to be proportional to the gas fraction with a constant rate coefficient $\nu$ of 0.4 Gyr$^{-1}$ . For the infall rate, we adopt the formula that is proportional to $t\exp(-t/\tau_{\rm in})$ with a timescale of infall $\tau_{\rm in}$=5 Gyr . The metallicity of an infall is assumed to be \[O/H\]=-4 with a SN-II-like elemental ratio. In this study, we treat the evolution of C, N, and O. From the implication of the nucleosynthesis ratios in SNe II from EMP stars (see Fig. 1), we assume that $\log$ C/O =-0.6 and $\log$ N/O=-1.5 for infalling gas. Owing to a very small amount of metals in the infalling gas, both stellar and gas abundances for \[O/H\]$>$-3 (i.e., $\log$ (O/H)$+12>$-5.7) are unaffected by the abundances of an infall.

For the nucleosynthesis yields of SNe II, we adopt the results by for C and O since their average nucleosynthesis C/O ratio is in good agreement with the observed C/O level. We discuss the N yield in Sect. 5. For the nucleosynthesis yields of type Ia SNe (SNe Ia), the O yield is taken from with zero yields for C and N. In addition, we assume metallicity-dependent C and N yields from AGB stars. As a basic yield, we adopt the results calculated with $Z$=0.001, $n_{\rm AGB}$=4, and $m_{\rm HBB}$=0.8 by , and we choose its metallicity dependence so as to follow the observed increasing feature seen in both C/O and N/O. These adopted nucleosynthesis yields are then combined with a power-law mass spectrum with a slope $x$ of -1.35, e.g., a Salpeter IMF, and the mass range (0.05 $M_{\odot}$, 40 $M_{\odot}$). The lower mass limit is deduced from the theoretical analysis of the mass-to-luminosity ratio in the solar neighborhood . For the occurrence frequency of SNe Ia, we assume that the fraction of the stars that eventually produce SNe Ia for $3-8$$M_{\odot}$ is 0.05 with the lifetime of SN Ia progenitors of 0.5$-$`<!-- -->`{=html}3 Gyr .

## Applications to other systems

Our prime purpose is to reproduce the chemical evolution of DLAs. Based on the knowledge acquired by analysing the C/O ratio in Sect. 2, we introduce the IMF with the upper mass end at $<$ 40 $M_{\odot}$ into the DLA case. In addition, we take into account the hypothesis that DLAs are metal-poor systems corresponding to dwarf galaxies or the outer disks of spiral galaxies . Accordingly, we assume a low SFR and a long timescale of infall for DLAs.

In the evolution of N/O, many extragalactic H II regions are likely to follow the path where chemical enrichment proceeds more efficiently than in the solar neighborhood and ends up with a more metal-rich ISM with $\log$(O/H)+12$>$`<!-- -->`{=html}9.0. In contrast to the DLA case, we assume a high SFR and a short timescale of infall for these metal-rich H II regions. The samples of HII regions with $\log$(O/H)+12$>$`<!-- -->`{=html}9.0 are indeed those in bright massive galaxies, mainly M31, M83, M51(an interacting galaxy), and NGC 4254, whose SFRs are expected to be higher than in the solar neighborhood. In addition, we introduce a top-heavy IMF, since the above choices do not agree sufficiently well with the observed trend of increasing N/O for these objects. In other words, since we find that the model with a high SFR and a normal IMF fails to reproduce the high N/O ratio observed in metal-rich H II regions, a top-heavy IMF is invoked (see Sect. 6).

In the end, we adopt IMFs that differ from Salpeter's for two objects to reduce or enhance the heavy-element yields ejected from SNe II. There are basically two approaches: varying the value of $m_u$ or changing the slope $x$. For the canonical case, $m_u$=40 $M_{\odot}$ and $x$=-1.35 are assumed. Since a steeper IMF, i.e., $x<$-1.35, does not give a satisfactory result because of its small effect on the elemental ratios for our adopted yields, $m_u$ is chosen as a variable parameter to reduce the yields. The restricted nucleosynthesis results lead to our choice of a flatter IMF, i.e., $x>$-1.35 to enhance the yields.

# C/O evolution

The predicted C/O evolution for the Galaxy is shown by the blue solid curve in Fig. 4. It clearly reproduces the observational measurements for $\log$ (O/H) + 12 $>$ 7.0, which include those of Galactic halo stars, extragalactic H II regions, and the Sun. The differences between the prediction and the observed data seen for $\log$ (O/H) + 12 $<$ 7.0 occur because in such a very metal-poor regime, stellar abundances (data of EMP stars) do not represent the abundance of the ISM (prediction) at each epoch when a star is formed, as already discussed.

The observed data of DLAs can be directly compared with the model results. The C/O ratio in DLAs is expected to be higher than that of the ISM in an early Galaxy. To change the relative abundances of C and O, we introduce an IMF with a variation in the upper mass limit, $m_u$. For two cases of $m_u$, 20$M_{\odot}$ and 25$M_{\odot}$, the calculations are performed with the choice of parameters ($\nu$, $t_{\rm in}$)=(0.3 Gr$^{-1}$, 7 Gyr). The results are shown by the red dashed and solid curves, respectively. As expected from the mass-dependence of the nucleosynthesis C/O ratio shown in the right panel of Fig. 2, the smaller $m_u$ corresponds to a higher C/O ratio in the ISM that is compatible with the observed level of DLAs. Here we note that a smaller $m_u$ is found to lower the C/O ratio more efficiently than a steeper slope of the IMF, which is another mechanism for reducing the number of very massive stars.

# N/O evolution

We continue our investigation of the N/O evolution, starting with the discussion of the N yield. In contrast to the case of the C yield, the nucleosynthesis results for N in massive stars vary significantly among different groups, thus we are far from converging on one reliable view. Therefore, we deduce the N yield as functions of the progenitor mass based on the theoretical argument about the observed N/O ratios of EMP stars, as well as those of DLAs combined with the knowledge previously acquired from the discussion of the C/O ratio. This deduction has two major steps:

\(i\) In contrast to the C/O feature, the N/O ratios of EMP stars apparently exhibit no clear trend. This can be interpreted as an outcome of (1) very massive stars, the first contributor of chemical enrichment, making a broad level of the observed N/O ratio, and (2) their prior contamination of the ISM with a high N/O ratio that prevents a lower N/O ratio from less massive stars from being established in stellar abundances. In other words, the N/O ratios of EMP stars roughly correspond to the nucleosynthesis N/O in very massive stars while the information on less massive stars is not established because of with a time delay their small amount of N production relative to more massive stars. Recent nucleosynthesis results support the idea that the N/O ratios in very massive stars broadly give the observed level of the N/O ratio for EMP stars.

\(ii\) After detecting the low N production in less massive stars and fixing the yield at the most massive mass, the small amount of N synthesized in less massive stars is not so critical. We can estimate the amount by referring to nucleosynthesis results and with supplementary information about the discussion of the C/O ratio: the abundances of DLAs are likely to reflect the nucleosynthesis products from less massive ($<$ $\sim$`<!-- -->`{=html}20-25$M_{\odot}$) stars. Thus, we assume that the nucleosynthesis products in 20-25 $M_{\odot}$ imply that log (N/O)$\sim$ -2.

Through these assessments, we finally obtain the N yield shown in the upper panel of Fig. 5.

Utilizing the deduced N yield, we first model the metal-rich systems aiming for the Galaxy and extragalactic H II regions, and show their results by blue and red curves, respectively, in the middle panel of Fig. 5. For the extragalactic H II regions, we set ($\nu$, $t_{\rm in}$)=(5 Gyr$^{-1}$, 1 Gyr), together with a flatter IMF of $x$=1.05. The choice of a top-heavy IMF, in addition to a high rate of star formation, is necessary to accelerate chemical enrichment during an early phase: the N ejected from AGB stars starts to occur on a very short timescale ($\sim 4\times 10^7$yr) because the major source of N, AGB stars, are massive (6-8 $M_{\odot}$). As an effect of a top-heavy IMF, the N/O ratio in an early phase becomes higher than the Galaxy case.

We emphasize that HII regions with a relatively low $\log$(N/O) for $\log$(O/H)+12.5ex8.0 are reproduced by the model with a flatter IMF. We know from observations that the more luminous disk galaxies are more metal-rich because they are redder . It is thus possible to claim that more luminous disk galaxies are likely to have a flatter IMF. Here, we propose that the N/O ratio in individual disk galaxies at a given metallicity is the more stringent parameter needed to deduce their IMFs. Investigation of the correlation between the N/O ratio and physical properties of disk galaxies such as the mean mass density and the total mass will surely provide a clearer understanding of the determinant factor for the IMF in disk galaxies. In addition, HII regions with a flatter IMF are predicted to exhibit a lower C/O ratio on the whole than those with a canonical IMF for the same reason, which will be confirmed by future observations.

The opposite situation corresponding to metal-poor DLAs is considered in the lower panel of Fig. 5. The N/O evolution is calculated with the models in which the upper mass in the IMF ends at 25 $M_{\odot}$ (red curve) and 20 $M_{\odot}$ (blue curve). Both models closely reproduce the overall chemical evolution of DLAs. Their most conspicuous property is represented by an early and steep increase in N/O owing to the ejection of N from AGB stars. This property is caused by the suppression of the ejection of O from SNe II because the IMF lacks very massive stars, while the contribution of N from AGB stars remains at the same level. In addition, we see a large scatter in the N/O ratio of DLAs, some of which deviate from our predicted paths. This may reflect the diverse chemical evolution of individual DLAs— some DLAs may have a canonical IMF, while some chemical enrichment proceeds with a higher SFR.

We add one comment. The low N abundance in DLAs is reminiscent of the Magellanic Clouds (MCs), which exhibit a severe deficiency of N in H II regions . However, an IMF without very massive stars, as predicted in metal-poor DLAs, is not likely to be representative of the MCs. Our prediction is that the upper-mass cut-off causes a low N/O but a high C/O, whereas the MCs also exhibit a deficiency of C abundance. Moreover, the low N/O occurs only in an early phase, before AGB stars begin to contribute to the N enrichment. Neither is a top-heavy IMF, which accelerates chemical enrichment and thereby results in the low N/O for young H II regions. propose that the origin of a low N/O in the MCs may be caused by an external factor, i.e., a dilution by a gas with a severely deficient N abundance such as high-velocity clouds, an idea that needs more research.

# Discussion

## Other evidence of a non-universal IMF

Although we first propose that the formation of very massive stars (masses higher than $\sim$ 20 $M_{\odot}$) is truncated in DLAs, observations have suggested that the formation of massive stars is truncated or significantly suppressed by the IMF with a steeper slope in some local galaxies. investigated the extinction-corrected flux ratio $F_{\rm H \alpha}/f_{\rm FUV}$ (where $F_{\rm H \alpha}$ and $f_{\rm FUV}$ are the observed H$\alpha$ flux and far-ultraviolet one, respectively) in various galaxies and found that the ratio is smaller for galaxies with lower surface-brightnesses in the $R$-band (see their Fig. 6). They considered the possibility that $F_{\rm H \alpha}$ comes mostly from massive O-type stars with masses higher than 20 $M_{\odot}$ that can ionize the ISM, while $f_{\rm FUV}$ comes from less massive stars. In the end, they concluded that the observed correlation results from there being a smaller number of massive O-type stars in lower surface-brightness galaxies. That led to their claim of a nonuniform IMF for galaxies in the local Universe. Furthermore, investigated the IMFs for a large number of galaxies with different magnitudes in a database of the Sloan Digital Sky Survey and found that the IMFs in bright galaxies such as the Milky Way are similar to the Salpeter IMF, whereas fainter galaxies tend to have steeper IMFs. Their results imply that fewer massive stars are formed in low-luminosity galaxies.

The possible truncation of the formation of massive stars in low-luminosity galaxies and LSB galaxies is discussed for other samples of galaxies . claimed that the flux ratio $F_{\rm H \alpha}/f_{\rm FUV}$ is barely related to galaxy properties such as stellar mass and effective surface brightness in late-type galaxies, and argued against the truncation of stars in these galaxies. In addition, the results by implied that a deficiency of H$\alpha$ possibly results from Poisson fluctuations in the formation of massive stars for a universal IMF. Thus, it is safe to say that it would be premature to claim that the truncated IMF is characteristic of low-luminosity galaxies and LSB galaxies.

From a theoretical aspect, the high-mass end ($m_{\rm u}$) has been found to depend on the mass of star clusters, thus on the physical properties of galaxies where the clusters are formed . If more massive clusters are less efficiently formed in the ISM with a lower gas density as in the outer parts of galaxies, $m_{\rm u}$ would be lower . This theoretical result appears to support a scenario where $m_{\rm u}$ is inclined to be lower in less massive and LSB galaxies of low mean mass-density and thereby the formation of massive star clusters is suppressed.

To explain the observed physical properties of massive galaxies, previous studies, both observational and theoretical, pointed out that the IMFs need to be top heavy . For example, the observed $\alpha$/Fe in typical $L_{\star}$ ellipticals can be reproduced by semi-analytic models based on the hierarchical galaxy formation scenario only if top-heavy IMFs are adopted . The ISM of these massive galaxies with top-heavy IMFs at high redshifts are predicted to exhibit low C/O and high N/O ratios, in contrast to the properties of high-$z$ DLAs. In this respect, massive galaxies are unlikely to be identified with DLAs.

Several theoretical works have suggested that the slope of IMFs in galaxies evolves with time in the sense that the slope becomes flatter, i.e., top-heavier in accordance with the growth of the galaxy . If we take into account the dependence of the IMFs on the physical properties of galaxies, we find that the formation of massive stars is correlated with the gas mass density of the ISM. We are then led to a picture where, as galaxy masses grow and surface mass densities become higher, galactic IMFs change from truncated IMFs into canonical ones. We recall, however, that the chemical abundances of the local halo and disk stars show no evidence of a variable IMF in the past.

Finally, there is strong evidence of a top-heavy IMF in our own Galaxy. The chemical properties of the Galactic bulge stars, especially their stellar metallicity distribution, are most consistent with a top-heavy IMF .

## Origin of the possible truncated IMF in DLAs

As discussed above, $m_{\rm u}$ can be smaller in dwarf galaxies and LSB galaxies. Our prediction is that a small $m_{\rm u}$ results in the relatively high C/O and low N/O ratios in the ISM. Therefore, this result is expected in these galaxies at an early epoch of their evolution. After all, the observed abundance pattern for high-z DLAs suggests that DLAs are gas-rich dwarf galaxies or LSBs. In addition, the gas-rich outer parts of high-$z$ disk galaxies represent another possible explanation of DLAs since the local outer disks may have lower $m_{\rm u}$ and thus lack H$\alpha$ regions . Future observations designed to measure the N/O and C/O ratios in the outer HII regions of high-$z$ disk galaxies will validate their connection with DLAs, as well as their potential to have truncated IMFs. we note that owing to a low $m_{\rm u}$ in gas-rich dwarfs and the outer part of disk galaxies, the formation of HII regions is significantly suppressed so that their HI-rich ISM can be identified with DLAs. The small $m_{\rm u}$ would increase the possibility of galactic ISM being detected as HI-rich DLAs.

For UV-selected galaxies, reveal the correlation of the N/O ratio with the luminosity. In their samples, galaxies exhibiting log N/O $<$ -1.6 are restricted to fainter ones with $M_{\rm B}>$-18, which correspond to dwarf galaxies. According to our present study, the correlation is caused by the IMF lacking massive stars in dwarf galaxies. In addition, the similarity of the N/O ratio in dwarf galaxies and DLAs supports the idea that DLAs are gas-rich parts of dwarf galaxies at high redshifts. DLAs can indeed be identified with dwarf galaxies .

# Conclusions

We have presented a new view of the chemical evolution of DLAs as well as extragalactic H II regions from our careful analysis of their C, N, and O abundances. We highlight two nucleosynthesis results: (i) the C/O ratio has a clear correlation with the progenitor mass in the sense that a less massive SN II yields a higher C/O, and (ii) the finding of a new channel for the very high N production in the most massive stars. These combined with the hypothesis that the abundances of EMP stars in the Galactic halo reflect those of individual SN remnants naturally lead to a theoretical interpretation for the abundances of DLAs. Metal-poor DLAs have puzzling properties, a relatively high C/O ratio and relatively low N/O ratio compared to EMP stars. Our claim is that these properties are compelling evidence that the IMF lacks very massive stars in metal-poor DLAs.

In contrast, we have found that some metal-rich extragalactic H II regions in local disk galaxies appear to have top-heavy IMFs. This suggests that the IMFs are likely to differ even among disk galaxies. There were earlier claims that high-$z$ early-type galaxies possess a top-heavy IMF . Our additional findings resulting from a detailed chemical analysis focusing on DLAs strengthen supporting evidence of a non-universal IMF from high to low redshifts in the Universe. However, on the other hand, there is extensive evidence of a universal IMF . Future work will definitely be required to validate the IMF variation. Finally, we emphasize that analysing the C/O and N/O abundances more precisely in galaxies, including our own, will help us to understand the IMF variation.

The authors wish to thank an anonymous referee for all valuable comments that helped improve the paper. TT is assisted in part by Grant-in-Aid for Scientific Research (21540246) of the Japanese Ministry of Education, Culture, Sports, Science and Technology.


abstract: The search for high energy neutrinos of astrophysical origin is being conducted today with two water/ice Cherenkov experiments. New instruments of higher performance are now in construction and more are in the R&D phase. No sources have been found to date. Upper limits on neutrino fluxes are approaching model predictions. Results are reported on the search for point sources, diffuse fluxes, gamma ray bursts, dark matter and other sources.
address: Department of Physics, University of Wisconsin, Madison, U.S.A.
author: Albrecht Karle
title: Observations of high energy neutrinos with water/ice neutrino telescopes.

# Introduction

High energy neutrino astronomy will offer a new view of the Universe. The observation of the highest energy cosmic rays and the detection of galactic and extragalactic sources of gamma rays are proof that very energetic sources of non-thermal radiation exist. The detection of sources of high energy neutrinos will provide insights into the origin of the highest energy radiation. Neutrinos have been detected from the Sun and from Supernova 1987a; both observations representing fundamental breakthroughs in astrophysics and physics. Neutrino telescopes now in operation and construction aim at higher energies. This report will focus on the energy range aimed for by under water/ice neutrino telescopes: $\approx 10^{11}$ eV to $\approx 10^{18}$eV. At the same time an increasing effort is underway to detect neutrinos at the highest energies to beyond $10^{20}$ eV. Markov and Zheleznikh suggested the detection of neutrinos in water via the process $\nu_l(\bar{\nu_l})+N \rightarrow l^{\pm} + X$ of upward or horizontal neutrinos interacting with a nucleon $N$ of the matter surrounding the detector. A significant fraction of the neutrino energy will be carried away by the produced lepton. For the important case of the muon the angle between the parent neutrino and the muon is less than $1^{\circ}$ at 1 TeV energy and it becomes very small at higher energies providing the basis for astronomy.

In the past decade high energy neutrino astronomy has evolved from exploratory experiments to the first generation of neutrino telescopes with substantial neutrino detection rates. The first attempt was made by the DUMAND project in the deep Pacific Ocean near Hawaii. This pioneering effort opened the road to the realization of the first working neutrino telescopes, AMANDA-B10 (1997, ), Baikal (NT-200 1998, ) and AMANDA-II (2000). Two telescopes are in construction, IceCube , at the South Pole, and ANTARES , in the Mediterranean Sea, and two projects are in the planning and prototyping phase in the Mediterranean Sea: NESTOR and NEMO. The sensitivity of some of the current experiments, Baikal and AMANDA-II, is already below some model predictions. This report will focus on the current status of observations.

# Detection principle and energy ranges

Water Cherenkov detectors detect the Cherenkov photons that are emitted by relativistic particles. Strings of optical sensors which are deployed in a transparent medium detect the Cherenkov photons of muon tracks or cascades. Measuring the photon arrival times with a resolution of a few nanoseconds allows an accurate reconstruction of the muon track.

Neutrino telescopes cover a wide energy range for two reasons. The first reason is that the neutrino-nucleon cross section increases substantially with energy. The second is that the muon range and energy loss increase substantially with energy. We single out muons because they allow the best angular resolution. The cascade channel (e's, $\tau$'s and neutral currents) however allow in many areas a comparable sensitivity. Figure shows the effective detection area for muon neutrinos in AMANDA-II and ANTARES for point source selection, for the Ultra High Energy (UHE) analysis of AMANDA B10.

# Detectors

## AMANDA

AMANDA is a neutrino telescope built and operated at the Geographic South Pole. The final detector configuration, called AMANDA-II, consists of 677 optical modules (OM) arranged on 19 vertical strings deployed at depths between 1300-2400 m. Most of the OMs are deployed at a depth of 1500-1950 m and arranged within a diameter of 200 m. All analog pulses of the 20 cm photomultipliers are transmitted to the surface on individual cables, with most of the channels using fiber optic cables. The analog transmission on fibers allows low dispersion (PMT pulse FWHM about 20 nsec), high time resolution and a high dynamic range while keeping the electronics in the ice at a minimum.

The IceCube neutrino observatory at the South Pole will consist of 4800 optical sensors, installed on 80 strings between the depths of 1450 m to 2450 m in the Antarctic Ice, and 320 sensors deployed in 160 ice Cherenkov tanks (referred to as IceTop) on the ice surface directly above the strings. Each sensor consists of a 25 cm photomultiplier tube, connected to a waveform-recording data acquisition circuit capable of resolving pulses at nsec time resolution and high dynamic range. In January 2005 76 such sensors were installed as a first part of the IceCube array. The collaboration reported first results which show that the detector system works as expected . The collaboration hopes to do first science with an array of 10 or more strings in 2006 and complete construction in 2010/11.

## Baikal

The Baikal neutrino telescope is located 30 km off the shore of Lake Baikal, Siberia, at a depth of 1.1 km. The configuration NT200 was commissioned in April, 1998. It consists of 192 optical modules (OMs), mounted on 8 vertical strings, which form a structure of 40 m diameter and 72 m height. Each OM contains a 37 cm diameter photomultiplier. For installation and maintenance, the Baikal collaboration takes advantage of the cold winters in Siberia, which cause the Lake to freeze over for several months allowing easy access to the deployment site. The Baikal Collaboration added three outer strings in 2005 to increase the sensitivity at high energies for cascades.

## ANTARES

Antares is a neutrino telescope under construction in the Mediteranean Sea at a depth of 2500 m near Toulon, France. The project aims to complete the full array, consisting of 12 lines with 75 photomultipliers each, by 2007. Simulations indicate a very good angular resolution of $0.3^{\circ}$ for 1 TeV $\mu's$ ($0.6^{\circ}$ degrees with respect to the neutrino direction). The neutrino effective area is expected to be similar to the AMANDA array as shown in fig. .

## Other projects and initiatives

NESTOR is a neutrino telescope with a proposed location near Pylos, on the Greek Ionian Sea coast. A calibration and engineering run of a test detector was carried out in 2003 . The detector was operated for more than one month and data was continuously transmitted to shore. With the data collected, the collaboration was able to reconstruct muon events at a rate consistent with predictions and produce a muon zenith angle distribution .

The initial goal of the NEMO project is research and development to lay the technological foundation for a future km scale detector in the Mediterranean. The collaboration has identified a possible site for a $1\,km^{3}•$-scale instrument at a depth of 3500 m, at a distance of about 80 km from the Sicilian coast. The NEMO architecture for a $1\,km^{3}•$ detector is based on 5832 optical sensors. Finally, it should be mentioned that the groups involved in the Meditarranean projects have formed a European consortium called KM3Net to further the development of a kilometer scale neutrino telescope in the Mediterranean Sea.

# Atmospheric muons and neutrinos

Cosmic ray muons form a background for neutrino telescopes. The background is rejected in two ways which illustrate two fundamental modes of operation:

1.  At low energies - below a few PeV: Reject downgoing muons by zenith angle and use upward muons as reliable neutrino messengers.

2.  At very high energies: Reject cosmic ray muon background by direction and energy. At energies above $\approx 10^{15-16}$ eV, (depending on the depth of the instrument) the down going background is negligible, especially at larger zenith angles. This allows to accept more downgoing signal events.

At the depth of AMANDA it is found that the down-going cosmic ray muon flux is about 5 orders of magnitude larger than the expected neutrino flux. At greater depth the muon flux is significantly smaller and vice versa. Neutrino telescopes can use the cosmic ray muon flux to study the detector response, possible systematic effects and also to perform physics measurements. The cosmic ray muon angular and vertical intensities were measured in AMANDA-II and they agree well with simulations and with other experimental measurements . AMANDA was able to use down-going muon events, which are in coincidence with SPASE surface array, to make a measurement of the cosmic ray mass composition as a function of the primary energy near the knee .

In order to measure the atmospheric neutrino component, a rejection factor of more than $\sim 10^6$ is needed to eliminate the cosmic ray muon background. The AMANDA collaboration records about 4 well reconstructed neutrino events per day in AMANDA-II. The preliminary unfolded $\nu_{\mu}$ sea-level energy spectrum was determined and is in agreement with the expectations up to above $10^4$ GeV . The observed spectrum is shown in figure .

# Search for sources of astrophysical neutrinos

## Diffuse fluxes

The search for diffuse astrophysical neutrino fluxes is in many respects the most challenging. The reason is that atmospheric neutrino background rejection is almost entirely based on energy reconstruction, and misreconstructed downgoing cosmic ray muons can complicate the analysis. A diffuse flux can only be measured if the spectrum is significantly harder than the steeply falling atmospheric neutrino spectrum.

Models of cosmic ray shock acceleration based on the Fermi mechanism naturally produce energy spectra with a power index $\gamma \approx 2$. $E^{-2}$-type spectra are often used as a reference flux. However, each specific model needs to be optimized and analyzed separately and the energy cut will move higher or lower depending on the hypothetical signal flux. The limits discussed here are not differential limits, instead they refer to an assumed $E^{-2}$ spectrum extending over the whole energy range. The accelerated protons interact with each other or with ambient photons. Pion decay yields a $\nu-$flux with roughly the same spectral index as the protons. At the production site, the $\nu-$flavors are distributed according to the ratios $\nu_{e} : \nu_{\mu} :\nu_{\tau} = 1 : 2 : 0$. Due to neutrino oscillations, an approximate equipartition of flavors is predicted to be observed at the detector. However, the assumption of a $\nu_{e} : \nu_{\mu} :\nu_{\tau} = 1 : 1:
1$ ratio is indeed a simplification, which is not accurate at energies above $\approx 10^{14}•\,eV$ .

We categorize the results in three groups and give the current results for Baikal and AMANDA. The reported results are compiled in the figure . AMANDA has reported results on all three channels, Baikal on the third one.

-   Through-going upward muon events:
    For the through-going upward muons AMANDA reports results in two analyses. In one analysis, an energy cut is applied and observed events above the energy cut are compared to the signal and background prediction. The optimization is done before the data is unblinded. The collaboration reported a sensitivity at the level of $\approx 10^{-7}•$GeV/(cm$^{2}$•s sr for a data set currently being analyzed. The other result is based on a neutrino energy spectrum analysis with a differential limit at about 30 TeV based on an energy unfolding method .

-   Contained cascade events:
    This analysis aims for the measurement of the diffuse flux of $\nu_{\e}$, $\nu_{\tau}$ and due to neutral current interactions. In order to determine a generic astrophysical flux limit the assumption is made that astrophysical neutrino fluxes appear in a $\nu_{e}:\nu_{\mu}:\nu_{\tau}$=1:1:1 ratio. While the directional resolution is relatively poor, a comparable flux limit is obtained due to the better energy resolution for cascades compared to $\nu_{\mu}$.

-   Not through-going muon and not contained cascade events:
    The Baikal collaboration was able to develop an analysis technique which allowed the effective detector volume for very bright cascade events to increase substantially over the geometric volume. Even though non contained events are accepted with a vertex up 200 m outside the instrument, the technique allows to reject downgoing background effectively by using a combination of energy and directionality of the timing pattern. The atmospheric neutrino background is easily rejected based on energy. An upper limit is reported based on 806 days of livetime. AMANDA developed a technique to detect non-contained events of very high energy. In this analysis the 2 billion cosmic ray muon events are rejected using the energy. This UHE analysis is sensitive to signal flux in the region from about 10 PeV to 1000 PeV.

## Search for point sources

The detection of astrophysical point sources of neutrinos is perhaps the most prominent and also name bearing goal of neutrino astronomy. Numerous sources of TeV energy gamma rays have been detected in the past decade. Neutrino telescopes not only overlap at low energies with the energy range of current ground based gamma ray telescopes, but also extend in sensitivity to much higher energies.

Results are available from the underground detectors MACRO and SuperKamiokande and also from Baikal and AMANDA.

None of the current observations have indicated any significant excess above the background of atmospheric neutrinos. Limits are reported for individual sources of special interest, which are also sources of gamma rays at TeV energies. Figure shows a number of individual limits for point sources from 4 years (807 days of livetime) of AMANDA in the northern hemisphere. The limts are at the level of E$^{2}\cdot6\cdot 10^{-7}$•GeV cm$^{-2}$sec$^{-2}•$ . In a few cases, for example the active galaxy Mrk 501, the limits on neutrino fluxes approach the observed gamma ray flux at TeV energy of these sources in their high states. Also shown are limits of the MACRO experiment from more than 10 years of data. The average sensitivity for point sources with an $E^{2}$-spectrum is shown for the lifetime of MACRO, for 4 years of AMANDA and for one year of IceCube.

AMANDA reported an upper limit from a search for neutrinos from the galactic plane in the Northern hemisphere . A larger km scale neutrino telescope will be needed to detect the a neutrino flux from the galactic plane. The Northern hemisphere large instrument will be ideally positioned to search for the higher flux in near the galactic center.

## Neutrinos from Gamma Ray Bursts, dark matter and other searches

Gamma-ray bursts (GRBs) are among the most energetic phenomena in the universe. The observation of high energy neutrinos in gamma ray bursts would confirm hadronic acceleration in the fireball; possibly revealing an acceleration mechanism for the highest energy cosmic rays. AMANDA reported searches for neutrinos from a total 451 GRBs, most of which were triggered by the BATSE telescope during 1997 to 2000. The short duration, typically 10 to 100 seconds, of the bursts makes background rejection easy. AMANDA observed zero neutrino events at an expected background of approximately 2 events for all bursts. Upper limits have been reported at the level of 3 $\cdot$ 10$^{-8}$ GeV•cm$^{-2}$•sr$^{-1}$•s$^{-1}•$ for 139 bursts from 2000 to 2003. This limit is about an order of magnitude above the Waxman-Bahcall prediction .

## Dark matter and other searches

Minimally Supersymmetric extensions to the Standard Model predict the existence of the neutralino, in the mass range GeV-TeV, which is a candidate for the cold dark matter. The particles will lose energy and become gravitationally trapped in the centers of object like the earth, the sun, or the center of the galaxy where they annihilate to produce neutrinos. These neutrinos may be detected in a high-energy neutrino telescope. The energy range of the predicted neutrino fluxes is from about 10$^{10}$•eV to a few times 10$^{12}$eV. AMANDA has performed a search for dark matter in the center of the Earth and the Sun . No excesses of events have been observed and upper limits were placed. Baikal reported an upper limit (see for example ) at a comparable level on the neutrino flux from WIMPs annihilating in the Earth. SuperKamiokande holds the lowest limit by a factor of 2 on the muon neutrino flux from the Sun at a level of 2000 km$^{-2}$•yr$^{-1}$ for neutralino masses around 10$^{12}$eV (based on a data set from many years). Limits from underground detectors like Baksan and MACRO are still comparable or lower, especially for lower neutralino masses around 10$^{11}$eV. ANTARES is expected to have a sensitivity of a factor 10 below the current limits; IceCube is expected to have a final sensitivity of about 2 orders of magnitude below the current AMANDA limit. Nothern hemisphere telescopes will be in a good position to also search for dark matter from the center of the Galaxy.

## Neutrinos of highest energies

New techniques are being explored to search for neutrinos in the energy range from $10^{18}$•eV to beyond $10^{21}$•eV. Some exploratory experiments are already underway. A recent overview of sensitivities of techniques and instruments can be found in ref. . Most of the current efforts have focused on the use of radio signals from highest energy events in ice. The application of the radio technique in salt domes is being investigated. Acoustic techniques are being explored to measure the acoustic pulse of GZK events. First data are available from a prototype version of ANITA, which has performed a first scan of the Antarctic ice sheet for neutrino radio pulses from a balloon . RICE uses radio sensors deployed in the deep ice near the South Pole. The techniques are still in rapid development and a substantial increase in sensitivity at highest energies can be expected in the next decade.

# Acknowledgments

This work was supported by the NSF under OPP-0337726 grant. The author thanks T. Montaruli, D. Boersma, G. Hill and S. Grullon for helpful comments.

# References


author: Krzysztof Chalupka
Computation and Neural Systems
California Insistute of Technology
E-mail: firstname.lastname@example.com Michael Dickinson
Biology
California Institute of Technology Pietro Perona
Electrical Engineering
California Institute of Technology
bibliography: GRM.bib
title: Generalized Regressive Motion: a Visual Cue to Collision

# Abstract

Brains and sensory systems evolved to guide motion. Central to this task is controlling the approach to stationary obstacles and detecting moving organisms. Looming has been proposed as the main monocular visual cue for detecting the approach of other animals and avoiding collisions with stationary obstacles. Elegant neural mechanisms for looming detection have been found in the brain of insects and vertebrates. However, looming has not been analyzed in the context of collisions between two moving animals. We propose an alternative strategy, *generalized regressive motion* (GRM), which is consistent with recently observed behavior in fruit flies. Geometric analysis proves that GRM is a reliable cue to collision among conspecifics, whereas agent-based modeling suggests that GRM is a better cue than looming as a means to detect approach, prevent collisions and maintain mobility.

# INTRODUCTION

Animals move to forage, approach potential mates, chase prey, escape from predators, and to maintain their position within a group. They must do so in environments that are often cluttered by rocks, plants and other moving animals. Whether the goal is making or avoiding contact, it is valuable to detect the proximity of both stationary and moving entities.

Looming, defined by Gibson as a visual pattern expanding symmetrically on the retina , is commonly believed to be a robust and reliable monocular visual cue to impending collision[^1]. It is often understood as a loosely-defined group of visual stimuli rather than a specific mechanism, with an underlying idea that an object approaching with a constant velocity produces expanding patterns on the observer's retina. When an animal is stationary, looming is a sufficient cue to detect approaching objects. When a moving animal is on a collision course with a stationary obstacle, time-to-collision can be estimated from looming patterns even when distance is unknown . It is generally accepted that looming is a cue used by various animals to avoid stationary obstacles, and elegant neural mechanisms for its detection have been unveiled. Experiments have revealed looming-sensitive neural pathways in many animals. The DCMD/LGMD neurons of the locust  as well as the pigeon nucleus rotundus  and the goldfish Mauthner cell  respond to divergence of image edges. Finally, the fruitfly uses looming-sensitive neurons during navigation (; see also Discussion).

Looming has been analyzed in the setting where either the animal or the obstacle is stationary. This is in contrast to regressive motion, used by the fruitfly to avoid collisions among multiple moving animals . Neural mechanisms behind regressive motion-driven behavior are unknown , but correlational motion detectors likely used by the fruitfly  can form a solid basis for regressive motion detection. The ecological usefulness of regressive motion has not been explored by . We build on the results of  by providing a theoretical and practical analysis of Generalized Regressive Motion (GRM) – visual stimulus similar to regressive motion but simpler to detect and more versatile. Whereas regressive motion occurs when there is clockwise motion in the left visual hemifield or counter-clockwise motion in the right visual hemifield, GRM occurs if there is clockwise motion in the left eye or counter-clockwise motion in the right eye, as shown in Fig. . Experiments by  admit the hypothesis that the fruitfly uses GRM, not pure regressive motion, for stopping. We show that GRM enjoys the advantages of both looming motion and regressive motion. Our contribution is threefold:

# GEOMETRY OF REGRESSIVE MOTION

Our geometric analysis of GRM is based on an abstract model of an agent: It is a point in the Euclidean plane, equipped with two "eyes" – centers of projection. Each agent has a well-defined orientation, which allows us to define its Contralateral Visual Angle (CVA). The CVA is the angle subtended by the nasal boundary of each eye, as in Fig.  (right). For now, we assume the distance between the eyes is zero and identify their position with the position of the agent (see Fig. ). This assumption is justified if the modeled animal's inter-eye distance is small compared to its typical distance from other animals. We nevertheless drop this assumption in simulations (described below), where we consider agents with two separate eyes and spatially extended bodies.

## GRM detection can prevent all collisions among moving conspecifics

Let a point (which could be stationary or in movement) project to azimuthal position $\phi$ on the observer agent's eye. Let $\phi=0$ be the direction in front of the agent, positive angles for the left side and negative angles for the right side, and restrict $\pi\in [-\pi,\pi)$. Denote the point's angular velocity as $\dot{\phi}$.

**Definition 1**. *A point projecting at $\phi$ is in *regressive motion* with respect to the observer if $\dot{\phi} \cdot \phi \leq 0$, and in progressive motion otherwise.*

Before proving our main theorem, we state an easy propsition whose proof we relegate to the Mathematical Appendix.

**Proposition 2**.  * Let the relative position and velocity of the observed object be $\mathbf{x} \text{ and } \mathbf{v}$ respectively. Then $\dot{\phi} = \frac{1}{\|\mathbf{x}\|^2}\langle \mathbf{v}^\perp,
  \mathbf{x} \rangle$. In particular, $\dot{\phi}$ scales as one over distance squared.*

The following theorem ensures that before any potential collision, one of the agents will perceive regressive motion.

**Theorem 3**. *Let f1 and f2 be two agents moving on straight, intersecting trajectories. If f1 reaches the trajectory intersection after f2, f1 perceives regressive motion at all times before f2 reaches the intersection and progressive motion afterwards. f2 perceives progressive motion before f1 reaches the intersection, and regressive motion afterward. *

*Proof.* Let two point-agents f1, f2 move on a flat uniform surface with constant velocities on intersecting (that is, non-parallel) trajectories. Align the reference frame's y-axis with the direction of f1's movement, and place the origin at the point at which the agents' trajectories cross. The situation is fully described by four parameters (see Fig. Left):

$v^1,v^2$

:   – the speed of f1 and f2 respectively,

$\Psi$

:   – the angle f2's velocity vector makes with f1's velocity (also called the *angle of approach*), and

$d$

:   – the y-coordinate of f1 at the moment when f2 reaches the origin.

First, we compute the angular position $\phi_{21}$ of f2 on f1's projection center and the angular velocity $\dot{\phi}_{21}$, at the moment when f1 is at distance $d+\epsilon$ from $0$.

For $\epsilon=0$, the positions of the two agents are respectively $$\begin{aligned}
    \mathbf{x}^1(d) &= (0, d),\\
    \mathbf{x}^2(d) &= (0, 0).

\end{aligned}$$ If $\epsilon \neq 0$, the time that passed since the original configuration is $\Delta t = \epsilon/
  v^1$, and since the velocity of f2 is $\mathbf{v}^2 = v^2(-\sin{\Psi}, \cos{\Psi})$, we have $$\begin{aligned}
    \mathbf{x}^1(d+\epsilon) &= (0, d+\epsilon),\\
    \mathbf{x}^2(d+\epsilon) &= \epsilon\frac{v^2}{v^1}(-\sin{\Psi}, \cos{\Psi}).

\end{aligned}$$ Define $\mathbf{x}^R, \mathbf{v}^R$ to be the relative position and velocity of f2 in f1's frame of reference. Then $$\begin{aligned}
    \mathbf{x}^R(d+\epsilon) &= \mathbf{x}^2(d+\epsilon) - \mathbf{x}^1(d+\epsilon)\\
    &= (-\frac{\epsilon v^2}{v^1}\sin{\Psi},
    \frac{\epsilon v^2}{v^1}\cos{\Psi}-(d+\epsilon)),\\
    \mathbf{v}^R(d+\epsilon) &= (-v^2\sin{\Psi}, v^2\cos{\Psi}-v^1).

\end{aligned}$$ From this we can directly compute $\phi_{21}$, and Proposition 2 enables us to compute the angular velocity of f2: $$\begin{aligned}
    \phi_{21} &= \arctan{(x^R_2/x^R_1)}-\pi/2,\\
    \dot{\phi}_{21} &= \frac{1}{D^2}\langle \mathbf{v}^{R\perp}, \mathbf{x}^R\rangle,

\end{aligned}$$ where $D^2=(x^R_1)^2+(x^R_2)^2$ is the distance between the two agents. Plugging in the values calculated above we get $$\begin{aligned}
    \phi_{21}(d,\epsilon) &=
    \arctan{\frac{\epsilon\frac{v^2}{v^1}\cos{\Psi}-(d+\epsilon)}{-\epsilon\frac{v^2}{v^1}\sin{\Psi}}}-\pi/2,\label{eq:theoryphi}\\
    \dot{\phi}_{21}(d,\epsilon) &= -\frac{1}{D^2}dv^2\sin{\Psi},\label{eq:theorydotphi}

\end{aligned}$$ with $D^2 =
  (\frac{v^2\epsilon}{v^1})^2+(d+\epsilon)^2-2\epsilon\frac{v^2(d+\epsilon)}{v^1}\cos{\Psi}.$ Now, assume that f2 arrives at the intersection first, that is $d<0$. From Equation it follows that $$\dot{\phi}_{21}(d,\epsilon) \geq 0 \iff \sin{\Psi} \geq 0.$$ But this implies that the denominator in the arctangent in Equation is nonnegative if and only if $\epsilon\leq0$, $$-\epsilon \frac{v_2}{v_1} \sin{\Psi} > 0 \iff \epsilon \leq 0,$$ A positive denominator restrict the range of the arctangent to $[-\pi/2,
    \pi/2]$, and thus $$\epsilon \leq 0 \iff \phi \in [-\pi, 0].$$

Thus, $$\dot{\phi}_{21}(d,\epsilon) \geq 0 \implies (\phi(d,\epsilon) \in [-\pi, 0] \iff \epsilon \leq 0).$$ Analogously, $$\dot{\phi}_{21}(d,\epsilon) \leq 0 \implies (\phi(d,\epsilon) \in
        [0, \pi] \iff \epsilon \leq 0).$$ This proves the first part of the theorem. The second part is proven exactly in the same way, but switching the reference frame to that of f2. ◻

The theorem easily generalizes to GRM:

**Definition 4**. *Let f2 project onto a projection center f1, with azimuthal position $\phi_{21}$. f1 perceives GRM if and only if $\left(\dot{\phi}_{21}>0 \text{ and } \phi_{21} \in [-\pi, \text{CVA}]\right)$ or $\left(\dot{\phi}_{21}<0 \text{ and } \phi_{21} \in [-\text{CVA},
    \pi]\right)$, where CVA is a fixed angle.*

**Theorem 5**. *Let f1 and f2 be two agents moving on straight, intersecting trajectories. If f1 reaches the trajectory intersection after f2, f1 perceives GRM at all times before f2 reaches the intersection.*

This follows from Theorem  because regressive motion implies GRM (see definitions above).

## GRM detection can prevent collisions with stationary objects

An agent using solely non-generalized regressive motion detection will always collide with stationary obstacles, which project expanding (progressive) patterns. However, GRM with CVA $> 0$ and appropriate motion thresholding provides a mechanism for stationary collision avoidance. Let T~GRM~ denote the smallest magnitude of GRM that causes the agent to stop. Intuitively, larger CVA's and smaller T~GRM~ provide better obstacle detection.

**Theorem 6**. *An agent on a collision course with a stationary object will perceive GRM before the collision, as long as T~GRM~$< \infty$ and CVA $>$ 0.*

*Proof.* We assume the agent is approaching an object such that the centerline of the agent $\phi=0$ does not point directly at a (non-smooth) corner of the object. We can then assume there is a neighborhood around the $\phi=0$ aziumuth that can be approximated as a wall segment. As in Fig.  (right), place the origin at the position of the agent as shown. The agent approaches the wall at an angle $0 < \alpha < \pi/2$ and with positive speed $v$, its velocity $\mathbf{v}=(v\sin{\alpha}, v\cos{\alpha})$. Consider an arbitrary wall-point $(x,y)$, marked red in the figure. The relative velocity of the point w. r. t. the agent is $-\mathbf{v}$, and its angular velocity on the agent's eye equals (by Proposition 2) $$\dot{\phi}(x,y,\alpha,v)=\frac{-v}{x^2+y^2}(x\cos{\alpha}-y\sin{\alpha}).\label{eq:wall_angvel}$$ We need to take into account very small $\theta$ and very large $T$. the following proposition implies the theorem: $$\forall_{\theta,T>0}\:\exists_{X>0} \text{ s.t. }
    \left\{\begin{array}{l}
      \|\dot{\phi}(X,y,\alpha,v)\|>T \text{  and  }\\
      y \tan{(\alpha-\theta)}>X>y\tan{\alpha}
    \end{array},\right.$$ where the last condition restricts the point we're searching for to be between the leftmost edge of perceived GRM and the center of the agent's visual field. Now, fix $X=y\tan{(\frac{\alpha-\theta}{k})}, k>1$. $k$ can always be chosen to make the point $(X,y)$ arbitrarily close to the centerline, and so contained in the smooth wall-like neighborhood on the obstacle. Clearly these $X,y$ satisfy the second condition above. We also have $$\dot{\phi}(X,y,\alpha,v)=\frac{-v\left(y\left[\tan{\frac{\alpha-\theta}{k}}\cos{\alpha}-\sin{\alpha}\right]\right)}{y^2(\tan{\frac{\alpha-\theta}{k}}+1)}.$$ Since $\alpha,\; \theta$ and $k$ are constant, this expression scales as $1/y$ and thus reaches arbitrarily large magnitudes as $y$ approaches $0$ – that is, as we place the observer closer and closer to the wall. ◻

The intuition behind this theorem is that a GRM detector detects any motion in a cone symmetrical about the center of the visual field. If the agent frontally approaches an object, the target is guaranteed to produce a strong enough signal within that cone at some positive time before the collision occurs.

## False alarms, safety and mobility

Theorems  and  give basis to the claim that GRM can be useful for collision avoidance. However, it can be argued that they are of limited practical use. One problematic area not explored by the theorems is that of false alarms. An agent using GRM as a stopping cue can stop unnecessarily in a variety of situations, some of which are shown in Fig. .

If perception of GRM of any magnitude greater than zero caused the agents to stop, they would be perfectly safe from collisions, but unable to move. In simulations described below we investigate the role of two parameters that enable GRM-based agents to trade-off mobility and safety from collisions: the CVA and tne threshold T~GRM~ on the magnitude of GRM that stops the agent. Varying these parameters in a population of agents changes the population's safety and mobility, where safety corresponds to the fraction of prevented collisions, and mobility to the fraction of useful stops. Formally, we can classify any encounter between agents f1 and f2 as

True Positive (TP)

:   : agent f1 stops due to perceived motion of agent f2, and *f1 would collide with f2* had both f1 and f2 continued to move with velocities they had at the moment of f1's stop, and f2 is outside of f1's collision radius.

True Negative (TN)

:   : agent f1 moves without stopping and doesn't collide with any entity.

False Positive (FP)

:   : agent f1 stops due to perceived motion of entity f2, and *f1 would not collide with f2* had both f1 and f2 continued to move at velocities they had at the moment of f1's stop. In addition, f2 is outside of f1's collision radius at the time of stopping.

False Negative (FN)

:   : agent f1 collides with agent f2.

We can then define mobility and safety as $$\begin{aligned}
 \text{mobility} &= \frac{TP}{TP+FP} \text{ and}\\
 \text{safety}  &= \frac{TP}{TP+FN},

\end{aligned}$$ Mobility is high if and only if the agent rarely stops without a good reason. Safety is a complementary measure that is high if and only if the agent avoided many out of all the potential collisions. In our view, any collision avoidance algorithm is useful inasmuch as it offers a range of good mobility-safety tradeoffs: It can be used to make mobile vehicles remain relatively safe (relative to other algorithms), as well as safe vehicles that retain good mobility.

# SIMULATIONS

A number of issues are not covered by our theory:

-   real agents have extended bodies, unlike the geometric points considered in Thm ,

-   real agents have visual systems with multiple centers of projection and

-   there is a tradeoff of safety and mobility to be explored.

Studying the tradeoff between mobility and safety is best done experimentally in a simulated environment. It is difficult to derive theoretical tradeoff curves given the statistical variability of the trajectories even in simple environments. Computational simulations are also a good tool for studying GRM agents with extended bodies and two separate eyes. Thus, to further study GRM we use agent-based modeling with populations of fly-like agents trying to avoid collisions. The agents use GRM- and looming-based algorithms for stopping. We compare the performance of GRM-based and looming-based collision avoidance using population safety and mobility as performance metrics. The details of the simulation setup are available in Methods below. Matlab code implementing the simulations is available online at `http://vision.caltech.edu/`$\sim$`kchalupk/cod` `e.html`.

## GRM detection offers good mobility and safety to a population of conspecifics

In each simulation, ten fly-like agents were placed in a toroidal arena and followed straight trajectories with constant speed. Each agent was equipped with a stopping mechanism triggered by the perception of GRM with specific CVA and T~GRM~(consistent across all the agents in a given simulation run). We performed 100 types of simulations, varying the CVA and T~GRM~ values systematically[^2]. After running multiple trials for each simulation type, we calculated the safety and mobility of the agent population in each case. Fig.  shows that varying the GRM parameters offers a wide variety of safety-mobility tradeoffs to the population.

## Looming detection offers poor mobility and safety to a population of conspecifics

As a point of reference we measured the usefulness of looming as a signal for collision avoidance. To this goal, we equipped each agent with both a GRM and a looming detector. We then performed a series of simulations varying CVA, T~GRM~ and T~LOOM~, where the latter is a looming threshold (Methods contains a detailed description of the stopping mechanism). Fig.  (left) shows safety and mobility in the simulations where only the looming signal was used for stopping (that is, T~GRM~ was very high). The figure shows that to achieve 95% safety, the agents had to stop unnecessarily 75% of the time. Fig.  (middle, right) shows full simulation data: each marker corresponds to one (CVA, T~GRM~, T~LOOM~) setting, and the three parameters vary independently. In Fig.  (middle) the value of T~GRM~ varies smoothly on the upper envelope of the scatter plot. It is very clear that choosing low T~GRM~ offers good safety but bad mobility, whereas higher T~GRM~ increases safety but decreases mobility. Fig.  (right) shows that the correspondence between T~LOOM~ and the safety-mobility tradeoff is much less clear. Whereas low values of T~LOOM~ decrease mobility of the agents, increasing the value past a certain point offers no additional safety-mobility advantages when GRM is also used for collision avoidance.

# METHODS: SIMULATION PARAMETERS

Each simulation puts ten fly-like agents, each defined by 14 visible points (see Fig. ), in a toroidal arena — a square of side length 50mm with opposing edges glued together.

The flies follow simple dynamics, described below. For each combination of $\text{CVA} \in \{0,10,...,90\}$ degrees and T~GRM~, T~LOOM~$\in \{0.1, 1, 2, 4, 6, 8, 10, 12, 14, 32\}$ rad/s, we ran 50 trials of 10000 time-steps (to a total of $1000\times 50$ simulations, 50s long each, of flies walking at realistic speeds). Movies S1-S4, available online at `http://vision.caltech.edu/`$\sim$`kchalupk/cod` `e.html`, show four example simulation runs resulting from different parameter settings (Movie Captions Appendix describes symbols used in the movies.)

```latex
\begin{table*}[!ht]\centering
   \caption{{\bf Parameters defining the fly model.}}
   \label{tab:parameters}
   \begin{tabular}{ccl}
     \toprule
     \textbf{Symbol} & \textbf{Default value} & \textbf{Meaning}\\
     \midrule
     $\Delta t$ & 0.005 s & Time-step in numerical simulations\\
     $R$ & 50 mm & Edge length of the walking arena (glued into a torus)\\
     $N$ & 10 & Number of flies in the arena\\
     $l$ & 2 mm & Length of an agent\\
     $d$ & 0.55 mm & distance between the eyes' centers\\
     $v_{\min}$ &10 mm/s & minimum agent speed \\
     $v_{\max}$ & 30 mm/s & maximum agent speed\\
     $v_i$ & $ \sim U(v_{\min}, v_{\max})$ & walking velocity of agent $i$\\
     $p_{01}$ & $\sim $0.8 & prob. stop-to-walk in a 1s time-interval\\
     $P_{01}$ & 0.008 & prob. stop-to-walk in one time-interval (i.e. $(1-p_{01})^{\Delta t} = (1-P_{01})$)\\
     T\textsubscript{LOOM}& 0$^\circ/s$ ($\dag$)& stopping threshold on looming motion\\
     T\textsubscript{GRM}& 0$^\circ/s$ ($\dag$)& stopping threshold on regressive motion\\
     $\text{CVA}$ & 30$^\circ$ ($\dag$) & the CVA\\
     $\theta_i$ & 120$^\circ$ & angle of ipsilateral visual field seen by each eye\\
     $\delta_\sigma$ & 30$^\circ$ &increment of standard deviation of agent reorientation motions $\sigma_{\vec{v}}$\\
     $\lambda_\sigma$ & 0.992 & decay constant for $\sigma_{\vec{v}}$ at each $\Delta t$ \\
     $n$ & 14 & number of points on each agent \\
     &($\dag$) & variable whose value is systematically explored in some experiments\\
     \bottomrule
   \end{tabular}
 \end{table*}
```

Numerical implementation:

:   the motion and control of each agent are computed at discrete time-intervals with constant time increments of $\Delta t$. For simplicity of notation where we write $t+1$ in the following, we mean $t + \Delta t$. The value of $\Delta t$ is given in Table  alongside all other simulation parameters.

Agent trajectories:

:   the $i$-th agent's trajectory is determined by walking speed $v_i$, initial position $\vec{x}_i^0$ and initial direction $\vec{v}_i$ with $\| \vec{v}_i\| = 1$, i.e. $\vec{x}_i(t)
       = \vec{x}_i^0 + t v_i \vec{v}_i$. Each agent has a different constant velocity chosen uniformly from the \[1-3\] cm/s interval.

Control:

:   while it is walking, each agent keeps walking at constant velocity until it is stopped by a GRM or looming percept caused by another agent (as explained below).

Spontaneous start:

:   when an agent is stationary it flips a coin at each time step: if the coin turns out to be heads, the agent starts moving in the current direction at its preferred velocity. If the coin is tails it stays put. The probability of obtaining heads in the interval of one second is denoted $p_{01}$; thus, the probability of obtaining tails in one time-step is $(1 -
       p_{01})^{\Delta t}$.

Orientation:

:   while an agent is walking, it keeps constant velocity and orientation (i.e. $\vec{v}_i (t+1) = \vec{v}_i(t)$). Upon stopping it samples a new orientation from a Gaussian pdf with standard deviation equal to the current value of a parameter $\sigma_i$ and centered at the current orientation. In addition, after reorientation, $\sigma_i$ increases by a fixed amount $\delta_\sigma$. $\sigma_i$ decays exponentially, with decay constant $\lambda_\sigma$, i.e. $\sigma_i(t+k\Delta t) =
       \sigma_i(t)\lambda_\sigma^k$. This mechanism, modeling basic neuronal sensitization, is a simple way to allow the flies to keep roughly straight trajectories when encountering transient obstacles (other moving flies), and avoid getting stuck around large static obstacles (groups of stopped flies).

GRM:

:   each eye sees GRM whenever the angular motion of any point projecting onto its retina is directed contralaterally, i.e. counterclockwise for the right eye and clockwise for the left eye. Each eye's visual field goes beyond the frontal direction to cover a given CVA. The CVA is a free parameter which we study to discover the best compromise between avoiding collisions and false alarms.

GRM stops:

:   agent boundaries are defined by 14 points visible to other flies, as shown in Fig. . Call the $j$-th point on the $i$-th agent $p_j^i$ and its azimuthal position on the observer's eye $\phi(p_j^i)$. Each agent measures the angular velocity of all the points on the other flies. If *any* GRM $\dot{\phi}(p_j^i)$ is detected, the agent compares the GRM magnitude $\|\dot{\phi}(p_j^i)\|$ to a threshold T~GRM~ and stops if $\|\dot{\phi}(p_j^i)\| >$ T~GRM~.

Looming motion:

:   Let $\|\dot{\phi}_L\|$ denote the largest magnitude of counter-clockwise motion that any point evokes in the left visual hemifield, and $\|\dot{\phi}_R\|$ the largest magnitude of clockwise motion evoked in the right visual hemifield. Then the strength of looming perceived by the agent equals $\omega_{LOOM} := \min{(\|\dot{\phi}_L\|,\|\dot{\phi}_R\|)}$.

Looming stops:

:   similarly to GRM stopping, the agent stops if $\omega_{LOOM} >\;\;$ T~LOOM~. This simple mechanism activates only if the agent can perceive points diverging at velocities larger than T~LOOM~.

All simulation parameter values are specified in Table . To summarize, agent $i$'s motion is governed by the following equations (see also the diagram in Fig. ).

In order to simplify the notation we omit index $i$ unless necessary, and write out the equations for a one-eyed agent; since both looming and GRM are monocular cues in our implementation, the extension to the two-eyed agent is trivial.

Control:

:   Let $z \in \{0, 1\}$ be the variable denoting whether an agent is stationary ($z=0$) or in motion ($z=1$). Let $u(t)$ be an i.i.d. random process with uniform probability density on $(0,
       1)$. Let $\|\dot{\phi}_{GRM}(t)\|$ be the largest observed magnitude of GRM at time t. Let $\omega_{LOOM}(t) := \min{(\|\dot{\phi}_L(t)\|,\|\dot{\phi}_R(t)\|)}$ denote the looming strength observed at time t.

    $$\begin{aligned}
         z(0) &= \begin{array}{l l} 1 & \quad \text{ initially the agent is set in motion } \end{array}\\
         z(t+1) & = \left\{
           \begin{array}{l l} 0 & \quad\text{ if }  z(t)=1 \text{ and}\\
                                & \quad\phantom{W}(\|\dot{\phi}_{GRM}(t)\| >  \text{T\textsubscript{GRM} or}\\
                                & \quad\phantom{W} \omega_{LOOM}(t) > \text{T\textsubscript{LOOM}})\\
                              1 & \quad \text{ if } z(t)=0\text{ and}\\
                                & \quad\phantom{W}\|\dot{\phi}_{GRM}(t)\| < \text{T\textsubscript{GRM} and}\\
                                & \quad\phantom{W} \omega_{LOOM}(t) < \text{T\textsubscript{LOOM} and }\\
                                & \quad\phantom{W} u(t) < P_{01} \\
                           z(t) & \quad\phantom{W} \text{ otherwise}
           \end{array} \right.

    \end{aligned}$$

Trajectory:

:   Call $\vec{v}$ the direction vector, i.e. $\|
       \vec{v} \| =1$. With a slight abuse of notation use $\vec{v}$ also for the angle of $\vec{v}$, i.e. write $\vec{v} \sim G( \cdot \,; \,
       \mu, \sigma)$ to indicate that the angle of $\vec{v}$ is drawn from a given Gaussian density mod $2 \pi$. $$\begin{aligned}
         x(0) &\sim  U(A)\footnote{}\\
         v(0) &\sim U(v_{\min}, v_{\max})\\
         \vec{v}(0) & \sim U(0, 2 \pi)\\
         x(t+ 1) &= x(t) + v(t)\vec{v}(t) \Delta t\\
         v(t+ 1) &= z(t) v\\
         \vec{v}(t+1) &\sim \left\{
           \begin{array}{l}
             G(\cdot \,; \, \vec{v}(t), \sigma_{\vec{v}}(t))\\
               \quad\text{if $z(t)-z(t+1) = 1$}\\
               \vec{v}(t) \text{otherwise}
           \end{array}
         \right.

    \end{aligned}$$

Direction change standard deviation:

:   $$\begin{aligned}
         \sigma_{\vec{v}}(0) & = & 0\\ \sigma_{\vec{v}}(t+1) & = &
         \lambda_\sigma \sigma_{\vec{v}}(t) + \delta_\sigma (1 - z(t))z(t-1) \\

    \end{aligned}$$

# DISCUSSION

Our results show that GRM constitutes a good cue for collision avoidance. Extending the analysis of , we showed both mathematically and in simulations that increasing the CVA to a non-zero value improves collision avoidance. In particular, it allows the detection of both stationary and moving objects on a collision course. In this respect, GRM can be viewed as a computationally efficient way to connect looming detection with regressive motion detection. We introduced safety and mobility as collision avoidance performance metrics and used looming as a reference point to show that GRM is a better cue for collision avoidance among conspecifics.

Whereas perception of regressive motion appears to influence the behavior of the fruit fly, the neural circuitry participating in this perception-action loop remains unknown. and showed that looming-sensitive neurons in *Drosophila* participate in a neural pathway that mediates escape behavior. However, escape is not the same as collision avoidance. We showed that looming might not be a practical collision-avoiding solution for groups of interacting animals as it can overly impede the mobility of a group. We argue that further research into the neural circuitry of GRM-based action in animals is an important future direction.

One reason why regerssive motion has remained a relatively obscure phenomenon might be the often overlooked difference between static and dynamic environments when testing collision avoidance algorithms. For example, constructed a robot guided by responses mimicking those of the locust looming-detection neurons. However, the robot's collision avoidance was only tested in an environment consisting of stationary obstacles. We have shown that in a more interactive environment, GRM has significant advantages over looming.

From the point of view of control and robotics, we studied what call *reciprocal collision avoidance*: avoiding crashes within a population of simple non-communicating robots all of which implement the same movement protocol. derive their elegant algorithm from first principles. Ours is biologically inspired and relies on a very simple visual cue that can be computed from the optical flow using Reichardt-detector-like circuits . While our simulations used agents that can immediately measure the position and velocity of any point on their retina, optic-flow-based motion detection has been successfully implemented on minimalistic hardware . We see GRM's largest potential in control of small swarming or flocking vehicles . More sophisticated algorithms for robot navigation are often based on Probabilistic Roadmaps  or Rapidly Exploring Random Trees . Such approaches allow a robot to avoid obstacles with lower failure rates than GRM, but require large computational power. Importance of low-complexity collision avoidance grows as fields such as drone flight control are rapidly developing . Alternatively, GRM could serve as a basis for a fallback collision prevention system for larger vehicles operating in highly dynamic environments.

# Acknowledgements

This work was supported by National Science Foundation grants 0914783 and 1216045, NASA Stennis grant NAS7.03001 and ONR MURI grant N00014-10-1-0933.

# MATHEMATICAL APPENDIX

**Proposition 2** *Let the relative position and velocity of the observed object be $\mathbf{x} \text{ and } \mathbf{v}$ respectively. Then $\dot{\phi} = \frac{1}{\|\mathbf{x}\|^2}\langle \mathbf{v}^\perp,
  \mathbf{x} \rangle$. In particular, angular velocity scales as one over distance squared.*

*Proof.* We wish to derive the angular velocity of a point in relative motion projecting onto an observer. Place the center of projection at the origin, and a particle moving with constant velocity $\mathbf{v}=(u,v)$ at position $\mathbf{x_0}=(x_0,y_0)$ at time $0$, as shown in Fig. .

The position at time $t$ equals $\mathbf{x}=\mathbf{x_0}+t\mathbf{v}$, and the azimuth of the particle $\phi$ is such that $$\tan{\phi} = \frac{\mathbf{y}}{\mathbf{x}} = \frac{y_0+vt}{x_0+ut}.$$ Taking the time derivative on both sides gives $$\frac{1}{\cos^2{\phi}}\dot{\phi} = \frac{v(x_0+tu)-u(y_0+tv)}{(x_0+ut)^2},$$ and thus $$\begin{aligned}
    \dot{\phi} &= \frac{cos^2{\phi}(v(x_0+tu)-u(y_0+tv))}{(x_0+tu)^2}\\
    &= \frac{(x_0+tu)^2(v(x_0+tu)-u(y_0+tv))}{D^2(x_0+tu)^2},\label{eq:dotphi}

\end{aligned}$$ where $D$ is the distance of the particle from the origin. Equation follows from the relation $\cos{\phi}=\frac{x_0+tu}{D}$ (see Fig. ). Simplifying the RHS yields $$\begin{aligned}
    \dot{\phi} &= \frac{v(x_0+tu)-u(y_0+tv)}{D^2}\\
    &= \frac{1}{\|\mathbf{x}\|^2}\langle \mathbf{v}^\perp, \mathbf{x} \rangle.\label{eq:angvel_derivation}

\end{aligned}$$ Since the position $\mathbf{x}$ moves on a line perpendicular to $\mathbf{v}^\perp$ we expect the dot product to be constant, and indeed it equals $vx_0-uy_0$. Hence angular velocity decays as one over distance squared. ◻

# MOVIE CAPTIONS APPENDIX

Available online are supplementary Movies S1-S4. The movies show full trials of our simulations for chosen parameter settings. In each movie, ten fly-like vehicles are visible, colored arbitrarily to make tracking the vehicles easy. Whenever two vehicles collide, their body size is temporarily increased. Whenever a vehicle stops, it is surrounded by a colored circle. A red circle means the stop is a False Positive. A green circle indicates a True Positive. In addition, a line segment is drawn from the stopping vehicle to the one (or more) causes of its stop.

[^1]: We will call 'collision' an event in which one moving agent comes into physical contact with another, or with a stationary obstacle. In a biological context, collisions can be harmful (as when a prey is caught by a predator, or a pathogen passes from one animal to another), or beneficial (as when one ant exchanges chemical information with a nest mate, or a predator succeeds in capturing a prey).

[^2]: The agents were also equipped with a looming detector, which in the simulations described in this section was set to be extremely insensitve.


abstract: The linear conductance spectrum of a metallic graphene junction formed by interconnecting two gapless graphene nanoribbons is calculated. A strong conductance suppression appears in the vicinity of the Dirac point. We found that such a conductance suppression arises from the antiresonance effect associated with the edge state localized at the zigzag-edged shoulder of the junction. The conductance valley due to the antiresonance is rather robust in the presence of the impurity and vacancy scattering. And the center of the conductance valley can be readily tuned by an electric field exerted on the wider nanoribbon.
address: Department of physics, Jilin University, Changchun 130023, China
author: Haidong Li; Lin Wang; Yisong Zheng
date: 2024-09-29
title: Suppressed conductance in a metallic graphene nano-junction

Graphene, an atomically thin layer of graphite, is regarded as a perspective base for the post-silicon electronics since its first experimental realization. Motivated by possible device applications, the electronic and transport properties of various graphene nano-structures have been studied both experimentally and theoretically, Among these structures, graphene nanoribbon(GNR) is the basic element to carry the current. Band structure calculations indicate that the zigzag-edged graphene nanoribbons are always metallic while the armchair-edged ones are either metallic or semiconducting, depending on the width of the armchair GNR. A graphene junction can be formed by interconnecting two semi-infinite GNRs with different widths. In such a graphene nanostructure, a traveling carrier is scattered by the junction interface, which causes a finite junction conductance. Recently, the conductance spectrum of a graphene metal-semiconductor junctions has been studied in details. It has been found that the presence of the lattice vacancy can efficiently enhance the junction conductance, because that a vacancy makes the coupling between the electron states of the two GNRs at the junction interface stronger.

Apart from the metal-semiconductor junction, a metallic graphene junction can be constructed by interconnecting two metallic armchair GNRs with different widths. In this letter, we investigate the conductance spectrum of such a metallic graphene junction, the lattice structure of which is depicted in Fig.(a). Unlike a metal-semiconductor junction, the gapless subband structures of the two metallic armchair nanoribbons shown in Fig.(b) indicates that electronic transmission through the metallic junction with an arbitrary energy is formally allowed. In particular, the electronic transmission via the linear subbands is reflectionless even in the presence of a long-range scattering potential, due to pseudospin conservation. Therefore, a plausible anticipation is that the electronic transmission probability near the Dirac point should be close to unity. However, our calculation gives the opposite result: the junction conductance at the Dirac point is equal to zero, and a conductance valley appears around this point. This means a strong conductance suppression in the vicinity of the Dirac point. Further analysis demonstrates that the conductance suppression arises from the antiresonance due to the existence of an edge state localized in the shoulder region of the junction.

To calculate the linear conductance of the metallic graphene junction $\mathcal{G}(E)$ as a function of the incident electron energy $E$, we adopt the Landauer-Bütikker formula in the discrete lattice representation. It gives $\mathcal{G}(E)={2e^2\over h}\text{Tr}(\Gamma_1G\Gamma_2G^\dag)$, where $G=[E+i0^+-H_d-\Sigma_1-\Sigma_2]^{-1}$ is a retarded Green function, and $H_d$ is the tight-binding Hamiltonian of the device region in the nearest-neighbor approximation. The contributions of the two semi-infinite leads are incorporated by the two self-energy terms $\Sigma_{1(2)}$ which are associated with the coupling functions $\Gamma_{1(2)}$ by $\Gamma_{1(2)}=i[\Sigma_{1(2)}-\Sigma_{1(2)}^\dag]$. The two self-energy terms are obviously the key quantities for calculating the conductance, which can be evaluated by the recursive method. In what follows we use the lattice constant $a$ and the hopping energy $t$ between the nearest neighbor atoms as the units of the length and energy.

The calculated conductance spectra($\mathcal{G}$ vs $E$) for some typical square junctions($\theta=90^\circ$) are shown in Fig.2(a). All these conductance spectra exhibit the staircase-like structures, which can be readily explained by matching of the subband structures of the two component ribbons as shown in Fig.(b). Another point to note in Fig.2(a) is that the conductance spectrum shows a notable suppression in the vicinity of the Dirac point when the difference of the widths of the two GNRs is larger than 3. In particular, a zero conductance occurs at the Dirac point. However, from Fig.1(b) we can readily find that the linear subbands of the two ribbons always match each other to provide an electron transmission mode, and hence a nonzero conductance at the Dirac point, in contradiction with the calculated zero conductance. As shown in Fig.2(b), we find that the conductance suppression near the Dirac point is tightly associated with the zigzag edge of the shoulder of the junction. On the contrary, when the edge of shoulder is of an armchair type($\theta=120^\circ$), the conductance spectrum no longer shows any suppression. From such a result we infer that the nature of the conductance suppression is the antiresonance effect, the detail of which is as follows. In the vicinity of the Dirac point, only the linear subbands are relevant to the electron transmission. Hence the two metallic GNRs can be viewed as single mode quantum wires coupling to each other directly. However, the zigzag-edged shoulder induces a localized edge state with the eigen-energy equal to the Dirac energy. Such a localized state couples to the linear subband of the wider GNR. Consequently, when the electronic transport is limited to the vicinity of the Dirac point, the graphene junction is equivalent to the T-shaped quantum dot structure as shown in Fig.2(c). The linear conductance of such a model has been extensively studied and can be expressed in terms of the model parameters $$\mathcal {G}(E)={2e^2\over
h}\frac{2\xi}{(1+\xi)^2}\frac{[E-\varepsilon]^2}
{(E-\varepsilon)^2+[\frac{\Gamma}{2(1+\xi)}]^2},\label{CF}$$ where $\Gamma=2\pi\rho_2 \tau^{2}$ and $\xi =\pi^2\rho_1\rho_2v^2$ with $\rho_{1(2)}$ being the electron density of the states in two leads. This expression presents a zero conductance at the quantum dot level $\varepsilon$, which is called the antiresonance effect. The antiresonance is in fact a result of quantum interference. The lateral quantum dot introduces new Feynman paths with a phase shift $\pi$. As a result, the destructive quantum interference occurs among electron Feynman paths. In the metallic graphene junction, the edge state attached to the zigzag-edged shoulder of the junction plays a role of laterally coupled quantum dot, which results in the antiresonance at the Dirac point.

When the widths of the two GNRs are fixed, the geometry of the junction can be changed by shifting downwards the narrower ribbon. In Fig.2(d) the conductance spectra are compared for differently shaped junctions. We can see that the width of conductance valley depends on the junction shape sensitively. To be more specific, in terms of the valley width the spectra are classified into three groups, each of which appears periodically whenever the narrower ribbon shifts downwards by three multiples of the lattice constant. This phenomenon can be explained with the help of the above quantum dot model. The parameter $v$ in the quantum dot model is a relevant quantity to the width of the antiresonance valley. From Fig.(a) we can see that $v$ is proportional to the product of the electron probability amplitudes of the A and B atoms interconnecting directly at the junction interface. And the probability amplitudes can be obtained by solving the Dirac equation . In such a way we work out the following relation $v\propto \sum_{j\in odd}\sin (2j\pi/3)\sin[2(j+n_0)\pi/3]$, where $j$ is any odd number within the range from 1 to $N_1$, and $n_0$ denotes the displacement of the narrower ribbon with respect to the upper edge of the wider ribbon. From this relation we can readily understand the periodic feature of the valley width of the conductance spectra shown in Fig.2(d).

The antiresonance picture of the conductance suppression in the metallic graphene junction is further demonstrated by the calculated spectra of the local density of states (LDOS) at some lattice points near the junction interface. From Fig.3 we can see that only for the lattice points at the zigzag edge of the shoulder of junction($\theta=90^\circ$ and $150^\circ$), the LDOS spectrum exhibits a very sharp peak at the Dirac point. This indicates the existence of a localized state in the shoulder region of the junction. An exception occurs for the junction with a very short shoulder($N_1=20$ and $N_2=23$ in Fig.3(a)). The LDOS spectrum of the lattice point at the zigzag-edged shoulder does not show a notable peak. This result is consistent with the previous work, which argued that a zigzag edge of the width smaller than three lattice constant can not induce any localized edged state.

We now proceed on to discuss the influence of the possible scatterers in an actual graphene junction on the antiresonance valley. At first, we consider the individual effect of an impurity appearing at distinct positions. In the numerical calculation, an impurity is simulated by the deviation of the on-site energy of an specific lattice point where the impurity appears. The calculated conductance spectra with an individual impurity positioned at different lattice points are compared in Fig.4(a). We can see that only when the impurity appears in the shoulder region, the variation of the conductance spectrum is notable. But the effect of such an impurity is not to destroy the antiresonance at the Dirac point. Instead it causes another conductance zero near the Dirac point. This result implies that the edge state is intricately affected by an impurity in the shoulder region. On the contrary, the impurity at other lattice points can not modify the antiresonance valley notably since it is irrelevant to the edge state. Fig.4(b) shows the conductance spectrum in the presence of many impurities distributed randomly in the device region with fluctuating strengths. We can see that the antiresonance is rather robust even if the impurity concentration and strength are nontrivially large. Fig.4(c) shows the effect of vacancies positioned uniformly at the zigzag edge of the shoulder of a square junction. A vacancy is simulated in the numerical calculation by simply cutting off a carbon atom of type A at the edge of the shoulder. If we use the notations $M_1$ and $M_2$ to denote the atom number of type A belonging to the edge of the shoulder and the number of the vacancy in this edge respectively, the ratio $r=M_2/M_1$ can be viewed as the concentration of vacancy. From Fig.4(c) we can see only when the concentration of the vacancies is about $r=1/3$, the zero conductance at the Dirac point can be completely eliminated. Higher or lower concentrations of vacancies can not destroy the conductance valley around the Dirac point. Our calculation also indicates that such a conclusion is independent of the size of the shoulder. This result implies the complicated effect of the edge defect on the edge state. Such an interesting topic is left for our study in the future. Finally, we can apply a step-like potential in the right hand side of the junction to tune the position of antiresonance, which can be simulated by shifting the on-site energy of all the lattice points of the wider GNR. The result is shown in Fig.4(d). We can see that the step-like potential simply shifts the antiresonance point, without drastically altering the lineshape of the conductance spectrum. What is noteworthy is that the conductance at the Dirac point can be easily tuned from zero to unity by a step-like potential. This suggests a possible device application of a nanoswitch based on such a metallic graphene junction.

This work was financially supported by the National Nature Science Foundation of China under Grant NNSFC10774055.