// 拟合结果位于 画图/拉曼结果拟合/250923
#page(flipped: true)[#figure({
  set text(size: 9pt);
  set par(justify: false);
  let m(n, content) = table.cell(colspan: n, content);
  let m2(content) = table.cell(colspan: 2, content);
  let m3(content) = table.cell(colspan: 3, content);
  let pm = sym.plus.minus;
  let m21(content) = table.cell(rowspan: 2, colspan: 1, content);
  let m22(content) = table.cell(rowspan: 2, colspan: 2, content);
  table(columns: 12, align: center + horizon,
    m2[*Notation of Mode*],
      [E#sub[2]-1], [E#sub[2]-2], [E#sub[1]-1], m2[A#sub[1]-1],
      [E#sub[1]-2], [E#sub[2]-3], [E#sub[2]-4], m2[A#sub[1]-2],
    table.cell(rowspan: 3)[*Frequency*],
      [Simulation (cm#super[-1])],
        // E2     E2        E1        A1          E1        E2        E2        A1
        [190.51], [197.84], [257.35], m2[591.90], [746.91], [756.25], [764.33], m2[812.87],
      [Experiment (cm#super[-1])],
        // E2    E2       E1       A1         E1           E2     E2           A1
        [195.5], [203.3], [269.7], m2[609.5], [Invisible], [776], [Invisible], m2[839],
      //           E2     E2     E1     A1       E1   E2     E2   A1
      [Error (%)], [2.6], [2.7], [4.6], m2[2.9], [-], [2.5], [-], m2[3.1],
    table.cell(rowspan: 4)[*Raman Tensor*],
      [Non-zero components],
        // E2          E2             E1        A1
        [xy, xx, -yy], [xy, xx, -yy], [xz, yz], [xx, yy], [zz],
        // E1     E2             E2             A1
        [xz, yz], [xy, xx, -yy], [xy, xx, -yy], [xx, yy], [zz],
      [Analysis result (a.u.)],
        // E2                          E2                      E1
        [$2epsilon_2-2zeta_2-4eta_2$], [$2epsilon_2-2zeta_2$], [$-2epsilon_1-2zeta_1$],
        // A1                             E1                       E2
        [$-2epsilon_5$], [$-2epsilon_6$], [$-2epsilon_1+2zeta_1$], [$bold(8a_2)+2epsilon_2-2zeta_2-4eta_2$],
        // E2                    A1
        [$-2epsilon_2+2zeta_2$], [$-2zeta_5$], [$-2zeta_6$],
      [Calculation result (a.u.)],
        // E2   E2      E1      A1              E1      E2       E2      A1
        [0.17], [1.13], [2.43], [2.83], [1.79], [0.09], [88.54], [0.50], [0.01], [1.78],
      [Experiment result #linebreak() (a.u.)],
        // E2   E2      E1      A1              E1           E2       E2           A1
        [0.15], [1.11], [0.68], [1.64], [3.15], [ ], [88.54], [Invisible], [Invisible], [2.97],
    table.cell(rowspan: 2)[*Raman Shift of Different Incidence Directions*],
      [Calculation result (cm#super[-1])],
        // E2    E2        E1      A1          E1      E2       E2      A1
        [0.263], [-0.221], m21[-], m2[-0.188], m21[-], [0.004], m21[-], m22[-],
      [Experiment result (cm#super[-1], average #pm standard deviation)],
        // E2              E2                  A1                    E2
        [0.228 #pm 0.051], [-0.260 #pm 0.025], m2[-0.279 #pm 0.016], [0 #linebreak() (used for calibration)],
  )},
  caption: [
    Negaligible-polarized Phonons at $Gamma$ Point.
    The calculated phonon frequencies had a slight underestimation of 2-5% comparing to experimental values,
    which might be attributed to the known tendency of the PBE functional underestimating interatomic forces (cite).
  ],
)<table-nopol>]

