Search

Although the unconventional $\pi$-magnetism at the zigzag edges of graphene holds promise for a wide array of applications, whether and to what degree it plays a role in their chemistry remains poorly understood. Here, we investigate the addition of a hydrogen atom $-$ the simplest yet the most experimentally relevant adsorbate $-$ to zigzag graphene nanoribbons (ZGNRs). We show that the $\pi$-magnetism governs the chemistry of ZGNRs, giving rise to a site-dependent reactivity of the carbon atoms and driving the hydrogenation process to the nanoribbon edges. Conversely, the chemisorbed hydrogen atom governs the $\pi$-magnetism of ZGNRs, acting as a spin-$\frac{1}{2}$ paramagnetic center in the otherwise antiferromagnetic ground state and spin-polarizing the charge carriers at the band extrema. Our findings establish a comprehensive picture of the peculiar interplay between chemistry and magnetism that emerges at the zigzag edges of graphene.

Twisted van der Waals multilayers are widely regarded as a rich platform to access novel electronic phases, thanks to the multiple degrees of freedom such as layer thickness and twist angle that allow control of their electronic and chemical properties. Here, we propose that the stacking domains that form naturally due to the relative twist between successive layers act as an additional "knob" for controlling the behavior of these systems, and report the emergence and engineering of stacking domain-dependent surface chemistry in twisted few-layer graphene. Using mid-infrared near-field optical microscopy and atomic force microscopy, we observe a selective adhesion of metallic nanoparticles and liquid water at the domains with rhombohedral stacking configurations of minimally twisted double bi- and tri-layer graphene. Furthermore, we demonstrate that the manipulation of nanoparticles located at certain stacking domains can locally reconfigure the moir\'e superlattice in their vicinity at the {\mu}m-scale. In addition, we report first-principles simulations of the energetics of adhesion of metal atoms and water molecules on the stacking domains in an attempt to elucidate the origin of the observed selective adhesion. Our findings establish a new approach to controlling moir\'e-assisted chemistry and nanoengineering., Comment: 11 pages, 3 figures

Vertical van der Waals heterostructures of semiconducting transition metal dichalcogenides realize moir\'e systems with rich correlated electron phases and moir\'e exciton phenomena. For material combinations with small lattice mismatch and twist angles as in MoSe$_2$-WSe$_2$, however, lattice reconstruction eliminates the canonical moir\'e pattern and instead gives rise to arrays of periodically reconstructed nanoscale domains and mesoscopically extended areas of one atomic registry. Here, we elucidate the role of atomic reconstruction in MoSe$_2$-WSe$_2$ heterostructures synthesized by chemical vapor deposition. With complementary imaging down to the atomic scale, simulations, and optical spectroscopy methods we identify the coexistence of moir\'e-type cores and extended moir\'e-free regions in heterostacks with parallel and antiparallel alignment. Our work highlights the potential of chemical vapor deposition for applications requiring laterally extended heterosystems of one atomic registry or exciton-confining heterostack arrays., Comment: 10 pages, 5 figures

The interlayer coupling in van der Waals heterostructures governs a variety of optical and electronic properties. The intrinsic dipole moment of Janus transition metal dichalcogenides (TMDs) offers a simple and versatile approach to tune the interlayer interactions. In this work, we demonstrate how the van der Waals interlayer coupling and charge transfer of Janus MoSSe/MoS2 heterobilayers can be tuned by the twist angle and interface composition. Specifically, the Janus heterostructures with a sulfur/sulfur (S/S) interface display stronger interlayer coupling than the heterostructures with a selenium/sulfur (Se/S) interface as shown by the low-frequency Raman modes. The differences in interlayer interactions are explained by the interlayer distance computed by density-functional theory (DFT). More intriguingly, the built-in electric field contributed by the charge density redistribution and interlayer coupling also play important roles in the interfacial charge transfer. Namely, the S/S and Se/S interfaces exhibit different levels of photoluminescence (PL) quenching of MoS2 A exciton, suggesting enhanced and reduced charge transfer at the S/S and Se/S interface, respectively. Our work demonstrates how the asymmetry of Janus TMDs can be used to tailor the interfacial interactions in van der Waals heterostructures.

