An Update to Isomers of Rydberg Excitations in Argon Clusters

This paper reports an improved Diatomic-In-Molecules (DIM) calculation for excited argon clusters that incorporates previously ignored strongly avoided crossings between 3p4s and 3p4p states to better understand the localization of excitation and the effect of diabatisation on cluster isomers.

The Gist

Imagine a group of Argon atoms hanging out together in a cluster. Usually, they are calm and quiet. But sometimes, one of them gets a little "excited" (like a person jumping up and down with energy). This paper is about figuring out exactly how that energy is shared among the group and what shape the group takes when this happens.

For a long time, scientists thought the excited energy was shared by a trio of atoms (a trimer) sitting right in the middle of the cluster. Think of it like a three-person huddle where everyone is holding hands and sharing a secret.

However, the authors of this paper found a problem with that old idea. They realized that the math they were using to predict this behavior was missing a crucial piece of the puzzle: a "traffic jam" in the energy levels.

Here is a breakdown of their work using simple analogies:

1. The Old Map vs. The New Map

• The Old Way (DIM method): Imagine trying to navigate a city using an old map that ignores a massive construction zone. The map told scientists that the excited energy was spread out over three atoms (a trimer).
• The Better Way (HPP method): A few years ago, the authors used a more detailed, high-tech GPS (called the HPP method). This GPS showed that the energy wasn't shared by three atoms; it was actually stuck on just two atoms (a dimer), like a pair of dancers spinning together while the rest of the crowd watches.
• The Problem: The high-tech GPS (HPP) is incredibly accurate but very slow and expensive to run. It's like having a super-precise but heavy tank that can't move fast enough to predict how the atoms will dance in real-time. The old map (DIM) is fast and light, but it was giving the wrong directions because it missed the "construction zone."

2. The "Traffic Jam" (Avoided Crossing)

The reason the old map was wrong is that two energy paths were trying to cross each other but couldn't quite do it. In physics, this is called an "avoided crossing."

• The Analogy: Imagine two cars on a highway trying to switch lanes. If they try to switch at the exact same spot, they crash. Instead, one car swerves up and the other swerves down to avoid the crash.
• The Mistake: The old math treated these two paths as if they were straight, separate lanes that never touched.
• The Fix: The authors realized they needed to account for that "swerve." They introduced a technique called Diabatisation. Think of this as drawing a new, smooth curve on the map that connects the two lanes correctly, acknowledging that they influence each other even when they don't crash.

3. The "Dummy" State

To fix the math without needing the super-slow, expensive GPS, the authors had to invent a "placeholder" or a "dummy" state.

• The Analogy: Imagine you are trying to balance a scale, but you don't know the weight of one of the objects. So, you put a "dummy" weight on the other side that you adjust until the scale balances perfectly.
• In this paper, they created a fake, made-up energy state (an ad hoc state) to help the math work out. It's not a "real" physical state they found, but it acts like a mathematical tool to make the equations behave correctly.

4. What They Found

When they used this new, improved "fast map" (Di-DIM) with the traffic jam fixed:

• The Shape Changed: They confirmed the old GPS finding: The excited energy lives on a pair of atoms (a dimer), not a trio.
• The Dance: The excited pair attaches itself to the rest of the cluster (the ground state atoms). It's like a glowing pair of dancers attaching themselves to a large group of people standing still.
• The Details: While the new map got the main shape right, it wasn't perfect.
  • The distance between the excited pair and the rest of the group was slightly shorter than the high-tech GPS predicted.
  • In some cases, the excited pair tilted slightly to the side (breaking symmetry), whereas the high-tech GPS showed them sitting perfectly straight. The authors admit this is because their "fast map" still misses some subtle forces (like polarization) that the "slow map" catches.

5. The Bottom Line

The authors successfully updated the "fast map" (DIM method) so that it now agrees with the "high-tech GPS" (HPP) on the most important fact: The excited energy in Argon clusters lives on a pair of atoms, not a trio.

They achieved this by fixing the "traffic jam" in the math using a clever trick with a "dummy" state. While their new map isn't 100% perfect on every tiny detail (like exact distances or tilting), it is now good enough to be used for fast, real-time simulations of how these excited atoms move and dance, which was the main goal of the study.

Technical Summary

Technical Summary: An Update to Isomers of Rydberg Excitations in Argon Clusters

