The emergence of chaos in an atom-field system is studied employing both semiclassical and numerical quantum techniques, taking advantage of the algebraic character of the Hamiltonian. A semiclassical Hamiltonian is obtained by considering the expectation value of the quantum Hamiltonian in Glauber (for the field) and Bloch (for the atoms) coherent states. Regular and chaotic regions are identified by looking at the Poincaré sections for different energies and parameter values. An analytical expression for the semiclassical energy density of states is obtained by integrating the available phase space, which provides an exact unfolding to extract the fluctuations in the level statistics. Quantum chaos is recognized in these fluctuations, as a function of the coupling strength, for different regions in the energy spectrum, evaluating the Anderson–Darling (A–D) parameter, which distinguishes the Wigner- or Poisson-like distributions. Peres lattices play a role similar to the Poincaré section for quantum states. They are calculated employing efficient numerical solutions and are a powerful visual tool to identify individual states belonging to a regular or chaotic region, classified by utilizing the Poincaré sections and the A–D parameter. Finally, the quantum Husimi function for selected excited states is shown to have a noticeable similitude with the Poincaré sections at the same energy.
The emergence of chaos in an atom-field system is studied employing both semiclassical and numerical quantum techniques, taking advantage of the algebraic character of the Hamiltonian. A semiclassical Hamiltonian is obtained by considering the expectation value of the quantum Hamiltonian in Glauber (for the field) and Bloch (for the atoms) coherent states. Regular and chaotic regions are identified by looking at the Poincaré sections for different energies and parameter values. An analytical expression for the semiclassical energy density of states is obtained by integrating the available phase space, which provides an exact unfolding to extract the fluctuations in the level statistics. Quantum chaos is recognized in these fluctuations, as a function of the coupling strength, for different regions in the energy spectrum, evaluating the Anderson–Darling (A–D) parameter, which distinguishes the Wigner- or Poisson-like distributions. Peres lattices play a role similar to the Poincaré section for quantum states. They are calculated employing efficient numerical solutions and are a powerful visual tool to identify individual states belonging to a regular or chaotic region, classified by utilizing the Poincaré sections and the A–D parameter. Finally, the quantum Husimi function for selected excited states is shown to have a noticeable similitude with the Poincaré sections at the same energy.
The nonintegrable Dicke model and its integrable approximation, the Tavis-Cummings model, are studied as functions of both the coupling constant and the excitation energy. The present contribution extends the analysis presented in the previous paper by focusing on the statistical properties of the quantum fluctuations in the energy spectrum and their relation with the excited-state quantum phase transitions. These properties are compared with the dynamics observed in the semiclassical versions of the models. The presence of chaos for different energies and coupling constants is exhibited, employing Poincar\'e sections and Peres lattices in the classical and quantum versions, respectively. A clear correspondence between the classical and quantum result is found for systems containing between N=80 and 200 atoms. A measure of the Wigner character of the energy spectrum for different couplings and energy intervals is also presented employing the statistical Anderson-Darling test. It is found that in the Dicke model, for any coupling, a low-energy regime with regular states is always present. The richness of the onset of chaos is discussed both for finite quantum systems and for the semiclassical limit, which is exact when the number of atoms in the system tends to infinite.
Abstract The interaction of a quantized electromagnetic field in a cavity with a set of two-level atoms inside it can be described with algebraic Hamiltonians of increasing complexity, from the Rabi to the Dicke models. Their algebraic character allows, through the use of coherent states, a semiclassical description in phase space, where the non-integrable Dicke model has regions associated with regular and chaotic motion. The appearance of classical chaos can be quantified calculating the largest Lyapunov exponent over the whole available phase space for a given energy. In the quantum regime, employing efficient diagonalization techniques, we are able to perform a detailed quantitative study of the regular and chaotic regions, where the quantum participation ratio ( P R ) of coherent states on the eigenenergy basis plays a role equivalent to the Lyapunov exponent. It is noted that, in the thermodynamic limit, dividing the participation ratio by the number of atoms leads to a positive value in chaotic regions, while it tends to zero in the regular ones.
Now, many different approaches have been presented to study the different semi-classical models derived from the Dicke Hamiltonian, which reflect a fact that the quantum-mechanical spin possesses no direct classical analog. The Hartree-Fock-type approximation is one of the widely used approaches, with which we derive the Heisenberg equations of motion for the system and replace the operators in these equations with the corresponding expectation values. In the present paper, we investigate the role of quantum phase transition in determining the chaotic property of the time-dependent driven Dicke model. The semi-classical Hamiltonian is derived by evaluating the expectation value of the Dicke Hamiltonian in a state, which is a product state of photonic and atomic coherent states. Based on the inverse of the relations between the position-momentum representation and the Bosonic creation-annihilation operators, the Hamiltonian is rewritten in the position-momentum representation and it undergoes a spontaneous symmetry-breaking phase transition, which is directly analogous to the quantum phase transition of the quantum system. In order to depict the Poincaré sections, which are used to analyze the trajectories through the four-dimensional phase space, we give the equations of motion of system from the derivatives of the semi-classical Hamiltonian for a variety of different parameters and initial conditions. According to the Dicke quantum phase transition observed from the experimental setup , we study the effect of a monochromatic non-adiabatic modulation of the atom-field coupling in Dicke model (i.e., the driven Dicke model) on the system chaos by adjusting the pump laser intensity. The change from periodic track to chaotic figure reflects the quantum properties of the system, especially the quantum phase transition point, which is a key position for people to analyse the shift from periodic orbit to chaos. In an undriven case, the system reduces to the standard Dicke model. We discover from the Poincaré sections that the system undergoes a change from the classical periodic orbit to a number of chaotic trajectories and in the superradiant phase area, the whole phase space is completely chaotic. When the static and driving coupling both exist, the system shows rich chaotic motion. The ground state properties are mainly determined by the static coupling, while the orbit of the system is adjusted by the driving coupling. If the static coupling is greater than the critical coupling, the system displays completely chaotic images in the Poincaré sections, and the periodic orbits in the chaos can also be adjusted by the strong driving coupling. While the static coupling is less than the critical coupling, the system can also show the chaotic images by adjusting the driving coupling strength and oscillation frequency. Moreover, by tuning the oscillation frequency, the Poincaré sections may change from the classical orbits to the chaos, and back to the classical orbits in the normal phase of the system.
In this thesis, we explore the Dicke model, well-known for its explanation of the superradiant
phase transition and signatures of Classical and Quantum Chaos. We
then propose a variational approach to determine the semi-classical (SC) equations
of motion, using coherent states. The SC dynamics exhibits signs of classical
chaos both with increasing coupling and as the energy of the system is gradually
increased, eventually leading to it completely engulfing the 2 dimensional phase
space of the pseudospin, which we examine through Poincare sections. We then
use numerical diagonalization to study the Quantum Dicke model - it is observed
that at and after exactly the critical value of the coupling, the system shows
a net macroscopic photon expectation value, and the pseudospin component in
the z direction starts to steadily decline - implying a net fraction of the atoms
are shifting to the excited state. We then study the dynamics of the Quantum
Dicke Hamiltonian using coherent states and compare the expectation values of
the dimensional with the classical trajectories from the semi-classical dynamics
for stable and chaotic orbits. We then introduce a linear quench which drives
the Dicke model across the Superradiant phase transition by changing the coupling
strength. The quenched dynamics is simulated using both semi-classical and
quantum dynamics - it is seen that dynamically varying quantities are periodically
oscillating around their instantaneous steady-state coupling values.
We study the nonintegrable Dicke model and its integrable approximation, the Tavis-Cummings model, as functions of both the coupling constant and the excitation energy. Excited-state quantum phase transitions (ESQPT) are found analyzing the density of states in the semiclassical limit and comparing it with numerical results for the quantum case in large Hilbert spaces, taking advantage of efficient methods recently developed. Two different ESQPTs are identified in both models, which are signaled as singularities in the semiclassical density of states; one static ESQPT occurs for any coupling, whereas a dynamic ESQPT is observed only in the superradiant phase. The role of the unstable fixed points of the Hamiltonian semiclassical flux in the occurrence of the ESQPTs is discussed and determined. Numerical evidence is provided that shows that the semiclassical results describe very well the tendency of the quantum energy spectrum for any coupling in both models. Therefore, the semiclassical density of states can be used to study the statistical properties of the fluctuation in the spectra, a study that is presented in a companion paper.
We present an analysis of chaos and regularity in the open Dicke model, when dissipation is due to cavity losses. Due to the infinite Liouville space of this model, we also introduce a criterion to numerically find a complex spectrum which approximately represents the system spectrum. The isolated Dicke model has a well-defined classical limit with two degrees of freedom. We select two case studies where the classical isolated system shows regularity and where chaos appears. To characterize the open system as regular or chaotic, we study regions of the complex spectrum taking windows over the absolute value of its eigenvalues. Our results for this infinite-dimensional system agree with the Grobe-Haake-Sommers (GHS) conjecture for Markovian dissipative open quantum systems, finding the expected 2D Poisson distribution for regular regimes, and the distribution of the Ginibre unitary ensemble (GinUE) for the chaotic ones, respectively.