Dynamic Optical Tuning of Interlayer Interactions in ... - ACS Publications

Nov 9, 2017 - Basic Energy Sciences, under Contract No. DE-SC0012509. D.J.S. and J.S. acknowledge support from the Department of. Energy, Office of Sc...
4 downloads 11 Views 942KB Size
Subscriber access provided by READING UNIV

Communication

Dynamic optical tuning of interlayer interactions in the transition metal dichalcogenides Ehren M Mannebach, Clara Nyby, Friederike Ernst, Yao Zhou, John Tolsma, Yao Li, MengJu Sher, I-Cheng Tung, Hua Zhou, Qi Zhang, Kyle Seyler, Genevieve Clark, Yu Lin, Diling Zhu, James M. Glownia, Michael Edmund Kozina, Sanghoon Song, Silke Nelson, Apurva Mehta, Yifei Yu, Anupum Pant, Ozgur Burak Aslan, Archana Raja, Yinsheng Guo, Anthony DiChiara, Wendy L Mao, Linyou Cao, Sefaattin Tongay, Jifeng Sun, David J. Singh, Tony F. Heinz, Xiaodong Xu, Allan H MacDonald, Evan J. Reed, Haidan Wen, and Aaron M. Lindenberg Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.7b03955 • Publication Date (Web): 09 Nov 2017 Downloaded from http://pubs.acs.org on November 10, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Reed, Evan; Stanford University, Wen, Haidan; Argonne National Laboratory, Lindenberg, Aaron; Stanford University, Department of Materials Science and Engineering; SLAC National Accelerator Laboratory

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Dynamic optical tuning of interlayer interactions in the transition metal dichalcogenides Ehren M. Mannebach1, Clara Nyby2, Friederike Ernst3,6, Yao Zhou1, John Tolsma4, Yao Li3, Meng-Ju Sher5, I-Cheng Tung7, Hua Zhou7, Qi Zhang7, Kyle L. Seyler8, Genevieve Clark8, Yu Lin5,9, Diling Zhu10, James M. Glownia10, Michael E. Kozina3,5, Sanghoon Song10, Silke Nelson10, Apurva Mehta10, Yifei Yu11, Anupum Pant12, Ozgur Burak Aslan3, Archana Raja3,6, Yinsheng Guo13, Anthony DiChiara7, Wendy Mao5,6,9, Linyou Cao11, Sefaattin Tongay12, Jifeng Sun14, David J. Singh14, Tony F. Heinz3,5,6, Xiaodong Xu8, Allan H. MacDonald4, Evan Reed1,6, Haidan Wen7, Aaron M. Lindenberg1,5,6*

1

Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305

2

Department of Chemistry, Stanford University, Stanford, CA 94305

3

Department of Applied Physics, Stanford University, Stanford, CA 94305

4

Department of Physics, University of Texas, Austin, Austin, TX 78712

5

Stanford Institute for Materials and Energy Sciences, SLAC National Accelerator Laboratory,

Menlo Park, CA 94025 6

PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA 94025

7

Advanced Photon Source, Argonne National Laboratory, Argonne, IL 64309

8

Department of Physics, University of Washington, Seattle, WA 98195

9

Department of Geological Sciences, Stanford University, Stanford, CA 94305

10

SLAC National Accelerator Laboratory, Menlo Park, CA 94025

1 ACS Paragon Plus Environment

Page 2 of 19

Page 3 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

11

Department of Materials Science and Engineering, North Carolina State University, Raleigh,

North Carolina, 27695 12

School for Engineering of Matter, Transport, and Energy, Arizona State University, Tempe,

AZ 85287 13

Department of Chemistry, Columbia University, New York, NY 10027

14

Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211-7010

Corresponding author: A.M. Lindenberg; [email protected]; phone: 650-725-2640

ABSTRACT:

Modulation of weak interlayer interactions between quasi-two-dimensional

