## Molecular influence on nuclear-quadrupole-coupling effects in laser induced alignment

Linda V. Thesing, 1, 2, 3 Andrey Yachmenev, 1, 2, a) Rosario González-Férez, 4, b) and Jochen Küpper 1, 2, 3, c)

- 1) Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
- 2) Center for Ultrafast Imaging, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- 3) Department of Physics, Universität Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany
- 4) Instituto Carlos I de Física Teórica y Computacional and Departamento de Física Atómica, Molecular y Nuclear, Universidad de Granada, 18071 Granada, Spain

(Dated: 2024-08-05)

We studied the effect of nuclear-quadrupole interactions on the field-free impulsive alignment of different asymmetric-top molecules. Our analysis is focused on the influence of the hyperfine- and rotational-energy-level structures. These depend on the number of nuclear spins, the rotational constants, and the symmetry of the tensors involved in the nuclear spin and external field interactions. Comparing the prototypical large-nuclearspin molecules iodobenzene, 1,2-diiodobenzene, 1,3-diiodobenzene, and 2,5-diiodobenzonitrile, we demonstrate that the magnitude of the hyperfine splittings compared to the rotational-energy splittings plays a crucial role in the spin-rotational dynamics after the laser pulse. Moreover, we point out that the impact of the quadrupole coupling on the rotational dynamics decreases when highly excited rotational states dominate the dynamics.

## I. INTRODUCTION

The coupling of rotational and nuclear spin angular momenta, known as hyperfine interaction, can have a significant influence on the rotational dynamics of molecules and is at the core of the nuclear magnetic resonance measurement. 1-9 There are several mechanisms of hyperfine interaction, the most prominent is the nuclearquadrupole coupling. This occurs in molecules containing nuclei with partly filled shells, resulting in non-spherical nuclear shapes and the appearance of an electric nuclear quadrupole moment. This quadrupole couples, electrostatically interacts, with the field gradient induced by other moving charged particles in the molecule, i. e., the nuclei and electrons.

under field-free conditions is often favorable. Field-free alignment relies on the preparation of rotational wave packets by means of intense optical laser pulses. 17 Since the wave packets are non-stationary states, they dephase and rephase periodically, changing the alignment with time, even after the pulsed excitation. Typical pulse shapes range from the short kick pulses 18-20 to pulses with long rising and short falling edges. 21-25 The timeevolution of the molecular alignment in field-free conditions is revival structure, ideally with nearly the same degree of alignment appearing every rotational period. The latter is, however, only valid for molecules with regular spacings between the rotational energy levels, such as rigid linear and symmetric top molecules.

Generally, nuclear-quadrupole coupling is strong in molecules containing heavy atoms, such as Se, Br, I, Fe, Au, or Pt, which are also found in many biologically relevant molecules. Due to their very large x-ray- and electron-scattering cross sections these atoms are often chemically attached to molecular compounds and used as marker atoms in diffractive imaging experiments. 10-13 Furthermore, in Coulomb-explosion imaging these heavy atoms often show good axial recoil enabling the observation of molecular orientation and internal structure during dynamical processes. 14-16

An indispensable component of high-resolution imaging experiments is the laser-induced alignment of molecules, enabling measurements to be performed in the moleculefixed frame without orientation averaging. To avoid perturbations due to the presence of external fields, alignment

a) Email: andrey.yachmenev@cfel.de

b) Email: rogonzal@ugr.es

c) Email: jochen.kuepper@cfel.de; website: https://www.controlledmolecule-imaging.org

The revivals of molecular alignment are highly sensitive to the molecule's rotational energy level structure. A detailed understanding of various molecular degrees of freedom that can couple to rotations is thus needed to accurately predict the alignment dynamics of molecules. It is well established that hyperfine interactions can have a significant impact on the rotational dynamics of molecules prepared in single rotational states. 9 Furthermore, experimental evidence of the influence of nuclear quadrupole effects on the impulsive alignment dynamics was published recently for I 2 molecules. 26 In computational studies of linear and asymmetric top molecules, we found a nontrivial dependence of the nuclear quadrupole effects on the intensity of the laser field. 27,28

Here, provide a deeper insight into the effect of the quadrupole coupling on the alignment for different molecular species. We compared the results across asymmetric top molecules with different number of nuclear quadrupolar nuclei, their different positions in the molecule, and different molecular symmetries. Our molecular set contains C 2 v -symmetric molecules with one and two iodine atoms, i. e., iodobenzene and the 1,2- and 1,3-isomers of diiodobenzene, besides the already investigated 1,4-isomer,

FIG. 1. Sketches of the molecules IB, DIBN, 1,3-DIB and 1,2DIB. For IB and 1,3-DIBN the most polarizable axis (MPA) is parallel to the a axis, for DIBN it forms an angle of 5 ◦ with the a axis and for 1,2-DIB the MPA is parallel to the b axis. The molecule-fixed z axis is indicated in the respective sketches.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-82065d63-2f1c-44b2-8abc-9f1c3df3dde1.jpg)

and the C -symmetric molecules diiodobenzonitrile and s 1,3-bromoiodobenzene. In addition, we consider a hypothetical modified single-spin 1,3-diiodobenzene. First, we summarize our theoretical approach in Section II. We then present impulsive alignment results for low- (Section III A) and high-intensity (Section III B) laser pulses for a rotational temperature of T rot = 0 K. Finally, we analyze the rotational dynamics of cold thermal ensembles and individual excited states.

## II. THEORETICAL DESCRIPTION

We studied the laser-induced rotational dynamics of asymmetric top molecules taking into account the nuclearquadrupole interactions. 27 The field-free spin-rotational Hamiltonian of the system is given by

$$\hat { H } _ { \text{mol} } = \hat { H } _ { \text{rot} } + \sum _ { l } \mathbf V ( l ) \cdot \mathbf Q ( l ), \quad \quad ( 1 ) \quad \text{ of a} \\ \text{ the} \quad \text{$cole}.$$

where ˆ H rot = AJ ˆ 2 a + BJ ˆ 2 b + CJ ˆ 2 c is the rigid-rotor Hamiltonian with the rotational constants A,B,C and the components ˆ J i , i = a, b, c , of the rotational angular momentum operator ˆ J . The second term in the Hamiltonian (1) describes the interaction of the quadrupole-moment tensor Q ( ) l of the l -th nucleus with the electric-field-gradient (EFG) tensor V ( ) l . The sum runs over all nuclei with significant quadrupole moments l in the molecule.

The interaction with the nonresonant laser field, linearly-polarized along the laboratory-fixed Z axis, can be written as

$$\hat { H } _ { \text{las} } ( t ) = - \frac { I ( t ) } { 2 \varepsilon _ { 0 } c } \alpha _ { Z Z }$$

