Sample Input for FROZSL: NaCl

The following is sample input for NaCl. We calculate the frequency of the transverse optical (TO) phonon mode at k=0. This is at the gamma point in the first Brillouin zone. In the input below, we indicate the gamma point with GM. Click on the start button below. This produces the output of frozsl_init, giving the frozen phonons that need to be calculated. From the output, we find that there is only one irrep involved: GM4-. There are two symmetry vectors, but since one combination is the acoustic mode at zero frequency, we find that we only need to calculate the energy of one frozen phonon.

In equilibrium, the atomic positions are at

 position of each atom (dimensionless)
    0.000000000000000   0.000000000000000   0.000000000000000
    0.500000000000000   0.000000000000000   0.000000000000000
 position of each atom (cartesian)
    0.000000000000000   0.000000000000000   0.000000000000000
    5.300000000000000   0.000000000000000   0.000000000000000
 Wyckoff positions
   Na  a  z =  0.00000  
   Cl  a  z =  0.50000

We have three choices above for expressing the atomic positions. The third choice gives the Wyckoff position in the isotropy subgroup 107 C4v-9 determined by the distortion produced by the frozen phonon. In the frozen phonon, the atomic positions are at

 atomic positions (dimensionless)
    0.018867924528302   0.000000000000000   0.000000000000000
    0.487763791681082   0.000000000000000   0.000000000000000
 atomic positions (cartesian)
    0.200000000000000   0.000000000000000   0.000000000000000
    5.170296191819464   0.000000000000000   0.000000000000000
 atomic positions (Wyckoff parameters in subgroup)
    0.00000   0.00000   0.07915
    0.00000   0.00000   0.44867

The energy of the frozen phonon is the difference in energy of these two structures. Using a first-principles energy calculation, we find the energy of this frozen phonon to be 0.83 millihartrees. This value has been entered into the energy input field below. Select "output from frozsl" below, and then click on the start button. Now we see the frequencies of the modes: the acoustic mode at zero frequency and the optical mode at 170/cm.

Output from frozsl_init
Output from frozsl (enter energies in box below)
Include effective charge tensors (enter tensors in box below)
Iterate, using output of frozsl as input to frozsl_init a second time

Input to frozsl_init:

Energy of each frozen mode (millihartrees)

Dielectric tensor and effective charge tensor for the first atom of each Wyckoff position.