Professor Gianluca Gorni,

Thank you for your response.

I have also recalculated and here are my findings:

1. I calculated U235NEW revision1 using both ver14 and ver11.2, but it did not converge, so I believe the solution is likely incorrect. (I apologize. My previous response was my mistake.)

2. I calculated U235OLD revision1 using ver11.2, and it converged, so I believe this is likely correct.

3. I calculated U235OLD revision1 using ver14, but I could not obtain the correct solution.

Thank you very much for performing various calculations for me.

Attachments:

U235NEW revision1.nb

U235OLD revision1.nb

Professor Gianluca Gorni,

I apologize.

I was mistaken. The question below was incorrect.

The correct question is as follows:

The attached U235OLD is the result calculated with Ver11.2, and the curve is mountain-shaped. The attached U235NEW is the result calculated with Ver14, and the curve is valley-shaped.

Since the coefficient of x² is 1 and the coefficient of x is approximately 92 (the atomic number of U), both calculations should be correct.

Could you please revise the calculation in the attached U235NEW to make the curve mountain-shaped?

I would appreciate your guidance.

Below is the incorrect question:

I have tried various calculations based on the code you provided.

As a result, when I used the same code to perform calculations and plot graphs in Version 11.2, as shown in the attached "Calculation of U235 (ver11.2)," the graph appeared mountain-shaped.

However, when I used the exact same code to perform calculations and plot graphs in Version 14, as shown in the attached "Calculation of U235 (ver14)," the graph appeared valley-shaped.

I would greatly appreciate it if you could kindly instruct me on how to make the graph appear mountain-shaped in Version 14, as it did in Version 11.2.

Thank you very much for your assistance.

Attachments:

U235NEW.nb

U235OLD.nb

Yes, here is the notebook with my code.

Attachments:

simple M14.nb

I just discovered that FindMinimum works fine after the change of variables k[n,m]->1/j[n,m]:

oldVars = Flatten[z];
newVars = oldVars /. k -> j;
changeOfVars = Thread[oldVars -> 1/newVars];
sol = FindMinimum[Total[w^2] /. changeOfVars, newVars]
sol /. (j[n__] -> v_) :> (k[n] -> 1/v)

We should check if the solution coincides with the one in version 9.

Thank you for your response. The solution for k in Version 9 has no particular significance. I apologize, but please look at the new attachment "U235". This η20 represents the distance (fm) between two fission fragments immediately after nuclear fission. For example, η20[atomic number 28, atomic number 64] indicates that the distance between these two fission fragments is 0.930472 fm. Additionally, this calculation can be performed in Mathematica versions 9 through 11.2. Correct solutions cannot be obtained in later versions.

Note: "fm" refers to femtometer, which is a unit of length commonly used in nuclear physics (1 fm = 10⁻¹⁵ meters).

Attachments:

U - 235.nb

This code performs important nuclear fission physics calculations. Specifically:

Purpose: - Uses nuclear fission yield data from JENDL-5 (Japanese Evaluated Nuclear Data Library version 5) - Calculates initial distance between two fission fragments at the moment of fission (effective fission distance) - Targets actinide nuclides (Np-237, Pu-239, Pu-240, U-235, U-238)

Issue: - Calculations work properly in Mathematica versions 9 to 11.2 - Cannot perform calculations in Mathematica 11.3 and later due to changes in FindMinimum specifications - Same issue exists in current Mathematica 14

Technical Question: How can we make FindMinimum function in the current version?

The goal is to make these important nuclear physics calculations executable in Mathematica 14.

Attachments:

Np - 237.nb

Pu - 239.nb

Pu-240.nb

U - 235.nb

U-238.nb