where I ( ) t is the intensity of the laser field and α ZZ is the corresponding element of the molecule's polarizability tensor in the laboratory-fixed frame.

The geometrical structures, principal axes of inertia, and most-polarizable (MPA) axes of iodobenzene (IB), diiodobenzonitrile (DIBN), 1,3-diiodobenzene (1,3-DIB) and 1,2-diiodobenzene (1,2-DIB) are shown in Fig. 1. All four molecules contain one or two 127 I nuclei, which produce considerably large hyperfine energy splittings of the rotational states due to iodine's large quadrupole moment eQ = -696 mb. 29 The IB molecule has only one iodine nucleus and the total spin angular momentum is ˆ = ˆ I I 1 , where I 1 = 5 2 / . For the molecules with two iodine nuclei ˆ = ˆ I I 1 + ˆ I 2 ; the associated quantum number can have values I = 0 , ..., 5 . The total angular momentum is given by ˆ = ˆ + ˆ F J I . The nitrogen nucleus in DIBN has an additional nonzero quadrupole moment eQ = 20 44 . mb, 29 which is, however, relatively small compared to iodine and thus neglected for the purpose of this study. This could also be achieved by the 15 N isotopologue with spin 1/2 and thus eQ = 0 . For 1,3-bromoiodobenzene (1,3-BIB) we assume two different quadrupole-coupled nuclei 127 I and 79 Br, the more abundant isotope of bromine with Q = 313 mb. 29 All important parameters for the molecules in this study, such as rotational constants, polarizability, and EFG tensor elements are listed in Appendix A.

The simulations of the spin-rotational dynamics of molecules subject to external electric fields were performed using the variational approach Richmol 27,30 and the same computational setup as described in our previous work. 28 The total Hamiltonian is composed of the sum ˆ H mol + ˆ H las and the total time-dependent wave function is expanded in the basis of eigenfunctions ∣ ∣ ∣ F, J ˜ ˜ K K a ˜ c , n, M 〉 of the fieldfree time-independent Hamiltonian ˆ H mol . The details of the functional form of these basis functions are described in Appendix B. The alignment of the molecules is quantified by the expectation value 〈 cos 2 θ 〉 , where θ is the angle between the laser polarization axis and the MPA of the molecule. In the present setup, the laser polarization is aligned along the laboratory-fixed Z axis and the MPA of all molecules is strictly or nearly parallel to one of the principal axes of inertia, see Fig. 1. This simplifies calculations of the expectation value 〈 cos 2 θ 〉 , since the θ angle can be associated with the corresponding Euler angle.

We compared the post-pulse molecular alignment with and without the nuclear-quadrupole-coupling term in (1) for initial rotational temperatures T rot = 0 0 1 , . , and 0 3 . K. We used the eigenstates of ˆ H mol and ˆ H rot as initial field-free molecular states in the calculations with and without quadrupole coupling, respectively.

We assumed equal populations for all hyperfine components of a single rotational state, mimicking experimental conditions in molecular beams, where nuclear-spin states typically do not cool and nuclear-spin temperatures correspond to the high-temperature limit before expansion. Thus, even for rotational temperatures T rot = 0 K the alignment was calculated by averaging over the results of independent calculations for every hyperfine component of the ground rotational state J = 0 . For the finitetemperature calculations further averaging was performed by Boltzmann weighting the rotational states.

## III. RESULTS AND DISCUSSION

## A. Low-intensity impulsive alignment

Previously, we found that the effect of the nuclearquadrupole coupling strongly depends on the laser intensity and is larger at a lower intensity. 27,28 Thus we begin by analyzing the influence of the quadrupole coupling on the impulsive alignment induced by a 1 ps (FWHM) Gaussian laser pulse with a moderate peak intensity I 0 = 10 11 W cm . Fig. 2 shows the temporal / 2 evolution of the post-pulse alignment dynamics, characterized by 〈 cos 2 θ 〉 ( ) t , computed for different molecules at T rot = 0 K. The results including nuclear-quadrupole coupling are plotted with thick lines and the results of rigidrotor calculations, i. e., neglecting the nuclear-quadrupole coupling, are plotted with thin black lines.

In the rigid-rotor case, the alignment revival structures in Fig. 2 look fairly regular for all molecules, which is the result of dephasing and rephasing of rather narrow rotational wave packets. The revivals change from one molecule to another due to their distinct rotationalenergy-level structures and, to some extent, due to slightly different interaction strengths with the laser field caused by the different polarizabilities. Due to the weak intensity of the field, the post-pulse dynamics of IB, DIBN and 1,3-DIB is dominated by the coherences between the lowenergy rotational states with K a = 0 . K a is the quantum number of the rotational angular momentum projection operator onto the molecule-fixed a axis. Since the energies of these states look fairly similar to the energy levels of a symmetric top, the revival patterns in Fig. 2 a-c repeat themselves after multiples of the respective rotational periods τ IB rot = 707 69 . ps, τ DIBN rot = 3365 7 . ps and τ 1,3-DIB rot = 2499 4 . ps, with τ rot = 1 ( / B + C ) . In the case of 1,2-DIB, the post-pulse dynamics (Fig. 2 d) does not show such a periodic behavior, due to the molecule's larger rotational asymmetry, see Table I.

The nuclear-quadrupole coupling increases the complexity of the post-pulse alignment dynamics. At early times after the laser pulse 〈 cos 2 θ 〉 does not show any significant deviation from the rigid-rotor dynamics as the quadrupole coupling is much weaker than the polarizability interaction with the laser field. However, the wavepackets quickly dephase, 28 characteristic of the non-

FIG. 2. Impulsive alignment induced by a Gaussian laser pulse with τ FWHM = 1 ps and I 0 = 10 11 W cm / 2 for (a) IB, (b) DIBN, (c) 1,3-DIB and (d) 1,2-DIB including (QC) and neglecting the nuclear-quadrupole coupling (noQC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-b1b7508b-6523-4db4-b093-dfbd094a4ce4.jpg)

regular hyperfine splittings of low-energy rotational states that dominate the dynamics. Each rotational state of the IB molecule is split into 6 hyperfine components, at most, while for the molecules with two iodine nuclei each state has up to 36 hyperfine components. The influence of the quadrupole coupling is noticeably stronger for DIBN and 1,3-DIB molecules, see Fig. 2 b, c. This is due to the rather small B and C rotational constants for these molecules, corresponding rotational-level spacings in DIBN and 1,3DIB at low rotational excitations that are comparable to or even smaller than the hyperfine splittings. This yields strong mixing of different rotational levels, even in the initial field-free states, giving rise to large dephasing effects. A similar behavior was observed previously for 1,4-DIB. 28 For IB and 1,2-DIB, the hyperfine splittings are much

FIG. 3. Impulsive-alignment dynamics for I 0 = 10 11 W cm / 2 with nuclear-quadrupole coupling (QC) for the initial states (a) ∣ ∣ 5 2 0 / , 00 , 1 2 / 〉 of IB and (b) | 0 0 , 00 , 0 ⟩ of 1,3-DIB and without the coupling (noQC). (c, d) Fast-Fourier transforms of 〈 cos 2 θ 〉 with and without quadrupole coupling.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-b498af18-66bf-4016-83b7-1df16a234c8f.jpg)

smaller than the rotational level spacings and thus the dephasing effect is weaker. Consequently, the alignment for these molecules oscillates at frequencies similar to the rigid-rotor results, see Fig. 2 a, d. The effect of the quadrupole coupling is, however, visible in the strongly modified amplitude of the degree of alignment.

In order to illustrate the interplay between the hyperfine-split and the pure rotational energy level spacings, we computed the fast-Fourier transforms (FFTs) of 〈 cos 2 θ 〉 ( ) t for IB and 1,3-DIB molecules, see Fig. 3. Here, the wavepackets obtained from the single initial hyperfine state ∣ ∣ ∣ F, J ˜ ˜ K K a ˜ c , n, M 〉 = 5 2 0 ∣ ∣ / , 00 , 1 1 , / 2 〉 for IB or | 0 0 , 00 , 1 0 , ⟩ for 1,3-DIB were analyzed.

the alignment than for the rigid-rotor, with the opposite effect observed for the main peak f ≈ 1 2 . GHz in 1,3-DIB, Fig. 3 d. The absolute value of the shift is similar for both molecules, however, relative to the peak frequency the shift in 1,3-DIB is about four times bigger than in IB. Thus, for 1,3-DIB, the deviations between the phases in the post-pulse alignment for the result with and without the quadrupole coupling occur on the timescale of the rotational revivals. In contrast for IB, the effect is slower than the molecule's longer rotational period. Furthermore, the larger number of hyperfine states in 1,3-DIB compared to IB led to an increased amount of irregularly spaced frequency contributions, giving rise to the non-periodic dynamics, see Fig. 3 b.

The Fourier transforms of the rigid-rotor results show only a few quasi-regularly-spaced peaks, reflecting the periodic behavior of 〈 cos 2 θ 〉 . For both molecules, the highest intensity peak corresponds to the coupling of the rigid-rotor states ∣ ∣ J K ,K a c 〉 = 0 | 00 , 0 ⟩ and | 2 02 , 0 ⟩ . When the quadrupole coupling is considered, the main rigidrotor peaks are split into several frequencies corresponding to transitions between different hyperfine states. For IB, Fig. 3 c, the main peak at f ≈ 4 2 . GHz is shifted toward higher frequency, which leads to a faster oscillation of

Other initial spin-rotational eigenstates of 1,3-DIB, which contribute to the alignment in Fig. 2 c, show analogous patterns, but with different frequency contributions. As a result, the revivals of 〈 cos 2 θ 〉 are on average almost entirely suppressed in 1,3-DIB. For both molecules, the frequency components below f ≈ 1 GHz originate from the energy gaps between different hyperfine states belonging to the same rotational level, i. e., to pure hyperfine-level population transfers similar to the transitions observed in NMR spectroscopy.

## 1. Influence of the number and strengths of nuclear spins

To gain more insight into the influence of the number of nuclear spins, we compared the alignment dynamics of 1,3DIB to that of a modified single-spin 1,3-DIB molecule, where we considered the quadrupole coupling of only one iodine nucleus with all other properties of 1,3-DIB unchanged. This allowed us to compare two molecular systems that differ only in the number of quadrupole-coupled nuclei. Furthermore, we investigated the alignment of 1,3-BIB, which can be considered an intermediate case between 1,3-DIB and the single-spin 1,3-DIB with respect to its rotational and hyperfine splittings.

The alignment of the single-spin 1,3-DIB, shown in Fig. 4 a, exhibits stronger or longer-lasting revivals than the regular 1,3-DIB with two coupled iodine nuclei, see Fig. 2 c.When comparing to a molecule with the same number of spins, IB, the impact of the quadrupole coupling in the single-spin 1,3-DIB is still more pronounced, Fig. 2 a, even though the quadrupole splittings for both molecules have similar values. However, the pure rotational energy spacings in IB are much larger than in 1,3-DIB. This highlights that it is the relative differences between the hyperfine and the pure rotational energy spacings in the molecule that determines the magnitude of the nuclearquadrupole effect in the alignment. This conclusion is further supported by the rotational dynamics of 1,2-DIB, which has two iodine nuclei, but much larger rotational energy splittings than 1,3-DIB. As a result, 〈 cos 2 θ 〉 is affected by the quadrupole interactions in a manner that qualitatively resembles the behavior of IB rather than 1,3-DIB. This can be seen, for example, when the peak alignment in 1,2-DIB increases after four nanoseconds

FIG. 4. Impulsive alignment induced by τ FWHM = 1 ps laser pulses with I 0 = 10 11 W cm / 2 for (a) 1,3-DIB neglecting the quadrupole coupling of the second iodine nucleus and (b) 1,3bromoiodobenzene including (QC) and neglecting the nuclear quadrupole coupling (noQC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-558a7b3e-b67a-4b00-8b26-dd1d2693ff9e.jpg)

in Fig. 2 d. The alignment for 1,3-BIB with and without quadrupole coupling is depicted in Fig. 4 b. This molecule has fewer hyperfine energy levels and slightly larger rotational-energy spacings than 1,3-DIB. As a result, the low-amplitude oscillations of the alignment have a more regular structure for 1,3-BIB, while the overall effect of the quadrupole coupling in 1,3-BIB is similar to that in 1,3-DIB, see Fig. 2 c.

## 2. Influence of the geometric symmetry

Apart from their rotational constants and the number of hyperfine levels, an important distinction of the molecules considered here lies in the symmetry of their EFG and polarizability tensors. In IB both tensors are diagonal in the coordinate system defined by the principal axes of inertia, while the EFG tensors of the other molecules and the polarizability tensor of DIBN have nonzero offdiagonal elements in the inertial frame. The diagonal elements of these tensors couple rotational states having the same symmetry in the D 2 rotation group. 31 The offdiagonal elements α ab and V ab couple rotational states with A ↔ B c symmetry as well as B a ↔ B b . For DIBN, and the two diiodobenzene molecules, however, a given spin-rotational eigenstate is a linear combination of basis states | F, J K K a c , I, M ⟩ with A and B c or B a and B b rotational symmetries. For an eigenstate ∣ ∣ ∣ F, J ˜ ˜ K K a ˜ c , n, M 〉 , only an approximate rotational symmetry can be assigned according to even or odd ˜ K a and ˜ K c . 31 Consequently, a weak laser-interaction coupling can be observed between the hyperfine states corresponding to rotational states which are not coupled by the laser field in the absence of the quadrupole interactions.

Without quadrupole interactions, the wave packets of IB, 1,3-DIB and 1,2-DIB consist exclusively of rotational states of the same symmetry in the D 2 group as the initial state, i. e., A symmetry for T rot = 0 K. For DIBN, states with B c symmetry also contribute due to the polarizability tensor's nonzero off-diagonal elements. When considering quadrupole coupling, nonzero populations of B c rotational states are also observed for 1,3-DIB and 1,2-DIB. These contributions are small, ≲ 5 % , for the initial states and laser intensity considered here and thus do not influence the alignment significantly. However, this aspect of the rotational dynamics is not captured if the nuclear-quadrupole interactions are neglected. In addition, for rotational levels which are strongly mixed by these interactions, a much larger impact on the population distribution is possible. This is, for instance, the case for the excited states | 2 12 , M J ⟩ and | 2 02 , M J ⟩ of 1,2-DIB, see Section III C.

## B. High-intensity impulsive alignment

Typical impulsive alignment experiments employ intenser laser pulses, 32-35 which excite a larger number of high-energy rotational states than discussed in the previous section. We investigated the strong field alignment using the same pulse shape as in the previous section but with the peak intensity increased tenfold, i. e., I 0 = 10 12 W cm . The temporal evolution of / 2 〈 cos 2 θ 〉 for all molecules except 1,2-DIB is characterized by J -type revivals 33,36 that appear at multiples of the rotational periods as well as integer fractions thereof, see Fig. 5. These beatings originate from the interference of highly excited rotational states with ∆ K a = 0 . The asymmetry splittings of these states prevent a complete rephasing of

FIG. 5. Impulsive alignment induced by a Gaussian laser pulse with τ FWHM = 1 ps and I 0 = 10 12 W cm / 2 for (a) IB, (b) DIBN, (c) 1,3-DIB and (d) 1,2-DIB including (QC) and neglecting the nuclear-quadrupole coupling (noQC).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-2d154bf6-55c1-4184-aeb6-15cd0f4530e4.jpg)

the wave packet and the peak alignment at the revivals decreases over time. 37

For 1,2-DIB, we observe clear revival features of a different type spaced by 304 ps in Fig. 5 d. The post-pulse dynamics of this molecule is dominated by the dephasing and rephasing of rotational states with ∆ J = ∆ K a = 2 and J ≈ K a . For large J , the energies of these states can be approximated by E J, K ( a ) ≈ BJ J ( +1)+( A -B K ) 2 a with B = 1 2 ( B + C ) . The relevant energy gaps are thus multiples of 4 A yielding the period of 1 4 / A ≈ 304 ps in 〈 cos 2 θ 〉 .

For the stronger laser pulse, the nuclear-quadrupole interactions had a much weaker impact on the post-pulse rotational dynamics. generally, it led to a loss of the peak alignment over time instead of qualitatively altering the dynamics, 28 see Fig. 5. An exception of this primarily destructive effect of the quadrupole coupling is found at certain alignment peaks of 1,2-DIB, e. g., t = 5 17 . ns, where 〈 cos 2 θ 〉 with quadrupole coupling is slightly larger than without. We attribute this constructive effect of the quadrupole coupling to the hyperfine splittings of some low-energy rotational states that cause a change in the amplitude of 〈 cos 2 θ 〉 similar to the weak laser intensity regime, Fig. 2 c. Overall, the weak impact of the quadrupole interactions is a result of increasingly uniform hyperfine splitting patterns at higher rotational excitations, as previously shown for I 2 and 1,4-DIB. 28

The FFTs of the post-pulse alignment for IB and 1,3DIB are shown in Fig. 6 c, d. The stronger-intensity pulse creates much broader wave packets that contain multiple frequency components. For both molecules the Fourier transforms calculated with and without quadrupole interactions are similar, especially for higher frequencies f that correspond to transitions between highly-excited rotational states. Deviations occur mainly in the lower frequency range, which contributes to a slight dephasing and the corresponding decrease of the alignment at longer timescales, noticed in Fig. 6 a, b. As for the weak laser intensity, the dephasing is more pronounced in 1,3-DIB due to the stronger mixing of less separated neighboring rotational states and generally larger number of the hyperfine components.

## C. Influence of the temperature

While single-state molecular ensembles can be achieved for small molecules, 35,38-41 this is generally not feasible for larger molecules. 41-43 We thus need to understand the effect of the nuclear-quadrupole coupling at finite rotational temperatures. We computed the post-pulse impulsive alignment of thermal ensembles of IB molecules for the two laser intensities used above and the rotational temperatures T rot = 0 1 . K and T rot = 0 3 . K. Such thermal ensembles are comparable 44 to state selected molecular ensembles 38,41,45 and represent conditions that are experimentally achievable.

The alignment of the thermal ensembles is shown in

FIG. 6. Impulsive alignment for I 0 = 10 12 W cm / 2 with (QC) and without (noQC) nuclear-quadrupole coupling for the initial states (a) ∣ ∣ 5 / , 2 0 00 , 1 2 / 〉 of IB and (b) | 0 0 , 00 , 0 ⟩ of 1,3-DIB. (c, d) Fast-Fourier transforms of 〈 cos 2 θ 〉 dynamics of the dynamics in (a) and (b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-07c8b707-526f-4bce-b31a-5c5c52ddcfd0.jpg)

Fig. 7. Without quadrupole interactions, 〈 cos 2 θ 〉 shows revival features with a decreased degree of alignment as compared to T rot = 0 K, which is due to the influence of 32 (130) initially populated excited states with up to J = 3 ( J = 5 ) for T rot = 0 1 . K (T rot = 0 3 . K). As for T rot = 0 K, the quadrupole-coupling effects are stronger at the weaker field intensity I 0 = 10 11 W cm , see Fig. 7 a, / 2 c. For this laser intensity, a periodic increase of the amplitude of 〈 cos 2 θ 〉 as observed for T rot = 0 K in Fig. 2 a is not present. At a stronger intensity I 0 = 10 12 W cm / 2 the effect of the quadrupole coupling is essentially negligible until the first full revival, see Fig. 7 b, d. Shortly thereafter the quadrupole interaction starts to manifest itself, reducing the overall degree of alignment.

I 0 = 10 12 W cm / 2 (Fig. 8 b) the effect is still more pronounced than for T rot = 0 K in Fig. 2 d. This can be attributed to the nuclear-quadrupole interactions strongly mixing the rotational states | 2 12 ⟩ and | 2 02 ⟩ , which are not coupled by the polarizability interaction. This mixing led to significant contributions of basis states | F, 2 02 , I, M ⟩ in the initial wave functions ∣ ∣ ∣ F, ˜ 2 ˜˜ 12 , n, M 〉 modifying the dynamics even before the laser field was applied. After the interaction with the laser pulse, we find significant populations of hyperfine states belonging to rotational states that are otherwise not excited when the quadrupole coupling is neglected.

