Abstract: Van der Waals quantum wells, naturally formed in two-dimensional layered materials with nanoscale thickness, possess many inherent advantages over conventional molecular beam epitaxy grown counterparts, and could bring up intriguing physics and applications. However, optical transitions originated from the series of quantized states in these emerging quantum wells are still elusive. Here, we show that multilayer black phosphorus appears to be an excellent candidate for van der Waals quantum wells with well-defined subbands and high optical quality. Using infrared absorption spectroscopy, we probe subband structures of multilayer black phosphorus with tens of atomic layers, revealing clear signatures for optical transitions with subband index as high as 10, far from what was attainable previously. Surprisingly, in addition to allowed transitions, an unexpected series of “forbidden” transitions is also evidently observed, which enables us to determine energy spacings separately for conduction and valence subbands. Furthermore, the linear tunability of subband spacings by temperature and strain is demonstrated. Our results are expected to facilitate potential applications for infrared optoelectronics based on tunable van der Waals quantum wells.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: Electrenes, an atomically thin form of layered electrides, are very recent members of the 2D materials family. In this work, we employed first principle calculations to determine stable, exfoliatable, and application-promising 2D electrene materials among possible M2X compounds, where M is a group II-A metal and X is a nonmetal element (C, N, P, As, and Sb). The promise of stable electrene compounds for battery applications is assessed via their exfoliation energy, adsorption properties, and migration energy barriers toward relevant Li, Na, K, and Ca atoms. Our calculations revealed five new stable electrene candidates in addition to previously known Ca2N and Sr2N. Among these seven dynamically stable electrenes, Ba2As, Ba2P, Ba2Sb, Ca2N, Sr2N, and Sr2P are found to be very promising for either K or Na ion batteries due to their extremely low migration energy barriers (5-16 meV), which roughly demonstrates 105 times higher mobility than graphene and two to four times higher mobility than other promising 2D materials such as MXene (Mo2C).

Abstract: High-performance thermoelectric materials are critical in recuperating the thermal losses in various machinery and promising in renewable energy applications. In this respect, the search for novel thermoelectric materials has attracted considerable attention. In particular, low dimensional materials have been proposed as potential candidates due to their unique and controllable thermal and electronic transport properties. The considerable potential of several two-dimensional materials as thermoelectric devices has already been uncovered and many new candidates that merit further research have been suggested. In this regard, we comprehensively investigate the thermoelectric coefficients and electronic fitness function (EFF) of a large family of structurally isotropic and anisotropic two-dimensional layered materials using density functional theory combined with semi-classical Boltzmann transport theory. With this high-throughput screening, we bring to light additional 2D crystals that haven't been previously classified as favorable TE materials. We predict that Pb2Se2, GeS2, As-2, NiS2, Hf2O6, Zr2O6, AsBrS, ISbTe, ISbSe, AsISe, and AsITe are promising isotropic thermoelectric materials due to their considerably high EFF values. In addition to these materials, Hf2Br4, Zr2Br4, Hf2Cl4, Zr2Cl4, Hf2O6, Zr(2)O(6)and Os(2)O(4)exhibit strong anisotropy and possess prominently high EFF values.

Abstract: Semiconductor monolayer transition metal dichalcogenides (TMDs) have brought a new paradigm by introducing optically addressable valley degree of freedom. Concomitantly, their high flexibility constitutes a unique platform that links optics to mechanics via valleytronics. With the intention to expedite the research in this direction, we investigated ten TMDs, namely MoS2, MoSe2, MoTe2, WS2, WSe2, WTe2, MoSSe, MoSeTe, WSSe, and WSeTe, which particularly includes their so-called janus types (JTMDs). First, we obtained their electronic band structures using regular and hybrid density functional theory (DFT) calculations in the presence of the spin-orbit coupling and biaxial or uniaxial strain. Our DFT results indicated that against the expectations based on their reported piezoelectric behavior, JTMDs typically interpolated between the standard band properties of the constituent TMDs without producing a novel feature. Next, by fitting to our DFT data we generated both spinless and spinful k center dot p parameter sets which are quite accurate over the K valley where the optical activity occurs. As an important application of this parametrization, we considered the circular and linear dichroism under strain. Among the studied (J)TMDs, WTe2 stood out with its largest linear dichroism under uniaxial strain because of its narrower band gap and large K valley uniaxial deformation potential. This led us to suggest WTe2 monolayer membranes for optical polarization-based strain measurements, or conversely, as strain tunable optical polarizers.