The relativistic charge carriers in monolayer graphene can be manipulated in manners akin to conventional optics (electron-optics): angle-dependent Klein tunneling collimates an electron beam (analogous to a laser), while a Veselago refraction process focuses it (analogous to an optical lens). Both processes have been previously investigated, but the collimation and focusing efficiency have been reported to be relatively low even in state-of-the-art ballistic pn-junction devices. These limitations prevented the realization of more advanced quantum devices based on electron-optical interference, while understanding of the underlying physics remains elusive. Here, we present a novel device architecture of a graphene microcavity defined by carefully-engineered local strain and electrostatic fields. We create a controlled electron-optic interference process at zero magnetic field as a consequence of consecutive Veselago refractions in the microcavity and provide direct experimental evidence through low-temperature electrical transport measurements. The experimentally observed first-, second-, and third-order interference peaks agree quantitatively with the Veselago physics in a microcavity. In addition, we demonstrate decoherence of the interference by an external magnetic field, as the cyclotron radius becomes comparable to the interference length scale. For its application in electron-optics, we utilize Veselago interference to localize uncollimated electrons and characterize its contribution in further improving collimation efficiency. Our work sheds new light on relativistic single-particle physics and provides important technical improvements toward next-generation quantum devices based on the coherent manipulation of electron momentum and trajectory.

The coupling between electrons and phonons in solids plays a central role in describing many phenomena, including superconductivity and thermoelecric transport. Calculations of this coupling are exceedingly demanding as they necessitate integrations over both the electron and phonon momenta, both of which span the Brillouin zone of the crystal, independently. We present here an ab initio method for efficiently calculating electron-phonon mediated transport properties by dramatically accelerating the computation of the double integrals with a dual interpolation technique that combines maximally localized Wannier functions with symmetry-adapted plane waves. The performance gain in relation to the current state-of-the-art Wannier-Fourier interpolation is approximately 2n_s \times M, where n_s is the number of crystal symmetry operations and M, a number in the range 5 - 60, governs the expansion in star functions. We demonstrate with several examples how our method performs some ab initio calculations involving electron-phonon interactions., Comment: 9 pages, 5 figures

We introduce twisted trilayer graphene (tTLG) with two independent twist angles as an ideal system for the precise tuning of the electronic interlayer coupling to maximize the effect of correlated behaviors. As established by experiment and theory in the related twisted bilayer graphene system, van Hove singularities (VHS) in the density of states can be used as a proxy of the tendency for correlated behaviors. To explore the evolution of VHS in the twist-angle phase space of tTLG, we present a general low-energy electronic structure model for any pair of twist angles. We show that the basis of the model has infinite dimensions even at a finite energy cutoff and that no Brillouin zone exists even in the continuum limit. Using this model, we demonstrate that the tTLG system exhibits a wide range of magic angles at which VHS merge and that the density of states has a sharp peak at the charge-neutrality point through two distinct mechanisms: the incommensurate perturbation of twisted bilayer graphene's flatbands or the equal hybridization between two bilayer moiré superlattices.

The kagome lattice based on 3d transition metals is a versatile platform for novel topological phases hosting symmetry-protected electronic excitations and exotic magnetic ground states. However, the paradigmatic states of the idealized two-dimensional (2D) kagome lattice - Dirac fermions and topological flat bands - have not been simultaneously observed, partly owing to the complex stacking structure of the kagome compounds studied to date. Here, we take the approach of examining FeSn, an antiferromagnetic single-layer kagome metal with spatially-decoupled kagome planes. Using polarization- and termination-dependent angle-resolved photoemission spectroscopy (ARPES), we detect the momentum-space signatures of coexisting flat bands and Dirac fermions in the vicinity of the Fermi energy. Intriguingly, when complemented with bulk-sensitive de Haas-van Alphen (dHvA) measurements, our data reveal an even richer electronic structure that exhibits robust surface Dirac fermions on specific crystalline terminations. Through band structure calculations and matrix element simulations, we demonstrate that the bulk Dirac bands arise from in-plane localized Fe-3d orbitals under kagome symmetry, while the surface state realizes a rare example of fully spin-polarized 2D Dirac fermions when combined with spin-layer locking in FeSn. These results highlight FeSn as a prototypical host for the emergent excitations of the kagome lattice. The prospect to harness these excitations for novel topological phases and spintronic devices is a frontier of great promise at the confluence of topology, magnetism, and strongly-correlated electron physics.

