text:
# Entropic-acoustic instability of shocked Bondi accretion I. What does perturbed Bondi accretion sound like ?

## 1 Introduction

The time variability of the emission from X-ray binaries may be related to hydrodynamic instabilities in the process of accretion onto a compact star. The simplest accretion flow is the stationary, spherically symmetric accretion of an adiabatic gas onto a compact object, first studied by Bondi (1952). The stability of Bondi accretion was established by Garlick (1979) and Petterson et al. (1980) in the newtonian limit, and by Moncrief (1980) in the framework of general relativity. Kovalenko & Eremin (1998) suggested that non radial perturbations might reach non linear amplitudes in the supersonic part of the flow if the size of the accretor is small enough. This effect was invisible in the 3-D numerical simulations of the Bondi flow by Ruffert (1994a), which confirmed its stability. By contrast, a strong instability appeared in the numerical simulations of Bondi-Hoyle-Lyttleton accretion when a bow shock is present (e.g. Ruffert 1994b for $`\gamma =5/3`$). The Bondi-Hoyle-Lyttleton accretion is named after the pioneering works of Hoyle & Lyttleton (1939) and Bondi & Hoyle (1944) concerning an accretor moving at supersonic speed with respect to the gas. Foglizzo & Tagger (2000) (hereafter FT2000) suggested that this instability might be explained by what they called the entropic-acoustic instability. This generic mechanism of instability in shocked converging flows is based on the cycle of entropic and acoustic waves in the subsonic region of the flow between the shock and the sonic surface surrounding the accretor. The advection of entropy perturbations towards the accretor produces outgoing acoustic waves which propagate towards the shock. Perturbed by these acoustic waves, the shock produces new entropy perturbations, thus closing the entropic-acoustic cycle. According to FT2000, this cycle is unstable (i.e. the amplitude of the new entropy perturbation exceeds the initial one) if the sound speed at the sonic point is much larger than the sound speed at the shock. This mechanism might play a destabilizing role in various astrophysical environments such as wind accretion or disc accretion. In this series of two papers, we describe the entropic-acoustic instability in the Bondi flow, including a spherical shock as an outer boundary (as in Foglizzo & Ruffert 1997). The study of this rather academic configuration serves two purposes:

(i) it is the first detailed description of the entropic-acoustic instability in a specific flow, beyond the general concepts and approximations of FT2000. The Bondi flow is simple enough to allow for an analytical treatment of 3D perturbations, including non radial perturbations of entropy and perturbations of vorticity.

(ii) understanding the instability of the Bondi flow with a spherical shock might be the first step towards a comprehension of the instability of non radial shocked accretion flows such as supersonic Bondi-Hoyle-Lyttleton accretion.
In this first paper, we investigate the key part of the entropic-acoustic instability: the excitation of outgoing acoustic waves by the advection of entropy and vorticity perturbations. As explained in FT2000, the efficiency of this process is directly linked to the efficiency of the refraction of acoustic waves, which is also studied. The efficiencies of these processes depend a priori on the frequency $`\omega `$ of the perturbation, its wavenumbers in spherical coordinates $`(l,m)`$, and the adiabatic index $`1<\gamma 5/3`$ of the gas. The accretor is a compact object such as a black hole, absorbing all the matter falling inside the sonic radius. The shock, which plays no role in this first paper, is introduced in the second paper of this series (Foglizzo 2001) in order to close the entropic-acoustic cycle and enable a description of the eigenmodes of the global instability.
The present paper is organized as follows. Linearized equations are established in Sect. 2. The refraction of acoustic waves is studied in Sect. 3. The excitation of acoustic waves by the advection of entropy and vorticity perturbations is examined in Sect. 4.

## 2 Linear perturbations of pressure, entropy and vorticity in the Bondi flow

### 2.1 Properties of the unperturbed flow

We consider the radial accretion flow of a gas with constant sound velocity $`c_{\mathrm{}}`$ and density $`\rho _{\mathrm{}}`$ at infinity, accelerated towards a point-like accretor of mass $`M`$, which is totally absorbing. The equations describing the unperturbed Bondi flow are recalled in Appendix A. In what follows, densities and velocities are normalized by $`c_{\mathrm{}}`$ and $`\rho _{\mathrm{}}`$, and distances are normalized to the Bondi radius $`GM/c_{\mathrm{}}^2`$. After this normalization, the unperturbed flow depends on a single parameter $`\gamma `$, the adiabatic index of the gas. The sound velocity $`c_\mathrm{s}`$ at the sonic radius $`r_\mathrm{s}`$ in the Bondi flow diverges for $`\gamma =5/3`$:

$`r_\mathrm{s}`$ $`=`$ $`{\displaystyle \frac{53\gamma }{4}},`$ (1)
$`c_\mathrm{s}`$ $`=`$ $`\left({\displaystyle \frac{2}{53\gamma }}\right)^{\frac{1}{2}}.`$ (2)
The ”natural” sonic frequency $`c_\mathrm{s}/r_\mathrm{s}`$ is directly related to the Keplerian frequency $`\omega _K`$ at the sonic radius:

$$\omega _K=2^{\frac{1}{2}}\frac{c_\mathrm{s}}{r_\mathrm{s}}.$$
(3)
### 2.2 Differential equation for small perturbations

The Euler equation, the equations of conservation of entropy and mass in the Bondi flow are linearized for small perturbations in Appendix B. After a Fourier transform in time, the perturbations are projected onto a basis of spherical harmonics $`Y_l^m`$, where $`l,m`$ are the numbers associated to the latitudinal ($`\theta `$) and longitudinal ($`\phi `$) angles. Kovalenko & Eremin (1998) studied the same set of equations in order to analyze the asymptotic behaviour of their solutions near the origin $`r=0`$. By contrast, we are interested in the effect of linear perturbations on the subsonic part of the flow ($`r>r_\mathrm{s}`$). A single differential equation of second order is obtained in Appendix B.2, with the following structure:

$$\frac{^2}{r^2}\frac{\delta p}{p}+a_1\frac{}{r}\frac{\delta p}{p}+a_0\frac{\delta p}{p}=b_0\delta S_R+b_1\delta K_R.$$
(4)
The coefficients $`a_0,a_1,b_0,b_1`$ are functions of $`r`$ which depend on the flow velocity and sound speed radial profiles in the unperturbed Bondi flow. Their rather lengthy expressions are given in Appendix B, Eqs. (77 to 80). The homogeneous equation corresponding to the left hand side of Eq. (4) describes the propagation of acoustic waves in the Bondi flow and their refraction by the sound speed gradients. The important feature of Eq. (4) is the presence of source terms on the right hand side:

(i) $`\delta S`$ is the perturbation of entropy. For the sake of simplicity, the ratio of the molecular weight $`\mu `$ to the gas constant $``$ is set to $`\mu /=1`$ throughout this paper, with no loss of generality.

(ii) $`\delta K`$ involves the curl of the vorticity vector $`\delta w\times \delta v`$ projected along the flow velocity:

$$\delta Kr^2v(\times \delta w)+l(l+1)c^2\frac{\delta S}{\gamma },$$
(5)
The radial dependence of the three components of the vorticity vector $`\delta w`$ is explicitely integrated in Appendix B.1 (Eqs. 61 to 63). Nevertheless, the vorticity contributes to the excitation of acoustic waves only through the radial part of its curl, as it appears in Eq. (5). In an adiabatic flow, both $`\delta S`$ and $`\delta K`$ are conserved when advected (see Appendix B.1). A more compact mathematical formulation is obtained in Appendix B.2.2 by writing the differential equation satisfied by the perturbation $`f`$ of the Bernoulli constant.

$`\left\{vc^2{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{v}}{\displaystyle \frac{}{r}}\right)+{\displaystyle \frac{\omega ^2\omega _l^2}{1^2}}\right\}\left(\mathrm{e}^{i\omega _R^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}f\right)`$
$`={\displaystyle \frac{i\omega vc^2}{\gamma }}\delta S_R{\displaystyle \frac{}{r}}\left[{\displaystyle \frac{1^2}{^2}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}}\right]`$
$`\delta K_R{\displaystyle \frac{c^2}{r^2}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}},`$ (6)
where the frequency $`\omega _l`$ is defined as:

$$\omega _l^2l(l+1)\frac{c^2v^2}{r^2}.$$
(7)
The structure of this differential equation is the same as Eq. (4), with much simpler coefficients. The right hand side of Eq. (6) contains two source terms proportional to $`\delta S_R`$ and $`\delta K_R`$. The pressure perturbation is related to $`f`$ by the following equation derived in Appendix B.1:

$$\frac{1}{\rho }\frac{\delta p}{t}=\frac{\mathrm{D}f}{\mathrm{D}t}.$$
(8)
## 3 Acoustic waves in the Bondi flow: the refraction cut-off

### 3.1 Acoustic refraction and regularity at the sonic point

The homogeneous equation associated with Eq. (6), describing the propagation of acoustic waves is:

$`\left\{vc^2{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{v}}{\displaystyle \frac{}{r}}\right)+{\displaystyle \frac{\omega ^2\omega _l^2}{1^2}}\right\}\left(\mathrm{e}^{i\omega _R^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}f\right)`$
$`=0.`$ (9)
This equation is of second order, which reflects the fact that pressure perturbations can be decomposed into ingoing and outgoing acoustic waves. Such an identification is possible only in the region of the flow where the lengthscale of the flow gradients is longer than the wavelength of the perturbation, where a WKB treatment of the wave propagation can be used (see Appendix C). This is always the case far from the accretor in the Bondi flow, since the density and sound speed are uniform at infinity. We may define $`\delta p^+`$ and $`\delta p^{}`$ as the solutions of the homogeneous equation propagating respectively inward and outward at infinity. Following FT2000, the index $`+`$ (resp. $``$) refers to waves propagating in the same (resp. opposite) direction as the flow. The acoustic flux carried by these acoustic waves is defined by Eq. (6) of FT2000:

$$F^\pm =\dot{M}_0c^2\frac{(1\pm )^2}{}\left|\frac{\delta p^\pm }{\gamma p}\right|^2.$$
(10)
Let us normalize $`\delta p_0^\pm `$ so that they carry the same acoustic flux at infinity. Any linear combination $`\delta p`$ of the two independent solutions $`\delta p_0^+,\delta p_0^{}`$, with complex coefficients $`\alpha ,\beta `$, is also a solution of Eq. (9):

$$\delta p\alpha \delta p_0^{}+\beta \delta p_0^+.$$
(11)
There is, however, only one ratio $`\alpha /\beta `$ such that $`\delta p`$ is regular at the sonic radius. Indeed, a local analysis of the solutions of Eq. (6) for $`\gamma <5/3`$, at the sonic point $`=1`$ using Frobenius series, reveals that there is only one solution crossing the sonic point regularly, the other exhibiting a logarithmic singularity (Bender & Orszag 1978, Chap. 3.3, p.68). This regularity condition reflects the fact that an incoming acoustic wave $`\beta \delta p_0^+`$ is only partially transmitted to the accretor: an outgoing acoustic wave of amplitude $`\alpha \delta p_0^{}`$ is refracted out. Using the same notations as FT2000, the fraction $`|\overline{}_s|^2`$ of refracted acoustic flux is deduced from Eqs. (10) and (11):

$`|\overline{}_s|^2`$ $``$ $`{\displaystyle \frac{F^{}}{F^+}},`$ (12)
$`=`$ $`\left|{\displaystyle \frac{\alpha }{\beta }}\right|^2.`$ (13)
The refraction coefficient is therefore directly related to the regularity of the solution at the sonic point.

### 3.2 Numerical calculation of the homogeneous solution

The differential equation (9) is integrated using a Runge-Kutta implicit method, from the sonic point to the WKB region. The refraction coefficient is then computed by identifying the ingoing and outgoing waves in the WKB region. If $`\gamma <5/3`$, a Frobenius expansion is used in the vicinity of $`r_\mathrm{s}`$ to start the integration away from the singularity. The calculation for $`\gamma =5/3`$ and $`l=0`$ is similar to the case $`\gamma <5/3`$, except that the sonic point is at $`r=0`$. If $`\gamma =5/3`$ and $`l1`$, the homogeneous solution presents an essential singularity at the sonic point, where one branch converges to zero, and the other diverges. Nevertheless, an asymptotic expansion of the converging solution (see Eqs. 120 and 121) allows a numerical integration from the vicinity of the origin towards infinity. As expected by Kovalenko & Eremin (1998), no overreflection occurs ($`|\overline{}_s|1`$). The typical frequency dependence of $`|\overline{}_s|^2`$ is shown in Fig. 1 for $`\gamma =1.33`$.

### 3.3 Cut-off frequency for the refraction of acoustic waves for $`\gamma <5/3`$

If $`\gamma <5/3`$, the accretor is surrounded by a supersonic region from which acoustic waves cannot escape. Only those waves with a high enough frequency may penetrate this region, because lower frequency waves are refracted out before the sonic radius. This leads us to expect the refraction $`|\overline{}_s|^2`$ to decrease to zero at high frequency, and be close to unity at low frequency (see Fig. 1). The refraction cut-off $`\omega _l^{\mathrm{cut}}`$ of a wave of latitudinal order $`l`$ is defined as the intermediate frequency corresponding to a refraction of half of the incoming acoustic flux:

$$|\overline{}_s|_{(\omega =\omega _l^{\mathrm{cut}})}^2\frac{1}{2}.$$
(14)
Bessel functions are used in Appendix D to approximate the homogeneous solution at high frequency, in order to obtain the leading order of the asymptotic behaviour of $`|\overline{}_s|`$. This latter decreases exponentially above $`\omega _l^{\mathrm{cut}}`$ (see Fig. 1):

$`|\overline{}_s|(\omega \omega _l^{\mathrm{cut}})`$ $``$ $`\mathrm{e}^{\frac{\omega }{\omega _{\mathrm{}}}},`$ (15)
$`\omega _{\mathrm{}}`$ $``$ $`{\displaystyle \frac{c_\mathrm{s}|\dot{}_\mathrm{s}|}{\pi }}={\displaystyle \frac{4}{(53\gamma )\pi }}.`$ (16)
#### 3.3.1 Non radial acoustic waves $`l1`$

