Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: Using the continuum model, we investigate the electronic properties of two types of bilayer graphene (BLG) quantum ring (QR) geometries: (i) An isolated BLG QR and (ii) a monolayer graphene (MLG) with a QR put on top of an infinite graphene sheet (hybrid BLG QR). Solving the Dirac-Weyl equation in the presence of a perpendicular magnetic field and applying the infinite mass boundary condition at the ring boundaries, we obtain analytical results for the energy levels and corresponding wave spinors for both structures. In the case of isolated BLG QR, we observe a sizable and magnetically tunable band gap which agrees with the tight-binding transport simulations. Our analytical results also show the intervalley symmetry EeK (m) = ???EK??? h (m) between the electron (e) and the hole (h) states (m is the angular momentum quantum number) for the energy spectrum of the isolated BLG QR. The presence of interface boundary in a hybrid BLG QR modifies drastically the energy levels as compared with that of an isolated BLG QR. Its energy levels are tunable from MLG dot to isolated BLG QR and to MLG Landau energy levels as the magnetic field is varied. Our predictions can be verified experimentally using different techniques such as by magnetotransport measurements.

Abstract: Using the continuum model, we investigate the confined states and the corresponding wave functions of ABC-stacked trilayer graphene (TLG) quantum dots (QDs). First, a general infinite-mass boundary condition is derived and applied to calculate the electron and hole energy levels of a circular QD in both the absence and presence of a perpendicular magnetic field. Our analytical results for the energy spectra agree with those obtained by using the tight-binding model, where a TLG QD is surrounded by a staggered potential. Our findings show that (i) the energy spectrum exhibits intervalley symmetry E-K(e)(m) = -E-K'(h)(m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number, (ii) the zero-energy Landau level (LL) is formed by the magnetic states with m <= 0 for both Dirac valleys, that is different from monolayer and bilayer graphene QD with infinite-mass potential in which only one of the cones contributes, and (iii) groups of three quantum Hall edge states in the tight-binding magnetic spectrum approach the zero LL, which results from the layer symmetry in TLG QDs.

Abstract: The electron energy spectrum in monolayer-bilayer-monolayer and in bilayer-monolayer-bilayer graphene structures is investigated and the effects of a perpendicular magnetic field and electric bias are studied. Different types of monolayer-bilayer interfaces are considered as zigzag (ZZ) or armchair (AC) junctions which modify considerably the bulk Landau levels (LLs) when the spectra are plotted as a function of the center coordinate of the cyclotron orbit. Far away from the two interfaces, one obtains the well-known LLs for extended monolayer or bilayer graphene. The LL structure changes significantly at the two interfaces or junctions where the valley degeneracy is lifted for both types of junctions, especially when the distance between them is approximately equal to the magnetic length. Varying the nonuniform bias and the width of this junction-to-junction region in either structure strongly influence the resulting spectra. Significant differences exist between ZZ and AC junctions in both structures. The densities of states (DOSs) for unbiased structures are symmetric in energy whereas those for biased structures are asymmetric. An external bias creates interface LLs in the gaps between the LLs of the unbiased system in which the DOS can be quite small. Such a pattern of LLs can be probed by scanning tunneling microscopy.

Abstract: Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalco-genides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory-based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber-type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectralenergy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.

Abstract: We theoretically examine behavior of superconductivity at parallel interfaces separating the domains of another dominant collective excitation, such as charge density waves or spin density waves. Due to their competitive coupling in a two-component Ginzburg-Landau model, suppression of the dominant order parameter at the interfacial planes allows for nucleation of the (hidden) superconducting order parameter at those planes. In such a case, we demonstrate how the number of the parallel interfacial planes and the distance between them are linked to the number and the size of the emerging superconducting gaps in the system, as well as the versatility and temperature evolution of the possible superconducting phases. These findings bear relevance to a broad selection of known layered superconducting materials, as well as to further design of artificial (e.g., oxide) superlattices, where the interplay between competing order parameters paves the way towards otherwise unattainable superconducting states, some with enhanced superconducting critical temperature.

Abstract: The Dzyaloshinskii-Moriya interaction (DMI) is a chiral interaction that favors formation of domain walls. Recent experiments and ab initio calculations show that there are multiple ways to modify the strength of the interfacially induced DMI in thin ferromagnetic films with perpendicular magnetic anisotropy. In this paper we reveal theoretically the effects of spatially varied DMI on the magnetic state in thin films. In such heterochiral 2D structures we report several emergent phenomena, ranging from the equilibrium spin canting at the interface between regions with different DMI, over particularly strong confinement of domain walls and skyrmions within high-DMI tracks, to advanced applications such as domain tailoring nearly at will, design of magnonic waveguides, and much improved skyrmion racetrack memory.

Abstract: The charge density is among the most fundamental solid state properties determining bonding, electrical characteristics, and adsorption or catalysis at surfaces. While atomic-scale charge densities have as yet been retrieved by solid state theory, we demonstrate both charge density and electric field mapping across a mono-/bilayer boundary in 2D MoS2 by momentum-resolved scanning transmission electron microscopy. Based on consistency of the four-dimensional experimental data, statistical parameter estimation and dynamical electron scattering simulations using strain-relaxed supercells, we are able to identify an AA-type bilayer stacking and charge depletion at the Mo-terminated layer edge.

Abstract: Starting from the simplified linear combination of atomic orbitals method in combination with first-principles calculations, we construct a tight-binding (TB) model in the two-centre approximation for borophene and hydrogenated borophene (borophane). The Slater and Koster approach is applied to calculate the TB Hamiltonian of these systems. We obtain expressions for the Hamiltonian and overlap matrix elements between different orbitals for the different atoms and present the SK coefficients in a nonorthogonal basis set. An anisotropic Dirac cone is found in the band structure of borophane. We derive a Dirac low-energy Hamiltonian and compare the Fermi velocities with that of graphene.

Abstract: Using the tight-binding approach coupled with mean-field Hubbard model, we theoretically study the effect of mechanical deformations on the magnetic properties of bilayer graphene (BLG) quantum dots (QDs). Results are obtained for AA-and AB(Bernal)-stacked BLG QDs, considering different geometries (hexagonal, triangular and square shapes) and edge types (armchair and zigzag edges). In the absence of strain, our results show that (i) the magnetization is affected by taking different dot sizes only for hexagonal BLG QDs with zigzag edges, exhibiting different critical Hubbard interactions, and (ii) the magnetization does not depend on the interlayer hopping energies, except for the geometries with zigzag edges and AA stacking. In the presence of in-plane and uniaxial strain, for all geometries we obtain two different magnetization regimes depending on the applied strain amplitude. The appearance of such different regimes is due to the breaking of layer and sublattice symmetries in BLG QDs.

Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: We model the fluorescence spectra of planar spin oscillators to find conditions that maximize spin resonance fluorescence. Spin oscillators perform Rabi oscillations under the effect of a periodic effective magnetic field caused by the winding motion of an electron in a gradient of magnetic field. We show that, despite the weak coupling of the spin magnetic dipole to the vacuum, spin oscillators excited by a direct current output a few nanowatts of microwave power, which is comparable to the best microwave sources. The large quantum efficiency relies on the combination of two effects. On the one hand, the spontaneous emission rate is enhanced by the synchronization of spin oscillators, which interact through the microwave field that they emit. On the other hand, the huge Rabi frequencies experienced by spin oscillators promote spins into upper levels of Zeeman transitions, from which a radiative cascade is triggered. We demonstrate different regimes of fluorescence which correspond to different values of the Rabi period relative to the spontaneous decay time and to the oscillator dwell time in the gradient of magnetic field. We investigate the device parameters which make these regimes experimentally accessible and find conditions that optimize microwave output. We find that microwave emission is centered around the cutoff frequency of spin oscillators. This has the advantage that the peak emission frequency may be tuned from zero continuously up to a few hundred gigahertz using an electrostatic gate. Quite remarkably for a spintronics effect, electrically induced spin resonance fluorescence does not require the injection of a spin polarized current. In fact, we show that microwave spectra are mostly independent of the incoming spin polarization except for magnetic waveguides which are shorter than a certain critical length, which we will specify.

Abstract: In this work we compute the torque field present in a ferromagnet in contact with a metallic nanowire when a skyrmion is present. If the nanowire is narrow enough, then the current is carried by a single conduction band. In this regime the classical torque model breaks down and we show that a skyrmion driven by spin transfer torque moves in a different direction than predicted by the classical model. However, the amount of charge current required to move a skyrmion with a certain velocity in the single-band regime is similar to a classical model of torque where it is implicitly assumed current transport by many conduction bands. The single-band regime is more efficient creating spin current from charge current because of the perfect polarization of the single band but is less efficient creating torque from spin current. Nevertheless, it is possible to take profit of the single-band regime to move skyrmions even with no net charge or spin current flowing between the device contacts. We have also been able to recover the classical limit considering an ensemble of only a few electronic states. In this limit we have discovered that electron diffusion needs to be considered even in ballistic nanowires due the effect of the skyrmion structure on the electron current.

Abstract: We examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.