We investigate the effects of lithium intercalation in twisted bilayers of graphene, using first-principles electronic structure calculations. To model this system we employ commensurate supercells that correspond to twist angles of 7.34$^\circ$ and 2.45$^\circ$. From the energetics of lithium absorption we demonstrate that for low Li concentration the intercalants cluster in the AA regions with double the density of a uniform distribution. The charge donated by the Li atoms to the graphene layers results in modifications to the band structure that can be qualitatively captured using a continuum model with modified interlayer couplings in a region of parameter space that has yet to be explored either experimentally or theoretically. Thus, the combination of intercalation and twisted layers simultaneously provides the means for spatial control over material properties and an additional knob with which to tune moir\'e physics in twisted bilayers of graphene, with potential applications ranging from energy storage and conversion to quantum information., Comment: 7 pages, 4 figures

The structural relaxation of multilayer graphene is essential in describing the interesting electronic properties induced by intentional misalignment of successive layers, including the recently reported superconductivity in twisted bilayer graphene. This is difficult to accomplish without an accurate interatomic potential. Here, we present a new, registry-dependent Kolmogorov-Crespi type interatomic potential to model interlayer interactions in multilayer graphene structures. It consists of two parts representing attractive interaction due to dispersion, and repulsive interaction due to anisotropic overlap of electronic orbitals. An important new feature is a dihedral-angle-dependent term that is added to the repulsive part in order to describe correctly several distinct stacking states that the original Kolmogorov-Crespi potential cannot distinguish. We refer to the new model as the Dihedral-angle-corrected Registry-dependent Interlayer Potential (DRIP). Computations for several test problems show that DRIP correctly reproduces the binding, sliding, and twisting energies and forces obtained from ab initio total-energy calculations based on density functional theory. We use the new potential to study the structural properties of a twisted graphene bilayer and the exfoliation of graphene from graphite. Our potential is available through the OpenKIM interatomic potential repository at https://openkim.org., Comment: Reparameterized against a different training set. The new parameterization gives better predictions for interlayer binding energy and c-axis elastic modulus

We study the quantum nonlinear Hall effect in two-dimensional materials with time-reversal symmetry. When only one mirror line exists, a transverse charge current occurs in second-order response to an external electric field, as a result of the Berry curvature dipole in momentum space. Candidate 2D materials to observe this effect are two-dimensional transition-metal dichalcogenides~(TMDCs). First we use an ab initio based tight-binding approach to demonstrate that monolayer $T_d$-stricture TMDCs exhibit a finite Berry curvature dipole. In the $1H$ and $1T'$ phase of TMDCs, we show the emergence of finite Berry curvature dipole with the application of strain and electrical displacement field respectively., Comment: 4 figures

