Search

We present an application of periodic coupled-cluster theory to the calculation of CO adsorption energies on the Pt(111) surface for different adsorption sites. The calculations employ a range of recently developed theoretical and computational methods. In particular, we use a recently introduced coupled-cluster ansatz, denoted as CCSD(cT), to compute correlation energies of the metallic Pt surface with and without adsorbed CO molecules. The convergence of Hartree-Fock adsorption energy contributions with respect to randomly shifted k -meshes is discussed. Recently introduced basis set incompleteness error corrections make it possible to achieve well-converged correlation energy contributions to the adsorption energies. We show that CCSD(cT) theory predicts the correct order of adsorption energies for the considered adsorption sites. Furthermore, we find that binding of the CO molecule to the top and fcc site is dominated by Hartree-Fock and correlation energy contributions, respectively.

We investigate the convergence of coupled-cluster (CC) correlation energies and related quantities with respect to the employed basis set size for the uniform electron gas (UEG) to gain a better understanding of the basis set incompleteness error (BSIE). To this end, coupled-cluster doubles (CCD) theory is applied to the three-dimensional UEG for a range of densities, basis set sizes, and electron numbers. We present a detailed analysis of individual diagrammatically decomposed contributions to the amplitudes at the level of CCD theory. In particular, we show that only two terms from the amplitude equations contribute to the asymptotic large-momentum behavior of the transition structure factor, corresponding to the cusp region at short interelectronic distances. However, due to the coupling present in the amplitude equations, all decomposed correlation energy contributions show the same asymptotic convergence behavior to the complete basis set limit. These findings provide an additional rationale for the success of a recently proposed correction to the BSIE of CC theory. Lastly, we examine the BSIE in the CCD plus perturbative triples [CCD(T)] method, as well as in the newly proposed CCD plus complete perturbative triples [CCD(cT)] method.

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Experimental setups for charge transport measurements are typically not compatible with the ultrahigh vacuum conditions for chemical doping, limiting the charge carrier density that can be investigated by transport methods. Field-effect methods, including dielectric gating and ionic liquid gating, achieve too low a carrier density to induce electronic phase transitions. To bridge this gap, we developed an integrated flip-chip method to dope graphene by alkali vapor in the diffusive regime, suitable for charge transport measurements at ultrahigh charge carrier density. We introduce a cesium droplet into a sealed cavity filled with inert gas to dope a monolayer graphene sample by the process of cesium atom diffusion, adsorption, and ionization at the graphene surface, with doping beyond an electron density of 4.7 × 10 14 cm -2 monitored by operando Hall measurement. The sealed assembly is stable against oxidation, enabling measurement of charge transport versus temperature and magnetic field. Cyclotron mass inversion is observed via the Hall effect, indicative of the change in Fermi surface geometry associated with the Liftshitz transition at the hyperbolic M point of monolayer graphene. The transparent quartz substrate also functions as an optical window, enabling nonresonant Raman scattering. Our findings show that chemical doping, hitherto restricted to ultrahigh vacuum, can be applied in a diffusive regime at ambient pressure in an inert gas environment and thus enable charge transport studies in standard cryogenic environments.

We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions: the treatment of the singularity at the Brillouin-zone center and discretization errors, which can cause severe convergence problems in anisotropic cells, necessary for the calculation of low-dimensional systems. We apply our strategy to the Hartree-Fock and coupled cluster (CC) theories and discuss the consequences of different sampling strategies on different theories. We show that sampling the Coulomb potential via the widely used probe-charge Ewald method is unsuitable for CC calculations in anisotropic cells. To demonstrate the applicability of our developed approach, we study two representative, low-dimensional use cases: the infinite carbon chain, for which we report the first periodic CCSD(T) potential energy surface, and a surface slab of lithium hydride, for which we demonstrate the impact of different sampling strategies for calculating surface energies. We find that our Coulomb sampling strategy serves as a vital solution, addressing the critical need for improved accuracy in plane-wave based CC calculations for low-dimensional systems., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.

Coupled-cluster theories can be used to compute ab initio electronic correlation energies of real materials with systematically improvable accuracy. However, the widely used coupled cluster singles and doubles plus perturbative triples [CCSD(T)] method is only applicable to insulating materials. For zero-gap materials the truncation of the underlying many-body perturbation expansion leads to an infrared catastrophe. Here, we present a novel perturbative triples formalism denoted as (cT) that yields convergent correlation energies in metallic systems. Furthermore, the computed correlation energies for the three-dimensional uniform electron gas at metallic densities are in good agreement with quantum Monte Carlo results. At the same time the newly proposed method retains all desirable properties of CCSD(T) such as its accuracy for insulating systems as well as its low computational cost compared to a full inclusion of the triples. This paves the way for ab initio calculations of real metals with chemical accuracy.

