“…ADE recasts the constitutive relations as where in the frequency domain P Metal , P Molecule , and M Molecule are defined as The polarization and magnetization fields are then inverse Fourier transformed, and all of the resulting derivatives are expanded with a central-difference approximation to generate the update equations, defined for chiral polarizations and chiral and achiral magnetizations as Here Δ t is the time step size, n is the current time step, i , j , and k are grid point coordinates in the x , y , and z directions, and E i +1/2, j +1/2, k +1/2 n , H i,j,k n –1/2 , and H i +1/2, j +1/2, k +1/2 n –1/2 are determined through spatially averaging the field components. The update equations for the nonchiral components of P̆ Metal and P̆ Molecule are defined in refs and , respectively. , …”