text:
# §1 Introduction

## §1 Introduction

Elliptic units are systems of units in abelian extensions of imaginary quadratic fields which are defined by evaluating some functions on CM elliptic curves at torsion points. These units play an essential role in Iwasawa theory of imaginary quadratic fields (\[dS\],\[Ru\] etc). In this paper we present certain systems of $`K_2`$ of function fields of CM elliptic curves. These systems are considered by Scholl (\[Sch\]) in a more general setting. We show that these systems have two similar properties to elliptic units.

First, in $`\mathrm{\S }2`$ we demonstrate that we can apply the theory of Coleman power series for $`K_2`$ of 2–dimensional local fields constructed by Fukaya (\[Fu\]) to our system. We show that the “Coleman power series” of our systems, which are elements in $`K_2`$ of some power series rings, are explicitly described. In $`\mathrm{\S }3`$ we modify our systems so that they are in the kernel of tame symbol maps. Using the result of $`\mathrm{\S }2`$ we will give a sufficient condition for the images of those systems to be non-zero under Chern class maps to Galois cohomology groups.

Second, in the appendix we show that our systems form Euler systems similar to those of elliptic units.

By \[Sch\] one can see that systems like ours also exist in higher $`K`$-groups of function fields of products of elliptic curves. This fact seems to suggest (to this author) that there is a further generalization of the theory of elliptic units and Coleman power series to higher $`K`$-groups.

The following question is asked by the referee:

Is there interpretation of the function which appears in corollary 3.4 in terms of L-function?

### Acknowledgment

The author would like to express his gratitude to Prof. K. Kato for valuable suggestions and to Dr.T. Fukaya for showing her paper at an early stage. He is also grateful to the referee for valuable comments on the first version of this paper. This work was carried out while the author stayed at Department of Mathematics, the University of Chicago. He is thankful to the people at Chicago for their hospitality.

## Notations

$`K`$: an imaginary quadratic field, 4

$`E`$: a CM elliptic curve, 4

$`F`$: an abelian extension of $`K`$, 4

$`\psi `$: the Hecke character of $`F`$ associated to $`E`$, 5

$`\varphi `$: a Hecke character of $`K`$ such that $`\psi =\varphi N_{F/K}`$, 5

$`𝔣`$: conductor of $`E`$, 5

$`p`$: a rational prime, 5

$`𝔭`$: a prime factor of $`p`$ in $`𝒪_K`$, 5

$`\pi `$: $`\varphi (𝔭)`$, 5

$`a`$: an integer relatively prime to $`6p𝔣`$, 5

$`g_a`$: a rational function on $`E`$, 5

$`v_n`$: a generator of $`E[\pi ^n]`$, 5

$`T_y`$: translation by a point $`y`$, 5

$`F_n`$: $`F(E[\pi ^n])`$, 5

$`\alpha _n`$: an element of $`K_2(F_n(E))`$, 5

$`W`$: the Witt ring of $`\overline{𝔽_p}`$, 6

$`𝔓`$: a prime of $`F`$ over $`p`$, 6

$`\widehat{}_W`$: the formal completion of a projective smooth model of $`E`$, 6

$`\theta `$: an isomorphism $`\widehat{𝔾_m}_W\widehat{}_W`$, 6

$`\mathrm{\Phi }`$: Frobenius of $`W`$, 6

$`t`$: a parameter on $`\widehat{}_W`$, 6

$`s`$: the standard parameter on $`\widehat{𝔾_m}_W`$, 6

$`(\zeta _n)_{n1}`$: a basis of $`_p(1)`$, 6

$`H`$: $`p`$-adic completion of $`W((s))`$, 6

$`H_n`$: $`H(s_n,\zeta _n)`$, 6

$`S`$: $`𝒪_H[[ϵ1]]`$, 6

$`S^{}`$: $`S[(ϵ1)^1]`$, 6

$`\nu ,\sigma ,\phi `$: endomorphisms of $`S`$, 6

$`h`$: an isomorphism $`K_2(S^{})_{}^{}{}_{}{}^{Nr_\phi =1}\underset{}{\mathrm{lim}}_nK_2(H_n)^{}`$, 7

$`G_k`$: $`Gal(k^{sep}/k)`$ for a field $`k`$, 7

$`H^i(k,M)`$: $`H^i(G_k,M)`$ for a $`G_k`$ module $`M`$, 7

$`\lambda ,`$: homomorphisms $`\underset{}{\mathrm{lim}}_nK_2(H_n)^{}\mathrm{\Omega }_{H_1}^1`$, 7

$`q`$ : $`s+1`$, 8

$`\beta _n`$: $`\theta _{}^{\mathrm{\Phi }^n}{}_{}{}^{}\alpha _n`$, 8

$`D`$: $`q\frac{d}{ds}`$, 9

$`F_n^{}`$: $`F(E[a\pi ^n])`$, 9

$`\alpha _n^\gamma `$: $`T_\gamma ^{}\alpha _n`$, 9

$`\alpha _n^{}`$: an element of $`K_2(F_n(E))`$, 9

$`\beta _n^{}`$: $`\theta _{}^{\mathrm{\Phi }^{n+1}}{}_{}{}^{}\alpha _n^{}`$, 10

$`g_a^\gamma `$: $`T_\gamma ^{}g_a`$, 11

## §2 norm-compatible systems in $`K_2`$

A general reference for this section is Chapter II of \[dS\]. We fix an imaginary quadratic field $`K`$. Let $`O_K`$ denote the ring of integers of $`K`$. Let F be an abelian extension of $`K`$ and let $`E`$ be an elliptic curve over $`F`$ which satisfies the following conditions:

(i) $`E`$ has complex multiplication by $`O_K`$, and

(ii) $`F(E_{tors})`$ is an abelian extension of $`K`$.

Here $`F(E_{tors})`$ denotes the field generated by coordinates of torsion points of a model of $`E`$ over $`F`$. If $`\psi `$ is the Hecke character of $`F`$ associated to $`E`$, then the condition (ii) is equivalent to the existence of a Hecke character $`\varphi `$ of $`K`$ such that $`\psi =\varphi N_{F/K}.`$

Let $`𝔣`$ be the conductor of $`E`$ and let $`p>3`$ be a rational prime which is relatively prime to $`𝔣`$ and which splits completely in $`F`$. Let $`p=𝔭\overline{𝔭}`$ be the decomposition of $`p`$ in $`O_K`$. We write $`\pi =\varphi (𝔭)`$. Note that $`𝔭=(\pi ).`$ We assume that $`Gal(F(E[p^n]))/F(O_K/(p^n))^{}`$ for all $`n1`$.

Let $`a`$ be an integer such that $`(a,\mathrm{\hspace{0.17em}\hspace{0.17em}6}p𝔣)=1`$. Then there is a function $`g_aF(E)`$, which is called the Kato-Siegel function in \[Sch\]. We now state the properties of $`g_a`$ which will be used in the subsequent discussions.

(a) div$`g_a=a^2(0)E[a]`$.

(b) For any isogeny $`j:EE`$ defined over $`\overline{F}`$ such that $`(\text{deg}j,a)=1`$, $`j_{}g_a=g_a`$.

(c) $`[a]_{}g_a=1`$.

(d) $`[1]^{}g_a=g_a`$. (This is a special case of (b)).

Here for an element $`u`$ of $`O_K`$, $`[u]`$ stands for an element of End($`E`$), which is justified by the condition (i). In our situation, the property (b) is equivalent to the following distribution relation

$$\underset{j(y)=x}{}g_a(y)=g_a(x).$$

We fix a system $`(v_n)_{n1}`$ of torsion points of $`E(\overline{F})`$ such that for each $`n2`$, $`v_n`$ is a generator of $`E[\pi ^n]`$ and that $`[\pi ]v_n=v_{n1}`$.

For a point $`y`$ on $`E`$ we denote by $`T_y:EE`$ the translation by $`y`$. For $`n1`$ consider the element

$$\alpha _n:=\{g_a,T_{v_n}^{}g_a\}K_2(F_n(E)).$$

Here $`F_n:=F(E[\pi ^n])`$ and $`F_n(E)`$ is the function field of $`E_{F_n}=E\times _FF_n`$. This element is considered by Scholl(\[Sch\], Appendix) in a more general setting. The following lemma claims that the system $`(\alpha _n)`$ forms an norm-compatible system.

###### Lemma 2.1

For $`n2`$, the equality

$$[\pi ]_{}N_{F_n(E)/F_{n1}(E)}\alpha _n=\alpha _{n1}$$

holds.

###### Demonstration Proof

This is another way of describing (EN1) in loc.cit. We will give a proof by calculation of symbols.

$$N_{F_n(E)/F_{n1}(E)}\alpha _n=\{g_a,\underset{vE[\pi ]}{}T_{v_n+v}^{}g_a\}=\{g_a,[\pi ]^{}(T_{v_{n1}}^{}g_a)\}$$

where the second equality holds by the distribution relation mentioned above, and

$$[\pi ]_{}\{g_a,[\pi ]^{}(T_{v_{n1}}^{}g_a)\}=\{[\pi ]_{}g_a,T_{v_{n1}}^{}g_a\}=\{g_a,T_{v_{n1}}^{}g_a\}.$$

We will apply the theory of Coleman power series for $`K_2`$ constructed by Fukaya \[Fu\] to this system. As at present her theory treats only the multiplicative group, we need to construct an isomorphism from $`E`$ to the multiplicative group. For this we will take the formal completion of $`E`$ along the $`0`$-section.

We will follow the argument of §4.3,Chapter II of \[dS\] rather closely. Let $`F^{}=F(E[\overline{\pi }^{\mathrm{}}])`$ and let $`R`$ be the completion of $`O_F^{}_{O_K}O_{K}^{}{}_{𝔭}{}^{}`$. Since the prime $`𝔭`$ finitely decomposes in $`F^{}`$, $`R`$ is a direct sum of local rings. Fix a direct summand of $`R`$ and denote it by $`W`$. Let $`𝔓`$ be the prime of $`F`$ over $`𝔭`$ which corresponds to $`W`$ and let $``$ be a smooth projective model of $`E`$ over $`O_{F}^{}{}_{𝔓}{}^{}=_p`$. Let $`\widehat{}_W`$ be the formal completion of $`\times __pW`$ along the $`0`$-section. Then by Proposition 4.3 in Chap.2 of \[dS\], there is an isomorphism

$$\theta :\widehat{𝔾_m}_W\widehat{}_W$$

which satisfies the equality $`[\pi ]\theta =\theta ^\mathrm{\Phi }[p]`$. Here $`\mathrm{\Phi }`$ is the Frobenius automorphism of $`W`$. We fix a parameter $`t`$ of $`\widehat{}_W`$ and let $`s`$ be the standard parameter of $`\widehat{𝔾_m}_W`$.

For $`n1`$, define $`\zeta _n\overline{\text{Frac}(W)}`$ by the equality $`t(v_n)=\theta ^{\mathrm{\Phi }^n}(\zeta _n1)`$. Then we see that the system $`(\zeta _n)_{n1}`$ is a basis of $`_p(1).`$

Now we will review the theory of Fukaya in the form we will use in the subsequent discussions.

Let $`H`$ be the $`p`$-adic completion of $`W((s))`$, the function field of $`\widehat{𝔾_m}_W`$. Consider a tower of field extensions $`(H_n)_{n1}`$ where $`H_n:=H(s_n,\zeta _n)`$ and $`(s_{n1}+1)=(s_n+1)^p`$ for $`n1`$. $`s_0=s`$ by convention.

We set $`S:=O_H[[ϵ1]]`$ to be the ring of formal power series and write $`S^{}=S[(ϵ1)^1]`$. Consider the ring homomorphisms $`\nu ,\sigma `$ and $`\phi :S^{}S^{}`$ defined as follows:

$`\nu `$ is defined by $`\nu (ϵ)=ϵ^p`$ and $`\nu |_{O_H}=`$ the identity. The homomorphism $`\sigma `$ is defined by $`\sigma |_{O_H}=[p]^{}\mathrm{\Phi }`$ and $`\sigma (ϵ)=ϵ.`$ The homomorphism $`\phi `$ is defined by the composition: $`\phi =\nu \sigma `$. The following theorem is a part of Theorem 1.5 of \[Fu\]. It asserts the existence of Coleman power series:

###### Theorem 2.2

There is a canonical isomorphism

$$h:K_2(S^{})_{}^{}{}_{}{}^{Nr_\phi =1}\underset{}{\mathrm{lim}}_nK_2(H_n)^{}$$

induced by

$$h_n:S^{}H_n,f(ϵ)(\sigma ^nf)(\zeta _n).$$

Here $`Nr_\phi `$ is the norm map on $`K_2(S^{})`$ induced by the finite morphism $`\phi `$ and $`^{Nr_\phi =1}`$ means the fixed part of $``$ under $`Nr_\phi `$. The map $`\sigma ^n:HH_n`$ is the composition of $`id\mathrm{\Phi }^n:HH`$ and the inclusion given by $`ss_n`$. The symbol $``$ means certain completion.

Also there is a fact which corresponds to the explicit reciprocity law. It asserts the coincidence with $``$ sign of the two homomorphisms $`\lambda `$ and $``$ defined as follows.

For a field $`k`$, we write $`G_k`$ for $`Gal(k^{sep}/k)`$. For a $`G_k`$ module $`M`$, we write $`H^i(k,M)`$ for $`H^i(G_k,M).`$

$`\lambda :\underset{}{\mathrm{lim}}_nK_2(H_n)^{}`$ $`\stackrel{cl}{}\underset{}{\mathrm{lim}}_nH^2(H_n,(/p^n)(2))`$
$`\stackrel{(\zeta _n)^1}{}\underset{}{\mathrm{lim}}_nH^2(H_n,(/p^n)(1))`$
$`\stackrel{Cor_{H_n/H_1}}{}\underset{}{\mathrm{lim}}_nH^2(H_1,(/p^n)(1))`$
$`=H^2(H_1,_p(1))`$
$`H^2(H_1,_p(1))\mathrm{\Omega }_{H_1}^1.`$
The last isomorphism is due to Hyodo \[Hy\]. Here $`\mathrm{\Omega }_{H_1}^1`$ is defined to be $`\mathrm{\Omega }_{O_{H_1}}`$ where $`\mathrm{\Omega }_{O_L}^1:=\underset{𝑛}{\underset{}{\mathrm{lim}}}\mathrm{\Omega }_{O_{H_1}/}^1/p^n\mathrm{\Omega }_{O_{H_1}/}^1`$.

The homomorphism $``$ is defined using the map $`h`$ in Theorem 2.2:

$`:\underset{}{\mathrm{lim}}_nK_2(H_n)^{}`$ $`\stackrel{h^1}{}K_2(S^{})^{Nr_\phi =1}`$
$`\stackrel{d\mathrm{log}}{}\mathrm{\Omega }_S^2(\mathrm{log})={\displaystyle \frac{S}{ϵ1}}d\mathrm{log}(q)d\mathrm{log}(ϵ)`$
$`\stackrel{(d\mathrm{log}(ϵ))^1}{}{\displaystyle \frac{S}{ϵ1}}d\mathrm{log}(q)`$
$`\mathrm{\Omega }_{H_1}^1`$
where $`q=s+1`$ and the last map is defined by

$$f(ϵ)d\mathrm{log}(q)\frac{1}{p}(\sigma ^1f)(\zeta _1)d\mathrm{log}(q).$$

The modules $`\mathrm{\Omega }_S^1(\mathrm{log})`$ and $`\mathrm{\Omega }_S^2(\mathrm{log})`$ are defined as follows:

Let

$$\mathrm{\Omega }_{S/}^1(\mathrm{log}):=(\mathrm{\Omega }_{S/}^1S_{}S_{}^{}{}_{}{}^{})/N$$

where $`N`$ is the $`S`$ submodule generated by elements $`(da,aa)`$ for $`aSS_{}^{}{}_{}{}^{}`$. The class $`(0,\mathrm{\hspace{0.17em}1}a)`$ is denoted by $`d\mathrm{log}(a).`$ The module $`\mathrm{\Omega }_S^2(\mathrm{log})`$ is defined by

$$\mathrm{\Omega }_S^2(\mathrm{log})=\underset{}{\mathrm{lim}}_n^2\mathrm{\Omega }_{S/}^1(\mathrm{log})/p^n^2\mathrm{\Omega }_{S/}^1(\mathrm{log}).$$

By proposition (19.3.5) of \[Gr\], $`S`$ is formally smooth over $``$, so by corollary (20.4.11) of loc.cit. $`\mathrm{\Omega }_{S/}^1(\mathrm{log})/p^n\mathrm{\Omega }_{S/}^1(\mathrm{log})`$ is a free $`S/p^nS`$ module on $`d\mathrm{log}(ϵ1)`$ and $`d\mathrm{log}(q)`$. Hence $`\mathrm{\Omega }_S^2(\mathrm{log})`$ is a free $`S`$ module on $`d\mathrm{log}(q)d\mathrm{log}(ϵ1).`$

The map $`d\mathrm{log}:K_2(S^{})\mathrm{\Omega }_S^2(\mathrm{log})`$ is described as $`\{f,g\}d\mathrm{log}(f)d\mathrm{log}(g)`$.

###### Theorem 2.3 (Theorem 2.4 in \cite{Fu})

$$\lambda =.$$

For $`n1`$, set $`\beta _n:=\theta _{}^{\mathrm{\Phi }^n}{}_{}{}^{}\alpha _nK_2(H_n).`$ Since $`[p]_{}\theta ^{}=\theta _{}^{\mathrm{\Phi }}{}_{}{}^{}[\pi ]_{}`$ it follows that $`(\beta _n)_n`$ is an element of $`\underset{}{\mathrm{lim}}_nK_2(H_n)`$ where projective limit is taken with respect to norm maps.

The following is our main result:

###### Theorem 2.4

We have the equality

$$h^1((\beta _n)_n)=B:=\{\theta ^{}g_a(s),\theta ^{}g_a(s[+]_{\widehat{𝔾_m}_W}ϵ1)\}.$$

###### Demonstration Proof

First of all, note that since $`g_a(t)t^{1a^2}W[[t]]^{}`$, $`\theta ^{}g_a(s)s^{1a^2}W[[s]]^{}O_H^{}`$. Since $`\theta ^{}g_a(s[+]_{\widehat{𝔾_m}_W}ϵ1)`$ is of the form $`\theta ^{}g_a(s)+(ϵ1)f(s,ϵ1)`$ with $`f(s,ϵ1)S`$, it belongs to $`S^{}`$. It can also be checked directly that $`h_n(B)=\beta _n`$ for $`n1.`$ What remains is to show that $`Nr_\phi B=B.`$