We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and, thus, reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes compared to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles, and perturbative triple particle-hole excitation operators., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Boron nitride nanotubes (BNNTs) are an emerging class of molecular container offering new functionalities and possibilities for studying molecules at the nanoscale. Herein, BNNTs are demonstrated as highly effective nanocontainers for polyoxometalate (POM) molecules. The encapsulation of POMs within BNNTs occurs spontaneously at room temperature from an aqueous solution, leading to the self-assembly of a POM@BNNT host-guest system. Analysis of the interactions between the host-nanotube and guest-molecule indicate that Lewis acid-base interactions between W═O groups of the POM (base) and B-atoms of the BNNT lattice (acid) likely play a major role in driving POM encapsulation, with photoactivated electron transfer from BNNTs to POMs in solution also contributing to the process. The transparent nature of the BNNT nanocontainer allows extensive investigation of the guest-molecules by photoluminescence, Raman, UV-vis absorption, and EPR spectroscopies. These studies revealed considerable energy and electron transfer processes between BNNTs and POMs, likely mediated via defect energy states of the BNNTs and resulting in the quenching of BNNT photoluminescence at room temperature, the emergence of new photoluminescence emissions at cryogenic temperatures (<100 K), a photochromic response, and paramagnetic signals from guest-POMs. These phenomena offer a fresh perspective on host-guest interactions at the nanoscale and open pathways for harvesting the functional properties of these hybrid systems.

A magnetic field modifies optical properties and provides valley splitting in a molybdenum disulfide (MoS 2 ) monolayer. Here we demonstrate a scalable approach to the epitaxial synthesis of MoS 2 monolayer on a magnetic graphene/Co system. Using spin- and angle-resolved photoemission spectroscopy we observe a magnetic proximity effect that causes a 20 meV spin-splitting at the Γ̅ point and canting of spins at the K̅ point in the valence band toward the in-plane direction of cobalt magnetization. Our density functional theory calculations reveal that the in-plane spin component at K̅ is localized on Co atoms in the valence band, while in the conduction band it is localized on the MoS 2 layer. The calculations also predict a 16 meV spin-splitting at the Γ̅ point and 8 meV K̅- K ' ¯ valley asymmetry for an out-of-plane magnetization. These findings suggest control over optical transitions in MoS 2 via Co magnetization. Our estimations show that the magnetic proximity effect is equivalent to the action of the magnetic field as large as 100 T.

We show that Cs intercalated bilayer graphene acts as a substrate for the growth of a strained Cs film hosting quantum well states with high electronic quality. The Cs film grows in an fcc phase with a substantially reduced lattice constant of 4.9 Å corresponding to a compressive strain of 11% compared to bulk Cs. We investigate its electronic structure using angle-resolved photoemission spectroscopy and show the coexistence of massless Dirac and massive Schrödinger charge carriers in two dimensions. Analysis of the electronic self-energy of the massive charge carriers reveals the crystallographic direction in which a two-dimensional Fermi gas is realized. Our work introduces the growth of strained metal quantum wells on intercalated Dirac matter.

A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride is reported employing an embedding approach for many-electron correlation methods. To this end, a plane-wave based implementation to obtain intrinsic atomic orbitals and Wannier functions for arbitrary localization potentials is presented. In our embedding scheme, the localized occupied orbitals allow for a separate treatment of short-range and long-range correlation contributions to the adsorption energy by a fragmentation of the simulation cell. In combination with unoccupied natural orbitals, the coupled cluster ansatz with single, double, and perturbative triple particle-hole excitation operators is used to capture the correlation in local fragments centered around the adsorption process. For the long-range correlation, a seamless embedding into the random phase approximation yields rapidly convergent adsorption energies with respect to the local fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors., (© 2021. The Author(s).)

We present a basis set correction scheme for the coupled-cluster singles and doubles (CCSD) method. The scheme is based on employing frozen natural orbitals (FNOs) and diagrammatically decomposed contributions to the electronic correlation energy, which dominate the basis set incompleteness error (BSIE). As recently discussed in the work of Irmler et al. [Phys. Rev. Lett. 123, 156401 (2019)], the BSIE of the CCSD correlation energy is dominated by the second-order Møller-Plesset (MP2) perturbation energy and the particle-particle ladder term. Here, we derive a simple approximation to the BSIE of the particle-particle ladder term that effectively corresponds to a rescaled pair-specific MP2 BSIE, where the scaling factor depends on the spatially averaged correlation hole depth of the coupled-cluster and first-order pair wavefunctions. The evaluation of the derived expressions is simple to implement in any existing code. We demonstrate the effectiveness of the method for the uniform electron gas. Furthermore, we apply the method to coupled-cluster theory calculations of atoms and molecules using FNOs. Employing the proposed correction and an increasing number of FNOs per occupied orbital, we demonstrate for a test set that rapidly convergent closed and open-shell reaction energies, atomization energies, electron affinities, and ionization potentials can be obtained. Moreover, we show that a similarly excellent trade-off between required virtual orbital basis set size and remaining BSIEs can be achieved for the perturbative triples contribution to the CCSD(T) energy employing FNOs and the (T*) approximation.

Lateral heterojunctions of atomically precise graphene nanoribbons (GNRs) hold promise for applications in nanotechnology, yet their charge transport and most of the spectroscopic properties have not been investigated. Here, we synthesize a monolayer of multiple aligned heterojunctions consisting of quasi-metallic and wide-bandgap GNRs, and report characterization by scanning tunneling microscopy, angle-resolved photoemission, Raman spectroscopy, and charge transport. Comprehensive transport measurements as a function of bias and gate voltages, channel length, and temperature reveal that charge transport is dictated by tunneling through the potential barriers formed by wide-bandgap GNR segments. The current-voltage characteristics are in agreement with calculations of tunneling conductance through asymmetric barriers. We fabricate a GNR heterojunctions based sensor and demonstrate greatly improved sensitivity to adsorbates compared to graphene based sensors. This is achieved via modulation of the GNR heterojunction tunneling barriers by adsorbates.

We present an implementation of the equation of motion coupled-cluster singles and doubles (EOM-CCSD) theory using periodic boundary conditions and a plane wave basis set. Our implementation of EOM-CCSD theory is applied to study F-centers in alkaline earth oxides employing a periodic supercell approach. The convergence of the calculated electronic excitation energies for neutral color centers in MgO, CaO, and SrO crystals with respect to the orbital basis set and system size is explored. We discuss extrapolation techniques that approximate excitation energies in the complete basis set limit and reduce finite size errors. Our findings demonstrate that EOM-CCSD theory can predict optical absorption energies of F-centers in good agreement with experiment. Furthermore, we discuss calculated emission energies corresponding to the decay from triplet to singlet states responsible for the photoluminescence properties. Our findings are compared to experimental and theoretical results available in the literature.

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.

Cluster superlattice membranes consist of a two-dimensional hexagonal lattice of similar-sized nanoclusters sandwiched between single-crystal graphene and an amorphous carbon matrix. The fabrication process involves three main steps, the templated self-organization of a metal cluster superlattice on epitaxial graphene on Ir(111), conformal embedding in an amorphous carbon matrix, and subsequent lift-off from the Ir(111) substrate. The mechanical stability provided by the carbon-graphene matrix makes the membrane stable as a free-standing material and enables transfer to other substrates. The fabrication procedure can be applied to a wide variety of cluster materials and cluster sizes from the single-atom limit to clusters of a few hundred atoms, as well as other two-dimensional layer/host matrix combinations. The versatility of the membrane composition, its mechanical stability, and the simplicity of the transfer procedure make cluster superlattice membranes a promising material in catalysis, magnetism, energy conversion, and optoelectronics.

We present laser-induced photothermal synthesis of atomically precise graphene nanoribbons (GNRs). The kinetics of photothermal bottom-up GNR growth are unravelled by in situ Raman spectroscopy carried out in ultrahigh vacuum. We photothermally drive the reaction steps by short periods of laser irradiation and subsequently analyze the Raman spectra of the reactants in the irradiated area. Growth kinetics of chevron GNRs (CGNRs) and seven atoms wide armchair GNRs (7-AGNRs) is investigated. The reaction rate constants for polymerization, cyclodehydrogenation, and interribbon fusion are experimentally determined. We find that the limiting rate constants for CGNR growth are several hundred times smaller than for 7-AGNR growth and that interribbon fusion is an important elementary reaction occurring during 7-AGNR growth. Our work highlights that photothermal synthesis and in situ Raman spectroscopy are a powerful tandem for the investigation of on-surface reactions.

A flat energy dispersion of electrons at the Fermi level of a material leads to instabilities in the electronic system and can drive phase transitions. Here we show that the flat band in graphene can be achieved by sandwiching a graphene monolayer by two cesium (Cs) layers. We investigate the flat band by a combination of angle-resolved photoemission spectroscopy experiment and the calculations. Our work highlights that charge transfer, zone folding of graphene bands, and the covalent bonding between C and Cs atoms are the origin of the flat energy band formation. Analysis of the Stoner criterion for the flat band suggests the presence of a ferromagnetic instability. The presented approach is an alternative route for obtaining flat band materials to twisting bilayer graphene which yields thermodynamically stable flat band materials in large areas.

We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the contribution of the ladder diagrams to the electronic correlation energy can be approximated in an effective manner using second-order perturbation theory. We employ the latter approximation to reduce the computational cost of coupled cluster theory calculations for insulators and semiconductors by 2 orders of magnitude without compromising accuracy.

For quasi-freestanding 2H-TaS 2 in monolayer thickness grown by in situ molecular beam epitaxy on graphene on Ir(111), we find unambiguous evidence for a charge density wave close to a 3 × 3 periodicity. Using scanning tunneling spectroscopy, we determine the magnitude of the partial charge density wave gap. Angle-resolved photoemission spectroscopy, complemented by scanning tunneling spectroscopy for the unoccupied states, makes a tight-binding fit for the band structure of the TaS 2 monolayer possible. As hybridization with substrate bands is absent, the fit yields a precise value for the doping of the TaS 2 layer. Additional Li doping shifts the charge density wave to a 2 × 2 periodicity. Unexpectedly, the bilayer of TaS 2 also displays a disordered 2 × 2 charge density wave. Calculations of the phonon dispersions based on a combination of density-functional theory, density-functional perturbation theory, and many-body perturbation theory enable us to provide phase diagrams for the TaS 2 charge density wave as functions of doping, hybridization, and interlayer potentials, and offer insight into how they affect lattice dynamics and stability. Our theoretical considerations are consistent with the experimental work presented and shed light on previous experimental and theoretical investigations of related systems.

We investigate the basis-set convergence of electronic correlation energies calculated using coupled cluster theory and a recently proposed finite basis-set correction technique. The correction is applied to atomic and molecular systems and is based on a diagrammatically decomposed coupled cluster singles and doubles (CCSD) correlation energy. Only the second-order energy and the particle-particle ladder term are corrected for their basis-set incompleteness error. We present absolute correlation energies and results for a large benchmark set. Our findings indicate that basis set reductions by two cardinal numbers are possible for atomization energies, ionization potentials, and electron affinities without compromising accuracy when compared to conventional CCSD calculations. In the case of reaction energies, we find that reductions by one cardinal number are possible compared to conventional CCSD calculations. The employed technique can readily be applied to other many-electron theories without the need for three- or four-electron integrals.

The random phase approximation (RPA) systematically overestimates the magnitude of the correlation energy and generally underestimates cohesive energies. This originates in part from the complete lack of exchange terms that would otherwise cancel Pauli exclusion principle violating (EPV) contributions. The uncanceled EPV contributions also manifest themselves in form of an unphysical negative pair density of spin-parallel electrons close to electron-electron coalescence. We follow considerations of many-body perturbation theory to propose an exchange correction that corrects the largest set of EPV contributions, while having the lowest possible computational complexity. The proposed method exchanges adjacent particle/hole pairs in the RPA diagrams, considerably improving the pair density of spin-parallel electrons close to coalescence in the uniform electron gas (UEG). The accuracy of the correlation energy is comparable to other variants of second-order screened exchange (SOSEX) corrections although it is slightly more accurate for the spin-polarized UEG. Its computational complexity scales as [Formula: see text] or [Formula: see text] in orbital space or real space, respectively. Its memory requirement scales as [Formula: see text].

Wet carbon interfaces are ubiquitous in the natural world and exhibit anomalous properties, which could be exploited by emerging technologies. However, progress is limited by lack of understanding at the molecular level. Remarkably, even for the most fundamental system (a single water molecule interacting with graphene), there is no consensus on the nature of the interaction. We tackle this by performing an extensive set of complementary state-of-the-art computer simulations on some of the world's largest supercomputers. From this effort a consensus on the water-graphene interaction strength has been obtained. Our results have significant impact for the physical understanding, as they indicate that the interaction is weaker than predicted previously. They also pave the way for more accurate and reliable studies of liquid water at carbon interfaces.

We report on a many-electron wavefunction theory study for the reaction energetics of hydrogen dissociation on the Si(100) surface. We demonstrate that quantum chemical wavefunction based methods using periodic boundary conditions can predict chemically accurate results for the activation barrier and the chemisorption energy in agreement with experimental findings. These highly accurate results for the reaction energetics enable a deeper understanding of the underlying physical mechanism and make it possible to benchmark widely used density functional theory methods.

