Degenerate Perturbation Theory Research Articles

A pendulum with a broken cylindrical symmetry: normal modes and beats

A pendulum, which has the cylindrical symmetry around its suspension wire, has two degenerate normal modes. We break this symmetry by putting a small magnet in the pendulum body and by running the pendulum in the presence of an inhomogeneous magnetic field. The force due to the magnetic field gradient modifies the equilibrium position and the spatial derivatives of this force define the normal modes of the pendulum with two different oscillation frequencies. We describe this system theoretically and we test these results by a series of experiments. This system provides a classical analogue of the degenerate perturbation theory in quantum mechanics.

Perturbation theory of nonlinear, non-self-adjoint eigenvalue problems: Semisimple eigenvalues

We discuss high-order adjoint based degenerate perturbation theory of semisimple eigenvalues resulting from nonlinear eigenvalue problems, which are typical of acoustic problems and vibrational analysis. The theory is discussed for a general matrix operator L with no further properties assumed. The unperturbed operator is thus in general non-normal, which requires the use of adjoint methods. The main challenge is to handle arbitrary order eigenvalue split for eigenvalue problems of arbitrary multiplicity: this is aided by a tree-structure formalism, which can be straightforwardly implemented numerically. We highlight the existence of a cross-branch solvability condition that needs to be satisfied when the perturbation causes the eigenvalues to split. After deriving the equations that govern the eigenproblem expansion order by order, we apply the method to three different test cases – which span eigenproblems of academic and industrial interest arising in structural mechanics and thermoacoustic stability – having different degrees of degeneracy and levels of complexity. We validate our algorithm by comparing the eigenvalues and eigenvectors reconstructed from perturbation theory up to 20th-order against exact solutions obtained using nonlinear eigenvalue solvers, and we discuss the advantages and limits of the method.

Low-temperature breakdown of many-body perturbation theory for thermodynamics

It is shown analytically and numerically that the finite-temperature many-body perturbation theory in the grand canonical ensemble has zero radius of convergence at zero temperature when the energy ordering or degree of degeneracy for the ground state changes with the perturbation strength. When the degeneracy of the reference state is partially or fully lifted at the first-order Hirschfelder-Certain degenerate perturbation theory, the grand potential and internal energy diverge as $T \to 0$. Contrary to earlier suggestions of renormalizability by the chemical potential $\mu$, this nonconvergence, first suspected by W. Kohn and J. M. Luttinger, is caused by the nonanalytic nature of the Boltzmann factor $e^{-E/k_\text{B}T}$ at $T=0$, also plaguing the canonical ensemble, which does not involve $\mu$. The finding reveals a fundamental flaw in perturbation theory, which is deeply rooted in the mathematical limitation of power-series expansions and is unlikely to be removed within its framework.

Complementary local-global approach for phonon mode connectivities

Sorting and assigning phonon branches (e.g., longitudinal acoustic) of phonon modes is important for characterizing the phonon bands of a crystal and the determination of phonon properties such as the Grüneisan parameter and group velocity. To do this, the phonon band indices (including the longitudinal and transverse acoustic) have to be assigned correctly to all phonon modes across a path or paths in the Brillouin zone. As our solution to this challenging problem, we propose a computationally efficient and robust two-stage hybrid method that combines two approaches with their own merits. The first is the perturbative approach in which we connect the modes using degenerate perturbation theory. In the second approach, we use numerical fitting based on least-squares fits to circumvent local connectivity errors at or near exact degenerate modes. The method can be easily generalized to other condensed matter problems involving Hermitian matrix operators such as electronic bands in tight-binding Hamiltonians or in a standard density-functional calculation, and photonic bands in photonic crystals.

ARC 3.0: An expanded Python toolbox for atomic physics calculations