Van der Waals (vdW) heterostructures are an emergent class of metamaterials comprised of vertically stacked two-dimensional (2D) building blocks, which provide us with a vast tool set to engineer their properties on top of the already rich tunability of 2D materials. One of the knobs, the twist angle between different layers, plays a crucial role in the ultimate electronic properties of a vdW heterostructure and does not have a direct analog in other systems such as MBE-grown semiconductor heterostructures. For small twist angles, the moir\'e pattern produced by the lattice misorientation creates a long-range modulation. So far, the study of the effect of twist angles in vdW heterostructures has been mostly concentrated in graphene/hexagonal boron nitride (h-BN) twisted structures, which exhibit relatively weak interlayer interaction due to the presence of a large bandgap in h-BN. Here we show that when two graphene sheets are twisted by an angle close to the theoretically predicted 'magic angle', the resulting flat band structure near charge neutrality gives rise to a strongly-correlated electronic system. These flat bands exhibit half-filling insulating phases at zero magnetic field, which we show to be a Mott-like insulator arising from electrons localized in the moir\'e superlattice. These unique properties of magic-angle twisted bilayer graphene (TwBLG) open up a new playground for exotic many-body quantum phases in a 2D platform made of pure carbon and without magnetic field. The easy accessibility of the flat bands, the electrical tunability, and the bandwidth tunability though twist angle may pave the way towards more exotic correlated systems, such as unconventional superconductors or quantum spin liquids., Comment: Main Text: 14 pages, 4 figures. Supplemental Information: 24 pages, 14 figures

Wave polarization contains valuable information for electromagnetic signal processing and the ability to manipulate it can be extremely useful in photonic devices. In this work, we propose designs comprised of one of the emerging and interesting two-dimensional media: Black Phosphorus. Due to substantial in-plane anisotropy, a single slab of Black Phosphorus can be very efficient for manipulating the polarization state of electromagnetic waves. We investigate Black Phosphorus slabs that filter the fields along one direction, or polarization axis rotation, or convert linear polarization to circular. These slabs can be employed as components in numerous mid-IR integrated devices.

We derive electronic tight-binding Hamiltonians for strained graphene, hexagonal boron nitride and transition metal dichalcogenides based on Wannier transformation of {\it ab initio} density functional theory calculations. Our microscopic models include strain effects to leading order that respect the hexagonal crystal symmetry and local crystal configuration, and are beyond the central force approximation which assumes only pair-wise distance dependence. Based on these models, we also derive and analyze the effective low-energy Hamiltonians. Our {\it ab initio} approaches complement the symmetry group representation construction for such effective low-energy Hamiltonians and provide the values of the coefficients for each symmetry-allowed term. These models are relevant for the design of electronic device applications, since they provide the framework for describing the coupling of electrons to other degrees of freedom including phonons, spin and the electromagnetic field. The models can also serve as the basis for exploring the physics of many-body systems of interesting quantum phases., Comment: 4 figures; the code for generating the strain TMDC Hamiltonians will be available on the research group website

The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nano-scale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domain. To address this problem, we developed a multi-scale mass transport computational framework based on the Lattice Boltzmann Method (LBM) that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

The negatively charged nitrogen-vacancy defect ($N-V^-$ center) in diamond is an important atomic-scale structure that can be used as a qubit in quantum computing and as a marker in biomedical applications. Its usefulness relies on the ability to optically excite electrons between well-defined gap states, which requires clear and detailed understanding of the relevant states and excitation processes. Here we show that by using hybrid density-functional-theory calculations in a large supercell we can reproduce the zero-phonon line and the Stokes and anti-Stokes shifts, yielding a complete picture of the spin-conserving excitation of this defect., Comment: 4 pages, 2 figures

Magnetic order in graphene-related structures can arise from size effects or from topological frustration. We introduce a rigorous classification scheme for the types of finite graphene structures (nano-flakes) which lead to large net spin or to antiferromagnetic coupling between groups of electron spins. Based on this scheme, we propose specific examples of structures that can serve as the fundamental (NOR and NAND) logic gates for the design of high-density ultra-fast spintronic devices. We demonstrate, using ab initio electronic structure calculations, that these gates can in principle operate at room temperature with very low and correctable error rates., Comment: Typo in title fixed

The nitrogen-vacancy center in diamond is a promising candidate for realizing the spin qubits concept in quantum information. Even though this defect is known for a long time, its electronic structure and other properties have not yet been explored in detail. We study the properties of the nitrogen-vacancy center in diamond through density functional theory within the local spin density approximation, using supercell calculations. While this theory is strictly applicable for ground state properties, we are able to give an estimate for the energy sequence of the excited states of this defect. We also calculate the hyperfine tensors in the ground state. The results clearly show that: (i) the spin density and the appropriate hyperfine constants are spread along a plane and unevenly distributed around the core of the defect; (ii) the measurable hyperfine constants can be found within about 7 A from the vacancy site. These results have important implications on the decoherence of the electron spin which is crucial in realizing the spin qubits in diamond., Comment: 28 pages, 7 figures, 2 tables