Abstract: The newest members of a two-dimensional material family, involving transition metal carbides and nitrides (called MXenes), have garnered increasing attention due to their tunable electronic and thermal properties depending on the chemical composition and functionalization. This flexibility can be exploited to fabricate efficient electrochemical energy storage (batteries) and energy conversion (thermoelectric) devices. In this study, we calculated the Seebeck coefficients and lattice thermal conductivity values of oxygen terminated M2CO2 (where M = Ti, Zr, Hf, Sc) monolayer MXene crystals in two different functionalization configurations (model-II (MD-II) and model-III (MD-III)), using density functional theory and Boltzmann transport theory. We estimated the thermoelectric figure-of-merit, zT, of these materials by two different approaches, as well. First of all, we found that the structural model (i.e. adsorption site of oxygen atom on the surface of MXene) has a paramount impact on the electronic and thermoelectric properties of MXene crystals, which can be exploited to engineer the thermoelectric properties of these materials. The lattice thermal conductivity kappa(l), Seebeck coefficient and zT values may vary by 40% depending on the structural model. The MD-III configuration always has the larger band gap, Seebeck coefficient and zT, and smaller kappa(l) as compared to the MD-II structure due to a larger band gap, highly flat valence band and reduced crystal symmetry in the former. The MD-III configuration of Ti2CO2 and Zr2CO2 has the lowest kappa(l) as compared to the same configuration of Hf2CO2 and Sc2CO2. Among all the considered structures, the MD-II configuration of Hf2CO2 has the highest kappa(l), and Ti2CO2 and Zr2CO2 in the MD-III configuration have the lowest kappa(l). For instance, while the band gap of the MD-II configuration of Ti2CO2 is 0.26 eV, it becomes 0.69 eV in MD-III. The zT(max) value may reach up to 1.1 depending on the structural model of MXene.

Abstract: Molybdenum carbide has various applications for which studying the material using classical molecular dynamics simulations would be valuable. Here, we develop an empirical potential within the Tersoff formalism using particle swarm optimization for the orthorhombic phase of Mo2C. The developed potential is shown to predict lattice constants, elastic properties, and equation of state results that are consistent with current and previously reported results from experiments and first principles calculations. We demonstrate the potential with simulations of indentation using multiple indenter sizes that load and unload in three different directions relative to the crystallographic lattice of orthorhombic Mo2C. Direction-dependent force-displacement trends are analyzed and explained in terms of the spatial distributions of stress and strain within the material during indentation. This study reveals the anisotropic elasticity of orthorhombic Mo2C and, more generally, provides researchers with a new empirical potential that can be used to explore the properties and behavior of the material going forward.

Abstract: We detail the effects of Cl and S functionalization on the superconducting properties of layered (bulk) and monolayer niobium carbide (Nb2C) MXene crystals, based on first-principles calculations combined with Eliashberg theory. For bulk layered Nb2CCl2, the calculated superconducting transition temperature (T-c) is in very good agreement with the recently measured value of 6 K. We show that T-c is enhanced to 10 K for monolayer Nb2CCl2, due to an increase in the density of states at the Fermi level, and the corresponding electron-phonon coupling. We further demonstrate feasible gate- and strain-induced enhancements of T-c for both bulk-layered and monolayer Nb2CCl2 crystals, resulting in T-c values of around 38 K. In the S-functionalized Nb2CCl2 crystals, our calculations reveal the importance of phonon softening in understanding their superconducting properties. Finally, we predict that Nb3C2S2 in bulk-layered and monolayer forms is also superconducting, with a T-c of around 28 K. Considering that Nb2C is not superconducting in pristine form, our findings promote functionalization as a pathway towards robust superconductivity in MXenes.

Abstract: Recent experiments demonstrated that proton transport through graphene electrodes can be accelerated by over an order of magnitude with low intensity illumination. Here we show that this photo-effect can be suppressed for a tuneable fraction of the infra-red spectrum by applying a voltage bias. Using photocurrent measurements and Raman spectroscopy, we show that such fraction can be selected by tuning the Fermi energy of electrons in graphene with a bias, a phenomenon controlled by Pauli blocking of photo-excited electrons. These findings demonstrate a dependence between graphene's electronic and proton transport properties and provide fundamental insights into molecularly thin electrode-electrolyte interfaces and their interaction with light. Recent experiments have shown that proton transport through graphene electrodes can be promoted by light, but the understanding of this phenomenon remains unclear. Here, the authors report the electrical tunability of this photo-effect, showing a connection between graphene electronic and proton transport properties.

Abstract: Defect-free graphene is impermeable to all atoms(1-5) and ions(6,7) under ambient conditions. Experiments that can resolve gas flows of a few atoms per hour through micrometre-sized membranes found that monocrystalline graphene is completely impermeable to helium, the smallest atom(2,5). Such membranes were also shown to be impermeable to all ions, including the smallest one, lithium(6,7). By contrast, graphene was reported to be highly permeable to protons, nuclei of hydrogen atoms(8,9). There is no consensus, however, either on the mechanism behind the unexpectedly high proton permeability(10-14) or even on whether it requires defects in graphene's crystal lattice(6,8,15-17). Here, using high-resolution scanning electrochemical cell microscopy, we show that, although proton permeation through mechanically exfoliated monolayers of graphene and hexagonal boron nitride cannot be attributed to any structural defects, nanoscale non-flatness of two-dimensional membranes greatly facilitates proton transport. The spatial distribution of proton currents visualized by scanning electrochemical cell microscopy reveals marked inhomogeneities that are strongly correlated with nanoscale wrinkles and other features where strain is accumulated. Our results highlight nanoscale morphology as an important parameter enabling proton transport through two-dimensional crystals, mostly considered and modelled as flat, and indicate that strain and curvature can be used as additional degrees of freedom to control the proton permeability of two-dimensional materials. A study using high-resolution scanning electrochemical cell microscopy attributes proton permeation through defect-free graphene and hexagonal boron nitride to transport across areas of the structure that are under strain.