atomic planes in the transition metal dichalcogenides (TMDCs) provides avenues for tuning their functional properties. Here we show that above-gap optical excitation in the TMDCs leads to an unexpected large-amplitude, ultrafast compressive force between the two-dimensional layers, as probed by in-situ measurements of the atomic layer spacing at femtosecond time resolution. We show that this compressive response arises from a dynamic modulation of the interlayer van der Waals interaction and that this represents the dominant light-induced stress at low excitation densities. A simple analytic model predicts the magnitude and carrier density dependence of the measured strains. This work establishes a new method for dynamic, non-equilibrium tuning of correlation-driven dispersive interactions and of the optomechanical functionality of TMDC quasi-two-dimensional materials.

2 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

KEYWORDS: 2D materials, interlayer van der Waals interactions, ultrafast, Casimir effect, femtosecond x-ray scattering

Transition metal dichalcogenide materials exhibit a host of novel optoelectronic functionalities in the monolayer limit, including the emergence of a direct band gap,1 piezoelectric response,2 nonlinear optical properties,3 and valley selection rules.4 In the few-tomultilayer limit, the relatively weak interlayer coupling provides new means of modulating and controlling their electronic band structure, optoelectronic properties, and atomic-scale structure through intercalation, pressure, temperature, and stacking of layers.5-8 Van der Waals forces that determine the interlayer spacing arise from correlated fluctuating dipoles coupled by electromagnetic fields between the atomic planes and are closely related to Casimir forces, in which dressed photon states lead to attractive forces driven by quantum fluctuations. A general theory by Lifshitz describes both effects.9,10 In particular, this interaction depends upon how fields are screened at each layer (or equivalently, on the frequency-dependent layer conductivity). A number of experimental and theoretical efforts in recent years have focused on attempts to tune and engineer the Casimir interaction through modification of material dielectric functions, with both attractive and repulsive interactions possible in the retarded and nonretarded limits.11-17 Mechanisms by which the Casimir interaction can be tuned18,19 and its behavior in the non-equilibrium or excited state limit, where direct comparison between experiment and theory is difficult to achieve,20 are topics of current interest. In this work, femtosecond-resolution x-ray scattering techniques are used to record optically-induced changes in interlayer spacing as a direct and in-situ probe of the effective time-dependent interlayer pressure. In contrast to expected thermally-induced expansive effects, we find that carrier

3 ACS Paragon Plus Environment

Page 4 of 19

Page 5 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

excitation leads to large-amplitude compressive strains, with this response representing the dominant contribution to the interlayer coupling at low light excitation densities. We show that this effect is observable within most of the prototypical TMDC materials known today, including MoS2, MoSe2, MoTe2, and WSe2, with peak observed pressures approaching 0.1 GPa. We further show that a simple analytical model based on a carrier-induced modulation of the dielectric function quantitatively matches the experimental trends and results and is supported by a more detailed many-body treatment. Our study thus demonstrates ultrafast all-optical modification of the near-field Casimir/van der Waals interaction and of the associated interlayer spacing in van der Waals materials. Measurements were carried out at the Linac Coherent Light Source, the Stanford Synchrotron Radiation Laboratory, and the Advanced Photon Source using focused hard x-ray beams to probe single TMDC domains/flakes. See the Supplementary Information for further details regarding the different experimental setups. Figure 1a shows a schematic of the pumpprobe geometry used in this study in which an above-bandgap optical pulse from a femtosecond laser excites a single-crystal TMDC flake which is then probed by a femtosecond x-ray pulse. Figure 1b shows a series of diffraction images from the (004) Bragg peak of MoS2 at several different pump-probe time delays for an absorbed pump fluence of 20 mJ/cm2 (see Supplementary Movie 1 for the full sequence). Analysis of the center-of-mass motion of the diffracted spot enables extraction of the interlayer d-spacing with quantitative details on this analysis in the Supplementary Information.21 The thickness of each flake was measured from the thickness fringes—seen as satellite peaks in the diffraction patterns in Fig. 1 and Fig. S7—and was confirmed by atomic force microscopy.

4 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 19

Figure 2a shows the extracted time-dependent interlayer strain along the c-axis, ηc, in a 50 nm MoS2 sample for absorbed fluences of 0.3 mJ/cm2 and 20 mJ/cm2. In these measurements, the strain is directly extracted from the shift in peak position:  = −

 

. Immediately after

