Primorial: asymptotics and bounds

  • Context: Graduate 
  • Thread starter Thread starter CRGreathouse
  • Start date Start date
  • Tags Tags
    Asymptotics Bounds
Click For Summary

Discussion Overview

The discussion revolves around approximations and bounds for the primorial function, defined as the product of all prime numbers less than or equal to a given number n. Participants explore asymptotic behaviors, comparisons to known functions, and the implications of various mathematical results related to the primorial function.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant notes that the primorial function n# is asymptotically similar to e^n, while also providing bounds such as π(n)! ≤ n# ≤ n! and questioning if e^n is consistently larger than n# for large n.
  • Another participant references Chris Caldwell's work, suggesting that the asymptotic expression n# ∼ e^n is related to the Prime Number Theorem (PNT) and introduces an asymptotic expression from Sloane's OEIS that appears to be larger than e^n.
  • A different participant mentions Pierre Dusart's work, which provides bounds for the primorial function and indicates that θ(x) < x for x < 10^11.
  • Concerns are raised about the equivalence of n# ∼ e^n and the implications for the logarithmic behavior of n#, with one participant asserting that this equivalence is false and that the difference between log(n#) and n grows arbitrarily large.
  • Another participant expresses skepticism about the asymptotic from Sloane's OEIS, suggesting it may not hold true and provides an upper bound for x# based on Dusart's results.
  • Participants acknowledge confusion between different notations and express a shared experience of making similar mistakes in their calculations.

Areas of Agreement / Disagreement

Participants express differing views on the validity of certain asymptotic expressions and their implications. There is no consensus on the equivalence of n# ∼ e^n, and multiple competing interpretations of the bounds and approximations for the primorial function are present.

Contextual Notes

Some participants highlight limitations in their understanding and calculations, particularly regarding the relationship between different mathematical notations and the implications of various bounds. The discussion reflects ongoing uncertainty and refinement of ideas.

CRGreathouse
Science Advisor
Homework Helper
Messages
2,832
Reaction score
0
I was wondering if anyone knew of good approximations of the primorial function for large numbers, or of reasonable bounds for it. By primorial I mean:

[tex]n\#=\prod_{p\le n}p[/tex] for p prime.

All I know is [tex]n\#\sim e^n[/tex] and the trivial [tex]\pi(n)!\le n\#\le n![/tex].

For small numbers (n < 1,000,000,000), e^n is larger than n#. Does this hold generally?

Obviously, the crude upper limit can be hacked slightly: take out the evens from the factorial and multiply by 2, take out the multiples of 3 and divide by 3, etc. In the extreme case this is just calculating the primorial itself, so not very useful. The lower bound can be improved to

[tex]p_n\#\le n!\prod^n_{k=5}\ln k+\ln\ln k-1[/tex]

using lower bounds on the primes, but this is unwieldy, and I'm pretty sure [tex]\lim_{n\rightarrow\infty}\frac{p_n\#}{n!\Pi}=0[/tex] where the product is as above. Using the best bounds out there (by Dusart) improves the result slightly at the cost of complicating it further.
 
Last edited:
Physics news on Phys.org
Well, I did some 'research' online:

Chris Caldwell's page says that [tex]n\#\sim e^n[/tex] is equivilent to the PNT. I knew the formula already, but at least this gives me a different way of thinking about it.

I found the asymptotic expression [tex]\exp((1 + o(1))\cdot n\cdot\log(n))=n^n(1+o(1))^{n\log(n)}[/tex] on Sloane's OEIS. This is larger than e^n, which already appears to be an upper bound... any ideas about this?
 
Last edited:
I forgot that Pierre Dusart worked on this already. He shows that [tex]\vartheta(x)<x[/tex] for x < 10^{11} (100 billion, better than I was able to calculate by a factor of 50) and gives some explicit bounds.
 
CRGreathouse said:
Chris Caldwell's page says that [tex]n\#\sim e^n[/tex] is equivilent to the PNT. I knew the formula already, but at least this gives me a different way of thinking about it.

This isn't true though. if n#~e^n then n#/e^n->1, which means log n#-n->0, i.e. [tex]\vartheta(n)-n\rightarrow 0[/tex]. This is false, the difference gets arbitrarily large.

What is true is log(n#)~n is equivalent to the prime number theorem.

That asymptotic from Sloanes appears to be for a(n), the nth primorial, a(1)=1,a(2)=2,a(3)=6,a(4)=30 etc. I don't believe this asymptotic for the same reasons, unless I misunderstand what they mean.
 
shmoe said:
This isn't true though. if n#~e^n then n#/e^n->1, which means log n#-n->0, i.e. [tex]\vartheta(n)-n\rightarrow 0[/tex]. This is false, the difference gets arbitrarily large.

:blushing: Sorry, I've been converting between the two much too readily. This is the second mistake I've made between theta and # in as many days (the last was writing one for the other!).

shmoe said:
That asymptotic from Sloanes appears to be for a(n), the nth primorial, a(1)=1,a(2)=2,a(3)=6,a(4)=30 etc. I don't believe this asymptotic for the same reasons, unless I misunderstand what they mean.

I don't trust it. I can get the upper bound

[tex]x\#<x^{\frac{x}{\ln x-1.1}}[/tex] for x >= 4

from Dusart's bounds on primes and the trivial observation that [tex]\vartheta(x)\le\pi(x)\ln x[/tex].
 
CRGreathouse said:
:blushing: Sorry, I've been converting between the two much too readily. This is the second mistake I've made between theta and # in as many days (the last was writing one for the other!).

I'd be lying if I said I've never made this mistake myself.


CRGreathouse said:
I don't trust it. I can get the upper bound

[tex]x\#<x^{\frac{x}{\ln x-1.1}}[/tex] for x >= 4

from Dusart's bounds on primes and the trivial observation that [tex]\vartheta(x)\le\pi(x)\ln x[/tex].

It looks fine if you take logs, log a(n)~n*log(n). Remember a(n)=p(n)#, where p(n) is the nth prime and p(n)~nlog(n).
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 33 ·
2
Replies
33
Views
9K
  • · Replies 150 ·
6
Replies
150
Views
21K
  • · Replies 125 ·
5
Replies
125
Views
20K
  • · Replies 4 ·
Replies
4
Views
3K