dif1d: provides analytical solution to 1d diffusion problem in a cylindrical geometry with inner sink.

This program computes the flux in one-dimensional diffusion problem in a cylindrical geometry, in application to studying the loss-cone problem of consumption of stars by a supermassive black hole. The evolution of the distribution function in the interval is governed by the equation:

with initial conditions , subject to the boundary conditions

The solution for the case is given by [Milosavljevic&Merritt 2003] in terms of series of Bessel functions (MM03, Eqs.24-28). The generalization to arbitrary is done by introducing the effective value of at which the distribution function would be zero. The general form of the boundary condition interpolates between empty ( ) and full ( ) regimes.

This method is also used to estimate the capture rates in axisymmetric galactic nuclei, as described by an effective one-dimensional approximation of [Vasiliev&Merritt 2013]. In this case, the meaning of is the effective capture boundary , defined in terms of flattening coefficient (Eq.72 of VM13).

In addition, the draining of loss region is accounted for with an approximate method described in section 6 of [Vasiliev&Merritt 2013]: initially, the draining flux is equal to the full loss cone capture rate, and declines with time so that the total captured mass from the draining region approaches Rdrain. At the same time, the flux due to relaxation is gradually scaled up from 0 to the flux in the time-dependent solution with initially empty loss region.

The methods implemented in this code are based on the following papers:

Milosavljevic, M., and Merritt, D., Long-Term Evolution of Massive Black Hole Binaries, Astrophysical Journal, 596, 860-878, 2003

Merritt, D. and Wang, J., Loss Cone Refilling Rates in Galactic Nuclei, Astrophysical Journal, 621, L101-L104, 2005

Vasiliev, E. and Merritt, D., The loss-cone problem in axisymmetric nuclei, Astrophysical Journal, 774, 87, 2013

Vasiliev, E. Rates of capture of stars by supermassive black holes in non-spherical galactic nuclei, Class.Quant.Grav., 2014