The Riemann zeta function and tuning: Difference between revisions
Wikispaces>genewardsmith **Imported revision 353276654 - Original comment: ** |
Wikispaces>genewardsmith **Imported revision 353296842 - Original comment: ** |
||
| Line 1: | Line 1: | ||
<h2>IMPORTED REVISION FROM WIKISPACES</h2> | <h2>IMPORTED REVISION FROM WIKISPACES</h2> | ||
This is an imported revision from Wikispaces. The revision metadata is included below for reference:<br> | This is an imported revision from Wikispaces. The revision metadata is included below for reference:<br> | ||
: This revision was by author [[User:genewardsmith|genewardsmith]] and made on <tt>2012-07- | : This revision was by author [[User:genewardsmith|genewardsmith]] and made on <tt>2012-07-16 01:32:48 UTC</tt>.<br> | ||
: The original revision id was <tt> | : The original revision id was <tt>353296842</tt>.<br> | ||
: The revision comment was: <tt></tt><br> | : The revision comment was: <tt></tt><br> | ||
The revision contents are below, presented both in the original Wikispaces Wikitext format, and in HTML exactly as Wikispaces rendered it.<br> | The revision contents are below, presented both in the original Wikispaces Wikitext format, and in HTML exactly as Wikispaces rendered it.<br> | ||
| Line 137: | Line 137: | ||
Removing 2 leads to increasing adjusted peak values corresponding to the division of 3 (the "tritave") into 4, 7, 9, 13, 15, 17, 26, 32, 39, 45, 52, 56, 71, 75, 88, 131, 245, 316 ... parts. A striking feature of this list is the appearance not only of [[13edt]], the [[Bohlen-Pierce]] division of the tritave, but the multiples 26, 39 and 52 also. | Removing 2 leads to increasing adjusted peak values corresponding to the division of 3 (the "tritave") into 4, 7, 9, 13, 15, 17, 26, 32, 39, 45, 52, 56, 71, 75, 88, 131, 245, 316 ... parts. A striking feature of this list is the appearance not only of [[13edt]], the [[Bohlen-Pierce]] division of the tritave, but the multiples 26, 39 and 52 also. | ||
=The Black Magic Formulas= | |||
When [[Gene Ward Smith|Gene Smith]] discovered these formulas in the 70s, he thought of them as "black magic" formulas not because of any aura of evil, but because they seemed mysteriously to give you something for next to nothing. They are based on Gram points and the Riemann-Siegel theta function θ(t). Recall that a Gram point is a points on the critical line where ζ(1/2 + ig) is real. This implies that exp(iθ(g)) is real, so that θ(g)/π is an integer. Theta has an [[http://en.wikipedia.org/wiki/Asymptotic_expansion|asymptotic expansion]] | |||
[[math]] | |||
\theta(t) \sim \frac{t}{2}\log \frac{t}{2\pi} - \frac{t}{2} - \frac{\pi}{8}+\frac{1}{48t}+ \frac{7}{5760t^3}+\cdots | |||
[[math]] | |||
=Computing zeta= | =Computing zeta= | ||
| Line 154: | Line 161: | ||
[[http://terrytao.wordpress.com/2009/07/12/selbergs-limit-theorem-for-the-riemann-zeta-function-on-the-critical-line/|Selberg's limit theorem]] by Terence Tao [[http://www.webcitation.org/5xrvgjW6T|Permalink]]</pre></div> | [[http://terrytao.wordpress.com/2009/07/12/selbergs-limit-theorem-for-the-riemann-zeta-function-on-the-critical-line/|Selberg's limit theorem]] by Terence Tao [[http://www.webcitation.org/5xrvgjW6T|Permalink]]</pre></div> | ||
<h4>Original HTML content:</h4> | <h4>Original HTML content:</h4> | ||
<div style="width:100%; max-height:400pt; overflow:auto; background-color:#f8f9fa; border: 1px solid #eaecf0; padding:0em"><pre style="margin:0px;border:none;background:none;word-wrap:break-word;width:200%;white-space: pre-wrap ! important" class="old-revision-html"><html><head><title>The Riemann Zeta Function and Tuning</title></head><body><!-- ws:start:WikiTextTocRule: | <div style="width:100%; max-height:400pt; overflow:auto; background-color:#f8f9fa; border: 1px solid #eaecf0; padding:0em"><pre style="margin:0px;border:none;background:none;word-wrap:break-word;width:200%;white-space: pre-wrap ! important" class="old-revision-html"><html><head><title>The Riemann Zeta Function and Tuning</title></head><body><!-- ws:start:WikiTextTocRule:31:&lt;img id=&quot;wikitext@@toc@@flat&quot; class=&quot;WikiMedia WikiMediaTocFlat&quot; title=&quot;Table of Contents&quot; src=&quot;/site/embedthumbnail/toc/flat?w=100&amp;h=16&quot;/&gt; --><!-- ws:end:WikiTextTocRule:31 --><!-- ws:start:WikiTextTocRule:32: --><a href="#Preliminaries">Preliminaries</a><!-- ws:end:WikiTextTocRule:32 --><!-- ws:start:WikiTextTocRule:33: --> | <a href="#Into the critical strip">Into the critical strip</a><!-- ws:end:WikiTextTocRule:33 --><!-- ws:start:WikiTextTocRule:34: --> | <a href="#The Z function">The Z function</a><!-- ws:end:WikiTextTocRule:34 --><!-- ws:start:WikiTextTocRule:35: --> | <a href="#Zeta EDO lists">Zeta EDO lists</a><!-- ws:end:WikiTextTocRule:35 --><!-- ws:start:WikiTextTocRule:36: --> | <a href="#Removing primes">Removing primes</a><!-- ws:end:WikiTextTocRule:36 --><!-- ws:start:WikiTextTocRule:37: --> | <a href="#The Black Magic Formulas">The Black Magic Formulas</a><!-- ws:end:WikiTextTocRule:37 --><!-- ws:start:WikiTextTocRule:38: --> | <a href="#Computing zeta">Computing zeta</a><!-- ws:end:WikiTextTocRule:38 --><!-- ws:start:WikiTextTocRule:39: --> | <a href="#Links">Links</a><!-- ws:end:WikiTextTocRule:39 --><!-- ws:start:WikiTextTocRule:40: --> | ||
<!-- ws:end:WikiTextTocRule: | <!-- ws:end:WikiTextTocRule:40 --><!-- ws:start:WikiTextHeadingRule:15:&lt;h1&gt; --><h1 id="toc0"><a name="Preliminaries"></a><!-- ws:end:WikiTextHeadingRule:15 -->Preliminaries</h1> | ||
Suppose x is a variable representing some equal division of the octave. For example, if x = 80, x reflects 80edo with a step size of 15 cents and with pure octaves. Suppose that x can also be continuous, so that it can also represent fractional or &quot;nonoctave&quot; divisions as well. The Bohlen-Pierce scale, 13 equal divisions of 3/1, is approximately 8.202 equal divisions of the &quot;octave&quot; (although the octave itself does not appear in this tuning), and would hence be represented by a value of x = 8.202.<br /> | Suppose x is a variable representing some equal division of the octave. For example, if x = 80, x reflects 80edo with a step size of 15 cents and with pure octaves. Suppose that x can also be continuous, so that it can also represent fractional or &quot;nonoctave&quot; divisions as well. The Bohlen-Pierce scale, 13 equal divisions of 3/1, is approximately 8.202 equal divisions of the &quot;octave&quot; (although the octave itself does not appear in this tuning), and would hence be represented by a value of x = 8.202.<br /> | ||
<br /> | <br /> | ||
| Line 213: | Line 220: | ||
so that we see that the absolute value of the zeta function serves to measure the error of an equal division.<br /> | so that we see that the absolute value of the zeta function serves to measure the error of an equal division.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:17:&lt;h1&gt; --><h1 id="toc1"><a name="Into the critical strip"></a><!-- ws:end:WikiTextHeadingRule:17 -->Into the critical strip</h1> | ||
So long as s is greater than or equal to one, the absolute value of the zeta function can be seen as an error measurement. However, the rationale for that view of things departs when s is less than one, particularly in the <a class="wiki_link_ext" href="http://mathworld.wolfram.com/CriticalStrip.html" rel="nofollow">critical strip</a>, when s lies between zero and one. As s approaches the value s=1/2 of the <a class="wiki_link_ext" href="http://mathworld.wolfram.com/CriticalLine.html" rel="nofollow">critical line</a>, the information content, so to speak, of the zeta function concerning higher primes increases and it behaves increasingly like a badness measure (or more correctly, since we have inverted it, like a goodness measure.) The quasi-symmetric <a class="wiki_link_ext" href="http://planetmath.org/encyclopedia/FunctionalEquationOfTheRiemannZetaFunction.html" rel="nofollow">functional equation</a> of the zeta function tells us that past the critical line the information content starts to decrease again, with 1-s and s having the same information content. Hence it is the zeta function between s=1/2 and s=1, and especially the zeta function along the critical line s=1/2, which is of the most interest.<br /> | So long as s is greater than or equal to one, the absolute value of the zeta function can be seen as an error measurement. However, the rationale for that view of things departs when s is less than one, particularly in the <a class="wiki_link_ext" href="http://mathworld.wolfram.com/CriticalStrip.html" rel="nofollow">critical strip</a>, when s lies between zero and one. As s approaches the value s=1/2 of the <a class="wiki_link_ext" href="http://mathworld.wolfram.com/CriticalLine.html" rel="nofollow">critical line</a>, the information content, so to speak, of the zeta function concerning higher primes increases and it behaves increasingly like a badness measure (or more correctly, since we have inverted it, like a goodness measure.) The quasi-symmetric <a class="wiki_link_ext" href="http://planetmath.org/encyclopedia/FunctionalEquationOfTheRiemannZetaFunction.html" rel="nofollow">functional equation</a> of the zeta function tells us that past the critical line the information content starts to decrease again, with 1-s and s having the same information content. Hence it is the zeta function between s=1/2 and s=1, and especially the zeta function along the critical line s=1/2, which is of the most interest.<br /> | ||
<br /> | <br /> | ||
| Line 220: | Line 227: | ||
Because the value of zeta increased continuously as it made its way from +infinity to the critical line, we might expect the values of zeta at these special Gram points to be relatively large. This would be especially true if -ζ'(z) is getting a boost from other small primes as it travels toward the Gram point. A complex formula due to <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann" rel="nofollow">Bernhard Riemann</a> which he failed to publish because it was so nasty becomes a bit simpler when used at a Gram point. It is named the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann%E2%80%93Siegel_formula" rel="nofollow">Riemann-Siegel formula</a> since <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Carl_Ludwig_Siegel" rel="nofollow">Carl Ludwig Siegel</a> went looking for it and was able to reconstruct it after rooting industriously around in Riemann's unpublished papers. From this formula, it is apparent that when x corresponds to a good edo, the value of ζ(1/2 + i g) at the corresponding Gram point should be especially large.<br /> | Because the value of zeta increased continuously as it made its way from +infinity to the critical line, we might expect the values of zeta at these special Gram points to be relatively large. This would be especially true if -ζ'(z) is getting a boost from other small primes as it travels toward the Gram point. A complex formula due to <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann" rel="nofollow">Bernhard Riemann</a> which he failed to publish because it was so nasty becomes a bit simpler when used at a Gram point. It is named the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann%E2%80%93Siegel_formula" rel="nofollow">Riemann-Siegel formula</a> since <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Carl_Ludwig_Siegel" rel="nofollow">Carl Ludwig Siegel</a> went looking for it and was able to reconstruct it after rooting industriously around in Riemann's unpublished papers. From this formula, it is apparent that when x corresponds to a good edo, the value of ζ(1/2 + i g) at the corresponding Gram point should be especially large.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:19:&lt;h1&gt; --><h1 id="toc2"><a name="The Z function"></a><!-- ws:end:WikiTextHeadingRule:19 -->The Z function</h1> | ||
The absolute value ζ(1/2 + i g) at a Gram point corresponding to an edo is near to a local maximum, but not actually at one. At the local maximum, of course, the partial derivative of ζ(1/2 + i t) with respect to t will be zero; however this does not mean its derivative there will be zero. In fact, the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann_hypothesis" rel="nofollow">Riemann hypothesis</a> is equivalent to the claim that all zeros of ζ'(s + i t) occur when s &gt; 1/2, which is where all known zeros lie. These do not have values of t corresponding to good edos. For this and other reasons, it is helpful to have a function which is real for values on the critical line but whose absolute value is the same as that of zeta. This is provided by the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Z_function" rel="nofollow">Z function</a>.<br /> | The absolute value ζ(1/2 + i g) at a Gram point corresponding to an edo is near to a local maximum, but not actually at one. At the local maximum, of course, the partial derivative of ζ(1/2 + i t) with respect to t will be zero; however this does not mean its derivative there will be zero. In fact, the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Riemann_hypothesis" rel="nofollow">Riemann hypothesis</a> is equivalent to the claim that all zeros of ζ'(s + i t) occur when s &gt; 1/2, which is where all known zeros lie. These do not have values of t corresponding to good edos. For this and other reasons, it is helpful to have a function which is real for values on the critical line but whose absolute value is the same as that of zeta. This is provided by the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Z_function" rel="nofollow">Z function</a>.<br /> | ||
<br /> | <br /> | ||
| Line 263: | Line 270: | ||
If you have access to <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Mathematica" rel="nofollow">Mathematica</a>, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematicia-generated plot of Z(2πx/ln(2)) in the region around 12edo:<br /> | If you have access to <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Mathematica" rel="nofollow">Mathematica</a>, which has Z, zeta and theta as a part of its suite of initially defined functions, you can do even better. Below is a Mathematicia-generated plot of Z(2πx/ln(2)) in the region around 12edo:<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextLocalImageRule: | <!-- ws:start:WikiTextLocalImageRule:41:&lt;img src=&quot;/file/view/plot12.png/219376858/plot12.png&quot; alt=&quot;&quot; title=&quot;&quot; /&gt; --><img src="/file/view/plot12.png/219376858/plot12.png" alt="plot12.png" title="plot12.png" /><!-- ws:end:WikiTextLocalImageRule:41 --><br /> | ||
<br /> | <br /> | ||
The peak around 12 is both higher and wider than the local maximums above 11 and 13, indicating its superiority as an edo. Note also that the peak occurs at a point slightly larger than 12; this indicates the octave is slightly compressed in the zeta tuning for 12. The size of a step in octaves is 1/x, and hence the size of the octave in the zeta peak value tuning for Nedo is N/x; if x is slightly larger than N as here with N=12, the size of the zeta tuned octave will be slightly less than a pure octave. Similarly, when the peak occurs with x less than N, we have stretched octaves.<br /> | The peak around 12 is both higher and wider than the local maximums above 11 and 13, indicating its superiority as an edo. Note also that the peak occurs at a point slightly larger than 12; this indicates the octave is slightly compressed in the zeta tuning for 12. The size of a step in octaves is 1/x, and hence the size of the octave in the zeta peak value tuning for Nedo is N/x; if x is slightly larger than N as here with N=12, the size of the zeta tuned octave will be slightly less than a pure octave. Similarly, when the peak occurs with x less than N, we have stretched octaves.<br /> | ||
| Line 269: | Line 276: | ||
For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how 270 completely dominates its neighbors:<br /> | For larger edos, the width of the peak narrows, but for strong edos the height more than compensates, measured in terms of the area under the peak (the absolute value of the integral of Z between two zeros.) Note how 270 completely dominates its neighbors:<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextLocalImageRule: | <!-- ws:start:WikiTextLocalImageRule:42:&lt;img src=&quot;/file/view/plot270.png/219383970/plot270.png&quot; alt=&quot;&quot; title=&quot;&quot; /&gt; --><img src="/file/view/plot270.png/219383970/plot270.png" alt="plot270.png" title="plot270.png" /><!-- ws:end:WikiTextLocalImageRule:42 --><br /> | ||
<br /> | <br /> | ||
Note that for one of its neighbors, 271, it isn't entirely clear which peak value corresponds to the line of real values from +infinity. This can be determined by looking at the absolute value of zeta along other s values, such as s=1 or s=3/4, and in this case the local minimum at 271.069 is the value in question. However, other peak values are not without their interest; the local maximum at 270.941, for instance, is associated to a different mapping for 3.<br /> | Note that for one of its neighbors, 271, it isn't entirely clear which peak value corresponds to the line of real values from +infinity. This can be determined by looking at the absolute value of zeta along other s values, such as s=1 or s=3/4, and in this case the local minimum at 271.069 is the value in question. However, other peak values are not without their interest; the local maximum at 270.941, for instance, is associated to a different mapping for 3.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:21:&lt;h1&gt; --><h1 id="toc3"><a name="Zeta EDO lists"></a><!-- ws:end:WikiTextHeadingRule:21 -->Zeta EDO lists</h1> | ||
If we examine the increasingly larger peak values of |Z(x)|, we find occur with values of x such that Z'(x) = 0 near to integers, so that there is a sequence of <a class="wiki_link" href="/edo">edo</a>s <a class="wiki_link" href="/1edo">1</a>, <a class="wiki_link" href="/2edo">2</a>, <a class="wiki_link" href="/3edo">3</a>, <a class="wiki_link" href="/4edo">4</a>, <a class="wiki_link" href="/5edo">5</a>, <a class="wiki_link" href="/7edo">7</a>, <a class="wiki_link" href="/10edo">10</a>, <a class="wiki_link" href="/12edo">12</a>, <a class="wiki_link" href="/19edo">19</a>, <a class="wiki_link" href="/22edo">22</a>, <a class="wiki_link" href="/27edo">27</a>, <a class="wiki_link" href="/31edo">31</a>, <a class="wiki_link" href="/41edo">41</a>, <a class="wiki_link" href="/53edo">53</a>, <a class="wiki_link" href="/72edo">72</a>, <a class="wiki_link" href="/99edo">99</a>, <a class="wiki_link" href="/118edo">118</a>, <a class="wiki_link" href="/130edo">130</a>, <a class="wiki_link" href="/152edo">152</a>, <a class="wiki_link" href="/171edo">171</a>, <a class="wiki_link" href="/217edo">217</a>, <a class="wiki_link" href="/224edo">224</a>, <a class="wiki_link" href="/270edo">270</a>, <a class="wiki_link" href="/342edo">342</a>, <a class="wiki_link" href="/422edo">422</a>, <a class="wiki_link" href="/441edo">441</a>, <a class="wiki_link" href="/494edo">494</a>, <a class="wiki_link" href="/742edo">742</a>, <a class="wiki_link" href="/764edo">764</a>, <a class="wiki_link" href="/935edo">935</a>, <a class="wiki_link" href="/954edo">954</a>, <a class="wiki_link" href="/1012edo">1012</a>, <a class="wiki_link" href="/1106edo">1106</a>, <a class="wiki_link" href="/1178edo">1178</a>, <a class="wiki_link" href="/1236edo">1236</a>, <a class="wiki_link" href="/1395edo">1395</a>, <a class="wiki_link" href="/1448edo">1448</a>, <a class="wiki_link" href="/1578edo">1578</a>, <a class="wiki_link" href="/2460edo">2460</a>, <a class="wiki_link" href="/2684edo">2684</a>, <a class="wiki_link" href="/3395edo">3395</a>, <a class="wiki_link" href="/5585edo">5585</a>, <a class="wiki_link" href="/6079edo">6079</a>, <a class="wiki_link" href="/7033edo">7033</a>, <a class="wiki_link" href="/8269edo">8269</a>, <a class="wiki_link" href="/8539edo">8539</a>, <a class="wiki_link" href="/11664edo">11664</a> ... of <em>zeta peak edos</em>. This is listed in the On-Line Encyclopedia of Integer Sequences as <a class="wiki_link_ext" href="http://oeis.org/A117536" rel="nofollow">sequence A117536</a>.<br /> | If we examine the increasingly larger peak values of |Z(x)|, we find occur with values of x such that Z'(x) = 0 near to integers, so that there is a sequence of <a class="wiki_link" href="/edo">edo</a>s <a class="wiki_link" href="/1edo">1</a>, <a class="wiki_link" href="/2edo">2</a>, <a class="wiki_link" href="/3edo">3</a>, <a class="wiki_link" href="/4edo">4</a>, <a class="wiki_link" href="/5edo">5</a>, <a class="wiki_link" href="/7edo">7</a>, <a class="wiki_link" href="/10edo">10</a>, <a class="wiki_link" href="/12edo">12</a>, <a class="wiki_link" href="/19edo">19</a>, <a class="wiki_link" href="/22edo">22</a>, <a class="wiki_link" href="/27edo">27</a>, <a class="wiki_link" href="/31edo">31</a>, <a class="wiki_link" href="/41edo">41</a>, <a class="wiki_link" href="/53edo">53</a>, <a class="wiki_link" href="/72edo">72</a>, <a class="wiki_link" href="/99edo">99</a>, <a class="wiki_link" href="/118edo">118</a>, <a class="wiki_link" href="/130edo">130</a>, <a class="wiki_link" href="/152edo">152</a>, <a class="wiki_link" href="/171edo">171</a>, <a class="wiki_link" href="/217edo">217</a>, <a class="wiki_link" href="/224edo">224</a>, <a class="wiki_link" href="/270edo">270</a>, <a class="wiki_link" href="/342edo">342</a>, <a class="wiki_link" href="/422edo">422</a>, <a class="wiki_link" href="/441edo">441</a>, <a class="wiki_link" href="/494edo">494</a>, <a class="wiki_link" href="/742edo">742</a>, <a class="wiki_link" href="/764edo">764</a>, <a class="wiki_link" href="/935edo">935</a>, <a class="wiki_link" href="/954edo">954</a>, <a class="wiki_link" href="/1012edo">1012</a>, <a class="wiki_link" href="/1106edo">1106</a>, <a class="wiki_link" href="/1178edo">1178</a>, <a class="wiki_link" href="/1236edo">1236</a>, <a class="wiki_link" href="/1395edo">1395</a>, <a class="wiki_link" href="/1448edo">1448</a>, <a class="wiki_link" href="/1578edo">1578</a>, <a class="wiki_link" href="/2460edo">2460</a>, <a class="wiki_link" href="/2684edo">2684</a>, <a class="wiki_link" href="/3395edo">3395</a>, <a class="wiki_link" href="/5585edo">5585</a>, <a class="wiki_link" href="/6079edo">6079</a>, <a class="wiki_link" href="/7033edo">7033</a>, <a class="wiki_link" href="/8269edo">8269</a>, <a class="wiki_link" href="/8539edo">8539</a>, <a class="wiki_link" href="/11664edo">11664</a> ... of <em>zeta peak edos</em>. This is listed in the On-Line Encyclopedia of Integer Sequences as <a class="wiki_link_ext" href="http://oeis.org/A117536" rel="nofollow">sequence A117536</a>.<br /> | ||
<br /> | <br /> | ||
| Line 280: | Line 287: | ||
Finally, taking the midpoints of the successively larger normalized gaps between the zeros of Z leads to a list of <em>zeta gap edos</em>. These are <a class="wiki_link" href="/2edo">2</a>, <a class="wiki_link" href="/3edo">3</a>, <a class="wiki_link" href="/5edo">5</a>, <a class="wiki_link" href="/7edo">7</a>, <a class="wiki_link" href="/12edo">12</a>, <a class="wiki_link" href="/19edo">19</a>, <a class="wiki_link" href="/31edo">31</a>, <a class="wiki_link" href="/46edo">46</a>, <a class="wiki_link" href="/53edo">53</a>, <a class="wiki_link" href="/72edo">72</a>, <a class="wiki_link" href="/270edo">270</a>, <a class="wiki_link" href="/311edo">311</a>, <a class="wiki_link" href="/954edo">954</a>, <a class="wiki_link" href="/1178edo">1178</a>, <a class="wiki_link" href="/1308edo">1308</a>, <a class="wiki_link" href="/1395edo">1395</a>, <a class="wiki_link" href="/1578edo">1578</a>, <a class="wiki_link" href="/3395edo">3395</a>, <a class="wiki_link" href="/4190edo">4190</a> ... Since the density of the zeros increases logarithmically, the normalization is to divide through by the log of the midpoint. These edos are listed in the OEIS as <a class="wiki_link_ext" href="http://oeis.org/A117537" rel="nofollow">sequence A117537</a>. The zeta gap edos seem to weight higher primes more heavily and have the advantage of being easy to compute from a table of zeros on the critical line.<br /> | Finally, taking the midpoints of the successively larger normalized gaps between the zeros of Z leads to a list of <em>zeta gap edos</em>. These are <a class="wiki_link" href="/2edo">2</a>, <a class="wiki_link" href="/3edo">3</a>, <a class="wiki_link" href="/5edo">5</a>, <a class="wiki_link" href="/7edo">7</a>, <a class="wiki_link" href="/12edo">12</a>, <a class="wiki_link" href="/19edo">19</a>, <a class="wiki_link" href="/31edo">31</a>, <a class="wiki_link" href="/46edo">46</a>, <a class="wiki_link" href="/53edo">53</a>, <a class="wiki_link" href="/72edo">72</a>, <a class="wiki_link" href="/270edo">270</a>, <a class="wiki_link" href="/311edo">311</a>, <a class="wiki_link" href="/954edo">954</a>, <a class="wiki_link" href="/1178edo">1178</a>, <a class="wiki_link" href="/1308edo">1308</a>, <a class="wiki_link" href="/1395edo">1395</a>, <a class="wiki_link" href="/1578edo">1578</a>, <a class="wiki_link" href="/3395edo">3395</a>, <a class="wiki_link" href="/4190edo">4190</a> ... Since the density of the zeros increases logarithmically, the normalization is to divide through by the log of the midpoint. These edos are listed in the OEIS as <a class="wiki_link_ext" href="http://oeis.org/A117537" rel="nofollow">sequence A117537</a>. The zeta gap edos seem to weight higher primes more heavily and have the advantage of being easy to compute from a table of zeros on the critical line.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:23:&lt;h1&gt; --><h1 id="toc4"><a name="Removing primes"></a><!-- ws:end:WikiTextHeadingRule:23 -->Removing primes</h1> | ||
The <a class="wiki_link_ext" href="http://mathworld.wolfram.com/EulerProduct.html" rel="nofollow">Euler product</a> for the Riemann zeta function is<br /> | The <a class="wiki_link_ext" href="http://mathworld.wolfram.com/EulerProduct.html" rel="nofollow">Euler product</a> for the Riemann zeta function is<br /> | ||
<br /> | <br /> | ||
| Line 301: | Line 308: | ||
Removing 2 leads to increasing adjusted peak values corresponding to the division of 3 (the &quot;tritave&quot;) into 4, 7, 9, 13, 15, 17, 26, 32, 39, 45, 52, 56, 71, 75, 88, 131, 245, 316 ... parts. A striking feature of this list is the appearance not only of <a class="wiki_link" href="/13edt">13edt</a>, the <a class="wiki_link" href="/Bohlen-Pierce">Bohlen-Pierce</a> division of the tritave, but the multiples 26, 39 and 52 also.<br /> | Removing 2 leads to increasing adjusted peak values corresponding to the division of 3 (the &quot;tritave&quot;) into 4, 7, 9, 13, 15, 17, 26, 32, 39, 45, 52, 56, 71, 75, 88, 131, 245, 316 ... parts. A striking feature of this list is the appearance not only of <a class="wiki_link" href="/13edt">13edt</a>, the <a class="wiki_link" href="/Bohlen-Pierce">Bohlen-Pierce</a> division of the tritave, but the multiples 26, 39 and 52 also.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:25:&lt;h1&gt; --><h1 id="toc5"><a name="The Black Magic Formulas"></a><!-- ws:end:WikiTextHeadingRule:25 -->The Black Magic Formulas</h1> | ||
When <a class="wiki_link" href="/Gene%20Ward%20Smith">Gene Smith</a> discovered these formulas in the 70s, he thought of them as &quot;black magic&quot; formulas not because of any aura of evil, but because they seemed mysteriously to give you something for next to nothing. They are based on Gram points and the Riemann-Siegel theta function θ(t). Recall that a Gram point is a points on the critical line where ζ(1/2 + ig) is real. This implies that exp(iθ(g)) is real, so that θ(g)/π is an integer. Theta has an <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Asymptotic_expansion" rel="nofollow">asymptotic expansion</a><br /> | |||
<!-- ws:start:WikiTextMathRule:13: | |||
[[math]]&lt;br/&gt; | |||
\theta(t) \sim \frac{t}{2}\log \frac{t}{2\pi} - \frac{t}{2} - \frac{\pi}{8}+\frac{1}{48t}+ \frac{7}{5760t^3}+\cdots&lt;br/&gt;[[math]] | |||
--><script type="math/tex">\theta(t) \sim \frac{t}{2}\log \frac{t}{2\pi} - \frac{t}{2} - \frac{\pi}{8}+\frac{1}{48t}+ \frac{7}{5760t^3}+\cdots</script><!-- ws:end:WikiTextMathRule:13 --><br /> | |||
<br /> | |||
<br /> | |||
<!-- ws:start:WikiTextHeadingRule:27:&lt;h1&gt; --><h1 id="toc6"><a name="Computing zeta"></a><!-- ws:end:WikiTextHeadingRule:27 -->Computing zeta</h1> | |||
There are various approaches to the question of computing the zeta function, but perhaps the simplest is the use of the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Dirichlet_eta_function" rel="nofollow">Dirichlet eta function</a> which was introduced to mathematics by <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Johann_Peter_Gustav_Lejeune_Dirichlet" rel="nofollow">Johann Peter Gustav Lejeune Dirichlet</a>, who despite his name was a German and the brother-in-law of <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Felix_Mendelssohn_Bartholdy" rel="nofollow">Felix Mendelssohn</a>.<br /> | There are various approaches to the question of computing the zeta function, but perhaps the simplest is the use of the <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Dirichlet_eta_function" rel="nofollow">Dirichlet eta function</a> which was introduced to mathematics by <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Johann_Peter_Gustav_Lejeune_Dirichlet" rel="nofollow">Johann Peter Gustav Lejeune Dirichlet</a>, who despite his name was a German and the brother-in-law of <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Felix_Mendelssohn_Bartholdy" rel="nofollow">Felix Mendelssohn</a>.<br /> | ||
<br /> | <br /> | ||
The zeta function has a <a class="wiki_link_ext" href="http://mathworld.wolfram.com/SimplePole.html" rel="nofollow">simple pole</a> at z=1 which forms a barrier against continuing it with its <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Euler_product" rel="nofollow">Euler product</a> or <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Dirichlet_series" rel="nofollow">Dirichlet series</a> representation. We could subtract off the pole, or multiply by a factor of (z-1), but at the expense of losing the character of a Dirichlet series or Euler product. A better method is to multiply by a factor of (1-2^(1-z)), leading to the eta function:<br /> | The zeta function has a <a class="wiki_link_ext" href="http://mathworld.wolfram.com/SimplePole.html" rel="nofollow">simple pole</a> at z=1 which forms a barrier against continuing it with its <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Euler_product" rel="nofollow">Euler product</a> or <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Dirichlet_series" rel="nofollow">Dirichlet series</a> representation. We could subtract off the pole, or multiply by a factor of (z-1), but at the expense of losing the character of a Dirichlet series or Euler product. A better method is to multiply by a factor of (1-2^(1-z)), leading to the eta function:<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextMathRule: | <!-- ws:start:WikiTextMathRule:14: | ||
[[math]]&lt;br/&gt; | [[math]]&lt;br/&gt; | ||
\eta(z) = (1-2^{1-z})\zeta(z) = \sum_{n=1}^\infty (-1)^{n-1} n^{-z}&lt;br /&gt; | \eta(z) = (1-2^{1-z})\zeta(z) = \sum_{n=1}^\infty (-1)^{n-1} n^{-z}&lt;br /&gt; | ||
= \frac{1}{1^z} - \frac{1}{2^z} + \frac{1}{3^z} - \frac{1}{4^z} + \cdots&lt;br/&gt;[[math]] | = \frac{1}{1^z} - \frac{1}{2^z} + \frac{1}{3^z} - \frac{1}{4^z} + \cdots&lt;br/&gt;[[math]] | ||
--><script type="math/tex">\eta(z) = (1-2^{1-z})\zeta(z) = \sum_{n=1}^\infty (-1)^{n-1} n^{-z} | --><script type="math/tex">\eta(z) = (1-2^{1-z})\zeta(z) = \sum_{n=1}^\infty (-1)^{n-1} n^{-z} | ||
= \frac{1}{1^z} - \frac{1}{2^z} + \frac{1}{3^z} - \frac{1}{4^z} + \cdots</script><!-- ws:end:WikiTextMathRule: | = \frac{1}{1^z} - \frac{1}{2^z} + \frac{1}{3^z} - \frac{1}{4^z} + \cdots</script><!-- ws:end:WikiTextMathRule:14 --><br /> | ||
<br /> | <br /> | ||
The Dirichlet series for the zeta function is absolutely convergent when s&gt;1, justifying the rearrangement of terms leading to the alternating series for eta, which converges conditionally in the critical strip. The extra factor introduces zeros of the eta function at the points 1 + 2πix/ln(2) corresponding to pure octave divisions along the line s=1, but no other zeros, and in particular none in the critical strip and along the critical line. The convergence of the alternating series can be greatly accelerated by applying <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Euler_summation" rel="nofollow">Euler summation</a>.<br /> | The Dirichlet series for the zeta function is absolutely convergent when s&gt;1, justifying the rearrangement of terms leading to the alternating series for eta, which converges conditionally in the critical strip. The extra factor introduces zeros of the eta function at the points 1 + 2πix/ln(2) corresponding to pure octave divisions along the line s=1, but no other zeros, and in particular none in the critical strip and along the critical line. The convergence of the alternating series can be greatly accelerated by applying <a class="wiki_link_ext" href="http://en.wikipedia.org/wiki/Euler_summation" rel="nofollow">Euler summation</a>.<br /> | ||
<br /> | <br /> | ||
<!-- ws:start:WikiTextHeadingRule: | <!-- ws:start:WikiTextHeadingRule:29:&lt;h1&gt; --><h1 id="toc7"><a name="Links"></a><!-- ws:end:WikiTextHeadingRule:29 -->Links</h1> | ||
<a class="wiki_link_ext" href="http://front.math.ucdavis.edu/0309.5433" rel="nofollow">X-Ray of Riemann zeta-function</a> by Juan Arias-de-Reyna<br /> | <a class="wiki_link_ext" href="http://front.math.ucdavis.edu/0309.5433" rel="nofollow">X-Ray of Riemann zeta-function</a> by Juan Arias-de-Reyna<br /> | ||
<a class="wiki_link_ext" href="http://terrytao.wordpress.com/2009/07/12/selbergs-limit-theorem-for-the-riemann-zeta-function-on-the-critical-line/" rel="nofollow">Selberg's limit theorem</a> by Terence Tao <a class="wiki_link_ext" href="http://www.webcitation.org/5xrvgjW6T" rel="nofollow">Permalink</a></body></html></pre></div> | <a class="wiki_link_ext" href="http://terrytao.wordpress.com/2009/07/12/selbergs-limit-theorem-for-the-riemann-zeta-function-on-the-critical-line/" rel="nofollow">Selberg's limit theorem</a> by Terence Tao <a class="wiki_link_ext" href="http://www.webcitation.org/5xrvgjW6T" rel="nofollow">Permalink</a></body></html></pre></div> | ||