<div dir="ltr">Hello all,<div><br></div><div>I am trying to getting phonon frequencies of biaxially strained silicon and germanium. Particularly, the lattice constants in the three dimensions are as such ax = ay != az in either case(!= refers to "not equal to"). </div><div><br></div><div>To do so I chose ibrav=10 which is an orthorhombic fcc structure, and I explicitly define A,B,C,cosAB,cosAC,cosBC. where A=B!=C and all cosines are 0 (simply because the angles remain at 90 degrees)</div><div><br></div><div>Now everything is OK, until I want to sample the 1BZ. For sampling in the 3 Cartesian directions, ie (100), (010), (001), the matdyn code seems to assume that the lattice constant in all three directions is the same (eg it assumes the value of "A"). Which in the case of the (001) direction is a problem, since the true lattice constant in that direction is C.  (C!=A)</div><div><br></div><div>Unfortunately I have not looked through enough of the code to distinguish a way around this without just compensating for this when I define each q-vector that I pass on to matdyn. However, this seems like a hassle in the case of a true orthorhombic structure where ax!=ay!=az. </div><div><br></div><div>Has anyone else experienced this and/or have a better alternative to simply defining scaled q-vectors?</div><div><br></div><div>Best,</div><div>Cameron</div><div><br></div><div><br></div></div>