Предположим, что момент инерции диполя конечный.
Я только хотел обратить Ваше внимание на этот вопрос. Может быть, Вам больше подойдёт модель, в которой диполь всегда ориентирован по электрическому полю.
Вероятно можно было бы ее вывести из действия для системы со связями, но тогда встает вопрос о слагаемом, отвечающим за взаимодействие с полем
Да, можно рассмотреть систему со связями, более того, это почти автоматически решает вопрос о слагаемых взаимодействия. Рассмотрим электрический диполь во внешнем ЭМ поле.
Сначала предположим, что диполь образован двумя точечными зарядами

и

(где

). Вектор от отрицательного к положительному заряду обозначим

. Радиус-вектор центра диполя

. Координаты зарядов

. Предполагать, что и масса сосредоточена в этих же точках, необязательно. Имеем голономную связь

. Слагаемые в лагранжиане, отвечающие за взаимодействие с внешним полем:

Пусть

мало. Разложим потенциалы в ряд в окрестности

до первой степени

. Ниже все полевые величины берутся в точке

. Используем тензорные обозначения для декартовых компонент. Индексы пробегают значения

.
![$\begin{array}{lclcl}L_{\text{int}}&=&q (A_i+\frac 12 A_{i,k}a_k)(v_i+\frac 12{\dot a}_i) &-&q(\varphi+\frac 12\varphi_{,i}a_i)\\[0.5ex]& -&q (A_i-\frac 12 A_{i,k}a_k)(v_i-\frac 12{\dot a}_i) &+&q(\varphi-\frac 12 \varphi_{,i}a_i)\\[0.5ex]&=&q (A_{i,k}a_k v_i+A_i {\dot a}_i -\varphi_{,i}a_i)\end{array}$ $\begin{array}{lclcl}L_{\text{int}}&=&q (A_i+\frac 12 A_{i,k}a_k)(v_i+\frac 12{\dot a}_i) &-&q(\varphi+\frac 12\varphi_{,i}a_i)\\[0.5ex]& -&q (A_i-\frac 12 A_{i,k}a_k)(v_i-\frac 12{\dot a}_i) &+&q(\varphi-\frac 12 \varphi_{,i}a_i)\\[0.5ex]&=&q (A_{i,k}a_k v_i+A_i {\dot a}_i -\varphi_{,i}a_i)\end{array}$](https://dxdy-01.korotkov.co.uk/f/0/1/a/01a994402bd5631edc96b92729c520b382.png)

— это декартовы компоненты трёхмерного вектора

, а не пространственные ковариантные компоненты 4-потенциала.
Введём вектор дипольного момента

, где

— единичный вектор, задающий ориентацию диполя.

Теперь можно устремить

так, чтобы

осталось конечным, на формулу это уже не повлияет.
(Оффтоп)
Я долго искал для дипольного момента подходящее обозначение, ведь буквы

и

уже заняты.
Чтобы лагранжиан приобрел благообразный вид, вычтем из него вот такую полную производную по времени (что не повлияет на вариации действия и уравнения движения):

Получим
![$\begin{array}{lcl}L_{\text{int}}&=&\mu\left((A_{i,k}-A_{k,i})v_i n_k+(-\varphi_{,k} - \frac{\partial A_k}{\partial t})n_k\right)\\[0.5ex]&=&\mu(\varepsilon_{ki\ell}B_\ell v_i n_k+E_k n_k)\;\;=\;\;\mu\mathbf n\cdot(\mathbf E+\mathbf v\times\mathbf B)\end{array}$ $\begin{array}{lcl}L_{\text{int}}&=&\mu\left((A_{i,k}-A_{k,i})v_i n_k+(-\varphi_{,k} - \frac{\partial A_k}{\partial t})n_k\right)\\[0.5ex]&=&\mu(\varepsilon_{ki\ell}B_\ell v_i n_k+E_k n_k)\;\;=\;\;\mu\mathbf n\cdot(\mathbf E+\mathbf v\times\mathbf B)\end{array}$](https://dxdy-01.korotkov.co.uk/f/8/8/c/88c2cad955cac5e6d94568111ae14dbb82.png)
Интересный результат. Скалярное произведение дипольного момента на силу Лоренца (действующую на единичный заряд).