Computational science and re-discovery: open-source implementations of ellipsoidal harmonics for problems in potential theory

We present two open-source (BSD) implementations of ellipsoidal harmonic expansions for solving problems of potential theory using separation of variables. Ellipsoidal harmonics are used surprisingly infrequently, considering their substantial value for problems ranging in scale from molecules to the entire solar system. In this article, we suggest two possible reasons for the paucity relative to spherical harmonics. The first is essentially historical---ellipsoidal harmonics developed during the late 19th century and early 20th, when it was found that only the lowest-order harmonics are expressible in closed form. Each higher-order term requires the solution of an eigenvalue problem, and tedious manual computation seems to have discouraged applications and theoretical studies. The second explanation is practical: even with modern computers and accurate eigenvalue algorithms, expansions in ellipsoidal harmonics are significantly more challenging to compute than those in Cartesian or spherical coordinates. The present implementations reduce the "barrier to entry" by providing an easy and free way for the community to begin using ellipsoidal harmonics in actual research. We demonstrate our implementation using the specific and physiologically crucial problem of how charged proteins interact with their environment, and ask: what other analytical tools await re-discovery in an era of inexpensive computation?Comment: 25 pages, 3 figure

Stationary Regime of Random Resistor Networks Under Biased Percolation

The state of a 2-D random resistor network, resulting from the simultaneous evolutions of two competing biased percolations, is studied in a wide range of bias values. Monte Carlo simulations show that when the external current $I$ is below the threshold value for electrical breakdown, the network reaches a steady state with a nonlinear current-voltage characteristic. The properties of this nonlinear regime are investigated as a function of different model parameters. A scaling relation is found between $/_0$ and $I/I_0$, where $$ is the average resistance, $_0$ the linear regime resistance and $I_0$ the threshold value for the onset of nonlinearity. The scaling exponent is found to be independent of the model parameters. A similar scaling behavior is also found for the relative variance of resistance fluctuations. These results compare well with resistance measurements in composite materials performed in the Joule regime up to breakdown.Comment: 9 pages, revtex, proceedings of the Merida Satellite Conference STATPHYS2