Twist-controlled moire superlattices (MS) have emerged as a versatile platform in which to realize artificial systems with complex electronic spectra. Bernal-stacked bilayer graphene (BLG) and hexagonal boron nitride (hBN) form an interesting example of the MS that has recently featured a set of unexpected behaviors, such as unconventional ferroelectricity and electronic ratchet effect. Yet, the understanding of the BLG/hBN MS electronic properties has, at present, remained fairly limited. Here we develop a multi-messenger approach that combines standard magnetotransport techniques with low-energy sub-THz excitation to get insights into the properties of this MS. We show that BLG/hBN lattice alignment results in the emergence of compensated semimetals at some integer fillings of the moire bands separated by van Hove singularities where Lifshitz transition occurs. A particularly pronounced semimetal develops when 8 electrons reside in the moire unit cell, where coexisting high-mobility electron and hole systems feature a strong magnetoresistance reaching 2350 % already at B=0.25 T. Next, by measuring the THz-driven Nernst effect in remote bands, we observe valley splitting, pointing to an orbital magnetization characterized by a strongly enhanced effective g-factor of 340. Last, using THz photoresistance measurements, we show that the high-temperature conductivity of the BLG/hBN MS is limited by electron-electron umklapp processes. Our multi-facet analysis introduces THz-driven magnetotransport as a convenient tool to probe the band structure and interaction effects in vdW materials and provides a comprehension of the BLG/hBN MS.

We describe an apparatus for investigating the excited-state dissociation dynamics of mass-selected ion-molecule clusters by mass-resolving and detecting photofragment-ions and neutrals, in coincidence, using an ultrafast laser operating at high repetition rates. The apparatus comprises a source that generates ion-molecule clusters, a time-of-flight spectrometer, and a mass filter that selects the desired anions, and a linear-plus-quadratic reflectron mass spectrometer that discriminates the fragment anions after the femtosecond laser excites the clusters. The fragment neutrals and anions are then captured by two channeltron detectors. The apparatus performance is tested by measuring the photofragments: I-, CF3I-, and neutrals from photoexcitation of the ion-molecule cluster CF3I·I- using femtosecond UV laser pulses with a wavelength of 266 nm. The experimental results are compared with our ground state and excited state electronic structure calculations as well as the existing results and calculations, with particular attention to the generation mechanism of the anion fragments and dissociation channels of the ion-molecule cluster CF3I·I- in the charge-transfer excited state.

In convergent beam electron diffraction (CBED) on two-dimensional (2D) materials, the intensity distributions within the individual CBED spots map the local atomic arrangements within the probed region. In this study we demonstrate that the average intensities within the CBED spots essentially depend on the scattering parameters by a single atom, thus, offering the possibility of the direct measurement of such parameters. Scattering potential of an individual atom can be approximated by the Gaussian function and its parameters, the standard deviation and the maximal phase shift, can be recovered from the ratio of the intensities in the zero- and higher orders spots in CBED pattern. In order to demonstrate this, we simulated CBED patterns and extracted the atomic scattering parameters from such patterns. In the simulated examples, no weak phase object approximation is applied and the proposed method provides accurate results even for materials with large phase shifts up to 1.5 rad, as for example tungsten. The scattering parameters recovered from the experimental CBED patterns of graphene and twisted graphene-hBN structure show good agreement with the theoretically obtained values.

We report the modification of a gas phase ultrafast electron diffraction (UED) instrument that enables experiments with both gas and condensed matter targets, where a time-resolved experiment with sub-picosecond resolution is demonstrated with solid state samples. The instrument relies on a hybrid DC-RF acceleration structure to deliver femtosecond electron pulses on the target, which is synchronized with femtosecond laser pulses. The laser pulses and electron pulses are used to excite the sample and to probe the structural dynamics, respectively. The new system is added with capabilities to perform transmission UED on thin solid samples. It allows for cooling samples to cryogenic temperatures and to carry out time-resolved measurements. We tested the cooling capability by recording diffraction patterns of temperature dependent charge density waves in 1T-TaS2. The time-resolved capability is experimentally verified by capturing the dynamics in photoexcited single-crystal gold.

