Rydberg Atom Ion Interaction ============================ Here we look at the interacton between a single ion and a Rydberg atom. More details regarding the theory can be found in `A. Duspayev et al., Phys. Rev. Research 3, 023114 (2021) `__ or `M. Deiß et al., Atoms 9, 2 (2021) `__. The inernuclear axis is aligned with the z-direction and the interaction is expanded in multipolar terms. With this choice of coordinate system, the interaction potential is rotationally symmetric around the quantization axis and only states of same :math:`m_J` couple with each other, reducing the number of basis states. We will reproduce here Fig. 2 (a) and (b) from Phys. Rev. Research 3, 023114 (2021), showing the influence of higher multipole order terms Import the Library ------------------ .. code:: ipython3 # We call an IPython magic function to make the output of plotting commands displayed inline. %matplotlib inline # Arrays import numpy as np # Plotting import matplotlib.pyplot as plt # Operating system interfaces import os, sys # pairinteraction :-) from pairinteraction import pireal as pi # Create cache for matrix elements if not os.path.exists("./cache"): os.makedirs("./cache") cache = pi.MatrixElementCache("./cache") Application: Potential Curves ----------------------------- The SystemOne class defines the Rydberg atom and the effect of the electric field of the ion can be included into the calculations. The charge of the ion in units of the elementary charge can be passed by SystemOne.setIonCharge(charge) and the distance between the ion and the Rydberg core with SystemOne.setRydIonDistance(distance) in units of micrometer. The orientatin is fixed such that the internuclear axis points along the z-axis. As an example we show how to calculate the energy of the Rydberg state for different internuclear distances and including a different number of multipole orders. Choosing a maximum multipole order of one is equivalent to calculating a StarkMap in a homogenious electric field. .. code:: ipython3 # Define the Rydberg state for which the interaction with the ion should be calculated state = pi.StateOne("Rb", 45, 1, 1.5, 0.5) # Setup system, considering states only with similar energy and principle quantum number. # (if a high accuracy is required, the energy and principle quantum number ranges must be increased) LowN = 42 HighN = 48 LowEnergy = -100 HighEnergy = 100 system = pi.SystemOne(state.getSpecies(), cache) system.restrictEnergy(state.getEnergy()+LowEnergy, state.getEnergy()+HighEnergy) system.restrictN(LowN, HighN) # Since the ion and Rydberg atom are placed on the z-axis, the magnetic momentum is conserved system.restrictM(0.5,0.5) # Define the charge of the ion in units of the elementary charge system.setIonCharge(1) # Define the maximum considered order of the multipole expansion. # 1: monopole-dipole, 2: monopole-quadrupole, ... system.setRydIonOrder(1) # Loop over different distances to the ion array_distances = np.linspace(1.5,1.1, 100) # um array_eigenvalues = [] array_overlaps = [] for distance in array_distances: # Set the ion Rydberg-atom distance in units of um system.setRydIonDistance(distance) # Diagonalize the system system.diagonalize() # Store the eigenenergies array_eigenvalues.append(system.getHamiltonian().diagonal()) # Store the overlap of the eigenstates with the defined state array_overlaps.append(system.getOverlap(state)) array_eigenvalues = np.ravel(array_eigenvalues) array_overlaps = np.ravel(array_overlaps) array_distances = np.repeat(array_distances, system.getNumBasisvectors()) array_eigenvalues = array_eigenvalues / 29.9792458 # Convert GHz into inverse cm wavenumber # Plot the interaction potential, the color code visualizes overlap of eigenstates with defined state plt.scatter(array_distances, array_eigenvalues, 8, array_overlaps) plt.xlabel("Distance (um)") plt.ylabel("Energy (1/cm)") plt.colorbar(label="Overlap with {}".format(state)) plt.title("Up to 1st order in the multipole expansion") plt.grid(); plt.ylim(-61.5,-61.2); plt.show(); # Calculate the same but including orders up to 6 system.setRydIonOrder(6) # Loop over different distances to the ion array_distances = np.linspace(1.5,1.1, 100) # um array_eigenvalues = [] array_overlaps = [] for distance in array_distances: # Set the ion Rydberg-atom distance in units of um system.setRydIonDistance(distance) # Diagonalize the system system.diagonalize() # Store the eigenenergies array_eigenvalues.append(system.getHamiltonian().diagonal()) # Store the overlap of the eigenstates with the defined state array_overlaps.append(system.getOverlap(state)) array_eigenvalues = np.ravel(array_eigenvalues) array_overlaps = np.ravel(array_overlaps) array_distances = np.repeat(array_distances, system.getNumBasisvectors()) array_eigenvalues = array_eigenvalues / 29.9792458 # Convert GHz into inverse cm wavenumber # Plot the interaction potential, the color code visualizes overlap of eigenstates with defined state plt.scatter(array_distances, array_eigenvalues, 8, array_overlaps) plt.xlabel("Distance (um)") plt.ylabel("Energy (1/cm)") plt.colorbar(label="Overlap with {}".format(state)) plt.title("Up to 6th order in the multipole expansion") plt.grid(); plt.ylim(-61.5,-61.2); plt.show(); .. image:: atom_ion_interaction_files/atom_ion_interaction_6_0.png .. image:: atom_ion_interaction_files/atom_ion_interaction_6_1.png