I'm not sure if this solution is available in closed form, but after drawing it out I do think there will be two unique solutions always. I unfortunately have no clue where to start.
Given:
- An ellipse with x radius $a$, y radius $b$ and centre point $(0, 0)$
- A point on that ellipse $P_1$
- A distance $d$
Find all points $P_{2j}$ (i.e. find their coordinates) on the ellipse which satisfy the following condition (where $C$ is the chord joining $P_1$ and $P_{2j}$):
- The distance between $C$ and a line which is tangent to the ellipse and parallel to $C$ (on the side of the smaller segment formed by $C$ i.e. the "nearer" side) is equal to $d$
In other words, the maximum distance between $C$, and the nearer side of the ellipse from $C$, must equal $d$.
The diagram I've drawn below is an example of this. $P_{21}$ and $P_{22}$ are the desired points, whereas $P_{23}$ is an example of an invalid point. (The distances are not fully accurate)
The reason I believe the solution is closed, is because as you sweep the point $P_2j$ around the ellipse, the distance (required to be $d$) increases, reaches a turning point, and then decreases.



An afterthought that looks me important (it's why I place it in front of my answer).
In fact your problem deals fundamentally with parallel (or offset) curves of an ellipse at distance $d$, as described here with parametric equations:
$$x=\left(a+\dfrac{bd}{\sqrt{a^2 \sin^2 t+b^2 \cos^2 t}}\right)\cos t, \ \ y=\left(b+\dfrac{ad}{\sqrt{a^2 \sin^2 t+b^2 \cos^2 t}}\right)\sin t$$
to which the two tangents issued from $P_1$ can/should be drawn. Getting the second point(s) $P_2$ (see my figure for notations) looks me accessible for example by using the method giving, for any point $(u,v)$ in the plane the closest point on the ellipse with equation $\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1$.
The big advantage of this point of view is that it gives a global understanding of the question.
Initial answer : Here is a methodology:
$$ u x + v y-p=0 \ \ \text{where} \ \ u^2+v^2=1, \tag{1}$$
(sometimes, one takes $u=\cos \theta, \ v=\sin \theta$ but working with angles in this context isn't necessary)
you can interpret $(u,v)$ as a unit normal vector to line $L$ and (the most important here) $p$ as the (shortest) distance of line $L$ from the origin (= the length of the shortest line segment $OH$ where $H$ belongs to the line). Therefore the ellipse's chord parallel to the tangent at distance $d$ from it has the following equation:
$$ u x + v y-(p\color{red}{-d})=0\tag{1'}$$
$$\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1 \tag{1''}$$
Let us extract from (1'') the following cartesian equation of its upper part is
$$y=f(x)=\dfrac{b}{a}\sqrt{a^2-x^2}\tag{2}$$
with $$f'(x)=\dfrac{-bx}{a}\dfrac{1}{\sqrt{a^2-x^2}}\tag{3}$$
$$y-f(x_0)=f'(x_0)(x-x_0)\tag{4}$$
Caution: in this perspective, $(x_0,y_0)$ aren't thought as the coordinates of $P_1$ (that I will call $(x_1,y_1)$) but the "kind of" midpoint on the future arc $P_1P_2$ where the tangent is to be taken (see figure).
which gives, by reporting (2) and (3) into (4):
$$y-\dfrac{b}{a}\sqrt{a^2-x_0^2}=\dfrac{-bx_0}{a}\dfrac{1}{\sqrt{a^2-x_0^2}}(x-x_0),\tag{4'}$$
$$(cx_0)x+y\sqrt{a^2-x_0^2}-ca^2=0 \ \text{where} \ c:=\frac{b}{a}\tag{4''}$$
then put the resulting equation under the form (1). Caution: obtaining a form $ux+vy+w=0$ isn't enough: you must have a $(u,v)$ with unit norm. If such is not the case, divide the equation by the norm $\sqrt{u^2+v^2}$
Then deduce the form (1') of the chord at distance $d$ from the tangent.
Express that $P_1(x_1,y_1)$ belongs to this chord by writing:
$$u x_1+ v y_1-(p-d)=0\tag{5}$$
This gives an equation constraining $x_0$ to take a certain value (in fact two values in general).
Here is a figure obtained with Geogebra: