Computational Physics - M. Jensen Episode 2 Part 8 docx

... . . . . . . . . . . 18 4 11 Monte Carlo methods in statistical physics 18 7 11 .1 Phase transitions in magnetic systems . . . . . . . . . . . . . . . . . . . . . . . 18 7 11 .1. 1 Theoretical background ... . 18 7 11 .1. 2 The Metropolis algorithm . . . . . . . . . . . . . . . . . . . . . . . . . 19 3 11 .2 Program example . . . . . . . . . . . . . . . . . . . . . . . . . . . . ....

Ngày tải lên: 07/08/2014, 12:22

... 0.935762E -1 3 -0 .34 21 3 4E-04 10 0 40.0 0. 424 835E -1 7 -0 .2 21 0 33E+ 01 127 50.0 0 .19 28 75E- 21 -0 .833851E+05 15 5 60.0 0.875651E -2 6 -0 .850381E+09 17 1 70.0 0.397545E-30 NaN 17 1 80.0 0 .18 0485E-34 NaN 17 1 90.0 ... -0 .3066 811 1E-04 10 0 40.000000 0. 424 83543E -1 7 -0 . 316 57 319 E+ 01 127 50.000000 0 .19 28 7498E- 21 0 .11 0 729 33E+05 15 5...

Ngày tải lên: 07/08/2014, 12:22

... equation that (3 .18 ) (3 .19 ) and (3. 20) These equations have the solution (3. 21) and (3. 22) yielding To determine , we require in the last equation that (3. 23) (3. 24) 42 CHAPTER 3. NUMERICAL DIFFERENTIATION Assume ... e.g., (3 .11 ) We could also define five-points formulae by expanding to two steps on each side of . Using a Taylor expansion around in a region we have (3...

Ngày tải lên: 07/08/2014, 12:22

... 20.0857 04 20.085539 20.085537 20.25 046 7 20.085537 4. 0 54. 643 6 64 54. 598605 54. 59 815 5 54. 59 815 1 54. 711 789 54. 59 815 0 5.0 14 8.536878 14 8. 41 4 396 14 8. 41 3 172 14 8. 41 3 1 61 150.635056 14 8. 41 3 159 Table 3 .1: Result ... 7.38905 610 0.00000003 3.0 20 .10 228 045 20.0897 217 6 20.08658307 20.08553692 0.00000009 4. 0 54. 643 66366 54. 60952560 54. 60099375 54. 5...

Ngày tải lên: 07/08/2014, 12:22

... the Gauss-Legendre method. 1 1.305 3.334 2 6 .74 7 7. 473 3 16 .030 10 .954 4 28.330 13 .463 5 42.556 14 .77 6 6 57. 444 14 .77 6 7 71 . 670 13 .463 8 83. 970 10 .954 9 93.253 7. 473 10 98.695 3.334 8.3.6 Applications ... Gauss-Legendre 10 1. 8 210 20 1. 214 025 0 .14 60448 20 0. 912 678 0.6098 97 0. 2 17 80 91 40 0. 478 456 0.333 71 4 0. 219 3834 10 0 0. 273 724 0.2 312...

Ngày tải lên: 07/08/2014, 12:22

... ran1 ran2 ran3 0. 0-0 .1 1 013 9 91 9 38 10 47 0. 1- 0 .2 10 02 10 09 10 40 10 30 0. 2-0 .3 989 999 10 30 993 0. 3-0 .4 939 960 10 23 937 0. 4-0 .5 10 38 10 01 1002 992 0. 5-0 .6 10 37 10 47 10 09 10 09 0. 6-0 .7 10 05 989 10 03 ... 1. 59064E-02 10 00 7 .83 486 E- 01 5 .14 102E-03 10 000 7 .85 488 E- 01 1.60 311 E-03 10 0000 7 .85 009E- 01 5. 087 45E-04...

Ngày tải lên: 07/08/2014, 12:22

... spin [ j ] spin [ j +1] ; } 10 .4. THE METROPOLIS ALGORITHM AND DETAILED BALANCE 18 1 Figure 10 .5: Probability distribution for one walker after 10 , 10 0 and 10 00 steps. 10 .5. PHYSICS PROJECT: SIMULATION ... statistical physics at a given temperature, we need a probability distribution (11 .2) 18 7 10 .3. MICROSCOPIC DERIVATION OF THE DIFFUSION EQUATION 17 7 with the obv...

Ngày tải lên: 07/08/2014, 12:22

