Some Unrelated Topics of Interest

A Publicon Demonstration Document

Michael Trott

Programming and Design Analyst

Wolfram Research, 100 Trade Center Dr, Champaign, IL 61820

info@wolfram.com

This paper offers a demonstration of the appearance and general document features of the Default style sheet. Although the topics in this paper were written specifically for Publicon, it is hoped they will nonetheless be regarded with interest by the physics community. © 2004 Wolfram Research, Inc.

High-precision Value for the Quartic Anharmonic Oscillator Ground State

Introduction

As is well known, only a very limited number of one-dimensional potentials allow for an exact solution of the Schrödinger equation1. This means that for many model potentials one has to resort to numerical solution methods. For judging their accuracy, reliability, and speed it is important to have high-precision values of certain nonexactly solvable potentials. The most investigated of such potentials is the quartic anharmonic oscillator [1-11][12-19] described by

-ψ _ k^''(z) + z^4 ψ _ k(z) = ε _ k ψ _ k(z)(1)

The eigenfunctions to the eigenvalues typeset structure decay exponentially for typeset structure.

The Hill Determinant Method

A classical method to solve Sturm-Liouville problems of type 1 is to calculate the eigenvalues of a truncated version of the corresponding Hill determinant. Using the harmonic oscillator basis typeset structure we write typeset structure where

-φ _ k^''(z) + z^2 φ _ k(z) = ϵ _ k φ _ k(z)(2)
φ _ n(z) = 1/(π^(1/2) 2^n n !)^(1/2) e^-z^2/2 H _ n(z)(3)

Forming the matrix elements typeset structure. For typeset structure we obtain