Abstract: Utilizing magnetic tunnel junctions (MTJs) for write/read and fast spin-orbit-torque (SOT)-driven domain-wall (DW) motion for propagation, enables non-volatile logic and majority operations, representing a breakthrough in the implementation of nanoscale DW logic devices. Recently, current-driven DW logic gates have been demonstrated via magnetic imaging, where the Dzyaloshinskii-Moriya interaction (DMI) induces chiral coupling between perpendicular magnetic anisotropy (PMA) regions via an in-plane (IP) oriented region. However, full electrical operation of nanoscale DW logic requires electrical write/read operations and a method to pattern PMA and IP regions compatible with the fabrication of PMA MTJs. Here, we study the use of a Hybrid Free Layer (HFL) concept to combine an MTJ stack with DW motion materials, and He+ ion irradiation to convert the stack from PMA to IP. First, we investigate the free layer thickness dependence of 100-nm diameter HFL-MTJ devices and find an optimal CoFeB thickness, from 7 to 10 angstrom, providing high tunneling magnetoresistance (TMR) readout and efficient spin-transfer torque (STT) writing. We then show that high DMI materials, like Pt/Co, can be integrated into an MTJ stack via interlayer exchange coupling with the CoFeB free layer. In this design, DMI values suitable for SOT-driven DW motion are measured by asymmetric bubble expansion. Finally, we demonstrate that He+ irradiation reliably converts the coupled free layers from PMA to IP. These findings offer a path toward the integration of fully electrically controlled DW logic circuits.

Abstract: Spin-orbit torques (SOT) allow ultrafast, energy-efficient toggling of magnetization state by an in-plane charge current for applications such as magnetic random-access memory (SOT-MRAM). Tailoring the SOT vector comprising of antidamping (T-AD) and fieldlike (T-FL) torques could lead to faster, more reliable, and low-power SOT-MRAM. Here, we establish a method to quantify the longitudinal (T-AD) and transverse (T-FL) components of the SOT vector and its efficiency chi(AD) and chi(FL), respectively, in nanoscale three-terminal SOT magnetic tunnel junctions (SOT-MTJ). Modulation of nucleation or switching field (B-SF) for magnetization reversal by SOT effective fields (B-SOT) leads to the modification of SOT-MTJ hysteresis loop behavior from which chi(AD) and chi(FL) are quantified. Surprisingly, in nanoscale W/CoFeB SOT-MTJ, we find chi(FL) to be (i) twice as large as chi(AD) and (ii) 6 times as large as chi(FL) in micrometer-sized W/CoFeB Hall-bar devices. Our quantification is supported by micromagnetic and macrospin simulations which reproduce experimental SOT-MTJ Stoner-Wohlfarth astroid behavior only for chi(FL) > chi(AD). Additionally, from the threshold current for current-induced magnetization switching with a transverse magnetic field, we show that in SOT-MTJ, T-FL plays a more prominent role in magnetization dynamics than T-AD. Due to SOT-MRAM geometry and nanodimensionality, the potential role of nonlocal spin Hall spin current accumulated adjacent to the SOT-MTJ in the mediation of T-FL and chi(FL) amplification merits to be explored.

Abstract: We show how snake-orbit states which run along a magnetic edge can be confined electrically. We consider a two-dimensional electron gas (2DEG) confined into a quantum wire, subjected to a strong perpendicular and steplike magnetic field B/ − B. Close to this magnetic step, new, spatially confined bound states arise as a result of the lateral confinement and the magnetic-field step. The number of states, with energy below the first Landau level, increases as B becomes stronger or as the wire width becomes larger. These bound states can be understood as an interference between two counter-propagating one-dimensional snake-orbit states.

Abstract: We present an ab initio study of the electronic structure and of the formation energies of various point defects in BaSnO3 and SrGeO3. We show that La and Y impurities substituting Ba or Sr are shallow donors with a preferred 1 + charge state. These defects have a low formation energy within all the suitable equilibrium growth conditions considered. Oxygen vacancies behave as shallow donors as well, preferring the 2 + charge state. Their formation energies, however, are higher in most growth conditions, indicating a limited contribution to conductivity. The calculated electron effective mass in BaSnO3, with a value of 0.21 me, and the very high mobility reported recently in La-doped BaSnO3 single-crystals, suggest that remarkably low scattering rates can be achieved in the latter. In the case of SrGeO3, our results point to carrier density and mobility values in the low range for typical polycrystalline TCOs, in line with experiment.