Search

We use photo-electron paramagnetic resonance (EPR) measurements and first-principles calculations to identify and explain the properties of carbon in AlN. We present clear evidence for carbon substitution on the nitrogen site (C$_{\rm N}$). We also clarify the origin of a widely observed EPR spectra in AlN that, although often attributed to a deep donor defect, we demonstrate is surprisingly due to C$_{\rm N}$. Finally, we show the presence of C$_{\rm N}$ is consistent with the absorption spectra at 4.7 eV observed in AlN., Comment: 5 pages, 3 figures

We find the recently developed strongly constrained and appropriately normed (SCAN) functional, now widely used in calculations of many materials, is not able to reliably describe the properties of deep defects and small polarons in a set of wide-bandgap semiconductors and insulators (ZnO, ZnSe, GaN, Ga$_2$O$_3$, and NaF). By comparing first-principles calculations using the SCAN functional against established experimental information and first-principles calculations using a hybrid functional, we find that the SCAN functional systematically underestimates the magnitude of the structural distortions at deep defects and tends to delocalize the charge density of these defect states., Comment: 10 pages, 8 figures

Chalcogen vacancies in the semiconducting monolayer transition-metal dichalcogenides (TMDs) have frequently been invoked to explain a wide range of phenomena, including both unintentional p-type and n-type conductivity, as well as sub-band gap defect levels measured via tunneling or optical spectroscopy. These conflicting interpretations of the deep versus shallow nature of the chalcogen vacancies are due in part to shortcomings in prior first-principles calculations of defects in the semiconducting two-dimensional (2D) TMDs that have been used to explain experimental observations. Here we report results of hybrid density functional calculations for the chalcogen vacancy in a series of monolayer TMDs, correctly referencing the thermodynamic charge transition levels to the fundamental band gap (as opposed to the optical band gap). We find that the chalcogen vacancies are deep acceptors and cannot lead to n-type or p-type conductivity. Both the (0/$-1$) and ($-$1/$-$2) transition levels occur in the gap, leading to paramagnetic charge states S=1/2 and S=1, respectively, in a collinear-spin representation. We discuss trends in terms of the band alignments between the TMDs, which can serve as a guide to future experimental studies of vacancy behavior.

The recent discovery of Ising superconductivity has garnered a lot of interest due in part to the resilience of these superconductors to large in-plane magnetic fields. In this Perspective we explain the basic concepts that define the behavior of Ising superconductors, provide an overview of the electronic structure and magnetic properties with a focus on NbSe$_2$, summarize key experimental observations that have been made in this class of superconductors, highlight the role that defects and proximity-induced effects at interfaces have on Ising superconductivity and finally discuss the prospects for observing Ising superconductivity in bulk materials., Comment: 10 pages, 8 figures

Low-dimensional materials have remarkable properties that are distinct from their bulk counterparts. A paradigmatic example is Ising superconductivity that occurs in monolayer materials such as NbSe2 which show a strong violation of the Pauli limit. In monolayers, this occurs due to a combination of broken inversion symmetry and spin-orbit coupling that locks the spins of the electrons out-of-plane. Bulk NbSe2 is centrosymmetric and is therefore not an Ising superconductor. We show that bulk misfit compound superconductors, (LaSe)1.14(NbSe2) and (LaSe)1.14(NbSe2)2, comprised of monolayers and bilayers of NbSe2, exhibit unexpected Ising protection with a Pauli-limit violation comparable to monolayer NbSe2, despite formally having inversion symmetry. We study these misfit compounds using complementary experimental methods in combination with first-principles calculations. We propose theoretical mechanisms of how the Ising protection can survive in bulk materials. We show how some of these mechanisms operate in these bulk compounds due to a concerted effect of charge-transfer, defects, reduction of interlayer hopping, and stacking. This highlights how Ising superconductivity can, unexpectedly, arise in bulk materials, and possibly enable the design of bulk superconductors that are resilient to magnetic fields., Comment: 14 pages, 4 figures, 4 extended data items

