-
Net and Hidden Spin-Valley Locking Enable Ultrahigh Hole Mobility in Covalent Bulk WN$_2$
Authors:
Rong-Tian Pang,
Zhongjuan Han,
Jiayi Gong,
Jiangang He,
Jin-Jian Zhou,
Yugui Yao
Abstract:
High carrier mobility at room temperature underpins high-performance electronics, yet high hole mobility remains rare in bulk semiconductors. Spin-valley locking can suppress intervalley scattering and enhance mobility, but it is limited to materials with broken inversion symmetry. Hidden spin polarization offers a possible route beyond this constraint, although whether its compensated spin textur…
▽ More
High carrier mobility at room temperature underpins high-performance electronics, yet high hole mobility remains rare in bulk semiconductors. Spin-valley locking can suppress intervalley scattering and enhance mobility, but it is limited to materials with broken inversion symmetry. Hidden spin polarization offers a possible route beyond this constraint, although whether its compensated spin textures could protect charge transport remains unclear. Using ab initio electron-phonon and transport calculations, we show that the two hexagonal phases of bulk WN$_2$ realize net and hidden spin-valley locking and exhibit ultrahigh room-temperature hole mobilities. In non-centrosymmetric $α$-WN$_2$, a large valley spin splitting produces net spin-valley locking that nearly eliminates phonon-mediated intervalley scattering. In centrosymmetric $β$-WN$_2$, hidden Zeeman-type spin polarization yields a compensated, sector-resolved spin texture that reverses between valleys and suppresses intervalley scattering as effectively as the net locking does. The stiff W-N/N-N covalent network further keeps the remaining intravalley scattering weak. Our results establish hidden spin polarization as an effective transport-protection mechanism and extend spin-valley engineering to centrosymmetric bulk semiconductors.
△ Less
Submitted 7 August, 2026;
originally announced August 2026.
-
Ultrahigh Intrinsic Hole Mobilities in $M$N$_2$ ($M$= Mo and W) at Room Temperature
Authors:
Zhongjuan Han,
Rong-Tian Pang,
Wu Xiong,
Zhonghao Xia,
Jin-Jian Zhou,
Jiangang He
Abstract:
High-mobility $p$-type semiconductors are essential for advanced electronic devices but remain scarce. Here, using a hierarchical screening framework that combines first-principles calculations with Boltzmann transport theory, we identify $M$N$_2$ ($M$= Mo and W) family as polar semiconductors with exceptionally high intrinsic hole mobilities. In particular, 1H-WN$_2$ exhibits a room-temperature h…
▽ More
High-mobility $p$-type semiconductors are essential for advanced electronic devices but remain scarce. Here, using a hierarchical screening framework that combines first-principles calculations with Boltzmann transport theory, we identify $M$N$_2$ ($M$= Mo and W) family as polar semiconductors with exceptionally high intrinsic hole mobilities. In particular, 1H-WN$_2$ exhibits a room-temperature hole mobility exceeding $10^{4}$~$\mathrm{cm^2\,V^{-1}\,s^{-1}}$. This exceptional transport performance arises from the synergistic suppression of polar-optical-phonon and acoustic-phonon scattering, together with a reduced intervalley-scattering phase space induced by spin--valley locking. These effects arise from anomalously small Born effective charges, strong covalent N--N bonds, and orbital hybridization between N-$2p_x$/$2p_y$ and W-$5d_{xy}$/$5d_{x^2-y^2}$ in the N$_2$-dimer-based structure. Our results establish MoN$_2$ and WN$_2$ as a promising class of high-mobility polar semiconductors and introduce a crystal-structure-based strategy for concurrently suppressing multiple electron--phonon scattering channels, thereby revising design principles for high-mobility materials.
△ Less
Submitted 1 August, 2026;
originally announced August 2026.
-
Optimization and variability can coexist
Authors:
Marianne Bauer,
William Bialek,
Chase Goddard,
Caroline M. Holmes,
Kamesh Krishnamurthy,
Stephanie E. Palmer,
Rich Pang,
David J. Schwab,
Lee Susman
Abstract:
Many biological systems perform close to their physical limits, but promoting this optimality to a general principle seems to require implausibly fine tuning of parameters. Using examples from a wide range of systems, we show that this intuition is wrong. Near an optimum, functional performance depends on parameters in a "sloppy'' way, with some combinations of parameters being only weakly constra…
▽ More
Many biological systems perform close to their physical limits, but promoting this optimality to a general principle seems to require implausibly fine tuning of parameters. Using examples from a wide range of systems, we show that this intuition is wrong. Near an optimum, functional performance depends on parameters in a "sloppy'' way, with some combinations of parameters being only weakly constrained. Absent any other constraints, this predicts that we should observe widely varying parameters, and we make this precise: the entropy in parameter space can be extensive even if performance on average is very close to optimal. This removes a major objection to optimization as a general principle, and rationalizes the observed variability.
△ Less
Submitted 29 May, 2025;
originally announced May 2025.
-
Exact minimax entropy models of large-scale neuronal activity
Authors:
Christopher W. Lynn,
Qiwei Yu,
Rich Pang,
Stephanie E. Palmer,
William Bialek
Abstract:
In the brain, fine-scale correlations combine to produce macroscopic patterns of activity. However, as experiments record from larger and larger populations, we approach a fundamental bottleneck: the number of correlations one would like to include in a model grows larger than the available data. In this undersampled regime, one must focus on a sparse subset of correlations; the optimal choice con…
▽ More
In the brain, fine-scale correlations combine to produce macroscopic patterns of activity. However, as experiments record from larger and larger populations, we approach a fundamental bottleneck: the number of correlations one would like to include in a model grows larger than the available data. In this undersampled regime, one must focus on a sparse subset of correlations; the optimal choice contains the maximum information about patterns of activity or, equivalently, minimizes the entropy of the inferred maximum entropy model. Applying this ``minimax entropy" principle is generally intractable, but here we present an exact and scalable solution for pairwise correlations that combine to form a tree (a network without loops). Applying our method to over one thousand neurons in the mouse hippocampus, we find that the optimal tree of correlations reduces our uncertainty about the population activity by 14% (over 50 times more than a random tree). Despite containing only 0.1% of all pairwise correlations, this minimax entropy model accurately predicts the observed large-scale synchrony in neural activity and becomes even more accurate as the population grows. The inferred Ising model is almost entirely ferromagnetic (with positive interactions) and exhibits signatures of thermodynamic criticality. These results suggest that a sparse backbone of excitatory interactions may play an important role in driving collective neuronal activity.
△ Less
Submitted 18 December, 2023;
originally announced February 2024.
-
Excitonic Instability in Ta2Pd3Te5 Monolayer
Authors:
Jingyu Yao,
Haohao Sheng,
Ruihan Zhang,
Rongtian Pang,
Jin-Jian Zhou,
Quansheng Wu,
Hongming Weng,
Xi Dai,
Zhong Fang,
Zhijun Wang
Abstract:
By systematic theoretical calculations, we have revealed an excitonic insulator (EI) in the Ta2Pd3Te5 monolayer. The bulk Ta2Pd3Te5 is a van der Waals (vdW) layered compound, whereas the vdW layer can be obtained through exfoliation or molecular-beam epitaxy. First-principles calculations show that the monolayer is a nearly zero-gap semiconductor with the modified Becke-Johnson functional. Due to…
▽ More
By systematic theoretical calculations, we have revealed an excitonic insulator (EI) in the Ta2Pd3Te5 monolayer. The bulk Ta2Pd3Te5 is a van der Waals (vdW) layered compound, whereas the vdW layer can be obtained through exfoliation or molecular-beam epitaxy. First-principles calculations show that the monolayer is a nearly zero-gap semiconductor with the modified Becke-Johnson functional. Due to the same symmetry of the band-edge states, the two-dimensional polarization $α_{2D}$ would be finite as the band gap goes to zero, allowing for an EI state in the compound. Using the first-principles many-body perturbation theory, the GW plus Bethe-Salpeter equation calculation reveals that the exciton binding energy is larger than the single-particle band gap, indicating the excitonic instability. The computed phonon spectrum suggests that the monolayer is dynamically stable without lattice distortion. Our findings suggest that the Ta2Pd3Te5 monolayer is an excitonic insulator without structural distortion.
△ Less
Submitted 23 August, 2024; v1 submitted 2 January, 2024;
originally announced January 2024.
-
Exactly solvable statistical physics models for large neuronal populations
Authors:
Christopher W. Lynn,
Qiwei Yu,
Rich Pang,
William Bialek,
Stephanie E. Palmer
Abstract:
Maximum entropy methods provide a principled path connecting measurements of neural activity directly to statistical physics models, and this approach has been successful for populations of $N\sim 100$ neurons. As $N$ increases in new experiments, we enter an undersampled regime where we have to choose which observables should be constrained in the maximum entropy construction. The best choice is…
▽ More
Maximum entropy methods provide a principled path connecting measurements of neural activity directly to statistical physics models, and this approach has been successful for populations of $N\sim 100$ neurons. As $N$ increases in new experiments, we enter an undersampled regime where we have to choose which observables should be constrained in the maximum entropy construction. The best choice is the one that provides the greatest reduction in entropy, defining a "minimax entropy" principle. This principle becomes tractable if we restrict attention to correlations among pairs of neurons that link together into a tree; we can find the best tree efficiently, and the underlying statistical physics models are exactly solved. We use this approach to analyze experiments on $N\sim 1500$ neurons in the mouse hippocampus, and show that the resulting model captures the distribution of synchronous activity in the network.
△ Less
Submitted 16 October, 2023;
originally announced October 2023.
-
Super-Exchange Theory for Polyvalent Anion Magnets
Authors:
Fang Zhang,
You-Chao Kong,
Rui Pang,
Liang Hu,
Peng-Lai Gong,
Xing-Qiang Shi,
Zi-Kang Tang
Abstract:
The Goodenough-Kanamori-Anderson (GKA) rules have been widely applied for explaining the magnetic properties induced by super-exchange interaction. As conclusions of the super-exchange theory, they reveal the antiferromagnetic (ferromagnetic) ordering along with bond angle of 180 degrees (90 degrees) in the cation-anion-cation interaction path, in which the theory sets a pre-condition that the ele…
▽ More
The Goodenough-Kanamori-Anderson (GKA) rules have been widely applied for explaining the magnetic properties induced by super-exchange interaction. As conclusions of the super-exchange theory, they reveal the antiferromagnetic (ferromagnetic) ordering along with bond angle of 180 degrees (90 degrees) in the cation-anion-cation interaction path, in which the theory sets a pre-condition that the electronic states of cations in all paths are identical. We observed that the GKA rules are in fact not universal and even invalid to materials containing anions with different valence states, for example, the layered CrOCl crystal (with two valence states of anions: O2- and Cl-). In this study, we propose an extended super-exchange theory (ESET) related to superposed electronic states of cation in a specific path. ESET is capable of predicting not only the sign and relative magnitude of magnetic exchange constants in different cation-anion-cation paths, but also the magnetic ground state. Through our proposed theory, we conclude that the magnetic ordering along with bond angle of 90 degrees in Cr-Cl-Cr path is moderately antiferromagnetic and of 180 degrees in Cr-O-Cr path is strongly ferromagnetic, which are opposite to the contents of GKA rules. Moreover, we clarify that monolayer CrOCl has antiferromagnetic ordering rather than ferromagnetic as reported recently. The reliability of ESET is verified via first-principles calculation and previous experimental report as well, and its universality is also demonstrated. Thus, our theory is powerful to predict the magnetic properties, which makes it possible to design new high Curie temperature two-dimensional semiconducting ferromagnets with polyvalent anion materials.
△ Less
Submitted 19 February, 2019;
originally announced February 2019.
-
Giant atomic magnetocrystalline anisotropy from degenerate orbitals around Fermi level
Authors:
Rui Pang,
Bei Deng,
Xingqiang Shi
Abstract:
Nano-structures with giant magnetocrystalline anisotropy energies (MAE) are desired in designing miniaturized magnetic storage and quantum computing devices. Through ab initio and model calculations, we propose that special p-element dimers and single-adatom on symmetry-matched substrates possess giant atomic MAE of 72-200 meV with room temperature structural stability. The huge MAE originates fro…
▽ More
Nano-structures with giant magnetocrystalline anisotropy energies (MAE) are desired in designing miniaturized magnetic storage and quantum computing devices. Through ab initio and model calculations, we propose that special p-element dimers and single-adatom on symmetry-matched substrates possess giant atomic MAE of 72-200 meV with room temperature structural stability. The huge MAE originates from degenerate orbitals around Fermi level. More importantly, we developed a simplified quantum mechanical model to understand the principle on how to obtain giant MAE for supported magnetic structures. These discoveries and mechanisms provide a paradigm to design giant atomic MAE in nanostructures.
△ Less
Submitted 27 December, 2016;
originally announced December 2016.
-
Fullerene antiferromagnetic reconstructed spinterface subsurface layer dominates multi-orbitals spin-splitting and large magnetic moment in C60
Authors:
Yangfan Shao,
Rui Pang,
Hui Pan,
Xingqiang Shi
Abstract:
The interfaces between organic molecules and metal surfaces with layered antiferromagnetic order have gained increasing interests in the field of antiferromagnetic spintronics. The C60 layered AFM spinterfaces have been studied for C60 bonded only to the outermost ferromagnetic layer. Using density functional theory calculations, here we demonstrate that C60 adsorption can reconstruct the layered…
▽ More
The interfaces between organic molecules and metal surfaces with layered antiferromagnetic order have gained increasing interests in the field of antiferromagnetic spintronics. The C60 layered AFM spinterfaces have been studied for C60 bonded only to the outermost ferromagnetic layer. Using density functional theory calculations, here we demonstrate that C60 adsorption can reconstruct the layered AFM Cr(001) surface so that C60 bonds to the top two Cr layers with opposite spin direction. Surface reconstruction drastically changes C60 s spintronic properties 1 the spin-split p-d hybridization involve multi-orbitals of C60 and metal double layers, 2 the subsurface layer dominates the C60 spin properties, and 3) reconstruction induces a large magnetic moment in C60 of 0.58 B, which is a synergetic effect of the top two layers as a result of a magnetic direct-exchange interaction. Understanding these complex spinterfaces phenomena is a crucial step for their device applications. The surface reconstruction can be realized by annealing at above room temperature in experiments.
△ Less
Submitted 26 October, 2016;
originally announced October 2016.
-
Manipulating Magnetism at Organic/Ferromagnetic Interfaces by Fullerene-Induced Surface Reconstruction
Authors:
R. Pang,
X. Shi,
M. A. van Hove
Abstract:
Fullerenes have several advantages as potential materials for organic spintronics. Through a theoretical first-principles study, we report that fullerene C$_{60}$ adsorption can induce a magnetic reconstruction in a Ni(111) surface and expose the merits of the reconstructed C$_{60}$/Ni(111) \emph{spinterface} for molecular spintronics applications. Surface reconstruction drastically modifies the m…
▽ More
Fullerenes have several advantages as potential materials for organic spintronics. Through a theoretical first-principles study, we report that fullerene C$_{60}$ adsorption can induce a magnetic reconstruction in a Ni(111) surface and expose the merits of the reconstructed C$_{60}$/Ni(111) \emph{spinterface} for molecular spintronics applications. Surface reconstruction drastically modifies the magnetic properties at both sides of the C$_{60}$/Ni interface. Three outstanding properties of the reconstructed structure are revealed, which originate from reconstruction enhanced spin-split $\mathrmπ$$-$d coupling between C$_{60}$ and Ni(111): 1) the C$_{60}$ spin polarization and conductance around the Fermi level are enhanced simultaneously, which can be important for read-head sensor miniaturization; 2) localized spin-polarized states appear in C$_{60}$ with a spin-filter functionality, and 3) magnetocrystalline anisotropic energy and exchange coupling in the outermost Ni layer are reduced enormously. Surface reconstruction can be realized simply by controlling the annealing temperature in experiments.
△ Less
Submitted 30 July, 2015;
originally announced July 2015.
-
Lithium Intercalation in Graphene/MoS2 Composites: First-Principles Insights
Authors:
Xiji Shao,
Kedong Wang,
Rui Pang,
Xingqiang Shi
Abstract:
As a storage material for Li-ion batteries, graphene/molybdenum disulfide (Gr/MoS2) composites have been intensively studied in experiments. But the relevant theoretical works from first-principles are lacking. In the current work, van-der-Waals-corrected density functional theory calculations are performed to investigate the interaction of Li in Gr/MoS2 composites. Three interesting features are…
▽ More
As a storage material for Li-ion batteries, graphene/molybdenum disulfide (Gr/MoS2) composites have been intensively studied in experiments. But the relevant theoretical works from first-principles are lacking. In the current work, van-der-Waals-corrected density functional theory calculations are performed to investigate the interaction of Li in Gr/MoS2 composites. Three interesting features are revealed for the intercalated Gr/Li(n)/MoS2 composites (n = 1 to 9). One is the reason for large Li storage capacity of Gr/MoS2: due to the binding energies per Li atom increase with the increasing number of intercalated Li atoms. Secondly, the band gap opening of Gr is found, and the band gap is enlarged with the increasing number of intercalated Li atoms, up to 160 meV with nine Li; hence these results suggest an efficient way to tune the band gap of graphene. Thirdly, the Dirac cone of Gr always preserve for different number of ionic bonded Li atoms.
△ Less
Submitted 21 July, 2015;
originally announced July 2015.
-
Molecular Precursors-Induced Surface Reconstruction at Graphene/Pt(111) Interfaces
Authors:
Qian Wang,
Rui Pang,
Xingqiang Shi
Abstract:
Inspired by experimental observations of Pt(111) surfaces reconstruction at the Pt/graphene (Gr) interfaces with ordered vacancy networks in the outermost Pt layer, the mechanism of the surface reconstruction is investigated by van-der-Waals-corrected density functional theory in combination with particle-swarm optimization algorithm and ab initio atomistic thermodynamics. Our global structural se…
▽ More
Inspired by experimental observations of Pt(111) surfaces reconstruction at the Pt/graphene (Gr) interfaces with ordered vacancy networks in the outermost Pt layer, the mechanism of the surface reconstruction is investigated by van-der-Waals-corrected density functional theory in combination with particle-swarm optimization algorithm and ab initio atomistic thermodynamics. Our global structural search finds a more stable reconstructed (Rec) structure than that was reported before. With correction for vacancy formation energy, we demonstrate that the experimental observed surface reconstruction occurred at the earlier stages of graphene formation: 1) reconstruction occurred when C60 adsorption (before decomposition to form graphene) for C60 as a molecular precursor, or 2) reconstruction occurred when there were (partial) hydrogens retain in the adsorbed carbon structures for C2H4 and C60H30 as precursors. The reason can be attributed to that the energy gain, from the strengthened Pt-C bonding for C of C60 or for C with partial H, compensates the energy cost of formation surface vacancies and makes the reconstruction feasible, especially at elevated temperatures. In the Rec structure, two Pt-C covalent bonds are formed per unit cell, which have a great impact on the adsorbed Gr electronic structures.
△ Less
Submitted 13 July, 2015;
originally announced July 2015.
-
Stability of Two-Dimensional Iron-Carbides Suspended across Graphene Pores: First-principles Particle Swarm Optimization
Authors:
Yangfan Shao,
Rui Pang,
Xingqiang Shi
Abstract:
Inspired by recent experimental realizations of two-dimensional (2D) metals and alloys, we theoretically investigate the stability and electronic properties of monolayer (ML) Fe-C compounds and pure Fe. According to our and others theoretical results, the experiment [Science 343, 1228 (2014)] proposed ML pure Fe square-lattices embedded in graphene (Gr) pores are energetically unstable compared to…
▽ More
Inspired by recent experimental realizations of two-dimensional (2D) metals and alloys, we theoretically investigate the stability and electronic properties of monolayer (ML) Fe-C compounds and pure Fe. According to our and others theoretical results, the experiment [Science 343, 1228 (2014)] proposed ML pure Fe square-lattices embedded in graphene (Gr) pores are energetically unstable compared to that of the Fe triangular-lattices in Gr. To solve the above contradiction, we search for the stable structures of ML Fe-C with various Fe to C ratios (as a generalization of ML Fe in Gr) using ab initio particle swarm optimization technique. A Fe1C1 square-lattice embedded in Gr is found. We propose and demonstrate that the square-lattices observed in the experiment were iron-carbides (Fe-C) but not pure Fe from the square-lattice shape, Fe-Fe lattice constant and energetic considerations. Note that the coexistence of C with Fe cannot be excluded from the experiment. More importantly, we find a lowest energy and dynamically stable structure, ML Fe2C2 with Fe atoms form distorted square lattices. High spin polarization around the Fermi level is predicted for different 2D Fe-C structures due to significant orbital hybridization between C and Fe.
△ Less
Submitted 5 September, 2015; v1 submitted 9 July, 2015;
originally announced July 2015.
-
Lattice dynamics in mono- and few-layer sheets of WS2 and WSe2
Authors:
Weijie Zhao,
Zohreh Ghorannevis,
Amara Kiran Kumar,
Jing Ren Pang,
Minglin Toh,
Xin Zhang,
Christian Kloc,
Ping Heng Tan,
Goki Eda
Abstract:
Thickness is one of the fundamental parameters that define the electronic, optical, and thermal properties of two-dimensional (2D) crystals. Phonons in molybdenum disulfide (MoS2) were recently found to exhibit unique thickness dependence due to interplay between short and long range interactions. Here we report Raman spectra of atomically thin sheets of WS2 and WSe2 in the mono- to few-layer thic…
▽ More
Thickness is one of the fundamental parameters that define the electronic, optical, and thermal properties of two-dimensional (2D) crystals. Phonons in molybdenum disulfide (MoS2) were recently found to exhibit unique thickness dependence due to interplay between short and long range interactions. Here we report Raman spectra of atomically thin sheets of WS2 and WSe2 in the mono- to few-layer thickness regime. We show that, similar to the case of MoS2, the characteristic and modes exhibit stiffening and softening with increasing number of layers, respectively, with a small shift of less than 3 cm-1 due to large mass of the atoms. Thickness dependence is also observed in a series of multiphonon bands arising from overtone, combination, and zone edge phonons, whose intensity exhibit significant enhancement in excitonic resonance conditions. Some of these multiphonon peaks are found to be absent only in monolayers. These features provide a unique fingerprint and rapid identification for monolayer flakes.
△ Less
Submitted 4 September, 2013; v1 submitted 3 April, 2013;
originally announced April 2013.