Publications

The valleytronic state found in group-VI transition-metal dichalcogenides such as MoS₂ has attracted immense interest since its valley degree of freedom could be used as an information carrier. However, valleytronic applications require spontaneous valley polarization. Such an electronic state is predicted to be accessible in a new ferroic family of materials, i.e., ferrovalley materials, which features the coexistence of spontaneous spin and valley polarization. Although many atomic monolayer materials with hexagonal lattices have been predicted to be ferrovalley materials, no bulk ferrovalley material candidates have been reported or proposed. Here, we show that a new non-centrosymmetric van der Waals (vdW) semiconductor Cr_0.32Ga_0.68Te_2.33, with intrinsic ferromagnetism, is a possible candidate for bulk ferrovalley material. This material exhibits several remarkable characteristics: (i) it forms a natural heterostructure between vdW gaps, a quasi-two-dimensional (2D) semiconducting Te layer with a honeycomb lattice stacked on the 2D ferromagnetic slab comprised of the (Cr, Ga)-Te layers, and (ii) the 2D Te honeycomb lattice yields a valley-like electronic structure near the Fermi level, which, in combination with inversion symmetry breaking, ferromagnetism, and strong spin-orbit coupling contributed by heavy Te element, creates a possible bulk spin-valley locked electronic state with valley polarization as suggested by our DFT calculations. Further, this material can also be easily exfoliated to 2D atomically thin layers. Therefore, this material offers a unique platform to explore the physics of valleytronic states with spontaneous spin and valley polarization in both bulk and 2D atomic crystals.

Exact density functionals for the exchange and correlation energies are\napproximated in practical calculations for the ground-state electronic\nstructure of a many-electron system. An important exact constraint for the\nconstruction of approximations is to recover the correct non-relativistic\nlarge-$Z$ expansions for the corresponding energies of neutral atoms with\natomic number $Z$ and electron number $N=Z$, which are correct to leading order\n($-0.221 Z^{5/3}$ and $-0.021 Z \\ln Z$ respectively) even in the lowest-rung or\nlocal density approximation. We find that hydrogenic densities lead to\n$E_x(N,Z) \\approx -0.354 N^{2/3} Z$ (as known before only for $Z \\gg N \\gg 1$)\nand $E_c \\approx -0.02 N \\ln N$. These asymptotic estimates are most correct\nfor atomic ions with large $N$ and $Z \\gg N$, but we find that they are\nqualitatively and semi-quantitatively correct even for small $N$ and for $N\n\\approx Z$. The large-$N$ asymptotic behavior of the energy is pre-figured in\nsmall-$N$ atoms and atomic ions, supporting the argument that widely-predictive\napproximate density functionals should be designed to recover the correct\nasymptotics. It is shown that the exact Kohn-Sham correlation energy, when\ncalculated from the pure ground-state wavefunction, should have no contribution\nproportional to $Z$ in the $Z\\to \\infty$ limit for any fixed $N$.\n

The recent development of accurate and efficient semilocal density functionals on the third rung of Jacob's ladder of density functional theory, such as the revised regularized strongly constrained and appropriately normed (r2SCAN) density functional, could enable rapid and highly reliable prediction of the elasticity and temperature dependence of thermophysical parameters of refractory elements and their intermetallic compounds using the quasi-harmonic approximation (QHA). Here, we present a comparative evaluation of equilibrium cell volumes, cohesive energy, mechanical moduli, and thermophysical properties (Debye temperature and thermal expansion coefficient) for 22 transition metals using semilocal density functionals, including the local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE) and PBEsol generalized gradient approximations (GGAs), and the r2SCAN meta-GGA. PBEsol and r2SCAN deliver the same level of accuracies for structural, mechanical, and thermophysical properties. PBE and r2SCAN perform better than LDA and PBEsol for calculating cohesive energies of transition metals. Among the tested density functionals, r2SCAN provides an overall well-balanced performance for reliably computing cell volumes, cohesive energies, mechanical properties, and thermophysical properties of various 3d, 4d, and 5d transition metals using QHA. Therefore, we recommend that r2SCAN could be employed as a workhorse method to evaluate thermophysical properties of transition metal compounds and alloys in high throughput workflows.

Immediately after the demonstration of the high-quality electronic properties in various two dimensional (2D) van der Waals (vdW) crystals fabricated with mechanical exfoliation, many methods have been reported to explore and control large scale fabrications. Comparing with recent advancements in fabricating 2D atomic layered crystals, large scale production of one dimensional (1D) nanowires with thickness approaching molecular or atomic level still remains stagnant. Here, we demonstrate the high yield production of a 1D vdW material, semiconducting Ta2Pd3Se8 nanowires, by means of liquid-phase exfoliation. The thinnest nanowire we have readily achieved is around 1 nm, corresponding to a bundle of one or two molecular ribbons. Transmission electron microscopy and transport measurements reveal the as-fabricated Ta2Pd3Se8 nanowires exhibit unexpected high crystallinity and chemical stability. Our low frequency Raman spectroscopy reveals clear evidence of the existing of weak inter-ribbon bindings. The fabricated nanowire transistors exhibit high switching performance and promising applications for photodetectors.

Birnessite, a layered-structure MnO₂, is an earth-abundant functional material with potential for various energy and environmental applications, such as water oxidation. An important feature of birnessite is the existence of Mn(III) within the MnO₂ layers, accompanied by interlayer charge-neutralizing cations. Using first-principles calculations, we reveal the nature of Mn(III) in birnessite with the concept of the small polaron, a special kind of point defect. Further taking into account the effect of the spatial distribution of Mn(III), we propose a theoretical model to explain the structure-performance dependence of birnessite as an oxygen evolution catalyst. We find an internal potential step which leads to the easy switching of the oxidation state between Mn(III) and Mn(IV) that is critical for enhancing the catalytic activity of birnessite. Finally, we conduct a series of comparative experiments which support our model.

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice I<i>h</i> at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.

The recently proposed rSCAN functional [J. Chem. Phys. 150, 161101 (2019)] is a regularized form of the SCAN functional [Phys. Rev. Lett. 115, 036402 (2015)] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.

We discuss the interplay between magnetic and structural degrees of freedom in elemental Mn. The equilibrium volume is shown to be sensitive to magnetic interactions between the Mn atoms. While the standard generalized gradient approximation underestimates the equilibrium volume, a more accurate treatment of the effects of electronic localization and magnetism is found to solve this longstanding problem. Our calculations also reveal the presence of a magnetic phase in strained $\ensuremath{\alpha}$-Mn that has been reported previously in experiments. This new phase of strained $\ensuremath{\alpha}$-Mn exhibits a noncollinear spin structure with large magnetic moments.