To allow for the optimisation of cell parameters based on atomic gradients from the embedded calculation, the molecular gradients are transformed into stress components. Since our model is based on the description of a single molecule embedded in an environment, defini-tions of stresses and strains use coordinates relative to the "central molecule". Contributions of all molecules m in the cluster are calculated and added up to