The recently discovered FeAs-based superconductors show intriguing behavior and unusual dynamics of electrons and holes which occupy the Fe $d$-orbitals and As $4s$ and $4p$ orbitals. Starting from the atomic limit, we carry out a strong coupling expansion to derive an effective hamiltonian that describes the electron and hole behavior. The hopping and the hybridization parameters between the Fe $d$ and As $s$ and $p$-orbitals are obtained by fitting the results of our density-functional-theory calculations to a tight-binding model with nearest-neighbor interactions and a minimal orbital basis. We find that the effective hamiltonian, in the strong on-site Coulomb repulsion limit, operates on three distinct sub-spaces coupled through Hund's rule. The three sub-spaces describe different components (or subsystems): (a) one spanned by the $d_{x^2-y^2}$ Fe orbital; (b) one spanned by the degenerate atomic Fe orbitals $d_{xz}$ and $d_{yz}$; and (c) one spanned by the atomic Fe orbitals $d_{xy}$ and $d_{z^2}$. Each of these Hamiltonians is an extended t-t'-J-J' model and is characterized by different coupling constants and filling factors. For the case of the undoped material the second subspace alone prefers a ground state characterized by a spin-density-wave order similar to that observed in recent experimental studies, while the other two subspaces prefer an antiferromagnetic order. We argue that the observed spin-density-wave order minimizes the ground state energy of the total hamiltonian., Comment: 11 pages, 6 figures

A multiscale approach is used to simulate the translocation of DNA through a nanopore. Within this scheme, the interactions of the molecule with the surrounding fluid (solvent) are explicitly taken into account. By generating polymers of various initial configurations and lengths we map the probability distibutions of the passage times of the DNA through the nanopore. A scaling law behavior for the most probable of these times with respect to length is derived, and shown to exhibit an exponent that is in a good agreement with the experimental findings. The essential features of the DNA dynamics as it passes through the pore are explored., Comment: 5 pages, 4 figures. to appear in International Journal of Modern Physics C

We employ a multiscale approach to model the translocation of biopolymers through nanometer size pores. Our computational scheme combines microscopic Molecular Dynamics (MD) with a mesoscopic Lattice Boltzmann (LB) method for the solvent dynamics, explicitly taking into account the interactions of the molecule with the surrounding fluid. We describe an efficient parallel implementation of the method which exhibits excellent scalability on the Blue Gene platform. We investigate both dynamical and statistical aspects of the translocation process by simulating polymers of various initial configurations and lengths. For a representative molecule size, we explore the effects of important parameters that enter in the simulation, paying particular attention to the strength of the molecule-solvent coupling and of the external electric field which drives the translocation process. Finally, we explore the connection between the generic polymers modeled in the simulation and DNA, for which interesting recent experimental results are available., Comment: 10 pages, 5 figures, 3 tables. to appear in International Journal for Multiscale Computational Engineering

We present a formalism for coupling a density functional theory-based quantum simulation to a classical simulation for the treatment of simple metallic systems. The formalism is applicable to multiscale simulations in which the part of the system requiring quantum-mechanical treatment is spatially confined to a small region. Such situations often arise in physical systems where chemical interactions in a small region can affect the macroscopic mechanical properties of a metal. We describe how this coupled treatment can be accomplished efficiently, and we present a coupled simulation for a bulk aluminum system., Comment: 15 pages, 7 figures

