Magnetoelastic Options to mcphas

-doeps	If the program mcphas is started with option -doeps and it finds elastic constants in the input file mcphas.j (note, that the elastic constants in the input file are normalised to the primitive crystallographic unit cell, units are meV / primitive crystallographic unit cell), it will use these and determine selfconsistently the strain $\epsilon$ by solving equations (136) and (134). Elastic energy and strain tensor are stored in results/mcphas.fum. In this way it is be possible to model Jahn Teller transitions, phase diagrams, magnetostriction, thermal expansion (magnetic part) and dynamics consistently based only on point charges and Born von Karman springs. If mcphas is used with option -doeps and it finds files mcphas.djdx, mcphas.djdy, mcphas.djdz with derivatives of the two ion interaction parameters with respect to spatial coordinates $x,y,z$ respectively, it performs a calculation of the corresponding correlation functions and computes the strain using equation (140) at each iteration. Similar, If mcphas is used with option -doeps and it finds files mcphas.djdeps1, mcphas.djdeps2,..., mcphas.djdeps6 with derivatives of the two ion interaction parameters with respect to strain components in Voigt notation $\epsilon_1,\dots,\epsilon_6$ respectively,it performs a calculation of the corresponding correlation functions and computes the strain using equation (140) at each iteration. Subsequently, making use of the Taylor expansion of the interaction constants (138) the computed strain is used in the next iteration for the the mean field Hamiltonian (mean field theory is used by McPhase to solve the Hamiltonian (137) and (134)). At the end of the iteration loop mean fields, moments $\langle I \rangle$ and a selfconsistent strain tensor $\epsilon$ is obtained.
-linepscf	The option -linepscf together with -doeps will trigger a calculation where equation (136) (and in case of exchange striction terms (140)) is used to calculate the strain $\epsilon$, however the last two terms in (134) are always evaluated for zero strain, i.e. the crystal field striction is not calculated selfconsistently, but only in the linear approximation assuming that the strain leads only to a negligible perturbation of the crystal field Hamiltonian.
-linepsjj	The option -linepsjj together with -doeps will trigger a calculation where equation (140) is used to calculate the strain $\epsilon$, yet the modification of the two ion interaction by the strain is assumed to be small and two ion interaction parameters (138) are used in the mean field loop always for zero strain, i.e. the exchange striction is not calculated selfconsistently, but only in the linear approximation assuming that the strain leads only to a negligible perturbation of the two ion interaction Hamiltonian.

The length change $\Delta L/L$ of a sample in a dilatometer experiment can be calculated from the strain tensor components using ²³

  $$\frac{\Delta L}{L}=\sum_{\alpha\beta} \epsilon_{\alpha\beta} \hat L_{\alpha} \hat L_{\beta}$$ (135)

where $\hat \mathbf L$ denotes the unit vector in the direction of measurement.