By enhancing the photoluminescence from aligned seven-atom wide armchair-edge graphene nanoribbons using plasmonic nanoantennas, we are able to observe blinking of the emission. The on- and off-times of the blinking follow power law statistics. In time-resolved spectra, we observe spectral diffusion. These findings together are a strong indication of the emission originating from a single quantum emitter. The room temperature photoluminescence displays a narrow spectral width of less than 50 meV, which is significantly smaller than the previously observed ensemble line width of 0.8 eV. From spectral time traces, we identify three optical transitions, which are energetically situated below the lowest bulk excitonic state E 11 of the nanoribbons. We attribute the emission to transitions involving Tamm states localized at the end of the nanoribbon. The photoluminescence from a single ribbon is strongly enhanced when its end is in the antenna hot spot resulting in the observed single molecule characteristics of the emission. Our findings illustrate the essential role of the end termination of graphene nanoribbons in light emission and allow us to construct a model for photoluminescence from nanoribbons.

The opening of a band gap in graphene nanoribbons induces novel optical and electronic properties, strongly enhancing their application potential in nanoscale devices. Knowledge of the optical excitations and associated relaxation dynamics are essential for developing and optimizing device designs and functionality. Here we report on the optical excitations and associated relaxation dynamics in surface aligned 7-atom wide armchair graphene nanoribbons as seen by time-resolved spontaneous Stokes and anti-Stokes Raman scattering spectroscopy. On the anti-Stokes side we observe an optically induced increase of the scattering intensity of the Raman active optical phonons which we assign to changes in the optical phonon populations. The optical phonon population decays with a lifetime of ∼2 ps, indicating an efficient optical-acoustic phonon cooling mechanism. On the Stokes side we observe a substantial decrease of the phonon peak intensities which we relate to the dynamics of the optically induced exciton population. The exciton population shows a multi-exponential relaxation on the hundreds of ps time scale and is independent of the excitation intensity, indicating that exciton-exciton annihilation processes are not important and the exsistence of dark and trapped exciton states. Our results shed light on the optically induced phonon and exciton dynamics in surface aligned armchair graphene nanoribbons and demonstrate that time-resolved spontaneous Raman scattering spectroscopy is a powerful method for exploring quasi-particle dynamics in low dimensional materials.

We employ ultra-high vacuum (UHV) Raman spectroscopy in tandem with angle-resolved photoemission (ARPES) to investigate the doping-dependent Raman spectrum of epitaxial graphene on Ir(111). The evolution of Raman spectra from pristine to heavily Cs doped graphene up to a carrier concentration of 4.4 × 10 14 cm -2 is investigated. At this doping, graphene is at the onset of the Lifshitz transition and renormalization effects reduce the electronic bandwidth. The optical transition at the saddle point in the Brillouin zone then becomes experimentally accessible by ultraviolet (UV) light excitation, which achieves resonance Raman conditions in close vicinity to the van Hove singularity in the joint density of states. The position of the Raman G band of fully doped graphene/Ir(111) shifts down by ∼60 cm -1 . The G band asymmetry of Cs doped epitaxial graphene assumes an unusual strong Fano asymmetry opposite to that of the G band of doped graphene on insulators. Our calculations can fully explain these observations by substrate dependent quantum interference effects in the scattering pathways for vibrational and electronic Raman scattering.

We investigate the electronic and vibrational properties of bottom-up synthesized aligned armchair graphene nanoribbons of N = 7 carbon atoms width periodically doped by substitutional boron atoms (B-7AGNRs). Using angle-resolved photoemission spectroscopy and density functional theory calculations, we find that the dopant-derived valence and conduction band states are notably hybridized with electronic states of Au substrate and spread in energy. The interaction with the substrate leaves the bands with pure carbon character rather unperturbed. This results in an identical effective mass of ≈0.2 m 0 for the next-highest valence band compared with pristine 7AGNRs. We probe the phonons of B-7AGNRs by ultrahigh-vacuum (UHV) Raman spectroscopy and reveal the existence of characteristic splitting and red shifts in Raman modes due to the presence of substitutional boron atoms. Comparing the Raman spectra for three visible lasers (red, green, and blue), we find that interaction with gold suppresses the Raman signal from B-7AGNRs and the energy of the green laser (2.33 eV) is closer to the resonant E 22 transition. The hybridized electronic structure of the B-7AGNR-Au interface is expected to improve electrical characteristics of contacts between graphene nanoribbon and Au. The Raman fingerprint allows the easy identification of B-7AGNRs, which is particularly useful for device fabrication.

