Diffractive imaging of a rotational wavepacket in nitrogen molecules with femtosecond megaelectronvolt electron pulses

Abstract

Imaging changes in molecular geometries on their natural femtosecond timescale with sub-Angström spatial precision is one of the critical challenges in the chemical sciences, as the nuclear geometry changes determine the molecular reactivity. For photoexcited molecules, the nuclear dynamics determine the photoenergy conversion path and efficiency. Here we report a gas-phase electron diffraction experiment using megaelectronvolt (MeV) electrons, where we captured the rotational wavepacket dynamics of nonadiabatically laser-aligned nitrogen molecules. We achieved a combination of 100 fs root-mean-squared temporal resolution and sub-Angstrom (0.76 Å) spatial resolution that makes it possible to resolve the position of the nuclei within the molecule. In addition, the diffraction patterns reveal the angular distribution of the molecules, which changes from prolate (aligned) to oblate (anti-aligned) in 300 fs. Our results demonstrate a significant and promising step towards making atomically resolved movies of molecular reactions.

ARTICLE

Received 3 Aug 2015 | Accepted 3 Mar 2016 | Published 5 Apr 2016

Jie Yang1, Markus Guehr2,3, Theodore Vecchione4, Matthew S. Robinson1, Renkai Li4, Nick Hartmann4, Xiaozhe Shen4, Ryan Coffee4, Jeff Corbett4, Alan Fry4, Kelly Gaffney4, Tais Gorkhover4, Carsten Hast4, Keith Jobe4, Igor Makasyuk4, Alexander Reid4, Joseph Robinson4, Sharon Vetter4, Fenglin Wang4, Stephen Weathersby4, Charles Yoneda4, Martin Centurion1 & Xijie Wang4

Imaging changes in molecular geometries on their natural femtosecond timescale with sub-Angstrm spatial precision is one of the critical challenges in the chemical sciences, as the nuclear geometry changes determine the molecular reactivity. For photoexcited molecules, the nuclear dynamics determine the photoenergy conversion path and efciency. Here we report a gas-phase electron diffraction experiment using megaelectronvolt (MeV) electrons, where we captured the rotational wavepacket dynamics of nonadiabatically laser-aligned nitrogen molecules. We achieved a combination of 100 fs root-mean-squared temporal resolution and sub-Angstrom (0.76 ) spatial resolution that makes it possible to resolve the position of the nuclei within the molecule. In addition, the diffraction patterns reveal the angular distribution of the molecules, which changes from prolate (aligned) to oblate (anti-aligned) in 300 fs. Our results demonstrate a signicant and promising step towards making atomically resolved movies of molecular reactions.

1 Department of Physics and Astronomy, University of Nebraska-Lincoln, 855 N 16th Street, Lincoln, Nebraska 68588, USA. 2 PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA. 3 Institute of Physics and Astronomy, Potsdam University, Potsdam 14476, Germany. 4 SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA. Correspondence and requests for materials should be addressed to M.G.(email: mailto:[email protected]

Web End [email protected] ) or to M.C. (email: mailto:[email protected]

Web End [email protected] ) or to X.W. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms11232 OPEN

Diffractive imaging of a rotational wavepacket in nitrogen molecules with femtosecond megaelectronvolt electron pulses

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232

The dynamical behaviour of molecules is governed by the complex interplay of a correlated system of nuclei and electrons which interact via Coulomb and exchange forces.

Determining how the individual nuclei within a molecule move relative to one another during a molecular transformation represents a key step to understanding chemical reactivity. Developments in the eld of femtochemistry enabled capturing the motion of nuclear wavepackets using purely spectroscopic measurements with femtosecond laser pulses18. For small molecules, structural information on the nuclear geometry can be indirectly inferred based on exact knowledge of spectral transitions involved in the probe process. For larger molecules, this inference becomes unfeasible. Diffraction techniques, which provide direct access to the position of each atom within a molecule, are a more effective approach. Over the last few decades, great strides have been taken to observe ultrafast dynamics with ultrashort X-ray pulses from synchrotrons and free electron lasers (FELs)912 and via diffraction of short electron pulses1316. With the advent of FELs, X-ray diffraction experiments have now reached a sub-100 fs temporal resolution. The spatial resolution is sufcient to observe larger-scale molecular processes like light-induced opening of a six-membered ring, but not to resolve individual atoms because of the limited wavelength available in FELs12. Laser-induced electron diffraction17,18 and photoelectron holography19 can also provide spatial and temporal resolution simultaneously, but require signicant theoretical input to retrieve the interatomic distances. Ultrafast electron diffraction (UED) from gas-phase molecules has achieved sub-Angstrom spatial resolution but the temporal resolution has not been sufcient to observe molecular geometry changes on the femtosecond timescale20,21. Here we report a crucial advance in gas-phase UED using multi-MeV relativistic electron pulses to achieve a temporal resolution of 100 fs root-mean-squared (RMS), or 230 fs full-width at half-maximum (FWHM) that opens the door to observing the motion of individual nuclei that result from structural changes during photochemical reactions of isolated molecules.

The majority of UED experiments have been performed on solid-state samples22,23, which are ideal to understand collective effects in condensed media such as superconductivity, heat transport and magnetism. Gaseous molecules are ideal to study prototypical processes in chemistry. They also provide a direct link between experiments and quantum chemical calculations, which can be performed on the highest level for isolated molecules2325. Early gas phase UED studies employed stroboscopic electron diffraction26. A breakthrough was achieved when picosecond electron beams became available in late 1980s (ref. 27). Zewail and co-workers were able to resolve non-equilibrium molecular structures28, transient molecular structures15 and radiationless dark structures16 using picosecond UED. However, in order to resolve molecular geometry changes in real time, a better temporal resolution is required. In the context of the crucial excited state photoisomerization reactions24, a 200-fs resolution is suitable for exploring the nuclear dynamics of the isolated azobenzene isomerization29. In addition, this temporal resolution is sufcient to explore photoprotection processes of nucleobases3033.

To improve UED temporal resolution, two major effects must be overcome: space-charge repulsion between electrons34 and velocity mismatch35 resulting from the electron pulse lagging behind light pulses used for the molecular excitation. Both of these limitations can be minimized using relativistic MeV electrons3645. The longitudinal space-charge pulse elongation is proportional to 1/b2g5, where b v/c, g (1 b2) 1/2 is