ARC 3.0 is a modular, object-oriented Python library combining data and algorithms to enable the calculation of a range of properties of alkali and divalent atoms. Building on the initial version of the ARC library (Šibalić et al., 2017), which focused on Rydberg states of alkali atoms, this major upgrade introduces support for divalent atoms. It also adds new methods for working with atom–surface interactions, for modelling ultracold atoms in optical lattices and for calculating valence electron wave functions and dynamic polarisabilities. Such calculations have applications in a variety of fields, e.g., in the quantum simulation of many-body physics, in atom-based sensing of DC and AC fields (including in microwave and THz metrology) and in the development of quantum gate protocols. ARC 3.0 comes with an extensive documentation including numerous examples. Its modular structure facilitates its application to a wide range of problems in atom-based quantum technologies. Program summaryProgram Title: ARC 3.0CPC Library link to program files:https://doi.org/10.17632/c4z4n2cdf7.1Licencing provisions: BSD-3-ClauseProgramming language: PythonExternal Routines: NumPy [1], SciPy [1], Matplotlib [2], SymPy [3], LmFit [4]Nature of problem: The calculation of atomic properties of alkali and divalent atoms including energies, Stark shifts and dipole–dipole interaction strengths using matrix elements evaluated through a variety of means.Solution method: Dipole matrix elements are calculated using an analytical semi-classical approximation or wave functions obtained by numerical integration of the radial Schrödinger equation for a one-electron model potential. Interaction energies and shifts due to external fields are calculated using second order degenerate perturbation theory or exact diagonalisation of the interaction Hamiltonian, yielding results valid even at large external fields or small interatomic separation.Additional comments including restrictions and unusual features: External electric and magnetic field must be parallel to the quantisation axis. The accuracy of short range (≲1μm) atom - atom interaction potentials is limited by the truncation of the basis. Only weak magnetic fields are supported as only linear Zeeman shifts are taken into account. Calculations for divalent atoms use a single-electron approximation and calculation of their wave functions is not supported.

The structure of normal modes in parallel ideal optical fibers with strong coupling

In this paper, we studied an effect of strong evanescent coupling on the structure of normal modes in a system of parallel ideal multimode optical fibers. Using the formalism of the degenerate perturbation theory and a scalar waveguide equation for this system, analytical expressions of higher-order supermodes and their propagation constants have been determined. We have shown that the structure of modes in the case of strong evanescent coupling coincides with the structure of normal modes for weakly coupled parallel fibers. We have demonstrated that in the presence of strong coupling, expressions for corrections to the scalar propagation constant are modified, deducing them analytically.

Degenerate perturbation theory in thermoacoustics: high-order sensitivities and exceptional points

Abstract

Bosonic representation of a Lipkin-Meshkov-Glick model with Markovian dissipation

We study the dynamics of a Lipkin-Meshkov-Glick model in the presence of Markovian dissipation, with a focus on late-time dynamics and the approach to thermal equilibrium. Making use of a vectorized bosonic representation of the corresponding Lindblad master equation, we use degenerate perturbation theory in the weak-dissipation limit to analytically obtain the eigenvalues and eigenvectors of the Liouvillian superoperator, which in turn give access to closed-form analytical expressions for the time evolution of the density operator and observables. Our approach is valid for large systems, but takes into account leading-order finite-size corrections to the infinite-system result. As an application, we show that the dissipative Lipkin-Meshkov-Glick model equilibrates by passing through a continuum of thermal states with damped oscillations superimposed, until finally reaching an equilibrium state with a temperature that in general differs from the bath temperature. We discuss limitations of our analytic techniques by comparing to exact numerical results.

Managing disease outbreaks: The importance of vector mobility and spatially heterogeneous control.

