Is there a closed-form expression for the term of $e^{t(c \hat{X} + d \hat{Y})}$ that is first-order in $d$, where $t$, $c$, and $d$ are scalars and $\hat{X}$ and $\hat{Y}$ are finite-dimensional linear operators? I guess that to calculate this, you would use the Zassenhaus formula to expand the exponential, then collect all the terms proportional to $d^1$. Unfortunately, so far as I can find the Zassenhaus formula is usually presented organized in powers of $t$, and each power of $t$ has a term proportional to $d^1$, so you would need to consider every term in the $t$ expansion just to get the first term in the $d$ expansion.
(A general expression for the Zassenhaus formula reorganized in powers of $d$ would be even better, but that's probably too much to ask for!)
I don't know about using the Zassenhaus formula, but I did obtain the following expression for the terms that are first-order in $d$: $$\sum_{n\geq1}\frac{c^{n-1}t^n}{n!}\sum\limits_{r,s\geq0,\ \ r+s=n-1}\hat{X}^r\hat{Y}\hat{X}^s. $$
Basically, expand $e^{t(c \hat{X} + d \hat{Y})}$ as a power series: $$e^{t(c \hat{X} + d \hat{Y})}=\sum_{n=0}^\infty \frac{t^n}{n!}(c \hat{X} + d \hat{Y})^n. $$ Then, for each term $(c \hat{X} + d \hat{Y})^n$, we get a single $d$ whenever we obtain one copy of $d\hat{Y}$ and $n-1$ copies of $c\hat{X}$. For example, when $n=2$ we get $$cd(\hat{X}\hat{Y}+\hat{Y}\hat{X}),$$ and for $n=3$ we get $$ c^2d(\hat{X}\hat{X}\hat{Y}+\hat{X}\hat{Y}\hat{X}+\hat{Y}\hat{X}\hat{X}). $$ Summing over all $n$, we obtain the result.