A combined experimental and theoretical study on realizing and using laser controlled torsion of molecules
C. B. Madsen∗ and L. B. Madsen Lundbeck Foundation Theoretical Center for Quantum System Research, Department of Physics and Astronomy, Aarhus University, 8000 Aarhus C, Denmark
S. S. Viftrup, M. P. Johansson, T. B. Poulsen, L. Holmegaard, V. Kumarappan, and K. A. Jørgensen Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
H. Stapelfeldt† Department of Chemistry and Interdisciplinary Nanoscience Center (iNANO), Aarhus University, 8000 Aarhus C, Denmark
It is demonstrated that strong laser pulses can introduce torsional motion in the axially chi- ral molecule 3,5-diflouro-3’,5’-dibromo-biphenyl (DFDBrBPh). A nanosecond laser pulse spatially aligns the stereogenic carbon-carbon (C-C) bond axis allowing a perpendicularly polarized, intense femtosecond pulse to initiate torsional motion accompanied by a rotation about the fixed axis. We monitor the induced motion by femtosecond time-resolved Coulomb explosion imaging. Our theo- retical analysis corroborates the experimental findings and on the basis of these results we discuss future applications of laser induced torsion, viz., time-resolved studies of de-racemization and laser controlled molecular junctions based on molecules with torsion.
I. INTRODUCTION tends beyond the external degrees of freedom and has also been demonstrated for the internal degrees of free- dom such as vibrational motion [8] in molecular hydro- Gaining control over the external and the internal de- gen. Notably, the electrical field from laser pulses can grees of freedom of molecules is of great interest for a modify energy potential barriers such that photoinduced range of areas in molecular science. One approach, dat- bond breakage of a small linear molecule is guided to ing as far back as the 1920s [1, 2], relies on the applica- yield a desired final product [9, 10]. tion of inhomogenous (static or low frequency) magnetic and electric fields. This has proven most useful for con- trolling the full three-dimensional motion of molecules, including control of the velocity distribution, bringing molecules to a stand still, storing them for extended pe- riods of times and trapping them in a confined volume (for a recent review see [3]). An alternative and comple- In the case of larger molecules control of the lowest fre- mentary approach relies on the use of non-resonant non- quency vibrational modes attracts special interest since ionizing laser fields typically supplied by pulsed lasers. some of these modes correspond to motion along well- Many studies during the last decade have shown that defined reaction coordinates separating two conforma- strong non-resonant laser fields can effectively manipu- tional minima (conformers). Although many molecules late the external degrees of freedom of isolated gas phase contain large number of conformers it is often just two molecules. The manipulation results from laser-induced conformers that dominate important chemical properties, forces and torques due to the interaction between the in- for instance chirality. A particularly important exam- duced dipole moment and the laser field itself. Examples ple is found in axially chiral molecules [11, 12] such as of manipulation include deflection [4], focusing [5] and biaryl systems. In these molecules rotation about a sin- slowing [6] of molecules through the dependence of the gle stereogenic carbon-carbon (C-C) bond axis changes non-resonant polarizability interaction on the intensity the molecule from one enantiomer into the opposite enan- distribution in a laser focus. Likewise, the dependence of tiomer (mirror image). Recently, we demonstrated that arXiv:0903.4463v1 [physics.atom-ph] 25 Mar 2009 the induced dipole interaction on molecular orientation the laser-induced nonresonant polarizability interaction has proven highly useful for controlling the rotation of a can also be used to influence the internal rotation of an variety of molecules [7]. In particular, the spatial orien- axially chiral molecule around the stereogenic C-C bond tation of molecules can be sharply confined with respect axis [13]. In particular, we showed that by fixing the C-C to axes that are fixed in the laboratory. bond axis of a substituted biphenyl molecule in space, us- Molecular manipulation by induced dipole forces ex- ing laser induced alignment by a long laser pulse, it was possible to initiate torsional motion of the two phenyl rings by a short laser pulse polarized perpendicular to the fixed axis. The purpose of the present work is to ex- ∗Corresponding author: [email protected] tend our recent paper [13]. In particular, we discuss our †Corresponding author: [email protected] theoretical modelling in detail. 2
II. OUTLINE OF THE STRATEGY FOR LASER CONTROLLED TORSION
A ns laser pulse aligns [7] the stereogenic axis of 3,5- −3 diflouro-3’,5’-dibromo-biphenyl (DFDBrBPh) along its x 10 polarization. This enables a much shorter fs pulse, which 2 we will refer to as the kick pulse, to both initiate tor- sional motion and set the molecule into rotation around 1 the fixed axis. The internal as well as external rotational 0 motion is monitored by fs time-resolved Coulomb explo- sion imaging. −1 The inset of Fig. 1 shows a model of the DFDBrBPh molecule with the stereogenic axis marked by red. Ac- −2 cording to our quantum chemical calculations, the laser- S R a a free torsional potential, illustrated by the red dotted −3 curve in Fig. 1, has minima at dihedral angles of hφdi = −4 Kick ◦ ±39 corresponding to the Ra and Sa enantiomeric forms, pulse and where the dihedral angle is the angle between the two −5 −2 phenyl rings. The twisted equilibrium shape is charac- 0 Time−dependent torsional potential (arb. units) 2 teristic of biphenyl compounds. The traditional view is 4 t (ps) 25 50 75 100 that this non-planarity results from competition between p −100 −75 −50 −25 0 stabilization via conjugation of the π-orbital systems and Dihedral angle (degrees) steric repulsion between the ortho-positioned atoms [14]. We might mention that alternative explanations have re- FIG. 1: (Color online) To illustrate the principle of laser- cently been discussed [15, 16, 17]. Our strategy for con- induced torsion and time-resolved de-racemization we show trolling φ relies on a transient modification of the field- the calculated kick pulse induced time-dependent torsional d potential as a function of time measured with respect to the free potential curve by a fs kick pulse (Fig. 1). The mod- center of the pulse and dihedral angle φd between the Br- ification, caused by the nonresonant polarizability inter- and F-phenyl planes. The asymmetry in the potential is ob- action [7, 9, 18] converts the initial stationary quantum tained by orienting the molecules (here with the Br-phenyl states, localized near the minima of the torsional poten- plane out of the paper), 3D aligning them, and by polarizing tial, into vibrational wave packets, i.e., coherent super- the kick pulse at an angle of 13◦ with respect to the second positions of several quantum states. The temporal evo- most polarizable axis (SMPA) (see Fig. 2). The red dotted lution of these wave packets gives rise to time-dependent curve at large times illustrates the laser-free time-independent torsional motion, and the grey and the black curves in torsional potential. For the Sa enantiomer, starting out with ◦ hφ ii = −39 , the time varying potential induces an oscilla- Fig. 1 show the calculated expectation values of φd for a d particular set of laser parameters. The stereogenic C-C tory motion (grey curve) corresponding to torsion confined within the initial well. By contrast, due to the induced asym- bond axis is fixed in the laboratory frame by adiabatic metry between the two wells, the initial R enantiomer is alignment utilizing a moderately intense, linearly polar- a traversing the central torsional barrier, and ends up as an Sa ized ns laser pulse [7, 19]. The ns pulse is intense enough enantiomer undergoing internal rotation (black curve). The to keep the C-C axis tightly confined, yet weak enough to kick pulse has an intensity of 1.2×1013 W/cm2 and a duration modify the torsional potential only slightly. The fs kick of 1.0 ps (FWHM). The torsional motion may be monitored pulse is applied with its polarization perpendicular to the by fs time-resolved Coulomb explosion imaging. The inset aligned C-C bond axis to ensure primarily influence on shows a model of the DFDBrBPh molecule with the stero- torsion while avoiding other possible vibrational motion. genic axis marked by red (grey). The DFDBrBPh compound is prepared as a racemate and the experiment is carried out on isolated molecules at rotational temperatures of a few Kelvin. Under these two phenyl planes at the time of the probe pulse. conditions we have equal numbers of molecules initially localized in the −39◦ or 39◦ conformation (see Fig. 1) and no thermally induced transitions between the two III. THEORETICAL METHODS conformations occur. The intense probe pulse, sent at A. Modelling the dynamics of the molecule time tp with respect to the kick pulse, removes several electrons from the molecules, thereby triggering Coulomb explosion into ionic fragments. In particular, the Br+ To account for the motion of the nuclei of the DFD- and F+ fragment ions recoil in the planes defined by the BrBPh molecule we make use of two sets of coordinate Br- and F-phenyl rings. By recording the velocities of systems (see Fig. 2): A molecular fixed (MF) frame at- both ion species with two-dimensional ion imaging [19], tached to the molecule and a laboratory fixed (LF) frame we determine the instantaneous orientation of each of the specified by the lasers. The MF coordinates are chosen 3
X obtain x 1 ∂2 1 ∂2 Br SMPA Hmol = − 2 + − 2 + Vtor(φd) 2I ∂Φ 2Irel ∂φd
= HΦ + Hφd . (2) F Here I = IBr + IF is the total moment of inertia of the molecule for rotation around the stereogenic axis and Y Irel = IBrIF/I is a relative moment of inertia for the two phenyl rings around the axis. A full rotation of ei- ther phenyl ring leaves us with the same molecule im- y plying 2π-periodic boundary conditions of the eigenfunc- tions of Hmol from Eq. (1), i.e., ψ(φBr +2πm, φF +2πn) = ψ(φBr, φF), with m and n integers. We shall assume that this property also applies to Φ and φd, so that we simply ˜ need to consider eigenfunctions ψ(Φ, φd) = ξ(Φ)χ(φd) of FIG. 2: Model of DFDBrBPh along with the molecular fixed Eq. (2) that separates into rotation of the molecule as (MF) xyz and the laboratory fixed XYZ coordinate axes. described by the 2π-periodic function ξ(Φ) and torsion The dotted line indicates the second most polarizable axis accounted for by the 2π-periodic function χ(φd). [Rig- ◦ (SMPA), which is located 11 from the Br-phenyl ring and orously, the bounds for Φ depend on φd [20]. At the 28◦ from the F-phenyl ring. time scales of interest, however, a molecule generally only makes a small fraction of a full rotation so that the effect of Φ-φd coupling is negligible.] The separation is physi- cally motivated by considering the energy scales related with the z axis along the stereogenic axis pointing from to rotation and torsion. The energy scale of the prior is the phenyl ring with the bromines towards the phenyl given by ~2/(2I) = 1.3 µeV. For torsion, on the other ring with the flourines, and the x axis is chosen along the hand, the relevant energy is determined by the torsional phenyl ring with the bromines. The LF coordinates are potential, and a harmonic approximation of the poten- chosen with the Z axis along the polarization direction tial near the minimum at 39◦ yields a frequency corre- of the ns pulse and the X axis along the kick pulse po- sponding to the energy 3.1 meV. Using that the period larization direction. In agreement with the experimental of motion is of the order Planck’s constant divided by observations (cf. Sec. V), we assume that the stereogenic the energy, we therefore see that the molecule rotates (Φ axis of the DFDBrBPh molecule is aligned along the Z change) on a nanosecond time scale, whereas the torsion axis. Also, we neglect all normal modes, except the low- (φd change) is of picosecond duration. est one, which corresponds to torsion (see Sec. III B for The small energy separation of the rotational levels details on quantum chemical calculations, justifying, e.g., compared to that of the torsional levels implies that many this assumption). In this simplified situation and in the rotational states will be occupied in thermal equillibrium. absence of the kick pulse the task is reduced to describ- Consequently, we expect the Φ dynamics to behave classi- ing the coupled rotations of the two (rigid) phenyl rings cally (see also Fig. 7), and we only treat the torsional dy- of the molecule as given in the LF frame by the Hamil- namics quantum mechanically. The torsion is described tonian [Atomic units, me = e = a0 = ~ = 1, are used by the stationary Schr¨odingerequation throughout unless indicated otherwise]
Hφd |χν i = Eν |χν i, (3)
1 ∂2 1 ∂2 with ν = 1, 2,... denoting the energy eigenstates. To H = − − + V (φ − φ ), (1) mol 2 2 tor Br F solve this equation we expand the Hamiltonian onto an 2IBr ∂φBr 2IF ∂φF orthonormal basis of 2π-periodic functions, which we di- agonalize. The eigenstates below the torsional barriers where φi, i = Br, F, is the angle of the i phenyl ring with exhibit an almost exact fourfold degeneracy, and linear respect to the kick pulse polarization axis, Ii the mo- combinations of the approximately degenerate states will ment of inertia for rotation of the i phenyl ring around also be stable at the time scales of the experiment. In the stereogenic axis (IBr = 8911925, IF = 1864705) and particular, we can produce such stable states that are Vtor(φBr − φF) is the laser-free torsional potential as ob- localized in the wells of the torsional potential. We will (i) tained from quantum chemical calculations (Fig. 3). By denote these states by |Lνmin i, with i = 1, 2, 3, 4 and νmin changing the coordinates to the dihedral angle φd = the smallest ν of the four degenerate states. In Fig. 3 φBr − φF between the two phenyl rings and the weighted we show the first four energy eigenstates and the corre- azimuthal angle Φ = (1 − η)φBr + ηφF, characterizing sponding localized states along with the torsional poten- the rotation of the molecule, with η = IF/(IF + IBr), we tial. The experiment is generally carried out on a gas of 4
a reduction of the dihedral angle, whereas an increase of the dihedral angle is resulting from the ⊥-geometry. Hence, the overall effect of the kick pulse will be to align the molecules into the k-geometry and drive a change of the dihedral angle. We proceed with the quantitative analysis of the field- induced dynamics. The state |χν i is no longer an eigen- state of the torsion when the kick pulse is applied, but will develop into a wave packet
Φ X Φ −iEν0 (t−t0) |χν i → |χν (t)i = cν0 (t)e |χν0 i (6) ν0 with the coefficients determined by
Φ X Φ −i(Eν −Eν0 )(t−t0) c˙ν0 (t) = −i cν (t)e hχν0 |Vkick(Φ, t)|χν i. FIG. 3: The field-free torsional states. (a) The four first (al- ν most) degenerate energy eigenstates lying 1.71 meV above (7) the minimum of the torsional potential. From linear combi- Once we find these new states of torsion the expecta- nations of these we obtain the corresponding four localized tion value of any operator O can be evaluated by tracing states shown in (b). The torsional potential is indicated with the product of the density matrix with the operator, i.e., a dashed line. hO(t)i = Tr[ρ(t)O]. In particular, we need the expecta- tion value of the dihedral angle and of the kick poten- tial from Eq. (5). From the latter we obtain the torque, molecules in a mixed state given by the density matrix which causes the molecule to rotate into the k-geometry. 4 If the molecule lies at an angle Φ it is exposed to a torque X X X ρ(t ) = P (i) |L(i) ihL(i) | + P |χ ihχ |, −∂hV (Φ, t)i/∂Φ directed along the Z axis and hence 0 νmin νmin νmin ν ν ν kick νmin i=1 ν achieves an angular acceleration given by (4) ∂hV (Φ, t)i IΦ¨ = − kick . (8) with the P ’s being weight factors that sum to unity, ∂Φ and where the latter sum includes the states above the torsional barrier. The kick pulse introduces dy- Along with the initial conditions given at Eq. (4), the namics that we determine from the alignment poten- Eqs. (7) and (8) provide a set of coupled differential equa- tions that may be integrated to obtain the coordinates tial created by the pulse. By definition F (t) = F0(t)Xˆ for the kick pulse, and we use a Gaussian envelope Φ(t) and hφd(t)i at time t. Rather than solving these 2 coupled equations, we, however, assume that the angle Φ F0(t) = F0 exp{− ln(2)[t/(τFWHM/2)] /2} corresponding to an intensity full width at half maximum (FWHM) of has the constant value Φ0 during the short time interval of the kick pulse, and we integrate Eq. (8) twice to arrive τFWHM. We calculate the polarizability tensor in the MF frame as a function of the dihedral angle (see Table I) and at consequently transform into the LF frame by the appli- Z ∞ 1 ∂ 0 0 cation of directional cosine matrices [21]. We then arrive Φ(t) = Φ0 − t dt hVkick(Φ, t )i . (9) I ∂Φ −∞ at the interaction Φ=Φ0 1 Consistently, we solve Eq. (7) with Φ fixed at Φ . V (Φ, φ , t) = − F 2(t)[α (φ ) cos2(Φ + ηφ ) 0 kick d 4 0 xx d d For the present work we start out with a distribution 2 + αyy(φd) sin (Φ + ηφd) of values for Φ0 in the interval [−π, π] along with an ini- tial density matrix, ρ(t ) and we propagate each member − 2α (φ ) cos(Φ + ηφ ) sin(Φ + ηφ )]. 0 xy d d d of the distribution according to Eqs. (7) and (9). This (5) procedure yields an ensemble of Φ(t)’s and hφd(t)i’s and For the succeeding analysis it is helpful to note a few we finally employ the relations φBr(t) = Φ(t) + ηhφd(t)i qualitative features of the potential and we refer to Fig. 2 and φF(t) = Φ(t) − (1 − η)hφd(t)i to find the angular for graphical guidance. For a fixed dihedral angle the distributions of bromines and flourines. potential is minimal when the second most polarizable axis (SMPA), which is located 11◦ from the Br-phenyl ring and 28◦ from the F-phenyl ring, is parallel to the X B. Quantum chemical calculations axis. We will denote this by the k-geometry. Conversely, the potential is maximal if the molecule is rotated 90◦ To determine the field-free torsional potential we ob- from the k-geometry, and we will denote this by the ⊥- tain the geometries of DFDBrBPh in various conforma- geometry. Next, in the k-geometry the potential favors tions at the density functional theory level [22, 23], with 5
TABLE I: The table lists the relevant polarizability components, αij, of DFDBrBPh in the MF frame as a function of the dihedral angle, φd. The components are π-periodic and fulfill αxx(φd) = αxx(π−φd), αyy(φd) = αyy(π−φd) and αxy(φd) = −αxy(π−φd). Also, αyx = αxy.
φd (degrees) 0 15 30 45 60 75 90 αxx (atomic units) 217.694 215.590 209.463 200.975 192.431 186.240 184.048 αyy (atomic units) 92.352 95.201 102.634 112.658 122.693 130.015 132.639 αxy (atomic units) 0.000 -9.488 -16.360 -18.810 -16.225 -9.295 0.000
B3LYP, Becke’s three-parameter hybrid exchange func- tional [24] in connection with the Lee-Yang-Parr corre- TABLE II: Vibrational frequencies and relative Raman cross sections of the eight lowest lying normal modes of DFD- lation [25] functional. We model the correlation of the BrBPh, computed at B3LYP/TZVPP level. uniform electron gas with the Vosko-Wilk-Nusair VWN5 formulation [26], and we use doubly polarized triple-zeta mode frequency (cm−1) cross section (%) quality basis-sets, TZVPP [27]. This level of theory has ν1 23.9 100.0 been shown to be suitable for the study of biphenyls [28]. ν2 55.7 21.4 To scan the potential energy surface of the relative tor- ν3 60.7 2.8 sional angle of the phenyl rings we optimize all coordi- ν4 111.9 5.1 ν 145.2 1.3 nates except the torsional angle. 5 ν 194.7 3.6 We do all electronic structure calculations with the 6 ν7 195.3 1.3 Turbomole 5.10 program suite [29]. In Table I we give ν8 213.4 0.8 the dynamic polarizabilities computed at a laser wave- length of 1064 nm. We note that the polarisabilities are quite stable with respect to choice of laser frequency. Using a value of 800 nm changes the individual compo- Na2CO3 (3 mL, degassed). We next heat the mixture ◦ nents only minutely (maximally by 2%). The differences to 90 C under argon for 5.5 h and cool it to room tem- between static and dynamic polarisabilities are equally perature. The mixture is diluted with water and ex- small. We calculate the vibrational frequencies analyt- tracted three times with Et2O. We dry the combined ically [30], and a normal mode analysis, using internal organic fractions with MgSO4, filter and concentrate to coordinates, shows that the lowest frequency vibration give approximately 1.4 g of a crude mixture. Based on corresponds almost purely to the torsion of the phenyl NMR we find that this mixture contains both the de- rings. We compute vibrational Raman cross-sections sired mono-Suzuki coupling product [31], double-coupled from derivatives of the polarizability tensor with respect products, and a homo-coupling product from the boronic to the normal modes of vibration. Table II shows the acid. Careful column chromatography on SiO2 eluting frequencies and Raman cross sections of the lowest vi- with n-hexane and c-hexane finally allows us to isolate brations of DFDBrBPh, computed at B3LYP/TZVPP 182 mg 3,5-dibromo-3’,5’-difluorobiphenyl (DFDBrBPh) level. The results provide the justification for assuming as a white amorphous solid (contains <3% of the homo- that the kick pulse primarily interacts with the lowest coupling product DFDBrBPh). NMR spectroscopic data 1 normal mode, since its Raman scattering cross section (see also Fig. 4): H NMR δ (400 MHz, CDCl3) 7.69 (t, J is over four times higher than the second lowest normal 1.7 Hz, 1H), 7.61 (d, J 1.7 Hz, 2H), 7.05 (m, 2H), 6.85 (m, 13 mode and 20 times higher than any other normal mode. 1H). C NMR δ (100 MHz, CDCl3) 163.5 (2C), 142.3, We note that the character of the second lowest (and sec- 141.5, 133.8, 128.9 (2C), 123.5 (2C), 110.1 (2C), 103.8. ond strongest Raman active) mode contains no torsional Elemental analysis, calculated: C 41.42%, H 1.74%, Br motion. 45.92%; found: C 41.20%, 1.66%, Br 45.90%.
IV. EXPERIMENTAL METHODS B. Setup of the laser experiment
A. Synthesis of 3,5-dibromo-3’,5’-difluorobiphenyl A schematic of the experimental setup is shown in Fig. 5. Less than hundred milligram of solid DFDBrBPh To a round-bottom flask equipped with a magnetic is heated to 130◦C and expanded in 90 bar of He into stirring bar under argon we add 3,5-difluorophenyl- vacuum using a pulsed Even-Lavie valve [19, 32]. The boronic acid (474 mg, 3.0 mmol), 1,3,5-tribromobenzene resulting supersonic molecular beam is skimmed and sub- (1134 mg, 3.6 mmol), and Pd(PPh3)4 (69.6 mg, 2 mol%). sequently crossed by three, focused, pulsed laser beams at The flask is then evacuated and backfilled with argon 90◦. The first laser beam originates from a Q-switched twice. To this mixture of solids we add toluene (18 Nd:YAG laser and consists of 9-ns-long pulses at 1064 mL, degassed), EtOH (6 mL, degassed), and 2M aq. nm. It is focused to a spotsize ω0 ∼ 35 µm yielding a 6
Computer-controlled CCD MCP + delay line phosphor 1064 nm, 9 ns, 1.5 J, 20 Hz (a) (b) Br Nd:YAG laser Br Br Pd(PPh3)4 F B(OH)2 Mu-metal + F Br F Br Br + Br Br o o F 0YAG 0 beam F F mirror splitter F+ 2.22 2.09 1.08 1.00 Grating 9.0 8.0 7.0 6.0 5.0 4.0 3.0 2.0 1.0 0.0 30 cm lens stretcher (c) (d) Br f1 (ppm) Ar-filled waveguide F Br Br + prism compressor F F F Br 800 nm, 130 fs, 2 mJ, 1 kHz F Even-Lavie valve F Ti:Sapphire laser F F FIG. 5: (Color online) Schematic representation of the ex- 1.00 0.05 perimental setup, showing the vacuum system and the three 180 160 140 120 100 80 60 40 20 0 -108.0 -109.0 -110.0 -111.0 -112.0 f1 (ppm) f1 (ppm) laser beam paths. The Nd:YAG (alignment) pulses and the probe pulses are polarized vertically, i.e., perpendicular to the FIG. 4: (Color online) (a), Synthesis of DFDBrBPh by Pd- detector plane, whereas the fs kick pulse is polarized in the catalyzed cross-coupling. (b)-(d), NMR-spectra of synthetic plane of the detector. DFDBrBPh. Integrated signals are indicated by green (grey) curves, and the corresponding numbers are situated directly (a) -0.87 ps 0.47 ps 1.47 ps 2.47 ps 5.14 ps 12.15 ps below the respective curves.
F + peak intensity of ∼ 7×1011 W/cm2. The purpose of these pulses is to align the stereogenic axis of the DFDBrBPh molecules. The second laser beam consists of 0.7-ps-long Br+ kick pulses (λ = 800 nm), obtained by passing part of the output from an amplified Ti:Sapphire fs laser system (b) (c) 27 0.50.5 27 through a grating stretcher. The spotsize of ω0 = 43 2626 12 2 > µm yields a peak intensity of ∼ 5 × 10 W/cm . These 0.40.4 > 2525 + kick pulses are used to induce torsion and rotation of the 0.30.3 F-kick.pol 2424 -kick. pol. φ + F < 2323 signal (arb. units) molecules. The third laser beam consists of 25-fs-long 0.20.2 < + pulses, which we achieve by passing another part of the F 2222 0.10.1 + 2121 output from the Ti:Sapphire system through a hollow F signal (arb. units) -180 -9090 180 0 27090 180 360 1.01.0 1.5 2.0 2.0 2.53.0 3.0 3.5 4.0 4.0 + (degrees) -kick.φF pol. + t (ps) wave guide compressor setup. The pulses have a spotsize F-kick.pol (degrees) ptp (ps) 14 2 ω0 = 25 µm and a peak intensity of ∼ 2 × 10 W/cm . These probe pulses Coulomb explode the molecules and FIG. 6: (Color online) (a) Ion images of F+ and Br+ frag- the produced ions are extracted by a weak static electric ments at probe times tp. The ns pulse is polarized per- field, in a velocity imaging geometry, and projected onto pendicularly to the image (detector) plane and the 5 × 1012 W/cm2, 0.7 ps (FWHM) kick pulse is polarized horizontally. a two-dimensional detector, consisting of a micro channel + (b) Angular distribution of the F ions, at tp =1.47 ps, ob- plate (MCP) detector backed by a phosphor screen. We + detect the F+ and Br+ ions separately by time gating tained by radially integrating the corresponding F ion im- age. The splitting of the pairwise peaks is twice the av- of the MCP detector around the respective arrival times + erage angle, hφF+- kick. poli, between the F ion recoil (and of the ions. We record the ion images on the phosphor thus the F-phenyl plane) and the kick pulse polarization. (c) screen by a CCD camera and determine the coordinates hφF+- kick. poli as a function of tp, for times where a clear four- of each individual ion hit. Figure 6(a) displays the ion peak structure is visible in the angular distributions. The images of F+ and Br+ at selected probe times. curve is a fit of the sum of a linear and a harmonic function to the experimental points (squares).
V. RESULTS AND DISCUSSION absence of ions in the innermost region is only compati- First, we establish that the ns pulse is able to hold ble with the C-C bond axis being aligned perpendicular the C-C bond axis of DFDBrBPh along its polarization to the detector plane and the F-phenyl planes uniformly + direction. At tp = −0.87 ps the F image, displayed in distributed around the C-C bond axis. The correspond- Fig. 6(a), is almost circularly symmetric and the small ing Br+ image is also essentially circularly symmetric al- deviation from circular symmetry is explained by noting though the hollow-out of the center is not as pronounced that the kick pulse has a finite value at −0.87 ps. The as for F+. 7
Next, we investigate the effect of the kick pulse. Al- 0.8 ready at tp = 0.47 ps the circular symmetry is broken and both F+ and Br+ ions start to localize around the 0.75 + polarization direction of the kick pulse. The F ions re- 0.7 main within the radial range of the ring structure, which shows that the kick pulse does not cause any significant 0.65 distortion of the C-C bond axis alignment but rather ini- 0.6 > Φ tiates an overall rotation of the molecule around this axis. 2 0.55 The latter is expected since the torque imparted by the kick pulse will force the SMPA (Fig. 2) to align along the (a) (a) No of 100 counts Prob. F 80 dens. 50 6 (b) 60 40 Br 0 20 5 R a -50 (c) 0.08 ps 1.47 ps 2.47 ps 42 4 40 -100 38 > (degrees) 36 (b) d (degrees) φ 3 < 100 34 d -1 0 1 2 34 φ S a tp (ps) 50 FIG. 8: (Color online) Angular distributions of (a) F-phenyl 2 and (b) Br-phenyl rings at tp = 0.08, 1.47 and 2.47 ps. (c) Expectation value of the dihedral angle for a molecule starting 0 out with the SMPA aligned along the kick pulse polarization. The kick pulse is as in Fig. 6. 1 -50 + fully consistent with the Br ion distribution, recorded -100 0 at 2.47 ps. -2 0 2 4 The theoretical value hφdi exhibits oscillations tp (ps) (Fig. 8(c)) with a period of ∼ 1.2 ps and amplitude of ∼ 2.45◦. The period agrees well with the experimental FIG. 9: (Color online) Time evolution of the dihedral angle value (∼ 1 ps), and smaller modulation in hφdi is ex- pected in the experiment (∼ 0.6◦) partly since here the for a molecule starting out as (a) an Ra or (b) an Sa enan- tiomer. Initially, the molecule is 3D oriented with the Br- SMPA is not pre-aligned and partly due to deviations phenyl end pointing out of the paper and the SMPA aligned from the stated laser peak intensity. The behavior is as- at an angle of 13◦ with respect to the kick pulse polariza- cribed to a wave packet of normal modes in the torsional tion. The kick pulse triggering the torsional motion has a double well potential (Fig. 3(b)) for a molecule starting peak intensity of 1.2 × 1013 W/cm2 and duration (FWHM) of out with the SMPA aligned. The qualitative agreement 1.0 ps. The torsional potential is scaled down by 1/4 rather of Figs. 6(b) and 8(c) corroborates the interpretation of than by increasing the kick strength. In practise reduction the kick pulse inducing time-dependent torsional motion. of the torsional barrier is possible, for instance, by modifying the aromatic rings or by using halogen substituted bipheny- lacetylene. VI. PERSPECTIVES We will now discuss some perspectives of the demon- and the kick pulse needs to be increased either through strated laser control of torsion. For one thing, we sug- higher intensity, a longer kick pulse or by trains of syn- gest a time-resolved study of de-racemization [34, 35, 36, chronized kick pulses [44, 45, 46]. 37, 38], where one enantiomer is steered into its mirror Assuming initial orientation and confinement of the form. To do this we break the inversion symmetry of the SMPA, we have calculated the φd dynamics for both con- C-C bond axis by orienting each molecule. This allows formations the DFDBrBPh molecule with a reduced tor- us to discriminate between the two enantiomeric forms. sional barrier (see caption to Fig. 9). The results are Theory [39, 40] and experiment [41, 42] show that orien- shown in Fig. 9, and clearly, the present method would tation can be added to 3D alignment by combining the allow for a time-resolved study of the transition from one ns alignment pulse with a static electric field. Next, the enantiomer into the other. A quantitative analysis of the excitation of torsional motion is optimized by eliminating efficiency of the process shows that after the pulse 99% initial rotation around the C-C bond through alignment of the molecules starting out as Ra has changed into Sa of the SMPA (Fig. 2) prior to the kick pulse using an el- enantiomers, whereas only 13% of the Sa enantiomers liptically rather than a linearly polarized ns pulse [7, 43]. changed into Ra. The inverse process causing an excess Finally, the interaction strength between the molecule of Ra enantiomers, is achieved simply by inverting the 9 orientation of the molecules. a fs kick pulse, polarized perpendicular to the fixed axis, We believe that the results presented in this paper also was applied. The torsion as well as the overall rotation of fuel the field of molecular junctions [47, 48, 49] with new the molecule was monitored with fs time resolution using possibilities. The conductivity of molecules like biphenyl a delayed probe pulse. To substantiate the observations, will depend on the dihedral angle [50], and as such this we developed an original theoretical model, which pro- type of molecules, if placed between two leads, can func- duces results in good agreement with the experiment. tion as a molecular junction [51, 52, 53, 54] that may The model divides the dynamics of the molecule into a be used to control the charge flow and act as, e.g., a the rotation with classical behavior and quantum me- switch [43]. Importantly, we can realize quantum control chanical torsional dynamics. As such the model provides of such a system by means of kick pulse schemes and in a transparent physical intepretation of essential features this way tailor very specific torsional wave packets [55] of laser controlled torsion of DFDBrBPh. Based on the to dictate the detailed time-dependent dihedral motion theoretical model, we further discussed the extension of and thereby the current flow through the molecular junc- the experiment to a setup that will facilitate a fs time- tion. Seen through the light of controllability, molecular resolved study of de-racemization. Finally, we pointed to junctions based on laser controlled torsion complements perspectives of using molecules like DFDBrBPh in molec- the available schemes, such as the mechanical break junc- ular junctions in order to realize laser controlled charge tions [56] and the previously suggested resonantly light flow. driven molecular junctions [57, 58, 59]. VIII. ACKNOWLEDGEMENTS VII. CONCLUSION This work was supported by the Carlsberg Founda- We presented a fs time-resolved study on laser con- tion, The Lundbeck Foundation, Danish Natural Re- trolled torsion of axially chiral molecules. Experimen- search Foundation, the Danish Natural Science Research tally, the symmetry axis of DFDBrBPh molecules was Council, and the Danish Research Agency (Grant no. held fixed in space by a long alignment laser pulse while 217-05-0081). [1] H. Kallmann and F. Reiche, Z. Phys. 6, 352 (1921). R. F. W. Bader, Chem. Eur. J. 9, 1940 (2003). [2] O. Stern, Z. Phys. A 39, 751 (1926). [16] J. Poater, M. Sol`a,and F. M. Bickelhaupt, Chem. Eur. [3] S. Y. T. van de Meerakker, H. L. Bethlem, and G. Meijer, J. 12, 2889 (2006). Nature Phys. 4, 595 (2008). [17] L. F. Pacios, Struct. Chem. 18, 785 (2007). [4] H. Stapelfeldt, Hirofumi Sakai, E. Constant, and [18] B. J. Sussman, J. G. Underwood, R. Lausten, M. Y. P. B. Corkum, Phys. Rev. Lett. 79, 2787 (1997). Ivanov, and A. Stolow, Phys. Rev. A 73, 053403 (2006). [5] Hoi Sung Chung, Bum Suk Zhao, Sung Hyup Lee, Sungu [19] V. Kumarappan, C. Z. Bisgaard, S. S. Viftrup, Hwang, Keunchang Cho Cho, Sang-Hee Shim, S. Wee L. Holmegaard, and H. Stapelfeldt, J. Chem. Phys. 125, Kyung Kang, and Doo Soo Chung, J. Chem. Phys. 114, 194309 (2006). 8293 (2001). [20] K. Hoki, D. Kr¨oner,and J. Manz, Chem. Phys. 267, 59 [6] R. Fulton, A. I. Bishop, M. N. Shneider, and P. F. Barker, (2001). Nature Phys. 2, 465 (2006), ISSN 1745-2473. [21] R. N. Zare, Angular momentum: understanding spatial [7] H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, aspects in chemistry and physics (John Wiley, New York, 543 (2003). 1988). [8] H. Niikura, P. B. Corkum, and D. Villeneuve, Phys. Rev. [22] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 Lett. 90, 203601 (2003). (1964). [9] B. J. Sussman, D. Townsend, M. Y. Ivanov, and [23] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). A. Stolow, Science 314, 278 (2006). [24] A. D. Becke, J. Chem. Phys. 98, 5648 (1993). [10] B. J. Sussman, M. Y. Ivanov, and A. Stolow, Phys. Rev. [25] C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 A 71, 051401(R) (2005). (1988). [11] E. L. Eliel and A. H. Wilen, Stereochemistry of Organic [26] S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, Compounds, Chapter 14 (Wiley, New York, 1994). 1200 (1980). [12] G. Bringmann, A. Mortimer, P. Keller, M. Gresser, [27] A. Sch¨afer,C. Huber, and R. Ahlrichs, J. Chem. Phys. J. Garner, and M. Breuning, Angew. Chem. Int. Ed. 44, 100, 5829 (1994). 5384 (2005). [28] M. P. Johansson and J. Olsen, J. Chem. Theory Comput. [13] C. B. Madsen, L. B. Madsen, S. S. Viftrup, M. P. Jo- 4, 1460 (2008). hansson, T. B. Poulsen, L. Holmegaard, V. Kumarap- [29] R. Ahlrichs, M. B¨ar,M. H¨aser,H. Horn, and C. K¨olmel, pan, K. A. Jørgensen, and H. Stapelfeldt, Phys. Rev. Chem. Phys. Lett. 162, 165 (1989). Lett. 102, 073007 (2009). [30] P. Deglmann, F. Furche, and R. Ahlrichs, Chem. Phys. [14] F. Grein, J. Phys. Chem. A 106, 3823 (2002). Lett. 362, 511 (2002). [15] C. F. Matta, J. Hernand´ez-Trujillo, T.-H. Tang, and [31] N. Miyaura and A. Suzuki, Chem. Rev. 95, 2457 (1995). 10 [32] U. Even, J. Jortner, D. Noy, N. Lavie, and C. Cossart- [46] K. F. Lee, I. V. Litvinyuk, P. W. Dooley, M. Spanner, Magos, J. Chem. Phys. 112, 8068 (2000). D. M. Villeneuve, and P. B. Corkum, J. Phys. B 37, L43 [33] S. Viftrup, V. Kumarappan, S. Trippel, H. Stapelfeldt, (2004). E. Hamilton, and T. Seideman, Phys. Rev. Lett. 99, [47] A. Nitzan and M. A. Ratner, Science 300, 1384 (2003). 143602 (2007). [48] Z. J. Donhauser et al. , Science 292, 2303 (2001). [34] K. Faber, Chem. Eur. J. 7, 5004 (2001). [49] J. Chen and M. Reed, Chem. Phys. 281, 127 (2002). [35] Y. Fujimura, L. Gonz´alez,K. Hoki, J. Manz, and Y. Oht- [50] J. M. Seminario, A. G. Zacarias, and J. M. Tour, suki, Chem. Phys. Lett. 306, 1 (1999). J. Am. Chem. Soc. 120, 3970 (1998). [36] M. Shapiro, E. Frishman, and P. Brumer, Phys. Rev. [51] M. C´ıˇzek,M.ˇ Thoss, and W. Domcke, Phys. Rev. B 70, Lett. 84, 1669 (2000). 125406 (2004). [37] D. Kr¨onerand B. Klaum¨unzer,Chem. Phys. 338, 268 [52] M. C´ıˇzek,M.ˇ Thoss, and W. Domcke, Czechoslovak Jour- (2007). nal of Physcis 55, 189 (2005). [38] D. Kr¨onerand B. Klaum¨unzer,Phys. Chem. Chem. Phys. [53] C. Benesch, M. C´ıˇzek,M.ˇ Thoss, and W. Domcke, Chem. 9, 5009 (2007). Phys. Lett. 430, 355 (2006). [39] B. Friedrich and D. Herschbach, J. Phys. Chem. A 103, [54] C. Benesch, M. C´ıˇzek,J.ˇ Klimeˇs,I. Kondov, M. Thoss, 10280 (1999). and W. Domcke, J. Phys. Chem. C 112, 9880 (2008). [40] B. Friedrich and D. Herschbach, J. Chem. Phys. 111, [55] J. Werschnik and E. K. U. Gross, J. Phys. B 40, R175 6157 (1999). (2007). [41] H. Tanji, S. Minemoto, and H. Sakai, Phys. Rev. A 72, [56] S. Wu, M. T. Gon´zalez,R. Huber, S. Grunder, M. Mayor, 063401 (2005). C. Sch¨onenberger, and M. Calame, Nature Nanotech. 3, [42] L. Holmegaard, J. H. Nielsen, I. Nevo, H. Stapelfeldt, 569 (2008). F. Filsinger, J. K¨upper, and G. Meijer, Phys. Rev. Lett. [57] S. Loudwig and H. Bayley, J. Am. Chem. Soc. 128, 12404 102, 023001 (2009). (2006). [43] S. Ramakrishna and T. Seideman, Phys. Rev. Lett. 99, [58] C. Zhang, M.-H. Du, H.-P-Cheng, X.-G. Zhang, A. E. 103001 (2007). Roitberg, and J. L. Krause, Phys. Rev. Lett. 92, 158301 [44] M. Leibscher, I. Sh. Averbukh, and H. Rabitz, Phys. Rev. (2004). Lett. 90, 213001 (2003). [59] M. del Valle, R. Guti´errez,C. Tejedor, and G. Cuniberti, [45] C. Z. Bisgaard, S. S. Viftrup, and H. Stapelfeldt, Phys. Nature Nanotech. 2, 176 (2007). Rev. A 73, 053410 (2006).