Management strategies for control of vector-borne diseases, for example Zika or dengue, include using larvicide and/or adulticide, either through large-scale application by truck or plane or through door-to-door efforts that require obtaining permission to access private property and spray yards. The efficacy of the latter strategy is highly dependent on the compliance of local residents. Here we develop a model for vector-borne disease transmission between mosquitoes and humans in a neighborhood setting, considering a network of houses connected via nearest-neighbor mosquito movement. We incorporate large-scale application of adulticide via aerial spraying through a uniform increase in vector death rates in all sites, and door-to-door application of larval source reduction and adulticide through a decrease in vector emergence rates and an increase in vector death rates in compliant sites only, where control efficacies are directly connected to real-world experimentally measurable control parameters, application frequencies, and control costs. To develop mechanistic insight into the influence of vector motion and compliance clustering on disease controllability, we determine the basic reproduction number R0 for the system, provide analytic results for the extreme cases of no mosquito movement, infinite hopping rates, and utilize degenerate perturbation theory for the case of slow but non-zero hopping rates. We then determine the application frequencies required for each strategy (alone and combined) in order to reduce R0 to unity, along with the associated costs. Cost-optimal strategies are found to depend strongly on mosquito hopping rates, levels of door-to-door compliance, and spatial clustering of compliant houses, and can include aerial spray alone, door-to-door treatment alone, or a combination of both. The optimization scheme developed here provides a flexible tool for disease management planners which translates modeling results into actionable control advice adaptable to system-specific details.

Non-degenerate fourth order perturbation theory

There exist several forms of perturbation theory in Quantum Mechanics, namely: non-degenerate perturbation theory, degenerate perturbation theory, time-dependent perturbation theory, time-constant perturbation theory and time-harmonic perturbation theory. This paper presents an extension of the perturbation theory for non-degenerate states. It has been observed that most texts and journal papers treat only the first –order and second-order non-degenerate perturbation theory. Therefore, this paper attempts to treat and present the third order and fourth order non-degenerate perturbation theory in quantum mechanics. It can thus be asserted that the higher the order of perturbations of quantum systems that we know, the more successful will be our efforts in suppressing or eliminating them.

Analytical double-unitary-transformation approach for strongly and periodically driven three-level systems

Floquet theory combined with the generalized Van Vleck nearly degenerate perturbation theory has been widely employed for studying various two-level systems that are driven by external fields via time-dependent longitudinal (i.e., diagonal) couplings. However, three-level systems strongly driven by time-dependent transverse (i.e., off-diagonal) couplings have rarely been investigated, due to the breakdown of the traditional rotating wave approximation. Meanwhile, the conventional perturbation theory is not directly applicable, since a small parameter for the perturbed part is no longer apparent. Here we develop a double-unitary-transformation approach to deal with periodically modulated and strongly driven systems, where the time-dependent Hamiltonian has large off-diagonal elements. The first unitary transformation converts strong off-diagonal elements to diagonal ones, and the second enables us to harness the generalized Van Vleck perturbation theory to deal with the transformed Floquet matrix and also allows us to reduce the infinite-dimensional Floquet Hamiltonian to a finite effective one. For a strongly modulated three-level system, with the combination of the Floquet theory and the transformed generalized Van Vleck perturbation theory, we obtain analytical results of the system, which agree well with exact numerical solutions. This method offers a useful tool to analytically study multilevel systems with strong transverse couplings.

The Coherent State of the Landau Hamiltonian and the Relativistic Corrections to the Zeeman Effect in He<sup>+</sup> Ions

We calculate the energy levels of He+ ion placed in a uniform magnetic field directed perpendicular to the direction of its center of mass (CM) velocity vector, correct to relative order . Our calculations are within the frame work of an approximately relativistic theory, correct to relative order , of a two-particle composite system bound by electromagnetic forces, and written in terms of the position, momentum and spin operators of the constituent particles as proposed by Krajcik and Foldy, and also by Close and Osborn. Since the He+ ion has a net electric charge, the total or the CM momentum is not conserved and a neat separation of the CM and the internal motion is not possible. What is new in our approach is that, for the basis states in a first order degenerate perturbation theory to calculate the effects of the external magnetic field, we use the direct product of the coherent state of the Landau Hamiltonian of the He+ ion in a uniform magnetic field and of the simultaneous eigenstate of the internal Hamiltonian h, j2, l2, s2 and jz, where j, l and s are the internal total, orbital and spin angular moments of the He+ ion. The coherent state is an excellent approximation to the expected classical circular motion of the center of mass (CM) of the He+ ion. In addition to the Z2 a2 corrections to the usual nonrelativistic results, including the small corrections due to the nuclear motion, we also obtain corrections which depend on the kinetic energy (ECM ) of the CM circular motion of the He+ ion, in a nontrivial way. Even though these corrections are proportional to , where M is the mass of the He+ ion, and are small for nonrelativistic CM motion, the results should be verifiable in careful experiments. Our results may also have application in astrophysical observations of the spectral lines of He+ ions in magnetized astrophysical objects.