The turning point of non radial perturbations appears explicitly on Eq. (9). It corresponds to the radius $`r_t`$ solution of $`\omega =\omega _l(r_t)`$, where $`\omega _l`$ is defined by Eq. (7) (see Fig. 2). For $`\gamma <5/3`$, this turning point is double when the frequency equals the maximum of the function $`\omega _l(r_t)`$. The refraction coefficient is then equal to $`1/2`$, according to Bender & Orszag (1978, Chap. 10.6 p. 524):

$`\omega _{l1}^{\mathrm{cut}}`$ $`=`$ $`\mathrm{Max}_{\left\{rr_\mathrm{s}\right\}}\left\{\omega _l\right\},`$ (17)
$`=`$ $`\left[{\displaystyle \frac{l(l+1)}{2}}\right]^{\frac{1}{2}}\omega _1^{\mathrm{cut}}.`$ (18)
This result was checked by measuring numerically the refraction of acoustic waves. Fig. 3 shows the dependence of $`\omega _{l1}^{\mathrm{cut}}`$ on the adiabatic index $`\gamma `$. $`\omega _{l1}^{\mathrm{cut}}`$ should scale as $`l^{1/2}(l+1)^{1/2}/(53\gamma )^{5/4}`$ according to the asymptotic behaviour of $`(1^2)^{1/2}c/rc_\mathrm{s}(\dot{}_\mathrm{s}/r_\mathrm{s})^{1/2}`$ when $`\gamma 5/3`$, deduced from Eqs. (54 to 56). A refined comparison with the results of numerical calculations of $`\omega _{l1}^{\mathrm{cut}}`$ leads to

$$\omega _{l1}^{\mathrm{cut}}l^{\frac{1}{2}}(l+1)^{\frac{1}{2}}\left(\frac{2}{53\gamma }\right)^{\frac{5}{4}},$$
(19)
over the range $`1.4<c_\mathrm{s}/c_{\mathrm{}}<100`$ considered numerically. The radius $`r_1`$ such that $`\omega _l(r_1)=\omega _l^{\mathrm{cut}}`$ is the minimum value of the turning radius $`r_t`$ (see Fig. 2). $`r_1(\gamma )`$ is independent of $`l,m`$ and lies just outside the sonic surface, and is well fitted by $`r_1\sqrt{2}r_\mathrm{s}`$ for $`\gamma `$ close to $`5/3`$.

#### 3.3.2 Radial acoustic waves $`l=0`$

$`\omega _0^{\mathrm{cut}}`$ is computed numerically in Fig. 3, and corresponds to a lower frequency than $`\omega _1^{\mathrm{cut}}`$. The value of $`\omega _0^{\mathrm{cut}}`$ can be estimated analytically for $`\gamma `$ close to $`5/3`$ from the criterion of validity of the WKB approximation of acoustic waves, derived in Appendix C:

$$\omega (1^2)c\left|\frac{\mathrm{log}v^2c^2}{r}\right|.$$
(20)
Since the acoustic flux is conserved in the region of validity of the WKB approximation, any refraction of acoustic waves must occur outside this region. We deduce that $`\omega _0^{\mathrm{cut}}`$ scales as the maximum of the function on the right hand side of Eq. (20). When $`\gamma `$ is close to $`5/3`$, this maximum scales as $`(1^2)c/r\dot{}_\mathrm{s}c_\mathrm{s}(53\gamma )^1`$. This is confirmed by numerical calculations of $`\omega _0^{\mathrm{cut}}`$ which indicate that in the range of sound velocities $`1.4<c_\mathrm{s}/c_{\mathrm{}}<100`$ studied numerically, the ratio $`\omega _0^{\mathrm{cut}}/\omega _{\mathrm{}}`$ is constant within the numerical errors:

$$\omega _0^{\mathrm{cut}}\frac{0.85}{53\gamma }.$$
(21)
The scaling of $`\omega _0^{\mathrm{cut}}`$ should be compared to the two other particular frequencies of the flow near the sonic point, $`\omega _{l1}^{\mathrm{cut}}`$ in Eq. (19) and the local Keplerian frequency $`\omega _K`$:

$$\omega _K\frac{8}{(53\gamma )^{\frac{3}{2}}},$$
(22)
and thus $`\omega _0^{\mathrm{cut}}\omega _1^{\mathrm{cut}}\omega _K`$ for $`\gamma `$ close to $`5/3`$. The values of $`\omega _0^{\mathrm{cut}},\omega _1^{\mathrm{cut}}`$ computed numerically are compared to $`c_\mathrm{s}/r_\mathrm{s}`$ and to the Keplerian frequency in Fig. 3.

### 3.4 Refraction of acoustic waves for $`\gamma =5/3`$

In a Bondi flow with $`\gamma =5/3`$, acoustic waves are perfectly refracted at any frequency:

$$|\overline{}_s|=1\mathrm{for}\gamma =\frac{5}{3}.$$
(23)
This can be deduced from the limit of flows with $`\gamma `$ close to $`5/3`$, because the cut-off frequency $`\omega _l^{\mathrm{cut}}`$ tends to infinity when $`\gamma 5/3`$. This is not surprising, since an ingoing acoustic wave of any frequency always ultimately meets ”a wall”, i.e. a subsonic region where the scale of the flow gradients are much shorter than the wavelength.
The WKB analysis at high frequency (Appendix C) indicates that ingoing and outgoing acoustic waves are coupled together at $`r_01/\omega `$ for $`l=0`$, and at their turning point $`r_\mathrm{t}1/\omega ^{\frac{4}{5}}`$ if $`l1`$. The continuity between the flow $`\gamma =5/3`$ and flows with $`\gamma `$ close to $`5/3`$ is checked by remarking that $`r_0r_\mathrm{s}`$ for $`\omega (53\gamma )^1`$ and $`r_tr_\mathrm{s}`$ for $`\omega (53\gamma )^{5/4}`$. Comparing this to $`\omega _l^{\mathrm{cut}}`$ in Eqs. (21) and (19), the cut-off frequency is simply interpreted as the maximum frequency such that the coupling of acoustic perturbations occurs in the subsonic flow.

## 4 Acoustic efficiency of entropy and vorticity perturbations

### 4.1 Analytical formulation

According to Eq. (6), it is natural to define two dimensionless complex acoustic efficiencies $`𝒬_S`$, $`𝒬_K`$ associated respectively to the two source terms $`\delta S_R`$ and $`\delta K_R`$. These definitions are made so that the acoustic flux $`F^{}`$ of an outgoing sound wave triggered by the perturbations of entropy and vorticity, measured far from the accretor is:

$$\frac{F^{}}{\dot{M}_0c_{\mathrm{}}^2}=\left|𝒬_S\delta S_R+𝒬_K\frac{\delta K_R}{c_{\mathrm{}}^2}\right|^2,$$
(24)
The integral expressions for $`𝒬_S`$ and $`𝒬_K`$ are derived in Appendix D:

$`𝒬_S`$ $`=`$ $`{\displaystyle \frac{i}{2\gamma \omega }}{\displaystyle _{r_\mathrm{s}}^{\mathrm{}}}\mathrm{e}^{^r\frac{i\omega \mathrm{d}r}{v(1^2)}}`$ (25)
$`{\displaystyle \frac{}{r}}\left[{\displaystyle \frac{v(1^2)^2}{^2}}{\displaystyle \frac{}{r}}\left(f_0\mathrm{e}^{i\omega ^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}\right)\right]\mathrm{d}r,`$
$`𝒬_K`$ $`=`$ $`{\displaystyle \frac{i}{2\omega }}{\displaystyle _{r_\mathrm{s}}^{\mathrm{}}}{\displaystyle \frac{f_0}{r^2v}}\mathrm{e}^{i\omega _R^r\frac{1+^2}{1^2}\frac{\mathrm{d}r}{v}}dr.`$ (26)
It is remarkable that both acoustic efficiencies are independent of the longitudinal number $`m`$. $`𝒬_S`$ and $`𝒬_K`$ depend only on $`\omega `$ and $`l`$.

### 4.2 Numerical calculations

The step of the numerical integration must be able to resolve both acoustic and entropy waves. Far from the accretor, the wavelength of the acoustic perturbation is asymptotically constant ($`2\pi c_{\mathrm{}}/\omega `$), whereas the wavelength of the entropy perturbation decreases to zero as $`2\pi v/\omega 1/r^2`$. It is therefore useful to transform the integrals in Eqs. (25) and (26) to accelerate their convergence. This is done in Appendix D by successive integrations by parts, thus obtaining Eqs. (129) and (142). If $`\gamma =5/3`$, the integrals are transformed near the origin in order to correctly treat the singularity.

### 4.3 Asymptotic expansion of $`|𝒬_S|`$ and $`|𝒬_K|`$ at high frequency for $`\gamma =5/3`$

Special properties of the Bondi flow $`\gamma =5/3`$ are related to the position of the sonic radius, at the central singularity. In particular, both cut-off frequencies $`\omega _0^{\mathrm{cut}}`$ and $`\omega _1^{\mathrm{cut}}`$ are infinite, which means that at any frequency, incoming sound waves are totally refracted outward:

$$|\overline{}_s|=1.$$
(27)
$`|𝒬_S|`$ is estimated in Appendix E.1 at high frequency for radial perturbations:

$$|𝒬_S|_{l=0}\frac{3}{5}\left(\frac{\omega }{6}\right)^{\frac{1}{3}}\mathrm{\Gamma }\left(\frac{2}{3}\right),$$
(28)
where the Gamma function satisfies $`\mathrm{\Gamma }(2/3)1.354`$.
The case of non radial perturbations is derived in Appendix E.2, where we define $`L^2l(l+1)`$:

$`|𝒬_S|_{l1}`$ $``$ $`{\displaystyle \frac{2}{5}}\left({\displaystyle \frac{\omega L^2}{2}}\right)^{\frac{1}{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right),`$ (29)
$`|𝒬_K|_{l1}`$ $``$ $`{\displaystyle \frac{2}{3}}\left({\displaystyle \frac{2}{\omega L^2}}\right)^{\frac{1}{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right).`$ (30)
The asymptotic expansions in Eqs. (28), (29) and (30), which involve no free scaling parameter, are successfully compared to numerical calculations of $`|𝒬_S|`$ and $`|𝒬_K|`$ in Figs. 6 and 7. The radial case of Eq. (28) is recovered by taking the continuous limit $`l0`$ in Eq. (29). Note that the asymptotic behaviour of the Bessel function $`K_{\frac{2}{3}}`$ for large arguments ($`L1`$),

$$K_{\frac{2}{3}}\left(\frac{2L}{3}\right)\left(\frac{3\pi }{4L}\right)^{\frac{1}{2}}\mathrm{e}^{\frac{2L}{3}},$$
(31)
makes the acoustic efficiencies decrease rapidly for high order $`l`$. Remembering that the frequency of the entropy perturbation is directly related to its radial wavelength ($`\lambda 2\pi v/\omega `$), our calculation of $`|𝒬_S|(\omega ,l)`$ indicates that the sound produced by the advection of an object of given mass in the Bondi flow depends a lot on its shape. Virtually no sound comes out if $`l`$ is large and $`\omega `$ is small (e.g. a stretched spaghetti). By contrast, the shape producing the highest acoustic flux corresponds to a small $`l`$ and a large $`\omega `$ (e.g. a wide thin saucer).
At high frequency, the coupling of entropy and vorticity perturbations to acoustic waves comes essentially from the region $`r_{\mathrm{eff}}1/\omega ^{2/3}`$ (i.e. the main contribution to the integrals in Eqs. 25 and 26, computed in Appendix E). The radial scaling at high frequency is therefore $`r_0r_\mathrm{t}r_{\mathrm{eff}}`$.

### 4.4 Case of accretion flows with $`\gamma `$ close to $`5/3`$

For $`\gamma `$ close to $`5/3`$ and at frequencies lower than the cut off frequency, we can use the approximation of $`|𝒬_S|`$ obtained for $`\gamma =5/3`$. According to Fig. 5, the maximum of $`|𝒬_S|`$ is reached at the frequency $`\omega _l^{\mathrm{max}}`$, with $`0.5<\omega _l^{\mathrm{max}}/\omega _l^{\mathrm{cut}}1`$. The asymptotic behaviour of the maximum value of $`|𝒬_S|`$ is therefore estimated by extrapolating Eqs. (28) and (29) to the cut-off frequency:

$`\mathrm{Max}|𝒬_S|_{l=0}`$ $``$ $`{\displaystyle \frac{1}{\gamma }}\left({\displaystyle \frac{\omega _0^{\mathrm{cut}}}{6}}\right)^{\frac{1}{3}}\mathrm{\Gamma }\left({\displaystyle \frac{2}{3}}\right),`$ (32)
$``$ $`r_\mathrm{s}^{\frac{1}{3}},`$ (33)
$`\mathrm{Max}|𝒬_S|_{l1}`$ $``$ $`{\displaystyle \frac{2^{\frac{1}{2}}}{3\gamma }}(\omega _1^{\mathrm{cut}})^{\frac{1}{3}}LK_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right),`$ (34)
$``$ $`LK_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right)r_\mathrm{s}^{\frac{5}{12}}.`$ (35)
It is interesting to note in Fig. 8 that the most efficient entropy perturbations at a given frequency when $`\gamma <5/3`$ are non radial ($`l=1`$), whereas radial ones are more efficient for $`\gamma =5/3`$ (Fig. 6). This is directly related to the fact that the critical refraction frequency is much higher for $`l=1`$ than for $`l=0`$ (Eqs. 19 and 21). This is illustrated in Fig. 9, where the acoustic efficiencies computed for $`c_\mathrm{s}/c_{\mathrm{}}=20`$ ($`\gamma =1.665`$) are very close to their values for $`\gamma =5/3`$, up to the refraction cut-off. In this figure, the efficiency $`|𝒬_S|`$ is higher for $`l=0`$ than for $`l=1`$ in the limited range of frequencies $`30<\omega <120`$, but the overall maximum is reached near $`\omega _1^{\mathrm{cut}}2300`$ for $`l=1`$.

