Computerized Nuclear Coordinate Calculations
A Perl program is available to calculate the position of one added nucleon, based on three or four existing protons and neutrons that were published for all elements. All 118 elements are thoroughly described and tabulated in my hardcover book:
Charge distributions on the nuclei
published in 2022 on Amazon. Also available as a 19 volume eBook on the Kindle platform.
(Hardcover $101.57, eBook volumes $3.44 to $33.44)
Free download file names
Uscan234-84.pl : this is the Perl program for uranium or any element.
fo-pn-33.txt : this file has 17,000 text lines describing the coordinates of every proton and neutron in all 118 elements. It uses the format for the Perl program. It was built directly from the 118 Tables in the digital book volumes.
Additional text data files needed to run the Perl program:
uranium234-baryons83.txt : make this file yourself by copying the 234 lines of text for uranium and pasting those lines into the file.
precise191-rev79.txt : This file should have the precision corrections for the rounded values of coordinates in the books' Tables. Here is the file content you can create:
3.41 3.414213562
4.82 4.8284271247
6.24 6.2426406871
7.65 7.6568542494
9.07 9.0710678118
10.48 10.48528137423
11.89 11.89949493661
13.31 13.31370849898
4.41 4.414213562
2.41 2.414213562
-3.41 -3.414213562
-4.82 -4.8284271247
-6.24 -6.2426406871
-7.65 -7.6568542494
-9.07 -9.0710678118
-10.48 -10.48528137423
-11.89 -11.89949493661
-13.31 -13.31370849898
-4.41 -4.414213562
-2.41 -2.414213562
optional file for interactive successive approximation subroutine, from a menu selection
goals-237.txt . This example is for francium, not uranium. Three baryons form the base, upon which, a sought-after baryon will have coordinates determined. The first three numbers are the baryon numbers of a base that is an equilateral triangle. The next three numbers are the "guess of the coordinates" of the sought-after baryon, which will start the successive approximation routine. After the first line of text, any comments will be ignored by the Perl program:
144 209 165 -4 2 4
Fr-227 goal find coordinates for Neutron #226
Output file:
output-237.txt : The Perl program appends and results that the menu choice will produce. Some Perl menu choices do not make that output file. The screen of a command prompt "terminal" gives all results for all menu selections. The terminal is like the old Microsoft DOS operating system that came before the Windows GUI.
How to use the Perl program, Part 1
I use Windows 10. You need to edit the Perl file so the file names are in a certain directory. You can customize which directory the files are in. Most of the edits are near the top of the Perl program, but some are elsewhere, so look for the path like C:/Perl/source_magnetism/rev17 to edit for accuracy on your file system. Here is what my "terminal" shows:
C:\Perl\source_magnetism\rev17>dir
Volume in drive C is Windows
Volume Serial Number is 68BA-F701
Directory of C:\Perl\source_magnetism\rev17
05/19/2022 08:07 PM <DIR> .
05/19/2022 08:07 PM <DIR> ..
04/20/2022 11:55 AM 67 goals-226.txt
04/19/2022 10:44 AM 410 precise191-rev79.txt
05/19/2022 09:07 AM 4,322 uranium234-baryons83.txt
05/19/2022 09:10 AM 30,189 UScan234-83.pl
05/19/2022 12:26 PM 30,226 UScan234-84.pl
05/19/2022 08:06 PM 30,260 UScan234-85.pl
05/22/2022 01:04 PM 30,213 UScan234-86.pl
7 File(s) 125,687 bytes
2 Dir(s) 747,424,788,480 bytes free
How to use the program, Part 2. Start the program to see the first menu:
C:\Perl\source_magnetism\rev17>UScan234-84.pl
the file loaded in as points, before Table is created.
Make a small table, baryon number and xyz before corrections for rounding in Table 58
Welcome to the NukeMaker Software Menu from rev17 RnScan220-81.pl for 234 baryons. Find Fr
Do you want to use 4 or 3 points to place another baryon or
will you enter 6 to use an arbitrary triangle ?
Notice that 234 baryons are used. That number changes for each element. Edit the Perl file to change the mass number A for the chemical element of your choice. That is the number of lines in the file with a name like uranium234-baryons83.txt :
my $A=234; # number of rows in Table 92
Use the menu item 3 for this example, for an equilateral triangle as a base. You need to choose which three baryon numbers will be the base for the fourth proposed baryon. See Volume 17 of my book, Table 92 for uranium. Charge distributions on the nuclei: Data set seventeen: U to Fm 1, Folmsbee, Alan - Amazon.com . The price is $33.44. Figure 1707 shows the plan for U-234 with numbered baryons. This can help, but an even easier exercise example is Figure 1711. This front view of the nucleus has some baryons numbered, some not numbered. For this example use the base triangle of 230 231 234: enter those numbers to see this screen. Menu item 3:
3
What is the first baryon number in an equilateral triangle?
28
What is the second baryon number?
31
What is the third baryon number?
32
The middle point will be calculated for an equilateral triangle.
The 20 precision corrections for trans-lanthanide element positions of baryons
3.41 3.414213562
4.82 4.8284271247 ...
two candidates
The top is 0 11.4283637147007 -0.666280163677192
The bot is 0 13.3134352000993 1.99961349701053
____________Candidate is vertex e. Next, find bottom length e-a
e-a length on EQ triangle
1.9996135156872
________________Candidate is vertex d. Next, find top d-b length
d-b length on EQ triangle for tetrahedron
1.99961351568793
------------Candidate is vertex e. Next, find bottom length e-c
e-c length on EQ triangle
1.99961351568793
------------Candidate is vertex d. Next, find top d-a length
d-a length on EQ triangle for tetrahedron
1.9996135156872
-------------Candidate is vertex e. Next, find bottom length e-b
e-b length on EQ triangle
1.99961351568793
--------------Candidate is vertex d. Next, find top d-c length
d-c length on EQ triangle for tetrahedron
1.99961351568793
How to use the program, Part 3: Understanding the terminal screen output
Before we get into interpretations, there is a link for the .xyz file format of the proton and neutron coordinates. That is not for this Perl program, but some errors are corrected in the Table format file that you already found, above. In addition, a PowerPoint file is there for the essay, "Quarks are on academic probation". Also see the Alan Folmsbee youtube channel. The files are at the bottom of this page:
Thoroughness (pyramidalcube.blogspot.com)
Understanding the terminal screen output
The output above is good because the last 6 numbers are nearly 2.000. Since the 6 numbers are about 1.99961351568793. The result are good, and two candidates are shown for the new baryon, above the diagnostic 6 numbers:
The top is 0 11.4283637147007 -0.666280163677192
The bot is 0 13.3134352000993 1.99961349701053
pick one xyz coordinate that is best.
How to use the program, Part 4: Understanding the terminal screen bad output
The end of the output might show six numbers like 2.21073361149907 .
If the six numbers are not almost 2.0000000000, some work needs to be done by editing the Perl source code. The radius of a proton or neutron is 1 unit of length. (That represents the true radius from 0.84fm to 0.95fm, depending on charge radius or gravitational radius). The distance between the centers of two neutrons is 2 units, if they are touching. Maybe the three base baryons have an error in their coordinates. Fix that. But if the coordinates are good, a scaling factor can be edited to make those six numbers all equal 2 units, plus or minus 0.01% .
# Xenon used menu item 4 and scale4: square, not rectangle scaling
# my $scale4 = 0.7071067811;
# my $scale4 = 1.0;
my $scale4 = 0.82;
my $scale3 = 0.61255000; # for menu item 3
# my $scale6 = 0.866;
my $scale6 = 0.701;
Edit $scale3 to be 1.000 or some number to make all 6 results equal to 2.000000 units of length. Run the program again.
How to use the program, Part 5: Neptunium-237 example, based on U-234 data file.
This example uses menu item 4, for a rectangular base. Figure 1716 shows neutron 236 was added onto U-234. The base for that neutron is four baryons numbered 151, 155, 171, 172. Although not all numbers are given, the author uses the Plan Figure 1707 for the right pyramid-2 to get numbers 171 172.
Start the Perl program U234-84.pl. Use menu choice 4, for a rectangular base...
two candidates
The top is 2.00515731798229 2.00088485587416 0
The bot is 4.40905624401771 2.41332870612584 0
____________Candidate is vertex e. Next, find bottom length e-a
e-a length on EQ triangle
1.99679993828213
________________Candidate is vertex d. Next, find top d-b length
d-b length on EQ triangle for tetrahedron
1.99679993828213
------------Candidate is vertex e. Next, find bottom length e-c
e-c length on EQ triangle
1.99679993828213
------------Candidate is vertex d. Next, find top d-a length
d-a length on EQ triangle for tetrahedron
1.99679993828213
-------------Candidate is vertex e. Next, find bottom length e-b
e-b length on EQ triangle
1.99679993828213
--------------Candidate is vertex d. Next, find top d-c length
d-c length on EQ triangle for tetrahedron
1.99679993828213
The result is good, because all 6 numbers at the end are about 2.000. Chose one candidate from the two candidates above the diagnostic 6 numbers.
How to use the program, Part 6: successive approximation
As an example, use Neptunium again, with neutron 236 again. But edit the "goal" file (goals-226.txt)to start a little off from the correct answer. See Table 93 of data set 17, Charge distributions on the nuclei: data set 17, U to Fm. Here is the correct answer from Part 5 of this tutorial:
The bot is 4.40905624401771 2.41332870612584 0
Edit the goal file to be a little wrong, only 3 points on the rectangular base are needed 171 155 151:
171 155 151 4.38 2.2 0
comments below the first line...
Np-227 goal find coordinates for Neutron #226
Edit the Perl source file to use the correct file names for the goals and outputs (without file types):
my $ascii_file3 = "goals-226";
my $ascii_file4 = "output-226"; # output file
Run the program using menu option 6 : UScan234-84.pl . The results are given and a new type of menu appears:
tape out 1 baryon base triangle data for 171 155 151
c-a b-a and c-b lengths of triangle
3.16227765983467 2.44948974235237 2
e-a, e-c, and e-b lengths to bottom.
2.26850369273455 2.11646930948387 1.77372743245793
d-a, d-b, and d-c lengths to top
2.26850369273455 1.77372743245793 2.11646930948387
iiT jjT and kkT are 4.54406057345961 2.85070534937402 0.333333333333333
iiB jjB and kkB are 1.73208180120706 2.36824606662598 0.333333333333333
Welcome to the Interactive Software Menu from Eu153Arti-66g.pl for 234 baryons.
Select a letter to adjust the position of a baryon added onto a base of three baryons.
mxyz ed aircf ls are the menu choices.
m show menu,x adjust x parameters, y and z
o omni mode controls xyz again, instead of 1 dimension
e calculate for bottom candidate
d calculate for top candidate
a again do the loop
i increase, r reduce in a dimension x or y or z to move the candidate
c f coarser or finer stepSize
l load file of base numbers and candidate baryon coordinate guess
s save file of state
use a much finer step size and start with y increases: f f y a
f
using finer stepSize
f
using finer stepSize
y
The y dimension parameters will be changeable
a
tape out 2 baryon base triangle variables for 171 155 151
e-a, e-c, and e-b bottom. 1.79500291363126 1.98864175099275 1.98864175099275 iiB jjB and kkB are 4.25 2.2343 0
Comments
Post a Comment