Ring dark solitons in three-dimensional Bose-Einstein condensates

In this work we present a systematic study of the three-dimensional extension of the ring dark soliton examining its existence, stability, and dynamics in isotropic harmonically trapped Bose-Einstein condensates. Detuning the chemical potential from the linear limit, the ring dark soliton becomes unstable immediately, but can be fully stabilized by an external cylindrical potential. The ring has a large number of unstable modes which are analyzed through spectral stability analysis. Furthermore, a few typical destabilization dynamical scenarios are revealed with a number of interesting vortical structures emerging such as the two or four coaxial parallel vortex rings. In the process of considering the stability of the structure, we also develop a modified version of the degenerate perturbation theory method for characterizing the spectra of the coherent structure. This semi-analytical method can be reliably applied to any soliton with a linear limit to explore its spectral properties near this limit. The good agreement of the resulting spectrum is illustrated via a comparison with the full numerical Bogolyubov-de Gennes spectrum. The application of the method to the two-component ring dark-bright soliton is also discussed.

Finite-temperature degenerate perturbation theory for bosons in optical lattices

Bosonic atoms confined in optical lattices can exist in two different phases, Mott-insulator and superfluid, depending on the strength of the system parameters, such as the on-site interaction between particles and the hopping parameter. This work is motivated by the fact that non-degenerate perturbation theory applied to the mean-field approximation of the Bose-Hubbard Hamiltonian at zero and finite temperature fails to give consistent results in the vicinity of the Mott-insulator-superfluid phase transition, e.g., the order parameter calculated via non-degenerate perturbation theory reveals an unphysical behavior between neighboring Mott lobes, which is an explicit consequence of degeneracy problems that artificially arise from such a treatment. Therefore, in order to fix this problem, we propose a finite-temperature degenerate perturbation theory approach based on a projection operator formalism which ends up solving such degeneracy problems in order to obtain physically consistent results for the order parameter near the phase transition.

Improving student understanding of fine structure corrections to the energy spectrum of the hydrogen atom

We discuss an investigation of student difficulties with fine structure corrections to the energy spectrum of the hydrogen atom in the context of degenerate perturbation theory. The investigation was carried out in advanced quantum mechanics courses by administering free-response and multiple-choice questions and conducting individual interviews with students. We find that students share many common difficulties related to these concepts. We describe how the research on student difficulties was used as a guide to develop and evaluate a Quantum Interactive Learning Tutorial (QuILT) which uses a guided inquiry-based approach to learning and strives to help students develop a functional understanding of concepts necessary for the fine structure corrections to the energy spectrum of the hydrogen atom in the context of degenerate perturbation theory. The QuILT strives to help students learn to focus on both the unperturbed and perturbation Hamiltonians to determine whether the basis is a good basis and then be able to find the first order corrections to the energy spectrum after ensuring that the basis chosen is a good basis and a good basis entails that the perturbation Hamiltonian be diagonal in each degenerate subspace of the unperturbed Hamiltonian.

Improving mean-field theory for bosons in optical lattices via degenerate perturbation theory

The objective of this paper is the theoretical description of the Mott-insulator to superfluid quantum phase transition of a Bose gas in an optical lattice. In former works the Rayleigh-Schr\"odinger perturbation theory was used within a mean-field approach, which yields partially non-physical results since the degeneracy between two adjacent Mott lobes is not taken into account. In order to correct such non-physical results we apply the Brillouin-Wigner perturbation theory to the mean-field approximation of the Bose-Hubbard model. Detailed explanations of how to use the Brillouin-Wigner theory are presented, including a graphical approach that allows to efficiently keep track of the respective analytic terms. To prove the validity of this computation, the results are compared with other works. Besides the analytic calculation of the phase boundary from Mott-insulator to superfluid phase, the condensate density is also determined by simultaneously solving two algebraic equations. The analytical and numerical results turn out to be physically meaningful and can cover a region of system parameters inaccessible until now. Our results are of particular interest provided an harmonic trap is added to the former calculations in an homogeneous system, in view of describing an experiment within the local density approximation. Thus, the paper represents an essential preparatory work for determining the experimentally observed wedding-cake structure of particle-density profile at both finite temperature and hopping.