Transition metal dichalcogenides offer unprecedented versatility to engineer 2D materials with tailored properties to explore novel structural and electronic phase transitions. In this work, we present the atomic-scale evolution of the electronic ground state of a monolayer of Nb$_{1-\delta}$Mo$_{\delta}$Se$_2$ across the entire alloy composition range (0 < ${\delta}$ < 1) using low-temperature (300 mK) scanning tunneling microscopy and spectroscopy (STM/STS). In particular, we investigate the atomic and electronic structure of this 2D alloy throughout the metal to semiconductor transition (monolayer NbSe$_2$ to MoSe$_2$). Our measurements let us extract the effective doping of Mo atoms, the bandgap evolution and the band shifts, which are monotonic with ${\delta}$. Furthermore, we demonstrate that collective electronic phases (charge density wave and superconductivity) are remarkably robust against disorder. We further show that the superconducting TC changes non-monotonically with doping. This contrasting behavior in the normal and superconducting state is explained using first-principles calculations. We show that Mo doping decreases the density of states at the Fermi level and the magnitude of pair-breaking spin fluctuations as a function of Mo content. Our results paint a detailed picture of the electronic structure evolution in 2D TMD alloys, which is of utmost relevance for future 2D materials design.

Characterizing and manipulating the circular polarization of light is central to numerous emerging technologies, including spintronics and quantum computing. Separately, monolayer tungsten disulfide (WS2) is a versatile material that has demonstrated promise in a variety of applications, including single photon emitters and valleytronics. Here, we demonstrate a method to tune the photoluminescence (PL) intensity (factor of x161), peak position (38.4meV range), circular polarization (39.4% range), and valley polarization of a Bi2Se3-WS2 2D heterostructure using a low-power laser (0.762uW) in ambient. Changes are spatially confined to the laser spot, enabling submicron (814nm) features, and are long-term stable (>334 days). PL and valley polarization changes can be controllably reversed through laser exposure in vacuum, allowing the material to be erased and reused. Atmospheric experiments and first-principles calculations indicate oxygen diffusion modulates the exciton radiative vs. non-radiative recombination pathways, where oxygen absorption leads to brightening, and desorption to darkening., Comment: Published: ACS Appl. Mater. Interfaces 2022

We conducted a tip-enhanced Raman scattering spectroscopy (TERS) and photoluminescence (PL) study of quasi-1D TaSe3 nanoribbons exfoliated onto gold substrates. At a selenium deficiency of ~0.25 (Se/Ta=2.75,), the nanoribbons exhibit a strong, broad PL peak centered around ~920 nm (1.35 eV), suggesting their semiconducting behavior. Such nanoribbons revealed a strong TERS response under 785-nm laser excitation, allowing for their nanoscale spectroscopic imaging. Nanoribbons with a smaller selenium deficiency of ~0.15 (Se/Ta=2.85) did not show any PL or TERS response. The confocal Raman spectra of these samples agree with the previously-reported spectra of metallic TaSe3. The differences in the optical response of the nanoribbons examined in this study suggest that even small variations in Se content can induce changes in electronic structure, causing samples to exhibit either metallic or semiconducting character. The temperature-dependent electrical measurements of devices fabricated with both types of materials corroborate these observations. The density-functional-theory calculations revealed that incorporation of an oxygen atom in a Se vacancy can result in band gap opening and thus enable the transition from a metal to a semiconductor. However, the predicted bandgap is substantially smaller than that derived from PL data. These results indicate that the properties of van der Waals materials can vary significantly depending on stoichiometry, defect types and concentration, and possibly environmental and substrate effects. In view of this finding, local probing of nanoribbon properties with TERS becomes essential to understanding such low-dimensional systems., Comment: 24 pages, 8 figures

We use hybrid density functional calculations to assess n-type doping in monoclinic (Al$_x$Ga$_{1-x}$)$_2$O$_3$ alloys. We focus on Si, the most promising donor dopant, and study the structural properties, formation energies and charge-state transition levels of its various configurations. We also explore the impact of C and H, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga$_2$O$_3$, Si$_{Ga}$ is an effective shallow donor, but in Al$_2O_3$ Si$_{Al}$ acts as a DX center with a (+/-) transition level in the band gap. Interstitial H acts as a shallow donor in Ga$_2$O$_3$, but behaves as a compensating acceptor in n-type Al$_2O_3$. Interpolation indicates that Si is an effective donor in (Al$_x$Ga$_{1-x}$)$_2$O$_3$ up to 70% Al, but it can be compensated by H already at 1% Al. We also assess the diffusivity of H and study complex formation. Si$_{cation}$-H complexes have relatively low binding energies. Substitutional C on a cation site acts as a shallow donor in Ga$_2$O$_3$, but can be stable in a negative charge state in (Al$_x$Ga$_{1-x}$)$_2$O$_3$ when x>5%. Substitutional C on an O site (C$_O$) always acts as an acceptor in n-type (Al$_x$Ga$_{1-x}$)$_2$O$_3$, but will incorporate only under relatively O-poor conditions. C$_O$-H complexes can actually incorporate more easily, explaining observations of C-related compensation in Ga$_2$O$_3$ grown by MOCVD. We also investigate C$_{cation}$-H complexes, finding they have high binding energies and act as compensating acceptors when x>56%; otherwise the H just passivates the unintentional C donors. C-H complex formation explains why MOCVD grown Ga$_2$O$_3$ can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of C. Our study highlights that, while Si is a suitable shallow donor in ALGO alloys, control of unintentional impurities is essential to avoid compensation., Comment: 17 pages, 13 figures

The twist angle between the monolayers in van der Waals heterostructures provides a new degree of freedom in tuning material properties. We compare the optical properties of WSe2 homobilayers with 2H and 3R stacking using photoluminescence, Raman spectroscopy, and reflectance contrast measurements under ambient and cryogenic temperatures. Clear stacking- dependent differences are evident for all temperatures, with both photoluminescence and reflectance contrast spectra exhibiting a blue shift in spectral features in 2H compared to 3R bilayers. Density functional theory (DFT) calculations elucidate the source of the variations and the fundamental differences between 2H and 3R stackings. DFT finds larger energies for both A and B excitonic features in 2H than in 3R, consistent with experimental results. In both stacking geometries, the intensity of the dominant A1g Raman mode exhibits significant changes as a function of laser excitation wavelength. These variations in intensity are intimately linked to the stacking- and temperature-dependent optical absorption through resonant enhancement effects. The strongest enhancement is achieved when the laser excitation coincides with the C excitonic feature, leading to the largest Raman intensity under 514 nm excitation in 2H stacking and at 520 nm in 3R stacked WSe2 bilayers.

NbSe$_{2}$ and NbS$_{2}$ are isostructural two-dimensional materials that exhibit contrasting superconducting properties when reduced to the single monolayer limit. Monolayer NbSe$_{2}$ is an Ising superconductor, while there have been no reports of superconductivity in monolayer NbS$_{2}$. NbS$_{x}$Se$_{2-x}$ alloys exhibit an intriguing non-monotonic dependence of the superconducting transition temperature with sulfur content, which has been interpreted as a manifestation of fractal superconductivity. However, several key questions about this result are not known: (1) Does the electronic structure of the alloy differ from the parent compounds, (2) Are spin fluctuations which have been shown to be prominent in monolayer NbSe$_{2}$ also present in the alloys? Using first-principles calculations, we show that the density of states at the Fermi level and the proximity to magnetism in NbS$_{x}$Se$_{2-x}$ alloys are both reduced compared to the parent compound; the former would decrease the transition temperature while the latter would increase it. We also show that Se vacancies, which are likely magnetic pair-breaking defects, may form in large concentrations in NbSe$_{2}$. Based on our results, we suggest an alternative explanation of the non-monotonic behavior the superconducting transition temperature in NbS$_{x}$Se$_{2-x}$ alloys, which does not require the conjecture of multifractality., Comment: [v1] had minor typos in Figure 1 and the main text that are now fixed. Please use the corrected version, [v2]

Hybrid Ising superconductor-ferromagnetic insulator heterostructures provide a unique opportunity to explore the interplay between proximity-induced magnetism, spin-orbit coupling and superconductivity. Here we use a combination of first-principles calculations of NbSe$_{2}$/CrBr$_{3}$ heterostructures and an analytical theory of Ising superconductivity to analyze the existing experiments and provide a complete explanation of highly nontrivial and largely counterintuitive effects: an increase in the magnitude of the superconducting gap accompanied by the broadening of the tunneling peaks; hysteretic behavior of the tunneling conductance that sets in $\approx 2$ K below $T_c$; and nematic symmetry breaking in the superconducting state. The microscopic reason in all three cases appears to be the interplay between the proximity-induced exchange splitting and intrinsic defects. Finally, we predict additional interesting effects that at the moment cannot be addressed experimentally: spin-filtering when tunneling across CrBr$_{3}$ and tunneling ``hot spots'' in momentum space that are anticorrelated with regions where the spin-orbit splitting is maximum.

Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. [Phys. Rev. B 90, 075202 (2014)] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture process. An approach for evaluating electron-phonon coupling within the projector augmented wave formalism is presented. We also show that the common procedure of replacing Dirac delta functions with Gaussians can introduce errors into the resulting capture rate, and implement an alternative scheme to properly account for vibrational broadening. Lastly, we assess the accuracy of using an analytic approximation to the Sommerfeld parameter by comparing with direct numerical evaluation.

A variety of experiments have been carried out to establish the origin of the chiral charge-density wave transition in 1T-TiSe$_2$, which in turn has led to contradictory conclusions on the origin of this transition. Some studies suggest the transition is a phonon-driven structural distortion while other studies suggest it is an excitonic insulator phase transition that is accompanied by a lattice distortion. First, we propose these interpretations can be reconciled if one analyzes the available experimental and theoretical data within a formal definition of what constitutes an excitonic insulator as initially proposed by Keldysh and Kopaev. Next, we present pump-probe measurements of circularly polarized optical transitions and first-principles calculations where we highlight the importance of accounting for structural distortions to explain the finite chirality of optical transitions in the CDW phase. We show that at the elevated electronic temperature that occurs upon photoexcitation, there is a non-centrosymmetric structure that is near-degenerate in energy with the centrosymmetric charge density wave structure, which explains the finite chirality of the optical transitions observed in the CDW phase of TiSe$_2$., Comment: 9 pages, 5 figures

Recent studies on superconductivity in NbSe$_2$ have demonstrated a large anisotropy in the superconducting critical field when the material is reduced to a single monolayer. Motivated by this recent discovery, we use density functional theory (DFT) calculations to quantitatively address the superconducting properties of bulk and monolayer NbSe$_2$. We demonstrate that NbSe$_2$ is close to a ferromagnetic instability, and analyze our results in the context of experimental measurements of the spin susceptibility in NbSe$_2$. We show how this magnetic instability, which is pronounced in a single monolayer, can enable sizeable singlet-triplet mixing of the superconducting order parameter, contrary to contemporary considerations of the pairing symmetry in monolayer NbSe$_2$, and discuss approaches as to how this degree of mixing can be addressed quantitatively within our DFT framework. Our calculations also enable a quantitative description of the large anisotropy of the superconducting critical field, using DFT calculations of monolayer NbSe$_2$ in the normal state, Comment: 13 pages, 6 figures

Using first-principles calculations we determine the role of compressive and tensile uniaxial and equibiaxial strain on the structural, electronic and magnetic properties of V$_2$O$_3$. We find that compressive strain increases the energy cost to transition from the high-temperature paramagnetic metallic phase to the low-temperature antiferromagnetic insulating phase. This shift in the energy difference can be explained by changes in the V-V bond lengths that are antiferromagnetically aligned in the low temperature structure. The insights that we have obtained provide a microscopic explanation for the shifts in the metal-insulator transition temperature that have been observed in experiments of V$_2$O$_3$ films grown on different substrates.

The vanadates VO$_2$ and V$_2$O$_3$ are prototypical examples of strongly correlated materials that exhibit a metal-insulator transition. While the phase transitions in these materials have been studied extensively, there is a limited understanding of how the properties of these materials are affected by the presence of defects and doping. In this study we investigate the impact of native point defects in the form of Frenkel defects on the structural, magnetic and electronic properties of VO$_2$ and V$_2$O$_3$, using first-principles calculations. In VO$_2$ the vanadium Frenkel pairs lead to a non-trivial insulating state. The unpaired vanadium interstitial bonds to a single dimer, which leads to a trimer that has one singlet state and one localized single-electron $S=1/2$ state. The unpaired broken dimer created by the vanadium vacancy also has a localized $S=1/2$ state. Thus, the insulating state is created by the singlet dimers, the trimer and the two localized $S=1/2$ states. Oxygen Frenkel pairs, on the other hand, lead to a metallic state in VO$_2$, but are expected to be present in much lower concentrations. In contrast, the Frenkel defects in V$_2$O$_3$ do not directly suppress the insulating character of the material. However, the disorder created by defects in V$_2$O$_3$ alters the local magnetic moments and in turn reduces the energy cost of a transition between the insulating and conducting phases of the material. We also find self-trapped small polarons in V$_2$O$_3$, which has implications for transport properties in the insulating phase., Comment: 10 pages, 6 figures

Deep level transient spectroscopy (DLTS) is used extensively to study defects in semiconductors. We demonstrate that great care should be exercised in interpreting activation energies extracted from DLTS as ionization energies. We show how first-principles calculations of thermodynamic transition levels, temperature effects of ionization energies, and nonradiative capture coefficients can be used to accurately determine actual activation energies that can be directly compared with DLTS. Our analysis is illustrated with hybrid functional calculations for two important defects in GaN that have similar thermodynamic transition levels, and shows that the activation energy extracted from DLTS includes a capture barrier that is temperature dependent, unique to each defect, and in some cases large in comparison to the ionization energy. By calculating quantities that can be directly compared with experiment, first-principles calculations thus offer powerful leverage in identifying the microscopic origin of defects detected in DLTS., Comment: Main text and Supplementary Material

The singular density of states and the two Fermi wavevectors resulting from a ring-shaped or "Mexican hat" valence band give rise to unique trends in the charged impurity scattering rates and charged impurity limited mobilities. Ring shaped valence bands are common features of many monolayer and few-layer two-dimensional materials including the III-VI materials GaS, GaSe, InS, and InSe. The wavevector dependence of the screening, calculated within the random phase approximation, is so strong that it is the dominant factor determining the overall trends of the scattering rates and mobilities with respect to temperature and hole density. Charged impurities placed on the substrate and in the 2D channel are considered. The different wavevector dependencies of the bare Coulomb potentials alter the temperature dependence of the mobilities. Moving the charged impurities 5 $\AA$ from the center of the channel to the substrate increases the mobility by an order of magnitude., Comment: 11 pages, 11 figures

Exciton condensation in an electron-hole bilayer system of monolayer transition metal dichalcogenides is analyzed at three different levels of theory to account for screening and quasiparticle renormalization. The large effective masses of the transition metal dichalcogenides place them in a strong coupling regime. In this regime, mean field (MF) theory with either an unscreened or screened interlayer interaction predicts a room temperature condensate. Interlayer and intralayer interactions renormalize the quasiparticle dispersion, and this effect is included in a GW approximation. The renormalization reverses the trends predicted from the unscreened or screened MF theories. In the strong coupling regime, intralayer interactions have a large impact on the magnitude of the order parameter and its functional dependencies on effective mass and carrier density., Comment: 7 pages, 7 figures

Interlayer electron transport through a graphene / hexagonal boron-nitride (h-BN) / graphene heterostructure is strongly affected by the misorientation angle $\theta$ of the h-BN with respect to the graphene layers with different physical mechanisms governing the transport in different regimes of angle, Fermi level, and bias. The different mechanisms and their resulting signatures in resistance and current are analyzed using two different models, a tight-binding, non-equilibrium Green function model and an effective continuum model, and the qualitative features resulting from the two different models compare well. In the large-angle regime ($\theta > 4^\circ$), the change in the effective h-BN bandgap seen by an electron at the $K$ point of the graphene causes the resistance to monotonically increase with angle by several orders of magnitude reaching a maximum at $\theta = 30^\circ$. It does not affect the peak-to-valley current ratios in devices that exhibit negative differential resistance. In the small-angle regime ($\theta < 4^\circ$), Umklapp processes open up new conductance channels that manifest themselves as non-monotonic features in a plot of resistance versus Fermi level that can serve as experimental signatures of this effect. For small angles and high bias, the Umklapp processes give rise to two new current peaks on either side of the direct tunneling peak.

