Energy of the 229Th nuclear clock transition.

Owing to its low excitation energy and long radiative lifetime, the first excited isomeric state of thorium-229, 229mTh, can be optically controlled by a laser1,2 and is an ideal candidate for the creation of a nuclear optical clock3, which is expected to complement and outperform current electronic-shell-based atomic clocks4. A nuclear clock will have various applications-such as in relativistic geodesy5, dark matter research6 and the observation of potential temporal variations of fundamental constants7-but its development has so far been impeded by the imprecise knowledge of the energy of 229mTh. Here we report a direct measurement of the transition energy of this isomeric state to the ground state with an uncertainty of 0.17 electronvolts (one standard deviation) using spectroscopy of the internal conversion electrons emitted in flight during the decay of neutral 229mTh atoms. The energy of the transition between the ground state and the first excited state corresponds to a wavelength of 149.7 ± 3.1 nanometres, which is accessible by laser spectroscopy through high-harmonic generation. Our method combines nuclear and atomic physics measurements to advance precision metrology, and our findings are expected to facilitate the application of high-resolution laser spectroscopy on nuclei and to enable the development of a nuclear optical clock of unprecedented accuracy.

Quantifying Strain in Moiré Superlattice.

In twisted van der Waals (vdW) bilayers, intrinsic strain associated with the moiré superlattice and unintentionally introduced uniaxial strain may be present simultaneously. Both strains are able to lift the degeneracy of the E2g phonon modes in Raman spectra. Because of the different rotation symmetry of the two types of strain, the corresponding Raman intensity exhibits a distinct polarization dependence. We compare a 2.5° twisted MoS2 bilayer, in which the maximal intrinsic moiré strain is anticipated, and a natural MoS2 bilayer with an intentionally introduced uniaxial strain. By analyzing the frequency shift of the E2g doublet and their polarization dependence, we can not only determine the direction of unintentional uniaxial strain in the twisted bilayer but also quantify both strain components. This simple strain characterization method based on far-field Raman spectra will facilitate the studies of electronic properties of moiré superlattices under the influence of combined intrinsic and external strains.

Photonic effects in the non-equilibrium optical response of two-dimensional semiconductors.

Transient absorption spectroscopy is a powerful tool to monitor the out-of-equilibrium optical response of photoexcited semiconductors. When this method is applied to two-dimensional semiconductors deposited on different substrates, the excited state optical properties are inferred from the pump-induced changes in the transmission/reflection of the probe, i.e., ΔT/T or ΔR/R. Transient optical spectra are often interpreted as the manifestation of the intrinsic optical response of the monolayer, including effects such as the reduction of the exciton oscillator strength, electron-phonon coupling or many-body interactions like bandgap renormalization, trion or biexciton formation. Here we scrutinize the assumption that one can determine the non-equilibrium optical response of the TMD without accounting for the substrate used in the experiment. We systematically investigate the effect of the substrate on the broadband transient optical response of monolayer MoS2 (1L-MoS2) by measuring ΔT/T and ΔR/R with different excitation photon energies. Employing the boundary conditions given by the Fresnel equations, we analyze the transient transmission/reflection spectra across the main excitonic resonances of 1L-MoS2. We show that pure interference effects induced by the different substrates explain the substantial differences (i.e., intensity, peak energy and exciton linewidth) observed in the transient spectra of the same monolayer. We thus demonstrate that the substrate strongly affects the magnitude of the exciton energy shift and the change of the oscillator strength in the transient optical spectra. By highlighting the key role played by the substrate, our results set the stage for a unified interpretation of the transient response of optoelectronic devices based on a broad class of TMDs.

Challenges of modeling nanostructured materials for photocatalytic water splitting.

