I am trying to determine some mechanical properties of Hf-Nb-Ta-Zr high entropy alloy. But I’m not being able to relax it properly & it’s causing absurd values of elastic moduli, double-necking during tensile tests etc. issues. Same codes work for single metals but I don’t know what to do for HEAs. The potential energy values do converge in thermo outputs… I can’t figure out what the exact issue is. The minimization part of the code is:
I have run the same code for for Al and Si; both yielded satisfactory results. Those two only required the direct minimization command. For the alloy, just the necessary create_box and create_atoms command were changed and it caused weird results even after nvt & npt. So the problem is either the potential or my relaxation procedure. Since the potential was procured from the reference of a journal paper, I think my code is the issue here
Yes, but it is more likely your code for creating the system is at fault for issues like you describe than the minimization protocol. This is supported by the fact, that you get expected results for single element systems.
Wow so that means the problem lies before minimization. I mistook the initialization as too simple to have a problem there, sorry. Can you please let me know where the issue is? This is the initial geometry creation upto minimization:
units metal
dimension 3
boundary p p p
atom_style atomic
variable latparam equal 3.63
variable fHf equal 25
variable fNb equal 25
variable fTa equal 25
variable fZr equal 25
----------------------- ATOM DEFINITION ----------------------------
lattice bcc ${latparam}
region HEA block 0 100 0 20 0 20
create_box 4 HEA
------------------ Create different types of Atoms --------------------
create_atoms 1 region HEA
set type 1 type/fraction 2 ((v_fNb+v_fTa+v_fZr)/100.0) 1738
set type 2 type/fraction 3 ((v_fTa+v_fZr)/(v_fNb+v_fTa+v_fZr)) 1755
set type 3 type/fraction 4 $(v_fZr/(v_fTa+v_fZr)) 6688
------------------------ FORCE FIELDS ------------------------------
pair_style eam/alloy
pair_coeff * * Hf-Nb-Ta-Zr.eam.alloy Hf Nb Ta Zr
neighbor 2.0 bin
neigh_modify delay 10 check yes