Consider a date sequence
\begin{align*}
0 \leq t_0 \leq T_s < T_e < T_p,
\end{align*}
where $t_0$ is the valuation date, $T_s$ is the Libor start date, $T_e$ is the Libor end date, and $T_p$ is the payment date. Let $\Delta_s^e = T_e-T_s$.
For $0\le t \le T_s$, define
\begin{align*}
L^e(t, T_s, T_e) = \frac{1}{\Delta_s^e}\bigg(\frac{P(t, T_s)}{P(t, T_e)}-1 \bigg),
\end{align*}
where $P(t, \mu)$ is the price at time $t$ of a zero-coupon bond with maturity $\mu$ and unit face value. Moreover, let $\Delta_e^p = T_p-T_e$.
For $0\le t \le T_e$, define
\begin{align*}
L^p(t, T_e, T_p) = \frac{1}{\Delta_e^p}\bigg(\frac{P(t, T_e)}{P(t, T_p)}-1 \bigg).
\end{align*}
Furthermore, let $Q^{T_e}$ denote the $T_e$-forward measure, and $Q^{T_p}$ denote the $T_p$-forward measure. We assume that, under $Q^{T_p}$, for $0\le t \le T_s$,
\begin{align*}
dL^e &= \mu(t) dt + \sigma_e L^e dW_t^e, \tag{1}\\
dL^p &= \sigma_p L^P\left(\rho dW_t^e + \sqrt{1-\rho^2} dW_t^p \right),
\end{align*}
where, $\sigma_e$, $\sigma_p$, and $\rho$ are constant, while $\{W_t^e, t \ge 0\}$ and $\{W_t^p, t \ge 0\}$ are two independent standard Brownian motions.
Let $E^{T_p}$ be the expectation operator under the $T_p$-forward measure $Q^{T_p}$.
We seek the value defined by
\begin{align*}
P(t_0, T_p)E^{T_p}\big(L^e(T_s, T_s, T_e) \mid \mathcal{F}_{t_0}\big).
\end{align*}
Note that, for $0 \le t \le T_e$,
\begin{align*}
\eta_t \equiv \frac{dQ^{T_e}}{dQ^{T_p}}\big|_{t} &= \frac{P(t, T_e)P(t_0, T_p)}{P(t_0, T_e)P(t, T_p)}\\
&= \frac{1+\Delta_e^pL^p(t, T_e, T_p) }{1+\Delta_e^pL^p(t_0, T_e, T_p) }.
\end{align*}
Then
\begin{align*}
d\eta_t &= \frac{\Delta_e^pL^p(t, T_e, T_p) }{1+\Delta_e^p L^p(t_0, T_e, T_p) }\sigma_p \left(\rho dW_t^e + \sqrt{1-\rho^2} dW_t^p \right)\\
&=\frac{\sigma_p\Delta_e^pL^p(t, T_e, T_p) }{1+\Delta_e^p L^p(t, T_e, T_p) }\eta_t\left(\rho dW_t^e + \sqrt{1-\rho^2} dW_t^p \right).
\end{align*}
Consequently,
\begin{align*}
\hat{W}_t^e = W_t^e - \int_0^t \frac{\rho\sigma_p\Delta_e^pL^p(s, T_e, T_p) }{1+\Delta_e^p L^p(s, T_e, T_p) }ds
\end{align*}
is a standard Brownian motion under the $T_e$-forward measure $Q^{T_e}$. Therefore,
\begin{align*}
\mu(t) &= -\frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t, T_e, T_p)}{1+\Delta_e^p L^p(t, T_e, T_p)}L^e(t, T_s, T_e)\\
&\approx -\frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)}{1+\Delta_e^p L^p(t_0, T_e, T_p)}L^e(t, T_s, T_e)
\end{align*}
in $(1)$ above, given that $L^e$ is a martingale under $T_e$-forward measure $Q^{T_e}$. That is, under the $T_p$-forward measure,
\begin{align*}
dL^e(t, T_s, T_e) &\approx L^e(t, T_s, T_e)\left(-\frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)}{1+\Delta_e^p L^p(t_0, T_e, T_p)} dt + \sigma_e dW_t^e\right),
\end{align*}
and
\begin{align*}
L^e(T_s, T_s, T_e)
&\approx L^e(t_0, T_s, T_e)\exp\bigg(-\frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)}{1+\Delta_e^p L^p(t_0, T_e, T_p)}(T_s-t_0) \\
&\qquad\qquad\qquad\qquad\qquad -\frac{\sigma^2}{2}(T_s-t_0) + \sigma \big(W_{T_s}^e -W_{t_0}^e\big) \bigg)
\end{align*}
Moreover,
\begin{align*}
&\ P(t_0, T_p)E^{T_p}\big(L^e(T_s, T_s, T_e) \mid \mathcal{F}_{t_0}\big) \\
\approx&\ P(t_0, T_p)L^e(t_0, T_s, T_e) \exp\left(-\frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)}{1+\Delta_e^p L^p(t_0, T_e, T_p)}(T_s-t_0) \right)\\
\approx&\ P(t_0, T_p)L^e(t_0, T_s, T_e)\left(1- \frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)}{1+\Delta_e^p L^p(t_0, T_e, T_p)}(T_s-t_0)\right).
\end{align*}
Here, the term
\begin{align*}
- \frac{\rho\sigma_e\sigma_p\Delta_e^pL^p(t_0, T_e, T_p)L^e(t_0, T_s, T_e)}{1+\Delta_e^p L^p(t_0, T_e, T_p)}(T_s-t_0)
\end{align*}
can be treated as the convexity adjustment.