Understanding the water splitting mechanism in photocatalysis is a rewarding goal as it will allow producing clean fuel for a sustainable life in the future. However, identifying the photocatalytic mechanisms by modeling photoactive nanoparticles requires sophisticated computational techniques based on multiscale modeling. In this review, we will survey the strengths and drawbacks of currently available theoretical methods at different length and accuracy scales. Understanding the surface-active site through Density Functional Theory (DFT) using new, more accurate exchange-correlation functionals plays a key role for surface engineering. Larger scale dynamics of the catalyst/electrolyte interface can be treated with Molecular Dynamics albeit there is a need for more generalizations of force fields. Monte Carlo and Continuum Modeling techniques are so far not the prominent path for modeling water splitting but interest is growing due to the lower computational cost and the feasibility to compare the modeling outcome directly to experimental data. The future challenges in modeling complex nano-photocatalysts involve combining different methods in a hierarchical way so that resources are spent wisely at each length scale, as well as accounting for excited states chemistry that is important for photocatalysis, a path that will bring devices closer to the theoretical limit of photocatalytic efficiency.

Phonon renormalization in reconstructed MoS2 moiré superlattices.

In moiré crystals formed by stacking van der Waals materials, surprisingly diverse correlated electronic phases and optical properties can be realized by a subtle change in the twist angle. Here, we discover that phonon spectra are also renormalized in MoS2 twisted bilayers, adding an insight to moiré physics. Over a range of small twist angles, the phonon spectra evolve rapidly owing to ultra-strong coupling between different phonon modes and atomic reconstructions of the moiré pattern. We develop a low-energy continuum model for phonons that overcomes the outstanding challenge of calculating the properties of large moiré supercells and successfully captures the essential experimental observations. Remarkably, simple optical spectroscopy experiments can provide information on strain and lattice distortions in moiré crystals with nanometre-size supercells. The model promotes a comprehensive and unified understanding of the structural, optical and electronic properties of moiré superlattices.

Dynamically encircling an exceptional point for asymmetric mode switching.

Physical systems with loss or gain have resonant modes that decay or grow exponentially with time. Whenever two such modes coalesce both in their resonant frequency and their rate of decay or growth, an 'exceptional point' occurs, giving rise to fascinating phenomena that defy our physical intuition. Particularly intriguing behaviour is predicted to appear when an exceptional point is encircled sufficiently slowly, such as a state-flip or the accumulation of a geometric phase. The topological structure of exceptional points has been experimentally explored, but a full dynamical encircling of such a point and the associated breakdown of adiabaticity have remained out of reach of measurement. Here we demonstrate that a dynamical encircling of an exceptional point is analogous to the scattering through a two-mode waveguide with suitably designed boundaries and losses. We present experimental results from a corresponding waveguide structure that steers incoming waves around an exceptional point during the transmission process. In this way, mode transitions are induced that transform this device into a robust and asymmetric switch between different waveguide modes. This work will enable the exploration of exceptional point physics in system control and state transfer schemes at the crossroads between fundamental research and practical applications.

Local embedding of coupled cluster theory into the random phase approximation using plane waves.

We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.

Band Nesting in Two-Dimensional Crystals: An Exceptionally Sensitive Probe of Strain.

Band nesting occurs when conduction and valence bands are approximately equispaced over regions in the Brillouin zone. In two-dimensional materials, band nesting results in singularities of the joint density of states and thus in a strongly enhanced optical response at resonant frequencies. We exploit the high sensitivity of such resonances to small changes in the band structure to sensitively probe strain in semiconducting transition metal dichalcogenides (TMDs). We measure and calculate the polarization-resolved optical second harmonic generation (SHG) at the band nesting energies and present the first measurements of the energy-dependent nonlinear photoelastic effect in atomically thin TMDs (MoS2, MoSe2, WS2, and WSe2) combined with a theoretical analysis of the underlying processes. Experiment and theory are found to be in good qualitative agreement displaying a strong energy dependence of the SHG, which can be exploited to achieve exceptionally strong modulation of the SHG under strain. We attribute this sensitivity to a redistribution of the joint density of states for the optical response in the band nesting region. We predict that this exceptional strain sensitivity is a general property of all 2D materials with band nesting.

Secondary Electron Emission by Plasmon-Induced Symmetry Breaking in Highly Oriented Pyrolytic Graphite.