excitation, instead of thermally expanding, the MoS2 layers initially and transiently compress. In Fig. 1b this effect appears at short times as a deflection of the diffraction peak center of mass to higher momentum transfer Qz at a probe delay time of ~10 ps. We observe a similar compressive effect in MoSe2, MoTe2, and WSe2 (Supplementary Information, Figure S8). In contrast, at low fluences, we find that this compressive response dominates over any thermal effect (Fig. 2a, bottom panel). Also observed is a coherent breathing mode response where the period of oscillation is determined by the acoustic round trip time in the flake.22 Under the assumption that the mechanically exfoliated domain is weakly bound to the surface (determined by the van der Waals gap between the flake and substrate), this period is simply 2d/v, where d is the overall thickness of the flake, and v is the longitudinal sound velocity in the cross-plane direction.23 In MoS2, we extract an interlayer speed of sound of ~2500 m/s, in agreement with previous measurements obtained from the bulk C33 elastic constant.24 We note that prior theoretical studies have pointed towards the possibility for compressive features in quasi-2D systems under direct vibrational excitation25 and previous experimental studies have indicated possibilities for compressive responses in other van der Waals bonded materials at high excitation fluences.26-28 The fluence dependence of the compressive effect (Figure 2b) is strongly sublinear, even down to the lowest measured pump fluences of ~5 µJ/cm2, corresponding to induced carrier densities of ~1018 cm-3 and extending over approximately three orders of magnitude in fluence. At higher pump fluences, the initial transient compressive feature is dominated by thermal expansion effects, which scale linearly with pump fluence but turn on after the initial 5 ACS Paragon Plus Environment

Page 7 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

compressive ones. This is consistent with an electronic, nearly instantaneous origin for the compressive stress, in contrast to a thermal stress which would develop on electron-phonon coupling time-scales of order a few picoseconds.29,30 Further support for an electronic origin of the compressive feature is found by comparing the lifetime of the compressive response at low fluences, a few hundred picoseconds in MoS2, to time-domain terahertz spectroscopic measurements of the lifetime of photoexcited carriers, where good agreement is obtained between the two (Figure 3a). Both the THz studies and the x-ray studies show a lifetime which depends on the carrier density, as would be expected for Auger-type relaxation processes (Supplementary Information Fig. S10). In contrast, multilayer ReS2—a direct band gap material where the free carrier lifetime is much shorter (~1.5 ps) such that the applied electronic stress turns off almost instantaneously—does not exhibit the nonlinear compressive response observed in the Mo- and W-dichalcogenides (Figure 3b). Taken together, these observations consistently imply a carrier-driven electronic origin for the compressive response. To obtain an understanding of the electronically-driven compressive strain associated with a modification of the interlayer bonding forces, we consider from a microscopic perspective how the attractive van der Waals forces responsible for this interaction depend on carrier density. Lifshitz developed a theory in which the force between two thick slabs separated by a distance l is related to their frequency-dependent dielectric functions. The Casimir interaction represents a special case of this theory valid for the limit of large l and perfect metal slabs. In the limit of small spacing between two plates, the effective pressure predicted by this theory can be written as

  − 1 ℏ ∆ =      , 1 8   + 1

6 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 19

where ε is the frequency dependent dielectric function evaluated at complex values, ℏ is the reduced Planck constant, and the integral is over all frequencies. 9,12 Although this model is not intended to be accurate for atomic scale separations, comparison with more detailed models based on electronic structure calculations suggests that it captures the essential physics of the response for some semiconducting 2D materials.31 As a simple approximation to simulate our experimental conditions, we assume a Drude model for the dielectric function,  = 1 + 

, where $% is the plasma frequency (proportional to the square root of the photo-excited

!!"#

carrier density n), and γ is the Drude damping term. The time-domain THz conductivity measurements support this assumed Drude-like response of the photoexcited electron gas.32 Previous studies investigating modifications to the Casimir force driven by doping or other modifications of the dielectric constants have used similar approaches.

11-13

In this

approximation, the integral in Equation 1 can be evaluated analytically. The Supplementary Information has further details on this calculation. In the limit $% ≫ ', which is relevant for our experimental conditions, this effective pressure can be written as