the Lorentz factor, and v and c are the speed of electrons and the speed of light in vacuum, respectively46. Concerning the

velocity mismatch, electrons with 3.7 MeV kinetic energy travels at v 0.993c. This results in only 5 fs delay with respect to an

optical pulse for a typical 200 mm interaction length.

Here we present the ultrafast laser-induced rotational dynamics of N2 molecules in the gas phase. We impulsively excited a rotational wavepacket in a N2 gas sample using a 35-fs FWHM laser pulse, spectrally centred at 800 nm, which interacts non-resonantly with the anisotropic molecular polarizability tensor. In the impulsive alignment regime, the alignment laser pulse duration is much shorter than the rotational period of the molecule, and the ensemble reaches its maximum degree of alignment after the interaction with the laser pulse. The phase evolution of the rotational wavepacket results in rotational revivals, with molecules alternating between aligned and anti-aligned directions47,48. Previous efforts to investigate rotational wavepackets with UED have captured dynamics on picosecond timescales49,50. The rotational dynamics of laser-aligned N2 has been previously observed with optical birefringence51, strong eld ionization52, high harmonic generation53,54 and Auger electron spectroscopy55. We have observed the temporal evolution of the full wavepacket revival with an 8.35 ps period by quantifying the anisotropy in the diffraction patterns. The temporal resolution was determined by a tting routine using the measured dynamics. Furthermore, we have retrieved molecular images of the aligned and anti-aligned molecular ensemble with atomic resolution.

ResultsExperimental layout and static diffraction. The experimental layout is shown in Fig. 1. The electron pulse (blue) is diffracted from the nitrogen gas jet (grey), which is introduced into the vacuum chamber using a pulsed nozzle (black). The diffraction pattern is captured by a phosphor screen and a detector. The 800-nm alignment laser pulse (red) is directed to the target and removed from the vacuum chamber by two holey mirrors at a 45 angle to the electron beam. The full setup is discussed in detail elsewhere56. The electron and the laser beam have a small angle of B5. This makes the velocity of the laser 0.996c along the electron beam direction, very close to 0.993c, the velocity of

fs Laser

N2 nozzle

Detector

3.7 MeV e

Figure 1 | Experimental layout. A sketch of the experimental setup.

A 3.7-MeV pulsed electron beam (blue) is directed towards a nitrogen gas jet (grey). The gas jet is introduced into the vacuum chamber using a pulsed nozzle (black). The pump laser pulse (red) is deected by two ring-shaped mirrors. The laser propagates at a small angle (B5) with respect to the electron beam as it traverses the target, and is then deected away from the detector by the second mirror. The electron diffraction pattern is recorded with a phosphor screen located 3.1 m downstream from the interaction region. The unscattered electron beam is transmitted through a hole in the phosphor screen.

2 NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232 ARTICLE

the electron beam. Therefore, the velocity mismatch is nearly eliminated. More details of the experiment can be found in the Experimental Setup section in the Methods.

The diffraction pattern is commonly expressed as a function of the momentum transfer

s

4pl siny=2 1 where l is the wavelength of the electron beam and y is the angle between the scattered and transmitted electrons. For a 3.7-MeV electron beam, l 0.30 pm. The total scattering intensity Itot is

the sum of the atomic scattering intensity Iat and the molecular scattering intensity Imol. Iat is dened as

Iat X

N

i1

fi s

j j2 2

where N is the number of atoms in the molecule and fi is the elastic scattering amplitude for the ith atom. For MeV electrons, fi can be calculated using the ELSEPA program57.

The structural information of the molecule is encoded in the molecular scattering intensity, Imol, given by

Imol X

N

i1

sinsrijsrij 3

where Zi is the scattering phase of the ith atom and rij is the distance between the ith and the jth atoms20.

The so-called modied diffraction intensity is dened as

sM s

s

Imols

Iats

4

The spatial resolution, d, is determined by the maximum measured s value smax using the formula

d 2p=smax 5 In Fig. 2a, the red curve shows an azimuthally averaged raw pattern after subtraction with a dark background pattern that is taken with the electron beam turned off. The azimuthal average of the subtracted background pattern is shown in the green curve in Fig. 2a, and the black dashed curve shows the experimental scattering background, determined using a standard tting procedure20. This tted background includes the atomic scattering Iat, scattering from background gas and other types of background scatterings. The details of the bond length

XN

j 6 i

fis

j j fjs

cosZ

i

determination are explained in the Bond Length Determination from Static Pattern section in the Methods. Figure 2b shows the azimuthally averaged experimental and theoretical modied diffraction intensity for static diffraction pattern. The bond length was determined to be 1.0730.027 , in agreement with the previously measured N2 bond length of 1.098 (ref. 58). The2.5% measurement uncertainty is due to the uncertainty in the calibration of the sample-to-detector distance and electron energy, as explained in the Experimental Setup section in the Methods. It should be noted that for time-resolved experiments usually the ground-state structure is known, in which case the static diffraction patterns can be used as a calibration. We show later that the interatomic distance can be determined very accurately from diffraction of transiently aligned molecules if the static diffraction is used as a calibration. In this experiment, scattering signal is available in the region between 3.5 and 12 1. The tting procedure relies on the zeros of sM(s) and this reduces the available data in Fig. 2 to s44.5 1. The measurement agrees well with the simulation up to sB12 1.

Temporal evolution of N2 alignment. For impulsive alignment, the full rotational revival is expected at t 1/2cB, where B is the

rotational constant and c is the speed of light in vacuum. For N2 molecules, the full revival is at 8.35 ps (B 1.998 cm 1). The

rapid evolution of the angular distribution can be used to determine the temporal resolution of the measurement technique.

Diffraction patterns from a molecular ensemble aligned with a polarization in the detector plane are not circularly symmetric, contrary to the static case in equation (3). The anisotropy, a(t), in a diffraction pattern can be used to trace the temporal evolution of alignment49. In addition, a(t) is a self-normalized parameter that is extracted directly from the diffraction patterns (see Data Processing of Diffraction Patterns and Anisotropy in the Methods). Figure 3 shows the temporal evolution of the simulated and experimentally measured a(t) values. The experimental data are recorded with 100 fs steps in the delay between laser and electron pulses, and at each point data are collected for 2 min (14,400 shots at a repetition rate of 120 Hz). The simulation is composed of two parts: an impulsive alignment simulation that calculates the angular distribution at different delay times59, followed by a diffraction pattern simulation based on the modelled angular distribution60. The anisotropy of the simulated patterns is calculated using the same method as for the experimental patterns. The details of the

Zj

a b

1

0.5

3,000

2,500

2,000

1,500

1,000

500

0 4 6 8 10 12 0

1

Counts

sM(s)

0

0.5

Simulation Data

5 10 15

s (1) s (1)

Figure 2 | Static N2 diffraction. (a) The red curve shows an azimuthally averaged raw experimental diffraction pattern after subtraction with a dark background pattern that is taken without the electron beam. The green curve shows the azimuthal average of the subtracted dark background pattern. The dashed black curve shows the tted background, including the atomic scattering Iat and other background scattering. The vertical axis is average detector counts per pixel per minute of exposure time, averaged over B100 min. The green curve varies between 0.8 and 2.5 counts. (b) The theoretical (blue) and experimental (red) modied diffraction intensity sM from N2 gas, which shows the enhanced diffraction rings. The experimental sM is calculated from the diffraction pattern in part (a).

NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232

[afii9851](a.u.)

1.041.02 1

0.98

0.960 50

RMS temporal resolution (fs)

100 150

Data Simulation

Anisotropy (a.u.)

0.61

0.6

0.59

0.58

0.57

2 0 2 4

Delay (ps)

6 8 10

0.45

0.41

0.37

0.32

0.27

<cos2 >

Figure 3 | Temporal evolution of the N2 rotational wavepacket. Anisotropy in the diffraction patterns from experimental data (red) and simulation (blue) versus time. Statistical error bars for the rst few points (before alignment) are shown to illustrate the uncertainty of this measurement. The right-hand side axis gives the degree of alignment ocos2a4 for the simulated curve. The simulation parameters of initial rotational temperature, alignment laser uence, temporal resolution and rescaling factor are obtained from a tting routine. Each data point is accumulated over 2 min. Reduced w2 error versus

RMS temporal resolution in the four-parameter t is shown in the inset.

simulation are explained in the Alignment Simulation and Diffraction Pattern Simulation section in the Methods. For diatomic molecules, the rotational states can be described by |J, M4, where J and M are the quantum numbers of the total and the z-component of the angular momentum, respectively. The angular distribution of |J, M4 is given by Laplaces spherical harmonics YmJ. Before the excitation, the ensemble is described by a

Boltzmann distribution of rotational states taking into account the exact nuclear statistics. Initially, the angular distribution is isotropic. For the initial rotational temperature of our sample of54 K, all values of J up to 10 are considerably populated before excitation. The laser pulse interacts with the polarizability tensor of the molecules, which is highest along the internuclear axis. Classically, it induces a dipole that minimizes its energy by aligning the molecular axis parallel to the laser eld. Quantum mechanically, the laser induces transitions of DJ 0, 2 and

DM 0. After excitation, the maximum value of J is about 14. The

phase of each rotational state evolves with time. When the rotational states are in phase, the angular distribution becomes aligned or anti-aligned with respect to the laser polarization axis.

The temporal evolution of the alignment has a rich structure that varies on a fast timescale. For example, at the half-revival (delay of 4 ps), the distribution changes from aligned to anti-aligned in 300 fs. The alignment peak corresponds to a prolate angular distribution, with the long axis along the direction of the laser polarization. The angular distribution during anti-alignment is oblate, with the molecules preferentially lying in a plane perpendicular to the laser polarization. The full revival at around 8 ps shows a similarly fast transition from oblate to prolate distribution. In between the revivals, the anisotropy of the diffraction pattern captures additional dynamics.

We have used the fast-changing distribution to characterize the temporal resolution of the measurement. The shape of the aligned molecular ensemble is determined by the initial temperature of the molecules and the uence of the alignment laser pulse61. In this experiment, the laser uence was measured to be 2.0 J cm 2 and the initial temperature was estimated to be 65 K, using a supersonic expansion model62. At these parameters, the wavepacket revivals are relatively sharp compared with the temporal resolution. The limited temporal resolution effectively blurs out the rotational dynamics and has a signicant effect on

the observed structure. We performed rst a two-parameter w2 tting by xing the laser uence and initial rotational temperature to the measured and calculated values, and then a full four-parameter t where all parameters were allowed to vary.

The two-parameter t, varying only the temporal resolution and a re-scaling factor that accounts for the spatial overlap between laser and electron pulse, returned a temporal resolution of 85 fs RMS (200 fs FWHM). The four-parameter t, which varies temporal resolution in addition to initial rotational temperature, laser uence and the re-scaling factor, accounts for uncertainties in the laser uence, initial temperature and the fraction of molecules excited by the laser. In the simulation, we assume that all excited molecules are exposed to the same laser intensity. This method achieved a best t resolution of 100 fs RMS (230 fs FWHM), with an initial temperature of 54 K and a laser uence of 1.8 J cm 2, which are comparable to our initial estimates for the two-parameter t. The re-scaling factor was 0.42, meaning that the best match to the experimental data is when the simulation assumes that 42% of the molecules in the diffraction volume are excited by the laser. The simulation shown in Fig. 3 shows the best t for all four parameters. The left axis in the gure shows the anisotropy in the diffraction pattern, and the axis on the right shows the degree of alignment quantied by the ocos2a4 value from simulation, where a is the angle between the molecular axis and the laser polarization direction. The reduced w2 error in the four-parameter t versus the RMS temporal resolution is shown as an inset in Fig. 3. We obtained a 0.96 reduced w2, indicating a good t. More details of the tting are explained in the Temporal Evolution Fitting section in the Methods.

The 100-fs (RMS) overall temporal resolution of this experiment is consistent with our expectations, based on the performance study of this machine56. A simulation showed that the electron bunch length was 70 fs RMS at the interaction region. Measurements of the phase and amplitude stability of the radio-frequency (RF) gun lead to an expected time of arrival jitter of 50 fs RMS56. The calculated overall temporal resolution was then 87 fs RMS, or 205 fs FWHM, which is close to the measured value.

Molecular images with different angular distribution. High-resolution molecular diffraction images were retrieved for

4 NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232 ARTICLE

prolate and oblate ensembles at the half revival. Diffraction patterns with adequate signal-to-noise ratio were recorded with 60 and 90 min of integration time for oblate and prolate distributions, respectively. We use diffraction-difference patterns to remove the experimental background and the diffraction signal from unexcited molecules. The diffraction intensity difference is given by DI(t) I(t) I(t 5 ps), where t 0 corresponds to

the maximum of the rst alignment peak after laser. Before t 0.4 ps when the pump laser arrives, the angular distribution

is isotropic. In the two-dimensional (2D) diffraction patterns, the

smax 8.3 1, corresponding to a spatial resolution of 0.76 .

This makes it possible to observe molecular structures with resolution better than the shortest possible bond lengths.

Figure 4 shows the experimental (left panels) and simulated (right panels) 2D diffraction patterns and their corresponding Fourier transforms for prolate and oblate distribution at the half revival. Figure 4a,b shows the experimental and simulated diffraction-difference patterns DI for the prolate distribution.

The diffraction pattern was captured at a time delay of 3.8 ps after the rst alignment peak (Fig. 3a). The diffraction-difference

a b

Experimental pattern (prolate)

FT of experimental pattern (prolate)

Simulated pattern (prolate)

0.5

0.5

0.5

0.5

8

8

4

4

1 )

1 )

s y(

0

0

arb. unit

s y(

0

0

arb. unit

4

4

8

1

0

1

8

8

4

0

4

8

8

4

0

4

8

sx (1)

sx (1)

c d

2

0

1

1

FT of simulated pattern (prolate)

1

0

1

1

r y()

arb. unit

r y()

2

0

1

arb. unit

1

0

1

1

0

1

rx ()

rx ()

e f

Experimental pattern (oblate)

8

8

0.5

0.5

Simulated pattern (oblate)

8

8

0.5

0.5

4

4

1 )

1 )

