(I have looked at other questions relating to Doppler Effect with acceleration but none seem to provide the formula in question.) Edit: I wish to avoid using the Einsteinian relativity model. So consider the pulses as sound pulses moving in a body of water at rest in the IRF of the source.
EXPERIMENT
A fixed source (at position $x=0$) emits brief spike pulses at a regular time interval between pulses with period $P_e$. The pulses travel away from the source at constant radial speed $s$ (positive). A distant receiver (at large $x+$) accelerates directly towards the source at constant rate of acceleration $a$ (negative). The receiver velocity $V$ (initially zero at some undefined time in the past) is always directed towards the source (the sign of $V$ is negative). The time interval between received pulses is $P$. Assume that the receiver remains distant from the source during the experiment (i.e. it doesn't arrive at or go past the source). Also the receiver velocity $V$ remains very small compared to $s$. Clearly the value of $P$ will decrease as time passes. The rate of change of $P$ with time is $\dot{P}$.
QUESTION
What is the formula for the rate of change of period in terms of $P_e, a, s, V$?
MY 1st ATTEMPT (shown to be incorrect by Sofia's answer).
The (wavelength) distance between pulses is $\lambda=P_e s$. The period between pulse receptions is $P=\lambda/(s-V)=P_es/(s-V)$. Consider the periods $P_1$,$P_2$ at times $t_1,t_2$ where $P_2=P_1+\Delta P$ and $t_2=t_1+\Delta t$. So $$ \Delta P = P_2-P_1 = P_e s \left( \frac{1}{s-V_2}-\frac{1}{s-V_1} \right) = P_e s \left( \frac{(s-V_1)- (s-V_2)}{s^2-sV_1-sV_2+V_1V_2} \right) $$
$$ = P_e s \left( \frac{V_2-V_1}{s^2-sV_1-s(V_1 + a\Delta t)+V_1(V_1 + a\Delta t)} \right) $$
$$ = P_e s \left( \frac{a\Delta t}{s^2 -2sV_1 - s a\Delta t + V_1^2 +V_1 a\Delta t} \right) $$
$$ = \frac{P_e s}{s^2} \left( \frac{a\Delta t}{1 -2V_1/s - a\Delta t/s + V_1^2/s^2 +V_1 a\Delta t/s^2} \right) $$ Thus $$ \frac{\Delta P}{\Delta t}= \frac{P_e }{s} \left( \frac{a}{1 -2V_1/s - a\Delta t/s + V_1^2/s^2 +V_1 a\Delta t/s^2} \right) $$
and as $\Delta t \rightarrow 0$ we get
$$ \dot P =\frac{\mathrm{d} P}{\mathrm{d} t}= \frac{P_e }{s} \left( \frac{a}{1 -2V/s + V^2/s^2 } \right). $$
Using $P = P_e s/(s-V) \longrightarrow P(s-V)/s = P_e$ we can express $\dot P$ in terms of $P$ rather than $P_e$ as follows
$$ \dot P = \frac { P(s-V) }{s^2} \left( \frac{a}{1 -2V/s + V^2/s^2 } \right) $$
and if $V << s$ then
$$ \dot P \approx \frac { P a }{s} . $$
(Thanks to Sofia for pointing out the error of this 1st solution).
My 2nd Attempt is contained in a separate Answer on this page.