Problem Statement
The characterization of the lowest energy isomers of excited argon clusters ($Ar^*_N$) has been a subject of debate. Early work by Naumkin et al. using the Diatomic-In-Molecules (DIM) method suggested that the excitation in the lowest energy isomers is localized over a trimer ($Ar^*_3$) with delocalized character. However, recent studies utilizing the Hole-Particle-Pseudopotential (HPP) method indicated that the excitation is actually localized on a dimer ($Ar^*_2$), forming an excimer attached to a ground-state cluster ($Ar^*_2-Ar_{N-2}$). While HPP provides accurate insights into spectroscopic properties and isomer geometries, it is computationally expensive and unsuitable for on-the-fly potential energy curve calculations required for quantum dynamics. Conversely, the DIM method is efficient for dynamics but previously failed to reproduce the HPP results, likely due to the neglect of strong avoided crossings between specific electronic states.

Methodology
This study presents an updated DIM framework, termed Di-DIM (Diabatised-DIM), designed to resolve the discrepancies between previous DIM calculations and HPP results. The core methodological advancement involves the explicit treatment of strong avoided crossings between the $3p4s$ and $3p4p$ $^1,3\Sigma$ states, which were previously ignored or treated as pure diabatic states.

1. Diabatisation Scheme: The authors implement a unitary transformation to decouple adiabatic states. They introduce an ad hoc diabatic state ($\beta$) to emulate the coupling of the $3p4s$ $\Sigma_g$ state with higher $3p4p$ states. This allows for the restoration of the correct diabatic character of the potential energy curves (PECs).
2. Hamiltonian Construction: The Di-DIM Hamiltonian incorporates spin-orbit couplings (SOC) and utilizes PECs derived from the HPP method for the $Ar^*_2$ dimer. The basis set is expanded to include $3p^54s$ and $3p^4p$ configurations, increasing the configuration space from 12 to 16 determinants.
3. Optimization: The lowest energy isomers for cluster sizes $N=2$ to $15$ and the magic number $N=55$ are determined using damped molecular dynamics. The electronic structure is analyzed by examining the hole distribution derived from the expansion coefficients of the atomic basis functions.

Key Contributions and Results

• Resolution of Trimer Discrepancy: The study demonstrates that without diabatisation, the DIM method yields a symmetric $D_{\infty h}$ trimer isomer with delocalized excitation, consistent with earlier literature. However, the Di-DIM method, by including the avoided crossing treatment, yields a minimum with $C_{\infty h}$ symmetry and a linear asymmetric geometry. This result aligns with the HPP finding that the excitation is localized on a dimer rather than a trimer.
• Localization of Excitation: For clusters $N \geq 3$, the Di-DIM method confirms that the excitation is localized on a protruding dimer ($Ar^*_2$), attached to a ground-state cluster of $N-2$ atoms. This contrasts with previous DIM results where 80% of the excitation was reported on the central atom of a trimer. In the Di-DIM model, the hole is equally distributed between the two outermost atoms of the excimer.
• Geometric and Energetic Comparisons:
  • The dissociation energies calculated with Di-DIM are lower than those from HPP and CASPT2 but significantly lower than previous non-diabatic DIM calculations.
  • The inter-atomic distance of the excimer in Di-DIM ($R_e = 4.52$ au) is shorter than in HPP ($6.10$ au) and CASPT2 ($5.90$ au).
  • For $N=4$, Di-DIM predicts a linear asymmetric isomer, differing from the symmetric $D_{\infty h}$ structure found in HPP, a discrepancy attributed to the lack of polarization effects in the DIM method.
  • For larger clusters (e.g., $N=9$), Di-DIM reproduces the general trend of an excimer attached to a ground-state cluster but observes symmetry breaking (tilting of the excimer) not present in the HPP model, likely due to the omission of rotational effects of P and D states.
• Cluster Size Trends: The study provides a comprehensive table of isomers up to $N=15$ and $N=55$. It observes that while the general motif of an excimer attached to a ground-state cluster holds, the specific geometries are a mixture of previous DIM and HPP predictions.

Significance and Limitations
The paper claims that the Di-DIM method, even with a crude approximation of the ad hoc state, successfully corrects the primary character of the lowest excited isomers to agree with the more rigorous HPP results. This validates the use of diabatisation to improve the transferability of the DIM method for studying excited state dynamics.

The authors remain modest regarding the quantitative accuracy of the method. They acknowledge that Di-DIM cannot fully replicate the dipole moments or the exact geometries (such as the specific symmetry breaking or bond lengths) observed in HPP due to the absence of polarization effects and the exclusion of higher excited states and their couplings. However, the method is deemed sufficient for studying the localization of excitation and the general dynamics of excited argon clusters, providing a computationally efficient alternative to HPP for on-the-fly calculations. The authors conclude that future work will investigate how the dynamics change if these specific couplings are ignored.