## 5 Conclusions

### 5.1 Summary

The acoustic response of the Bondi flow to pressure, entropy and vorticity perturbations has been studied thoroughly. Let us summarize the main results:

(i) the critical refraction frequency for acoustic waves $`\omega _l^{\mathrm{cut}}`$ is quite different from the ”natural” frequency $`c_\mathrm{s}/r_\mathrm{s}`$ one might think of at first glance. Moreover $`\omega _l^{\mathrm{cut}}`$ is also very different for radial and non radial perturbations. The two fundamental frequencies $`\omega _0^{\mathrm{cut}},\omega _1^{\mathrm{cut}}`$ of this flow near the sonic point were accurately determined (Eqs. 19 and 21), and are ordered as follows: $`\omega _0^{\mathrm{cut}}<\omega _1^{\mathrm{cut}}<\omega _K`$, with

$`\omega _0^{\mathrm{cut}}`$ $``$ $`\mathrm{Max}\left\{(1^2){\displaystyle \frac{c}{r}}\right\}\left({\displaystyle \frac{}{\mathrm{log}r}}\right)_\mathrm{s}{\displaystyle \frac{c_\mathrm{s}}{r_\mathrm{s}}},`$ (36)
$`\omega _1^{\mathrm{cut}}`$ $``$ $`\mathrm{Max}\left\{(1^2)^{\frac{1}{2}}{\displaystyle \frac{c}{r}}\right\}\left({\displaystyle \frac{}{\mathrm{log}r}}\right)_\mathrm{s}^{\frac{1}{2}}{\displaystyle \frac{c_\mathrm{s}}{r_\mathrm{s}}}.`$ (37)
(ii) two quantities conserved in the linear approximation of the perturbed Bondi flow are source terms for acoustic waves: the entropy perturbation $`\delta S`$, and a quantity $`\delta K`$ related to both vorticity and entropy perturbations defined by Eq. (5).

(iii) the acoustic efficiencies $`𝒬_S`$ and $`𝒬_K`$ of entropy and vorticity perturbations are independent of the longitudinal number $`m`$.

(iv) $`|𝒬_S|`$ and $`|𝒬_K|`$ are highest for $`l=1`$ perturbations if $`\gamma <5/3`$.

(v) $`|𝒬_K|`$ is highest at relatively low frequency, and is always rather moderate ($`|𝒬_K|<1`$).

(vi) $`|𝒬_S|`$ is highest at high frequency near the refraction cut-off, and can be much larger than one if $`\gamma `$ is close to $`5/3`$, i.e. if $`c_\mathrm{s}/c_{\mathrm{}}1`$. This confirms the physical argument of FT2000 which stressed the importance of the enthalpy increase between the outer part of the flow and the sonic point.

(vii) asymptotic scalings of $`|𝒬_S|,|𝒬_K|`$ were obtained at high frequency for $`\gamma =5/3`$ (Eqs. 28 to 30), in excellent agreement with numerical calculations.
The analytical expressions obtained for $`\omega _l^{\mathrm{cut}}(\gamma )`$, and for $`|𝒬_S|`$ and $`|𝒬_K|`$ at high frequency for $`\gamma =5/3`$, are accurate enough (see Figs. 3, 6 and 7) to be used as tests of the accuracy of 3D hydrodynamical codes used for Bondi accretion, and particularly their limitations at high frequency.

These results suggest that if a shock were present as outer boundary, the entropic-acoustic cycle would be unstable with respect to high frequency waves if the sound speed at the sonic point is high enough,i.e. if $`\gamma `$ is close to $`5/3`$. Indeed, an entropy perturbation $`\delta S_1`$ advected towards the accretor triggers an outgoing acoustic flux $`F^{}=\dot{M}_0c^2|𝒬_S|^2|\delta S|^2`$. These acoustic waves propagate against the flow until they reach the shock, where they produce new entropy perturbations $`\delta S_2`$ (see FT2000), with $`\delta S_2>\delta S_1`$ since $`|𝒬_S|1`$. The eigenmodes corresponding to this unstable cycle will be investigated in the second paper of this series (Foglizzo 2001). The fact that the most unstable simulations of Bondi-Hoyle-Lyttleton accretion also correspond to $`\gamma =5/3`$ (Ruffert & Arnett 1994, Ruffert 1994b) already encourages us to look for an extrapolation of these results to the case of non-radial shocked accretion flows.

### 5.2 Discussion

#### 5.2.1 Upper bound for $`|𝒬_S|^2`$ imposed by the Schwarzschild radius

In the flow $`\gamma =5/3`$, the fact that the outgoing acoustic flux produced by entropy perturbations diverges at high frequency (Eqs. 28 and 29) is rather surprising, and directly related to the hypothesis of a point-like newtonian accretor ($`r_\mathrm{s}=0`$). In this section we show that $`|𝒬_S|^2`$ is naturally bounded if the finite size of the accretor is taken into account. Simple arguments are used in order to roughly estimate this upper bound in the case of an accreting black hole, without computing the exact relativistic corrections near its horizon.
Let us denote by $`r_{\mathrm{Sch}}=2GM/c_{\mathrm{light}}^2`$ the Schwarzschild radius of the black hole, where $`c_{\mathrm{light}}`$ is the speed of light. With a sound speed at infinity $`c_{\mathrm{}}c_{\mathrm{light}}`$, Eq. (1) guarantees that relativistic effects are negligible except for $`\gamma `$ close to $`5/3`$. Indeed Petrich et al. (1989) checked numerically that relativistic corrections are small for $`\gamma =1.1`$ and $`\gamma =4/3`$. If $`\gamma =5/3`$, the sonic point is bound to lie at a few Schwarzschild radii:

$$\eta \frac{r_\mathrm{s}}{r_{\mathrm{Sch}}}>1.$$
(38)
Our description of non radial acoustic waves in the flow $`\gamma =5/3`$ is not affected by relativistic corrections as long as their turning point $`r_t`$ is far enough from the horizon. Since $`r_t(L/\omega )^{4/5}`$ at high frequency (Eq. 122), a natural relativistic cut-off is introduced for $`r_t\eta r_{\mathrm{Sch}}`$:

$$\omega _{l1}^{\mathrm{cut}}\frac{0.25}{\eta ^{\frac{5}{4}}}\left[\frac{l(l+1)}{2}\right]^{\frac{1}{2}}\left(\frac{c_{\mathrm{}}}{c_{\mathrm{light}}}\right)^{\frac{1}{2}}\frac{c^3}{GM}.$$
(39)
The efficiency $`|𝒬_S|^2`$ of the generation of non radial acoustic waves close to this cut-off frequency is deduced from Eq. (29):

$$|𝒬_S|^2\frac{0.1}{\eta ^{\frac{5}{6}}}\left(\frac{L^2}{2}\right)^{\frac{7}{6}}\mathrm{e}^{\frac{4L}{3}}\left(\frac{c_{\mathrm{light}}}{c_{\mathrm{}}}\right)^{\frac{5}{3}}.$$
(40)
The numerical values of the cut-off frequency and maximum acoustic efficiency of the entropy perturbations $`l=1`$ are thus typically:

$`\omega _1^{\mathrm{cut}}`$ $``$ $`{\displaystyle \frac{100}{\eta ^{\frac{5}{4}}}}\left({\displaystyle \frac{c_{\mathrm{}}}{100\mathrm{km}/\mathrm{s}}}\right)^{\frac{1}{2}}\left({\displaystyle \frac{M}{10M_{\mathrm{sol}}}}\right)^1\mathrm{Hz},`$ (41)
$`|𝒬_S|_{\mathrm{max}}^2`$ $``$ $`{\displaystyle \frac{10^4}{\eta ^{\frac{5}{6}}}}\left({\displaystyle \frac{c_{\mathrm{}}}{100\mathrm{km}/\mathrm{s}}}\right)^{\frac{5}{3}}.`$ (42)
Therefore relativistic effects set an upper bound on the apparent divergence of $`|𝒬_S|`$ at high frequency for $`\gamma =5/3`$. Eq. (42) indicates that entropy perturbations may excite low degree acoustic waves very efficiently, up to non linear amplitudes if the component $`l=1`$ of the entropy perturbations exceeds a few percent. This efficiency is much smaller for inhomogeneities with a small angular scale, since $`|𝒬_S|`$ decreases exponentially when $`l`$ increases (Eq. 40).

#### 5.2.2 Non adiabatic limitations

Non adiabatic processes are discussed by Chang & Ostriker (1985) in the context of Bondi accretion. The efficiency $`|𝒬_S|`$ would be strongly affected by a strong thermal conduction which smoothes out entropy perturbations. This effect would be more pronounced for short wavelength perturbations, which are also the most efficient for $`\gamma `$ close to $`5/3`$. Without going into the details of such processes, our understanding of the acoustic efficiency allows the following remark: if the inner region of the flow were dominated by non adiabatic processes, the optimal frequency for $`|𝒬_S|`$ would be at most reduced to the maximum frequency of acoustic waves refracted inside the adiabatic region of the flow. Even when there is no shock, the mechanism of excitation of acoustic waves from the advection of entropy perturbations opens interesting perspectives to the problem of time delay in Cygnus X-1, approached numerically by Manmoto et al. (1996). Interpreting their results in terms of acoustic efficiency of entropy perturbations suggest that this physical process survives the inclusion of rotation, radiative cooling and viscous heating. The study of effects deserve a careful analysis which is beyond the scope of the present paper.

###### Acknowledgements.

The author acknowledges stimulating discussions with Michel Tagger. The Runge-Kutta algorithm used in this work for the numerical calculations was kindly provided by Roland Lehoucq.

## Appendix A Description of the unperturbed Bondi flow

### A.1 $`\gamma <5/3`$

The Bernoulli equation and the conservation of mass are written after normalizing velocities to the sound speed $`c_{\mathrm{}}`$ at infinity, and distances to the Bondi radius $`GM/c_{\mathrm{}}^2`$

$`{\displaystyle \frac{v^2}{2}}+{\displaystyle \frac{c^21}{\gamma 1}}`$ $`=`$ $`{\displaystyle \frac{1}{r}},`$ (43)
$`r^2vc^{\frac{2}{\gamma 1}}`$ $`=`$ $`\left({\displaystyle \frac{1}{2}}\right)^{\frac{\gamma +1}{2(\gamma 1)}}\left({\displaystyle \frac{4}{53\gamma }}\right)^{\frac{53\gamma }{2(\gamma 1)}}.`$ (44)
By taking the radial derivative of these equation, we obtain:

$`{\displaystyle \frac{\mathrm{log}c}{\mathrm{log}r}}`$ $`=`$ $`{\displaystyle \frac{\gamma 1}{2}}{\displaystyle \frac{^2}{1^2}}\left(2{\displaystyle \frac{1}{rv^2}}\right),`$ (45)
$`{\displaystyle \frac{\mathrm{log}v}{\mathrm{log}r}}`$ $`=`$ $`{\displaystyle \frac{1}{1^2}}\left({\displaystyle \frac{1}{rc^2}}2\right).`$ (46)
At the sonic point, $`r_\mathrm{s},c_\mathrm{s}`$ are defined by Eqs. (1) and 2). A Taylor expansion of Eqs. (45) and (46) leads to:

$`{\displaystyle \frac{\mathrm{log}v}{\mathrm{log}r}}(r_\mathrm{s})`$ $`=`$ $`{\displaystyle \frac{2}{\gamma +1}}\left[\gamma 1+\left({\displaystyle \frac{53\gamma }{2}}\right)^{\frac{1}{2}}\right],`$ (47)
$`{\displaystyle \frac{\mathrm{log}c}{\mathrm{log}r}}(r_\mathrm{s})`$ $`=`$ $`{\displaystyle \frac{\gamma 1}{\gamma +1}}\left[2\left({\displaystyle \frac{53\gamma }{2}}\right)^{\frac{1}{2}}\right],`$ (48)
$`{\displaystyle \frac{\mathrm{log}}{\mathrm{log}r}}(r_\mathrm{s})`$ $`=`$ $`\left({\displaystyle \frac{53\gamma }{2}}\right)^{\frac{1}{2}}.`$ (49)
### A.2 $`\gamma =5/3`$

In accretion flows where $`\gamma =5/3`$, Eqs. (43) and (44) become:

$`{\displaystyle \frac{v^2}{2}}+{\displaystyle \frac{3}{2}}(c^21)`$ $`=`$ $`{\displaystyle \frac{1}{r}},`$ (50)
$`r^2vc^3`$ $`=`$ $`{\displaystyle \frac{1}{4}}.`$ (51)
$`r,v,c`$ can be expressed as explicit functions of the variable $`x`$ defined by:

$`x`$ $``$ $`{\displaystyle \frac{1}{^{\frac{1}{2}}}},`$ (52)
$`{\displaystyle \frac{r}{x}}`$ $`=`$ $`{\displaystyle \frac{1^2}{2}},`$ (53)
$`r`$ $`=`$ $`{\displaystyle \frac{(x1)^2(9x^4+2x^2+1)}{6x^3(3x^22x+1)}},`$ (54)
$`v^2`$ $`=`$ $`{\displaystyle \frac{3}{(x1)^2}}{\displaystyle \frac{3x^22x+1}{9x^4+2x^2+1}},`$ (55)
$`c^2`$ $`=`$ $`{\displaystyle \frac{3x^4}{(x1)^2}}{\displaystyle \frac{3x^22x+1}{9x^4+2x^2+1}}.`$ (56)
## Appendix B Linearized equations for perturbations

### B.1 Second order differential system

The Euler equation is written as follows:

$$\frac{v}{t}+w\times v+\left(\frac{v^2}{2}+\frac{c^2}{\gamma 1}\frac{GM}{r}\right)=c^2\frac{S}{\gamma }.$$
(57)
Projecting the Euler equation onto the flow velocity, we obtain the equation of evolution of the Bernoulli constant:

$$\left(\frac{}{t}+v\right)\left(\frac{v^2}{2}+\frac{c^2}{\gamma 1}\frac{GM}{r}\right)=\frac{1}{\rho }\frac{p}{t}.$$
(58)
By combining the curl of Euler equation and the mass conservation, we obtain the equation for the evolution of the vorticity $`w`$:

$$\frac{}{t}\frac{w}{\rho }+(v)\frac{w}{\rho }=\left(\frac{w}{\rho }\right)v+\frac{1}{\rho }c^2\times \frac{S}{\gamma }.$$
(59)
We make a Fourier transform in time of the linearized equations describing the radial flow perturbed in spherical coordinated $`(r,\theta ,\phi )`$. The linearized equation of entropy conservation can be directly integrated:

$$\delta S=\delta S_R\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}.$$
(60)
The vorticity equation (59) can also be integrated when linearized:

$`w_r`$ $`=`$ $`\left({\displaystyle \frac{R}{r}}\right)^2(w_r)_R\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (61)
$`w_\theta `$ $`=`$ $`{\displaystyle \frac{1}{rv}}\left[Rv_R(w_\theta )_R{\displaystyle \frac{c^2c_R^2}{\mathrm{sin}\theta }}{\displaystyle \frac{}{\phi }}{\displaystyle \frac{\delta S_R}{\gamma }}\right]\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (62)
$`w_\phi `$ $`=`$ $`{\displaystyle \frac{1}{rv}}\left[Rv_R(w_\phi )_R+(c^2c_R^2){\displaystyle \frac{}{\theta }}{\displaystyle \frac{\delta S_R}{\gamma }}\right]\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}.`$ (63)
Thus the product $`r^2w_r`$ is conserved during advection, as remarked by Kovalenko & Eremin (1998). In order to write the linearized Euler equation in the simplest form, let us define the two functions $`f,g`$ as follows:

$`f`$ $``$ $`v\delta v_r+{\displaystyle \frac{2}{\gamma 1}}c\delta c,`$ (64)
$`g`$ $``$ $`{\displaystyle \frac{\delta v_r}{v}}+{\displaystyle \frac{2}{\gamma 1}}{\displaystyle \frac{\delta c}{c}}.`$ (65)
$`f`$ is the perturbation of the Bernoulli constant, which is directly related to the pressure variations according to Eq. (58). $`g`$ is related to the perturbation of the mass accretion rate and entropy through:

$$g=\frac{\delta \dot{M}}{\dot{M}}+\delta S.$$
(66)
By linearization of Eq. (57), we express the velocity components $`(\delta v_\theta ,\delta v_\phi )`$ as follows:

$`{\displaystyle \frac{\delta v_\theta }{v}}`$ $`=`$ $`{\displaystyle \frac{w_\phi }{i\omega }}+{\displaystyle \frac{1}{i\omega rv}}{\displaystyle \frac{}{\theta }}f{\displaystyle \frac{c^2}{i\omega rv}}{\displaystyle \frac{}{\theta }}{\displaystyle \frac{\delta S_R}{\gamma }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (67)
$`{\displaystyle \frac{\delta v_\phi }{v}}`$ $`=`$ $`{\displaystyle \frac{w_\theta }{i\omega }}+{\displaystyle \frac{1}{i\omega rv\mathrm{sin}\theta }}\left[{\displaystyle \frac{}{\phi }}fc^2{\displaystyle \frac{}{\phi }}{\displaystyle \frac{\delta S_R}{\gamma }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}\right].`$ (68)
The radial part of the Euler equation, together with the continuity equation, lead to the following differential system:

$`v{\displaystyle \frac{f}{r}}+{\displaystyle \frac{i\omega ^2f}{1^2}}`$ $`=`$ $`{\displaystyle \frac{i\omega v^2g}{1^2}}+i\omega c^2{\displaystyle \frac{\delta S_R}{\gamma }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (69)
$`v{\displaystyle \frac{g}{r}}+{\displaystyle \frac{i\omega ^2g}{1^2}}`$ $`=`$ $`{\displaystyle \frac{i\omega f}{c^2(1^2)}}+{\displaystyle \frac{i}{\omega }}\mathrm{\Delta }_{\theta ,\phi }f+{\displaystyle \frac{i\delta K_R}{r^2\omega }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$
where $`\mathrm{\Delta }_{\theta ,\phi }`$ is the non radial part of the Laplacian in spherical coordinates. The constant $`\delta K_R`$ is defined as follows:

$$\delta Kr^2\left[v(\times w)c^2\mathrm{\Delta }_{\theta ,\phi }\frac{\delta S}{\gamma }\right].$$
(71)
We deduce from Eqs. (62-63) that $`\delta K`$ is conserved when advected:

$$\left\{\frac{}{t}+v\frac{}{r}\right\}\delta K=0.$$
(72)
$`\delta K`$ can be written independently of the system of coordinates, as the first order perturbation of the quantity $`K`$ defined by:

$$K\frac{\dot{M}}{4\pi }\frac{𝐯}{\rho v}.\left[\times 𝐰c^2\times \left(\frac{𝐯}{v^2}\times \frac{S}{\gamma }\right)\right],$$
(73)
but the conservation of $`K`$ is established only in the linear approximation. By projecting $`f,g,\delta S_R`$ and the radial component of the curl of the vorticity $`(\times w_R)_r`$ onto the spherical harmonics $`Y_l^m(\theta ,\phi )`$, which are the eigenvectors of the Laplacian, we obtain:

$`v{\displaystyle \frac{f}{r}}+{\displaystyle \frac{i\omega ^2f}{1^2}}`$ $`=`$ $`{\displaystyle \frac{i\omega v^2g}{1^2}}+i\omega c^2{\displaystyle \frac{\delta S_R}{\gamma }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (74)
$`v{\displaystyle \frac{g}{r}}+{\displaystyle \frac{i\omega ^2g}{1^2}}`$ $`=`$ $`{\displaystyle \frac{i\omega f}{c^2(1^2)}}{\displaystyle \frac{iL^2}{\omega r^2}}f+{\displaystyle \frac{i\delta K_R}{r^2\omega }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}},`$ (75)
where we have defined $`L^2l(l+1)`$ for the sake of the simplicity of the equations.

### B.2 Second order differential equation

#### B.2.1 Pressure perturbation

The most natural function to describe the behaviour of acoustic waves is the pressure perturbation $`\delta p/p`$, which is related to the sound speed perturbation through:

$$\frac{\delta p}{\gamma p}=\frac{2}{\gamma 1}\frac{\delta c}{c}\frac{\delta S_R}{\gamma }\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}.$$
(76)
Using this equation in Eqs. (64-65) enables us to write $`(f,g)`$ as functions of $`\delta v_r`$ and $`\delta p`$ in Eqs. (74-75), and thus obtain Eq. (4) with:

$`a_1`$ $``$ $`{\displaystyle \frac{\mathrm{log}}{r}}{\displaystyle \frac{(c^2v^2)^2}{v\mathrm{\Delta }}}+{\displaystyle \frac{2i\omega v}{c^2v^2}},`$ (77)
$`a_0`$ $``$ $`{\displaystyle \frac{1}{c^2v^2}}\{\omega ^2L^2{\displaystyle \frac{c^2}{r^2}}{\displaystyle \frac{1}{v}}{\displaystyle \frac{}{r}}\left({\displaystyle \frac{c^4}{v}}{\displaystyle \frac{^2}{r}}\right)+`$ (78)
$`2i\omega {\displaystyle \frac{v}{r}}+c^2{\displaystyle \frac{\mathrm{log}\mathrm{\Delta }}{r}}[{\displaystyle \frac{\mathrm{log}^2}{r}}(1+^2){\displaystyle \frac{i\omega }{v}}]\},`$
$`b_0`$ $``$ $`{\displaystyle \frac{i\mathrm{\Delta }}{vc^2(1^2)}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}{\displaystyle \frac{}{r}}\left[{\displaystyle \frac{c^2}{\mathrm{\Delta }}}\left(\omega iv{\displaystyle \frac{\mathrm{log}}{r}}{\displaystyle \frac{1}{^2}}\right)\right],`$ (79)
$`b_1`$ $``$ $`{\displaystyle \frac{i\gamma \mathrm{\Delta }}{vc^2(1^2)}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v}}{\displaystyle \frac{}{r}}{\displaystyle \frac{v^2}{\omega r^2\mathrm{\Delta }}},`$ (80)
$`\mathrm{\Delta }`$ $``$ $`\omega ^2+2i\omega {\displaystyle \frac{v}{r}}+L^2{\displaystyle \frac{v^2}{r^2}}.`$ (81)
The function $`\mathrm{\Delta }`$ is introduced by our choice of writing the differential equation satisfied by the pressure perturbation, which is justified from the physical point of view. Note that $`\mathrm{\Delta }`$ never vanishes as long as the real part of the frequency is different from zero.

#### B.2.2 A more compact mathematical formulation

It is convenient to define ($`\stackrel{~}{f},\stackrel{~}{g}`$) as:

$`\stackrel{~}{f}`$ $``$ $`\mathrm{e}^{i\omega _R^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}f,`$ (82)
$`\stackrel{~}{g}`$ $``$ $`\mathrm{e}^{i\omega _R^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}g.`$ (83)
The differential system Eqs. (74-75) is then simpler:

$`{\displaystyle \frac{\stackrel{~}{f}}{r}}`$ $`=`$ $`{\displaystyle \frac{i\omega v\stackrel{~}{g}}{1^2}}+i\omega {\displaystyle \frac{c^2}{v}}{\displaystyle \frac{\delta S_R}{\gamma }}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}},`$ (84)
$`{\displaystyle \frac{\stackrel{~}{g}}{r}}`$ $`=`$ $`{\displaystyle \frac{i\stackrel{~}{f}}{\omega v}}\left[{\displaystyle \frac{\omega ^2}{c^2(1^2)}}{\displaystyle \frac{L^2}{r^2}}\right]+{\displaystyle \frac{i\delta K_R}{r^2\omega v}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}}.`$ (85)
The homogeneous differential equation ( 9) can be transformed into a more compact form using the new variable $`X`$:

$`{\displaystyle \frac{\mathrm{d}X}{\mathrm{d}r}}`$ $``$ $`{\displaystyle \frac{v}{1^2}},`$ (86)
$`W`$ $``$ $`{\displaystyle \frac{1}{v^2c^2}}(\omega ^2\omega _l^2),`$ (87)
where $`\omega _l`$ is defined by Eq. (7).

$`{\displaystyle \frac{^2\stackrel{~}{f}}{X^2}}+W\stackrel{~}{f}=`$
$`{\displaystyle \frac{1^2}{v}}\mathrm{e}^{i\omega {\scriptscriptstyle {\scriptscriptstyle \frac{\mathrm{d}X}{v^2}}}}\left\{{\displaystyle \frac{\omega }{^2}}{\displaystyle \frac{\delta S_R}{\gamma }}\left({\displaystyle \frac{\omega }{v}}+i{\displaystyle \frac{\mathrm{log}^2}{r}}\right)+{\displaystyle \frac{\delta K_R}{vr^2}}\right\}.`$
(88)
With this new variable, the sonic point $`r_\mathrm{s}`$ corresponds to $`X+\mathrm{}`$, while the spatial infinity corresponds to $`X=0`$.

## Appendix C Approximations of the homogeneous solution

### C.1 WKB approximation far from the accretor

The general solution $`f`$ of the homogeneous equation associated to Eq. (88) is a linear combination of ingoing ($`f_+`$) and outgoing ($`f_{}`$) waves, which can be approximated by the WKB method:

$$\stackrel{~}{f}^\pm \frac{\omega ^{\frac{1}{2}}}{W^{\frac{1}{4}}}\mathrm{exp}\left(\pm i^r\frac{vW^{\frac{1}{2}}}{1^2}dr\right).$$
(89)
This approximation is asymptotically valid for

$$\frac{\mathrm{log}W}{X}W^{\frac{1}{2}},$$
(90)
which is satisfied at high frequency or far from the accretor. The Wronskien $`𝒲`$ associated to the couple ($`\stackrel{~}{f}^+,\stackrel{~}{f}^{}`$) of solutions is:

$`𝒲`$ $``$ $`\stackrel{~}{f}^+{\displaystyle \frac{\stackrel{~}{f}^{}}{r}}\stackrel{~}{f}^{}{\displaystyle \frac{\stackrel{~}{f}^+}{r}},`$ (91)
$`=`$ $`{\displaystyle \frac{2i\omega v}{1^2}}.`$ (92)
The pressure perturbation $`\mathrm{\Delta }p`$ associated to a perturbation $`f`$ is deduced from Eq. (58):

$$\frac{\mathrm{\Delta }p}{\gamma p}=\frac{1}{c^2}\mathrm{e}^{i\omega ^r\frac{^2}{1^2}\frac{\mathrm{d}r}{v}}\left(\frac{\stackrel{~}{f}}{1^2}+\frac{iv}{\omega }\frac{\stackrel{~}{f}}{r}\right).$$
(93)
The normalization factor $`\omega ^{\frac{1}{2}}`$ in Eq. (89) is such that the pressure perturbations $`\mathrm{\Delta }p_\pm `$ associated to these solutions by Eq. (58) have the following asymptotic behaviour:

$$\frac{\mathrm{\Delta }p^\pm }{\gamma p}\frac{c_{\mathrm{}}}{c}\frac{^{\frac{1}{2}}}{1\pm }\mathrm{e}^{i\omega \frac{r}{c_{\mathrm{}}}}.$$
(94)
Both carry the same acoustic flux $`F^\pm =\dot{M}_0c_{\mathrm{}}^2`$ according to Eq. (10).

### C.2 Definition of the refraction coefficient $`\overline{}_s`$ and estimate at high frequency

#### C.2.1 $`\gamma <5/3`$

A Frobenius expansion in the vicinity of the sonic point ($`=1`$) shows that any solution of the homogeneous equation can be projected onto the couple $`f_1,f_2`$ of solutions, where $`f_1(r)`$ is regular and $`f_2(r)`$ is singular:

$`f_1(r)`$ $``$ $`1+𝒪(r),`$ (95)
$`f_2(r)`$ $``$ $`\left[1+𝒪(r)\right]\mathrm{exp}\left\{{\displaystyle \frac{i\omega }{2c_\mathrm{s}\dot{}_\mathrm{s}}}\mathrm{log}(rr_\mathrm{s})\right\}.`$ (96)
This defines a unique coefficient $`\overline{}_s`$ such that the solution

$$f_0f^++\overline{}_sf^{},$$
(97)
is regular at the sonic point. The coefficient $`|\overline{}_s|^2`$ is the refraction coefficient of sound waves. The linearity of the equations guarantees of course that the solution $`\mathrm{\Delta }p_0\mathrm{\Delta }p^++\overline{}_s\mathrm{\Delta }p^{}`$ is also regular.
An expansion of $`W`$ in the vicinity of the sonic point leads to:

$`W`$ $`=`$ $`q_0+q_1(rr_\mathrm{s})+𝒪(rr_\mathrm{s})^2,`$ (98)
$`q_0`$ $``$ $`{\displaystyle \frac{\omega ^2}{c_\mathrm{s}^4}},`$ (99)
$`q_1`$ $``$ $`{\displaystyle \frac{\omega ^2}{c_\mathrm{s}^4}}{\displaystyle \frac{\mathrm{log}c^2v^2}{r}}+{\displaystyle \frac{2L^2\dot{}_\mathrm{s}}{c_\mathrm{s}^2r_\mathrm{s}^2}},`$ (100)
$`rr_\mathrm{s}`$ $``$ $`\mathrm{exp}\left({\displaystyle \frac{2\dot{}_\mathrm{s}}{c_\mathrm{s}}}X\right).`$ (101)
The WKB approximation of $`\stackrel{~}{f}^\pm `$ is valid at high frequency $`\omega \omega _l^{\mathrm{cut}}`$, for $`rr_\mathrm{s}1`$:

$$\stackrel{~}{f}^\pm c_\mathrm{s}\mathrm{exp}\pm i\left[\frac{\omega }{2c_\mathrm{s}\dot{}_\mathrm{s}}\mathrm{log}(rr_\mathrm{s})+𝒪(1)\right].$$
(102)
The asymptotic solution to the differential equation (9) can also be approximated by a Bessel function in the vicinity of the sonic point, for $`rr_\mathrm{s}1`$. The regularity of $`f_0`$ at the sonic radius implies:

$`\stackrel{~}{f}_0`$ $``$ $`\lambda J_\nu \left[{\displaystyle \frac{c_\mathrm{s}q_1^{\frac{1}{2}}}{\dot{}_\mathrm{s}}}(rr_\mathrm{s})^{\frac{1}{2}}\right],`$ (103)
$`\nu `$ $``$ $`{\displaystyle \frac{ic_\mathrm{s}q_0^{\frac{1}{2}}}{\dot{}_\mathrm{s}}}={\displaystyle \frac{i\omega }{c_\mathrm{s}\dot{}_\mathrm{s}}}.`$ (104)
The Bessel and WKB approximations are both valid in the region where $`\omega ^1rr_\mathrm{s}1`$. In this region, the matching of the two approximations gives the normalization constant $`\lambda `$ and the refraction coefficient $`\overline{}_s`$. The asymptotic behavior of the Bessel function, when both the imaginary order $`\nu `$ and the argument tend to infinity, is obtained from Watson (1952, Chap. 8.6, p.262). The Bessel function $`J_\nu (z)`$ is written in terms of the Hankel functions $`H_\nu ^{(1)}(z)`$ and $`H_\nu ^{(2)}(z)`$. Let $`x,z`$ be real numbers, with $`|x|,z1`$, such that the order $`\nu `$ is purely imaginary $`\nu iz\mathrm{sinh}x`$:

$`J_\nu (z)`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left[H_\nu ^{(1)}(z)+H_\nu ^{(2)}(z)\right],`$ (105)
$`H_\nu ^{(1)}(z)`$ $``$ $`\left[{\displaystyle \frac{2\mathrm{exp}(\pi z\mathrm{sinh}x)}{\pi z\mathrm{cosh}x}}\right]^{\frac{1}{2}}\mathrm{e}^{+i\left[z(\mathrm{cosh}xx\mathrm{sinh}x)\frac{\pi }{4}\right]},`$ (106)
$`H_\nu ^{(2)}(z)`$ $``$ $`\left[{\displaystyle \frac{2\mathrm{exp}(\pi z\mathrm{sinh}x)}{\pi z\mathrm{cosh}x}}\right]^{\frac{1}{2}}\mathrm{e}^{i\left[z(\mathrm{cosh}xx\mathrm{sinh}x)\frac{\pi }{4}\right]}.`$ (107)
Applying these asymptotic expansions to Eq. (103), with

$`z`$ $``$ $`{\displaystyle \frac{c_\mathrm{s}q_1^{\frac{1}{2}}}{\dot{}_\mathrm{s}}}(rr_\mathrm{s})^{\frac{1}{2}}>0,`$ (108)
$`\mathrm{sinh}x`$ $`=`$ $`{\displaystyle \frac{\nu }{iz}}\left|(rr_\mathrm{s}){\displaystyle \frac{\mathrm{log}c^2v^2}{r}}\right|^{\frac{1}{2}}<0,`$ (109)
we obtain $`\lambda `$ and $`|f_0(r_\mathrm{s})|`$ to first order at high frequency:

$`\lambda `$ $`=`$ $`\left|{\displaystyle \frac{2\pi c_\mathrm{s}\omega }{\dot{}_\mathrm{s}}}\right|^{\frac{1}{2}}\mathrm{e}^{\frac{\omega \pi }{2c_\mathrm{s}\dot{}_\mathrm{s}}},`$ (110)
$`|f_0(r_\mathrm{s})|`$ $`=`$ $`c_\mathrm{s}.`$ (111)
The exponential decrease of $`|\overline{}_s|`$ at high frequency is deduced from Eqs. (106) and (107), by identifying the Hankel functions with the ingoing and outgoing acoustic waves:

$$|\overline{}_s|\mathrm{e}^{\frac{\omega \pi }{c_\mathrm{s}\dot{}_\mathrm{s}}}.$$
(112)
#### C.2.2 $`\gamma =5/3`$, $`l=0`$

The WKB approximation is valid in the range $`r1/\omega `$.

$$\stackrel{~}{f}^\pm c^{\frac{1}{2}}\mathrm{exp}\pm i\omega _0^r\frac{}{v(1^2)}dr,$$
(113)
where the normalization is chosen in agreement with Eq. (89). In the region $`r1/\omega `$ where the WKB approximation is not valid, we approximate the solution by the Bessel function $`J_0`$, which is regular at the origin:

$`{\displaystyle \frac{^2\stackrel{~}{f}_0}{r^2}}`$ $`+`$ $`{\displaystyle \frac{1}{r}}\left[1+𝒪(r^{\frac{1}{2}})\right]{\displaystyle \frac{\stackrel{~}{f}_0}{r}}+{\displaystyle \frac{\omega ^2}{8}}\stackrel{~}{f}_0\left[1+𝒪(r^{\frac{1}{2}})\right]=0,`$ (114)
$`\stackrel{~}{f}_0`$ $``$ $`\stackrel{~}{f}_0(0)J_0\left({\displaystyle \frac{\omega r}{2^{\frac{3}{2}}}}\right).`$ (115)
The other Bessel solution, $`Y_0`$, is singular with a logarithmic divergence, as in the case of $`\gamma <5/3`$. The Bessel approximation for $`r1/\omega `$ matches the WKB approximation:

$`\stackrel{~}{f}_0(0)J_0\left({\displaystyle \frac{\omega r}{2^{\frac{3}{2}}}}\right)`$ $``$ $`\stackrel{~}{f}_0(0)\left({\displaystyle \frac{2^{\frac{5}{2}}}{\pi \omega r}}\right)^{\frac{1}{2}}\mathrm{cos}\left({\displaystyle \frac{\omega r}{2^{\frac{3}{2}}}}{\displaystyle \frac{\pi }{4}}\right),`$ (116)
$``$ $`{\displaystyle \frac{1}{2^{\frac{1}{2}}\gamma r^{\frac{1}{2}}}}\left(\mathrm{e}^{i\frac{\omega r}{2^{\frac{3}{2}}}}+\overline{}_s\mathrm{e}^{+i\frac{\omega r}{2^{\frac{3}{2}}}}\right).`$ (117)
From this we deduce that:

$`\stackrel{~}{f}_0(0)`$ $`=`$ $`𝒪\left(\omega ^{\frac{1}{2}}\right),`$ (118)
$`|\overline{}_s|`$ $`=`$ $`1.`$ (119)
#### C.2.3 $`\gamma =5/3`$, $`l1`$

The case of non radial perturbations needs to be treated separately, because the behaviour of the solution is quite different near the sonic point. If $`l1`$, there is always a turning point beyond which the solution is evanescent, such that it decays exponentially on the sonic point. The leading behaviour of $`f_0,g_0`$ on this essential singularity is the following, within a multiplicative constant:

$`f_0`$ $``$ $`r^{\frac{1}{8}}\mathrm{exp}{\displaystyle \frac{2L}{r^{\frac{1}{4}}}},`$ (120)
$`g_0`$ $``$ $`2^{\frac{3}{2}}{\displaystyle \frac{iL}{\omega r^{\frac{1}{8}}}}\mathrm{exp}{\displaystyle \frac{2L}{r^{\frac{1}{4}}}}.`$ (121)
The solutions of the other branch diverge exponentially on the sonic point: as in the case of a regular sonic point, there is a unique combination of ingoing and outgoing waves such that the solution does not diverge at the sonic point. The turning point $`r_\mathrm{t}`$ corresponds to the radius where $`W=0`$ in Eq. (88):

$$r_\mathrm{t}\left(\frac{2L^2}{\omega ^2}\right)^{\frac{2}{5}}.$$
(122)
The WKB approximation is valid for $`r_\mathrm{t}r1`$:

$`\stackrel{~}{f}_\pm `$ $``$ $`{\displaystyle \frac{c^{\frac{1}{2}}}{\left[1\frac{L^2c^2}{r^2\omega ^2}(1^2)\right]^{\frac{1}{4}}}}\times `$ (123)
$`\mathrm{exp}\pm i\omega {\displaystyle _{r_0}^r}{\displaystyle \frac{}{v(1^2)}}\left[1{\displaystyle \frac{L^2c^2}{r^2\omega ^2}}(1^2)\right]^{\frac{1}{2}}dr.`$
## Appendix D Acoustic efficiencies $`𝒬_S,𝒬_K`$ of entropy and vorticity perturbations

### D.1 Acoustic efficiency $`𝒬_K`$

The outgoing pressure perturbation produced by the advection of a perturbation $`\delta K_R`$ is computed by solving Eq. (88), imposing the regularity at the sonic point and the absence of an incoming sound wave from infinity. The general solution of Eq. (88) with $`\delta S_R=0`$ can be written as:

$`\stackrel{~}{f}`$ $`=`$ $`{\displaystyle \frac{i\delta K_R}{2\omega }}[\stackrel{~}{f}^{}\{\beta +{\displaystyle _{r_0}^r}{\displaystyle \frac{\stackrel{~}{f}^+}{r^2v}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\}`$ (124)
$`\stackrel{~}{f}^+\{\alpha +{\displaystyle _{r_0}^r}{\displaystyle \frac{\stackrel{~}{f}^{}}{r^2v}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\}],`$
$`{\displaystyle \frac{\delta p}{\gamma p}}`$ $`=`$ $`{\displaystyle \frac{i\delta K_R}{2\omega }}[{\displaystyle \frac{\mathrm{\Delta }p^{}}{\gamma p}}\{\beta +{\displaystyle _{r_0}^r}{\displaystyle \frac{\stackrel{~}{f}^+}{r^2v}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\}`$ (125)
$`{\displaystyle \frac{\mathrm{\Delta }p^+}{\gamma p}}\{\alpha +{\displaystyle _{r_0}^r}{\displaystyle \frac{\stackrel{~}{f}^{}}{r^2v}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\}].`$
According to Eqs. (94) and (95-96), the integrals involved in Eq. (125) are well defined both at infinity and at the sonic radius. The condition that no sound wave comes from infinity leads to:

$$\alpha =_{r_0}^{\mathrm{}}\frac{\stackrel{~}{f}^{}}{r^2v}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}.$$
(126)
The regularity at the sonic radius can be written:

$$\beta =\overline{}_s\alpha _{r_0}^{r_\mathrm{s}}\frac{\stackrel{~}{f}_0}{r^2v}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}.$$
(127)
Thus the efficiency of sound wave emission by the advection of vorticity perturbations is:

$$𝒬_K=\frac{i}{2\omega }_{r_\mathrm{s}}^{\mathrm{}}\mathrm{e}^{i\omega _R^r\frac{1+^2}{1^2}\frac{\mathrm{d}r}{v}}\frac{f_0}{r^2v}dr.$$
(128)
This integral converges as $`1/r^3`$ at infinity. It can be transformed using integrations by parts together with the homogeneous differential system satisfied by $`f_0,g_0`$, so that it converges rapidly at infinity, thus allowing an easier numerical computation. The following expression converges as $`1/r^5`$:

$`𝒬_K={\displaystyle \frac{1}{2\omega ^2}}`$ $`{\displaystyle _{r_\mathrm{s}}^{\mathrm{}}}\mathrm{e}^{i\omega _R^r\frac{1+^2}{1^2}\frac{\mathrm{d}r}{v}}`$ (129)
$`\times \left\{{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{r^2}}\right)f_0+{\displaystyle \frac{i\omega v}{r^2}}g_0\right\}\mathrm{d}r.`$
### D.2 Acoustic efficiency $`𝒬_S`$

The outgoing pressure perturbation produced by an incoming entropy perturbation $`\delta S_R`$ is also computed by solving Eq. (88), imposing the regularity at the sonic point and the absence of an incoming sound wave from infinity. The general solution of Eq. (88) with $`\delta K_R=0`$ can be written as:

$`\stackrel{~}{f}={\displaystyle \frac{\delta S_R}{2\gamma }}\times `$
$`[\stackrel{~}{f}^{}\{\beta +{\displaystyle _{r_0}^r}\stackrel{~}{f}^+{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{^2}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\right)\}`$
$`\stackrel{~}{f}^+\{\alpha +{\displaystyle _{r_0}^r}\stackrel{~}{f}^{}{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{^2}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\right)\}],`$ (130)
$`{\displaystyle \frac{\delta p}{\gamma p}}={\displaystyle \frac{\delta S_R}{2\gamma }}\times `$
$`[{\displaystyle \frac{\mathrm{\Delta }p^{}}{\gamma p}}\{\beta +{\displaystyle _{r_0}^r}\stackrel{~}{f}^+{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{^2}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\right)\}`$
$`{\displaystyle \frac{\mathrm{\Delta }p^+}{\gamma p}}\{\alpha +{\displaystyle _{r_0}^r}\stackrel{~}{f}^{}{\displaystyle \frac{}{r}}\left({\displaystyle \frac{1^2}{^2}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\right)\}].`$ (131)
Since the integrals of Eq. (131) converge at the sonic point, the regularity of the solution at the sonic point can be written:

$$\beta =\alpha \overline{}_s_{r_0}^r\stackrel{~}{f}_0\frac{}{r}\left(\frac{1^2}{^2}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\right).$$
(132)
In Eq. (131) the integrals diverge at infinity. Two integrations by parts leads to define the two functions $`F_\pm ,H_\pm `$ as follows:

$`F^\pm `$ $``$ $`c^2{\displaystyle \frac{(1^2)^2}{^2}}{\displaystyle \frac{\mathrm{\Delta }p^\pm }{\gamma p}}\mathrm{e}^{i\omega ^r\frac{1+^2}{1^2}\frac{\mathrm{d}r}{v}},`$ (133)
$`H^\pm `$ $``$ $`\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}{\displaystyle \frac{}{r}}\left[{\displaystyle \frac{v(1^2)^2}{i\omega ^2}}{\displaystyle \frac{\stackrel{~}{f}^\pm }{r}}\right].`$ (134)
The pressure perturbation associated to the entropy perturbation $`\delta S_R`$ is:

$`{\displaystyle \frac{\delta p}{\gamma p}}`$ $`=`$ $`{\displaystyle \frac{\delta S_R}{2\gamma }}[{\displaystyle \frac{\mathrm{\Delta }p^{}}{\gamma p}}\{\beta +[F^+(r_0)+{\displaystyle _{r_0}^r}H^+\mathrm{d}r]\}`$ (135)
$`+`$ $`{\displaystyle \frac{\mathrm{\Delta }p^+}{\gamma p}}\{\alpha [F^{}(r_0)+{\displaystyle _{r_0}^r}H^{}\mathrm{d}r]\}].`$
Far from the accretor, the pressure perturbation can be projected onto the ingoing and outgoing sound waves. $`\alpha `$ is chosen such that there is no sound wave coming from infinity:

$$\alpha =F^{}(r_0)+_{r_0}^{\mathrm{}}H^{}dr.$$
(136)
Using Eqs. (132) and (136), and the fact that $`F_0(r_\mathrm{s})=0`$, we take the limit $`r\mathrm{}`$, $`r_0r_\mathrm{s}`$ in Eq. (135) and obtain:

$$𝒬_S=\frac{i}{2\gamma \omega }_{r_\mathrm{s}}^{\mathrm{}}\mathrm{e}^{i\omega ^r\frac{\mathrm{d}r}{v(1^2)}}\frac{}{r}\left[\frac{v(1^2)^2}{^2}\frac{\stackrel{~}{f}_0}{r}\right]dr.$$
(137)
This integral can be written with the following format, using the homogeneous differential system (Eqs. 84, 85) satisfied by $`f_0,g_0`$:

$$𝒬_S=\frac{1}{2\gamma }_{r_\mathrm{s}}^{\mathrm{}}\mathrm{e}^{i\omega ^r\frac{1+^2}{1^2}\frac{\mathrm{d}r}{v}}(A_kf_0+B_kg_0)dr,$$
(138)
where the functions $`A_k(r),B_k(r)`$ depend only on the unperturbed flow, and the integral converges as $`1/r^k`$ far from the accretor:

$`A_1`$ $``$ $`{\displaystyle \frac{i\omega }{v}}\left(1{\displaystyle \frac{\omega _l^2}{\omega ^2}}\right),`$ (139)
$`B_1`$ $``$ $`{\displaystyle \frac{}{r}}(c^2v^2).`$ (140)
Two successive integrations by parts make the convergence of the integral much faster, thus allowing an easier numerical calculation:

$`A_5`$ $``$ $`{\displaystyle \frac{}{r}}\left[(1^2)\left(1{\displaystyle \frac{\omega _L^2}{\omega ^2}}\right)\right]`$ (141)
$`+`$ $`{\displaystyle \frac{1}{c^2}}\left(1{\displaystyle \frac{\omega _L^2}{\omega ^2}}\right){\displaystyle \frac{}{r}}(c^2v^2){\displaystyle \frac{i\omega v}{c^2}}\left(1{\displaystyle \frac{\omega _L^2}{\omega ^2}}\right)^2,`$
$`B_5`$ $``$ $`{\displaystyle \frac{}{r}}\left\{(1^2)\left[{\displaystyle \frac{v}{i\omega }}{\displaystyle \frac{}{r}}(c^2v^2)v^2\left(1{\displaystyle \frac{\omega _L^2}{\omega ^2}}\right)\right]\right\}.`$ (142)
## Appendix E Asymptotic estimate of $`|𝒬_S|,|𝒬_K|`$ at high frequency for $`\gamma =5/3`$

Eq.(137) can be estimated in the limit $`\omega c_{\mathrm{}}^3/GM`$, when $`\gamma =5/3`$. We treat the case of radial ($`l=0`$) and non radial ($`l1`$) perturbations separately.

### E.1 Radial perturbations $`l=0`$

We use a local expansion of $``$ where $`r1`$:

$$^2=14r^{\frac{1}{2}}+𝒪(r).$$
(143)
Writing $`\stackrel{~}{f}_0=\stackrel{~}{f}_{}+\overline{}_s\stackrel{~}{f}_+`$, we obtain for $`r_11/\omega `$ and $`r_21`$:

$`{\displaystyle _{r_1}^{r_2}}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}}{\displaystyle \frac{}{r}}\left({\displaystyle \frac{v(1^2)^2}{^2}}{\displaystyle \frac{\stackrel{~}{f}_0}{r}}\right)dr`$
$`{\displaystyle \frac{6^{\frac{2}{3}}}{5}}\omega ^{\frac{1}{3}}\overline{}_s{\displaystyle _{u_1}^{u_2}}{\displaystyle \frac{\mathrm{e}^{iu}}{u^{\frac{1}{3}}}}du,`$
$`u{\displaystyle \frac{2^{\frac{1}{2}}r^{\frac{3}{2}}\omega }{3}}.`$ (144)
The range of validity $`1/\omega r1`$ implies that $`u_11/\omega ^{1/2}`$ and $`u_2\omega `$. The contribution of the region $`r_2r\mathrm{}`$ is negligible.

$$_{r_1}^{r_2}\mathrm{e}^{i\omega _R^r\frac{\mathrm{d}r}{v(1^2)}}\frac{}{r}\left(\frac{v(1^2)^2}{^2}\frac{\stackrel{~}{f}_0}{r}\right)dr\omega ^{\frac{1}{3}}.$$
(145)
The contribution of the region $`rr_1`$ is estimated using the asymptotic behaviour of the Bessel function $`J_1`$:

$`{\displaystyle \frac{}{z}}\left({\displaystyle \frac{1^2}{^2}}{\displaystyle \frac{\stackrel{~}{f}_0}{z}}\right)2^{\frac{11}{4}}{\displaystyle \frac{\stackrel{~}{f}_0(0)}{\omega ^{\frac{1}{2}}}}{\displaystyle \frac{}{z}}\left[z^{\frac{1}{2}}J_1(z)\right],`$ (146)
$`{\displaystyle _0^{z_1}}\mathrm{e}^{iz}{\displaystyle \frac{}{z}}\left({\displaystyle \frac{1^2}{^2}}{\displaystyle \frac{\stackrel{~}{f}}{z}}\right)dz{\displaystyle \frac{\stackrel{~}{f}_0(0)}{\omega ^{\frac{1}{2}}}}{\displaystyle \frac{2^{\frac{9}{4}}}{\pi ^{\frac{1}{2}}}}\mathrm{e}^{i\frac{\pi }{4}}z_1`$ (147)
Thus we deduce the acoustic efficiency at high frequency:

$$𝒬_Si\mathrm{e}^{iz_R}\frac{\overline{}_s}{\gamma }\left(\frac{\omega }{6}\right)^{\frac{1}{3}}_0^{\mathrm{}}\frac{\mathrm{e}^{iu}}{u^{\frac{1}{3}}}du.$$
(148)
The main contribution to the integral comes from the region where $`u1`$, i.e. $`r_{\mathrm{eff}}(6/\omega )^{2/3}/2`$. Eq. (148) is rewritten in Eq. (28) using the Gamma function.

### E.2 Non radial perturbations $`l1`$

As in the case of radial perturbations, the main contribution to the integral comes from the region accessible to the WKB approximation. The outgoing wave is the main contributor to the integral (137). Using the asymptotic development of $`v,c`$ for $`r1`$ in Appendix. A.2,

$`{\displaystyle _{r_0}^r}\left\{1\left[1{\displaystyle \frac{L^2c^2}{r^2\omega ^2}}(1^2)\right]^{\frac{1}{2}}\right\}{\displaystyle \frac{\mathrm{d}r}{v(1^2)}}`$
$`{\displaystyle \frac{1}{3}}\left(2^{\frac{1}{2}}r^{\frac{3}{2}}\omega +{\displaystyle \frac{L^2}{2^{\frac{1}{2}}r^{\frac{3}{2}}\omega }}\right).`$ (149)
With $`u2^{\frac{1}{2}}r^{\frac{3}{2}}\omega /3`$, the acoustic efficiencies at high frequency are deduced from Eqs. (123), (128), (137) and (149):

$`𝒬_S`$ $``$ $`i{\displaystyle \frac{\overline{}_s}{\gamma }}\left({\displaystyle \frac{\omega }{6}}\right)^{\frac{1}{3}}{\displaystyle _0^{\mathrm{}}}\mathrm{exp}i\left(u{\displaystyle \frac{L^2}{9u}}\right){\displaystyle \frac{\mathrm{d}u}{u^{\frac{1}{3}}}},`$ (150)
$`𝒬_K`$ $``$ $`i{\displaystyle \frac{2^{\frac{1}{3}}\overline{}_s}{3^{\frac{5}{3}}\omega ^{\frac{1}{3}}}}{\displaystyle _0^{\mathrm{}}}\mathrm{exp}i\left(u{\displaystyle \frac{L^2}{9u}}\right){\displaystyle \frac{\mathrm{d}u}{u^{\frac{5}{3}}}}.`$ (151)
Using the new variable

$$t\frac{3}{2L}\left[u\frac{L^2}{9u}\right],$$
(152)
the integrals in Eqs. (150-151) can be written using $`K_{\frac{2}{3}}`$, the modified Bessel function of order $`2/3`$ (Gradshteyn & Ryzhik 1980, p. 430):

$`{\displaystyle _0^{\mathrm{}}}\mathrm{exp}i\left(u{\displaystyle \frac{L^2}{9u}}\right){\displaystyle \frac{\mathrm{d}u}{u^{\frac{1}{3}}}}`$ $`=`$ $`2\left({\displaystyle \frac{L}{3}}\right)^{\frac{2}{3}}\mathrm{e}^{i\frac{\pi }{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right),`$ (153)
$`{\displaystyle _0^{\mathrm{}}}\mathrm{exp}i\left(u{\displaystyle \frac{L^2}{9u}}\right){\displaystyle \frac{\mathrm{d}u}{u^{\frac{5}{3}}}}`$ $`=`$ $`2\left({\displaystyle \frac{3}{L}}\right)^{\frac{2}{3}}\mathrm{e}^{+i\frac{\pi }{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right).`$ (154)
Thus we obtain:

$`𝒬_S`$ $``$ $`{\displaystyle \frac{2i}{3\gamma }}\overline{}_s\left({\displaystyle \frac{\omega L^2}{2}}\right)^{\frac{1}{3}}\mathrm{e}^{\frac{i\pi }{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right),`$ (155)
$`𝒬_K`$ $``$ $`{\displaystyle \frac{2i}{3}}\overline{}_s\left({\displaystyle \frac{2}{\omega L^2}}\right)^{\frac{1}{3}}\mathrm{e}^{+\frac{i\pi }{3}}K_{\frac{2}{3}}\left({\displaystyle \frac{2L}{3}}\right).`$ (156)
Here again, the main contribution to the integrals comes from the region $`r_{\mathrm{eff}}(6/\omega )^{2/3}/2`$. By contrast with $`|𝒬_S|`$, $`|𝒬_K|`$ decreases at high frequency.


format: text

text:
# References

In this article considered long-range interactions of the ball under assumption that the unit of the its volume interact with probe particle as $`A/r^n`$. There are many examples of the such singular potentials: e.g. Van-der-Vaals interactions. Also must be noticed that in some of the theories with higher dimensions (see e.g. , and references therein) appear the following corrections to Newton potential :

$$U(r)=\frac{Gm_1m_2}{r}(1+A/r^n)$$
(1)
Previously has been considered the following corrections to Newton potential which also appear in some of multidimensional theories (see references in ):

$$U(r)=\frac{Gm_1m_2}{r}(1+se^{Mr})$$
(2)
In number of papers (see references in ) has been calculated the potential of the ball with homogeneous density under assumption that any particle of the ball interact with probe particle in accordance with formula (2). In  has been also published the experimental restriction on parameters $`s`$ and$`M`$. Analogously, using data of the experiments described in  in principle possible to obtain some restriction on $`A`$ and $`n`$.

We present also model-independent corrections to the energy levels of the electron in the Coulomb field of the nuclei from regularized version of the potential $`A/r^n`$.

We obtain for potential of ball with constant density :

$$U_n(r)=d^3r_1\frac{A}{|\stackrel{}{r}\stackrel{}{r^{}}|^n}=\frac{2\pi A}{r}r^{}𝑑r^{}(\frac{|rr^{}|^{(2n)}}{2n}\frac{(r+r^{})^{(2n)}}{2n})$$
(3)
At $`r>R`$ we have after integration:

$$U_n(r)=\frac{2\pi A}{r}\frac{1}{n2}[\frac{(rR)^{4n}}{4n}\frac{r(rR)^{3n}}{3n}\frac{(r+R)^{4n}}{4n}+\frac{r(r+R)^{3n}}{3n}]$$
(4)
From (17) it is seen that potential is singular at $`rR`$ ($`U(r)|rR|^{(4n)}`$.

From this formula it is seen that cases $`n=2,3,4`$ must be considered separately.We obtain:

$$U_2(r)=\frac{2\pi A}{r}(rR+\frac{1}{2}(R^2r^2)log(\frac{r+R}{|rR|}),$$
(5)
$$U_3(r)=\frac{2\pi A}{r}[rlog(\frac{r+R}{rR})2R],$$
(6)
$$U_4(r)=\frac{\pi A}{r}[log(\frac{(rR)}{r+R})+\frac{2rR}{r^2R^2}].$$
(7)
During derivation of $`U_2(r)`$ has been taken into account potential $`U_n(r)`$ at $`r<R`$ (see formula (11) below). From (4) for $`n=5,6,7`$ we obtain the following simplifications:

$$U_5(r)=\frac{4\pi AR^3}{3r}\frac{1}{(r^2R^2)^2}$$
(8)
$$U_6(r)=\frac{4\pi AR^3}{3}\frac{1}{(r^2R^2)^3}$$
(9)
$$U_7(r)=\frac{4\pi AR^3}{15r}\frac{(5r^2+R^2)}{(r^2R^2)^4}$$
(10)
At $`n=1`$ we obtain as it must be $`U_1(r)=4\pi R^3A/3r`$ . At $`r<R`$ and $`n<3`$ we have:

$$U_n(r)=\frac{2\pi A}{r}\frac{1}{n2}[\frac{(Rr)^{4n}}{4n}+\frac{r(Rr)^{3n}}{3n}\frac{(r+R)^{4n}}{4n}+\frac{r(r+R)^{3n}}{3n}]$$
(11)
At $`n3`$ integral $`r^{}𝑑r^{}\frac{1}{(rr^{})^{(2n)}}`$ is divergent if $`r<R`$.

It mean that inside ball ($`r<R`$) and $`n3`$ we must regularised interaction $`A/r^n`$. We choose the following regularizations:

$$V_n(r)=\frac{A}{(r^2+a^2)^{n/2}}$$
(12)
$$V_n(r)=\frac{A}{(r+a)^n},$$
(13)
$$V_n(r)=\frac{A}{r(r+r_0)^n}.$$
(14)
where$`r_0,a<<R`$. For example $`a`$ may be a size of the molecules (atoms) which contained in this ball.The last potential transformed into Coulomb potential at $`r<<r_0`$.

Using regularized potential (12) we obtain for $`n=5`$ the following result:

$$V_5(r)=\frac{2\pi A}{5r}[f_{}+f_++\frac{r}{a^2}((Rr)f_{}+(R+r)f_+)]$$
(15)
where $`f_\pm =((R\pm r)^2+a^2)^{\frac{1}{2}}`$ At $`r>R`$ and $`|rR|>>a`$ we again obtain $`A/r^n`$, and at $`r<R`$ ,$`r,|rR|>>a`$ we have:

$$V_5(r)=\frac{4\pi A}{5}[\frac{1}{a^2}\frac{1}{R^2r^2}].$$
(16)
Now we consider in model independent way (the nature of correction may be any including cases which has been considered above) corrections to the energy levels of the electron in the Coulomb field of the nuclei from regularised version of the potential $`A/r^n`$ expressed by formula (14):

$$\delta E=d^3r|\psi |^2A/(r(r+r_0)^n)4Aa_B^3r_0^{2n}\frac{2n3}{(n2)(n1)}$$
(17)
Here has been taken into account that at $`r_0<<a_B`$ in integral we can put $`\psi ^2=1/\pi a_B^3`$ at $`l=0`$. Analogously, in case of potential (13) we obtain:

$$\delta E=d^3r|\psi |^2A/((r+r_0)^n)8Aa_B^3r_0^{3n}\frac{1}{(n3)(n2)(n1)}.$$
(18)
This formulas are valid if $`n>3`$. At $`n=3`$ by consideration of two range in integral ($`0<r<L`$, $`L<r<\mathrm{}`$ where $`r_0<<L<<a_B`$ ) e.g. from (22) we obtain:

$$\delta E=d^3r|\psi |^2A/((r+r_0)^3)4\pi A(1/32+log(a_B/2r_0)C)$$
(19)
where $`C`$-is Eiler constant.

The author express his sincere gratitude to G.V.Grigoryan for helpful discussions.


text:
# On a character sum problem of H. Cohn

## 1. Introduction

Let $`p`$ be prime and let $`F_{p^k}`$ be the finite field with $`p^k`$ elements. Let $`f:F_{p^k}^\times 𝐂`$ be a nontrivial multiplicative character, and extend $`f`$ to a function on $`F_{p^k}`$ by letting $`f(0)=0`$. It is then easy to see that the following holds:

(1.1)
$$\underset{xF_{p^k}}{}f(x)\overline{f(x+h)}=\{\begin{array}{cc}1\hfill & \text{if }h0\hfill \\ p^k1\hfill & \text{if }h=0\hfill \end{array}$$

Cohn asked (see p. 202 in ) if the converse is true in the following sense: if a function $`f:F_{p^k}𝐂`$ satisfies

(1.2)
$$f(0)=0,f(1)=1,\text{and }|f(x)|=1\text{ for }x0$$

and equation 1.1, does it follow that $`f`$ is a multiplicative character?

The problem has recently received some attention. In , Choi and Siu proved that the converse is not true for $`k>1`$. One of the arguments given is quite pretty, and proceeds as follows: Let $`\lambda `$ be a linear automorphism of $`F_{p^k}`$ so that $`\lambda (1)=1`$. If $`f`$ satisfies 1.1 and 1.2, so does $`f`$ composed with $`\lambda `$. Now, if $`f`$ is an injective multiplicative character then the converse being true implies that $`f`$ composed with $`\lambda `$ must be an injective multiplicative character. On the other hand, a simple counting argument shows that the number of possible $`\lambda `$’s is greater than the number of injective characters.

However, the case $`k=1`$ remains unresolved. In , Biro proved that there are only finitely many functions satisfying equation 1.1 and 1.2 for each $`p`$. Biro also solved the following “characteristic $`p`$” version of the problem (, Theorem 2):

###### Theorem (Biro).

Let $`p`$ be a prime, let $`F_p`$ be the finite field with $`p`$ elements, and $`FF_p`$ any field of characteristic $`p`$. Assume that there is given an $`a_iF`$ for every $`iF_p`$ such that $`a_0=0,a_1=1,a_i0`$ for $`i0`$, and

$$\underset{iF_p^\times }{}\frac{a_{i+j}}{a_i}=1$$

for every $`jF_p^\times `$. Then $`a_i=i^A`$ for every $`iF_p`$ with some $`1Ap2`$.

Using this Biro deduces that the converse holds for functions taking values in $`\{1,0,1\}`$.111There appears to be several independent proofs of this result, see the introduction in . In fact, if $`m`$ is coprime to $`p`$, then the case of the nonzero values of $`f`$ being $`m`$-th roots of unity can be deduced in a similar way: Let $`𝔒`$ be the ring of integers in $`𝐐(e^{2\pi i/m})`$, and let $`𝔓𝔒`$ be a prime ideal lying above $`p`$. The result then follows from the theorem by letting $`F=𝔒/𝔓`$ and noting that $`m`$-th roots of unity are distinct modulo $`p`$. (Since $`|f(x)|=1`$ for $`x0`$ we have $`\overline{f(x)}=1/f(x)`$.)

The aim of this paper is to show that the converse is true for the case $`k=1`$, under the additional assumption that the nonzero values of $`f:F_p𝐂`$ are $`m`$-th roots of unity, including the case $`p|m`$. We begin by giving a proof that does not depend on Biro’s result for the case $`(m,p)=1`$, and we then show how to modify the argument for the general case.

Acknowledgements: I would like to thank Ernest Croot, Andrew Granville, Robert Rumely, and Mark Watkins for helpful and stimulating discussions. I would also like to thank the referee for several suggestions on how to improve the exposition, and for pointing out that the case $`p|m`$ can be deduced independently of Biro’s theorem.

## 2. Preliminaries

In what follows we assume that $`p`$ is odd since the case $`p=2`$ is trivial.

We will use the following conventions: if a function $`f`$ takes values in $`𝐂`$ and $`\sigma \mathrm{Aut}(𝐂/𝐐)`$, then we let $`f^\sigma `$ be the function defined by $`f^\sigma (x)=\sigma (f(x))`$. We regard $`\psi (x)=e^{2\pi ix/p}`$ as a nontrivial additive character of $`F_p`$. For an integer $`t`$, $`\psi _t`$ will denote the character $`\psi _t(x)=\psi (tx)`$. By $`\zeta _m`$ we denote the $`m`$-th root of unity $`\zeta _m=e^{2\pi i/m}`$.

Let $`m`$ be even and large enough so that all nonzero values of $`f`$ are $`m`$-th roots of unity, and write $`m=np^k`$, where $`(n,p)=1`$. Let $`K=𝐐(\zeta _n)`$, $`L=K(\zeta _p,\zeta _{p^k})`$, and let $`G=\mathrm{Gal}(L/𝐐)`$, $`H=\mathrm{Gal}(L/K)`$ denote the Galois groups of $`L/𝐐`$ and $`L/K`$. By $`𝔒_K`$ and $`𝔒_L`$ we will denote the ring of integers in $`K`$ respectively $`L`$.

The “Gauss sum”

$$G(f,\psi )=\underset{x=0}{\overset{p1}{}}f(x)\psi (x)$$

is clearly an algebraic integer. As in the case of classical Gauss sums, the absolute value of $`G(f,\psi )`$ can easily be determined:

###### Lemma 1.

If $`f`$ satisfies 1.1, then

$$|G(f,\psi _t)|=\{\begin{array}{cc}\sqrt{p}\hfill & \text{if }t0modp\text{,}\hfill \\ 0\hfill & \text{if }t0modp\text{.}\hfill \end{array}$$

###### Proof.

We have

$$|G(f,\psi _t)|^2=\underset{x,yF_p}{}f(x)\overline{f(y)}\psi (t(xy))=\underset{x,hF_p}{}f(x)\overline{f(x+h)}\psi (th)$$

$$=\psi (0)\underset{xF_p}{}f(x)\overline{f(x)}+\underset{hF_p^\times }{}\psi (th)\underset{xF_p}{}f(x)\overline{f(x+h)}$$

$$=p1\underset{hF_p^\times }{}\psi (th)=\{\begin{array}{cc}p\hfill & \text{if }t0modp\text{,}\hfill \\ 0\hfill & \text{if }t0modp\text{,}\hfill \end{array}$$

The action of complex conjugation on $`K`$ is given by an element in $`G`$, and since $`G`$ is abelian, equation 1.1 is $`G`$-invariant. I.e., if $`f`$ satisfies 1.2, so does $`f^\sigma `$ for all $`\sigma G`$. But if $`\sigma G`$ then $`\sigma (G(f,\psi ))=G(f^\sigma ,\psi _t)`$, where $`\sigma (\zeta _p)=\zeta _p^t`$. Since $`f^\sigma `$ also satisfies 1.1, we find that $`|G(f^\sigma ,\psi _t)|=p^{1/2}`$, and hence the $`𝐐`$-norm of $`G(f,\psi )`$ is a power of $`p`$. The factorization of the principal ideal $`G(f,\psi )𝔒_L`$ thus consists only of prime ideals $`𝔓_L|p`$.

It is well known that $`𝐐(\zeta _{p^k})/𝐐`$ is totally ramified over $`p`$, and that $`𝐐(\zeta _n)/𝐐`$ does not ramify at $`p`$ if $`(n,p)=1`$. Comparing ramification indices gives that if $`𝔓_K`$ is a prime ideal in $`𝔒_K`$ that divides $`p`$, then $`𝔓_K`$ is totally ramified in $`L`$. In particular, if $`𝔓_L`$ is any prime ideal in the ring of integers in $`𝔒_L`$ that lies above $`p`$, then $`\sigma (𝔓_L)=𝔓_L`$ for all $`\sigma H`$.

Let $`l=\mathrm{max}(1,k)`$. Then $`H`$ consists of elements $`\sigma _t`$ such that

$`\sigma _t(\zeta _{p^l})=\zeta _{p^l}^t,`$ $`\sigma _t(\zeta _n)=\zeta _n.`$
Choose $`t`$ so that $`\sigma _t`$ generates $`H`$. Applying $`\sigma _t`$ to the principal ideal

$$G(f,\psi )𝔒_L=\underset{𝔓_L|p}{}𝔓_{L}^{}{}_{}{}^{\eta (𝔓_L)}$$

we find that

$$\sigma _t(G(f,\psi )𝔒_L)=\sigma _t(\underset{𝔓_L|p}{}𝔓_{L}^{}{}_{}{}^{\eta (𝔓_L)})=\underset{𝔓_L|p}{}𝔓_{L}^{}{}_{}{}^{\eta (𝔓_L)}=G(f,\psi )𝔒_L$$

and hence $`\sigma _t(G(f,\psi ))=uG(f,\psi )`$ for some unit $`u`$.

Since the absolute value of any complex embedding of $`G(f,\psi )`$ equals $`\sqrt{p}`$, we find that all conjugates of $`u=\sigma (G(f,\psi ))/G(f,\psi )`$ has absolute value one. Hence $`u`$ is in fact a root of unity, and there are integers $`a,b`$ such that

(2.1)
$$\sigma _t(G(f,\psi ))=\zeta _{p^l}^a\zeta _n^bG(f,\psi ).$$

## 3. The case $`(m,p)=1`$

Since $`f`$ is fixed by $`H`$ we find that $`\sigma _t(G(f,\psi ))=G(f,\psi _t)`$, and equation 2.1 can, after the change of variable $`xt^1x`$, be written as

(3.1)
$$\underset{x=1}{\overset{p1}{}}f(x)\psi (x)=\zeta _p^a\zeta _n^b\underset{x=1}{\overset{p1}{}}f(t^1x)\psi (x).$$

###### Lemma 2.

If $`f`$ takes values in $`n`$-th roots of unity for $`x0modp`$ and equation 3.1 holds then $`a0modp`$.

###### Proof.

From 3.1 we obtain that

(3.2)
$$\underset{i=1}{\overset{p1}{}}A_i\zeta _p^i=\underset{i=0}{\overset{p1}{}}B_i\zeta _p^i$$

where $`A_i=f(i)`$ and $`B_i=\zeta _n^bf(t^1(i+a))`$. (Note that $`B_{pa}=0`$.) Since $`1=_{i=1}^{p1}\zeta _p^i`$ we may rewrite 3.2 as

(3.3)
$$\underset{i=1}{\overset{p1}{}}A_i\zeta _p^i=\underset{i=1}{\overset{p1}{}}(B_iB_0)\zeta _p^i.$$

The elements $`\{\zeta _p,\zeta _p^2,\zeta _p^3,\mathrm{}\zeta _p^{p1}\}`$ are linearly independent over $`K`$, hence $`A_i=B_iB_0`$. From lemma 1 we have $`_{x=0}^{p1}f(x)=0`$, which implies that $`_{i=1}^{p1}A_i=0`$, as well as $`_{i=0}^{p1}B_i=0`$. Therefore,

$$0=\underset{i=1}{\overset{p1}{}}A_i=\underset{i=1}{\overset{p1}{}}(B_iB_0)=\underset{i=0}{\overset{p1}{}}B_ipB_0=pB_0.$$

But $`B_0=\zeta _n^bf(t^1(0+a))`$ which is nonzero unless $`a0modp`$. ∎

Thus

(3.4)
$$\underset{x=1}{\overset{p1}{}}f(x)\psi (x)=\zeta _n^b\underset{x=1}{\overset{p1}{}}f(t^1x)\psi (x)$$

and the linear independence of $`\{\zeta _p,\zeta _p^2,\zeta _p^3,\mathrm{}\zeta _p^{p1}\}`$ over $`K`$ implies that

$$f(t^1x)=f(x)\zeta _n^b$$

for all $`x0`$. Thus

$$f(t^k)=f(t^{(k1)})\zeta _n^b=\mathrm{}=f(1)\zeta _n^{kb}=\zeta _n^{kb}.$$

Taking $`k=p1`$ we find that $`\zeta _n^b`$ is a $`(p1)`$-th root of unity, and that $`f`$ is a multiplicative character.

## 4. The general case

In this case $`m=np^k`$ where $`(n,p)=1`$ and $`k>0`$. We will need the following:

###### Lemma 3.

If $`a_iK`$ and $`_{i=0}^{p^k1}a_i\zeta _{p^k}^iK(\zeta _p)`$ then

(4.1)
$$\underset{i=0}{\overset{p^k1}{}}a_i\zeta _{p^k}^i=\underset{j=0}{\overset{p1}{}}a_{p^{k1}j}\zeta _p^j$$

###### Proof.

We may assume that $`k>1`$. The minimal polynomial for $`\zeta _{p^k}`$ (over $`K`$ as well as over $`𝐐`$) is given by

$$\frac{x^{p^k}1}{x^{p^{k1}}1}=1+x^{p^{k1}}+x^{2p^{k1}}+\mathrm{}+x^{(p1)p^{k1}}.$$

Hence, by letting $`\stackrel{~}{i}[0,p^{k1}1]`$ be a representative of $`i`$ modulo $`p^{k1}`$, we can rewrite the left hand side of equation 4.1 as

$$\underset{i=0}{\overset{(p1)p^{k1}1}{}}(a_ia_{(p1)p^{k1}+\stackrel{~}{i}})\zeta _{p^k}^i$$

with no further relations among the $`\zeta _{p^k}^i`$’s, and thus

$$\underset{i=0}{\overset{(p1)p^{k1}1}{}}(a_ia_{(p1)p^{k1}+\stackrel{~}{i}})\zeta _{p^k}^iK(\zeta _p)$$

if and only if $`a_ia_{(p1)p^{k1}+\stackrel{~}{i}}=0`$ for all $`i`$ not congruent to zero modulo $`p^{k1}`$. ∎

Recall from equation 2.1 (note that $`l=k`$ since $`k1`$) that

$$\sigma _t(G(f,\psi ))=\zeta _{p^k}^a\zeta _n^bG(f,\psi ).$$

Let $`\stackrel{~}{G}=\zeta _{p^k}^sG(f,\psi )`$ where $`\sigma _t(\zeta _{p^k}^s)/\zeta _{p^k}^s=\zeta _{p^k}^a`$. (Such an $`s`$ exists as $`\sigma _t(\zeta _{p^k}^s)/\zeta _{p^k}^s=\zeta _{p^k}^{(t1)s}`$, and $`t1modp`$ since $`\sigma _t`$ generates $`H`$.) We then have

$$\sigma _t(\stackrel{~}{G})=\sigma _t(\zeta _{p^k}^sG(f,\psi ))$$

$$=\sigma _t(\zeta _{p^k}^s)\sigma _t(G(f,\psi ))=\sigma _t(\zeta _{p^k}^s)\zeta _{p^k}^a\zeta _n^bG(f,\psi ))=\zeta _n^b\stackrel{~}{G}.$$

The following lemma shows that $`\stackrel{~}{G}`$ must transform by a nontrivial $`n`$-th root of unity:

###### Lemma 4.

There is no integer $`s`$ such that $`\zeta _{p^k}^sG(f,\psi )K`$.

###### Proof.

We first assume that $`\zeta _{p^k}^s=1`$. Let $`G(f,\psi )𝔒_L=_{𝔓_L|p}𝔓_{L}^{}{}_{}{}^{\eta (𝔓_L)}`$ be the factorization of the principal ideal $`G(f,\psi )𝔒_L`$. Since $`p`$ does not ramify in $`K`$, we have $`p𝔒_K=_{𝔓_K|p}𝔓_K`$, and hence $`p𝔒_L=_{𝔓_L|p}𝔓_{L}^{}{}_{}{}^{e}`$ where $`e`$ is the ramification index of $`𝔓_K`$ in $`L`$.

Since $`\psi (x)=\zeta _p^x`$ is congruent to $`1`$ modulo $`𝔓_L`$ for all $`x`$, we find that

$$G(f,\psi )=\underset{x=0}{\overset{p1}{}}f(x)\psi (x)\underset{x=1}{\overset{p1}{}}f(x)mod𝔓_L$$

for all $`𝔓_L|p`$. Now, since $`f(0)=0`$, we have $`_{x=1}^{p1}f(x)=G(f,\psi _0)`$ and by lemma 1, $`G(f,\psi _0)=0`$. Thus $`G(f,\psi )𝔓_L`$ for all $`𝔓_L|p`$, i.e., $`\eta (𝔓_L)>0`$ for all $`𝔓_L|p`$. But if $`G(f,\psi )K`$ then $`e|\eta (𝔓_L)`$ for all $`𝔓_L|p`$, and since complex conjugation permutes the set of primes of $`𝔒_L`$ that lies above $`p`$, and

$$p=G(f,\psi )\overline{G(f,\psi )},$$

we get that $`𝔓_{L}^{}{}_{}{}^{2e}|p𝔒_L`$ for all $`𝔓_L`$, contradicting that the ramification index is $`e`$.

For the general case, the previous argument carries through by noting that $`\zeta _{p^k}^s`$ is a unit (and thus multiplication of $`G(f,\psi )`$ by $`\zeta _{p^k}^s`$ does not change the ideal factorization) and that $`G(f,\psi )𝔓_L`$ if and only if $`\zeta _{p^k}G(f,\psi )𝔓_L`$. ∎

Since $`\sigma _t`$ has order $`p^{k1}(p1)`$ and $`(n,p)=1`$ we find that $`\zeta _n^b`$ must be a nontrivial $`(p1)`$-th root of unity. Hence there exists a nontrivial multiplicative character $`\chi `$ of $`F_p^\times `$ such that $`\chi (t^1)=\zeta _n^b`$. But $`\sigma _t(G(\chi ,\psi ))=G(\chi ,\psi _t)=\chi (t^1)G(\chi ,\psi )`$ and thus

$$\delta =\frac{\stackrel{~}{G}}{G(\chi ,\psi )}$$

is $`\sigma _t`$-invariant and hence an element of $`K`$. Moreover, $`|\delta |=1`$ (for all complex embeddings) since $`|\stackrel{~}{G}|=|G(\chi ,\psi )|=p^{1/2}`$.

Write $`f(x)=f_1(x)f_2(x)`$ where $`f_1(x)`$ takes values in $`p^k`$-th roots of unity and $`f_2(x)`$ takes values in $`n`$-th roots of unity. We will show that $`f_1(x)`$ must be constant.

###### Lemma 5.

Let

$$a_i=\underset{x:\zeta _{p^k}^sf_1(x)\psi (x)=\zeta _{p^k}^i}{}f_2(x)$$

If

(4.2)
$$\zeta _{p^k}^s\underset{x=1}{\overset{p1}{}}f(x)\psi (x)=\delta \underset{x=1}{\overset{p1}{}}\chi (x)\psi (x),$$

then $`|a_i|=0`$ unless $`i=p^{k1}j`$ for $`j=1,2,\mathrm{},p1`$, in which case $`|a_i|=1`$. In particular, $`\zeta _{p^k}^sf_1(x)\psi (x)`$ ranges over all nontrivial $`p`$-th roots of unity.

###### Proof.

Collecting terms in 4.2 according to the values of $`\zeta _{p^k}^sf_1(x)\psi (x)`$, we obtain

(4.3)
$$\underset{i=0}{\overset{p^k1}{}}a_i\zeta _{p^k}^i=\delta \underset{i=1}{\overset{p1}{}}\chi (i)\zeta _p^iK(\zeta _p).$$

Clearly $`a_iK`$ and $`a_i0`$ for at most $`p1`$ values of $`i`$. Letting $`A_i=a_{p^{k1}i}`$ we may, by lemma 3, write equation 4.3 as

$$\underset{i=0}{\overset{p1}{}}A_i\zeta _p^i=\delta \underset{i=1}{\overset{p1}{}}\chi (i)\zeta _p^i.$$

Since $`1=_{i=1}^{p1}\zeta _p^i`$ we get that

$$\underset{i=1}{\overset{p1}{}}(A_iA_0)\zeta _p^i=\underset{i=0}{\overset{p1}{}}A_i\zeta _p^i=\delta \underset{i=1}{\overset{p1}{}}\chi (i)\zeta _p^i$$

and hence $`A_iA_0=\delta \chi (i)`$ for all $`i`$.

Since $`a_i0`$ for at most $`p1`$ values of $`i`$, $`A_00`$ implies that $`A_j=0`$ for some $`j0`$, and thus $`|A_0|=|\delta \chi (j)A_j|=1`$. Since

$$0=\delta \underset{i=1}{\overset{p1}{}}\chi (i)=\underset{i=1}{\overset{p1}{}}(A_iA_0)=\underset{i=0}{\overset{p1}{}}A_ipA_0,$$

we find that $`|_{i=0}^{p1}A_i|=p|A_0|=p`$. On the other hand, $`|_{i=0}^{p1}A_i|_{x=1}^{p1}|f_2(x)|=p1`$. Thus $`A_0=0`$, and it follows that $`A_i=\delta \chi (i)`$ for $`i0`$. In other words, $`a_{p^{k1}j}=A_j=\delta \chi (j)`$ for $`j=1,2,\mathrm{},p1`$, and since there are at most $`p1`$ nonzero values among the $`a_i`$’s, the remaining ones must all be equal to zero. ∎

Now, the lemma gives that $`\zeta _{p^k}^sf_1(1)\psi (1)=\zeta _{p^k}^s\zeta _p`$ is a $`p`$-th root of unity, hence $`p^{k1}`$ must divide $`s`$, and the nonzero values of $`f_1(x)\psi (x)`$ are thus distinct $`p`$-th roots of unity. Replacing $`\psi `$ by $`\psi _r`$, for $`r0modp`$, in the previous argument gives that $`f_1(x)\psi (rx)`$ also ranges over distinct $`p`$-th roots of unity. On the other hand, if $`f_1(x)`$ is not constant, then there exists $`r0modp`$ such that the set $`\{f_1(x)\psi _r(x)\}_{x=1}^{p1}`$ contains strictly less than $`p1`$ elements. (If $`f_1(x_1)f_1(x_2)`$, write $`f_1(x_1)=\zeta _p^{y_1},f_1(x_2)=\zeta _p^{y_2}`$ and take $`r(y_2y_1)(x_2x_1)^1modp`$.) Hence $`f_1(x)`$ must be constant, and since $`f_1(1)=1`$, we find that the nonzero values of $`f(x)`$ are in fact $`n`$-th roots of unity. The result has thus been reduced to the case $`(m,p)=1`$.