h _ (m, n) = {                0                                                                ...                                             m, n                   m, n + 2         m, n + 4  else(4)

A rough estimation shows that one obtains about 0.2 digits per harmonic oscillator state. So by taking into account the first 500 eigenstates and carrying out the calculation with about five thousand digits one obtains about 120 reliable digits for typeset structure. (This calculation takes about 20 minutes on a 2000 vintage workstation using Mathematica 4 [20].)

ε _ 0 = 1.060362090484182899647046016692663545515208728528977933216245241695943563044344421126896299134671703510546244358582525580982763829 ...

The Hill determinant approach allows in addition to the calculation of the eigenvalues, the calculation of the eigenvectors. The following graphic visualizes the matrix of eigenvectors of typeset structure. The graphic shows that the lowest eigenfunctions are quite similar to the harmonic oscillator eigenfunctions. Higher states are complicated mixtures of harmonic oscillator states. The overall “checkerboard”-like structure results from the fact that the contribution of the antisymmetric (symmetric) harmonic oscillator states to the symmetric (antisymmetric) anharmonic oscillator states is identical zero. The very high states are dominated by truncation effects and do not correctly mimic the anharmonic oscillator states.

[Graphics:HTMLFiles/DefaultSample_14.gif]

Figure 1: The matrix of eigenvectors of typeset structure.

The New Algorithm

To get a very high-precision approximation of  

-ψ^''(x) + z^4 ψ(x) = λ ψ(x)(5)

we start with the series expansion

ψ(x) = y _ n(x) = Underoverscript[∑, k = 0, arg3] a _ k(λ) x^k(6)

For the ground state we choose (ignoring normalization) typeset structure. For “suitable chosen” typeset structure we then find high-precision approximations for the zeros of typeset structure and typeset structure. These zeros then bound typeset structure from below and above.   

Using the differential equation one obtains the following recursion relation for the typeset structure:

a _ k(λ)    = (a _ (m - 6)(λ) - λ a _ (m - 2)(λ))/(m^2 - m)(7)

For large typeset structure (typeset structure) we want the function typeset structure to vanish as typeset structure. For a λ smaller than the smallest possible λ, the function typeset structure will not have a zero, but the function typeset structure will have a zero for a certain typeset structure. For a λ larger than the smallest possible λ, the function typeset structure will have a zero, but the function typeset structure will not have a zero for a certain typeset structure. This allows to find a bounding interval for typeset structure. The next two graphics show typeset structure and typeset structure for 10 equidistant values for λ from the interval typeset structure to visualize this bounding process. (For more details, see [21].)

[Graphics:HTMLFiles/DefaultSample_39.gif]

Figure 2: typeset structure and typeset structure for 10 equidistant values for λ from the interval typeset structure.

It is straightforward to implement the calculation of the bounding interval for typeset structure in Mathematica in a three-line program (see [22]). Using FindRoot we calculate high-precision values for the zeros of typeset structure and typeset structure.

λBounds[n_, ξ_, opts___] :=
   Function[f, λ /. FindRoot[f[n, λ, ξ] == 0,
       {λ, 106/100, 107/100}, opts]] /@ {y, yPrime};

The calculation of typeset structure and typeset structure is also straightforward based on a recursive calculations of the typeset structure.

y[n_, λ_Real, ξ_] := Module[{a6, a4, a2, ak, σ},
   {a6, a4, a2} = {1, -λ/2, λ^2/24};
    σ = a6 + a4*ξ^2 + a2*ξ^4;
    Do[ak = a6 - λ*(a2/(k*(k - 1)));
      {a6, a4, a2} = {a4, a2, ak};
      σ = σ + ak*ξ^k, {k, 6, n, 2}]; σ]

y[n_, λ_Real, ξ_] := Module[{a6, a4, a2, ak, σ},
   {a6, a4, a2} = {1, -λ/2, λ^2/24};
    σ = a6 + 2*ξ*a4 + 4*ξ^3*a2;
    Do[ak = a6 - λ*(a2/(k*(k - 1)));
      {a6, a4, a2} = {a4, a2, ak};
      σ = σ + k*ak*ξ^(k - 1), {k, 6, n, 2}]; σ]

Calculating now λBounds[16000, 16, startingValues, WorkingPrecision->6000, AccuracyGoal->600, MaxIterations->100] (where startingValues has been obtained from a call to λBounds recursively, one gets in a few minutes a 1184 digit approximation to the ground state energy of the quartic anharmonic oscillator.

ε _ 0 = 1.06036209048418289964704601669266354551520872852897793321624524169594356304434 ... 5479719467596466796613467688586437952654519560568286715958338884743467012042420714919290048732 ...

Statistical analysis of the number does not show any regularity.

Summary

A power series based approach to the high-precision calculation of the ground state of the anharmonic oscillator was presented. Mathematica code to carry out the calculation, as well as results were given. The method can straightforwardly be used to calculate ten thousands of digits of the quartic anharmonic, as well as other anharmonic oscillators. Work concerning the application of the method to higher states is in progress.

All calculations and visualizations have been carried out in Mathematica 4.

Local Density of States for the Harmonic Oscillator

The next graphic shows the local density of states typeset structure for the harmonic oscillator [21].

[Graphics:HTMLFiles/DefaultSample_51.gif]

Figure 3: The local density of states typeset structure for the harmonic oscillator.

Acknowledgments

The author would like to thank André Kuzniarek for making a prerelease version of the Publicon typesetting system. This work was supported by Wolfram Research, Inc.

Are Brillouin Zones of High Order Fractal?

Introduction

Brillouin zones are among the most popular objects a solid state physicist deals with [23-26]. Despite the fundamental importance for the explanation of most properties of crystalline solids, Brillouin zones as an own subject have rarely been investigated (the only ones we are aware of are [27-30]). For electronic properties, mostly the low order Brillouin zones matter, as an ownstanding subject the high order Brillouin zones are interesting. Mathematically, the typeset structureth Brillouin zone is the set of points that a line to them crosses exactly typeset structure bisector planes. In computational geometry a typeset structureth order Brillouin zones is also called typeset structureth degree Voronoi region or nth nearest point Voronoi diagrams. The most important fact for high order Brillouin zones is that their shape approaches that of a thin spherical shell and their volume is a constant. Here, for the first time we report on some computational results of higher order Brillouin zones. All calculations and visualizations were done with Mathematica 4 [20].

Recursive definition of Brillouin zones: Given a lattice typeset structure in typeset structurewith lattice points typeset structure (typeset structure being a multiindex) the first Brillouin zone typeset structure is the closure of the set of all points typeset structure such that typeset structuretypeset structure for all typeset structure!=typeset structure. The typeset structureth order Brillouin zones is the closure of the set of all points typeset structure such that typeset structuretypeset structure for all typeset structure!=typeset structure and typeset structure.

2D Hexagonal Lattice

Figure 4 shows the first twenty Brillouin zones of a 2D hexagonal lattice. It is interesting to observe that the first, third, and fourth Brillouin zones have the shape of an hexagon. For higher orders the shape becomes much more complicated.

[Graphics:HTMLFiles/DefaultSample_74.gif]

Figure 4: The first twenty Brillouin zones of a 2D hexagonal lattice.

Figure 5 shows the 200th Brillouin zone in one symmetry unit (inside an angle of typeset structure). One sees many small and a few quite large Landsberg zones. The distribution typeset structure of the area of the Landsberg zones in the limit typeset structure might be an interesting subject to study.

[Graphics:HTMLFiles/DefaultSample_78.gif]

Figure 5: The 200th Brillouin zone in one symmetry unit (inside an angle of typeset structure).

A good numerical fit to number of faces (line segments) typeset structure of the typeset structureth Brillouin zone is typeset structure. Figure 6 shows the circumference of the Brillouin zones normalized to the circumference of a circle with the same radius. Does the ratio approach a finite value in the limit typeset structure?

[Graphics:HTMLFiles/DefaultSample_84.gif]

Figure 6: The circumference of the Brillouin zones normalized to the circumference of a circle with the same radius.

3D Cubic Lattices

In [31] we gave a complete implementation for the effective calculation of higher order Brillouin zones of the three cubic lattice in typeset structure. Figure 7 shows the (outside of) 15th Brillouin zone for the simple cubic, Figure 8 shows the 18th for the face-centered cubic and Figure 9 shows the 10th for the body-centered cubic lattice. The higher order Brillouin zones show quite complicated behavior. Large faces alternate with small ones, relatively plane regions alternate with quite structured ones. The appearance of the typeset structureth Brillouin zone is typically completely independent of the appearance of the typeset structureth Brillouin zone.

[Graphics:HTMLFiles/DefaultSample_88.gif]

Figure 7: The (outside of) 15th Brillouin zone for the simple cubic.

[Graphics:HTMLFiles/DefaultSample_89.gif]

Figure 8: The (outside of) 18th Brillouin zone for the face-centered cubic.

[Graphics:HTMLFiles/DefaultSample_90.gif]

Figure 9: The (outside of) 10th Brillouin zone for the body-centered cubic lattice.

Using this implementation we analysed various properties of higher order Brillouin zones. Figure 10 shows the area of the Brillouin zones normalized to the area of a sphere of the same volume. The order of the three curves from the bottom is sc, fcc, bcc.

[Graphics:HTMLFiles/DefaultSample_91.gif]

Figure 10: The area of the Brillouin zones normalized to the area of a sphere of the same volume.

The last result to be given here is the number of faces of the Brillouin zones. By a face we mean any connected planar part of the outside facing side of a Brillouin zone (point contacts separate faces). Table 1 gives the results for the first 25 Brillouin zones.

Table 1. The results for the first 25 Brillouin zones.

nscfccbcc
 
161412
2127248
3729630

Table 2. The results for the first 25 Brillouin zones.

nsc
 
16

Table 3. The results for the first 25 Brillouin zones.

nscfcc
 
1614
21272

Table 4. The results for the first 25 Brillouin zones.

nscfccbccbcc
 
16141212
212724848
372963030

Summary

Preliminary results about some computational results about higher order Brillouin zones have been presented. Further work is in progress and will be published elsewhere.

All calculations and visualizations have been carried out in Mathematica 4.

Acknowledgments

The author would like to take the opportunity to thank Amy Young for assistance in typesetting this paper.

Appendix: Index Gymnastics

The technique of extracting the content from geometric (tensor) equations by working in component notation and rearranging indices as required2. Index gymnastics is a fundamental component of special and general relativity [32]. Examples of index gymnastics include

S^αβ _ γ = g^βμ S^α _ μγ(A.1)
S^α _ μγ = g _ μβ S^αβ _ γ(A.2)
A^2 = A^α A _ α(A.3)
g _ αβ g^βγ = δ _ α^γ(A.4)
N^α _ β^(, γ) = N^α _ (β, μ) g^μγ(A.5)
(R _ α M _ β) _ (, γ) = R _ (α, γ) M _ β + R _ α M _ (β, γ)(A.6)
F _ [αβ] = 1/2 (F _ αβ - F _ βα)(A.7)
F _ (αβ) = 1/2 (F _ αβ + F _ βα)(A.8)

[32], where typeset structure is the metric tensor, typeset structure is the Kronecker delta, typeset structure is a comma derivative, typeset structure is the antisymmetric tensor part, and typeset structure is the symmetric tensor part3.

Notes

1 This footnote was not added by the author, but was placed for the purpose of demonstration.

2 From MathWorld, http://mathworld.wolfram.com

3 See also: Index, Index Lowering, Index Raising, Tensor.

References

[1]  Tymczak, C. J., Japaridze, G. S., Handy, C. R., & Wang, X. Q. (1998, March 14). New perspective on inner product quantization. Phys. Rev. Lett., 80(17), 3673-3677.

[2]  Graffi, S., & Grecchi, V. (1973, November 15). Rayleigh-Ritz Method, Secular Determinant, and Anharmonic Oscillators. Phys. Rev. D, 8(10), 3487-3492.

[3]  Sergeev, A. V., & Goodson, D. Z. (1998, March 14). Bloch electrons in a magnetic field. J. Phys., A(31), 4301-4301.

[4]  Suzuki, J. (1999, March 14). A new ambient pressure organic superconductor based on BEDT-TTF with 10.4-K higher than 10K. J. Phys., A(32), 183-187.

[5]  Chebotarev, L. V. (1997, Dec 12). Parabolic connection formulae in quantum mechanics. Ann. Phys., 255(2), 305-332.

[6]  Shanley, P. E. (1988, March 14). Large-order analysis of the convergent renormalized strong-coupling perturbation theory for the quartic anharmonic oscillator. Ann. Phys., 186(3), 292-325.

[7]  Hatsuda, T., Kunihiro, T., & Tanaka, T. (1997, April 28). Optimized Perturbation Theory for Wave Functions of Quantum Systems. Phys. Rev. Lett., 78(17), 3229-3230.

[8]  Weniger, E. J. (1996, September 30). Construction of the Strong Coupling Expansion for the Ground State Energy of the Quartic, Sextic, and Octic Anharmonic Oscillator via a Renormalized Strong Coupling Expansion. Phys. Rev. Lett., 77(14), 2859-2862.

[9]  Ho, K. C., et al. (1996, March 12). Study of quantum anharmonic oscillators by state-dependent diagonalization. Phys. Rev. A., 53(3), 1280-1284.

[10]  Čizek, J., Weniger, E. J., Bracken, P., & Špiro, V. (1996, Dec 12). Effective characteristic polynomials and two-point Padé approximants as summation techniques for the strongly divergent perturbation expansions of the ground state energies of anharmonic oscillators. Phys.Rev., 53(3), 2925-2939.

[11]  Bay, K., & Lay, W. (1997, Dec 12). The spectrum of the quartic oscillator. J. Math. Phys., 38(5), 2127-2131.

[12]  Voros, A. (1999, March 14). Magnetic field dependence of the cyclotron effective mass in the organic superconductor. J. Phys., 32(A), 5993-6000.

[13]  Khan, P. B., & Zarmi, Y. (1999, Dec 12). Properties of new organic conductors. J. Math. Phys., 40(1), 4658-4659.

[14]  Antonsen, F. (1997, Dec 12). Quantum theory in curved spacetime using the Wigner function. Phys. Rev. D., 56(2), 920-935.

[15]  Skála, L., Čizek, J., Kapsa, V., & Weniger, E. J. (1997, December). Large-order analysis of the convergent renormalized strong-coupling perturbation theory for the quartic anharmonic oscillator. Phys. Rev. A, 56(6), 4471-4476.

[16]  Meißner, H., & Steinborn, E. O. (1997, March 14). Large-order analysis of the convergent renormalized strong-coupling perturbation theory for the quartic anharmonic oscillator. Phys. Rev., 56(A), 4471-4476.

[17]  Delabaere, E., & Pham, F. (1997, Dec 12). Unfolding the quartic oscillator. Ann. Phys, 261(2), 180-218.

[18]  Fernández, F., & Guardiola, R. (1993, Dec 12). Accurate eigenvalues and eigenfunctions for quantum-mechanical anharmonic oscillators. J. Phys. A., 26(23), 7169-7180.

[19]  Vinette, F., & Čizek, J. (1991, March 14). Spinless fermions on frustrated lattices in a magnetic field. J. Math. Phys., 26(3), 3392-3396.

[20]  Wolfram, S. (1999). The Mathematica Book. Champaign: Cambridge University Press and Wolfram Media.

[21]  Trott, M. (2000). The Mathematica GuideBook: Mathematics in Mathematica. New York: Springer-Verlag.

[22]  Trott, M. (2000). The Mathematica GuideBook: Programming in Mathematica. New York: Springer-Verlag.

[23]  Brillouin, L. (1946). Wave Propogation in Periodic Structures. New York: McGraw-Hill.

[24]  Ashcroft, N. W., & Mermin, N. D. (1976). Solid State Physics. Philadelphia: Saunders College.

[25]  Crackwell, A. P., & Wong, K. C. (1973). The Fermi Surface. Oxford: Clarendon Press.

[26]  Landsberg, P. T. (1969). Solid State Theory. London: Wiley.

[27]  Bieberbach, L. (1939, Dec 12). Über die Inhaltsgleichheit der Brillouinschen Zonen. Monatsh. Math. Phys., 48(1), 509-515.

[28]  Jones, G. A. (1984, Dec 12). An algorithm for machine calculation of complex Fourier series. Bull. London Math. Soc., 16(1), 241-242.

[29]  Skriganov, M. M. (1984, March 14). Mulliken-Wolfsberg-Helmholtz band structure of di-tetramethyltetraselenafulvalene-X typeset structure: Role of the basis set. Zap. Nautn. Sem.Leningrad. Otdel. Mat. Inst. Steklov., 134(3), 206-210.

[30]  Veerman, J. J., et al. (2000, March 14). On Brillouin zones. Comm. Math. Phys., 212(3), 725-744.

[31]  Trott, M. (2000). The Mathematica GuideBook: Graphics in Mathematica. New York: Springer-Verlag.

[32]  Misner, C. W., Thorne, K. S., & Wheeler, J. A. (1973). Gravitation. San Francisco: W. H. Freeman.