... two-dimensional lattice we obtain the following partition function (11 .22 ) and mean energy (11 .23 ) 20 4 CHAPTER 12. QUANTUM MONTE CARLO METHODS which can be rewritten as ( 12. 5) In general, the integrals involved ... the one-dimensional Ising model with spins with free ends (FE) and periodic boundary conditions (PBC). Number spins up Degeneracy Energy (FE) Energy (PBC) Magnetization 2...

Ngày tải lên: 07/08/2014, 12:22

... 7.71 826 E-04 1 .20 000E+00 1.08667E+00 2. 233 89E-01 1.49462E- 03 1 .30 000E+00 1.17168E+00 4.78446E-01 2. 18 734 E- 03 1.40000E+00 1 .26 37 4E+00 8.55 524 E-01 2. 924 93E- 03 1.50000E+00 1 .38 897E+00 1 .30 720 E+00 3. 61553E- 03 ... of 2 was used in order to obtain an acceptance of . 5.00000E-01 2. 06479E+00 5.78 739 E+00 7.60749E- 03 6.00000E-01 1.50495E+00 2. 32 7 82E+00 4....

Ngày tải lên: 07/08/2014, 12:22

... 9.898985E-01 2. 949 052E+00 4. 86 622 3E+00 6.739916E+00 8.56 844 2E+00 100 9.9 748 93E-01 2. 98 744 2E+00 4. 96 727 7E+00 6.936913E+00 8.89 628 2E+00 20 0 9.993715E-01 2. 996864E+00 4. 991877E+00 6.9 843 35E+00 8.9 743 01E+00 40 0 ... the potential is symmetric with nm eV nm [eV] 86 42 0 -2 -4 -6 -8 0 -1 0 -2 0 -3 0 -4 0 -5 0 -6 0 Figure 12. 5: Plot of f...

Ngày tải lên: 07/08/2014, 12:22

... second-order equation goes typically like (14 .2) A well-known second-order equation is Newton’s second law (14.3) where is the force constant. ODE depend only on one variable, whereas 25 5 26 2 CHAPTER ... derivative in our Taylor expansion. We have then (14 .22 ) The second derivative can be rewritten as (14 .23 ) and we can rewrite Eq. (14.14) as (14 .24 ) 25 0 CHAPTER 13. EIGENSYS...

Ngày tải lên: 07/08/2014, 12:22

... errmax ; double yout [ 2] , y_h [ 2] , y_m [2] , y1 [ 2] , y2 [ 2 ] , delt a [ 2] , ys cal [ 2 ]; const double eps =1.0e 6; const double s a fet y =0.9; const double er rcon =6. 0e 4; const double ... phi and v / yout [0]= yin [ 0 ] +1.0 / 6 .0 ( k1 [0] +2 k2 [0] +2 k3 [0]+ k4 [ 0]) ; yout [1]= yin [ 1 ] +1.0 / 6 .0 ( k1 [1] +2 k2 [1] +2 k3 [1]+ k4 [ 1]) ; } 28 0 CHAPTER 1...

Ngày tải lên: 07/08/2014, 12:22

... following condition at (15 .22 ) We can modify this expression by normalizing the function . Then Eq. (15 .22 ) becomes at (15 .23 ) For an arbitary value of the eigenvalue Eq. (15 .22 ) will not be satisfied. ... simplified as (16 .26 ) and (16 . 27 ) The one-dimensional diffusion equation can then be rewritten in its discretized version as (16 .28 ) Defining results in the explicit scheme (16...

Ngày tải lên: 07/08/2014, 12:22

... Possible boundary and initial conditions with are (16. 58) 3 12 CHAPTER 16. PARTIAL DIFFERENTIAL EQUATIONS programs/chap16/program2.f90 ! Program to solve the 2 dim Laplace equation usi ng i t e r a t i ... as (16.47) or in matrix-vector form as (16. 48) where the vector is the same as defined in the implicit case while the matrix is (16.49) 16 .2. 5 Non-linear terms and implementation...

Ngày tải lên: 07/08/2014, 12:22

... Computational Physics (Wiley, 199 7). [7] E.W. Schmid, G. Spitz and W. Lösch, Theoretische Physik mit dem Personal Computer, (Springer, 198 7). [8] J .M. Thijssen, Computational Physics, (Springer, 199 9). [9] ... (Springer, 199 3). [4] S.E. Koonin and D. Meredith, Computational Physics, (Addison Wesley, 199 0). [5] F.S. Acton, Numerical Methods that work, (Harper and Row, 197 0)....

Ngày tải lên: 07/08/2014, 12:22