We report on the experimental demonstration and electrical characterization of N = 7 armchair graphene nanoribbon (7-AGNR) field effect transistors. The back-gated transistors are fabricated from atomically precise and highly aligned 7-AGNRs, synthesized with a bottom-up approach. The large area transfer process holds the promise of scalable device fabrication with atomically precise nanoribbons. The channels of the FETs are approximately 30 times longer than the average nanoribbon length of 30 nm to 40 nm. The density of the GNRs is high, so that transport can be assumed well-above the percolation threshold. The long channel transistors exhibit a maximum I ON / I OFF current ratio of 87.5.

We report on the synthesis and spectroscopic characterization of alkali metal intercalated ZrSe 2 single crystals. ZrSe 2 is produced by chemical vapour transport and then Li is intercalated. Intercalation is performed from the liquid phase (via butyllithium) and from the vapour phase. Raman spectroscopy of intercalated ZrSe 2 reveals phonon energy shifts of the Raman active A 1g and E g phonon modes, the disappearance of two-phonon modes and new low wavenumber Raman modes. Angle-resolved photoemission spectroscopy is used to perform a mapping of the Fermi surface revealing an electron concentration of 4.7 × 10 14 cm -2 . We also perform vapour phase intercalation of K and Cs into ZrSe 2 and observe similar changes in the Raman modes as for the Li case.

In recent years, solid-state spin systems have emerged as promising candidates for quantum information processing. Prominent examples are the nitrogen-vacancy (NV) center in diamond, phosphorus dopants in silicon (Si:P), rare-earth ions in solids, and V Si -centers in silicon-carbide. The Si:P system has demonstrated that its nuclear spins can yield exceedingly long spin coherence times by eliminating the electron spin of the dopant. For NV centers, however, a proper charge state for storage of nuclear spin qubit coherence has not been identified yet. Here, we identify and characterize the positively charged NV center as an electron-spin-less and optically inactive state by utilizing the nuclear spin qubit as a probe. We control the electronic charge and spin utilizing nanometer scale gate electrodes. We achieve a lengthening of the nuclear spin coherence times by a factor of 4. Surprisingly, the new charge state allows switching of the optical response of single nodes facilitating full individual addressability.

Molecular adsorption on surfaces plays an important part in catalysis, corrosion, desalination, and various other processes that are relevant to industry and in nature. As a complement to experiments, accurate adsorption energies can be obtained using various sophisticated electronic structure methods that can now be applied to periodic systems. The adsorption energy of water on boron nitride substrates, going from zero to 2-dimensional periodicity, is particularly interesting as it calls for an accurate treatment of polarizable electrostatics and dispersion interactions, as well as posing a practical challenge to experiments and electronic structure methods. Here, we present reference adsorption energies, static polarizabilities, and dynamic polarizabilities, for water on BN substrates of varying size and dimension. Adsorption energies are computed with coupled cluster theory, fixed-node quantum Monte Carlo (FNQMC), the random phase approximation, and second order Møller-Plesset theory. These wavefunction based correlated methods are found to agree in molecular as well as periodic systems. The best estimate of the water/h-BN adsorption energy is -107±7 meV from FNQMC. In addition, the water adsorption energy on the BN substrates could be expected to grow monotonically with the size of the substrate due to increased dispersion interactions, but interestingly, this is not the case here. This peculiar finding is explained using the static polarizabilities and molecular dispersion coefficients of the systems, as computed from time-dependent density functional theory (DFT). Dynamic as well as static polarizabilities are found to be highly anisotropic in these systems. In addition, the many-body dispersion method in DFT emerges as a particularly useful estimation of finite size effects for other expensive, many-body wavefunction based methods.

We demonstrate the alignment-preserving transfer of parallel graphene nanoribbons (GNRs) onto insulating substrates. The photophysics of such samples is characterized by polarized Raman and photoluminescence (PL) spectroscopies. The Raman scattered light and the PL are polarized along the GNR axis. The Raman cross section as a function of excitation energy has distinct excitonic peaks associated with transitions between the one-dimensional parabolic subbands. We find that the PL of GNRs is intrinsically low but can be strongly enhanced by blue laser irradiation in ambient conditions or hydrogenation in ultrahigh vacuum. These functionalization routes cause the formation of sp 3 defects in GNRs. We demonstrate the laser writing of luminescent patterns in GNR films for maskless lithography by the controlled generation of defects. Our findings set the stage for further exploration of the optical properties of GNRs on insulating substrates and in device geometries.

We present a comprehensive benchmark study of the adsorption energy of a single water molecule on the (001) LiH surface using periodic coupled cluster and quantum Monte Carlo theories. We benchmark and compare different implementations of quantum chemical wave function based theories in order to verify the reliability of the predicted adsorption energies and the employed approximations. Furthermore we compare the predicted adsorption energies to those obtained employing widely used van der Waals density-functionals. Our findings show that quantum chemical approaches are becoming a robust and reliable tool for condensed phase electronic structure calculations, providing an additional tool that can also help in potentially improving currently available van der Waals density-functionals.