Nonlinear mechanics of solids is an exciting field that encompasses both beautiful mathematics, such as the emergence of instabilities and the formation of complex patterns, as well as multiple applications. Two-dimensional crystals and van der Waals (vdW) heterostructures allow revisiting this field on the atomic level, allowing much finer control over the parameters and offering atomistic interpretation of experimental observations. In this work, we consider the formation of instabilities consisting of radially oriented wrinkles around mono- and few-layer "bubbles" in two-dimensional vdW heterostructures. Interestingly, the shape and wavelength of the wrinkles depend not only on the thickness of the two-dimensional crystal forming the bubble, but also on the atomistic structure of the interface between the bubble and the substrate, which can be controlled by their relative orientation. We argue that the periodic nature of these patterns emanates from an energetic balance between the resistance of the top membrane to bending, which favors large wavelength of wrinkles, and the membrane-substrate vdW attraction, which favors small wrinkle amplitude. Employing the classical "Winkler foundation" model of elasticity theory, we show that the number of radial wrinkles conveys a valuable relationship between the bending rigidity of the top membrane and the strength of the vdW interaction. Armed with this relationship, we use our data to demonstrate a nontrivial dependence of the bending rigidity on the number of layers in the top membrane, which shows two different regimes driven by slippage between the layers, and a high sensitivity of the vdW force to the alignment between the substrate and the membrane.

Monolayer hexagonal boron nitride (hBN) is a 2D crystal with a unique combination of properties. It is also predicted to be piezoelectric, yet experimental evidence has been lacking. In article 1905504, Francisco Guinea, Laura Fumagalli, Konstantin S. Novoselov, and co-workers use electrostatic force microscopy and directly visualize the piezoelectric field in monolayer hBN in nonhomogeneous strain regions around bubbles and creases.

The convergent beam electron diffraction (CBED) patterns of twisted bilayer samples exhibit interference patterns in their CBED spots. Such interference patterns can be treated as off-axis holograms and the phase of the scattered waves, meaning the interlayer distance can be reconstructed. A detailed protocol of the reconstruction procedure is provided in this study. In addition, we derive an exact formula for reconstructing the interlayer distance from the recovered phase distribution, which takes into account the different chemical compositions of the individual monolayers. It is shown that one interference fringe in a CBED spot is sufficient to reconstruct the distance between the layers, which can be practical for imaging samples with a relatively small twist angle or when probing small sample regions. The quality of the reconstructed interlayer distance is studied as a function of the twist angle. At smaller twist angles, the reconstructed interlayer distance distribution is more precise and artefact free. At larger twist angles, artefacts due to the moiré structure appear in the reconstruction. A method for the reconstruction of the average interlayer distance is presented. As for resolution, the interlayer distance can be reconstructed by the holographic approach at an accuracy of ±0.5 Å, which is a few hundred times better than the intrinsic z-resolution of diffraction limited resolution, as expressed through the spread of the measured k-values. Moreover, we show that holographic CBED imaging can detect variations as small as 0.1 Å in the interlayer distance, though the quantitative reconstruction of such variations suffers from large errors.