$`Nr_\phi B=`$ $`Nr_\sigma Nr_\nu B`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,{\displaystyle \underset{\zeta ^p=1}{}}\theta ^{}g_a(s[+]_{\widehat{𝔾_m}_W}\zeta ϵ1)\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,{\displaystyle \underset{\zeta ^p=1}{}}\theta ^{}g_a(s[+]_{\widehat{𝔾_m}_W}ϵ1[+]_{\widehat{𝔾_m}_W}\zeta 1)\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,{\displaystyle \underset{\zeta ^p=1}{}}g_a(\theta (s[+]_{\widehat{𝔾_m}_W}ϵ1)[+]_{\widehat{}_W}\theta (\zeta 1))\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,{\displaystyle \underset{vE[\pi ]}{}}g_a(\theta (s[+]_{\widehat{𝔾_m}_W}ϵ1)[+]_{\widehat{}_W}v)\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,g_a([\pi ](\theta (s[+]_{\widehat{𝔾_m}_W}ϵ1)))\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,g_a(\theta ^\mathrm{\Phi }([p](s[+]_{\widehat{𝔾_m}_W}ϵ1)))\}`$
$`=`$ $`Nr_\sigma \{\theta ^{}g_a,[p]^{}(\theta _{}^{\mathrm{\Phi }}{}_{}{}^{}g_a)(s[+]_{\widehat{𝔾_m}_W}ϵ1)\}.`$
Since $`\sigma =[p]^{}\mathrm{\Phi }`$, $`Nr_\sigma =\mathrm{\Phi }^1[p]_{}.`$ From this and the equality $`[p]_{}(\theta ^{}g_a)=\theta _{}^{\mathrm{\Phi }}{}_{}{}^{}[\pi ]_{}g_a=\theta _{}^{\mathrm{\Phi }}{}_{}{}^{}g_a`$ it follows that $`Nr_\phi B=B.`$

Let $`D`$ be a differential operator on $`H`$ defined by $`D=q\frac{d}{ds}`$. Then by straightforward calculation we obtain the following corollary.

###### Corollary 2.5

$$((\beta _n)_n)=\frac{1}{p}\mathrm{\Phi }^1(D\mathrm{log}(\theta ^{}g_a)(s)\times D\mathrm{log}(\theta ^{}g_a)(s[+]_{\widehat{𝔾_m}_W}\zeta _11))d\mathrm{log}q.$$

## §3. Elements in the kernel of tame symbol map

For $`n1,`$ let $`F_n^{}:=F(E[a\pi ^n])`$. For a point $`\gamma E[a]`$, let $`\alpha _n^\gamma :=T_\gamma ^{}\alpha _nK_2(F_n^{}(E))`$. We set

$$\alpha _n^{}:=N_{F_n^{}/F_n}[\pi ]_{}\underset{\gamma E[a]}{}\alpha _n^\gamma .$$

We can check that for $`n2`$,

$$[\pi ]_{}N_{F_n(E)/F_{n1}(E)}\alpha _n^{}=\alpha _{n1}^{}$$

in a similar way as in the case of $`\alpha _n`$. We have the following exact sequence (\[Qu\]):

$$\underset{xE(\overline{F})}{}K_2(\kappa _x)K_2(E_{F_n})\stackrel{𝑟}{}K_2(F_n(E))\stackrel{\underset{xE(\overline{F})}{}t_x}{}\underset{xE(\overline{F})}{}\kappa _x^{}\stackrel{𝑁}{}F_n^{}$$

where $`r`$ is the restriction, $`N`$ is the norm and for $`xE(\overline{F})`$, $`t_x`$ is defined by

$$\{f,g\}(1)^{ord_xford_xg}\frac{f^{ord_xg}}{g^{ord_xf}}(x).$$

By this exact sequence we see that $`[\pi ^n]^{}[\pi ^n]_{}\underset{\gamma E[a]}{}\alpha _n^\gamma =\underset{uE[a\pi ^n]}{}T_u^{}\alpha _n`$ is in the kernel of$`\underset{xE(\overline{F})}{}t_x`$, and since the map $`[\pi ^n]:EE`$ is etale, it follows that $`\underset{xE(\overline{F})}{}t_x([\pi ^{n1}]_{}\alpha _n^{})=0.`$

Hence there is an element $`𝒜K_2(E_{F_n})`$ such that $`r(𝒜)=[\pi ^{n1}]_{}\alpha _n^{}.`$ Since $`K_2`$ of number fields is torsion (\[Ga\]), $`𝒜`$ is well defined up torsion.

###### Lemma 3.1

Let $`V_p(E)(1):=T_p(E)(1)`$ and set $`F_{\mathrm{}}:=F(E[\pi ^{\mathrm{}}]).`$ Then

$$(V_p(E)(1)/T_p(E)(1))^{G_F_{\mathrm{}}}=0.$$

###### Demonstration Proof

There is a decomposition $`T_p(E)(1)=T_𝔭(E)(1)T_{\overline{𝔭}}(E)(1)`$, so we will treat each factor separately. The factor $`T_𝔭(E)(1)`$ is isomorphic to $`_p(1)`$ as a $`G_F_{\mathrm{}}`$ module. By corollary 1.7 in Chap.II of \[dS\], the fields $`F_{\mathrm{}}`$ and $`F(E[\overline{\pi }^{\mathrm{}}])`$ are linearly disjoint over $`F`$. Since $`_p(1)T_𝔭(E)T_{\overline{𝔭}}(E)`$, the fields $`F_{\mathrm{}}`$ and the field $`F(\zeta _p^{\mathrm{}})`$, generated by $`p`$-power roots of unity, are also linearly disjoint over $`F`$. Hence the assertion holds for this factor.

As for $`T_{\overline{𝔭}}(E)(1)`$, it is isomorphic to $`_p(2)`$ as a $`G_F_{\mathrm{}}`$ module. Since $`p>3`$, the assertion holds for this factor as well. ∎

By the lemma, $`H^1(F_n,T_p(E)(1))`$ has no torsion and we have a well defined class $`\stackrel{~}{cl}([\pi ^{n1}]_{}\alpha _n^{}):=cl(𝒜)H^1(F_n,T_p(E)(1)).`$

We set $`\beta _n^{}:=\theta _{}^{\mathrm{\Phi }^{n+1}}{}_{}{}^{}\alpha _n^{}K_2(H_n).`$ Since all the points in $`E[a]`$ are rational over Frac$`W`$ and $`\theta _{}^{\mathrm{\Phi }}{}_{}{}^{}[\pi ]_{}=[p]_{}\theta ^{}`$, we have the equality

$$\beta _n^{}=[F_0^{}:F][p]_{}\underset{\gamma E[a]}{}\theta _{}^{\mathrm{\Phi }^n}{}_{}{}^{}\alpha _n^\gamma .$$

One can check that $`(\beta _n^{})_n`$ form an element of $`\underset{}{\mathrm{lim}}_nK_2(H_n)`$ in a similar way to the case of $`\beta _n`$.

It can also be shown that $`\theta ^{}([\pi ^{n1}]_{}\alpha _n^{})=[p^{n1}]_{}\beta _n^{}`$. By (5.10) in \[Fu\], we have a commutative diagram

$$\begin{array}{ccccc}\underset{}{\mathrm{lim}}_nK_2(H_n)^{}& \stackrel{cl}{}& \underset{}{\mathrm{lim}}_nH^2(H_n,/p^n(2))& \stackrel{C_1}{}& H^2(H_1,_p(1))\\ \left[p^{n1}\right]_{}& & \text{Cor}_{H_n/H_1\left(\zeta _n\right)}& & \text{id}& & \\ \underset{}{\mathrm{lim}}_nK_2(H_1(\zeta _n))^{}& \stackrel{cl}{}& \underset{}{\mathrm{lim}}_nH^2(H_1(\zeta _n),/p^n(2))& \stackrel{C_2}{}& H^2(H_1,_p(1))\end{array}$$

Here $`C_1=\text{Cor}_{H_n/H_1}(\zeta _n)^1`$ and $`C_2=\text{Cor}_{H_1(\zeta _n)/H_1}(\zeta _n)^1`$. By this diagram it follows that if $`((\beta _n^{})_n)0`$, then $`cl([p^{n1}]_{}\beta _n^{})0`$ for sufficiently large $`n`$.

It implies that $`0cl([\pi ^{n1}]_{}\alpha _n^{})H^2(F_{n}^{}{}_{𝔓}{}^{}(E),/p^n(2))`$ where $`F_{n}^{}{}_{𝔓}{}^{}`$ is completion of $`F_n`$ by the prime over $`𝔓.`$

Consider the restriction of $`\stackrel{~}{cl}([\pi ^{n1}]_{}\alpha _n^{})`$ to $`H^1(F_{n}^{}{}_{𝔓}{}^{},T_p(E)(1))`$. We need a splitting $`H^1(F_{n}^{}{}_{𝔓}{}^{},T_p(E)(1))H^2(E_{F_{n}^{}{}_{𝔓}{}^{}},/p^n(2))`$ and the following lemma enables us to define such a splitting.

###### Lemma 3.2

$`\mathrm{}H^2(F_{n}^{}{}_{𝔓}{}^{},/p^n(2))=\mathrm{}H^0(F_{n}^{}{}_{𝔓}{}^{},/p^n(1))`$ is bounded with respect to $`n`$.

###### Demonstration Proof

Considering the Weil pairing $`E[\pi ^n]\times E[\overline{\pi }^n]\mu _{p^n}`$, we see that if $`F_{n}^{}{}_{𝔓}{}^{}F_𝔓(\mu _{p^m})`$ for a number $`mn`$, then $`F_{n}^{}{}_{𝔓}{}^{}F_𝔓(E[\overline{\pi }^m]).`$ Since $`𝔓`$ is finitely decomposed in $`F(E[\overline{\pi }^{\mathrm{}}])`$, the residue field of $`F_𝔓(E[\overline{\pi }^{\mathrm{}}])`$is $`\overline{𝔽_p}`$. As $`F_{n}^{}{}_{𝔓}{}^{}`$ is totally ramified over $`F_𝔓`$, it follows that $`m`$ is bounded. ∎

By this lemma we have a number $`M`$ which is a power of $`p`$ and a map

$$u_n:H^1(F_{n}^{}{}_{𝔓}{}^{},T_p(E)(1))H^2(E_{F_{n}^{}{}_{𝔓}{}^{}},/p^n(2))$$

for each $`n1`$ such that

$$u_n(\stackrel{~}{cl}([\pi ^{n1}]_{}\alpha _n^{})=Mcl([\pi ^{n1}]_{}\alpha _n^{}).$$

Noting that the action of the Frobenius $`\mathrm{\Phi }`$ on $`E[a]`$ is by $`\pi `$, one can show the following in a similar way to the proof of Theorem 2.4.

###### Theorem 3.3

For $`\gamma E[a]`$, set $`g_a^\gamma :=T_\gamma ^{}g_a.`$ Then we have the equality

$$h^1((\beta _n^{}))=B^{}:=[F_0^{}:F][p]_{}\underset{\gamma E[a]}{}\{\theta ^{}g_a^\gamma ,\theta ^{}g_a^\gamma (s[+]_{\widehat{𝔾_m}_W}ϵ1)\}.$$

Since $`d\mathrm{log}`$ and the pull back $`[p]^{}`$ commutes on $`S^{}`$, it also commutes on $`K_2(S^{})`$. As

$`[p]^{}B^{}=`$ $`[F_0^{}:F]{\displaystyle \underset{\zeta ^p=1}{}}{\displaystyle \underset{\gamma E[a]}{}}`$
$`\{\theta ^{}g_a^\gamma (s[+]_{\widehat{𝔾_m}_W}\zeta 1),\theta ^{}g_a^\gamma (s[+]_{\widehat{𝔾_m}_W}\zeta 1[+]_{\widehat{𝔾_m}_W}ϵ1)\}`$
we see that

$`[p]^{}((\beta _n^{}))`$
$`=`$ $`{\displaystyle \frac{[F_0^{}:F]}{p}}\mathrm{\Phi }^1({\displaystyle \underset{\zeta ^p=1}{}}{\displaystyle \underset{\gamma E[a]}{}}D\mathrm{log}\theta ^{}g_a^\gamma (s[+]_{\widehat{𝔾_m}_W}\zeta 1)`$
$`\times D\mathrm{log}\theta ^{}g_a^\gamma (s[+]_{\widehat{𝔾_m}_W}\zeta 1[+]_{\widehat{𝔾_m}_W}\zeta _11))d\mathrm{log}q.`$
We write the last term as $`\frac{[F_0^{}:F]}{p}\mathrm{\Phi }^1G(s)d\mathrm{log}q.`$

Let $`\lambda :\widehat{E}\widehat{𝔾}_a`$ be a logarithm defined over $`F`$ and let $`z`$ be the parameter of $`\widehat{𝔾}_a`$ such that $`z=\lambda (t)=t+\mathrm{}`$. We have the equality

$$\mathrm{\Omega }_p\frac{d}{dz}\lambda _{}\theta _{}=\lambda _{}\theta _{}D$$

where $`\mathrm{\Omega }_pW^{}`$ is such that $`t=\theta (s)=\mathrm{\Omega }_ps+\mathrm{}`$. Let $`\xi :E()/L`$ be the inverse of the analytic uniformization associated to $`dz`$. Here $`L`$ is the period lattice of $`dz`$. For a point $`vE[\pi ]`$, $`\frac{d}{dz}\lambda _{}g_a^\gamma (z[+]_{\widehat{𝔾}_a}\lambda (v))`$ is the Taylor expansion of $`\frac{d}{dz}\xi _{}g_a^\gamma (z+\xi (v))`$ at $`z=0.`$ Hence by looking at $`G(s)`$ we obtain the following corollary.

###### Corollary 3.4

If the function

$$\underset{vE[\pi ]}{}\underset{\gamma E[a]}{}\frac{d}{dz}\mathrm{log}(\xi _{}g_a^\gamma )(z+\xi (v))\times \frac{d}{dz}\mathrm{log}(\xi _{}g_a^\gamma )(z+\xi (v+v_1))$$

on $`/L`$ is nonzero, then $`\stackrel{~}{cl}([\pi ^{n1}]_{}\alpha _n^{})`$ is nonzero for sufficiently large $`n`$.

## Appendix: Euler Systems

Let $`𝒮`$ be the set of ideals of $`O_K`$ which are divisible only by primes which split completely in $`F`$ and which are relatively prime to $`𝔣a`$. For each $`𝔪𝒮`$, Let $`F_𝔪:=F(E[𝔪]).`$ Fix a system of torsion points $`(v_𝔪)_{𝔪𝒮}`$ such that

(1) $`v_𝔪`$ is a generator of $`E[𝔪]`$ for each $`𝔪𝒮.`$

(2) If $`𝔪_2=𝔪_1𝔫`$ with $`𝔪_2𝒮`$, then $`v_{𝔪_1}=\varphi (𝔫)v_{𝔪_2}`$.

We set $`\alpha _𝔪:=\{g_a,T_{v_𝔪}^{}g_a\}K_2(F_𝔪(E))`$ and set $`\alpha _𝔪^{}:=[\varphi (𝔪)]_{}\underset{\gamma E[a]}{}T_\gamma ^{}\alpha _𝔪.`$ Then it is seen that $`\alpha _𝔪^{}\mathrm{\Gamma }(E_{F_𝔪},𝒦_2)`$ in the similar way to the argument in the beginning of §3. So that there is a class $`\stackrel{~}{cl}(\alpha _𝔪^{})H^1(F_𝔪,T_p(E)(1)).`$ In the similar way to the proof of (EN1) in (\[Sch\], Appendix) one obtains the following result:

###### Proposition A.1

Let $`𝔪𝒮`$ and assume that $`𝔩|𝔪`$ for a prime $`𝔩`$. Then the following equality holds:

$$N_{F_𝔪/F_{𝔪/𝔩}}\alpha _𝔪^{}=\{\begin{array}{cc}& \alpha _{𝔪/𝔩}^{}𝔩^2|𝔪\hfill \\ & (1[\varphi (𝔩)]_{}Fr_𝔩^1)\alpha _{𝔪/𝔩}^{}𝔩^2𝔪.\hfill \end{array}$$

The classes $`\stackrel{~}{cl}(\alpha _𝔪^{})`$ satisfy the same relation as $`\alpha _𝔪^{}`$. This implies the classes $`\stackrel{~}{cl}(\alpha _𝔪^{})`$ form an Euler system.


format: text

text:
# Wide field interferometric imaging with single-mode fibers

## 1 Introduction

In the visible and near-infrared, single-mode optical fibers allow efficient transport of coherent light over large distances: the transmission losses inside the fiber are typically 1dB/km in the near-infrared. The throughput of a fiber interferometer can be very high, since the number of optical elements before the fiber is small, and there is no need for numerous beam steering mirrors. Moreover, the beam combination can be done using fiber optics couplers, and does not require high numbers of glass to air interfaces. A second interesting property of single mode optical fibers is the spatial filtering of the incoming wavefronts. Spatial filtering with single-mode optical fibers is superior to pinhole filtering (Keen et al. (2001)) and is used in interferometry to restore coherence of the telescope’s beams. Single-mode fibers only allow the propagation of a fundamental mode ($`LP_{01}`$). The loss of light coherence of the entrance pupil due to wavefront aberrations is then traded for a loss of coupling efficiency (Coude du Foresto et al. (2000)). The object visibility measurements can be corrected for this loss of coupling efficiency by recording in real time the amount of light coupled in each fiber (real-time self calibration of the interferometer), which greatly improves the accuracy of the object visibility measurement. Thanks to these two properties, fiber optics offer a reliable and simple solution to transport and combine coherent beams in interferometers with few apertures, and the recent progress in integrated optics will soon allow complex beam combiners to be built for larger interferometer arrays (Malbet et al (1999); Berger et al. (1999, 2001)).

Unfortunately, the restoration of coherence by spatial filtering comes at the expense of a very small field of view (FOV), which is the size of the diffraction spot of the largest telescope in the array. This is a serious limitation for the new generation of interferometers (Keck, VLTI, OHANA) which make use of large AO-equipped telescopes: the FOV is very small for such large telescopes. Moreover, these interferometers make it possible to observe fainter non-stellar sources such as AGNs and YSOs, for which the required FOV is often larger than for stellar observations. In this work, it is shown in §2 and §3 that spatial filtering reduces the field of view and that wavefront aberrations corrupt the object visibility measurements. This later effect is quantified and discussed in §4 and §5. The use of multiple fibers per telescope offers an attractive solution to these limitations (§6) and greatly improves the wide field imaging capabilities of space (§7) and ground-based (§8) interferometers.

## 2 Coupling efficiency and field of view

### 2.1 Coupling efficiency

Coupling of starlight into a single-mode fiber is achieved by positioning the fiber head in the focal plane and centering it on the PSF of the star. The coupling efficiency, which is the fraction of the starlight coupled into the fiber, is given by

$$\rho =\frac{\left|E_{fp}E_{01}^{}𝑑A\right|^2}{\left|E_{fp}\right|^2𝑑A\times \left|E_{01}\right|^2𝑑A},$$
(1)
where $`E_{fp}`$ is the incident electric field at the head of the fiber, in the focal plane, and $`E_{01}`$ is the electric field of the fundamental ($`LP_{01}`$) mode of propagation in the fiber. $`E_{fp}`$ is obtained by Fourier transform of the electric field in the entrance pupil, and $`E_{01}`$ is very well approximated by a Gaussian of constant phase across the width of the fiber. The coupling efficiency between the $`LP_{01}`$ mode and the Airy function is 0.82 (Shaklan & Roddier (1988); Ruilier (1998)). Therefore, under optimal conditions, no more than 82% of the incoming starlight can be coupled into the single-mode fiber. Fresnel reflection by the surface of the fiber head reduces this coupling efficiency to 78 % unless an anti-reflective coating is used on the fiber head.

### 2.2 Field of view of the fiber

The coupling efficiency decreases rapidly with the distance of the point source to the optical axis of the telescope. Using equation (1), Shaklan and Roddier (1988) have shown that the coupling efficiency, in the case of a circular aperture and a fundamental mode approximated by a Gaussian, is

$`\rho (\alpha ,d,f,\omega ,\lambda )={\displaystyle \frac{8}{\omega ^2}}e^{2\left(\frac{|\alpha |f}{\omega }\right)^2}`$ (2)
$`\times \left[{\displaystyle e^{\left(\frac{r}{\omega }\right)^2}I_0(2r|\alpha |f/\omega ^2)J_1(\pi dr/\lambda f)𝑑r}\right]^2.`$
Where $`f`$ is the telescope focal length, $`d`$ is the telescope diameter, $`\lambda `$ is the wavelength, $`\omega `$ is the $`1/e`$ width of the fundamental mode, $`\alpha `$ is the angular position of the point source relative to the optical axis, $`J_1`$ is the first order Bessel function and $`I_0`$ is the zero order modified Bessel function. We note that this expression is only a function of $`X_1=|\alpha |/(\lambda /d)`$ and $`X_2=(f\lambda )/(\omega d)`$ :

$`\rho (X_1,X_2)=8e^{2(X_1X_2)^2}`$ (3)
$`\times \left[{\displaystyle e^{r^2}I_0(2X_1X_2r)J_1(\frac{\pi r}{X_2})𝑑r}\right]^2`$
where $`r`$ is now dimensionless. $`\rho (0,X_2)`$ is maximum for $`X_2=1.401`$, and this value will be used in this work. We can then express $`\rho `$ as a function of $`X_1`$ only, which leads to :

$`\rho (\alpha )=8e^{3.923\left(\frac{|\alpha |}{\lambda /d}\right)^2}`$ (4)
$`\times \left[{\displaystyle e^{r^2}I_0(2.802\frac{|\alpha |}{\lambda /d}r)J_1(\frac{\pi r}{1.402})𝑑r}\right]^2.`$
The relation between the amplitude $`E`$ of the electric field coupled in the fiber and the luminosity $``$ of a point source which position is $`\alpha `$ is then

$$E(\alpha )=\sqrt{}\frac{d\sqrt{\pi }}{2}\sqrt{\rho (\alpha )}.$$
(5)
Figure 1 shows the coupling efficiency as a function of $`|\alpha |`$, in units of $`\lambda /d`$, for an unobstructed circular pupil and a circular pupil with an central obscuration of diameter $`0.4\times d`$. In each case, $`f`$ has been chosen to maximize the coupling efficiency for $`|\alpha |=0`$. At $`|\alpha |=\lambda /d`$, the coupling efficiency is 10 times smaller than it is on-axis for a non-obstructed circular aperture, and $`\rho (0.58\lambda /d)=0.5\rho (0)`$.

Therefore, the field of view of an interferometer using one single-mode fiber per aperture to transport the beam is limited to $`\lambda /d`$. This very limited field of view (FOV) is a serious constraint for the observation of many objects, such as active galactic nuclei (AGNs) and young stellar objects (YSOs). Aperture synthesis using either the Earth rotation for ground-based interferometers or a rotation of the array for space interferometers could allow reconstruction of images with a much larger FOV than permitted by the use of single-mode fibers.

A point of concern is the use of single-mode fibers on an interferometer whose telescopes have different diameters and/or pupil shapes. The coupling efficiency, as projected on the sky, is then different for each telescope, and this effect has to be taken into account for image reconstruction of objects larger than $`\lambda /d_{max}`$ ($`d_{max}`$ is the diameter of the largest aperture in the array). Residual wavefront aberrations also produce a similar discrepancy between the fiber coupling efficiencies. These effects are studied in the next 2 sections.

## 3 Observation of extended sources with a fiber interferometer

### 3.1 Instrumental response to a point source

A general description of fiber interferometers is given in Coude du Foresto et al. (1997), and we only recall here some of the basic concepts for a 2 telescope interferometer (Fig. 2).

A fiber interferometer cophased (zero optical pathlength difference (OPD) from the source to the detector) for a direction $`r_0`$ on the sky couples the electric field of the observed source on the telescopes’ focal planes into the fibers. With a 2 telescope interferometer observing an unresolved point source at a direction $`r`$, the electric fields $`E_1`$ and $`E_2`$ in the fibers, as a function of time $`t`$, are given by

$$𝐄_𝐢(t)=E_ie^{j(\alpha k_i+\omega t)}$$
(6)
for $`i=1,2`$, with $`\omega =2\pi c/\lambda `$, $`\alpha =rr_0`$ and $`k_1`$,$`k_2`$ are the positions of the 2 telescopes. The beam combiner (a fiber coupler in this case, noted X in Fig. 2) sums the two electric fields, with a time-dependent phase term $`\psi (t)`$ between the two, and the electric field at its output is

$`𝐄(t)`$ $`=`$ $`𝐄_\mathrm{𝟏}(t)+𝐄_\mathrm{𝟐}(t)e^{j\psi (t)}`$ (7)
$`=`$ $`E_1e^{j(\alpha k_1+\omega t)}+E_2e^{j(\alpha k_2+\omega t+\psi (t))}.`$
The intensity measured at the output of the beam combiner is

$$I(t)=𝐄(t)𝐄^{}(t)$$
(8)
$$I(t)=E_1^2+E_2^2+2E_1E_2cos(\alpha k+\psi (t)),$$
(9)
with the baseline vector $`k=k_1k_2`$. If a 50/50 fiber splitter is used, there are in fact two interferometric outputs (Fig. 2) :

$$I_1(t)=\frac{E_1^2+E_2^2}{2}+E_1E_2cos(\alpha k+\psi (t))$$
(10)
and

$$I_2(t)=\frac{E_1^2+E_2^2}{2}E_1E_2cos(\alpha k+\psi (t)).$$
(11)
The OPD between the two beams is usually modulated ($`\psi (t)`$ is not constant), and the temporal modulation of $`I`$ obtained by this changing OPD leads to the measurement of $`E_1^2+E_2^2`$ (mean level), $`E_1E_2`$ (fringe amplitude) and $`\varphi =\alpha k`$ (fringe phase).For example, most interferometers introduce a time-dependent OPD variation to scan the fringe packet and extract the fringe parameters in each scan.

The phase of the interferometric signal, $`\varphi =\alpha k`$, can only be referenced through the simultaneous observation of the object and a reference star or by using phase closure techniques, which requires at least 3 apertures.

### 3.2 Instrumental response to an extended source

The expression of the interferometric signal for an extended source, of light distribution $`(\alpha )`$, is obtained by integration of equation (9) over $`\alpha `$ :

$$I(t)=_\alpha E_1^2(\alpha )+E_2^2(\alpha )+2E_1(\alpha )E_2(\alpha )cos(\alpha k+\psi (t))d\alpha .$$
(12)
From equation (5),

$$E_i(\alpha )=\sqrt{(\alpha )}\frac{d_i\sqrt{\pi }}{2}\sqrt{\rho _i(\alpha )}$$
(13)
for $`i=1,2`$, where $`\rho _1`$ and $`\rho _2`$ are the coupling efficiencies functions for telescopes 1 and 2 respectively (given by equation (4) for non-obstructed circular pupils without wavefront aberrations). Therefore,

$`I(t)={\displaystyle _\alpha }(\alpha )[{\displaystyle \frac{d_1^2\pi }{4}}\rho _1(\alpha )+{\displaystyle \frac{d_2^2\pi }{4}}\rho _2(\alpha )`$ (14)
$`+2{\displaystyle \frac{d_1d_2\pi }{4}}\sqrt{\rho _1(\alpha )\rho _2(\alpha )}cos(\alpha k+\psi (t))d\alpha ]`$
### 3.3 Measurements of the object’s Fourier transform without wavefront aberrations

Several cases are briefly discussed in the absence of wavefront aberrations :

1. $`d_1=d_2=d`$, $`\rho (\alpha )=\rho _0`$
The coupling of the fiber is considered constant across the field. This corresponds to the small field approximation of fiber interferometers, a valid approximation when the object is unresolved by the individual apertures. Equation (14) becomes

$$I(t)=\rho _0\frac{d^2\pi }{2}\left[(\alpha )𝑑\alpha +e\left[e^{j\psi (t)}\stackrel{~}{}(k)\right]\right]$$
(15)
where $`\stackrel{~}{}`$ is the Fourier transform of $``$. The interferometer measures one component of the Fourier transform of the source’s light distribution per baseline.
2. $`d_1=d_2=d`$
Because the two apertures have the same diameter, the coupling efficiency is also the same for the two apertures for all values of $`\alpha `$. Equation (14) becomes

$$I(t)=\frac{d^2\pi }{2}\left[_a(\alpha )𝑑\alpha +e\left[e^{j\psi (t)}\stackrel{~}{}_a(k)\right]\right]$$
(16)
where

$$_a(\alpha )=\rho (\alpha )\times (\alpha ).$$
(17)
The interferometer measures here one Fourier component (per baseline) of an apodized image of the source. The apodization profile is the transmission profile of the fiber on the sky.
3. $`d_1d_2`$
The coupling efficiencies of the two fibers on the sky are different. Equation (14) becomes

$`I(t)={\displaystyle \frac{d_1^2\pi }{4}}\left[{\displaystyle _{a1}(\alpha )𝑑\alpha }\right]`$ (18)
$`+{\displaystyle \frac{d_2^2\pi }{4}}\left[{\displaystyle _{a2}(\alpha )𝑑\alpha }\right]+{\displaystyle \frac{d_1d_2\pi }{2}}e\left[e^{j\psi (t)}\stackrel{~}{}_{a12}(k)\right]`$
where

$$_{ai}(\alpha )=\rho _i(\alpha )\times (\alpha ),$$
(19)
for $`i=1,2`$, and

$$_{a12}(\alpha )=\sqrt{\rho _1(\alpha )\rho _2(\alpha )}\times (\alpha ).$$
(20)
The interferometer measures a Fourier component of the product of the source and an apodization mask. The apodization mask is the square root of the product of the two fiber apodization profiles on the sky.

### 3.4 Effect of time-variable wavefront perturbations

Time-variable wavefront aberrations change the fiber coupling efficiencies of each fiber for an on-axis point source. More seriously, the on-sky maps of the fiber coupling efficiency are changing in structure. In the following discussion, we note $`\rho _1(\alpha ,t)`$ and $`\rho _2(\alpha ,t)`$ the coupling efficiencies of the two fibers for the sky position $`\alpha `$ at time $`t`$.

#### 3.4.1 Observation of an object of small size

Here, we consider an object which is not resolved by the individual apertures : the size of the object is much smaller than the diffraction limit of each telescope. The coupling efficiency is then constant across the object for each aperture. Because $`\rho _1`$ and $`\rho _2`$ are not functions of $`\alpha `$, equation (14) becomes :

$`I(t)={\displaystyle \frac{d_1^2\pi }{4}}\rho _1(t){\displaystyle (\alpha )𝑑\alpha }+{\displaystyle \frac{d_2^2\pi }{4}}\rho _2(t){\displaystyle (\alpha )𝑑\alpha }`$ (21)
$`+{\displaystyle \frac{d_1d_2\pi }{2}}\sqrt{\rho _1(t)\rho _2(t)}e\left[e^{j\psi (t)}\stackrel{~}{}(k)\right].`$
The complex visibility of the observed object, $`V(k)`$, can be expressed as

$$V(k)=\frac{\stackrel{~}{}(k)}{\stackrel{~}{}(0)}.$$
(22)
The photometric outputs (Fig. 2), $`P_1(t)`$ and $`P_2(t)`$, are measuring respectively $`\rho _1(t)`$ and $`\rho _2(t)`$:

$$P_i(t)=\beta _i\frac{d_i^2\pi }{4}\rho _i(t)(\alpha )𝑑\alpha ,$$
(23)
for $`i=1,2`$, where $`\beta _1`$ and $`\beta _2`$ are the fraction of fiber flux being sent into the photometric outputs by the fiber splitter $`Y_1`$ and $`Y_2`$ respectively. $`Q_1`$ and $`Q_2`$ represent the flux in each fiber to the combiner (fiber splitter X).

$$Q_i(t)=\frac{1\beta _i}{\beta _i}\times P_i(t)$$
(24)
for $`i=1,2`$.

$$I(t)=Q_1(t)+Q_2(t)+\sqrt{Q_1(t)\times Q_2(t)}\times e\left[e^{j\psi (t)}V(k)\right].$$
(25)
From equations (24) and (25), it appears that $`|V(k)|`$ can be computed from the photometric and interferometric output(s) $`P_1(t)`$, $`P_2(t)`$ and $`I(t)`$. Precise measurement of $`\psi (t)`$ is required to compute the phase of $`V(k)`$. This self-referencing technique is used on the fiber-fed FLUOR beam combination unit of the IOTA interferometer.

#### 3.4.2 Observation of an extended object

If the size of the object is a significant fraction of the diffraction limit of the telescopes, $`\rho `$ is not uniform across the object: for a given telescope, the fiber coupling efficiency is a function of both the angular vector $`\alpha `$ and of time. Figure 3 shows a typical instantaneous coupling efficiency $`\rho `$ as a function of position on the sky $`\alpha `$ for two values of atmospheric turbulence. The photometric outputs are:

$$P_i(t)=\beta _i\frac{d_i^2\pi }{4}\rho _i(\alpha ,t)(\alpha )𝑑\alpha ,$$
(26)
for $`i=1,2`$. The interferometric output is:

$$I(t)=Q_1(t)+Q_2(t)+\sqrt{Q_1(t)Q_2(t)}e\left[e^{j\psi (t)}\frac{\stackrel{~}{}_{a12}(k)}{\stackrel{~}{}_{a12}(0)}\right]$$
(27)
where

$$_{a12}(\alpha )=I(\alpha )\sqrt{\rho _1(\alpha ,t)\rho _2(\alpha ,t)}.$$
(28)
In this case, the knowledge of $`I(t)`$, $`P_1(t)`$ and $`P_2(t)`$ is not sufficient to compute the visibility $`|V(k)|`$ of the object: $`\rho _1(\alpha ,t)`$ and $`\rho _2(\alpha ,t)`$ also need to be known to accurately interpret the observed fringe visibility. To illustrate this effect, the simple example of the observation of a double star observed with two identical telescopes is considered. The wavefront of each telescope is affected by a random tip-tilt, due to atmospheric turbulence or telescope pointing errors. $`\alpha _1`$ and $`\alpha _2`$ are the positions of the two stars, of identical brightness and both unresolved by the interferometer. Three cases for an observation at time $`t_0`$ are considered:
(1) Telescope 1 is pointing at star 1 and telescope 2 is pointing at star 2 ($`\rho _1(\alpha _1,t_0)=1`$, $`\rho _1(\alpha _2,t_0)=0`$, $`\rho _2(\alpha _1,t_0)=0`$ and $`\rho _2(\alpha _2,t_0)=1`$).
(2) Telescope 1 and 2 are both pointing halfway between star 1 and 2 ($`\rho _1(\alpha _1,t_0)=\rho _1(\alpha _2,t_0)=\rho _2(\alpha _1,t_0)=\rho _2(\alpha _2,t_0)=0.5`$).
(3) Telescope 1 and 2 are both pointing at star 1 (($`\rho _1(\alpha _1,t_0)=\rho _2(\alpha _1,t_0)=1`$, $`\rho _1(\alpha _2,t_0)=\rho _2(\alpha _2,t_0)=0`$).
$`P_1(t_0)`$ and $`P_2(t_0)`$ are identical in each case : it is impossible to tell between the 3 cases from the values of the photometric signals. In case 1 there is no coherence between the light in the two fibers and $`_{a12}=0`$ and the measured fringe visibility is 0, while in case 2, $`_{a12}=0.5\times `$ and the measured fringe visibility corresponds to the separation of the 2 stars. Finally, in case 3, the observed visibility is 1, because the interferometer behaves just as if there was only one star. During real observations, the atmospheric seeing and pointing errors of the telescopes would induce variations into the measured visibilities which cannot be corrected by measurements of $`P_1`$ and $`P_2`$.

## 4 Effect of atmospheric turbulence on object visibility measurements.

In this section, the amplitude of the effect presented above is estimated for the observation of various astrophysical sources. In these numerical simulations, only the effect of atmospheric turbulence has been simulated, for various values of $`r_0`$, the Fried parameter. This is adequate for ground-based interferometers without Adaptive Optics (such as CHARA, PTI, IOTA and NPOI). The effect of atmospheric turbulence on ground-based interferometers with adaptive optics (such as Keck, OHANA and VLTI) can still be obtained from this study by adopting a $`r_0`$ value which best represents the corrected wavefront: the value of $`r_0`$ that would yield the full width at half maximum (FWHM) of the corrected long-exposure PSF if there was no AO system. For each telescope+AO system, this derived value of $`r_0`$ will decrease as a function of the brightness of the source used for wavefront sensing.

To reproduce the effect of atmospheric turbulence, a simple 2-telescope fiber-fed interferometer has been simulated. Two independent time series of 2000 wavefronts ($`d/r_0=1`$) were first generated, one for each telescope. These wavefronts are then scaled to the appropriate value of $`d/r_0`$ for the simulation. For each step, the simulator, written in C, first computes the on-sky fiber coupling efficiency map for each telescope. These two maps are then multiplied by the intensity map of the source to yield the values of the photometric outputs (equation 26). The square root of the product of these two maps is multiplied by the intensity map of the source (equation 28), and the result is then Fourier transformed to obtain the uncorrected fringe visibility, which is then corrected using the values of the photometric outputs (equation 27). Although the code can work with different telescope pupils (shape, diameter, central obstruction etc…), the results presented in this work were obtained with two circular pupils (no central obstruction) of identical diameter.

### 4.1 One example : the observation of a double star

Mass estimates of stars rely upon precise measurements of orbits of double stars. Thanks to interferometers, it is possible to make such measurements on otherwise unresolved binaries (spectroscopic binaries).

In this simulation, one of the stars is on the optical axis of each telescope. As can be seen on Fig. 4, when the second star is close to the edge of the field of view of the fiber ($`\alpha /(\lambda /d)`$ is close to 1), the measured visibility is affected by the vignetting introduced by the fiber. One can correct for this effect because the electric field $`E_{01}`$ of the fundamental mode of the fiber is known. Atmospheric turbulence will introduce a variation (as a function of time) of the measured visibilities. One might think that it is possible to correct for this effect by averaging a large number of measurements. Figure 5 shows that such an average (here, an average of 2000 visibility measurements, each one of those has been previously been corrected according to the values of the photometric outputs $`P_1`$ and $`P_2`$) is significantly different from the measurement that would be obtained without atmospheric turbulence. The difference is especially large when the visibility of the object is low (the interferometer is resolving the double star): the error can then be larger than the measured value of the visibility. Figure 6 shows how these errors affect the measure of the system separation.

For $`d/r_0`$ values larger than 1, the errors are almost independent of $`d/r_0`$. This is an expected result since the error is then a function of the ratio between the size of the object (separation between the 2 stars in this example) and the size of a speckle: this ratio defines the gradient of fiber coupling efficiency across the object, which is the cause of this effect. Figure 7 shows how those errors vary as a function of $`d/r_0`$ for $`\alpha /(\lambda /d)=0.22`$ and $`B/d=5.92`$. Errors for other values of $`\alpha /(\lambda /d)`$ and $`B/d`$ vary in a similar way. It is found that the standard deviation between measurements is very close to the systematic error introduced by the atmospheric turbulence. This is true for all points except the ones close to $`\alpha /(\lambda /d)=0.7`$, for which, although the standard deviation is relatively high, the average visibility measurement is not significantly affected by atmospheric turbulence. This “equilibrium” value of $`\alpha /(\lambda /d)`$ is a function of $`d/r_0`$ and is therefore not constant during an observation: this property cannot be used to hope to obtain exact visibility measurements. It is found that, in order to significantly decrease the visibility measurement error, $`d/r_0`$ needs to be smaller than 1. Existing Adaptive Optics systems deliver Strehl ratios of typically 30% to 70% on bright sources in the near-infrared. This degree of correction, as suggested in Fig. 7, can only reduce the visibility measurement error by a factor of 2.5 at best. The adaptive optics correction is lower for fainter sources and in many cases, Adaptive Optics does not significantly reduce the visibility measurement error. Therefore, the results presented in this study for $`d/r_0=1`$ are also valid for most ground-based interferometers, with or without adaptive optics. High-performance ($`Strehl>0.9`$) could however greatly improve the quality of the measurements.

### 4.2 Other types of sources : Stellar diameter measurement, AGNs and YSOs

Existing interferometers can measure stellar visibilities with an accuracy of about one percent. If the star is partially resolved by a telescope of the interferometer, it is however often impossible to reach this accuracy. We consider the observation of a star of diameter $`d_{star}`$ with a 2-telescope interferometer. The error on the measurement of the fringe visibility due to the apparent size of the star and atmospheric turbulence have been computed with the same code as for the double star (averaging of 2000 independant photometry-corrected visibility measurements). Figure 8 shows how these errors affect the precision of the stellar diameter measurement (with $`d/r_0=1`$). The errors are generally smaller than in the case of the double star because the flux is more concentrated in the central region. Moreover, a very simple model (uniform disk) is adopted for the star and only one simple quantity (the diameter) is estimated. The errors tend to be large (more than a few percent) if the visibility is measured close to the first null. For $`B/d>2`$, the accuracy on the stellar diameter is generally better than $`1\%`$ if the measure is done at baselines shorter than the resolving baseline (for which the visibility is null).

The effect of atmospheric turbulence on the fringe visibility measurements of AGNs and YSOs is hard to quantify precisely due to the lack of a simple and accurate model of light intensity distribution for those sources. However, the two examples studied above (double stars and stellar disks) provide us with some general estimate of the measurement errors for various sources and some conclusions can be drawn:

* Atmospheric turbulence seriously alters visibility measurements of sources that are not very concentrated: stellar diameters measurements are less corrupted than double star separation measurements.
* The error on the measured physical quantity is larger when the visibility of the object is small. Measurements near or beyond the first zero of the visibility curve should be avoided.
* This effect still allows measurements of simple physical quantities (stellar diameters) but is likely to seriously alter more evolved measurements (small variations of stellar diameters, limb darkening, aperture synthesis imaging).

## 5 Discussion

Imaging of extended sources with fiber-fed interferometers is made very delicate because of two effects.

* The field of view of the single mode fiber is small ($`\lambda /d`$).
* When the source size is a significant fraction of the field of view of the single mode fiber, the visibility measurements are corrupted by atmospheric turbulence.

##### Field of view of the fiber:

To increase the field of view of the single mode fiber, the telescope diameter could be reduced, at the cost of lowering the flux in the fibers. This solution is of limited use because it seriously affects the limiting magnitude of the interferometer. Another obvious solution is to increase the operating wavelength.

##### Visibility measurement errors due to atmospheric turbulence:

This effect has been studied in §4. One might think that because these errors can be estimated, the visibility measurements could be corrected for, using this work. Unfortunately, during actual observations, $`d/r_0`$ does not stay constant and additional effects (dome seeing, mirror seeing, shape of the primary mirror etc.) prevent such correction: the wavefront perturbation key parameters are not constant in time. This effect is also a function of the morphology of the source, which is generally poorly known. However, several solutions to this problem can be implemented without much complexity :

* Better statistical analysis of the measurements Measurements could be “weighted” according to the values of the photometric outputs. For example, visibility measurements are often very corrupted when at least one of the photometric outputs is low.
* Reducing the size of the apertures Figures 6 and 8 show that reducing $`d`$ (the point of measure in these figures is then moved to the upper left, along a line of constant visibility, parallel to the dashed line) decreases the visibility measurement error. This technique is only efficient for bright (the apertures can be stopped down) sources.
* Monitoring the PSF Real-time imaging of the two PSFs on the fiber heads can yield to an estimate of $`\rho _1(\alpha ,t)`$ and $`\rho _2(\alpha ,t)`$ (equations 26 to 28).

The statistical analysis of the measurements can help with the standard deviation between measurements, but most of the systematic visibility measurement errors cannot be corrected for. Reducing the size of the aperture affects the limiting magnitude of the interferometer and degrades the precision of the measurements because of photon noise. Monitoring of the PSF seems to be the most promising technique and can be done by imaging the image reflected by the fiber head with a high-speed camera. Alternatively, efficient photon-counting detectors, such as APDs, could be used, since only a small number of “pixels” are required. This last technique is efficient for sources smaller than $`\lambda /d`$, because large sources will tend to “blur” the PSF (the information that leads to estimates of $`\rho _1`$ and $`\rho _2`$ is lost).

One elegant solution to both effects is to use several fibers per telescope. By segmenting an aperture into subapertures and using one fiber per subaperture, the field of view of the fibers is increased and the visibility measurement errors due to atmospheric turbulence are reduced. This approach offers the advantage of reducing the size of the apertures (see above) without significant loss of light.

## 6 The multifiber apertures

### 6.1 The principle

The use of one fiber per aperture limits the FOV to the size of a coherent domain in the focal plane, $`\lambda /d`$. Increasing the FOV can be done by increasing the size of this coherent domain (reducing $`d`$) or by coupling $`N`$ coherent domains in $`N`$ single-mode fibers. There are therefore two solutions to increase the FOV.

* Reducing the size of the pupil to be coupled in a fiber. To avoid losing most of the signal, the pupil is segmented into N subpupils and the light of each subpupil is coupled to a fiber.
* Paving the focal plane with an array of fibers. Each coherent domain is coupled to a fiber. Alternatively, this technique could be implemented sequentially, by mapping the field of view with numerous pointings of the array.

Paving the focal plane with fibers would increase the FOV but poses problems for beam recombination. We have therefore chosen to use the first solution. A possible implementation of this scheme is shown on Fig. 9: a lenslet array in a pupil plane is used to couple each subaperture to its corresponding fiber. This is very similar to the optical scheme used to couple subapertures to multi-mode fibers in curvature adaptive optics systems.

All subapertures of the interferometer have the same physical size. As a result, the number of fibers on a given aperture is proportional to its area.

### 6.2 Field of view

$`N`$ is the number of apertures in the array, $`d_i,i=1\mathrm{}N`$, the diameter of each of those apertures, $`d_{sa}`$ the diameter of each subaperture and $`M_i,i=1\mathrm{}N`$, the number of fibers for each aperture. Since $`d_{sa}`$ should be the same for all subapertures,

$$M_i=\left(\frac{d_i}{d_{sa}}\right)^2.$$
(29)
The field of view of the interferometer is

$$FOV=\frac{\lambda }{d_{sa}}=\sqrt{M_i}\times \frac{\lambda }{d_i}$$
(30)
for $`i=1\mathrm{}N`$. At each telescope, the FOV is multiplied by the square root of the number of fibers used. Paving the focal plane with optical fibers would have resulted in the same increase of FOV.

### 6.3 Coupling efficiency

Figure 10 shows four possible subaperture geometries. We compute the coupling efficiency for each of those geometries, neglecting the edge effects ($`M_i1`$).

In geometry (a) and (b), the coupling efficiency is the product of the filling factor of the disks across the pupil and the coupling efficiency of a circular pupil into a fiber:

$$\rho _a=\frac{\pi }{4}\times 0.82=0.64$$
(31)
$$\rho _b=\frac{\pi }{4sin(\frac{\pi }{3})}\times 0.82=0.74.$$
(32)
For geometries (c), and (d), we have computed the optimal coupling efficiency for square and hexagonal apertures using equation (1).

$$\rho _c=0.793$$
(33)
$$\rho _d=0.815$$
(34)
A tight hexagonal paving of the pupil plane allows a coupling efficiency $`\rho _d`$ almost (0.5 % difference) as good as a direct coupling of the entire circular pupil into a single fiber. For apertures with a central obstruction of diameter $`a\times d`$, equation (2) becomes

$`\rho (\alpha )={\displaystyle \frac{8}{\omega ^2}}e^{2\left(\frac{\alpha f}{\omega }\right)^2}\times `$ (35)
$`\left[{\displaystyle e^{\left(\frac{r}{\omega }\right)^2}I_0(2r\alpha f/\omega ^2)\left[J_1(\pi dr/\lambda f)aJ_1(a\pi dr/\lambda f)\right]𝑑r}\right]^2.`$
Using this equation, we can compute, for each value of $`a`$, the combination of $`\omega `$, $`f`$ and $`d`$ which yields the best coupling efficiency. As can be seen on Fig. 11, the coupling efficiency decreases rapidly with the size of the central obscuration. For a central obstruction of 0.4, the coupling efficiency is 56%. However, if the pupil is subdivided into hexagonals subapertures, each coupled to a fiber, the overall coupling efficiency is $`\rho _4=0.815`$.
Under diffraction-limited conditions, the division of the pupil in several hexagonal subpupils does not result in a loss of coupling efficiency. In the case of a circular pupil with a central obstruction, there is a gain of coupling efficiency.

### 6.4 Beam combining

The high number of optical fibers in this concept is a strong constraint for the design of a beam combiner. Two possible schemes have been identified.

* Reconstruction of the entrance pupils. The most intuitive approach is to optically reconstruct the incoming wavefront (in fact, some approximation of it) for each pupil of the interferometer, as seen on Fig. 12. Wide-field beam recombination techniques can then be used as if each entire telescope beam had been transported. In this scheme, the fibers (several per telescope) are only used to transport the wavefront from the telescopes to the beam combiner. One attractive approach for beam combination is a Fizeau-type beam combiner (shown on Fig. 12), where the beams are combined in a common focal plane. The pupil of the beam combiner can be identical to the entrance pupil, it can be a densified (Labeyrie (1996)) version of it, or it can be a fixed pupil. In the first 2 cases, the beam combiner pupil needs to be rearranged to follow the geometry of the entrance pupil of the array (as seen by the source).
* Use of integrated optics components. A more promising technique is the use of a large number of fiber couplers to determine the complex amplitude of each subpupil in the array, relative to a common (for the whole array) phase. Phase closure techniques could be used between subapertures of different telescopes if more than 2 telescopes are used. Integrated optics components (Malbet et al (1999); Berger et al. (1999, 2001)) can include a high number of fiber couplers in a small volume and seem to be a very promising technology for this beam combination scheme, illustrated in Fig. 13. In this case, the beam combiner considers the fibers as the coherent beams from a high number of sub-telescopes (as many as there are fibers). The beams of sub-telescopes belonging to the same telescope are already cophased thanks to adaptive optics. It would be very interesting to determine what is the optimal set of beam couplings for wide-field imaging with a limited number of photons.

### 6.5 Photometric calibration

In this concept, the coupling field of view of an individual fiber is increased by dividing the entrance pupil into subapertures, each of those being coupled to a fiber. When this coupling field of view is significantly larger than the observed source, the effect presented in §3 and §4 is greatly reduced. If this coupling field of view is significantly larger than the atmospheric seeing, the flux in each fiber is constant in time and there is no need for photometric calibration of the interferometer.

In practice, this second requirement is hard to meet because of the very large number of fibers required. Therefore, photometric calibration of the interferometer would still greatly improve the fringe visibility measurement accuracy, as it currently does for monofiber interferometers (FLUOR instrument on the IOTA interferometer). This photometric calibration of the interferometer would be done by spitting the light of each fiber before the beam combination, yielding one photometric output per subaperture. Such a calibration unit could be integrated directly into the integrated optics processor presented in §6.4. It is very interesting to note that in this scheme, the different subapertures can be considered as different telescopes of an interferometric array, regardless of their belonging to the same main aperture.

## 7 Use on a space interferometer

### 7.1 (u,v) plane coverage and field of view.

Guyon & Roddier (2001) have demonstrated that with only six 2-meter apertures, a full (u,v) plane coverage up to a 60m baseline is accessible with the rotation of an array. Wide field of view imaging is possible with such arrays if the phase information across the pupil is preserved: collapsing each pupil on a single mode fiber would reduce the field of view. To exploit this full (u,v) plane coverage and reconstruct wide field images, the flux received by the telescope over a large solid angle needs to be brought to the beam combiner. We consider such a optimized array geometry and review the characteristic field of view quantities for a fiber interferometer using rotational aperture synthesis. In this section, we consider a array of $`N`$ identical telescopes: $`d_0=d_1=\mathrm{}=d_N=d`$ and $`M_0=M_1=\mathrm{}=M_N=M`$.

* Coupling field of view
Using the results of §6.2,

$$FOV_c=\frac{\lambda }{d_{sa}}=\sqrt{M}\times \frac{\lambda }{d},$$
(36)
with $`d`$, the diameter of each telescope, $`d_{sa}`$ the diameter of each subaperture and $`M`$, the number of fibers for each aperture.
* Fourier field of view
We first consider an array of N telescopes of infinitely small diameter. There are $`N(N1)+1`$ (u,v) points measured per snapshot exposure (two symmetrical points per baseline plus the origin). In the case of a rotating array, the (u,v)-plane coverage of the final reconstructed image is a series of $`N(N1)/2`$ circles (excluding the origin). This correctly describes the (u,v) plane coverage if each aperture’s flux is coupled into a fiber (the information about the size of the pupil is not preserved). The maximum distance between two consecutive baseline lengths measured is

$$d_c=e(N)\frac{B}{N(N1)/2}$$
(37)
where $`e(N)`$ is the efficiency of the (u,v) plane coverage of the array (Guyon & Roddier (2001)) and $`B`$ is the longest baseline of the array. This efficiency decreases from 1 to 0.869 as $`N`$ increases from 4 to 10. In this simple estimate, we will consider $`e(N)=1`$, which is a valid approximation for $`N<10`$. In a (u,v)-plane coverage obtimized array, the diameter or the individual apertures is sufficient to fill the “gap” between consecutive baselines: $`d=d_c`$ (Guyon & Roddier (2001)).

Reconstructing an image requires one to interpolate the (u,v) plane between the measured spatial frequencies: this is equivalent to a convolution in the (u,v) plane by a kernel of characteristic size equal to the distance between the known (u,v) points. In this case, this distance is $`d_c`$. Therefore, with only one fiber per aperture, this imposes a limitation $`FOV_f`$ on the field of view of the reconstructed image:

$$FOV_f=\frac{\lambda }{d_c}=\frac{\lambda N(N1)}{2B}.$$
(38)
If no interpolation is done in the (u,v) plane, there is no “windowing” in the image plane, but the PSF has strong secondary peaks separated by $`FOV_f`$: the clean field of view would therefore be $`FOV_f`$. In the case of a subdivision of each pupil into $`M`$ hexagonal subpupils ($`M`$ fibers per aperture), equation (37) becomes

$$d_c<\frac{e(N)\frac{B}{N(N1)/2}}{\sqrt{M}}$$
(39)
and equation (38) becomes

$$FOV_f>\sqrt{M}\frac{\lambda N(N1)}{2B}.$$
(40)

From equations (36) and (40), and noting that, for an array with a full (u,v) plane coverage ($`d=d_c`$), with $`e(N)1`$, equation (37) yields $`N^2d2B`$,

$$FOV_f>FOV_c.$$
(41)
Therefore, on the reconstructed image, the field of view will be limited by the coupling of the light into the fibers, and

$$FOV=\sqrt{M}\times \frac{\lambda }{d}.$$
(42)
Equation (42) states that the number of (u,v) points measured is sufficient to yield a clean PSF without strong secondary peaks inside the field of view permitted by the fiber coupling.

### 7.2 Spatial filtering by the fibers

Single-mode fibers have been proposed for use on space interferometers because of their spatial filtering capabilities. The wavefront errors at the telescope’s pupil result in a loss of coupling efficiency into the fiber but there is no loss of coherence inside the fiber: it only allows one coherent mode to be coupled. Because image reconstruction and nulling are more tolerant to amplitude variations than phase variation, the constraints on the wavefront quality are reduced by the use of such fibers.

In this concept, the wavefront over each subaperture is cleaned by this spatial filtering and the coupling efficiency decreases with the amount of wavefront distortion over this subaperture. The spatial filtering is done at a higher spatial frequency than if only one fiber was used per aperture. The signals in the fibers are perfectly cophased only if the wavefront is flat at spatial frequencies lower than $`1/d_{sa}`$: it is therefore essential to accurately control the quality of the wavefront at low spatial frequencies, either before injection into the fibers (adaptive optics) or after (active delay lines to cophase the fibers). Spatial filtering transforms the wavefront distortions of higher spatial frequencies into coupling efficiency losses, which can also affect the the visibility measurements. This effect can be corrected for by photometric calibration of the interferometer, as discussed in §6.5.

## 8 Use on a ground-based interferometer : OHANA

### 8.1 Presentation of the OHANA project

OHANA (Optical Hawaiian Array for Nanoradian Astronomy) makes use of the existing optical and near-infrared telescopes on the summit of Mauna Kea to build a fiber interferometer (Mariotti et al. (1996, 1998)). Optical fibers are preferred to classical beam propagation for the simplicity of the beam collection at the telescope foci and the easy transport of coherent light over hundreds of meters. The spatial filtering is also very desirable since the Strehl ratios delivered by the adaptive optics systems on the OHANA telescopes, in the near-infrared, are typically about 50% or less.

### 8.2 Fiber coupling Field of View

The OHANA array, because it makes use of pre-existing telescopes, is made of apertures of various sizes, from the 10m-class telescopes (Keck I, Keck II, Subaru and Gemini) to the 4m-class telescopes (CFHT, IRTF and UKIRT). The field of view of the interferometer is the smallest FOV seen by the fibers. If only one fiber per telescope is implemented, the field of view is limited to 20 milli-arc-second at $`\lambda =1\mu m`$ and 40 milli-arc-second at $`\lambda =2\mu m`$.
A more serious problem for the imaging of extended objects is the matching of the field of view for the different telescopes studied in §3.3: with only one fiber per telescope, the field of view coupled into the fibers varies by a factor of 2 from telescope to telescope (factor of 4 in area). This is a serious problem when the observed source overfills the smallest field of view. Another obstacle to high accuracy object visibility measurements is the effect of residual wavefront aberrations, studied in §4.

The use of several fibers per aperture, as presented in §6, would solve these two problems and allow observations of extended objects. Table 2 lists, for each telescope of the OHANA array, the approximate number of fibers required for an observation with a specific FOV.

### 8.3 (u,v) plane coverage

The instantaneous (u,v) plane coverage of the OHANA array is very sparse due to the large baselines and relatively small number of apertures. However, aperture supersynthesis greatly improves this situation: the measured (u,v) points lie along arcs that are drawn as the object moves across the sky.

The size of the telescopes of the array (up to 10 meters) is not negligible when compared to the distance between them (up to 800m). Therefore, a significant gain in (u,v) plane coverage can be achieved by preserving the phase information across the pupil. The points measured in the (u,v) plane are then thick arcs. Figure 14 shows the (u,v) coverage accessible to OHANA when observing sources at various declinations, when their elevation is more than 20 degrees above the horizon. Although there is a significant gain in (u,v) plane coverage, the number of telescopes and their diameters are still too small to fill most of the frequency gaps between the arcs. The Fourier clean field of view, which is given by the spacing between the arcs, is about $`\lambda /20m`$ for OHANA, which is about twice as small as the field of view of the fibers of the big telescopes of the array. Therefore, in this case, the use of several fibers per aperture is mostly motivated by the need to reduce the turbulence-induced visibility measurement errors in the array and increase the coupling FOV: the size of the gaps in the (u,v) plane coverage is not significantly reduced by using several fibers per telescope.

## 9 Conclusion

Wide field imaging with fiber-fed interferometers (and more generally with interferometers which use spatial filtering) is made very difficult because of the small field of view of the fibers (or the pinhole). This work has also demonstrated that the spatial filtering of partially corrected turbulent wavefronts limits the fringe visibility measurement accuracy. This effect is especially serious for the next generation of interferometers (Keck, VLTI, OHANA) which will make use of large (10m-class) telescopes with adaptive optics. For such projects, this work has shown that photometry-calibrated fringe visibility measurement errors will commonly be of the order 5% or more for extended objects (AGNs and YSOs) if only one fiber per telescope is used. The use of several fibers per aperture is an attractive solution to both improve the photometry-corrected visibility measurement accuracy and extend the field of view of these interferometers, and it can also allow fiber-based wide field interferometric imagers to operate from space.

###### Acknowledgements.

The author is very grateful to Olivier Lai, Guy Perrin and Katherine Roth for constructive discussions about this work and suggestions to improve the manuscript.