s y(

0

0

arb. unit

s y(

0

0

arb. unit

4

4

8

8

8

4

0

4

8

4

0

4

sx (1)

sx (1)

g h

FT of experimental pattern (oblate)

1

0

1

1

2

0

1

FT of simulated pattern (oblate)

1

0

1

1

r y()

arb. unit

r y()

2

0

1

arb. unit

1

0

1

1

0

1

rx ()

rx ()

Figure 4 | 2D N2 diffraction patterns at half revival. (a) Experimentally measured and (b) simulated diffraction-difference patterns of the prolate distribution. Images shown in c,d are Fourier transforms of a,b, respectively. The Fourier transform of the diffraction-difference patterns show the changes in the angular distribution of the molecules. The positive regions (red colour) indicate where the population has increased and the negative regions (blue colour) indicate where the population has decreased. (e) Experimentally measured and (f) simulated diffraction-difference pattern of the oblate distribution. Images shown in g,h are Fourier transforms of e,f, respectively. In patterns (a,e), the data inside the black circles are missing due to the beam stop. They are obtained by extrapolating the pattern and letting the counts smoothly go to zero towards the centre. For illustrative purpose, angular distributions are shown on the side of panel (c,d,g,h) for visual guidance. In these angular distributions, the colour code indicates polar angle.

NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232

pattern is anisotropic as a result of molecular alignment. The experimental pattern shows excellent agreement with the simulation for a range of s 3.58.3 1.

Figure 4c,d shows the Fourier transforms of the difference signals in Fig. 4a,b, respectively. The Fourier transform of the diffraction pattern displays the autocorrelation of the molecular structure, convolved with the angular distribution and projected onto the detector plane. The centre of the Fourier transforms goes to zero because they are generated from the difference of two diffraction patterns. For a diatomic molecule, the autocorrelation is directly related to the molecular image. For more complex molecules, an image of the structure can be reconstructed using phase retrieval algorithms60. In the autocorrelation functions depicted in Fig. 4c,d, the positive regions indicate an increase in population, and the negative regions indicate a decrease. Specically, Fig. 4c,d indicates that the population of molecules that are lying perpendicular to the laser polarization has decreased, whereas the population parallel to the polarization (vertical in Fig. 4) has increased, that is, more molecules are aligned along the vertical direction. Similarly, Fig. 4e,f shows the measured and simulated diffraction pattern for the oblate distribution at half revival, corresponding to a time delay of 4.1 ps. Figure 4g,h are the Fourier transforms of Fig. 4e,f, which show that the molecular ensemble is aligned in the horizontal plane.

Bond length measurement from diffraction patterns. We have used a tting method to extract the N2 bond length from diffraction patterns of aligned and anti-aligned molecules, using the static diffraction pattern as a calibration. The details of the tting are given in the section Bond Length Fitting from Aligned Patterns in the Methods. The extracted bond length is1.0910.036 for the prolate distribution and 1.0960.056 for the oblate distribution, in good agreement with 1.098 , the bond length of the ground-state N2. The small uncertainties indicate that changes in interatomic distances could be measured very accurately with this method. The precision of determining the bond length, 0.036 and 0.056 , should not be confused with the 0.76 spatial resolution. The spatial resolution gives the capability to distinguish two bond lengths that are very close to each other, whereas the precision gives how accurate a single bond length can be determined. In real space, the resolution is determined by the width of the peak, whereas the precision is determined by how accurate one can nd the centre of the peak. Generally, for well-separated peaks, the centre can be determined to much higher accuracy than its width.

Spatial resolution of the molecular images. The spatial resolution of the 2D diffraction patterns shown in Fig. 4 can be determined in two different ways. First we can use equation (6), with smax 8.3 1 we get the spatial resolution d 0.76 .

We can also determine spatial resolution directly from the autocorrelation images (Fig. 4c,g). For example, in Fig. 4c, the spatial resolution can be determined by converting the image into polar coordinates o(r, y), then using a Gaussian function to t to along the r dimension. The FWHM of the Gaussian t was0.76 using this method, consistent with the spatial resolution obtained using equation (5) and smax 8.3 1.

Angular distribution based on 2D images of aligned molecules. The angular distribution of the molecules can be extracted from the patterns in Fig. 4c,g. The resulting distributions for prolate and oblate molecular ensembles are shown in Fig. 5a,b, respectively (see Angular Distribution section in the Methods). The prolate distribution peaks at a 0 and 180, in the direction

of the laser polarization. The oblate distribution peaks at a 90,

in the direction perpendicular to the laser polarization. The degree of alignment is commonly measured by the quantity

cos2a

R

p0 f a cos2a sina da

R

p0 f a

a

0.02

0.015

0.01

0.005

0.015

0.01

0.005

00 45 90 135 180

Population (a.u.) Population (a.u.)

0 0 45 90

[afii9825] ()

[afii9825] ()

135 180

b

Figure 5 | Angular distributions. (a) Prolate angular distribution at the revival: experimental (red), simulated (black) and simulated convolved with 100 fs RMS temporal resolution (blue). The effect of the temporal resolution is to reduce the measured degree of alignment at the peak, due to time averaging. The measurement agrees very well with the simulation after the convolution. The ocos2a4 values for red, black and blue curves are 0.41, 0.45 and 0.42, respectively. (b) Oblate angular distribution:

Experimental (red), simulated (black) and simulated convolved with 100 fs RMS temporal resolution (blue). The ocos2a4 values for red, black and blue curves are 0.28, 0.25 and 0.28, respectively. The experimental curves are measured from Fig. 4c,g, respectively.

6

where f(a) is the angular distribution and a is the angle between the molecular axis and the laser polarization. A value of the alignment parameter ocos2a4 1 corresponds to perfect

alignment, whereas random orientation gives ocos2a4 1/3.

Any value between 1/3 and 1 indicates alignment, whereas any value below 1/3 indicates anti-alignment. In Fig. 5, three curves are displayed: f(a) extracted from data (red), f(a) obtained from simulation using the uence and temperature extracted from the tting routine (black), and the simulation results convolved with the 100 fs RMS temporal resolution (blue). The temporal resolution is included by convolving a Gaussian pulse with the simulated temporal evolution of the angular distribution, and it has a signicant effect because the distribution changes from prolate to oblate in 300 fs. The extracted angular distributions are in good agreement with the simulation for both prolate and

sin a da

6 NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232 ARTICLE

oblate distributions, after taking into account the effect of the temporal resolution.

The measured angular distribution can also be used to provide additional conrmation of the temporal resolution. Using the laser parameters and the initial temperature as determined by the tting from the data in Fig. 3, we ran an optimization to nd the temporal convolution that results in a best t between the measured and simulated angular distributions. From this, we obtain a resolution of 107 fs RMS for the prolate angular distribution and 102 fs RMS for the oblate distribution, in good agreement with the value of 100 fs RMS obtained before.

DiscussionIn summary, we have shown that MeV UED can simultaneously reach a temporal resolution of 100 fs RMS and spatial resolution of 0.76 , which allows us to characterize the ultrafast evolution of a rotational wavepacket in N2, measuring the angular distribution and the acquired molecular images with atomic resolution.

This achievement opens the door to a new class of experiments where changes in molecular geometry during a chemical reaction can be followed in space and time. The current results demonstrate not only excellent temporal resolution, but also that sufcient signal-to-noise ratio for sub-Angstrom spatial resolution can be acquired with MeV electrons. Previous picosecond UED experiments with similar spatial resolution have successfully retrieved the structure of transient states15. Nitrogen has a low atomic number, Z, and therefore a low scattering cross-section, so we expect that the method will be successful for observing molecular dynamics in a large class of molecules. RF-compressed electron pulses63,64 and compact UED guns6567 have achieved similar temporal resolution in condensed matter experiments, but have not yet been successfully applied to gas-phase experiments. For future developments, a number of upgrades can potentially improve both spatial and temporal resolution signicantly. For example, a RF compression cavity could be used to compress the electron beam longitudinally68, which can potentially lead to a temporal resolution on the order of 10 fs together with orders-of-magnitude increase in charge per pulse. Time jitter and drift could be addressed by time-stamping techniques similar to those developed for X-Ray free electron lasers (XFELs)69. Spatial resolution is limited by the signal-to-noise ratio (SNR) in the scattering at larger angles, so RF compression would also improve the spatial resolution by increasing the number of electrons per pulse. Complementary metal-oxide semiconductor (CMOS) active pixel sensor could potentially achieve single electron, single shot detection and thus eliminate detector noise70, which would also improve the SNR.

Methods

Experimental setup. The electron gun used in this experiment is a replica of the photo-injector used at the LCLS facility at SLAC National Accelerator Laboratory. A 15 fC per pulse charge is generated at the photocathode and 6 fC per pulse is delivered on target. A 50-nm-thick Si3N4 membrane is used to separate the sample chamber and preserve the high vacuum in the electron gun. Roughly 60% of the charge is lost after the Si3N4 membrane and a 200-mm diameter collimator. The electron beam and pump laser are generated with a repetition rate of 120 Hz.

The N2 gas is delivered into the chamber through a pulsed valve at a 120-Hz repetition rate. The backing pressure is 0.7 bar, and the nozzle orice is 100 mm in diameter. During operation, the chamber pressure is B6 10 5 torr. The

interaction region is roughly 350 mm away from the nozzle exit. The width of the laser beam, the electron beam and the gas jet are all around 200 mm FWHM at the interaction region. The sample-to-detector distance is calibrated with diffraction from a single crystal gold sample. There is a 2.5% uncertainty in the distance calibration, due to the quality of the gold sample and quality of the electron beam after the Si3N4 membrane. The Si3N4 window has a strong diffraction and makes the emittance of the beam much worse. The electrons scattered inelastically from Si3N4 also diffract from the gold sample, making the Bragg spots larger. Spatial alignment of the gas jet, the electron beam and the pump laser is obtained by

positioning the focus of the laser approximately underneath the tip of the nozzle, and maximizing for the electron beam deection by the plasma. Time zero can then be found, to within B200 fs, by adjusting the delay between the pump and the probe beams and observing plasma deection effects on the unscattered beam20. The laser focus is adjusted to 50 mm FWHM spot size for plasma lensing, giving a peak intensity of 8 1014 W cm 2. For the alignment experiments, the lens is

driven 4 mm out of the focus to reach a spot size of 200 mm in the interaction region. The peak intensity for the experiment is 5 1013 W cm 2; only very faint

plasma is observed in the interaction region for this intensity.

The detector comprises a phosphor screen, a ring-shaped mirror at a 45 angle to the beam, an f/0.85 lens and an electron multiply charge-coupled device. A 4-mm diameter hole is drilled in the centre of the phosphor screen to allow the unscattered electron beam to pass through. Data in the region so3.5 1 are not captured due to the hole in the phosphor screen. This can be improved by replacing the phosphor screen by a similar one with a smaller hole. We have seen that by replacing the phosphor screen the missing region can be reduced to 1.6 1.

The beam divergence s is 28 mrad RMS. Using the denition l/2ps, we get a transverse coherence length of 1.7 nm.

Bond length determination from static pattern. The following steps are used to determine the bond length from the static N2 diffraction pattern20. (i)Starting from an initial guess of the bond length the modied diffraction intensity sM is calculated. (ii) A series of zero points is determined from the simulated sM.(iii) The background is determined by tting an exponentially decaying function at the zero points. (iv) The experimental sM is obtained by subtracting the background. (v) The experimental sM is compared with the simulated sM and the error is dened as the square of the difference. (vi) Steps iiv are repeated for different bond lengths, until the minimum error is found.

Data processing of diffraction patterns and anisotropy. The experimental pattern is symmetrized over four quadrants, and the central region (so3.5 1) is extrapolated from existing data by letting the pattern smoothly go to zero towards the centre.

The diffraction pattern of the prolate distribution is integrated over 90 min, and that of the oblate distribution is integrated over 60 min. The anisotropy is calculated by dividing the total counts in a horizontal cone by the total counts in a vertical cone for each diffraction pattern. The region used for this calculation is between s 3 1 and s 4.5 1. The horizontal cone has a half angle of 35

and the vertical cone has a half angle of 55.

Alignment simulation. The simulation of impulsive alignment is calculated using a linear rigid rotor interacting with a non-resonant pulse described by the time-dependent Schrodinger equation59. The temporal resolution is implemented by a convolution with a Gaussian function in time.

Diffraction pattern simulation. The simulations of diffraction patterns of a given angular distribution are calculated using an incoherent weighted sum of diffraction patterns from single molecules with different orientations. In the simulations, the atoms are assumed to be stationary at their equilibrium positions without vibrations60.

Temporal evolution tting. The temporal evolution tting shown in Fig. 3a is obtained by a w2 tting of the simulation and data. The simulated anisotropy is obtained by an alignment simulation that gives angular distribution as a function of time, followed by a diffraction simulation according to the angular distribution. The anisotropy is calculated using the same method as for data, and the temporal resolution is obtained by a convolution of the calculated anisotropy with a Gaussian beam in time. The reduced w2 error, plotted in the Fig. 3 inset, is dened asw2red 1w P

O E

s , where w is the number of degrees of freedom, O and E are the data and the best t, s is the s.e., respectively. Here s is calculated using the standard deviation of 13 patterns before t 0.

Bond length tting from aligned patterns. The optimal bond length is found by comparing the experimental diffraction pattern with many simulated patterns with different bond lengths using the least-square error. For the experimental patterns, a vertical cone with a half opening angle of 30 is selected, and this part is averaged radially to generate a one-dimensional curve. For the simulated patterns, diffraction patterns with different bond lengths and the same angular distribution are simulated, and the same region is selected to generate the one-dimensional curve. The prolate/ oblate angular distributions are obtained by the tting in Fig. 3. To get statistics, ten independent data set (6 min data in each set) were used in the ts separately, and the mean and standard deviation of the ten tting results are taken as the nal result and the standard error. In this tting, we use the static pattern (shown in Fig. 2) as a calibration of the sample-detector distance and electron beam energy.

Angular distribution. The angular distribution in Fig. 5a,b are extracted from Fig. 4c,g, respectively. The procedure includes three steps: (i) extracting the

NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232

28. Ruan, C. Y. et al. Ultrafast diffraction and structural dynamics: the nature of complex molecules far from equilibrium. Proc. Natl Acad. Sci. USA 98, 71177122 (2001).

29. Schultz, T. et al. Mechanism and dynamics of azobenzene photoisomerization.J. Am. Chem. Soc. 125, 80988099 (2003).30. Crespo-Hernandez, C., Cohen, B., Hare, P. & Kohler, B. Ultrafast excited-state dynamics in nucleic acids. Chem. Rev. 104, 19772019 (2004).

31. McFarland, B. K. et al. Ultrafast X-ray Auger probing of photoexcited molecular dynamics. Nat. Commun. 5, 4235 (2014).

32. Middleton, C. T. et al. DNA excited-state dynamics: from single bases to the double helix. Annu. Rev. Phys. Chem. 60, 217239 (2009).

33. Schreier, W. J. et al. Thymine dimerization in DNA is an ultrafast photoreaction. Science 315, 625629 (2007).

34. Siwick, B. J., Dwyer, J. R., Jordan, R. E. & Miller, R. J. D. Ultrafast electron optics: propagation dynamics of femtosecond electron packets. J. Appl. Phys. 92, 16431648 (2002).

35. Dantus, M., Kim, S. B., Williamson, J. C. & Zewail, A. H. Ultrafast electron diffraction. 5. Experimental time resolution and applications. J. Phys. Chem. 98, 27822796 (1994).

36. Hastings, J. B. et al. Ultrafast time-resolved electron diffraction with megavolt electron beams. Appl. Phys. Lett. 89, 184109 (2006).

37. Li, R. et al. Experimental demonstration of high quality MeV ultrafast electron diffraction. Rev. Sci. Instrum 80, 083303 (2009).

38. Zhu, P. et al. Femtosecond time-resolved MeV electron diffraction. N. J. Phys. 17, 063004 (2015).

39. Wang, X., Qiu, X. & Ben-Zvi, I. Experimental observation of high-brightness microbunching in a photocathode rf electron gun. Phys. Rev. E 54, R3121R3124 (1996).

40. Wang, X. J., Wu, Z. & Ihee, H. Femto-seconds electron beam diffraction using photocathode RF gun. Proc. 2003 Part Accel. Conf. 1, 420422 (2003).

41. Wang, X. J., Xiang, D., Kim, T. K. & Ihee, H. Potential of femtosecond electron diffraction using near-relativistic electrons from a photocathode RF electron gun. J. Korean Phys. Soc. 48, 390396 (2006).

42. Musumeci, P., Moody, J. T. & Scoby, C. M. Relativistic electron diffraction at the UCLA Pegasus photoinjector laboratory. Ultramicroscopy 108, 14501453 (2008).

43. MuroOka, Y. et al. Transmission-electron diffraction by MeV electron pulses. Appl. Phys. Lett. 98, 20092012 (2011).

44. Fu, F. et al. High quality single shot ultrafast MeV electron diffraction from a photocathode radio-frequency gun. Rev. Sci. Instrum 85, 083701 (2014).45. Manz, S. et al. Mapping atomic motions with ultrabright electrons: towards fundamental limits in space-time resolution. Faraday Discuss. 177, 467491 (2015).

46. Reiser, M. Theory and Design of Charged Particle Beams (Wiley, 1994).47. Stapelfeldt, H. & Seideman, T. Colloquium: aligning molecules with strong laser pulses. Rev. Mod. Phys. 75, 543557 (2003).

48. Rosca-Pruna, F. & Vrakking, M. J. Experimental observation of revival structures in picosecond laser-induced alignment of I2. Phys. Rev. Lett. 87, 153902 (2001).

49. Hensley, C. J., Yang, J. & Centurion, M. Imaging of isolated molecules with ultrafast electron pulses. Phys. Rev. Lett. 109, 133202 (2012).

50. Reckenthaeler, P. et al. Time-resolved electron diffraction from selectively aligned molecules. Phys. Rev. Lett. 102, 213001 (2009).

51. Chen, Y.-H., Varma, S. & Milchberg, H. M. Space- and time-resolved measurement of rotational wave packet revivals of linear gas molecules using single-shot supercontinuum spectral interferometry. J. Opt. Soc. Am. B 25, B122 (2008).

52. Litvinyuk, I. V. et al. Alignment-dependent strong eld ionization of molecules. Phys. Rev. Lett. 90, 233003 (2003).

53. Itatani, J. et al. Controlling high harmonic generation with molecular wave packets. Phys. Rev. Lett. 94, 123902 (2005).

54. McFarland, B. K., Farrell, J. P., Bucksbaum, P. H. & Ghr, M. High harmonic generation from multiple orbitals in N2. Science 322, 12321235 (2008).

55. Cryan, J. P. et al. Auger electron angular distribution of double core-hole states in the molecular reference frame. Phys. Rev. Lett. 105, 083004 (2010).56. Weathersby, S. P. et al. Mega-electron-volt ultrafast electron diffraction at SLAC National Accelerator Laboratory. Rev. Sci. Instrum. 86, 073702 (2015).

57. Salvat, F., Jablonski, A. & Powell, C. J. ElsepaDirac partial-wave calculation of elastic scattering of electrons and positrons by atoms, positive ions and molecules. Comput. Phys. Commun. 165, 157190 (2005).

58. Huber, K. P. & Herzberg, G. in Molecular Spectra and Molecular StructureIV. Constants of Diatomic Molecules (Springer, 1979).59. Ortigoso, J., Rodriguez, M., Gupta, M. & Friedrich, B. Time evolution of pendular states created by the interaction of molecular polarizability with a pulsed nonresonant laser eld. J. Chem. Phys. 110, 38703875 (1999).

60. Yang, J., Makhija, V., Kumarappan, V. & Centurion, M. Reconstruction of three-dimensional molecular structure from diffraction of laser-aligned molecules. Struct. Dyn. 1, 044101 (2014).

difference angular distribution from the 2D difference pattern, (ii) adding a baseline that accounts for the angular distribution of the reference pattern to the difference in angular distribution and (iii) normalizing the angular distribution. Step i is implemented by rst converting Fig. 4c,g to polar coordinates, then integrating over the radial coordinate within the FWHM of the peak. In step ii, the baseline is simulated using equations (2) and (3), corrected by the rescaling parameter that accounts for the spatial overlap between laser and electrons(0.42, obtained from the tting for Fig. 3). In step iii, the normalization used is R

p0 f a

sina da 1.

References

1. Pedersen, S., Herek, J. L. & Zewail, A. H. The validity of the diradical hypothesis: direct femtoscond studies of the transition-state structures. Science 266, 13591364 (1994).

2. Schoenlein, R. W., Peteanu, L. A., Mathies, R. A. & Shank, C. V. The rst step in vision: femtosecond isomerization of rhodopsin. Science 254, 412415 (1991).

3. Mathies, R. A., Brito Cruz, C. H., Pollard, W. T. & Shank, C. V. Direct observation of the femtosecond excited-state cis-trans isomerization in bacteriorhodopsin. Science 240, 777779 (1988).

4. Blanchet, V., Zgierski, M. Z., Seideman, T. & Stolow, A. Discerning vibronic molecular dynamics using time-resolved photoelectron spectroscopy. Nature 401, 5254 (1999).

5. Bisgaard, C. Z. et al. Time-resolved molecular frame dynamics of xed-in-space CS2 molecules. Science 323, 14641468 (2009).

6. Deb, S. & Weber, P. M. The ultrafast pathway of photon-induced electrocyclic ring-opening reactions: the case of 1,3-cyclohexadiene. Annu. Rev. Phys. Chem. 62, 1939 (2011).

7. Gessner, O. et al. Femtosecond multidimensional imaging of a molecular dissociation. Science 311, 219222 (2006).

8. Milne, C. J., Penfold, T. J. & Chergui, M. Recent experimental and theoretical developments in time-resolved X-ray spectroscopies. Coord. Chem. Rev. 277, 4468 (2014).

9. Kpper, J. et al. X-ray diffraction from isolated and strongly aligned gas-phase molecules with a free-electron laser. Phys. Rev. Lett. 112, 83002 (2014).

10. Ihee, H. et al. Ultrafast x-ray diffraction of transient molecular structures in solution. Science 309, 12231227 (2005).

11. Siders, C. W. et al. Detection of nonthermal melting by ultrafast X-ray diffraction. Science 286, 13401342 (1999).

12. Minitti, M. P. et al. Imaging molecular motion: femtosecond X-ray scattering of an electrocyclic chemical reaction. Phys. Rev. Lett. 114, 255501 (2015).

13. Gao, M. et al. Mapping molecular motions leading to charge delocalization with ultrabright electrons. Nature 496, 343346 (2013).

14. Siwick, B. J., Dwyer, J. R., Jordan, R. E. & Miller, R. J. D. An atomic-level view of melting using femtosecond electron diffraction. Science 302, 13821385 (2003).

15. Ihee, H. et al. Direct imaging of transient molecular structures with ultrafast diffraction. Science 291, 458462 (2001).

16. Srinivasan, R., Feenstra, J. S., Park, S. T., Xu, S. & Zewail, A. H. Dark structures in molecular radiationless transitions determined by ultrafast diffraction. Science 307, 558563 (2005).

17. Blaga, C. I. et al. Imaging ultrafast molecular dynamics with laser-induced electron diffraction. Nature 483, 194197 (2012).

18. Meckel, M. et al. Laser-induced electron tunneling and diffraction. Science 320, 14781482 (2008).

19. Krasniqi, F. et al. Imaging molecules from within: ultrafast angstrm-scale structure determination of molecules via photoelectron holography using free-electron lasers. Phys. Rev. A At. Mol. Opt. Phys. 81, 111 (2010).

20. Srinivasan, R., Lobastov, V. A., Ruan, C.-Y. & Zewail, A. H. Ultrafast electron diffraction (UED): a new development for the 4D determination of transient molecular structures. Helv. Chim. Acta 86, 17631838 (2003).

21. Yang, J., Beck, J., Uiterwaal, C. J. & Centurion, M. Imaging of alignment and structural changes of carbon disulde molecules using ultrafast electron diffraction. Nat. Commun. 6, 8172 (2015).

22. Sciaini, G. & Miller, R. J. D. Femtosecond electron diffraction: heralding the era of atomically resolved dynamics. Rep. Prog. Phys. 74, 096101 (2011).23. Miller, R. J. D. Mapping atomic motions with ultrabright electrons: The Chemists Gedanken Experiment Enters the Lab Frame. Annu. Rev. Phys. Chem. 65, 583604 (2014).

24. Levine, B. G. & Martnez, T. J. Isomerization through conical intersections. Annu. Rev. Phys. Chem. 58, 613634 (2007).

25. Kochman, M. A., Tajti, A., Morrison, C. A. & Miller, R. J. D. Early events in the nonadiabatic relaxation dynamics of 4-(N, N -dimethylamino)benzonitrile.J. Chem. Theory Comput. 11, 11181128 (2015).26. Ischenko, A. A., Ewbank, J. D. & Lothar, S. Structural kinetics by stroboscopic gas electron diffraction Part 1. Time-dependent molecular intensities of dissociative states. J. Mol. Struct. 320, 147158 (1994).

27. Elsayed-Ali, H. E. & Mourou, G. a. Picosecond reection high-energy electron diffraction. Appl. Phys. Lett. 52, 103104 (1988).

8 NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11232 ARTICLE

61. Holmegaard, L. et al. Control of rotational wave-packet dynamics in asymmetric top molecules. Phys. Rev. A 75, 051403 (2007).

62. Hagena, O. F. Nucleation and growth of clusters in expanding nozzle ows. Surface Sci. Lett. 106, 101116 (1981).

63. Chatelain, R. P., Morrison, V. R., Godbout, C. & Siwick, B. J. Ultrafast electron diffraction with radio-frequency compressed electron pulses. Appl. Phys. Lett. 101, 26 (2012).

64. Van Oudheusden, T. et al. Compression of subrelativistic space-charge-dominated electron bunches for single-shot femtosecond electron diffraction. Phys. Rev. Lett. 105, 264801 (2010).

65. Waldecker, L., Bertoni, R. & Ernstorfer, R. Compact femtosecond electron diffractometer with 100 keV electron bunches approaching the single-electron pulse duration limit. J. Appl. Phys. 117, 1310981901 (2015).

66. Gerbig, C., Senftleben, A., Morgenstern, S., Sarpe, C. & Baumert, T. Spatio-temporal resolution studies on a highly compact ultrafast electron diffractometer. N. J. Phys. 17, 043050 (2015).

67. Sciaini, G. et al. Electronic acceleration of atomic motions and disordering in bismuth. Nature 458, 5659 (2009).

68. Li, R. K., Musumeci, P., Bender, H. A., Wilcox, N. S. & Wu, M. Imaging single electrons to enable the generation of ultrashort beams for single-shot femtosecond relativistic electron diffraction. J. Appl. Phys. 110, 074512 (2011).

69. Beye, M. et al. X-ray pulse preserving single-shot optical cross-correlation method for improved experimental temporal resolution. Appl. Phys. Lett. 100, 15 (2012).

70. Battaglia, M. et al. Characterisation of a CMOS active pixel sensor for use in the TEAM microscope. Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 622, 669677 (2010).

Acknowledgements

We thank SLAC management for the strong support. The technical support by SLAC Accelerator Directorate, Technology Innovation Directorate, LCLS Laser Science and Technology division and Test Facilities Department is gratefully acknowledged.

This work was supported in part by the US Department of Energy (DOE) Contract No. DE-AC02-76SF00515, DOE Ofce of Basic Energy Sciences Scientic User Facilities

Division, the SLAC UED/UEM Initiative Program Development Fund and by the AMOS program within the Chemical Sciences, Geosciences, and Biosciences Division of the Ofce of Basic Energy Sciences, Ofce of Science, US Department of Energy. J.Y. and M.C. were partially supported by the US Department of Energy Ofce of Science, Ofce of Basic Energy Sciences under Award Number DE-SC0003931. M.S.R. was supported by the National Science Foundation EPSCoR RII Track-2CA Award No. IIA-1430519.

Author contributions

J.Y., M.G., T.V., M.S.R., R.L., X.S., T.G., F.W., S.W. and X.W. carried out the experiments. N.H., R. C., J.C., I.M., S.V. and A.F. developed the laser system. M.G. and J.Y. constructed the setup for gas phase experiments. C.H., K.J., A.R. and C.Y. helped on experimental setup. J.Y. performed the data analysis and simulations. The experiment was conceived by M.G., M.C. and X.W. The manuscript was prepared by J.Y., M.C., M.S.R. and M.G. with discussion and improvements from all authors. M.C. and X.W. supervised the work.

Additional information

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Yang, J. et al. Diffractive imaging of a rotational wavepacket in nitrogen molecules with femtosecond megaelectronvolt electron pulses. Nat. Commun. 7:11232 doi: 10.1038/ncomms11232 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 7:11232 | DOI: 10.1038/ncomms11232 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9