Interlayer misorientation in transition metal dichalcogenides alters the interlayer distance, the electronic band structure, and the vibrational modes, but, its effect on the interlayer resistance is not known. This work analyzes the coherent interlayer resistance of misoriented 2H-MoS2 for low energy electrons and holes as a function of the misorientation angle. The electronic interlayer resistance monotonically increases with the supercell lattice constant by several orders of magnitude similar to that of misoriented bilayer graphene. The large hole coupling gives low interlayer hole resistance that weakly depends on the misorientation angle. Interlayer rotation between an n-type region and a p-type region will suppress the electron current with little effect on the hole current. We estimate numerical bounds and explain the results in terms of the orbital composition of the bands at high symmetry points. Density functional theory calculations provide the interlayer coupling used in both a tunneling Hamiltonian and a non-equilibrium Green function calculation of the resistivity., Comment: 6 pages, 3 figures

The valence band of a variety of few-layer, two-dimensional materials consists of a ring of states in the Brillouin zone. The energy-momentum relation has the form of a `Mexican hat' or a Rashba dispersion. The two-dimensional density of states is singular at or near the band edge, and the band-edge density of modes turns on nearly abruptly as a step function. The large band-edge density of modes enhances the Seebeck coefficient, the power factor, and the thermoelectric figure of merit ZT. Electronic and thermoelectric properties are determined from ab initio calculations for few-layer III-VI materials GaS, GaSe, InS, InSe, for Bi$_{2}$Se$_{3}$, for monolayer Bi, and for bilayer graphene as a function of vertical field. The effect of interlayer coupling on these properties in few-layer III-VI materials and Bi$_{2}$Se$_{3}$ is described. Analytical models provide insight into the layer dependent trends that are relatively consistent for all of these few-layer materials. Vertically biased bilayer graphene could serve as an experimental test-bed for measuring these effects.

As the only non-carbon elemental layered allotrope, few-layer black phosphorus or phosphorene has emerged as a novel two-dimensional (2D) semiconductor with both high bulk mobility and a band gap. Here we report fabrication and transport measurements of phosphorene-hexagonal BN (hBN) heterostructures with one-dimensional (1D) edge contacts. These transistors are stable in ambient conditions for >300 hours, and display ambipolar behavior, a gate-dependent metal-insulator transition, and mobility up to 4000 $cm^2$/Vs. At low temperatures, we observe gate-tunable Shubnikov de Haas (SdH) magneto-oscillations and Zeeman splitting in magnetic field with an estimated g-factor ~2. The cyclotron mass of few-layer phosphorene holes is determined to increase from 0.25 to 0.31 $m_e$ as the Fermi level moves towards the valence band edge. Our results underscore the potential of few-layer phosphorene (FLP) as both a platform for novel 2D physics and an electronic material for semiconductor applications., Comment: minor correction of typos, equations and references

The electronic and thermoelectric properties of one to four monolayers of MoS$_{2}$, MoSe$_{2}$, WS$_{2}$, and WSe$_{2}$ are calculated. For few layer thicknesses,the near degeneracies of the conduction band $K$ and $\Sigma$ valleys and the valence band $\Gamma$ and $K$ valleys enhance the n-type and p-type thermoelectric performance. The interlayer hybridization and energy level splitting determine how the number of modes within $k_BT$ of a valley minimum changes with layer thickness. In all cases, the maximum ZT coincides with the greatest near-degeneracy within $k_BT$ of the band edge that results in the sharpest turn-on of the density of modes. The thickness at which this maximum occurs is, in general, not a monolayer. The transition from few layers to bulk is discussed. Effective masses, energy gaps, power-factors, and ZT values are tabulated for all materials and layer thicknesses.

A number of the charge-density-wave materials reveal a transition to the macroscopic quantum state around 200 K. We used graphene-like mechanical exfoliation of titanium diselenide crystals to prepare a set of films with different thicknesses. The transition temperature to the charge-density-wave state was determined via modification of Raman spectra of titanium diselenide films. It was established that the transition temperature can increase from its bulk value to ~240 K as the thickness of the van-der-Waals films reduces to the nanometer range. The obtained results are important for the proposed applications of such materials in the collective-state information processing, which require room-temperature operation., Comment: 13 pages