James Jehiel Ramsey @ 04:37 PM on 27 Jun, 2012

I’ve been trying to test the parameters from Jancu et al., PRB 57, 6493 (1998) for use in NEMO5, and I’m running into a problem. The results are not agreeing with results of previously written Matlab/Octave scripts that determine band edges in bulk GaAs as a function of strain, and the Matlab/Octave results agree with published results. (Compare the Octave results in comp_NEMO_octave-Jancu-biax.eps in the attached archive with the dotted lines in Fig. 1 of the paper by Boykin et al. in PRB 66, 125207 (2002).) With the GaAs parameters from Boykin et al., however, I do get agreement between NEMO5 and the scripts.

I’ve checked the obvious stuff like typos in the parameters, and so far, the problem doesn’t seem to be that simple. Another possible issue is that parameters from Jancu et al. aren’t meant to be used with the scheme from Boykin et al. for strain-dependent shifts of diagonal parameters, but I thought that simply having zero values for the strain_constant_C_* parameters should make these shifts equal to zero. The strain_exponent_V_* parameters are of course set to the values from Jancu et al.

Obviously, I’m missing something, but I’m not sure what.

bulkTests.zip