We present first-principles calculations in the framework of density-functional theory and the pseudopotential approach, aiming to model the intermediate stages of the reduction of NO in the presence of MoO$_3$(010). In particular, we study the formation of dinitrosyl, which proves to be an important intermediate stage in the catalytic reduction. We find that the replacement of an oxygen of MoO$_3$ by NO is energetically favorable, and that the system lowers further its energy by the formation of (NO)$_2$. Moreover, the geometry and charge distribution for the adsorbed dinitrosyl indicates a metal-oxide mediated coupling between the two nitrogen and the two oxygen atoms. We discuss the mechanisms for the dinitrosyl formation and the role of the oxide in the reaction., Comment: 6 pages, 4 figs, RevTeX. To be published in J. Chem. Phys

The interaction of vacancy with dislocations in Al is studied using the Semidiscrete Variational Peierls-Nabarro model with ab initio determined gamma surface. For the first time, we confirm theoretically the so-called vacancy lubrication effect on dislocation motion in Al, a discovery that can settle a long-standing controversy in dislocation theory for fcc metals. We provide insights on the lubrication effect by exploring the connection between dislocation mobility and its core width. We predict an increased dislocation splitting in the presence of vacancy. We find that on average there is a weak repulsion between vacancies and dislocations which is independent of dislocation character., Comment: 6 pages, 2 figures

We investigate the effects of geometrical micro-irregularities on the conversion efficiency of reactive flows in narrow channels of millimetric size. Three-dimensional simulations, based upon a Lattice-Boltzmann-Lax-Wendroff code, indicate that periodic micro-barriers may have an appreciable effect on the effective reaction efficiency of the device. Once extrapolated to macroscopic scales, these effects can result in a sizeable increase of the overall reaction efficiency., Comment: 12 pages, 12 figures

The role of geometrical micro-barriers on the conversion efficiency of reactive flows in narrow three-dimensional channels of millimetric size is investigated. Using a Lattice-Boltzmann-Lax-Wendroff code, we show that micro-barriers have an appreciable effect on the effective reaction efficiency of the device. If extrapolated to macroscopic scales, these effects can result in a sizeable increase of the overall reaction efficiency., Comment: 5 pages, 7 figures

We calculate the free energies of unstable stacking fault (USF) configurations on the glide and shuffle slip planes in silicon as a function of temperature, using the recently developed Environment Dependent Interatomic Potential (EDIP). We employ the molecular dynamics (MD) adiabatic switching method with appropriate periodic boundary conditions and restrictions to atomic motion that guarantee stability and include volume relaxation of the USF configurations perpendicular to the slip plane. Our MD results using the EDIP model agree fairly well with earlier first-principles estimates for the transition from shuffle to glide plane dominance as a function of temperature. We use these results to make contact to brittle-ductile transition models., Comment: 6 pages revtex, 4 figs, 16 refs, to appear in Phys. Rev. B

We present a study of the equilibrium properties of $sp$-bonded solids within the pseudopotential approach, employing recently proposed generalized gradient approximation (GGA) exchange correlation functionals. We analyze the effects of the gradient corrections on the behavior of the pseudopotentials and discuss possible approaches for constructing pseudopotentials self-consistently in the context of gradient corrected functionals. The calculated equilibrium properties of solids using the GGA functionals are compared to the ones obtained through the local density approximation (LDA) and to experimental data. A significant improvement over the LDA results is achieved with the use of the GGA functionals for cohesive energies. For the lattice constant, the same accuracy as in LDA can be obtained when the nonlinear coupling between core and valence electrons introduced by the exchange correlation functionals is properly taken into account. However, GGA functionals give bulk moduli that are too small compared to experiment., Comment: 15 pages, latex, no figures

We present a method for electronic structure calculations that retains all of the advantages of real space and addresses the inherent inefficiency of a regular grid, which has equal precision everywhere. The computations are carried out on a {\em regular} mesh in {\em curvilinear space}, which allows natural and efficient decomposition on parallel computers, and effective use of iterative numerical methods. A novel feature is the use of error analysis to optimize the curvilinear grid for highly inhomogeneous electronic distributions. We report accurate all-electron calculations for H$_2$, O, and O$_{2}$., Comment: 4 pages, 2 Postscript figures