∆( ≈



%$32√2 

=

ℏ-

32.2/0 

√( 2

where m is the effective mass of the electron gas, and ϵ0 is the vacuum permittivity. Evaluation of this expression based on the Mo-to-Mo layer separation = 6.2 Å gives pressures of order 10 MPa at a carrier density of n = 1021 cm-3. To compare with the experimental data, we make use of the known absorption coefficients of the materials to convert the applied pump fluence into a carrier concentration. We then calculate the strain induced by the electronic force using the interlayer bulk modulus. In this approximation, only nearest layer interactions are considered, a 7 ACS Paragon Plus Environment

Page 9 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

reasonable approximation given the 1/l3 scaling of the force. As can be seen in Fig. 2b (dashed line) for MoS2, good quantitative agreement in both magnitude and carrier-density dependence is achieved despite the simplicity of this model. Moreover, the n1/2 scaling explains why the electronic effect dominates over the usual thermal response (scaling linearly with n) at low carrier densities. More detailed first principles models, including both interband and intraband excitations and the full dielectric function, give similar behavior with carrier density, as described in the Supplementary Information. This suitably defined Lifshitz model has been found to give vdW interaction energy to within 8–20% of the adiabatic-connection fluctuationdissipation theorem within random phase approximation (ACFDT-RPA) calculations for a variety of layered semiconducting heterostructures.33 DFT-based approaches modeling the vdW interaction in MoSe2 using the method of reference 33 further support the arguments made here, showing that the attractive contribution arising from the electronic polarizability dominates over expansive effects that would arise from excitation into antibonding orbitals.34 Finally, a manybody calculation described in the Supplementary Information shows from a different perspective how quantum fluctuations between optically-excited electrons and holes lead to compressive effects. It also quantifies the role of near field effects neglected in the Lifshitz approximation that are expected to become important when the layer separation is small compared to the carrier Fermi wavelength. This more detailed treatment (see in particular equation S19) yields Equation 2 above at small Fermi wavelengths (up to a factor of 2/π) and further confirms the basic physical model presented above. Further support for this model can be obtained by considering the balance between the electronic and thermal contributions. One can write the net light-induced, time-dependent interlayer stress as the sum of two terms representing the electronic (σe) and thermal (σth) 8 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 19

pressures, with 45 = −6√( and 478 = 9 6 ∆: = 9 6 ℏ$ − ;< (/9. Here, α is the prefactor in Equation 2, C33 is the interlayer elastic constant, αc is the c-axis thermal expansion coefficient, ∆T is the photo-induced temperature jump, ℏ$ − ;< is the photon energy minus the band gap (i.e. the excess photon energy), and C is the specific heat. Assuming the sample is homogeneously pumped, the time-dependent stress at short times can then be written as 4> ∝ −6 √(@> + A(@> 1 − -

B

7 CD

, 3

where A = 9 6 ℏ$ − ;< ⁄9 , H(t) is the Heaviside step function representing an instantaneous turn-on of the free carrier population, and τep is the effective electron-phonon coupling time which gives rise to the observed delayed onset of the thermal response. This simple model accurately predicts the experimentally observed cross-over fluence, FT, at which the thermal and electronic effects become comparable in magnitude, determined quantitatively in terms of the G ℏ

above coefficients by F = H

GI

~ 100 µJ/cm2, where 6J is the optical absorption coefficient.

Following this time-dependent stress, the film responds by launching strain waves from the free and buried interfaces of the form K, > =

L

MN

4> − K/O.35 A more detailed model presented in

the Supplementary Information predicts how these strain waves are reflected at both the free and buried interfaces, with the observed period 2d/v signifying a quasi-free surface at the substrate interface, in agreement with recent studies of exfoliated flakes van der Waals bonded to a substrate.36 The results of this second model are shown in Figure 4 with good agreement to experimental results. We note that additional modifications to the repulsive part of the interlayer potential and to interlayer orbital hybridization may also be occurring under carrier excitation; separation of

9 ACS Paragon Plus Environment

Page 11 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

these effects from the modulations discussed above is the subject of future work. Carrier/exciton diffusion processes are not likely to play a significant role in the observed response given the quasi-homogeneous excitation conditions used here. Additionally, one can rule out the role of radiation pressure effects or other related stresses which only act during the presence of the incident electromagnetic wave, in contrast to the above considered carrier-induced stresses which stay on for time-scales corresponding to recombination/trapping time-scales. In summary, we demonstrate that light can be used to tune interlayer van der Waals/Casimir interactions in quasi-two-dimensional transition metal dichalcogenides on ultrashort time scales. This arises from a strong carrier-driven modulation in the polarizability, which leads to attractive interlayer bonding, increasing this interaction strength. The square rootlike scaling with fluence corresponds to an effectively divergent response function (with increasing slope at low carrier densities) and potentially enables new opportunities for nanooptomechanical devices driven by low light excitation densities.25 For example, in a TMDC structure with volume of order a few nm3, a single absorbed photon is predicted to generate an effective compressive strain of ~0.1–1% and enable new optical means for modifying band structures and associated electronic properties. In addition to the pulsed excitation used here, an additional direction for exploration may involve photoexcitation in the continuous wave limit for tuning 2D layer interactions. Optical control of the van der Waals interaction may also open up new means of optically modulating exciton transfer processes governed by resonant energy transfer.37 Finally, measurements in a tailored heterostructure architecture with stacked layers of different dielectric functions may enable the observation of the converse effect, namely carrierinduced repulsive interlayer interactions, opening new means for dynamic optical tuning of nanoscale stiction effects.

10 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 19

Methods. Samples were prepared via mechanical exfoliation from bulk crystals onto c-axis sapphire substrates. The thickness of each flake was measured from the thickness fringes—seen as satellite peaks in the diffraction patterns in Fig. 1 and Fig. S7—and was confirmed by atomic force microscopy. The thicknesses of the flakes investigated were between 30–150 nm. For all experiments, the samples were pumped above their optical band gaps. The repetition rates of the pump lasers in the time-resolved diffraction measurements were matched to the repetition rate of the relevant x-ray source; the THz transmission measurements were performed on a separate apparatus with a 960 Hz laser chopped at 480 Hz. The optical pump parameters were 40 fs pulses at 650 nm (1.9 eV) from an optical parametric amplifier driven by an amplified a Ti:sapphire laser and at the Linac Coherent Light Source (LCLS), 2 ps pulses at 400 nm (3.1 eV) from a frequency-doubled Ti:sapphire laser at the Advanced Photon Source (APS), 500 fs pulses at 515 nm (2.4 eV) from a frequency-doubled Yb-doped fiber laser at the Stanford Synchrotron Radiation Lightsource (SSRL), and 100 fs pulses at 400 nm from a frequency-doubled Ti:sapphire laser for the THz transmission measurements. The polarization of the pump light was in the plane of the x-ray scattering. The pump beam was roughly collinear with the probe beam for the LCLS, SSRL, and THz transmission measurements and was perpendicular to the probe beam for the APS experiments. Time-resolved diffraction experiments were carried out at the X-Ray Pump-Probe (XPP) instrument at LCLS, at beamline 14-ID-B at APS, and at beamline 10-2 at SSRL. All measurements were performed in a vertically-scattering geometry. The LCLS experiments used 9.5 keV photons at a repetition rate of 120 Hz and the Cornell-SLAC Pixel Array Detector (CS-

11 ACS Paragon Plus Environment

Page 13 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

PAD). The APS experiments used 12 keV photons at a repetition rate of 1 kHz and a gated Pilatus 100K detector. The SSRL experiments used 15 keV photons at a repetition rate of 1.28 MHz and a gated Pilatus 100K detector. All x-ray probes were monochromatized to a bandwidth of a few eV and focused onto the sample to probe a single domain TMDC flake. THz pulses were generated via a two-color laser-induced gas plasma. The THz beam transmitted through the samples was detected by electro-optic sampling in a 1 mm-thick ZnTe crystal. A lock-in amplifier was used to record the differential signal changes of the transmitted THz beam induced by photoexcitation.

These measurements average over many flakes in

contrast to the x-ray probes.

ASSOCIATED CONTENT Supporting Information: First principles modeling of the vdW interaction; electron-hole fluid superlattice model of compressive strain; Analytical Lifshitz model; Extraction of strain from time-dependent diffraction images; Response of other TMDCs; THz conductivity measurements; Acoustic model. AUTHOR INFORMATION Corresponding Author *Email: [email protected] Author Contributions E.M., C.N., F.E, M.S., I.T., H.Z., Q.Z., K.S., Y.L., D.Z., J.G., M.K., S.S., S.N., A.M., A.D., W.M., T.H., X.X., H.W. and A.L. carried out the experiments. Y.Z., J.T., Y.L., A.H.M., E.R., J.S., D.J.S. and A.L. carried out the theoretical work. K.S., G.C., Y.Y., A.P., O.A., A.R., Y.G., L.C., S.T., and X.X. were responsible for sample growth. A.L. conceived the experiment. Notes The authors declare no competing financial interest. ACKNOWLEDGMENTS

12 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 19

This work was supported by the Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division under contract DE-AC02-76SF00515. Use of the Linac Coherent Light Source (LCLS), SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515. Use of the Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC0276SF00515. This research used resources of the Advanced Photon Source, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. Use of BioCARS was also supported by the National Institute of General Medical Sciences of the National Institutes of Health under grant number R24GM111072. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The time-resolved set-up at Sector 14 was funded in part through a collaboration with Philip Anfinrud (NIH/NIDDK). C.N. acknowledges support from the NSF through a Graduate Research Fellowship (DGE-114747). F.E. gratefully acknowledges Grant LPDS 2013-13 from the German National Academy of Sciences Leopoldina. I.T., Q.Z., K.S., G.C., X. X., H.W. acknowledge the support from the Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-SC0012509. D.J.S. and J.S. acknowledge support from the Department of Energy, Office of Science, Computational Materials Science Program, MAGICS Center, award DE-SC0014607. REFERENCES (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14)

Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Phys. Rev. Lett. 2010, 105, 136805. Zhu, H. Nat. Nanotech 2014, 10, 151. Li, Y.; Rao, Y.; Mak, K. F.; You, Y.; Wang, S.; Dean, C. R.; Heinz, T. F. Nano Lett. 2013, 13, 3329–3333. Mak, K. F.; He, K.; Shan, J.; Heinz, T. F. Nat. Nanotech 2012, 7, 494–498. Nayak, A. P.; Bhattacharyya, S.; Zhu, J.; Liu, J.; Wu, X.; Pandey, T.; Jin, C.; Singh, A. K.; Lin, J.-F.; Akinwande, D. Nat. Commun. 2014, 5, 3731. Zhang, H.; Yuan, H.; Wang, S.; Lin, Y.; Zeng, Q.; Xu, G.; Liu, Z.; Solanki, G. K.; Patel, K. D.; Cui, Y.; Hwang, H. Y.; Mao, W. L.; Zhao, Z. Nat. Commun. 2015, 6, 7312. Xiong, F.; Wang, H.; Liu, X.; Sun, J.; Brongersma, M.; Pop, E.; Cui, Y. Nano Lett. 2015, 15, 6777–6784. Tongay, S.; Zhou, J.; Ataca, C.; Lo, K.; Matthews, T. S.; Li, J.; Grossman, J. C.; Wu, J. Nano Lett. 2012, 12, 5576–5580. Lifshitz, E. M. Soviet Physics JETP 1956, 2, 73. Rodriguez, A. W.; Capasso, F.; Johnson, S. G. Nature Photon. 2011, 5, 211. Torricelli, G.; van Zwol, P. J.; Shpak, O.; Palasantzas, G.; Svetovoy, V. B.; Binns, C.; Kooi, B. J.; Jost, P.; Wuttig, M. Adv. Funct. Mater. 2012, 22, 3729–3736. Galkina, E. G.; Ivanov, B. A.; Savel’ev, S.; Yampol’skii, V. A.; Nori, F. Phys. Rev. B 2009, 80, 125119–11. Boström, M.; Sernelius, B. E. Phys. Rev. A 2012, 85, 012508–5. Chen, F.; Klimchitskaya, G. L.; Mostepanenko, V. M.; Mohideen, U. Phys. Rev. B 2007, 76, 035338–15. 13 ACS Paragon Plus Environment

Page 15 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

(15) (16) (17) (18) (19) (20) (21)

(22) (23) (24) (25) (26) (27) (28)

(29)

(30) (31) (32) (33) (34) (35) (36) (37)

Arnold, W.; Hunklinger, S.; Dransfeld, K. Phys. Rev. B 1979. Tkatchenko, A. Adv. Funct. Mater. 2014, 25, 2054–2061. Woods, L. M.; Dalvit, D. A. R.; Tkatchenko, A.; Rodriguez-Lopez, P.; Rodriguez, A. W.; Podgornik, R. Rev. Mod. Phys. 2016, 88, 045003–045048. Iannuzzi, D.; Lisanti, M. Proc. Natl. Acad. Sci. 2004, 101, 4019. Klimchitskaya, G. L.; Mohideen, U.; Mostepanenko, V. M. Rev. Mod. Phys. 2009, 81, 1827–1885. Cohen, A. E.; Mukamel, S. Phys. Rev. Lett. 2003, 91 (23), 233202–233204. Kozina, M.; Hu, T.; Wittenberg, J. S.; Szilagyi, E.; Trigo, M.; MIller, T. A.; Uher, C.; Damodaran, A.; Martin, L.; Mehta, A.; Corbett, J.; Safranek, J.; Reis, D. A.; Lindenberg, A. M. Struct. Dyn. 2014, 1, 034301. Ge, S.; Liu, X.; Qiao, X.; Wang, Q.; Xu, Z.; Qiu, J.; Tan, P.-H.; Zhao, J.; Sun, D. Sci. Rep. 2014, 4, 5722. Antonelli, G. A.; Perrin, B.; Daly, B. C.; Cahill, D. G. MRS Bull. 2006, 31, 607. Peelaers, H.; Van de Walle, C. G. J. Phys. Chem. C 2014, 118, 12073. Miyamoto, Y.; Zhang, H.; Miyazaki, T.; Rubio, A. Phys. Rev. Lett. 2015, 114, 116102– 116105. Carbone, F.; Baum, P.; Rudolf, P.; Zewail, A. H. Phys. Rev. Lett. 2008, 100, 035501– 035504. Raman, R. K.; Murooka, Y.; Ruan, C.-Y.; Yang, T.; Berber, S.; Tománek, D. Phys. Rev. Lett. 2008, 101, 077401–077404. Ferguson, K. R.; Bucher, M.; Gorkhover, T.; Boutet, S.; Fukuzawa, H.; Koglin, J. E.; Kumagai, Y.; Lutman, A.; Marinelli, A.; Messerschmidt, M.; Nagaya, K.; Turner, J.; Ueda, K.; Williams, G. J.; Bucksbaum, P. H.; Bostedt, C. Sci. Adv. 2016, 2, e1500837– e1500837. Mannebach, E. M.; Li, R.; Duerloo, K.-A.; Nyby, C.; Zalden, P.; Vecchione, T.; Ernst, F.; Reid, A. H.; Chase, T.; Shen, X.; Weathersby, S.; Hast, C.; Hettel, R.; Coffee, R.; Hartmann, N.; Fry, A. R.; Yu, Y.; Cao, L.; Heinz, T. F.; Reed, E. J.; Dürr, H. A.; Wang, X.; Lindenberg, A. M. Nano Lett. 2015, 15, 6889–6895. Waldecker, L.; Bertoni, R.; Hübener, H.; Brumme, T.; Vasileiadis, T.; Zahn, D.; Rubio, A.; Ernstorfer, R. Phys. Rev. Lett. 2017, 119, 036803–036805. Zhou, Y.; Pellouchoud, L. A.; Reed, E. J. 2016. Strait, J. H.; Nene, P.; Rana, F. Phys. Rev. B 2014, 90, 245402–245409. Zhou, Y.; Pellouchoud, L. A.; Reed, E. J. 2D Mater. 2017, 4, 025005–025014. Tkatchenko, A.; Scheffler, M. Phys. Rev. Lett. 2009, 102, 073005–4. Nicoul, M.; Shymanovich, U.; Tarasevitch, A.; Linde, von der, D.; Sokolowski-Tinten, K. Appl. Phys. Lett. 2011, 98, 191902. Beardsley, R.; Akimov, A. V.; Greener, J. D. G.; Mudd, G. W.; Sandeep, S.; Kudrynskyi, Z. R.; Kovalyuk, Z. D.; Patanè, A.; Kent, A. J. Sci. Rep. 2016, 6, 26970. Rowland, C. E.; Fedin, I.; Zhang, H.; Gray, S. K.; Govorov, A. O.; Talapin, D. V.; Schaller, R. D. Nat. Mater. 2015, 14, 484–489.

14 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 19

Figure 1. Schematics of the time-resolved experiment and data collection. (A) Schematic showing the optical pump/x-ray probe experimental setup wherein a TMDC flake is probed by an x-ray pulse (black) following ultrafast above-bandgap optical excitation (red). Changing the carrier density via optical excitation increases the magnitude of the attractive interlayer van der Waals forces. (B) Diffraction images of the (004) Bragg spot of MoS2 at selected probe delay times for an absorbed pump fluence of 20 mJ/cm2. The vertical direction corresponds to the momentum transfer Qz in units of Å-1. (C) The thickness of the sample can be extracted from the spacing of the thickness fringes that lie along the Qz direction:  = 2 /ΔQR .

15 ACS Paragon Plus Environment

Page 17 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 2. The photo-induced strain effects in exfoliated MoS2. (A) Interlayer strain in MoS2 following above-bandgap optical pumping. The bottom panel shows the gray shaded region of the top panel. In MoS2 (as well as MoSe2, MoTe2, and WSe2, see Fig. S3), the material initially compresses (ηc, σc < 0) following photoexcitation before thermally expanding. Echoes of the compressive response can be seen every ~40 ps in MoS2, corresponding to an acoustic breathing of the 50 nm-thick flake. The right axis shows the effective stress needed to achieve the experimentally-measured strains. The inset shows a representative histogram of the spread of strains measured at each delay point, showing the effective sensitivity of the measurement. Error bars on the measurement are indicated by the shaded regions. (B) The fluence dependence of the compressive strain in MoS2 shows a nonlinear response, scaling roughly as the square root of the carrier density/pump fluence. The black points are experimental values, and the red dashed line 16 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 19

is the result of the analytical model of the compressive strain described by Equation 1. See the Supplementary Information for discussion of the calculation of the error bars.

Figure 3. Comparison of the free-carrier lifetime with the lifetime of the optically-induced compressive effect. Normalized transient terahertz (THz) transmission following optical pumping (top) and transient interlayer strain (bottom) in MoS2 (A) and ReS2 (B); note the difference in time scales. Decreased THz transmission corresponds to an increase in the density of free carriers. The time-resolution of the x-ray measurements in (A, bottom) is ~50 ps. The inset in (B, bottom) shows the longer-time response in ReS2 at a higher pump fluence. In the top panels, the open circles are data points, and solid lines are exponential fits to the data. In MoS2, the lifetime of the compressive strain follows the response time seen in the THz transmission measurements determined by the free carrier population. In ReS2, the free carrier lifetime is not long enough to induce a measurable compressive strain; only a thermal expansive strain is seen.

17 ACS Paragon Plus Environment

Page 19 of 19

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 4. Comparison of the simulated average strain with the measured strain in MoS2. Simulated strain profiles in a 50 nm MoS2 flake with carrier densities of 2×1020 cm-3 (A) and 1022 cm-3 (B), corresponding to the pump fluences in Figure 2a. Blue shading corresponds to compression, while red corresponds to rarefaction. (C) Time-dependent average strains in the MoS2 flake, where the left and right panels correspond to the low and high carrier densities, respectively, comparing experiment to the theoretical Lifshitz model described in the text.

For TOC only:

18 ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

19 ACS Paragon Plus Environment

Page 20 of 19