Example: TD-CDFT for MoS2 Monolayer (NewResponse)¶

Download NewResp_2DMoS2Restart.run

This example demonstrates how to calculate the frequency-dependent dielectric function with the help of the NewResponse implementation for a two-dimensional system. (see NewResponse) Furthermore, the general setup to run the TD-CDFT section as a restart calculation is presented as well. This allows for splitting of the frequency range into several parts, which can then be calculated in separate calculation without the overhead of evaluating the groundstate properties for each of them! Hence, it is a trivial parallelization possibility.

Groundstate Optimization

$ADFBIN/band <<EOF
  GeometryFile MoS2_2D_1L.xyz

  Symmetry NoSym

  NumericalQuality basic
  Relativistic ZORA

  KSpace
    Grid 5 5
  End

  BasisDefaults
    BasisType DZP
    Core      Large
  End

  End input
EOF
mv RUNKF restart.runkf

TD-CDFT+Restart Calculation

$ADFBIN/band <<EOF
  GeometryFile MoS2_2D_1L.xyz

  Symmetry NoSym

  NumericalQuality basic
  Relativistic ZORA

  KSpace
    Grid 5 5
  End

  BasisDefaults
    BasisType DZP
    Core      Large
  End

  Restart
    SCF
    File restart.runkf
  End

  NewResponse
    nFreq        11
    FreqLow      2.0
    FreqHigh     3.0
    ActiveESpace 10.0
  End

  NewResponseSCF
    nCycle       30
  End

  End input
EOF
mv RUNKF prop.runkf

The results are accessible via the standard output or via the prop.runkf file. For the latter, one can use the ADFreport command $ADFBIN/adfreport prop.runkf RESPDIELRE and $ADFBIN/adfreport prop.runkf RESPDIELIM to print the components of the dielectric function for the real (RESPDIELRE) and imaginary (RESPDIELIM) part separately. In the following tables, only the diagonal components are presented:

The more convenient option is to plot the spectral data directly with the help of ADFspectra. Just type:

$ADFBIN/adfspectra prop.runkf