Two-particle spectroscopy with correlated electron pairs is used to establish the causal link between the secondary electron spectrum, the (π+σ) plasmon peak, and the unoccupied band structure of highly oriented pyrolytic graphite. The plasmon spectrum is resolved with respect to the involved interband transitions and clearly exhibits final state effects, in particular due to the energy gap between the interlayer resonances along the ΓA direction. The corresponding final state effects can also be identified in the secondary electron spectrum. Interpretation of the results is performed on the basis of density-functional theory and tight-binding calculations. Excitation of the plasmon perturbs the symmetry of the system and leads to hybridization of the interlayer resonances with atomlike σ^{*} bands along the ΓA direction. These hybrid states have a high density of states as well as sufficient mobility along the graphite c axis leading to the sharp ∼3 eV resonance in the spectrum of emitted secondary electrons reported throughout the literature.

Localized Intervalley Defect Excitons as Single-Photon Emitters in WSe_{2}.

Single-photon emitters play a key role in present and emerging quantum technologies. Several recent measurements have established monolayer WSe_{2} as a promising candidate for a reliable single-photon source. The origin and underlying microscopic processes have remained, however, largely elusive. We present a multiscale tight-binding simulation for the optical spectra of WSe_{2} under nonuniform strain and in the presence of point defects employing the Bethe-Salpeter equation. Strain locally shifts excitonic energy levels into the band gap where they overlap with localized intragap defect states. The resulting hybridization allows for efficient filling and subsequent radiative decay of the defect states. We identify intervalley defect excitonic states as the likely candidate for antibunched single-photon emission. This proposed scenario is shown to account for a large variety of experimental observations including brightness, radiative transition rates, the variation of the excitonic energy with applied magnetic and electric fields as well as the variation of the polarization of the emitted photon with the magnetic field.

Topologically Nontrivial Valley States in Bilayer Graphene Quantum Point Contacts.

We present measurements of quantized conductance in electrostatically induced quantum point contacts in bilayer graphene. The application of a perpendicular magnetic field leads to an intricate pattern of lifted and restored degeneracies with increasing field: at zero magnetic field the degeneracy of quantized one-dimensional subbands is four, because of a twofold spin and a twofold valley degeneracy. By switching on the magnetic field, the valley degeneracy is lifted. Because of the Berry curvature, states from different valleys split linearly in magnetic field. In the quantum Hall regime fourfold degenerate conductance plateaus reemerge. During the adiabatic transition to the quantum Hall regime, levels from one valley shift by two in quantum number with respect to the other valley, forming an interweaving pattern that can be reproduced by numerical calculations.

Electrostatically Confined Monolayer Graphene Quantum Dots with Orbital and Valley Splittings.

The electrostatic confinement of massless charge carriers is hampered by Klein tunneling. Circumventing this problem in graphene mainly relies on carving out nanostructures or applying electric displacement fields to open a band gap in bilayer graphene. So far, these approaches suffer from edge disorder or insufficiently controlled localization of electrons. Here we realize an alternative strategy in monolayer graphene, by combining a homogeneous magnetic field and electrostatic confinement. Using the tip of a scanning tunneling microscope, we induce a confining potential in the Landau gaps of bulk graphene without the need for physical edges. Gating the localized states toward the Fermi energy leads to regular charging sequences with more than 40 Coulomb peaks exhibiting typical addition energies of 7-20 meV. Orbital splittings of 4-10 meV and a valley splitting of about 3 meV for the first orbital state can be deduced. These experimental observations are quantitatively reproduced by tight binding calculations, which include the interactions of the graphene with the aligned hexagonal boron nitride substrate. The demonstrated confinement approach appears suitable to create quantum dots with well-defined wave function properties beyond the reach of traditional techniques.

Two-Color Coherent Control of Femtosecond Above-Threshold Photoemission from a Tungsten Nanotip.

We demonstrate coherent control of multiphoton and above-threshold photoemission from a single solid-state nanoemitter driven by a fundamental and a weak second harmonic laser pulse. Depending on the relative phase of the two pulses, electron emission is modulated with a contrast of the oscillating current signal of up to 94%. Electron spectra reveal that all observed photon orders are affected simultaneously and similarly. We confirm that photoemission takes place within 10 fs. Accompanying simulations indicate that the current modulation with its large contrast results from two interfering quantum pathways leading to electron emission.

