The calculation is described by “quatermion” which is frequently used in computer graphics. Quatermion provides the simplest means for coordinate translation and reduces the calculation cost. A quatemion is defined as follows. Display FormulaHere, is a point or vector in 3-D coordinates. The conjugation of is described as . Although the real number is arbitrary, is the simplest description. For the translation to be applied, two conditions must be satisfied. The first condition is that the oriented vector in the interaction coordinates must fit that in the simulation coordinates. The second is that the incident point and vector in the interaction coordinates must also fit those in the simulation coordinates. For the first condition, we use Eq. (8). is translated to by the following process: Display Formulawhere is the vector in the direction of orientation. The vectors and are the unit vectors in the - and -directions in pciMC simulation coordinates. After the incident point and the vector in the photon-RBC interaction coordinates have been selected by Eq. (9), they are rotated around the -axis so that becomes . The translated incident vector and the rotation are described by the following equations: Display Formulawhere is randomly selected to be either positive or negative so that the photon azimuthal angle is uniformly distributed. Next, for the second condition, the translated point and vector are additionally rotated so that becomes . The rotation is described as follows: Display FormulaFinally, using these equations, the scattered point and the vector are calculated as follows. Display Formulawhere is the incident point prior to the scattering event, is the scale parameter to fit the scale of the interaction coordinates with those of the simulation coordinates.