Example: Output Distributions

RustBCA includes a feature to output reflected ion and sputtered atom energy-angle distributions, as well as implanted x-, y-, and z-distributions. To enable this feature, build with --features distributions. Optionally, build with --features distributions,no_list_output to also disable the particle list output.

In order to use this, you must add the following to the [options] block of your input file:

energy_min = (float)
energy_max = (float)
energy_num = (int)
angle_min = (float)
angle_max = (float)
angle_num = (int)
x_min = (float)
y_min = (float)
z_min = (float)
x_max = (float)
y_max = (float)
z_max = (float)
x_num = (int)
y_num = (int)
z_num = (int)

Running cargo run --release --features distributions 1D examples/layered_geometry_1D.toml produces the following:

[energies]
v = 1
dim = [20]
data = [0.0, 105.26315789473684, 210.52631578947367, 315.7894736842105, 421.05263157894734, 526.3157894736842, 631.578947368421, 736.8421052631578, 842.1052631578947, 947.3684210526316, 1052.6315789473683, 1157.8947368421052, 1263.157894736842, 1368.421052631579, 1473.6842105263156, 1578.9473684210525, 1684.2105263157894, 1789.4736842105262, 1894.7368421052631, 1999.9999999999998]

[angles]
v = 1
dim = [20]
data = [0.0, 4.7368421052631575, 9.473684210526315, 14.210526315789473, 18.94736842105263, 23.684210526315788, 28.421052631578945, 33.1578947368421, 37.89473684210526, 42.63157894736842, 47.368421052631575, 52.10526315789473, 56.84210526315789, 61.57894736842105, 66.3157894736842, 71.05263157894737, 75.78947368421052, 80.52631578947367, 85.26315789473684, 90.0]

[x_range]
v = 1
dim = [20]
data = [0.0, 0.005263157894736842, 0.010526315789473684, 0.015789473684210527, 0.021052631578947368, 0.02631578947368421, 0.031578947368421054, 0.03684210526315789, 0.042105263157894736, 0.04736842105263158, 0.05263157894736842, 0.05789473684210526, 0.06315789473684211, 0.06842105263157895, 0.07368421052631578, 0.07894736842105263, 0.08421052631578947, 0.08947368421052632, 0.09473684210526316, 0.09999999999999999]

[y_range]
v = 1
dim = [20]
data = [-0.1, -0.08947368421052632, -0.07894736842105264, -0.06842105263157895, -0.05789473684210527, -0.04736842105263159, -0.0368421052631579, -0.026315789473684223, -0.015789473684210534, -0.005263157894736845, 0.005263157894736831, 0.01578947368421052, 0.02631578947368421, 0.0368421052631579, 0.04736842105263156, 0.05789473684210525, 0.06842105263157894, 0.07894736842105263, 0.08947368421052632, 0.09999999999999998]

[z_range]
v = 1
dim = [20]
data = [-0.1, -0.08947368421052632, -0.07894736842105264, -0.06842105263157895, -0.05789473684210527, -0.04736842105263159, -0.0368421052631579, -0.026315789473684223, -0.015789473684210534, -0.005263157894736845, 0.005263157894736831, 0.01578947368421052, 0.02631578947368421, 0.0368421052631579, 0.04736842105263156, 0.05789473684210525, 0.06842105263157894, 0.07894736842105263, 0.08947368421052632, 0.09999999999999998]

[reflected_ead]
v = 1
dim = [20, 20]
data = [11, 25, 43, 45, 56, 69, 83, 72, 74, 83, 67, 49, 44, 22, 7, 7, 0, 0, 0, 0, 9, 23, 41, 50, 68, 65, 83, 94, 69, 80, 73, 78, 61, 39, 23, 9, 3, 0, 0, 0, 3, 27, 35, 52, 69, 84, 87, 82, 92, 82, 86, 70, 72, 47, 24, 12, 5, 0, 0, 0, 13, 28, 43, 54, 73, 103, 95, 93, 93, 95, 106, 88, 70, 56, 50, 24, 8, 2, 0, 0, 9, 24, 45, 73, 77, 105, 104, 99, 106, 91, 107, 85, 78, 52, 41, 31, 6, 3, 0, 0, 11, 40, 48, 59, 103, 104, 100, 96, 100, 116, 98, 85, 74, 57, 37, 44, 15, 2, 0, 0, 19, 37, 54, 77, 84, 117, 107, 102, 119, 125, 97, 100, 76, 65, 55, 30, 13, 0, 0, 0, 16, 21, 39, 48, 76, 90, 102, 128, 112, 103, 114, 90, 74, 60, 53, 29, 24, 0, 0, 0, 9, 23, 20, 39, 58, 59, 74, 87, 82, 97, 90, 84, 72, 77, 48, 26, 24, 2, 0, 0, 8, 10, 29, 40, 48, 72, 57, 63, 86, 92, 83, 71, 66, 61, 48, 40, 17, 2, 0, 0, 2, 14, 22, 39, 43, 46, 53, 81, 63, 68, 69, 79, 65, 54, 50, 31, 18, 5, 0, 0, 5, 8, 18, 22, 30, 38, 35, 51, 67, 54, 56, 53, 63, 49, 41, 31, 17, 4, 0, 0, 2, 5, 11, 17, 24, 32, 27, 51, 57, 39, 58, 56, 60, 32, 41, 27, 16, 9, 0, 0, 2, 3, 10, 10, 18, 17, 34, 38, 36, 34, 45, 54, 46, 35, 35, 33, 13, 2, 0, 0, 0, 0, 0, 0, 0, 3, 1, 5, 15, 23, 26, 43, 34, 45, 46, 27, 27, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 3, 5, 20, 31, 31, 17, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

[sputtered_ead]
v = 1
dim = [20, 20]
data = [119, 334, 488, 656, 795, 847, 912, 927, 918, 919, 832, 739, 617, 521, 367, 232, 164, 95, 42, 0, 4, 13, 14, 32, 27, 30, 40, 37, 40, 43, 34, 35, 24, 14, 15, 7, 1, 0, 0, 0, 2, 5, 9, 5, 13, 7, 9, 10, 14, 9, 10, 8, 14, 6, 2, 2, 0, 0, 0, 0, 0, 0, 1, 3, 4, 6, 6, 4, 5, 2, 7, 3, 3, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 2, 4, 5, 3, 4, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 3, 2, 0, 1, 3, 1, 0, 3, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

[implanted_x]
v = 1
dim = [20]
data = [9826, 17963, 14187, 14036, 12028, 8905, 5504, 2722, 1097, 409, 112, 26, 5, 2, 0, 0, 0, 0, 0, 0]

[implanted_y]
v = 1
dim = [20]
data = [0, 0, 0, 0, 2, 88, 1268, 7390, 20149, 28901, 20149, 7452, 1324, 98, 1, 0, 0, 0, 0, 0]

[implanted_z]
v = 1
dim = [20]
data = [0, 0, 0, 0, 1, 88, 1275, 7379, 20142, 28948, 20181, 7480, 1248, 79, 1, 0, 0, 0, 0, 0]