Introduction
Here is an interesting, but slightly detached from the previous two articles, result ...
Let's look at the following two sequences:
an=√pn+1−√pnbn=ln(1+√pn+11+√pn)1+√pn+1
Here is a short Python code to visualise the sequences:
import math
import matplotlib.pyplot as plt
primes = []
def isPrime(n):
l = int(math.sqrt(n)) + 1
for i in xrange(2,l):
if (n % i) == 0:
return False
return True
def calculateLog(sqrt_p1, sqrt_p2):
sqrt_1_p2 = 1.0 + sqrt_p2
r = math.log(sqrt_1_p2/(1.0 + sqrt_p1))
return r * sqrt_1_p2
N = 1500000
print "populate primes ..."
for i in xrange(2, N):
if isPrime(i):
primes.append(i);
sqrt_diff = []
diff = []
log_calcs = []
x = []
for i in xrange(1, len(primes)):
sqrt_p2 = math.sqrt(primes[i])
sqrt_p1 = math.sqrt(primes[i-1])
sqrt_diff.append(sqrt_p2 - sqrt_p1)
diff.append(primes[i] - primes[i-1])
log_calcs.append(calculateLog(sqrt_p1, sqrt_p2))
x.append(i)
for i in xrange(len(sqrt_diff)):
print sqrt_diff[i]," = s(",primes[i+1],") - s(",primes[i],")"
plt.subplot(311)
plt.plot(x, sqrt_diff)
plt.subplot(312)
plt.plot(x, log_calcs)
plt.subplot(313)
plt.hist(diff, 1000)
plt.show()
And here is how both sequences look like (an the first and bn the second): 
Quite asymptotic, aren't they? Indeed they are ...
Lemma 3.
√pn+1−√pn≤ln(1+√pn+11+√pn)1+√pn+1≤(1+√pn+11+√pn)⋅(√pn+1−√pn)
Let's look at this function f6(x)=√pn+11+x⋅√pn+1. Obviously ln(1+x⋅√pn+1)′=f6(x) As a result
∫1√pnpn+1f6(x)dx=ln(1+x⋅√pn+1)|1√pnpn+1=ln(1+√pn+11+√pn)
According to Mean Value Theorem, ∃μ∈(√pnpn+1,1) such that:
∫1√pnpn+1f6(x)dx=f6(μ)⋅(1−√pnpn+1)
Putting all together:
ln(1+√pn+11+√pn)=√pn+1−√pn1+μ⋅√pn+1
Because:
√pnpn+1<μ<1⇒1+√pn<1+μ⋅√pn+1<1+√pn+1
And we get:
√pn+1−√pn1+√pn+1≤ln(1+√pn+11+√pn)≤√pn+1−√pn1+√pn
which proves this lemma.
Noting Δn=√pn+1−√pn, this becomes:
Δn1+√pn+1≤ln(1+Δn1+√pn)≤Δn1+√pn
or
1+√pn1+√pn+1≤ln(1+Δn1+√pn)1+√pnΔn≤1
Is this result of any use? I don't know yet, but it looks like:
(1+Δn1+√pn)1+√pnΔn→e,n→∞