Werden die Werte für `ytick` explizit angegeben, dann werden auch die `yminorticks` ebenfalls angezeigt:
angezeigt. Eine Höhenangabe ist hier nicht erforderlich.
    \documentclass{standalone}
    \usepackage{pgfplots}
    
    \pgfplotsset{compat=newest}
    \usepgfplotslibrary{dateplot}
    
    \begin{document}
    
    \begin{tikzpicture}
      \begin{semilogyaxis}[
          clip=false,
          width = \textwidth,
          axis on top = true,
          xmin = 2016-01-01,
          xmax = 2021-12-31,
          xtick = {2016-01-01, 2017-01-01, 2018-01-01, 2019-01-01, 2020-01-01, 2021-01-01},
          xticklabel style ={/pgf/number format/1000 sep=,
                             xshift = 25pt},
          ymin = 1E-14,
          ymax = 1,
          ylabel = {Wahrscheinlichkeit},
          yminorticks = true,
          ytick = {1E-13, 1E-12, 1E-11, 1E-10, 1E-09, 1E-08, 1E-07, 1E-06, 1E-05, 1E-04, 1E-03, 1E-02, 1E-01, 1}, % <=== Siehe hier
          date coordinates in = x,
          xticklabel=\year,
          date ZERO = 2016-01-01
      ]
      \addplot[blue, only marks, mark size = 0.5] coordinates {
        (2016-03-04, 8.27e-11)
        (2017-03-03, 3.30e-10)
        (2018-01-15, 2.22e-10)
        (2019-01-11, 2.47e-10)
        (2020-01-13, 3.69e-10)
        (2021-01-26, 2.55e-10)
      };
    
      \end{semilogyaxis}
      \end{tikzpicture}
    
    \end{document}