Embedded correlated wavefunction schemes: theory and applications.

Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFT's ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.

Implementation of density functional embedding theory within the projector-augmented-wave method and applications to semiconductor defect states.

We report a new implementation of the density functional embedding theory (DFET) in the VASP code, using the projector-augmented-wave (PAW) formalism. Newly developed algorithms allow us to efficiently perform optimized effective potential optimizations within PAW. The new algorithm generates robust and physically correct embedding potentials, as we verified using several test systems including a covalently bound molecule, a metal surface, and bulk semiconductors. We show that with the resulting embedding potential, embedded cluster models can reproduce the electronic structure of point defects in bulk semiconductors, thereby demonstrating the validity of DFET in semiconductors for the first time. Compared to our previous version, the new implementation of DFET within VASP affords use of all features of VASP (e.g., a systematic PAW library, a wide selection of functionals, a more flexible choice of U correction formalisms, and faster computational speed) with DFET. Furthermore, our results are fairly robust with respect to both plane-wave and Gaussian type orbital basis sets in the embedded cluster calculations. This suggests that the density functional embedding method is potentially an accurate and efficient way to study properties of isolated defects in semiconductors.

Photovoltaic effect in an electrically tunable van der Waals heterojunction.

Semiconductor heterostructures form the cornerstone of many electronic and optoelectronic devices and are traditionally fabricated using epitaxial growth techniques. More recently, heterostructures have also been obtained by vertical stacking of two-dimensional crystals, such as graphene and related two-dimensional materials. These layered designer materials are held together by van der Waals forces and contain atomically sharp interfaces. Here, we report on a type-II van der Waals heterojunction made of molybdenum disulfide and tungsten diselenide monolayers. The junction is electrically tunable, and under appropriate gate bias an atomically thin diode is realized. Upon optical illumination, charge transfer occurs across the planar interface and the device exhibits a photovoltaic effect. Advances in large-scale production of two-dimensional crystals could thus lead to a new photovoltaic solar technology.

Time-dependent potential-functional embedding theory.

We introduce a time-dependent potential-functional embedding theory (TD-PFET), in which atoms are grouped into subsystems. In TD-PFET, subsystems can be propagated by different suitable time-dependent quantum mechanical methods and their interactions can be treated in a seamless, first-principles manner. TD-PFET is formulated based on the time-dependent quantum mechanics variational principle. The action of the total quantum system is written as a functional of the time-dependent embedding potential, i.e., a potential-functional formulation. By exploiting the Runge-Gross theorem, we prove the uniqueness of the time-dependent embedding potential under the constraint that all subsystems share a common embedding potential. We derive the integral equation that such an embedding potential needs to satisfy. As proof-of-principle, we demonstrate TD-PFET for a Na4 cluster, in which each Na atom is treated as one subsystem and propagated by time-dependent Kohn-Sham density functional theory (TDDFT) using the adiabatic local density approximation (ALDA). Our results agree well with a direct TDDFT calculation on the whole Na4 cluster using ALDA. We envision that TD-PFET will ultimately be useful for studying ultrafast quantum dynamics in condensed matter, where key regions are solved by highly accurate time-dependent quantum mechanics methods, and unimportant regions are solved by faster, less accurate methods.

Analysis of and remedies for unphysical ground states of the multireference averaged coupled-pair functional.

A Multireference Configuration Interaction (MRCI) wavefunction includes both static and dynamic electron correlation. MRCI's well-known flaw, a lack of size extensivity, can be ameliorated with the Multireference Averaged Coupled-Pair Functional (MRACPF). However, the original MRACPF is frequently unstable, sometimes producing unphysical results. The more Multireference Averaged Quadratic Coupled-Cluster and MRACPF2 functionals also occasionally exhibit unphysical behavior. We find that these instabilities are avoided crossings with unphysical solutions to the MRACPF equations. We present two approaches to avoid the undesirable unphysical solutions.

Size-extensivity-corrected multireference configuration interaction schemes to accurately predict bond dissociation energies of oxygenated hydrocarbons.

Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs.