A systematic, rigorous, and complete investigation of the Bloch equations in time-harmonic driving classical field is performed. Our treatment is unique in that it takes full advantage of the partial fraction decomposition over real number field, which makes it possible to find and classify all analytic solutions. Torrey's analytic solution in the form of exponentially damped harmonic oscillations [Phys. Rev. {\bf 76}, 1059 (1949)] is found to dominate the parameter space, which justifies its use at numerous occasions in magnetic resonance and in quantum optics of atoms, molecules, and quantum dots. The unorthodox solutions of the Bloch equations, which do not have the form of exponentially damped harmonic oscillations, are confined to rather small detunings $δ^2\lesssim (γ-γ_t)^2/27$ and small field strengths $Ω^2\lesssim 8 (γ-γ_t)^2/27$, where $γ$ and $γ_t$ describe decay rates of the excited state (the total population relaxation rate) and of the coherence, respectively. The unorthodox solutions being readily accessible experimentally are characterized by rather featureless time dependence.
A systematic, rigorous, and complete investigation of the Bloch equations in time-harmonic driving classical field is performed. Our treatment is unique in that it takes full advantage of the partial fraction decomposition over real number field, which makes it possible to find and classify all analytic solutions. Torrey's analytic solution in the form of exponentially damped harmonic oscillations [Phys. Rev. {\bf 76}, 1059 (1949)] is found to dominate the parameter space, which justifies its use at numerous occasions in magnetic resonance and in quantum optics of atoms, molecules, and quantum dots. The unorthodox solutions of the Bloch equations, which do not have the form of exponentially damped harmonic oscillations, are confined to rather small detunings δ2≲(γ−γt)2/27 and small field strengths Ω2≲8(γ−γt)2/27, where γ and γt describe decay rates of the excited state (the total population relaxation rate) and of the coherence, respectively. The unorthodox solutions being readily accessible experimentally are characterized by rather featureless time dependence.
The Bloch equation is the fundamental dynamical model applicable to arbitrary two-level systems. Analytical solutions to date are incomplete for a number of reasons that motivate further investigation. The solution obtained here for the propagator, which generates the time evolution of the system and embodies all the system dynamics, is compact and completely general. The parameter space that results in division by zero in previous treatments is explicitly defined and accommodated in the solution. Polynomial roots required for the solution are expressed in terms of a single real root obtained using simple functional forms. A simple graphical rendition of this root is developed that clarifies and characterizes its dependence on the physical parameters of the problem. As a result, the explicit time dependence of the system as a function of its physical parameters is immediately evident. Several intuitive models of system dynamics are also developed. In particular, the Bloch equation is separable in the proper coordinate system, written as the sum of a relaxation operator and either a null operator or a commuting rotation. The propagator thus drives either pure relaxation or relaxation followed by a rotation. The paper provides a basis for increased physical insight into the Bloch equation and its widespread applications.
We present a comprehensive analytical study of the optical Bloch equations for a closed two-level atom. We obtain a map in phase space of the Rabi frequency and the laser frequency detuning, discriminating the characteristics of solutions to the optical Bloch equations separated by a boundary in phase space. Outside the boundary, oscillatory behavior is shown, while simple decaying behavior is shown inside the boundary. The behavior of the solutions is also found to be dependent on the relative magnitude of the decay rates of the optical coherence and the excited state population.
Aspects of coherence and decoherence are analyzed within the optical Bloch equations. By rewriting the analytic solution in an alternate form, we are able to emphasize a number of unusual features: (a) despite the Markovian nature of the bath, coherence at long times can be retained; (b) the long-time asymptotic degree of coherence in the system is intertwined with the asymptotic difference in level populations; (c) the traditional population-relaxation and decoherence times, T1 and T2, lose their meaning when the system is in the presence of an external field, and are replaced by more general overall time scales; (d) increasing the field strength, quantified by the Rabi frequency Omega, increases the rate of decoherence rather than reducing it, as one might expect; and (e) maximum asymptotic coherence is reached when the system parameters satisfy Omega2=1(T1T2).
The Bloch equation and its variants constitute the fundamental dynamical model for arbitrary two-level systems. Many important processes, including those in more complicated systems, can be modeled and understood through the two-level approximation. It is therefore of widespread relevance, especially as it relates to understanding dissipative processes in current cutting-edge applications of quantum mechanics. Although the Bloch equation has been the subject of considerable analysis in the 70 years since its inception, there is still, perhaps surprisingly, significant work that can be done. This paper extends the scope of previous analyses. It provides a framework for more fully understanding the dynamics of dissipative two-level systems. A solution is derived that is compact, tractable, and completely general, in contrast to previous results. Any solution of the Bloch equation depends on three roots of a cubic polynomial that are crucial to the time dependence of the system. The roots are typically only sketched out qualitatively, with no indication of their dependence on the physical parameters of the problem. Degenerate roots, which modify the solutions, have been ignored altogether. Here the roots are obtained explicitly in terms of a single real-valued root that is expressed as a simple function of the system parameters. For the conventional Bloch equation, a simple graphical representation of this root is presented that makes evident the explicit time dependence of the system for each point in the parameter space. Several intuitive, visual models of system dynamics are developed. A Euclidean coordinate system is identified in which any generalized Bloch equation is separable, i.e., the sum of commuting rotation and relaxation operators. The time evolution in this frame is simply a rotation followed by relaxation at modified rates that play a role similar to the standard longitudinal and transverse rates. These rates are functions of the applied field, which provides information towards control of the dissipative process. The Bloch equation also describes a system of three coupled harmonic oscillators, providing additional perspective on dissipative systems.
This work presents a systematic numerical study of the dissipative Jaynes-Cummings model, consisting of a two-level system (qubit) coupled to a truncated harmonic oscillator, both subject to Markovian dissipation. The model is standard in quantum optics and open quantum systems. The master equation is solved using the Lindblad formalism with BDF integration. The study focuses on two standard observables: the qubit excited-state population P_exc(t) and the qubit Bloch vector length |r(t)|. A three-dimensional parameter scan is performed over the coupling g, the qubit dephasing rate gamma_phi, and the oscillator relaxation rate gamma_a, giving 216 cells in total. The parameters used are omega_0 = 6 Hz, omega_1 = 40 Hz, giving a detuning delta = 34 Hz, and thermal occupation n_th = 0.05. MAIN FINDINGS: The single-exponential description of the excited-state population is valid only at intermediate coupling (g approximately 4 to 8 Hz). Outside this region, the dynamics exhibits a multi-timescale structure. The decay time T is non-monotonic in g, with a minimum at g approximately 7 Hz. The location of the minimum depends only weakly on gamma_a, following a power law with exponent 0.127. At strong coupling (g greater than 15 Hz), a two-exponential model is strongly preferred by the Akaike Information Criterion in 43 out of 45 cells. The two timescales are separated by a factor of approximately 150. The fast timescale T_1 is best approximated by the Purcell formula (median ratio of numerical to theoretical value 0.68). The Fermi golden rule overestimates the decay rate by a factor of approximately 4.5. Neither the dressed-state correction nor the Bloch-Redfield calculation improves on the Purcell formula. The Bloch vector length does not follow any simple functional form (single exponential, stretched exponential, power law, or double exponential). All models have mean R-squared in the range 0.31 to 0.36. Its time derivative is, however, almost perfectly correlated with that of the population (mean correlation 0.91). The results are converged with respect to the Fock-space truncation at N_0 = 12, with identical observables for N_0 between 12 and 20. SCOPE: This is a purely numerical study of a well-defined quantum-optical model. It contains no empirical data, no claims about biological or cognitive systems, and no speculative ontology. The model, the numerical method, and the observables are standard. The results include several negative findings (the failure of perturbative formulas, the lack of a simple functional form for the Bloch vector) which are reported honestly and which identify the regime where non-perturbative methods are required. The purpose of this release is to provide reproducible numerical data and a benchmark for future analytical work on the dissipative Jaynes-Cummings model in the strong-coupling regime. CONTENTS: The LaTeX technical note (approximately 20 pages), the Python scripts for the numerical simulations, and the CSV files with the fit results, regime classification, and theoretical comparison. REFERENCES: E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963).G. Lindblad, Commun. Math. Phys. 48, 119 (1976).H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press (2002).C. W. Gardiner and P. Zoller, Quantum Noise, 3rd ed., Springer (2004).E. M. Purcell, Phys. Rev. 69, 681 (1946).C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom-Photon Interactions, Wiley (1992).W. P. Schleich, Quantum Optics in Phase Space, Wiley-VCH (2001).