No Arabic abstract
Twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure--belonging to symmetry group P3m1--and an out-of-plane intrinsic electric polarization are shown to be metastable by three independendent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2$_1$ (or Pmmn for PbS, PbSe, and PbTe) phase to estimate {em energy barriers} $E_B$ that must be overcome during such structural transformation. Using the curvature of the local minima and $E_B$ as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (nevertheless, and with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe, and can be considered stable for practical purposes only in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature {em ab initio} molecular dynamics for selected compounds. The magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.
We study the injection current response tensor (also known as circular photogalvanic effect or ballistic current) in ferrolectric monolayer GeS, GeSe, SnS, and SnSe. We find that the injection current is perpendicular to the spontaneous in-plane polarization and could reach peak (bulk) values of the order of $10^{10}$A/V$^{2}$s in the visible spectrum. The magnitude of the injection current is the largest reported in the literature to date for a two dimensional material. To rationalize the large injection current, we correlate the injection current spectrum with the joint density of states, electric polarization, strain, etc. We find that various factors such as anisotropy, in-plane polarization and wave function delocalization are important in determining the injection current tensor in these materials. We also find that compression along the polar axis can increase the injection current (or change its sign), and hence strain can be an effective control knob for their nonlinear optical response. Conversely, the injection current can be a sensitive probe of the crystal structure.
We survey the state-of-the-art knowledge of ferroelectric and ferroelastic group-IV monochalcogenide monolayers. These semiconductors feature remarkable structural and mechanical properties, such as a switchable in-plane spontaneous polarization, soft elastic constants, structural degeneracies, and thermally-driven two-dimensional structural transformations. Additionally, these 2D materials also display selective valley excitations, valley Hall effects, and persistent spin helix behavior. After a description of their Raman spectra, a discussion of optical properties arising from their lack of centrosymmetry---such as an unusually strong second-harmonic intensity, large bulk photovoltaic effects, photostriction, and tunable exciton binding energies---is provided as well. The physical properties observed in these materials originate from (correlate with) their intrinsic and switchable electric polarization, and the physical behavior hereby reviewed could be of use in non-volatile memory, valleytronic, spintronic, and optoelectronic devices: these 2D multiferroics enrich and diversify the 2D materials toolbox.
We performed density functional theory calculations with self-consistent van der Waals corrected exchange-correlation (XC) functionals to capture the structure of black phosphorus and twelve monochalcogenide monolayers and find the following results: (a) The in-plane unit cell changes its area in going from the bulk to a monolayer. The change of in-plane distances implies that bonds weaker than covalent or ionic ones are at work within the monolayers themselves. This observation is relevant for the prediction of the critical temperature $T_c$. (b) There is a hierarchy of independent parameters that uniquely define a ground state ferroelectric unit cell (and square and rectangular paraelectric unit cells as well): only 5 optimizable parameters are needed to establish the unit cell vectors and the four basis vectors of the ferroelectric ground state unit cell, while square and rectangular paraelectric structures are defined by only 3 or 2 independent optimizable variables, respectively. (c) The reduced number of independent structural variables correlates with larger elastic energy barriers on a rectangular paraelectric unit cell when compared to the elastic energy barrier of a square paraelectric structure. This implies that $T_c$ obtained on a structure that keeps the lattice parameters fixed (for example, using an NVT ensemble) should be larger than the transition temperature on a structure that is allowed to change in-plane lattice vectors (for example, using the NPT ensemble). (d) The dissociation energy (bulk cleavage energy) of these materials is similar to the energy required to exfoliate graphite and MoS$_2$. (e) There exists a linear relation among the square paraelectric unit cell lattice parameter and the lattice parameters of the rectangular ferroelectric ground state unit cell. These results highlight the subtle atomistic structure of these novel 2D ferroelectrics.
Topological semimetals in ferromagnetic materials have attracted enormous attention due to the potential applications in spintronics. Using the first-principles density functional theory together with an effective lattice model, here we present a new family of topological semimetals with a fully spin-polarized nodal loop in alkaline-metal monochalcogenide $MX$ ($M$ = Li, Na, K, Rb, Cs; $X$ = S, Se, Te) monolayers. The half-metallic ferromagnetism can be established in $MX$ monolayers, in which one nodal loop formed by two crossing bands with the same spin components is found at the Fermi energy. This nodal loop half-metal survives even when considering the spin-orbit coupling owing to the symmetry protection provided by the $mathcal{M}_{z}$ mirror plane. The quantum anomalous Hall state and Weyl-like semimetal in this system can be also achieved by rotating the spin from the out-of-plane to the in-plane direction. The $MX$ monolayers hosting rich topological phases thus offer an excellent materials platform for realizing the advanced spintronics concepts.
We present a three-band tight-binding (TB) model for describing the low-energy physics in monolayers of group-VIB transition metal dichalcogenides $MX_2$ ($M$=Mo, W; $X$=S, Se, Te). As the conduction and valence band edges are predominantly contributed by the $d_{z^{2}}$, $d_{xy}$, and $d_{x^{2}-y^{2}}$ orbitals of $M$ atoms, the TB model is constructed using these three orbitals based on the symmetries of the monolayers. Parameters of the TB model are fitted from the first-principles energy bands for all $MX_2$ monolayers. The TB model involving only the nearest-neighbor $M$-$M$ hoppings is sufficient to capture the band-edge properties in the $pm K$ valleys, including the energy dispersions as well as the Berry curvatures. The TB model involving up to the third-nearest-neighbor $M$-$M$ hoppings can well reproduce the energy bands in the entire Brillouin zone. Spin-orbit coupling in valence bands is well accounted for by including the on-site spin-orbit interactions of $M$ atoms. The conduction band also exhibits a small valley-dependent spin splitting which has an overall sign difference between Mo$X_{2}$ and W$X_{2}$. We discuss the origins of these corrections to the three-band model. The three-band TB model developed here is efficient to account for low-energy physics in $MX_2$ monolayers, and its simplicity can be particularly useful in the study of many-body physics and physics of edge states.