  The goal of this paper is to study the existence of peak solutions for the
following fractional Schr\"{o}dinger-Poisson system:
  \begin{eqnarray*}
  \left\{ \arraycolsep=1.5pt
  \begin{array}{ll}
  \varepsilon^{2s}(-\Delta)^{s}u+u+\phi u=u^p,\ \ \ &\ \mbox{in}\
\Omega,\\[2mm]
  (-\Delta)^{s}\phi=u^2,\ \ \ &\ \mbox{in}\ \Omega,\\[2mm]
  u=\phi=0,\ \ \ \ &\ \mbox{in}\ \mathbb{R}^N\setminus \Omega,
  \end{array}
  \right.
  \end{eqnarray*} where $s\in(0,1)$, $N>2s$, $p\in (1,\frac{N+2s}{N-2s})$,
$\Omega$ is a bounded domain in $\mathbb{R}^N$ with Lipschitz boundary, and
$(-\Delta)^{s}$ is the fractional Laplacian operator, $\varepsilon$ is a small
positive parameter. By using the Lyapunov-Schmidt reduction method, we
construct a single peak solution $(u_\varepsilon,\phi_\varepsilon)$ such that
the peak of $u_\varepsilon$ is in the domain but near the boundary. In order to
characterize the boundary concentration of solutions, which concentrates at an
approximate distance $\varepsilon^{2/3}$ away from the boundary
$\partial\Omega$ as $\varepsilon$ tends to 0, some new estimates and analytic
technique are used.
