

A147815


An integer differential Schroedinger Hamiltonian for primes as a triangular sequence: g(n)]= K=integer<=n: Ep(n,K) = (K n Prime[n]  2 K n Prime[1 + n] + K n Prime[2 + n] + Prime[n] g(n))/(n*Prime[n]).


0



2, 3, 3, 3, 3, 5, 11, 17, 23, 7, 1, 9, 17, 25, 11, 21, 31, 41, 51, 61, 13, 1, 11, 23, 35, 47, 59, 17, 31, 45, 59, 73, 87, 101, 115, 19, 35, 51, 67, 83, 99, 115, 131, 147, 23, 13, 49, 85, 121, 157, 193, 229, 265, 301, 29, 69, 109, 149, 189, 229, 269
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


COMMENTS

The equation derivation goes like this: If you have a function of integers n : f[n] then the derivative equivalents are: f[n] f[n+1]f[n] f[n+2]2*f[n+1]+f[n] f[n+3]3*f[n+2]+2*f[n+1]f[n] Something like: d[n_,k_]=Sum[(1)^m*Binomial[n,m]*f[n+m],{m,0,k}] So the second order differential equation equivalent would look like: a*(f[n+2]2*f[n+1]+f[n])+b*(f[n+1]f[n])+c*f[n]=0.
I ran across this in an old problems book by Ulam where they were using the "new" computers ( late 50's / early 60's) to simulate different differential equations. So quantum mechanic is based on Schroedinger's equation: Hamiltonian : Kinetic energy term: K =( hbar2/(2*m))* d2/dx2. Potential energy term: V=g[x]/x. The total Hamiltonian energy ( a constant with quantum levels of n): E K*phi(x)+V*phi(x)=E*phi(x) Putting that into f[n] terms: ( hbar2/(2*m))*(f[n+2]2*f[n+1]+f[n])+g[n]*f[n]/n=E[n]*f[n].
Making ( hbar2/(2*m)) a constant K Clear[f, E0, K, n, g] Solve[K*(f[n + 2]  2*f[n + 1] + f[n]) + g[n ]*f[n]/n  E0*f[n] == 0, E0] E0=(K n f[n]  2 K n f[1 + n] + K n f[2 + n] + f[n] g[n])/(n*f[n]) So a quantum prime second order differential Hamiltonian might be: Ep[n_] = (K n Prime[n]  2 K n Prime[1 + n] + K n Prime[2 + n] + Prime[n] g[n])/(n*Prime[n]).


LINKS

Table of n, a(n) for n=1..61.


FORMULA

g(n)]= K=integer<=n: Ep(n,K) = (K n Prime[n]  2 K n Prime[1 + n] + K n Prime[2 + n] + Prime[n] \ g(n))/(n*Prime[n]); t(n,m)=n*Prime[n]*Ep(n,m).


EXAMPLE

{2, 3},
{3, 3, 3},
{5, 11, 17, 23},
{7, 1, 9, 17, 25},
{11, 21, 31, 41, 51, 61},
{13, 1, 11, 23, 35, 47, 59},
{17, 31, 45, 59, 73, 87, 101, 115},
{19, 35, 51, 67, 83, 99, 115,131, 147},
{23, 13, 49, 85, 121, 157, 193, 229, 265, 301},
{29, 69, 109, 149, 189, 229, 269, 309, 349, 389, 429},
...


MATHEMATICA

Clear[K, g, Ep]; g[n_] = 1; Ep[n_, K_] = (K n Prime[n]  2 K n Prime[1 + n] + K n Prime[2 + n] + Prime[n] g[n])/(n*Prime[n]); Table[n*Prime[n]*Ep[n, m], {n, 1, 10}, {m, 0, n}]; Flatten[%]


CROSSREFS

Sequence in context: A129263 A035367 A042959 * A227246 A200924 A111913
Adjacent sequences: A147812 A147813 A147814 * A147816 A147817 A147818


KEYWORD

sign,tabf,less


AUTHOR

Roger L. Bagula, Nov 13 2008


STATUS

approved