We study a tensor hypercontraction decomposition of the Coulomb integrals of periodic systems where the integrals are factorized into a contraction of six matrices of which only two are distinct. We find that the Coulomb integrals can be well approximated in this form already with small matrices compared to the number of real space grid points. The cost of computing the matrices scales as O(N 4 ) using a regularized form of the alternating least squares algorithm. The studied factorization of the Coulomb integrals can be exploited to reduce the scaling of the computational cost of expensive tensor contractions appearing in the amplitude equations of coupled cluster methods with respect to system size. We apply the developed methodologies to calculate the adsorption energy of a single water molecule on a hexagonal boron nitride monolayer in a plane wave basis set and periodic boundary conditions.

The explicitly correlated approach is one of the most important breakthroughs in ab initio electronic structure theory, providing arguably the most compact, accurate, and efficient ansatz for describing the correlated motion of electrons. Since Hylleraas first used an explicitly correlated wave function for the He atom in 1929, numerous attempts have been made to tackle the significant challenges involved in constructing practical explicitly correlated methods that are applicable to larger systems. These include identifying suitable mathematical forms of a correlated wave function and an efficient evaluation of many-electron integrals. R12 theory, which employs the resolution of the identity approximation, emerged in 1985, followed by the introduction of novel correlation factors and wave function ansätze, leading to the establishment of F12 theory in the 2000s. Rapid progress in recent years has significantly extended the application range of explicitly correlated theory, offering the potential of an accurate wave-function treatment of complex systems such as photosystems and semiconductors. This perspective surveys explicitly correlated electronic structure theory, with an emphasis on recent stochastic and deterministic approaches that hold significant promise for applications to large and complex systems including solids.

Electronically driven nematic order is often considered as an essential ingredient of high-temperature superconductivity. Its elusive nature in iron-based superconductors resulted in a controversy not only as regards its origin but also as to the degree of its influence on the electronic structure even in the simplest representative material FeSe. Here we utilized angle-resolved photoemission spectroscopy and density functional theory calculations to study the influence of the nematic order on the electronic structure of FeSe and determine its exact energy and momentum scales. Our results strongly suggest that the nematicity in FeSe is electronically driven, we resolve the recent controversy and provide the necessary quantitative experimental basis for a successful theory of superconductivity in iron-based materials which takes into account both, spin-orbit interaction and electronic nematicity.

We present a method to correct for finite size errors in coupled cluster theory calculations of solids. The outlined technique shares similarities with electronic structure factor interpolation methods used in quantum Monte Carlo calculations. However, our approach does not require the calculation of density matrices. Furthermore we show that the proposed finite size corrections achieve chemical accuracy in the convergence of second-order Møller-Plesset perturbation and coupled cluster singles and doubles correlation energies per atom for insulating solids with two atomic unit cells using 2 × 2 × 2 and 3 × 3 × 3 k-point meshes only.

We present a simple, robust, and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cutoff radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artifacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid, and water adsorption on a LiH surface, at the level of second-order Møller-Plesset perturbation theory.

The full exploration of the potential, which graphene offers to nanoelectronics requires its integration into semiconductor technology. So far the real-world applications are limited by the ability to concomitantly achieve large single-crystalline domains on dielectrics and semiconductors and to tailor the interfaces between them. Here we show a new direct bottom-up method for the fabrication of high-quality atomically precise interfaces between 2D materials, like graphene and hexagonal boron nitride (hBN), and classical semiconductor via Ge intercalation. Using angle-resolved photoemission spectroscopy and complementary DFT modelling we observed for the first time that epitaxially grown graphene with the Ge monolayer underneath demonstrates Dirac Fermions unaffected by the substrate as well as an unperturbed electronic band structure of hBN. This approach provides the intrinsic relativistic 2D electron gas towards integration in semiconductor technology. Hence, these new interfaces are a promising path for the integration of graphene and hBN into state-of-the-art semiconductor technology.

We employ Hartree-Fock, second-order Møller-Plesset perturbation, coupled cluster singles and doubles (CCSD) as well as CCSD plus perturbative triples (CCSD(T)) theory to study the pressure induced transition from the rocksalt to the cesium chloride crystal structure in LiH. We show that the calculated transition pressure converges rapidly in this series of increasingly accurate many-electron wave function based theories. Using CCSD(T) theory, we predict a transition pressure for the structural phase transition in the LiH crystal of 340 GPa. Furthermore, we investigate the potential energy surface for this transition in the parameter space of the Buerger path.