To account for the nuclear-quadrupole effect in thermal ensembles of different molecular species, we additionally computed the alignment dynamics of the excited rotational state | 2 12 ⟩ of 1,2-DIB, see Fig. 8. We compared the average rigid-rotor alignment over all possible M J values of the initial state to the average result obtained for all initial hyperfine states ∣ ∣ ∣ F, ˜ 2 ˜˜ 12 , n, M 〉 . As for the rotational ground state, the effect of the quadrupole coupling is stronger for weaker intensity laser field. For

We thus conclude that for finite-temperature molecular beams the influence of the nuclear-quadrupole interaction generally decreases with increasing the laser intensity. This is further strengthened by our previous analysis of the excited state rotational dynamics of I 2 and 1,4-DIB. 28 However, the hyperfine interactions can have a significant impact on the dynamics of some excited initial states and cannot be neglected even on a time scale where they do not affect the alignment for T rot = 0 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-2665a6fd-4dfe-4f2f-b577-b5a318368dff.jpg)

FIG. 7. Impulsive alignment of a thermal ensemble of IB molecules induced by laser pulses with τ FWHM = 1 ps including and neglecting the quadrupole coupling (QC and noQC) for (a, b) T rot = 0 1 . K and (c, d) T rot = 0 3 . K. (a, c) show the results for I 0 = 10 11 W cm / 2 and (b, d) for I 0 = 10 12 W cm . / 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01125v1-figure-06fa3340-686d-44fb-9c48-d86184d11081.jpg)

strong nuclear-quadrupole interactions. We focused on differences in the hyperfine interactions due to different chemical structures of the molecules. We showed that the interplay of hyperfine energy shifts and rotational energy level spacings is a deciding factor for the extent to which the coupling affects the post-pulse alignment. For all molecules analyzed in the present work and in previous studies 27,28 we observed that the influence of the quadrupole coupling becomes small when the field-dressed dynamics is dominated by excited rotational states. Due to the variety of molecules considered, we expect that this observation is general and that pushing population into highly-excited rotational states is advantageous for minimizing the dephasing due to nuclear-quadrupole coupling.

FIG. 8. Impulsive alignment averaged over the initial states | 2 12 , M J ⟩ of 12-DIB induced by 1 ps laser pulses with (a) I 0 = 10 11 W cm , (b) / 2 I 0 = 10 12 W cm . The alignment is / 2 shown including (QC) and neglecting the nuclear quadrupole coupling (noQC).

## IV. SUMMARY AND CONCLUSIONS

We presented a computational analysis of the post-pulse alignment dynamics of asymmetric top molecules with

We point out that while nuclear-quadrupole interactions pose difficulties in achieving strong field-free alignment, their effect on the rotational dynamics might be exploited to probe nuclear or electronic excitations which alter the nuclear spins and quadrupole momenta as well as electric field gradients. The sensitivity of the alignment to the hyperfine energy level structures of low-energy rotational states may thus be viewed as an opportunity rather than a disadvantage.

## V. ACKNOWLEDGEMENTS

This work was supported by Deutsches ElektronenSynchtrotron DESY, a member of the Helmholtz Association (HGF), including the maxwell compute-cluster operated by DESY, by the Deutsche Forschungsgemeinschaft (DFG) through the priority program 'Quantum Dynamics in Tailored Intense Fields' (QUTIF, SPP1840, KU 1527/3, YA 610/1), and through the clusters of excellence 'Center for Ultrafast Imaging' (CUI, EXC 1074, ID 194651731) and 'Advanced Imaging of Matter' (AIM, EXC 2056, ID 390715994). R.G.F. gratefully acknowledges financial support by the Spanish Project PID2020113390GB-I00 (MICIN), and by the Andalusian research group FQM-207.

## Appendix A: Molecular parameters

The rotational constants, asymmetry parameters κ = (2 B -A -C / A ) ( -C ) and elements of the polarizability and EFG tensors of the iodine and bromine nuclei in the principle axes of inertia frame are listed in Table I. Except for IB, the rotational constants were computed from the equilibrium geometry obtained using density functional theory (DFT) with the B3LYP functional and the def2-QZVPP basis set. 47,48 The effective core potential def2-ECP 49 was used for the iodine atoms. The electric field gradient and polarizability tensors of DIBN, 1,3-DIB and 1,2-DIB were calculated at the DFT/B3LYP level of theory using the all-electron scalar relativistic DouglasKroll-Hess Hamiltonian 50 with the DKH-def2-TZVP basis set. 51,52 All electronic structure calculations were carried out with the quantum-chemistry package ORCA. 53,54 For IB, we used the experimental nuclear-quadrupole coupling constants 46 χ aa = -1892 039 . MHz, χ bb = 978 816 . MHz, χ cc = 913 222 . MHz instead of a calculated electric field gradient and quadrupole moment. The molecule-fixed frame ( x, y, z ) is defined by the principle axes of inertia so that the z axis is parallel to the most polarizable axis or, in the case of DIBN, forms the smallest possible angle ( ≈ 5 ◦ ) with it. The z axis is thus chosen parallel to the a axis for IB, DIBN, 1,3-DIB and 1,3-bromoiodobenzene and to the b axis for 1,2-DIB.

## Appendix B: Matrix representation of the Hamiltonian

To obtain the eigenstates of H mol , we first solved the time-independent Schrödinger equation for H rot yielding the rigid-rotor energy levels and eigenstates | J K K a c , M J ⟩ . Here, M J is the eigenvalue of ˆ J Z and K a,c are associated with the projections ˆ J a,c and are good quantum numbers in the prolate and oblate symmetric top cases, respectively. The matrix representation of H mol was then constructed in the coupled basis 55,56 | F, J K K a c , I, M ⟩ , M being the eigenvalue of ˆ F Z . The field-free energies and eigenstates

∣ ∣ ∣ F, J ˜ ˜ K K a ˜ c , n, M 〉 were computed by numerically solving the associated time-independent Schrödinger equation. Here, n = 1 2 , , ... is an index labeling states with the same F and the approximate quantum numbers ˜ J , ˜ K a and ˜ K c were assigned according to the coupled basis state | F, J K K a c , I, M ⟩ with the largest contribution to the given eigenstate. For IB, the nuclear spin quantum number I = I 1 remains a good quantum number. The matrix elements of the polarizability element α ZZ are known analytically 55 in the basis | F, J K K a c , I, M ⟩ and can be transformed easily into the field-free eigenbasis. 27

- 1 R. Altkorn, R. N. Zare, and C. H. Greene, 'Depolarization of optically prepared molecules by two randomly oriented spins,' Mol. Phys. 55 , 1-9 (1985).
- 2 C. Yan and A. C. Kummel, 'Effect of hyperfine depolarization upon creation and detection of alignment in free-jet expansions via selective photodissociation,' J. Chem. Phys. 98 , 6869-6882 (1993).
- 3 T. A. Cool and N. Hemmi, 'Hyperfine polarization quantum beats in cyanogen,' J. Chem. Phys. 103 , 3357-3366 (1995).
- 4
- 6 D. Sofikitis, L. Rubio-Lago, M. R. Martin, D. J. Ankeny Brown, N. C.-M. Bartlett, A. J. Alexander, R. N. Zare, and T. P. Rakitzis, 'Optical control of ground-state atomic orbital alignment: Cl ( 2 P 3 / 2 ) atoms from HCl ( ν = 2 , j = 1) photodissociation,' J. Chem. Phys. 127 , 144307 (2007).
- J. Zhang, C. W. Riehn, M. Dulligan, and C. Wittig, 'An experimental study of HF photodissociation: Spin-orbit branching ratio and infrared alignment,' J. Chem. Phys. 104 , 7027-7035 (1996). 5 E. R. Wouters, L. D. Siebbeles, K. L. Reid, B. Buijsse, and W. J. van der Zande, 'Observation of fine structure and hyperfine structure depolarization in the photofragment anisotropy in triplet H ,' Chem. Phys. 2 218 , 309-323 (1997).
- 7 N. C.-M. Bartlett, D. J. Miller, R. N. Zare, A. J. Alexander, D. Sofikitis, and T. P. Rakitzis, 'Time-dependent depolarization of aligned hd molecules,' Phys. Chem. Chem. Phys. 11 , 142-147 (2009).
- 8 N. C.-M. Bartlett, J. Jankunas, R. N. Zare, and J. A. Harrison, 'Time-dependent depolarization of aligned D 2 caused by hyperfine coupling,' Phys. Chem. Chem. Phys. 12 , 15689-15694 (2010).
- 9 K. Grygoryeva, J. Rakovský, O. Votava, and M. Fárník, 'Imaging of rotational wave-function in photodissociation of rovibrationally excited HCl molecules,' J. Chem. Phys. 147 , 013901 (2017).
- 10 J. C. H. Spence and H. N. Chapman, 'The birth of a new field,' Phil. Trans. R. Soc. B 369 , 20130309-20130309 (2014).
- 11 J. Küpper, S. Stern, L. Holmegaard, F. Filsinger, A. Rouzée, A. Rudenko, P. Johnsson, A. V. Martin, M. Adolph, A. Aquila, S. Bajt, A. Barty, C. Bostedt, J. Bozek, C. Caleman, R. Coffee, N. Coppola, T. Delmas, S. Epp, B. Erk, L. Foucar, T. Gorkhover, L. Gumprecht, A. Hartmann, R. Hartmann, G. Hauser, P. Holl, A. Hömke, N. Kimmel, F. Krasniqi, K.-U. Kühnel, J. Maurer, M. Messerschmidt, R. Moshammer, C. Reich, B. Rudek, R. Santra, I. Schlichting, C. Schmidt, S. Schorb, J. Schulz, H. Soltau, J. C. H. Spence, D. Starodub, L. Strüder, J. Thøgersen, M. J. J. Vrakking, G. Weidenspointner, T. A. White, C. Wunderer, G. Meijer, J. Ullrich, H. Stapelfeldt, D. Rolles, and H. N. Chapman, 'X-ray diffraction from isolated and strongly aligned gas-phase molecules with a free-electron laser,' Phys. Rev. Lett. 112 , 083002 (2014), arXiv:1307.4577 [physics].
- 12 J. Yang, X. Zhu, T. J. A. Wolf, Z. Li, J. P. F. Nunes, R. Coffee, J. P. Cryan, M. Gühr, K. Hegazy, T. F. Heinz, K. Jobe, R. Li, X. Shen, T. Veccione, S. Weathersby, K. J. Wilkin, C. Yoneda, Q. Zheng, T. J. Martínez, M. Centurion, and X. Wang, 'Imaging CF I conical intersection and photodissociation dynamics with 3 ultrafast electron diffraction,' Science 361 , 64-67 (2018).
- 13 M. S. Hunter, C. H. Yoon, H. DeMirci, R. G. Sierra, E. H. Dao, R. Ahmadi, F. Aksit, A. L. Aquila, H. Ciftci, S. Guillet, M. J. Hayes, T. J. Lane, M. Liang, U. Lundström, J. E. Koglin,

TABLE I. Rotational constants A B C , , , asymmetry parameter κ , elements of the polarizability tensors α ij and electric field gradient tensors V ij ( ) l of the l th nucleus in the coordinate system defined by the principle axes of inertia. In the case of 1,3-BIB, V ij (1) and V ij (2) refer to the EFG of the iodine and bromine nucleus, respectively. For IB, the rotational constants from ref. 46 and the polarizability tensor from ref. 36 were used.

|                 | IB a             | DIBN       | 1,3-DIB    | 1,2-DIB    | 1,3-BIB    |
|-----------------|------------------|------------|------------|------------|------------|
| A (MHz)         | 5669 . 131       | 1608 . 324 | 1859 . 912 | 822 . 0850 | 2040 . 244 |
| B (MHz)         | 750 . 414 293    | 155 . 402  | 210 . 776  | 522 . 1492 | 281 . 807  |
| C (MHz)         | 662 . 636 146    | 141 . 709  | 189 . 321  | 319 . 3276 | 247 . 607  |
| κ               | - 0 . 965        | - 0 . 981  | - 0 . 974  | - 0 . 193  | - 0 . 962  |
| α aa (Å 3 )     | 21 . 5           | 34 . 506   | 28 . 847   | 22 . 072   | 25 . 774   |
| α bb (Å 3 )     | 15 . 3           | 21 . 662   | 19 . 789   | 24 . 830   | 18 . 234   |
| α cc (Å 3 )     | 10 . 2           | 12 . 123   | 11 . 292   | 11 . 103   | 10 . 192   |
| α ab (Å 3 )     | 0                | - 1 . 209  | 0          | 0          | - 0 . 419  |
| V aa (1) (a.u.) | - 1892 . 039 /eQ | 12 . 327   | 7 . 34     | - 0 . 30   | 8 . 3735   |
| V bb (1) (a.u.) | 978 . 816 /eQ    | - 6 . 6802 | - 1 . 73   | 5 . 88     | - 2 . 7590 |
| V cc (1) (a.u.) | 913 . 222 /eQ    | - 5 . 6472 | - 5 . 61   | - 5 . 57   | - 5 . 6146 |
| V ab (1) (a.u.) | 0                | 1 . 5092   | - 7 . 85   | 8 . 84     | - 7 . 1745 |
| V aa (2) (a.u.) | 0                | 12 . 155   | 7 . 34     | - 0 . 30   | 4 . 1656   |
| V bb (2) (a.u.) | 0                | - 6 . 4317 | - 1 . 73   | 5 . 88     | - 0 . 4784 |
| V cc (2) (a.u.) | 0                | - 5 . 7230 | - 5 . 61   | - 5 . 57   | - 3 . 6872 |
| V ab (2) (a.u.) | 0                | 0 . 570 71 | 7 . 85     | - 8 . 84   | 5 . 6800   |

- a In the calculations for IB we used experimentally determined quadrupole coupling constants [46, Table 6].
- P. Mgbam, Y. Rao, L. Zhang, S. Wakatsuki, J. M. Holton, and S. Boutet, 'Selenium single-wavelength anomalous diffraction de novo phasing using an x-ray-free electron laser,' Nat. Commun. 7 , 13388 (2016).
- 14 A. Barty, J. Küpper, and H. N. Chapman, 'Molecular imaging using x-ray free-electron lasers,' Annu. Rev. Phys. Chem. 64 , 415-435 (2013).
- 'Switched wave packets with spectrally truncated chirped pulses,' J. Chem. Phys. 148 , 221105 (2018), arXiv:1803.03953 [physics]. 26 E. F. Thomas, A. A. Søndergaard, B. Shepperson, N. E. Henriksen, and H. Stapelfeldt, 'Hyperfine-structure-induced depolarization of impulsively aligned I 2 molecules,' Phys. Rev. Lett. 120 , 163202 (2018), arXiv:1804.04416 [physics].
- 15 L. Christensen, J. H. Nielsen, C. B. Brandt, C. B. Madsen, L. B. Madsen, C. S. Slater, A. Lauer, M. Brouard, M. P. Johansson, B. Shepperson, and H. Stapelfeldt, 'Dynamic stark control of torsional motion by a pair of laser pulses,' Phys. Rev. Lett. 113 , 073005 (2014).
- 16 R. J. D. Miller, 'Mapping atomic motions with ultrabright electrons: The chemists' gedanken experiment enters the lab frame,' Annu. Rev. Phys. Chem. 65 , 583-604 (2014).
- 17 H. Stapelfeldt and T. Seideman, 'Colloquium: Aligning molecules with strong laser pulses,' Rev. Mod. Phys. 75 , 543-557 (2003).
- 18 T. Seideman, 'Revival structure of aligned rotational wave packets,' Phys. Rev. Lett. 83 , 4971-4974 (1999).
- 19 F. Rosca-Pruna and M. J. J. Vrakking, 'Revival structures in picosecond laser-induced alignemnt of I 2 molecules. II. Numerical modeling,' J. Chem. Phys. 116 , 6579 (2002).
- 20 E. Hamilton, T. Seideman, T. Ejdrup, M. D. Poulsen, C. Z. Bisgaard, S. S. Viftrup, and H. Stapelfeldt, 'Alignment of symmetric top molecules by short laser pulses,' Phys. Rev. A 72 , 043402 (2005).
- 21 T. Seideman, 'On the dynamics of rotationally broad, spatially aligned wave packets,' J. Chem. Phys. 115 , 5965 (2001).
- 22 J. Underwood, M. Spanner, M. Ivanov, J. Mottershead, B. Sussman, and A. Stolow, 'Switched wave packets: A route to nonperturbative quantum control,' Phys. Rev. Lett. 90 , 223001 (2003).
- 23 J. Underwood, B. Sussman, and A. Stolow, 'Field-free three dimensional molecular axis alignment,' Phys. Rev. Lett. 94 , 143002 (2005).
- 24 A. Goban, S. Minemoto, and H. Sakai, 'Laser-field-free molecular orientation,' Phys. Rev. Lett. 101 , 013001 (2008).
- 25 A. Chatterley, E. T. Karamatskos, C. Schouder, L. Christiansen, A. V. Jörgensen, T. Mullins, J. Küpper, and H. Stapelfeldt,
- 27 A. Yachmenev, L. V. Thesing, and J. Küpper, 'Laser-induced dynamics of molecules with strong nuclear quadrupole coupling,' J. Chem. Phys. 151 , 244118 (2019), arXiv:1910.13275 [physics].
- 28 L. V. Thesing, A. Yachmenev, R. González-Férez, and J. Küpper, 'The effect of nuclear-quadrupole coupling in the laser-induced alignment of molecules,' J. Phys. Chem. A 124 , 2225-2230 (2020), arXiv:1911.12270 [physics].
- 29 P. Pyykkö, 'Year-2008 nuclear quadrupole moments,' Mol. Phys. 106 , 1965-1974 (2008).
- 30 A. Owens and A. Yachmenev, 'RichMol: A general variational approach for rovibrational molecular dynamics in external electric fields,' J. Chem. Phys. 148 , 124102 (2018), arXiv:1802.07603 [physics].
- 31 W. Gordy and R. L. Cook, Microwave Molecular Spectra , 3rd ed. (John Wiley &amp; Sons, New York, NY, USA, 1984).
- 32 F. Rosca-Pruna and M. J. J. Vrakking, 'Experimental observation of revival structures in picosecond laser-induced alignment of I 2 ,' Phys. Rev. Lett. 87 , 153902 (2001).
- 33 L. Holmegaard, S. S. Viftrup, V. Kumarappan, C. Z. Bisgaard, H. Stapelfeldt, E. Hamilton, and T. Seideman, 'Control of rotational wave-packet dynamics in asymmetric top molecules,' Phys. Rev. A 75 , 051403 (2007).
- 34 S. Trippel, T. Mullins, N. L. M. Müller, J. S. Kienitz, J. J. Omiste, H. Stapelfeldt, R. González-Férez, and J. Küpper, 'Strongly driven quantum pendulum of the carbonyl sulfide molecule,' Phys. Rev. A 89 , 051401(R) (2014), arXiv:1401.6897 [quant-ph].
- 35 E. T. Karamatskos, S. Raabe, T. Mullins, A. Trabattoni, P. Stammer, G. Goldsztejn, R. R. Johansen, K. Długołęcki, H. Stapelfeldt, M. J. J. Vrakking, S. Trippel, A. Rouzée, and J. Küpper, 'Molecular movie of ultrafast coherent rotational dynamics of OCS,' Nat. Commun. 10 , 3364 (2019), arXiv:1807.01034 [physics].

- 36 M. D. Poulsen, E. Peronne, H. Stapelfeldt, C. Z. Bisgaard, S. Viftrup, E. Hamilton, and T. Seideman, 'Nonadiabatic alignment of asymmetric top molecules: Rotational revivals,' J. Chem. Phys. 121 , 783-791 (2004).
- 37 A. Rouzee, S. Guerin, V. Boudon, B. Lavorel, and O. Faucher, 'Field-free one-dimensional alignment of ethylene molecule,' Phys. Rev. A 73 , 033418-9 (2006).
- 38 Y.-P. Chang, D. A. Horke, S. Trippel, and J. Küpper, 'Spatiallycontrolled complex molecules and their applications,' Int. Rev. Phys. Chem. 34 , 557-590 (2015), arXiv:1505.05632 [physics].
- 39 J. H. Nielsen, P. Simesen, C. Z. Bisgaard, H. Stapelfeldt, F. Filsinger, B. Friedrich, G. Meijer, and J. Küpper, 'Starkselected beam of ground-state OCS molecules characterized by revivals of impulsive alignment,' Phys. Chem. Chem. Phys. 13 , 18971-18975 (2011), arXiv:1105.2413 [physics].
- 40 D. A. Horke, Y.-P. Chang, K. Długołęcki, and J. Küpper, 'Separating para and ortho water,' Angew. Chem. Int. Ed. 53 , 1196511968 (2014), arXiv:1407.2056 [physics].
- 41 S. Y. T. van de Meerakker, H. L. Bethlem, N. Vanhaecke, and G. Meijer, 'Manipulation and control of molecular beams,' Chem. Rev. 112 , 4828-4878 (2012).
- 42 K. Wohlfart, F. Grätz, F. Filsinger, H. Haak, G. Meijer, and J. Küpper, 'Alternating-gradient focusing and deceleration of large molecules,' Phys. Rev. A 77 , 031404(R) (2008), arXiv:0803.0650 [physics].
- 43 H. L. Bethlem, M. R. Tarbutt, J. Küpper, D. Carty, K. Wohlfart, E. A. Hinds, and G. Meijer, 'Alternating gradient focusing and deceleration of polar molecules,' J. Phys. B 39 , R263-R291 (2006), arXiv:0604020 [physics].
- 44 L. V. Thesing, J. Küpper, and R. González-Férez, 'Timedependent analysis of the mixed-field orientation of molecules without rotational symmetry,' J. Chem. Phys. 146 , 244304 (2017), arXiv:1705.03225 [physics].
- 45 F. Filsinger, G. Meijer, H. Stapelfeldt, H. N. Chapman, and J. Küpper, 'State- and conformer-selected beams of aligned and oriented molecules for ultrafast diffraction studies,' Phys. Chem. Chem. Phys. 13 , 2076-2087 (2011), arXiv:1009.0871 [physics].
- 46 J. L. Neill, S. T. Shipman, L. Alvarez-Valtierra, A. Lesarri, Z. Kisiel, and B. H. Pate, 'Rotational spectroscopy of iodobenzene and iodobenzene-neon with a direct digital 2-8ghz chirped-pulse fourier transform microwave spectrometer,' J. Mol. Spectrosc. 269 , 21 - 29 (2011).
- 47 F. Weigend, F. Furche, and R. Ahlrichs, 'Gaussian basis sets of quadruple zeta valence quality for atoms H-Kr,' J. Chem. Phys. 119 , 12753-12762 (2003).
- 48 F. Weigend and R. Ahlrichs, 'Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy,' Phys. Chem. Chem. Phys. 7 , 3297 (2005).
- 49 K. A. Peterson, D. Figgen, E. Goll, H. Stoll, and M. Dolg, 'Systematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16-18 elements,' J. Chem. Phys. 119 , 11113-11123 (2003).
- 50 F. Neese, A. Wolf, T. Fleig, M. Reiher, and B. A. Hess, 'Calculation of electric-field gradients based on higher-order generalized Douglas-Kroll transformations,' J. Chem. Phys. 122 , 204107 (2005).
- 51 F. E. Jorge, A. C. Neto, G. G. Camiletti, and S. F. Machado, 'Contracted gaussian basis sets for Douglas-Kroll-Hess calculations: Estimating scalar relativistic effects of some atomic and molecular properties,' J. Chem. Phys. 130 , 064108 (2009).
- 52 C. Campos and F. Jorge, 'Triple zeta quality basis sets for atoms Rb through Xe: application in CCSD(T) atomic and molecular property calculations,' Mol. Phys. 111 , 167-173 (2012).
- 53 F. Neese, 'The ORCA program system,' Wiley Interdiscip. Rev. Comput. Mol. Sci. 2 , 73-78 (2011).
- 54 F. Neese, 'Software update: the ORCA program system, version 4.0,' Wiley Interdiscip. Rev. Comput. Mol. Sci. 8 , e1327 (2017).
- 55 R. L. Cook and F. C. de Lucia, 'Application of the theory of irreducible tensor operators to molecular hyperfine structure,' Am. J. Phys. 39 , 1433-1454 (1971).
- 56 A. Yachmenev and J. Küpper, 'Communication: General variational approach to nuclear-quadrupole coupling in rovibrational spectra of polyatomic molecules,' J. Chem. Phys. 147 , 141101 (2017), arXiv:1709.08558 [physics].
- 57 P. M. Felker, 'Rotational coherence spectroscopy: studies of the geometries of large gas-phase species by picosecond time-domain methods,' J. Phys. Chem. 96 , 7844-7857 (1992).
- 58 P. M. Felker, J. S. Baskin, and A. H. Zewail, 'Rephasing of collisionless molecular rotational coherence in large molecules,' J. Phys. Chem. 90 , 724-728 (1986).
- 59 C. Riehn, 'High-resolution pump-probe rotational coherence spectroscopy - rotational constants and structure of ground and electronically excited states of large molecular systems,' Chem. Phys. 283 , 297-329 (2002).
Linda V. Thesing Andrey Yachmenev Rosario González-Férez Jochen Küpper
Molecular influence on nuclear-quadrupole-coupling effects in laser induced alignment
We studied the effect of nuclear-quadrupole interactions on the field-free
impulsive alignment of different asymmetric-top molecules. Our analysis is
focused on the influence of the hyperfine- and rotational-energy-level
structures. These depend on the number of nuclear spins, the rotational
constants, and the symmetry of the tensors involved in the nuclear spin and
external field interactions. Comparing the prototypical large-nuclear-spin
molecules iodobenzene, 1,2-diiodobenzene, 1,3-diiodobenzene, and
2,5-diiodobenzonitrile, we demonstrate that the magnitude of the hyperfine
splittings compared to the rotational-energy splittings plays a crucial role in
the spin-rotational dynamics after the laser pulse. Moreover, we point out that
the impact of the quadrupole coupling on the rotational dynamics decreases when
highly excited rotational states dominate the dynamics.