Abstract:Lipkin-Meshkov-Glick model is paradigmatic to study quantum phase transition in equilibrium or non-equilibrium systems and entanglement dynamics for a variety of disciplines. In thermodynamics limit, quantum fluctuations are negligible, its semiclassical dynamics in presence of only one nonlinear couplings, as a good benchmark to study quantum fluctuation in finite-size system, can be well obtained in terms of Jacobi elliptic functions. In this work, we extend this semiclassical analysis into the regime where both nonlinear interactions are present, and successfully obtain its exact solutions of semiclassical equations by constructing an auxiliary function that is a linear combination of the $y$ and $z$ component of the classical spin in thermodynamic limit. Taking implementation of Lipkin-Meshkov-Glick model in a Bose-Einstein condensate setup as an example, we figure out all classical dynamical modes, specially find out mesoscopic self-trapping mode in population and phase-difference space even persists in presence of both nonlinear couplings. Our results would be useful to analyze dynamical phase transitions and entanglement dynamics of Lipkin-Meshkov-Glick model in presence of both nonlinear couplings.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the semiclassical dynamics of the Lipkin - Meshkov - Glick (LMG) model in the presence of two nonlinear couplings. Specifically, the author extended the previous semiclassical analysis method which only considered a single nonlinear coupling, and successfully obtained the exact solution of the LMG model containing two nonlinear coupling terms. By constructing an auxiliary function, which is a linear combination of the y - component and z - component of the classical spin in the thermodynamic limit, the author was able to describe the dynamic behavior of the system and found that the mesoscopic self - trapping mode still exists in the particle number and phase difference space even when the two nonlinear couplings exist simultaneously.
### Background of the Paper and Problem Definition
The Lipkin - Meshkov - Glick (LMG) model was initially introduced in nuclear physics to describe an N - spin system with infinite - range interactions. This model can also be regarded as a special limit of the XY Heisenberg model in condensed - matter physics, and has been engineered on multiple platforms in recent years, such as Bose - Einstein condensates (BEC) in double - well traps, two - component BECs, quantum gases in cavity quantum electrodynamics (cavity QED), thermal Bose gases, and superconducting circuit QED, to explore its potential applications in quantum information and precision measurement.
Although the exact solution of the integrable LMG model has been obtained, many studies are still dedicated to discovering unexplored singular phenomena by adding new aspects, such as chaotic dynamics caused by non - equilibrium quantum fluctuations. For the general form of the LMG model:
\[
\hat{H}=h \hat{J}_{x}+\frac{1}{N}\left(g_{1} \hat{J}_{y}^{2}+g_{2} \hat{J}_{z}^{2}\right)
\]
where \(\hat{J}_{\alpha}=\sum_{j = 1}^{N}\frac{1}{2}\sigma_{\alpha}^{j}\) (\(\alpha=x, y, z\)) are collective spin operators, satisfying the \(SU(2)\) algebra, \([\hat{J}_{x}, \hat{J}_{y}]=i \hat{J}_{z}\) (\(\hbar\equiv1\) in this paper), \(h\) is the magnetic field in the x - direction, and \(g_{1}, g_{2}\) are the nonlinear interaction strengths of the spin in the y(z) - direction. Since \(\hat{J}^{2}=\sum_{\alpha = x, y, z}\hat{J}_{\alpha}^{2}=J(J + 1)\) is a conserved quantity, most attention is focused on the part with the maximum angular momentum \(J = N/2\).
When \(g_{1}=0\) or \(g_{2}=0\), the model has \(Z_{2}\) symmetry and is widely used to generate spin squeezing and multi - body entanglement. However, when \(g_{1}\neq0\) and \(g_{2}\neq0\), the model has \(Z_{2}\times Z_{2}\) symmetry, but its analytical solution for semiclassical dynamics is still lacking so far.
### Research Objectives
The main objective of this paper is to obtain the exact solution of the LMG model containing two nonlinear coupling terms by constructing an auxiliary function, which is a linear combination of the \(\langle\hat{J}_{y}\rangle\) and \(\langle\hat{J}_{z}\rangle\) components. Specifically, the author hopes:
1. **Analytical Solution**: Analytically solve the semiclassical dynamic equations by constructing an auxiliary function.
2. **Dynamic Modes**: Explore the dynamic modes under different parameters, especially the mesoscopic self - trapping mode.
3. **Frequency Analysis**: Give the evolution frequency of the system and compare it with the numerical results.
### Main Contributions
1. **Analytical Solution**: Successfully obtained the analytical solution of the semiclassical dynamics of the LMG model containing two nonlinear coupling terms.
2. **Dynamic Modes**: Discovered that even when the two nonlinear couplings exist simultaneously...