Dear @vincentm
In RESCU+, we use something likerlx = Relax.from_totalenergy(calc, fixatoms=[1,2,3,4,5])
to fix certain atoms. Could we fix certain axis as well? For example, if I only want to fix the c direction for [1,2,3,4,5], what should I do?
Thank you so much for your help!