The availability of accessible fabrication methods based on deterministic transfer of atomically thin crystals has been essential for the rapid expansion of research into van der Waals heterostructures. An inherent issue of these techniques is the deformation of the polymer carrier film during the transfer, which can lead to highly nonuniform strain induced in the transferred two-dimensional material. Here, using a combination of optical spectroscopy, atomic force, and Kelvin probe force microscopy, we show that the presence of nanometer scale wrinkles formed due to transfer-induced stress relaxation can lead to strong changes in the optical properties of MoSe2/WSe2 heterostructures and the emergence of linearly polarized interlayer exciton photoluminescence. We attribute these changes to local breaking of crystal symmetry in the nanowrinkles, which act as efficient accumulation centers for interlayer excitons due to the strain-induced interlayer band gap reduction. Surface potential images of the rippled heterobilayer samples acquired using Kelvin probe force microscopy reveal variations of the local work function consistent with strain-induced band gap modulation, while the potential offset observed at the ridges of the wrinkles shows a clear correlation with the value of the tensile strain estimated from the wrinkle geometry. Our findings highlight the important role of the residual strain in defining optical properties of van der Waals heterostructures and suggest effective approaches for interlayer exciton manipulation by local strain engineering.

Convergent beam electron diffraction is routinely applied for studying deformation and local strain in thick crystals by matching the crystal structure to the observed intensity distributions. Recently, it has been demonstrated that CBED can be applied for imaging two-dimensional (2D) crystals where a direct reconstruction is possible and three-dimensional crystal deformations at a nanometre resolution can be retrieved. Here, we demonstrate that second-order effects allow for further information to be obtained regarding stacking arrangements between the crystals. Such effects are especially pronounced in samples consisting of multiple layers of 2D crystals. We show, using simulations and experiments, that twisted multilayer samples exhibit extra modulations of interference fringes in CBED patterns, i. e., a CBED moiré. A simple and robust method for the evaluation of the composition and the number of layers from a single-shot CBED pattern is demonstrated.

When two-dimensional (2D) atomic crystals are brought into close proximity to form a van der Waals heterostructure, neighbouring crystals may influence each other's properties. Of particular interest is when the two crystals closely match and a moiré pattern forms, resulting in modified electronic and excitonic spectra, crystal reconstruction, and more. Thus, moiré patterns are a viable tool for controlling the properties of 2D materials. However, the difference in periodicity of the two crystals limits the reconstruction and, thus, is a barrier to the low-energy regime. Here, we present a route to spectrum reconstruction at all energies. By using graphene which is aligned to two hexagonal boron nitride layers, one can make electrons scatter in the differential moiré pattern which results in spectral changes at arbitrarily low energies. Further, we demonstrate that the strength of this potential relies crucially on the atomic reconstruction of graphene within the differential moiré super cell.

2D hexagonal boron nitride (hBN) is a wide-bandgap van der Waals crystal with a unique combination of properties, including exceptional strength, large oxidation resistance at high temperatures, and optical functionalities. Furthermore, in recent years hBN crystals have become the material of choice for encapsulating other 2D crystals in a variety of technological applications, from optoelectronic and tunneling devices to composites. Monolayer hBN, which has no center of symmetry, is predicted to exhibit piezoelectric properties, yet experimental evidence is lacking. Here, by using electrostatic force microscopy, this effect is observed as a strain-induced change in the local electric field around bubbles and creases, in agreement with theoretical calculations. No piezoelectricity is found in bilayer and bulk hBN, where the center of symmetry is restored. These results add piezoelectricity to the known properties of monolayer hBN, which makes it a desirable candidate for novel electromechanical and stretchable optoelectronic devices, and pave a way to control the local electric field and carrier concentration in van der Waals heterostructures via strain. The experimental approach used here also shows a way to investigate the piezoelectric properties of other materials on the nanoscale by using electrostatic scanning probe techniques.

Van der Waals heterostructures have been lately intensively studied because they offer a large variety of properties that can be controlled by selecting 2D materials and their sequence in the stack. The exact arrangement of the layers as well as the exact arrangement of the atoms within the layers, both are important for the properties of the resulting device. However, it is very difficult to control and characterize the exact position of the atoms and the layers in such heterostructures, in particular, along the vertical (z) dimension. Recently it has been demonstrated that convergent beam electron diffraction (CBED) allows quantitative three-dimensional mapping of atomic positions in three-dimensional materials from a single CBED pattern. In this study we investigate CBED in more detail by simulating and performing various CBED regimes, with convergent and divergent wavefronts, on a somewhat simplified system: a two-dimensional (2D) monolayer crystal. In CBED, each CBED spot is in fact an in-line hologram of the sample, where in-line holography is known to exhibit high intensity contrast in detection of weak phase objects that are not detectable in conventional in-focus imaging mode. Adsorbates exhibit strong intensity contrast in the zero and higher order CBED spots, whereas lattice deformation such as strain or rippling cause noticeable intensity contrast only in the first and higher order CBED spots. The individual CBED spots can thus be reconstructed as typical in-line holograms, and a resolution of 2.13 Å can in principle be achieved in the reconstructions. We provide simulated and experimental examples of CBED of a 2D monolayer crystal. The simulations show that individual CBED spots can be treated as in-line holograms and sample distributions such as adsorbates, can be reconstructed. Individual atoms can be reconstructed from a single CBED pattern provided the later exhibits high-order CBED spots. The experimental results were obtained in a transmission electron microscope (TEM) at 80 keV on free-standing monolayer hBN containing adsorbates. Examples of reconstructions obtained from experimental CBED patterns at a resolution of 2.7 Å are shown. CBED technique can be potentially useful for imaging individual biological macromolecules, because it provides a relatively high resolution and does not require additional scanning procedure or multiple image acquisitions and therefore allows minimizing the radiation damage.

The van der Waals heterostructures, which explore the synergetic properties of 2D materials when assembled into 3D stacks, have already brought to life a number of exciting phenomena and electronic devices. Still, the interaction between the layers in such assembly, possible surface reconstruction, and intrinsic and extrinsic defects are very difficult to characterize by any method, because of the single-atomic nature of the crystals involved. Here we present a convergent beam electron holographic technique which allows imaging of the stacking order in such heterostructures. Based on the interference of electron waves scattered on different crystals in the stack, this approach allows one to reconstruct the relative rotation, stretching, and out-of-plane corrugation of the layers with atomic precision. Being holographic in nature, our approach allows extraction of quantitative information about the 3D structure of the typical defects from a single image covering thousands of square nanometers. Furthermore, qualitative information about the defects in the stack can be extracted from the convergent diffraction patterns even without reconstruction, simply by comparing the patterns in different diffraction spots. We expect that convergent beam electron holography will be widely used to study the properties of van der Waals heterostructures.

Resonance magnetic tunneling in heterostructures formed by graphene single sheets separated by a hexagonal boron nitride barrier and bounded by two gates has been investigated in a strong magnetic field, which has allowed observing transitions between spin- and valley-split Landau levels with various indices belonging to different graphene sheets. An unexpected increase with the temperature in the interlayer tunneling conductance owing to transitions between the Landau levels in strong magnetic fields cannot be explained by existing theories.

The mechanism of graphene growth on tantalum has been studied. It has been demonstrated that graphene growth occurs through the formation of tantalum carbide, which exhibits catalytic properties towards cracking of hydrocarbons. We found that excessive carbon, not involved in carbide formation, forms self-limiting layer of graphene on the surface. The layer of tantalum carbide and graphene provides an efficient protective layer against chemical corrosion and high temperature oxidation.

Two-dimensional materials are characterised by a number of unique physical properties which can potentially make them useful to a wide diversity of applications. In particular, the large thermal conductivity of graphene and hexagonal boron nitride (hBN) has already been acknowledged and these materials have been suggested as novel core materials for thermal management in electronics. However, it was not clear if mass produced flakes of hBN would allow one to achieve an industriallyrelevant value of thermal conductivity. Here we demonstrate that laminates of hBN exhibit thermal conductivity of up to 20W/m.K, which is significantly larger than that currently used in thermal management. Wealso show that the thermal conductivity of laminates increases with the increasing volumetric mass density, which creates a way of fine tuning its thermal properties.