We introduce a novel and efficient explicitly correlated implementation of second-order perturbation theory for solids. The required three-electron integrals are computed directly using a plane wave basis set. We parametrize the employed correlation factors using results previously obtained for a finite uniform electron gas simulation cell. We demonstrate for a range of solids that basis set converged correlation energies, equilibrium volumes, and bulk moduli can be obtained efficiently in this theory using a few ten orbitals per atom. To stretch the capabilities of this novel method we compute the Schottky defect formation energy in MgO, studying systems with 54 atoms in the supercell. We verify the accuracy of the calculated formation energies using the more accurate coupled cluster singles and doubles theory. Furthermore, we discuss other potential applications for the derived and implemented expressions such as an occupied orbital only correlation energy functional.

With the discovery and first characterization of graphene, its potential for spintronic applications was recognized immediately. Since then, an active field of research has developed trying to overcome the practical hurdles. One of the most severe challenges is to find appropriate interfaces between graphene and ferromagnetic layers, which are granting efficient injection of spin-polarized electrons. Here, we show that graphene grown under appropriate conditions on Co(0001) demonstrates perfect structural properties and simultaneously exhibits highly spin-polarized charge carriers. The latter was conclusively proven by observation of a single-spin Dirac cone near the Fermi level. This was accomplished experimentally using spin- and angle-resolved photoelectron spectroscopy, and theoretically with density functional calculations. Our results demonstrate that the graphene/Co(0001) system represents an interesting candidate for applications in devices using the spin degree of freedom.

Oxygen reduction reaction (ORR) plays a key role in lithium-air batteries (LABs) that attract great attention thanks to their high theoretical specific energy several times exceeding that of lithium-ion batteries. Because of their high surface area, high electric conductivity, and low specific weight, various carbons are often materials of choice for applications as the LAB cathode. Unfortunately, the possibility of practical application of such batteries is still under question as the sustainable operation of LABs with carbon cathodes is not demonstrated yet and the cyclability is quite poor, which is usually associated with oxygen reduced species side reactions. However, the mechanisms of carbon reactivity toward these species are still unclear. Here, we report a direct in situ X-ray photoelectron spectroscopy study of oxygen reduction by lithiated graphene and graphene-based materials. Although lithium peroxide (Li2O2) and lithium oxide (Li2O) reactions with carbon are thermodynamically favorable, neither of them was found to react even at elevated temperatures. As lithium superoxide is not stable at room temperature, potassium superoxide (KO2) prepared in situ was used instead to test the reactivity of graphene with superoxide species. In contrast to Li2O2 and Li2O, KO2 was demonstrated to be strongly reactive.

Many propositions have been already put forth for the practical use of N-graphene in various devices, such as batteries, sensors, ultracapacitors, and next generation electronics. However, the chemistry of nitrogen imperfections in this material still remains an enigma. Here we demonstrate a method to handle N-impurities in graphene, which allows efficient conversion of pyridinic N to graphitic N and therefore precise tuning of the charge carrier concentration. By applying photoemission spectroscopy and density functional calculations, we show that the electron doping effect of graphitic N is strongly suppressed by pyridinic N. As the latter is converted into the graphitic configuration, the efficiency of doping rises up to half of electron charge per N atom.

The ionization potential is a fundamental key quantity with great relevance to diverse material properties. We find that state of the art methods based on density functional theory and simple diagrammatic approaches as commonly taken in the GW approximation predict the ionization potentials of semiconductors and insulators unsatisfactorily. Good agreement between theory and experiment is obtained only when diagrams resulting from the antisymmetry of the many-electron wave function are taken into account via vertex corrections in the self-energy. The present approach describes both localized and delocalized states accurately, making it ideally suited for a wide class of materials and processes.

Electron-phonon coupling and the emergence of superconductivity in intercalated graphite have been studied extensively. Yet, phonon-mediated superconductivity has never been observed in the 2D equivalent of these materials, doped monolayer graphene. Here we perform angle-resolved photoemission spectroscopy to try to find an electron donor for graphene that is capable of inducing strong electron-phonon coupling and superconductivity. We examine the electron donor species Cs, Rb, K, Na, Li, Ca and for each we determine the full electronic band structure, the Eliashberg function and the superconducting critical temperature Tc from the spectral function. An unexpected low-energy peak appears for all dopants with an energy and intensity that depend on the dopant atom. We show that this peak is the result of a dopant-related vibration. The low energy and high intensity of this peak are crucially important for achieving superconductivity, with Ca being the most promising candidate for realizing superconductivity in graphene.