I’ve now added a new module which allows coupling between the magnetic order and lattice strains.
Currently the module allows computations with uniform external stresses, as well as non-uniform strains (elastostatics). This module is now being developed to solve the elastic equation of motion with implementation details to follow in a publication.
As an example, below we have a poly-crystalline magnetic thin film, starting from no applied stress at zero field:
Now we apply an in-plane stress (e.g. as might be obtained from coupling to a piezo-electric substrate). As a result the magnetisation tends to rotate towards one of the crystallite cubic axes: