##
1Collisions of excited Na atoms with N_{2}, C_{2}H_{2} and N_{2}O molecules. Surface-hopping calculation of the Na(3p) fine structure population
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj1

1

Collisions of excited Na atoms with N

_{2}, C_{2}H_{2}and N_{2}O molecules. Surface-hopping calculation of the Na(3p) fine structure population
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj1

2

The results of optical collision experiments: fine structure population ratios and polarization alignment angles obtained from scattering of Na atoms which are optically 3s → 3p excited during the collision with N

_{2}, C_{2}H_{2}and N_{2}O molecules, have been simulated numerically and analysed using quantum chemically computed potential energy surfaces, classical trajectory calculations and the semiclassical surface hopping method combined with a specially developed diabatization procedure.
Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj2

## Introduction

3

Optical collisions are scattering processes where the projectile is optically excited during its interaction with the target.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac1

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac2

5

The experiments can briefly be explained as follows.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac2

6

The collisions take place in a crossed beam arrangement.

Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met1

7

The collision complexes are excited with polarized laser light and a detection laser is used to separate the excited Na atoms from the background.

Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met2

8

The excitation laser is detuned by a few 100 cm

^{−1}from the resonant 3s → 3p transition of Na, as a rule, but not always to higher frequency.
Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp1

9

Therefore optical transitions occur only at those geometries, the so-called Condon configurations, where the energy difference in the collision complex corresponds just to the detuned laser frequency.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac3

10

Excitation probability as a function of detuning and polarization angle is recorded for a few selected scattering angles.

Type: Experiment |
Advantage: None |
Novelty: None |
ConceptID: Exp2

11

From flight time measurement the velocity of the Na atoms after collision is obtained; this allows to reconstruct the kinematics of the observed scattering processes.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met3

12

The goal of the experiments is to test the interaction potentials of excited Na with the target molecules and to obtain “snapshots” of the collision complexes, in particular of the direction of the 3p orbital after the process of optical excitation.

Type: Goal |
Advantage: None |
Novelty: None |
ConceptID: Goa1

13

The results are presented as polar diagrams of the probability as a function of an orientation angle; they are analysed with the aid of numerical simulations using quantum chemically calculated potential surfaces and classical trajectories.

^{2,3,6}
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met4

14

A related experiment consists in measuring the ratio of excited Na atoms leaving the complex in the lower 3p

_{1/2}fine structure state compared to the sum 3p_{1/2}+ 3p_{3/2}, the two levels being separated by spin–orbit interaction.^{3,5,7}
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met5

15

Under the conditions the experiments have been carried out, the non-zero results for this ratio can be explained only by assuming non-adiabatic coupling effects between different excited electronic states.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod1

16

Therefore precise knowledge of the interaction between the potentials is needed for a detailed understanding of the underlying process.

Type: Motivation |
Advantage: None |
Novelty: None |
ConceptID: Mot1

17

It is the aim of the present work to explain the observed fine structure ratios as functions of the kinetic energy after the collision numerically by combining classical trajectory methods with a semiclassical surface hopping procedure.

Type: Goal |
Advantage: None |
Novelty: None |
ConceptID: Goa2

## Optical collisions

18

Upon collision of an excited Na(3p) atom with a target molecule the degeneracy of the atomic 3p levels is removed to an extent depending on the separation and arrangement of the collision complex.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac4

19

In the case of a linear closed-shell target molecule as treated here three excited potential energy surfaces can be distinguished, two of

^{2}A′ (electronic) symmetry and one^{2}A″.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2

20

The 3p charge distribution is no longer spherically symmetric as in the unperturbed atom.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2

21

In each of the excited components the 3p orbital is described by a real wave function and accordingly the excited Na atom has a quadrupole moment (of about −16.8 au).

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2

22

The interaction of Na(3p) with the target is determined by the interaction of this quadrupole with the molecular multipole moments as well as by the dispersion effect.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2

23

Quantum chemically the latter is an effect of electron correlation whereas the multipole interactions can be treated already on the Hartree–Fock (HF) level.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod2

24

It is important to have the correct orientation of the 3p atomic orbitals already in low order of the theoretical treatment.

Type: Motivation |
Advantage: None |
Novelty: None |
ConceptID: Mot2

25

In our calculations we therefore choose a slightly distorted (fixed) geometry for the target molecules such that the HF multipoles correctly reproduce the experimental moments.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met6

26

Alternatively the molecular moments obtained from a correlated wave function could be used; we hope that the experimental moments also account for the average over the configurations of the rotating and vibrating target molecule.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met6

27

In a previous paper

^{6}we have explained how a qualitative understanding of the excited potential surfaces can be obtained if the most important contributions to the multipole interactions are discussed.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac5

28

For every geometrical arrangement of the collision complex each one of the three 3p surfaces is characterized by the spatial orientation of the corresponding 3p orbital.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod3

29

Whereas in the

^{2}A″ state the 3p atomic orbital is perpendicular to the plane defined by the linear target and the Na atom, the two^{2}A′ states have in-plane 3p orbitals.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod3

30

The expectation value of the 3p orbital direction after excitation leads to the so-called alignment angle, defined with respect to the center-of-mass (CM) trajectory of the incoming particle.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod3

31

With quantum chemistry the 3p direction angles can easily be obtained by calculating the 3s → 3p transition moment vectors; except for very close separations they are mutually orthogonal.

Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met7

32

For the descriptions of the experiments we also need the potential energy surface of the 3s ground state.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod4

33

In this state the Na atom is approximately spherical; the surface is nearly structureless down to rather short distances where the Pauli repulsion sets in.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod4

34

Furthermore on the three excited surfaces the spin–orbit interaction of the Na(3p) electron must be taken into account which contributes small corrections to the energies determined by the interaction of Na(3p) and the target molecule, in the order of the spin–orbit coupling constant of the Na(3p) state (11.5 cm

^{−1}).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod4

35

Thus numerical simulation of the scattering experiments is carried out by running classical trajectories on the surfaces calculated with

*ab-initio*quantum chemistry including electron correlation and the spin–orbit corrections.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met8

36

The polarization of the laser radiation is taken into account by projecting onto the 3s → 3p transition dipole vector.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met8

37

For this projection only the orientation of the transition vector is needed along the trajectories, not its value.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met8

38

Therefore our method to obtain the correct multipole interaction on the HF level is very useful and saving computer time.

Type: Method |
Advantage: Yes |
Novelty: New |
ConceptID: Met8

39

The three potential surfaces derived from the Na(3p) components may penetrate each other at certain geometrical arrangements,

*i.e.*change their energetic ordering.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5

40

This situation favours non-adiabatic couplings.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5

41

Penetration of the two surfaces of

^{2}A′ symmetry leads to conical intersections.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5

42

The

^{2}A′ surfaces of the Na(3p) + HCCH system have such intersection points on the molecular axis, about 8*a*_{0}off the position of the H atoms.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5

43

In the case of the N

_{2}, N_{2}O (and CO_{2}) target molecules, whose quadrupole moments are negative in contrast to HCCH, the penetration occurs on rings around the molecular axis with radii between 10*a*_{0}(CO_{2}) and 16.5*a*_{0}(N_{2}).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod5

44

Evidently the probability that a Na trajectory passes close to an intersection point such that non-adiabatic processes can be expected is much higher for a ring-shaped intersection line than in the case of only two axial intersection points.

Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con1

## Theory

45

In this section we shall investigate the dynamics of the collision process taking place on the quantum chemically calculated potential energy surfaces.

Type: Object |
Advantage: None |
Novelty: New |
ConceptID: Obj3

46

The dynamics of nuclear motion is treated according to the concept of classical trajectories.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod6

47

The relative velocity of the colliding particles at large distance and the rotational temperature of the molecule are specified by the experimental setup.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met9

48

The impact parameter

*b*and the initial rotational state of the molecule are determined by a random number generator for every trajectory.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met10

49

The bond distances and angles of the target molecule are kept fixed at linear geometry, thus the surfaces become two dimensional and only five coordinates remain describing the dynamics of the system in the center of mass system.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met11

50

At a distance depending on the detuning of the laser frequency the collision complex is excited to one of the states of the 3p-manifold.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met12

### Quantum chemical calculations

51

The potential energy and transition dipole moment surfaces of the collision complexes Na + N

_{2}and Na + C_{2}H_{2}have been presented previously.^{2,6}
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac6

52

Graphical representations of their transition vectors to the uppermost surfaces and of the Condon configurations for several detuning frequencies are shown below (Figs. 1 and 2).

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7

53

The value of the transition moments is practically constant for all transitions and every configuration.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7

54

The direction of these vectors characterizes the surfaces and is important for the analysis.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7

55

The sudden change of orientation in the figures indicates the position of the conical intersection where two surfaces penetrate each other.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod7

56

The quadrupole moments of the molecules were computed in HF approximation leading to the following results: 5.44 au for C

_{2}H_{2}and −0.93 au for N_{2}.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met13

57

The quantum chemical calculations for the system Na + N

_{2}O have been performed in a similar way.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met14

58

The surfaces are functions of only two geometrical parameters, the center of mass distance

*R*between the atom and the target molecule and the angle*ϑ*between*R*and the molecular axis.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod8

59

The fixed bond lengths are 2.18

*a*_{0}for NN and 2.17*a*_{0}for NO.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod9

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod9

61

0.167 D) and a quadrupole moment of −2.85 au.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod9

62

For the calculations a 7s7p4d3f ANO (atomic natural orbitals) basis set from Widmark

*et al.*has been used^{10}for the target molecule and a selfmade 8s6p4d3f basis for the Na atom.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod10

63

The dipole vector surfaces have been calculated in HF approximation; the uppermost dipole vector surface is shown in Fig. 3.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod11

64

For the interaction energies we used the MCCEPA (multiconfiguration coupled electron pair approximation) program of Fink and Staemmler

^{11}and single determinant (HF) reference functions.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met15

65

In order to justify the approximation for the orientation of the transition dipoles as introduced in section 2 we have carried out test calculations for the Na + N

_{2}and Na + N_{2}O systems, evaluating the transition dipole vectors with correlated wave functions.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met16

66

It turned out that although the values of the moments may change by up to 10 or 20% (mainly due to the reduced weight of the reference configurations in the correlated function) the orientation angles remain constant within 1° even for geometries close to a conical intersection where the electronic states are considerably mixed.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res1

67

The potential and transition dipole surfaces were interpolated by two-dimensional spline functions.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met17

68

The derivatives of the potential resulted from local expansions for every grid point and subsequent interpolation by spline functions.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met18

69

We found the achieved smoothness sufficient for our Monte Carlo based simulation method.

Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con2

70

The surfaces can be obtained from the authors on request.

Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs1

71

In order to simulate the population of the fine structure states the adiabatic potential surfaces have to be modified, so that they include the spin–orbit interaction.

Type: Motivation |
Advantage: None |
Novelty: None |
ConceptID: Mot3

72

This can be achieved approximately by diagonalizing the matrix for every point of the 3p-manifold.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met19

73

In the asymptotic region the three adiabatic energies

*e*_{1},*e*_{2}and*e*_{3}are degenerate and the correct splitting *a*with the spin–orbit coupling constant*a*= 11.5 cm^{−1}results for Na(3p).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod12

74

Usually with positive detuning the optical excitation of the collision complex leads exclusively to the population of the highest of the 3p states.

Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met20

75

Therefore under condition of strictly adiabatic dynamics all sodium atoms would be in the

^{2}P_{3/2}state after the collision.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod13

76

The involved surfaces however show conical intersections which lead to non-adiabatic transitions and to the population of the lower

^{2}P_{1/2}state.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod14

77

In the asymptotic limit of infinite relative velocity of the colliding particles the

^{2}P_{1/2}/^{2}P_{3/2}ratio is given by projection of the molecular state after excitation onto the final atomlike states^{2}P_{J}.^{5}
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met21

78

Thus for high collision energy the relative

^{2}P_{1/2}population is expected to approach 33%.
Type: Hypothesis |
Advantage: None |
Novelty: None |
ConceptID: Hyp1

### The semiclassical trajectory surface-hopping method

79

We use the semiclassical trajectory surface-hopping method developed by Preston and Tully

^{12–15}which will be reviewed shortly in the following.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met22

80

For every trajectory the probability of being on the excited surface number

*j*is given by the squared absolute value of a time dependent coefficient*c*_{j}(*t*).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod15

81

The equation for the time derivative

*ċ*_{j}(*t*) can be obtained by inserting the expansion for the electronic wave function*Ψ*^{el}(*r*,*R*(*t*)) into the time-dependent Schrödinger equation.^{14,16}
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod15

82

The nuclear coordinates symbolized by

*R*(*t*) are assumed to be functions of time*t*obeying classical equations of motion;*r*stands for the coordinates of the electrons.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod15

83

In general the resulting expression for

*ċ*_{j}(*t*) contains matrix elements 〈*ψ*_{j}|*H*^{el}|*ψ*_{k}〉 of the electronic Hamiltonian as well as matrix elements of the nuclear coordinate derivatives 〈*ψ*_{j}|∂/∂*R*|*ψ*_{k}〉.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod15

84

Conventional quantum chemical programs yield adiabatic wave functions

*ψ*^{a}that diagonalize the electronic Hamiltonian within the space of the electronic wave functions.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met23

85

The calculation of matrix elements of ∂/∂

*R*with these wave functions (neglected in the Born–Oppenheimer approximation) is very time consuming.
Type: Method |
Advantage: No |
Novelty: Old |
ConceptID: Met23

86

On the other hand these matrix elements become large in the region of a conical intersection or avoided crossing and are responsible for the radiationless transitions between different (adiabatic) electronic states.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac7

87

A unitary transformation leads to a (quasi)-diabatic basis set

*ψ*^{d}of electronic functions,^{17–20}which change as smoothly as possible with the nuclear coordinates and where the matrix elements of ∂/∂*R*are often assumed to be negligibly small.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met24

88

In the case of only two involved states

*ψ*d1 and*ψ*d2 this transformation is determined by the definition of an angle*ϕ*describing the mixing of the adiabatic states and depending on the nuclear coordinates.^{21}
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod16

89

In the diabatic picture the matrix of the electronic Hamiltonian is no longer diagonal, but the matrix elements

*H*_{12}= 〈*ψ*d1|*H*^{el}|*ψ*d2〉 can easily be calculated^{14,21}as soon as the mixing angle is known.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod16

90

The expression for the off-diagonal element reads

*H*_{12}= (*E*_{1}−*E*_{2})sin*ϕ*cos*ϕ*, where*E*_{1}and*E*_{2}are the adiabatic energies; the corresponding diagonal expressions are*H*_{11}=*E*_{1}cos^{2}*ϕ*+*E*_{2}sin^{2}*ϕ*and*H*_{22}=*E*_{1}sin^{2}*ϕ*+*E*_{2}cos^{2}*ϕ*.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod16

91

With

*a*_{j}and*b*_{j}we denote the complex time-dependent coefficients of eqn. (2) in the diabatic and adiabatic basis respectively.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod17

92

Using the mixing angle

*ϕ*the coefficients of the two different pictures can immediately be transformed into each other.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod17

93

Including the partial derivatives with respect to the nuclear coordinates the following system of differential equations is obtained

^{16,21}where the matrix elements*H*_{jk}(*R*(*τ*)) and 〈∂/∂*R*〉_{jk}(*R*(*τ*)) are implicit functions of the runtime*τ*on the trajectory.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod17

94

These semiclassical equations and the classical equations of motion are solved simultaneously by numerical integration.

Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met25

95

After transformation to (quasi)diabatic electronic wave functions the matrix elements 〈∂/∂

*R*〉_{jk}may be neglected for intermediate relative velocities*Ṙ*.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met26

96

However in our examples of optical collisions the correct

^{2}P_{1/2}population for the asymptotic case of high velocity can only be obtained by maintaining the full operator.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met27

97

A simple approximation for 〈∂/∂

*R*〉_{jk}will be introduced in section 4.1.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met27

98

The transition to the uppermost 3p surface takes place when a trajectory crosses the line (Condon configuration) corresponding to the appropriate energy difference defined by the detuning.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod18

99

At this moment the time-dependent coefficients in the adiabatic picture are set to

*b*_{1}= 1 and*b*_{2}= 0.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod18

100

Every trajectory is computed twice, corresponding to the two possibilities of excitation on the incoming and on the outgoing part of the path.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met28

101

We denote these two cases with type 1 and type 2 respectively.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met28

102

The (complex valued) diabatic coefficients

*a*_{j}are computed by numerical integration using eqn. (3) and transformed to*b*_{1}=*a*_{1}cos*ϕ*+*a*_{2}sin*ϕ*and*b*_{2}= −*a*_{1}sin*ϕ*+*a*_{2}cos*ϕ*after every time step.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met29

103

When the probability |

*b*_{i}|^{2}falls below 0.5 a random decision is made whether to hop to the other surface or not, and then the coefficients are reset accordingly.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met29

104

In the case of a hop the momentum of nuclear motion has to be corrected in order to ensure energy conservation.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met29

105

Following this procedure several hops between the two in plane states including the spin–orbit correction (1) are possible per trajectory.

Type: Method |
Advantage: Yes |
Novelty: New |
ConceptID: Met29

106

As is to be expected most hops occur not too far from the conical intersection.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res2

107

The number of trajectories ending in the two states

^{2}P_{1/2}and^{2}P_{3/2}have to be weighted with the values of the impact parameter before being interpreted as probability.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met30

### Diabatization

108

In the neighbourhood of a conical intersection the adiabatic wave functions are rapidly changing.

Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs2

109

This is illustrated

*e.g.*by the behaviour of the orientation of the transition moment vector to the higher^{2}A′ state (shown in Figs. 1–3) which in fact maps the orientation of the Na(3p) orbital in the final state.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs2

110

Due to the orthogonality constraint the 3p orbital in the lower excited

^{2}A′ state (not shown) has an orientation rotated by 90°.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res3

111

Diabatic wave functions are obtained by a transformation of the mixed adiabatic states (in our case a 2 × 2 rotation) such that a suitable molecular property behaves as smoothly as possible.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met31

112

For that we compare the orientation of the 3p orbitals with the orientation of a classical 3p quadrupole

*Q*in the field of the electric multipoles of the target molecule.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met31

113

The classical quadrupole–dipole term of the interaction potential

*e.g.*reads where*μ*is the dipole moment of the molecule.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod19

114

The position of the sodium atom relative to the molecule is given in polar coordinates (

*R*,*ϑ*) and*α*is the direction angle of the Na quadrupole (with*α*= 0 defined as the direction of the vector connecting the molecule and the sodium atom).
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod20

115

Analogous expressions hold for the quadrupole-quadrupole interaction and for higher interaction terms; in general the potential

*V*(*R*,*ϑ*,*α*) is a sum of such terms.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod21

116

For energy calculations this approach would not be sufficient because the polarization effects are missing.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res4

117

But as we need only orientation angles and since polarizabilities are essentially isotropic and do not affect the 3p orientation our approximation seems to be justified.

Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con3

118

In the classical picture we cannot distinguish two interacting electronic states.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod22

119

Hence there is no conical intersection, but there are two mutually orthogonal angles

*α*_{1}and*α*_{2}corresponding to minimal and maximal interaction energy.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod22

120

The angles

*α*_{i}can be determined by the condition The direction*α*_{1}as a function of*R*and*ϑ*is also shown in Figs. 1 to 3.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met32

121

In the asymptotic region far away from the conical intersection the angles

*α*_{i}coincide with the two directions of the transition dipole moments to the ground state, but they behave smoothly also near the intersection points.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res5

122

Therefore we force the 3p orbitals of the diabatic wave functions to coincide with these classical angles everywhere.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met33

123

In other words the deviation of the transition dipole vector direction to the highest 3p surface from the angle

*α*_{1}can serve as a “diabatization angle”*ϕ*defining the transformation to the (quasi)-diabatic states.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con4

124

Note that the classical orientation angles are only weakly dependent on the exact values of the molecular multipoles.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac8

125

In all interaction terms (

*e.g.*eqn. (4)) the Na quadrupole moment is only a factor which does not affect the*α*values of stationary energy.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac8

126

Furthermore it is rather the ratio of the molecular multipole moments than their absolute value which determines the orientation angles.

Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac9

127

The function

*α*(*R*,*ϑ*) has been approximated using the values of the molecular dipole and quadrupole moments given in section 3.1.
Type: Method |
Advantage: None |
Novelty: Old |
ConceptID: Met34

128

In addition in order to improve the agreement of the two angles in the areas where there is no mixing of the adiabatic states we have adjusted in a semiempirical way the hexadecupole terms for all three investigated molecules.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met35

129

The Figs. 1 to 3 show the directions of

*α*_{1}and of the transition dipole vectors for the three molecules.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs3

130

The mixing angle

*ϕ*has been determined as difference between these two directions.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met36

131

Thus in principle the transformation can directly be read from the figures.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res6

132

In our examples with only two involved states a necessary condition for a considerable change of probability in (3) is that the absolute values of

*H*_{12}and thus of sin*ϕ*and cos*ϕ*are not too small, and this is the case whenever the two adiabatic states are considerably mixed.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res7

133

It can easily be seen in Figs. 1 to 3 that the necessary mixing angles around 45° are found only in rather restricted regions close to the conical intersections.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res8

## Results and discussion

### Population of the Na(3p) fine structure states

134

Fig. 4 shows the resulting fine structure populations for optical collisions of Na with the molecules N

_{2}and N_{2}O for the rotational temperature 100 K and at a detuning value of 120 cm^{−1}.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs4

135

To facilitate comparison of our trajectory simulations with the experimental results

^{3}(and Fig. 5), measured at a laboratory scattering angle of 18°, the abscissa is chosen to be the relative kinetic energy after collision.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met37

136

Our

^{2}P_{1/2}fractions represented by solid lines were computed neglecting the matrix elements of ∂/∂*R*in the diabatic picture.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met38

137

Instead in order to ensure the correct asymptotic probability (33%) for high relative velocity a simple correction was applied: the matrix element 〈∂/∂

*R*〉(*R*,*ϑ*) in (3) was assumed to be proportional to the radial derivative of the mixing angle*ϕ*.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met38

138

The proportionality constant

*f*in*Ṙ*〈∂/∂*R*〉_{12}≈*fṘ*∂*ϕ*(*R*,*ϑ*)/∂*R*was adjusted individually for N_{2}and N_{2}O at high velocities.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met38

139

Maintaining the same constant in the full range of relative motion leads to the dashed curves in Fig. 4.

Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs5

140

The error bars in the figures were determined by separate statistical test runs and do not represent any systematic error of our calculation method.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res9

141

A qualitative discussion of the results has to consider the position of the conical intersection as well as the Na-molecule distance at which the optical excitation occurs.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met39

142

For Na + N

_{2}the location of the conical intersection is definitely the most favorable for a high^{2}P_{1/2}fraction.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res10

143

Our calculated percentage, which is slightly decreasing in the investigated range of collision energy, is in excellent agreement with measurement.

^{3}
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con5

144

The relative contribution of type 2 trajectories on the

^{2}*P*_{1/2}surface decreases from about 50% to about 25% with increasing relative velocity.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs6

145

For the Na + N

_{2}O system the 120 cm^{−1}Condon line is situated at a radial distance which lies just inside the conical intersection.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs7

146

Thus less than 10% of the

^{2}P_{1/2}trajectories are of type 2 except for very low collisional energy.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res11

147

The experimental population of the lower spin–orbit state is reproduced satisfactorily by our computations for this collision complex too.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res12

148

Only excitations to the uppermost 3p surface were taken into account in the calculations for Na + N

_{2}and Na + N_{2}O. In the case of Na + C_{2}H_{2}corresponding simulations lead to a population of the^{2}P_{1/2}state which is considerably much lower than according to the experiment.^{3}
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res13

149

The conical intersection is located in collinear geometry for C

_{2}H_{2}, what means that much fewer trajectories pass through this region of strong mixing than for N_{2}and N_{2}O, where the non-adiabatic coupling is strong at T-shaped geometry.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res14

150

Further the Condon configuration for the detuning by 120 cm

^{−1}is situated clearly outside of the conical intersection region on the C_{2}H_{2}surface.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res15

151

For that reason no non-adiabatic transitions to the

^{2}P_{1/2}state can result from type 2 trajectories for C_{2}H_{2}.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con6

152

The calculated

^{2}P_{1/2}populations for Na + C_{2}H_{2}at a detuning value of 120 cm^{−1}and for the rotational temperature 100 K are shown in Fig. 6.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs8

153

The low probability of non-adiabatic transitions in the case of Na + C

_{2}H_{2}represented by the dashed curves is not able to explain the rather high measured^{2}P_{1/2}fraction.^{3}
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res16

154

But in contrast to the other two investigated molecules the energy differences between the 3p surfaces are much smaller, so that excitations of the lower two 3p states are possible.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res17

155

At a detuning of 120 cm

^{−1}the number of excitations to the uppermost surface is about three times higher than to each one of the lower surfaces.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res18

156

After excitation to the lower one of the two mixed A′ surfaces most trajectories will contribute to the

^{2}P_{1/2}population because of the low non-adiabatic transition likelihood.
Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod23

157

The third state (A″) with the

^{2}P_{3/2}asymptotic limit was assumed not to mix with the other two states.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met40

158

The results represented by solid lines in Fig. 6 were obtained including excitations to all three surfaces and are in good agreement with experiment.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res19

159

A test calculation for Na + N

_{2}O showed that less than 15% of the excitations take place to the lower A′ state and many of these trajectories jump back to the uppermost surface.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res20

160

So it seems reasonable to consider only excitations to the highest 3p surface except for the system Na + C

_{2}H_{2}.
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con7

### Polarization and alignment

161

As mentioned in section 2 the intensity of the Na(3s → 3p) transition depends on the polarization direction angle

*β*of the excitation laser.
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac10

162

Table 1 shows classical trajectory results for the scattering angle

*Θ*_{CM}= 40° and several values of the detuning frequency.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs9

163

A qualitative discussion of this method can be found in .

^{ref. 2}
Type: Background |
Advantage: None |
Novelty: None |
ConceptID: Bac11

164

*β*is measured relative to the incidence direction of the Na atom.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met41

165

For example

*β*=*Θ*_{CM}denotes the situation where polarization is parallel to the Na momentum vector after collision.
Type: Observation |
Advantage: None |
Novelty: None |
ConceptID: Obs10

166

For all three examined molecules the weight of the type 1 trajectories is dominant, leading to an optimal polarization direction in backward direction (

*β*> 90°).
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res21

167

Figs. 1, 2 and 3 are useful for a qualitative understanding: For large detuning frequencies the transition dipole vector points more or less in radial direction towards the center of mass of the molecule for all possible Condon points.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res22

168

The distance between the condon points for type 1 and type 2 trajectories decreases with increasing detuning, resulting in alignment angles near 90° and large alignment ratios.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res23

169

For 120 cm

^{−1}detuning frequency the direction of the transition dipole vector shows more variation for Na + N_{2}O and Na + C_{2}H_{2}, explaining the rather small alignment ratios.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res24

170

The calculated alignment results may be compared with measurements for Na + N

_{2}and Na + C_{2}H_{2}^{2,3}(and Table 1).
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met42

171

While the experimental data are reproduced satisfactorily for the N

_{2}system the rotation of the alignment angle and the decrease of the alignment ratio at small detuning frequency are clearly underestimated for C_{2}H_{2}by our calculations.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res25

172

In the case of C

_{2}H_{2}excitations to all three surfaces of the 3*p*-manifold were included.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met42

## Summary and conclusion

173

The presented results for the population of the Na(3p) fine structure states after an optical collision with a linear molecule are in good agreement with the experiment.

^{3}
Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con8

174

Our simulations are based on quantum chemically calculated potential energy surfaces and dipole moment functions.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met43

175

Nuclear dynamics was treated in the framework of the semiclassical surface hopping method.

Type: Model |
Advantage: None |
Novelty: None |
ConceptID: Mod17

176

In order to test the effect of molecular vibration we introduced a harmonic potential for N

_{2}allowing the exchange of vibrational energy with the other degrees of freedom.
Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met44

177

Neither for the population of the Na(3p) fine structure states nor for the alignment direction and ratio a significant change within the statistical uncertainty could be found.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res26

178

We therefore conclude that our assumption of fixed molecular geometry is a reasonable approximation.

Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con9

179

The analysis of the results for the

^{2}*P*_{1/2}fraction after collision shows a sensitive dependence on the details of the energy surfaces, the non-adiabatic mixing of states and especially on the location of the conical intersections.
Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res27

180

The experimentally observed dependence of the Na(3s → 3p) transition intensity on the polarization direction angle of the excitation laser may be understood qualitatively from the properties of the potential and dipole moment surfaces.

Type: Method |
Advantage: None |
Novelty: New |
ConceptID: Met45

181

The deviations between experiment and simulation are largest for small detuning frequencies.

Type: Result |
Advantage: None |
Novelty: None |
ConceptID: Res28

182

The success of our simulations of the optical collision experiments is a strong indication for the correct understanding of the underlying physical processes.

Type: Conclusion |
Advantage: None |
Novelty: None |
ConceptID: Con10