--- title: Two-Dimensional Doi–Onsager Model url: https://www.emergentmind.com/topics/two-dimensional-doi-onsager-model type: topic --- # Two-Dimensional Doi–Onsager Model The two-dimensional Doi–Onsager model is a mean-field theory for orientational ordering of rod-like molecules when the orientation is reduced to a single angular variable on the circle. In its planar equilibrium form, the unknown is a probability measure or density on \(S^1\), subject to head–tail symmetry, and the competition between entropic spreading and an interaction kernel determines whether the isotropic state remains stable or nonuniform stationary states emerge. Recent phase-transition analysis for the kernel \(W_{DO}(\theta)=-|\sin(2\pi\theta)|\) establishes that the critical coupling is exactly \(K_c=K_\#=3\pi/4\) and that the transition is continuous, while earlier nonlinear-functional-analytic work on the stationary equation with general even periodic kernels described uniqueness of the isotropic state in a low-parameter regime and the local bifurcation structure of nontrivial branches [2604.16288]. ## 1. Mean-field formulation on the circle In the phase-transition formulation, the state space is the set of probability measures \(q\in \mathcal P(T)\) on the unit circle \[ T=\left[-\frac12,\frac12\right) \] with periodic boundary conditions. The free energy is \[ \mathcal F_K(q)=\int_T \log q\,dq(\theta)-K\iint_{T\times T}W(\theta-\theta')\,dq(\theta)\,dq(\theta'), \] where \(K\ge 0\) is the coupling strength and \(W\) is an even interaction potential. For the two-dimensional Doi–Onsager model, the interaction is \[ W_{DO}(\theta)=-|\sin(2\pi\theta)|. \] The uniform distribution \[ q_u\equiv 1 \] is the isotropic state. At \(K=0\), the free energy reduces to relative entropy with respect to \(q_u\), so the isotropic state is uniquely minimizing. As \(K\) increases, the interaction term can destabilize \(q_u\), and the resulting change in the structure of minimizers is the phase transition [2604.16288]. A closely related stationary formulation writes the unknown as an orientation distribution function \(f(\theta)\ge 0\) on \([0,2\pi)\), normalized by \[ \int_0^{2\pi}f(\theta)\,d\theta=1, \] and constrained by the reversal symmetry \[ f(\theta)=f(\theta+\pi). \] In that formulation, the constant distribution \[ f_0(\theta)=\frac1{2\pi} \] represents the isotropic phase, and equilibrium states are stationary solutions of a nonlinear self-consistency equation [1503.01477]. ## 2. Symmetry, stationary equations, and isotropic equilibrium The head–tail symmetry of rods is encoded by \(\pi\)-periodicity. In the stationary formulation with convolution kernel \(K(\theta-\theta')\), the interaction potential is \[ U(f)(\theta)=\lambda\int_0^{2\pi}K(\theta-\theta')f(\theta')\,d\theta', \] where \(\lambda>0\) is the concentration/inverse temperature parameter, and the kernel satisfies \[ K(\theta+\pi)=K(\theta),\qquad K(-\theta)=K(\theta). \] Stationary solutions satisfy the Doi–Onsager self-consistent equation \[ f(\theta)=\frac{e^{-U(f)(\theta)}}{\int_0^{2\pi}e^{-U(f)(\varphi)}\,d\varphi}, \qquad f(\theta)=f(\theta+\pi), \] so equilibria are fixed points of a nonlinear integral operator [1503.01477]. For functional-analytic treatment, the stationary problem can be rewritten in terms of the centered potential \[ V(\theta)=U(\theta)-\lambda\bar K,\qquad \bar K=\frac1{2\pi}\int_0^{2\pi}K(\theta)\,d\theta, \] which satisfies \[ V(\theta)=\lambda\,\frac{\displaystyle\int_0^{2\pi}\widetilde K(\theta-\theta')e^{-V(\theta')}\,d\theta'}{\displaystyle\int_0^{2\pi}e^{-V(\theta)}\,d\theta}, \qquad \widetilde K(\theta)=K(\theta)-\bar K, \] together with \[ \int_0^{2\pi}V(\theta)\,d\theta=0,\qquad V(\theta)=V(\theta+\pi)=V(2\pi-\theta). \] The natural Hilbert space is \[ H=\left\{V\in H^1([0,2\pi]) : V(\theta)=V(\theta+\pi)\text{ a.e.},\ \int_0^{2\pi}V=0,\ V(\theta)=V(2\pi-\theta)\text{ a.e.}\right\}, \] and the equation becomes the fixed-point problem \(V=\lambda T(V)\) with compact nonlinear operator \[ T(V)(\theta)= \frac{\displaystyle\int_0^{2\pi}\widetilde K(\theta-\theta')e^{-V(\theta')}\,d\theta'} {\displaystyle\int_0^{2\pi}e^{-V(\theta)}\,d\theta}. \] This representation isolates the isotropic state as the trivial solution \(V\equiv 0\) and makes degree-theoretic and bifurcation arguments available [1503.01477]. ## 3. Fourier structure and mode selection The two-dimensional model is especially tractable in Fourier variables. For the phase-transition study, the Doi–Onsager kernel has the explicit expansion \[ W_{DO}(\theta)=-\frac{2}{\pi}+\frac{4}{\pi}\sum_{\ell=1}^\infty \frac{\cos(4\pi \ell \theta)}{4\ell^2-1}. \] This shows that the interaction is \(\frac12\)-periodic, so only the even Fourier modes are active. Its nonzero Fourier coefficients are \[ W_{DO}(2\ell)=\frac{2}{\pi(4\ell^2-1)},\qquad \ell\ge 1, \] and in particular \[ 2W_{DO}(2)=\frac{4}{3\pi}. \] The dominant attractive mode is therefore the \(k=2\) mode, which determines the linear stability threshold of the isotropic state [2604.16288]. In the stationary-equation literature with general \(\pi\)-periodic kernels, one writes \[ K(\theta)=\sum_{m=0}^\infty k_m\cos(2m\theta), \] with coefficients \(k_m\) controlling uniqueness and bifurcation. For the classical Onsager kernel in that normalization, \[ K(\theta)=|\sin\theta|, \qquad k_m=-\frac{4}{4m^2-1},\qquad m\ge 1. \] The corresponding local bifurcation points are \[ \lambda_m=-\frac{2}{k_m}=\frac{4m^2-1}{2},\qquad m=1,2,3,\dots, \] and all bifurcations are supercritical. The first bifurcating branch occurs at \[ \lambda_1=\frac32, \] is stable, and all higher bifurcating branches are unstable [1503.01477]. These two Fourier descriptions are written in different normalizations, but both emphasize the same structural fact: the planar Doi–Onsager problem is mode-selective, and the first loss of isotropic stability is governed by a distinguished low Fourier mode. ## 4. Critical coupling and the continuous phase transition The phase-transition threshold \(K_c\) is defined by the minimization problem for \(\mathcal F_K\): if \(KK_c\), there exists a nonuniform minimizer. The transition is continuous when the unique minimizer at criticality is \(q_u\), equivalently when minimizers for \(K\downarrow K_c\) converge back to the isotropic state. In general, \[ K_c\le K_\#, \] where \(K_\#\) is the linear stability threshold determined by the second variation at \(q_u\) [2604.16288]. For a general interaction, the second variation at the uniform state is \[ \frac{\delta^2 \mathcal F_K}{\delta q^2}(q_u)[\varphi,\varphi] = 2\sum_{k=1}^\infty \bigl(1-2KW(k)\bigr)|\varphi(k)|^2. \] Hence \(q_u\) is linearly stable as long as \[ K\le \frac{1}{2\max_k W(k)}. \] For the Doi–Onsager kernel, the maximal attractive coefficient occurs at \(k=2\), so \[ K_\#=\frac{1}{2W_{DO}(2)}=\frac{3\pi}{4}. \] The main result for the two-dimensional Doi–Onsager model is that the phase transition occurs exactly at this linear threshold and is continuous: \[ K_c=K_\#=\frac{3\pi}{4}. \] At criticality, the isotropic state is still the unique global minimizer, so nonuniform minimizers do not pre-exist below the stability threshold [2604.16288]. In the stationary formulation with parameter \(\lambda\), a related but differently normalized picture appears. If \[ \lambda\sum_{m=1}^\infty |k_m|<1, \] then the only stationary solution is the isotropic one, \[ f(\theta)\equiv \frac1{2\pi}. \] For the classical Onsager kernel, since \[ \sum_{m=1}^\infty |k_m|=2, \] this gives the uniqueness regime \[ \lambda<\frac12. \] Because the first bifurcation for that kernel occurs at \(\lambda=3/2\), there is a gap between the low-\(\lambda\) uniqueness theorem and the first local bifurcation point. The phase-transition analysis at \(K_c=3\pi/4\) resolves an analogous threshold question in the mean-field minimization setting by proving exact coincidence of onset and loss of linear stability [1503.01477]. ## 5. Proof methods: coercivity, degree theory, and local bifurcation The recent proof of continuity for the planar Doi–Onsager transition is based on a sharp coercivity estimate derived from the constrained Lebedev–Milin inequality. In the normalized \(\frac{1}{n+1}\)-periodic setting, the key estimate is \[ H(q\mid q_u)\ge (n+1)\sum_{k=1}^\infty |k|^{-1}|q(k)|^2, \] with equality only for the one-parameter family \[ q(\theta)=\frac{1-c^2}{1+c^2-2c\cos\bigl(2\pi(n+1)(\theta-\theta_0)\bigr)}, \qquad c\in[0,1),\ \theta_0\in T. \] This yields the decomposition \[ \mathcal F_K(q)-\mathcal F_K(q_u) = \Bigl(H(q\mid q_u)-(n+1)\sum_{k\ge1}|k|^{-1}|q(k)|^2\Bigr) + \sum_{k\ge1}\bigl((n+1)|k|^{-1}-2KW(k)\bigr)|q(k)|^2, \] so nonnegativity follows from the entropy inequality and a Fourier decay condition. In the normalized theorem, if \(W\) is \(\frac{1}{n+1}\)-periodic, normalized by \(2W(n+1)=1\), and satisfies \[ 2W(k)\le \frac{n+1}{k}, \qquad \forall k\ge 1, \] then \[ K_c=K_\#=1, \] the phase transition is continuous, and if \(K\le \tfrac12\), then \(q_u\) is the unique critical point. For the Doi–Onsager kernel, the verification reduces to the explicit coefficient inequality \[ W_{DO}(2\ell)< \frac{W_{DO}(2)}{\ell},\qquad \ell\ge 2, \] equivalently \[ \frac{1}{4\ell^2-1}<\frac{1}{3\ell},\qquad \ell\ge 2, \] after rescaling by \(2W_{DO}(2)\) [2604.16288]. Earlier stationary analysis used a different toolkit. The low-parameter uniqueness theorem was obtained by applying Leray–Schauder degree theory to \[ F(V):=V-\lambda T(V). \] The proof combines compactness of \(T:H\to H\), an a priori bound for solutions, homotopy invariance of degree, \[ \deg(I-\lambda T,\Omega,0)=\deg(I,\Omega,0)=1, \] and the fact that every solution has index \(1\). The local bifurcation analysis assumes an ordering of negative Fourier coefficients, written as \[ k_11\); for \(m>1\) all bifurcated solutions are unstable, while for \(m=1\) they are stable if \(2k_2<1\) and unstable if \(2k_2>1\) [1503.01477]. Together, these results separate two analytical layers of the planar Doi–Onsager problem. Degree theory and local bifurcation describe existence and branch structure of stationary solutions, whereas the sharp coercivity argument identifies the exact mean-field phase-transition threshold and proves continuity at criticality. ## 6. Scope, variants, and the meaning of “two-dimensional” The expression “two-dimensional Doi–Onsager model” is not uniform across the literature. In the planar equilibrium and stationary works discussed above, two-dimensional means that the orientation of a rod is described by a single angle on \(S^1\), or equivalently by a probability measure on the circle, with \(\pi\)-periodicity representing head–tail symmetry [2604.16288]. By contrast, in the small Deborah number analysis without hydrodynamics, the physical space dimension may be \(d=2\) or \(3\), but the orientation variable remains \[ m\in S^2, \] and the model is a Smoluchowski-type equation on \(\mathbb R^d\times S^2\), not a planar orientation model on \(S^1\) [1707.05962]. The same distinction appears in kinetic-to-continuum derivations. The Bingham-closure derivation of a Q-tensor equation and the subsequent small Deborah number limit to the Ericksen–Leslie system are written in the standard 3D rod-like setting with orientational distribution \(f(x,m,t)\), \(m\in S^2\), symmetric traceless \(3\times 3\) Q-tensors, and 3D tensor identities. That framework is part of Doi–Onsager theory, but it is not a literal reduction to planar orientation space [1305.4721]. A common misconception is therefore to identify every “2D Doi–Onsager” reference with the same mathematical model. The planar circle model, the spatially two-dimensional Smoluchowski equation with \(S^2\)-valued orientations, and the Q-tensor or director limits derived from the 3D kinetic theory belong to the same Doi–Onsager hierarchy but address different variables, symmetries, and asymptotic regimes. The planar model distinguished by the kernel \(W_{DO}(\theta)=-|\sin(2\pi\theta)|\) is the one for which the exact threshold \[ K_c=K_\#=\frac{3\pi}{4} \] and the continuity of the phase transition are now rigorously established [2604.16288]. Source: https://www.emergentmind.com/topics/two-dimensional-doi-onsager-model