Jahn-Teller effect in density-functional theory

Jahn and Teller used degenerate perturbation theory to prove that open shell molecules with symmetry Jahn-Teller distort. In that vein, we have developed a perturbative approach to computing Jahn-Teller distortions in Kohn-Sham density-functional theory, starting from a symmetrized molecule with an electronic shell closed by fractional occupation numbers. The resulting highly symmetric state is an energy extremum but not necessarily a minimum. Using second-order perturbation theory, we find the changes in geometry that occur when the symmetry of the electron density is broken to form a state with integer occupation numbers. This methodology allows us to retain many of the computational benefits of working in higher symmetry. As a demonstration, we solve the resulting equations for ten electrons in a superatomlike harmonic-oscillator potential.

A simple mechanism for unconventional superconductivity in a repulsive fermion model

Motivated by a scarcity of simple and analytically tractable models of superconductivity from strong repulsive interactions, we introduce a simple tight-binding lattice model of fermions with repulsive interactions that exhibits unconventional superconductivity (beyond BCS theory). The model resembles an idealized trilayer. The Cooper pair consists of electrons on opposite sides of the dielectric, which mediates the attraction. In the strong coupling limit, we use degenerate perturbation theory to show that the model reduces to a superconducting hard-core Bose-Hubbard model. Above the superconducting critical temperature, an analog of pseudo-gap physics results where the fermions remain Cooper paired with a large single-particle energy gap.

Fermi-Fermi crossover in the ground state of one-dimensional few-body systems with anomalous three-body interactions

In one spatial dimension, quantum systems with an attractive three-body contact interaction exhibit a scale anomaly. In this work, we examine the few-body sector for up to six particles. We study those systems with a self-consistent, non-perturbative, iterative method, in the subspace of zero total momentum. Exploiting the structure of the contact interaction, the method reduces the complexity of obtaining the wavefunction by three powers of the dimension of the Hilbert space. We present results on the energy, and momentum and spatial structure, as well as Tan's contact. We find a Fermi-Fermi crossover interpolating between large, weakly bound trimers and compact, deeply bound trimers: at weak coupling, the behavior is captured by degenerate perturbation theory; at strong coupling, the system is governed by an effective theory of heavy trimers (plus free particles in the case of asymmetric systems). Additionally, we find that there is no trimer-trimer attraction and therefore no six-body bound state.

Dissipative quantum mechanics beyond the Bloch-Redfield approximation: A consistent weak-coupling expansion of the Ohmic spin boson model at arbitrary bias

We study the time dynamics of the ohmic spin boson model at arbitrary bias $\epsilon$ and small coupling $\alpha$ to the bosonic bath. Using perturbation theory and the real-time renormalization group (RG) method we present a consistent zero-temperature weak-coupling expansion for the time evolution of the reduced density matrix one order beyond the Bloch-Redfield solution. We develop a renormalized perturbation theory and present an analytical solution covering the whole range from small to large times, including further results for exponentially small or large times. Resumming all secular terms in all orders of perturbation theory we find exponential decay for all terms of the time evolution. We determine the preexponential functions and find slowly varying logarithmic terms with the renormalized Rabi frequency $\Omega$ as energy scale together with strongly varying parts falling off asymptocially as $1/t$ in leading order, in contrast to the unbiased case. Resumming all logarithmic terms in all orders of perturbation theory via real-time RG we find the correct renormalized tunneling and a power-law behaviour for the oscillating modes with exponent crossing over from $2\alpha$ for exponentially small times to a bias-dependent value $2\alpha \epsilon^2/\Omega^2$ for exponentially large times. Furthermore, we present a degenerate perturbation theory to calculate consistently the purely decaying mode one order beyond Bloch-Redfield.