2026-11-10 | | Total: 71
Phase identification is a computational task involving the classification of material structure using the positions of atoms. Using data from Molecular Dynamics (MD) simulations, our task is to identify the material phase of a set of atoms, based on how the simulations predict their positioning relative to each other. By leveraging existing research on mechanistic neural network theory, particularly Neural Tangent Kernels (NTK), we develop a new and efficient method for classifying simulated material samples. Our method involves defining a set of "atom detectors" that are given positions around the origin in cartesian coordinates. The sum of Gaussians centered at each atom is calculated at each atom detector, and is used as the input to a neural network. We present theoretical and experimental evidence to show that this method is highly effective and efficient at the task. Our method obtains better or similar accuracy scores to commonly used methods while having a major computational speedup.
At finite temperature and pressure, the equilibrium structure of a crystal minimizes the Gibbs free energy rather than the static Born-Oppenheimer energy. We present WARP (Wyckoff-constrained Anharmonic Relaxations at Pressure and temperature), a symmetry-constrained, variable-cell relaxation based on the stochastic temperature-dependent effective potential method, including anharmonicity and zero-point motion. Every position update is projected onto Wyckoff positions and every strain onto symmetric strains, so a relaxation can raise the symmetry but never lower it. For aluminum, the thermal expansion follows measurements at cryogenic temperatures, including the quantum plateau below 50 K. For lanthanum hydrides at megabar pressures, equations of state agree with phonon-based and path-integral calculations, and LaH$_{10}$ relaxes from its rhombohedral static minimum to the observed cubic phase, also with a machine-learned potential fine-tuned on 150 supercells. For brucite, Mg(OD)$_2$, lattice parameters, bond distances and thermal expansion agree with neutron diffraction, and the O-D isotope effect is captured.
The epitaxial growth of monolayer graphene by Si sublimation of 6H-SiC(0001) is one of the most effective ways of producing large-scale monocrystalline graphene samples of high quality. However, graphene is left interacting with a reconstructed carbon buffer layer known as zerolayer graphene (ZLG) and the SiC substrate underneath. Beyond the mere creation of Dirac cone replicas, the effect of this interaction on photoemission measurements is poorly understood. We utilized synchrotron radiation for photon-energy-dependent angle-resolved photoemission spectroscopy (ARPES) to investigate the complex final-state interactions of photoelectrons within the heterostructure of monolayer graphene, ZLG, and SiC. We identify multiple resonances in the photoemission intensity of the Dirac cone as a function of photon energy that originate from interlayer interference, confirmed with a multiple scattering model. We also demonstrate for the first time that the relative intensity of the replica bands induced by the $(6\sqrt{3} \times 6\sqrt{3})\text{R}30^\circ$ reconstruction scales with the electron kinetic energy. Finally, we resolved four nondispersing electronic states, which originate from the ZLG. Our findings show that the three-dimensional (3D) structure of SiC and ZLG manifests itself in the photoemission from two-dimensional (2D) epitaxial graphene. By revealing how 3D substrate architecture dictates the photoemission of a 2D layer, this work establishes a crucial framework for investigating the electronic structure of complex heterostructures via ARPES.
Phase diagrams are a constant work in progress as Gibbs-energy models are fitted, revised, and disputed in the light of an endless stream of experimental and computational data. Yet, phase diagrams are still largely communicated as static figures, in much the same way as a century ago. We present GibbScope, a web-based framework for interactive calculation and visualization of binary phase diagrams. Written in Rust and compiled to WebAssembly, the new CALPHAD engine recomputes a complete phase diagram in the browser in a fraction of a second with zero installation on the user or server side. The high speed allows new ways of interacting with phase diagrams and underlying assessments. Validation against well-established CALPHAD software shows agreement to within 0.01 K for invariant temperatures and within 1e-6 for tie-line compositions for most databases. GibbScope offers interactive CALPHAD parameter adjustments, overlays of different assessments and point-by-point misfit to experimental data, which makes the diagram a live, working object. In addition, GibbScope provides convenient access to 2745 thermodynamic databases from the SGTE collection (Hallstedt, 2025) and their live adjustment. For over 1200 of these databases, all phases are modeled in full with the current implementation.
We use machine-learned interatomic potentials combined with higher-order interatomic force constant models to show that finite-range polar correlations in PbTe are driven by long-range harmonic interactions rather than lattice anharmonicity. Our molecular dynamics simulations reveal a spatial polar correlation range exceeding 50 Å at low temperatures and decreasing with increasing temperature. Selective interatomic force constant mapping demonstrates that crystal anharmonicity usually associated with the formation of polar nanoregions suppresses polar correlations. Instead, long-range harmonic interactions extending to 13 Å are necessary to reproduce the spatial coherence of local dipoles. In reciprocal space, these extended harmonic interactions correspond to high curvature of the transverse optical (TO) phonon branch at the $Γ$ point. Comparison with isostructural NaCl connects these extended interactions to local underscreening of predominantly Coulombic forces. Our results provide a key microscopic description of finite-range polar correlations and a framework for understanding local symmetry breaking in centrosymmetric materials.
Predicting magnetic ordering temperatures requires accurate exchange interactions and an appropriate treatment of spin thermodynamics. Here we assess these requirements through complementary quantum Monte Carlo benchmarks. For five insulating antiferromagnets, we compare ferromagnetic-antiferromagnetic energy differences from diffusion Monte Carlo (DMC) and density functional theory with references reconstructed from inelastic neutron scattering (INS) exchanges. DMC underestimates these references but achieves a mean absolute error comparable to the best-performing functionals tested. Separately, we directly test the previously proposed (S+1)/S rescaling of classical Monte Carlo temperatures against quantum stochastic series expansion simulations for 17 compounds using corresponding INS-derived Heisenberg Hamiltonians. Rescaling reduces the mean absolute classical-quantum temperature difference from 81.7 to 8.2 K, with a mean absolute percentage difference of 5.9 %. These results establish a direct quantum benchmark for the rescaling within the tested models and emphasize that improvements in electronic magnetic energy scales should be accompanied by attention to quantum spin thermodynamics.
Van de Waals materials possess great potential as core active materials for novel electronics, optoelectronics, and sensing applications. In particular, transition metal dichalcogenides (TMDs) and other low-dimensional materials, such as hexagonal boron nitride and binary oxides exhibit memristive behavior, making them attractive candidates as building blocks for resistive memories and neuromorphic computing devices. In this work, we investigate two-terminal lateral memristors based on mechanically exfoliated few layer $MoS_2$ flakes. We study the influence of the Schottky barrier height (SBH) as a function of the contact area. We report different SBHs from 60 meV up to 208 meV in the same $MoS_2$ flake and analyze the effect of the SBH in the memristive-like $J-V_{DS}$ curves and the observed non-volatile/volatile resistance switching. Our results highlight the importance of contact engineering to obtain different memristive properties in $MoS_2$-based devices.
In the iron oxide Fe$_3$BO$_6$ weak ferromagnetism stems from competing contributions from two-site Dzyaloshinskii-Moriya interactions (DMI) and from on-site anisotropies at the two $4c$ and $8d$ types of nonequivalent Fe$^{3+}$ sites. As a result, spin reorientation (SR) transition $Γ_2$ $\rightarrow$ $Γ_4$ in Fe$_3$BO$_6$ uniquely comprises 90$^\circ$ spin switching and a 35% change of weak ferromagnetic moment. We explore ultrafast kinetics of such a transition as a proxy for weak ferromagnetism contributions modified by femtosecond laser pulses. Using magneto-optical pump-probe technique we reveal that the SR transition is triggered by an abrupt decrease of the weak ferromagnetic moment in the $Γ_2$ phase, and relate it to the laser-induced modification of the DMI and anisotropy contributions to the spin canting. By tuning the pump pulse phonon energy between different optical absorption ranges, we show that laser-driven impact on spin canting possesses a nonthermal character and is achieved through the charge-transfer transitions $2p\rightarrow3d$ within Fe$^{3+}$O$^{2-}_6$ octahedra rather than by localized $3d\rightarrow3d$ transitions at Fe$^{3+}$. The decreased weak ferromagnetic moment controls the threshold and the latency of the SR transition, as well as a final state. The latter is achieved via the precessional switching to the $Γ_4$ phase activated by the laser-induced heating. Final state is found to be a coexistence of $Γ_4$ and $Γ_2$ phase with modified spin canting persisting even at hundreds of picoseconds.
In this study, we investigate the construction of efficient model potentials for perovskite heterostructures by transferring relevant dynamic quantities from bulk systems and short-period superlattices. Through first-principles calculations, we provide a comprehensive analysis and comparison of the dynamical properties of different perovskite-based superlattice structures, including (PbTiO3 )n /(SrTiO3 )m superlattices, (ATiO3 )n /(AO) 2 superlattices, and An+1 Tin O 3n+1 Ruddlesden-Popper phases, where A can be Sr, Ba, or Pb. We demonstrate that as the thickness of these structures increases, the dynamical properties of the perovskite blocks converge rapidly towards their bulk counterparts. We discuss this convergence in terms of phonon density of states in reciprocal space and inter-atomic force constants in real space. Our findings suggest a straightforward strategy for developing transferable model potentials that accurately capture the dynamics of layered oxide compounds with progressively thicker layers.
We present a direct-space formulation and implementation of quantum density functional theory (DFT) embedding for molecular active subsystems embedded in frozen environments (hosts) of arbitrary dimensionality, from finite molecular hosts to extended systems. The host is computed once without the embedded molecule (the guest) and is represented only by its electron density fitted in an atom-centered auxiliary Gaussian basis together with its nuclear charges with the periodic electrostatics separated into an explicit near-field and a continuous fast-multipole far-field operator. The method allows the removal of host atoms from this representation to create vacancies or other voids without requiring a supercell creation. Analytic nuclear gradients are also implemented, enabling structural relaxation of the active subsystem. We demonstrate the versatility of the approach through applications to adsorbates on periodic surfaces and molecules embedded in bulk solvent environments. For CO on MgO(001), we show that DFT-in-DFT converges rapidly and reproduces the periodic interaction energy within 0.44 kJ/mol while reducing the cost by a factor of 27. Real-time time-dependent embedded DFT (RT-TDDFET) of methylenecyclopropene in periodic water is showcased to produce a ~0.6 eV solvatochromic blue shift, while strong-field propagation of CO2 on hexagonal boron nitride (hBN) yields substrate-induced effects at the cost of a gas-phase propagation. Coupling the embedding scheme to excited-state correlated wave function theory methods such as CC2 yields a solvent shift within 0.01 eV of supermolecular CC2 references, while reducing the computational cost by over 500 times.
The properties of two-dimensional (2D) materials are governed not only by their intrinsic electronic structure but also in particular by their surface chemistry and interactions with the environment. The surfaces of non-van der Waals (non-vdW) 2D materials are intrinsically reactive and may require chemical stabilization under realistic experimental conditions. Here, we perform a systematic first-principles study of water passivation of non-vdW 2D materials and demonstrate that molecular water adsorption is energetically and dynamically favored for a wide range of systems. Rather than dissociating, water molecules frequently remain intact and form adsorbed ordered surface overlayers. These overlayers can adopt a kagome lattice geometry, leading to the emergence of flat electronic bands directly at the valence-band edge. Analysis of the electronic structure reveals a strong intermolecular coupling between neighboring water lone pairs within this kagome topology, while the locked molecular orientations simultaneously produce an electrostatically frustrated surface potential. These findings establish water not only as a stabilizing agent but also as an active element for electronic structure engineering of non-vdW 2D surfaces.
Y6-based homojunction organic photovoltaic devices have recently achieved exceptional power conversion efficiencies, surpassing those of conventional organic photovoltaic materials. Although unusually long exciton diffusion lengths and efficient exciton dissociation have been reported, the fundamental mechanism underlying this outstanding performance remains under debate. Here, we report that structurally distinct aggregates in Y6 thin films exhibit complementary charge-transfer (CT) characters that enable directional energy transfer and enhanced exciton delocalization. Using ultrafast spectroscopy, we identify highly efficient Förster resonance energy transfer (FRET) between Y6 aggregates, with an ensemble-averaged efficiency of 65% occurring on a sub-picosecond timescale (<1 ps). This aggregate-mediated energy-transfer pathway drives long-range energy migration in Y6 thin films and accounts for the remarkable external quantum efficiencies observed in homojunction devices, even in the absence of conventional donor-acceptor bulk heterojunction morphologies and despite substantial fractions of weak-CT components. Here, using Y6 as a model non-fullerene acceptor (NFA) system, we show that structurally distinct aggregates coexisting in NFA thin films possess complementary CT characters, which are expected to enable directional energy transfer and enhanced exciton delocalization. We establish a general multi-component FRET framework that quantitatively links molecular aggregation, CT character, and excited-state dynamics in systems with multiple packing motifs. This work provides insight into the development of next-generation high-performance organic electronic materials and devices.
Atomically narrow magnetic domain walls can arise in ultrathin magnetic films, van der Waals magnets, and systems with reduced effective exchange stiffness due to competing exchange interactions. As their width approaches the atomic lattice spacing, the continuum approximation underlying micromagnetics becomes increasingly unreliable. We investigate the structure and stability of such domain walls in a one-dimensional nearest-neighbor J-K-D spin-lattice model incorporating Heisenberg exchange, uniaxial anisotropy, and Dzyaloshinskii-Moriya interaction (DMI). In the absence of DMI, we demonstrate analytically a continuous transition to an atomically sharp easy-axis wall at the critical exchange-to-anisotropy ratio J/K = 3/2, characterized by the softening of a localized Hessian mode. Atomistic spin-dynamics simulations confirm this transition and reveal that fitting discrete wall profiles with the conventional micromagnetic solution, as commonly done in experimental analyses, underestimates the exchange interaction and overestimates the anisotropy. Finite DMI removes the perfectly collinear wall as a stationary configuration, replacing the transition by a smooth crossover toward a strongly localized but canted wall that approaches the collinear easy-axis configuration only asymptotically in the strong-anisotropy limit. We derive an analytical expression for the DMI-induced canting in this regime and demonstrate its strongly enhanced canting response near the critical ratio. These findings establish the limitations of micromagnetic descriptions of atomically narrow domain walls and suggest opportunities for controlling their structure through anisotropy and DMI.
X-ray free-electron lasers (XFELs) provide the temporal resolution needed to extend total scattering and pair distribution function (PDF) analysis to investigate structural dynamics on femtosecond-to-picosecond timescales. We recently demonstrated that high-quality total scattering data and PDFs can be obtained from individual XFEL pulses; here, we take the next step towards time-resolved studies of structural dynamics in solution. Using the Materials Imaging and Dynamics (MID) instrument at the European XFEL, we obtain total scattering data up to $Q_{\max} = 16.4\,\mathring{\mathrm{A}}^{-1}$ from both conventional capillary samples and, critically, from samples delivered by a liquid jet. Despite the large air and solvent backgrounds inherent to this setup, PDFs suitable for quantitative structural analysis are recovered from liquid dispersions and molecular clusters in solution. We further demonstrate pump-probe total scattering from $\mathrm{CsPbBr}_3$ nanocrystals in a liquid jet, resolving a fluence-dependent photoinduced isotropic lattice strain. Additionally, empirical potential structure refinement (EPSR) of XFEL total scattering data from an aqueous Keggin cluster reveals site-specific hydration environments and enables three-dimensional mapping of the surrounding solvent structure. Together, these advances open the door to tracking not only ultrafast structural changes in solution, but also the evolving solvent environments that govern molecular behaviour and reactivity.
The layered van der Waals semiconductor CrSBr combines intrinsic magnetic order with semiconducting and optical properties, making it a promising material for multifunctional electronic and photonic applications. Understanding how defects and structural transformations affect its properties under experimentally relevant conditions, however, remains challenging for first-principles simulations due to the relevant length and time scales. Here, we develop and benchmark machine-learning interatomic potentials (MLIPs) for the Cr-S-Br system, enabling large-scale atomistic simulations of CrSBr with accuracy approaching that of density functional theory (DFT). The charge-aware variant qNEP achieves near-DFT accuracy in predicting lattice constants, formation enthalpies, vacancy energetics, and the anisotropic pressure response of o-CrSBr. Nudged elastic band calculations reveal that pre-existing S vacancies lower the Cr Frenkel-pair migration barrier by up to 1.1 eV, providing an atomistic mechanism for beam-induced Cr redistribution and structural transformations. Displacement threshold calculations further show that Br is the most susceptible species to electron-beam-induced displacement, identifying the relative ease of removing different atomic species. Together, these results establish a validated computational framework for investigating atomic-scale disorder and structural transformations in CrSBr, providing a basis for defect engineering in materials for future magneto-optoelectronic applications.
Polarization rotation is a fundamental mechanism in ferroelectric materials, yet how such rotation arises microscopically in fluorite ferroelectrics and gives rise to distinct domain walls (DW) remains unresolved. Phonon modes offer a microscopic route to such structural transformations, as geometric orientations can modify their frequencies, dispersion relations, and vibrational eigenstates when their characteristic wavelengths become comparable to relevant physical dimensions. Here, by combining condensed phonon mode decomposition with first-principles calculations, we investigate the origin of 90 polarization rotation in the ferroelectric phase of HfO2. Using the cubic parent phase as reference, we show that polarization rotation originates from multiple symmetry equivalent phonon modes of distinct eigenstates at different wavevectors in the Brillouin zone, whose coupling and cooperative reorientation continuously evolve to generate distinct 90 Pca21 variants and switching pathways. Among the 90 variants, the lowest-energy switching pathway through the Pmn21 exhibits an energy barrier comparable to, but slightly lower than, that of the 180 polarization reversal pathway through the P42/nmc phase, indicating that 90 polarization rotation is an energetically accessible alternative switching route. We further analyze 90 DW and find that frustrated phonon modes emerge at the wall center, where the coexistence and competition of condensed modes promote local reconfiguration of the lattice distortions, contributing to the elevated DW energies. Importantly, experimental observations with STEM imaging provide direct evidence of atomic mode distortions in stable 90 DW, consistent with the theoretically identified low-energy DW. These findings provide the fundamental structural mechanism underlying 90 switching dynamics, competing domain-wall physics, and macroscopic ferroelectric behavior.
Finite-temperature phonons are usually fitted to first-principles forces in randomly displaced supercells, and each force calculation can be expensive. We introduce special optimal thermal cells (SOTCs), in which the displacements are designed rather than drawn at random. Our design, SOTC-D, starts from a harmonic trial model. It gives every supercell phonon mode its quantum-thermal amplitude, chooses the sign of each mode so that the force constants are best determined, and pairs each configuration with its inverse so that cubic forces cancel. Updating the trial leads to the self-consistent harmonic solution. Cells that match displacement correlations instead give good thermodynamics but not better dispersions. For eleven crystals, eight force calculations reproduce the effective dispersions within 0.03-0.31 THz, and SOTC-D is the most accurate of five designs in 9 of them. After one trial update, it beats random sampling in all crystals at 4-16 calculations and two-iteration stochastic sampling (sTDEP) in all at eight (supercell basis). On model potentials, the self-consistent form needs 2-16 times fewer calculations than the stochastic self-consistent harmonic approximation in 13 of 14 cases. With the HSE06 hybrid functional, it is ahead of pooled sTDEP at 4-16 calculations in Si and NaCl, and it gives the thermal expansion of Al within 3% of experiment at 300K. The method has limits: in Al and Si cells of more than 100 atoms, sTDEP is more accurate beyond 4-16 calculations (not so in NaCl up to 16); for dynamically stabilized bcc Zr, the outcome depends on the seed; and above 400 K, the free-energy curvature is not resolved.
Dislocation pattern formation is accompanied by the emergence of lattice rotations and concomitant misorientations across dislocation walls (cell boundaries). Recent advances in mesoscopic modelling of dislocation patterns enable for the first time a systematic simulation study of this phenomenon. We present simulation results for lattice rotation and misorientation patterns in crystals with face-centered cubic and rocksalt lattice structure. Statistical analysis of the orientation changes across cell walls reveals a transition associated with a symmetry-breaking bifurcation in slip system activity: Below the bifurcation, statistical tests indicate the presence of a single population of misoriented walls whereas, above the bifurcation, two sub-populations of 'incidental' and 'geometrically necessary' dislocation walls are present. Our data analytic approach demonstrates that this distinction does not require prior image segmentation and emphasizes the statistical robustness of the distinction between two types of dislocation walls.
Iodine-vacancy migration and association are important factors in the structural stability of halide perovskites. Although neighboring vacancies have been suggested to form energetically favorable pairs, their stability and reconfiguration under finite-temperature lattice deformation remain unclear. Here, we combine machine-learning molecular dynamics (MLMD) and density functional theory (DFT) to investigate neighboring iodine vacancies in tilted CsPbI3 under isotropic and biaxial strain. The MLMD trajectories reveal a nonmonotonic dependence of short-time hopping fractions on strain, with strong isotropic compression favoring point-like separation and tension promoting pair-preserving relay hopping. Notably, biaxial compression produces distinct orientation-dependent responses: vacancy pairs initially aligned with the compressed plane preferentially undergo separation, whereas those perpendicular to the plane predominantly exhibit relay hopping. Complementary static calculations using DFT and the SevenNet-0 potential yield repulsive interactions for all configurations examined. These findings highlight the importance of finite-temperature dynamics and strain-dependent local geometry in understanding vacancy-pair stability and reconfiguration in halide perovskites.
We investigate hydrogen diffusion in tungsten under thermal gradients in the presence of trapping sites using multiscale simulations. Using molecular dynamics simulations, we find that the heat of transport, which quantifies the coupling between temperature gradients and hydrogen transport, for hydrogen in tungsten is positive (thermophobic) at lower temperatures but changes sign near $1000$ K. We investigate the origin of the heat of transport and find that harmonic transition state theory with elastic thermal-expansion corrections is insufficient to explain the observed behaviour, consistent with either a breakdown of local thermodynamic equilibrium or non-negligible anharmonic contributions. We further extend previous studies of hydrogen thermodiffusion by including the effect of trapping sites in both simulations and analytical models. We show that when trapped hydrogen accounts for a significant fraction of the total hydrogen concentration, the heat of transport inferred from measurements of the total concentration is dominated by the temperature dependence of the trapped hydrogen population and no longer reflects the gradient in the mobile hydrogen concentration. Consequently, experimentally measured heats of transport cannot in general be used directly in continuum-scale models that explicitly account for trapping and detrapping. In such models, the heat of transport should be based on mobile hydrogen concentrations only to avoid incorporating trapping effects twice.
Current machine-learned interatomic potentials offer a choice between small, fast models and large, expressive networks. This study demonstrates how an intermediate potential can approach the limit of transferable accuracy at the cost of the cheapest existing formalisms. Specifically, we model He-bearing W-Ta-Hf-Cr-V-Ti alloys with an atomic cluster expansion that learns the radial components of each polynomial as independent linear transforms and embeds numerous densities with learned exponents. While (near-)training loss appears largely a function of model size, the error on a distinct validation set is close to that of leading potentials with ten times more parameters. Taking advantage of the new structure, we develop optimized implementations that run more than an order of magnitude faster than the larger models and many times faster than smaller models with equivalent polynomials. The potential is trained on systematically enumerated small structures but generalizes to point and planar defects and can be extended to even more complex chemistries. As a demonstration, we simulate radiation-damage cascades in billion-atom systems for hundreds of thousands of timesteps.
Chemical separation processes represent one of the largest contributors to industrial energy consumption, underscoring the need for the most efficient separation strategies. Many separation processes rely on selective transport through porous media. In this context, two-dimensional (2D) materials offer a promising route to atomically precise porosity beyond the reach of traditional top-down fabrication, with the potential to enhance separation performance by increasing flow. Here we report measurements of gas transport through a one-molecule-thick layer of suspended, nanoporous poly(carboxyfantrip), a 2D polymer that has identical, nanometer-sized pores that cover nearly a quarter of the surface. This combination of ultimate thinness with high areal pore density gives normalized gas flows (permeances) that approach the effusion limit, three to four orders of magnitude above gas-sieving graphene membranes. We find that gases spanning nearly two orders of magnitude in molecular weight (H$_2$, He, N$_2$, CO$_2$, and SF$_6$) obey Graham's law and approach the ideal limit of non-interacting ballistic transport, despite strong, species-dependent gas-pore interactions. Simulations reproduce the measured permeances within experimental uncertainty and provide an atomistic picture of the transport process. This work demonstrates the potential of synthetic crystalline 2D polymers for chemical separation, water filtration, and quantum sensing.
NiI$_2$ is a van der Waals semiconductor characterized by complex electronic and magnetic phases, holding promise towards functional applications. For this purpose, a detailed understanding of its band structure and its excitations is essential to elucidate the electronic and transport properties. We present a comprehensive description of the band structure and of the lowest valence-to-conduction band transition in bulk NiI$_2$, by means of time- and angle-resolved photoemission spectroscopy supported by first-principles calculations. We describe the orbital character of Ni and I manifolds, and identify the main momentum-resolved features of the valence band and the shallow minima of the conduction band. With this study, we provide a foundational understanding of the electronic structure of NiI$_2$ and its single-particle interband electronic excitation, a benchmark for theoretical models, and a framework for future investigations on this material.
Combining high-resolution wide-angle neutron powder diffraction with magnetic-superspace symmetry analysis, we refine the crystal structure and present, to our knowledge, the first systematic refinement of cubic $P2_13$ MnSi's zero-field magnetic structure across multiple magnetic satellite families. We compare single-arm rhombohedral 3+1D and cubic four-arm 3+4D magnetic superspace-group (MSSG) descriptions with propagation vectors along $\langle111\rangle$. Simultaneous extinction of the two satellite groups around $\{100\}$ dictates a transverse circular-helix model. This suppression and refinements of 18 separately varied satellite-group intensities support this ideal limit with only one adjustable magnetic amplitude. In the ideal single-arm model with MSSG $R3.1'(0,0,g)ts$, the moment has constant magnitude 0.535(3) $μ_B$ per Mn. An explicit four-arm 3+4D cubic superspace construction presented here combines four interpenetrating helices and, for the illustrated phase choice, has moments ranging from zero to 0.874(5) $μ_B$ per Mn, with the same root-mean-square moment as the 3+1D structure. Their identical powder profiles prevent the present measurements from distinguishing them. The average crystal structure remains cubic within experimental resolution, providing structural motivation for the cubic magnetic model. The mean lattice parameter, determined with formal relative precision $σ_a/a \lesssim 10^{-5}$, exhibits nonmonotonic temperature dependence below the Néel temperature.
sdos is an open Python package that computes the density of states (DOS) and related spectral quantities of very large sparse Hermitian matrices A, and of generalized problems A psi = lambda B psi, from Chebyshev moments estimated with random vectors. The density is reconstructed by the kernel polynomial method (KPM, the default) or by maximum entropy (Maxent); the same interface gives band energies, Fermi levels, local and projected densities, harmonic thermodynamics, exact HOMO-LUMO gaps and eigenpairs in a spectral window. Interfaces cover electronic (empirical tight-binding, DFTB, SIESTA) and vibrational (matrix-free LAMMPS Hessians) problems. I checked most ingredients against dense diagonalization or an independent code; the exceptions are listed in Sec. V. On a workstation (4 threads, 19 GB) the DOS of an 884,736-atom amorphous silicon model (a 6^3 tiling of a 4096-atom cell; 3.5 million orbitals; the orthogonal empirical tight-binding model of Kwon et al.) takes 71 s (M=100, R=2, no exact moments; 383 s with the default exact moments), and a matrix-free 2,097,152-atom crystal 10 min in 1.8 GB. A 17^3 replica of a 216-atom amorphous cell (1,061,208 atoms) agrees with exact Bloch spectra: electronic DOS in 229 s, matrix-free ACE vibrational DOS in 2.2 h on one GPU. Shift-invert, incomplete-Cholesky setup, exact moments of high order and Maxent at high moment counts do not scale. The underlying methods are due to Silver, Roder, Kress, Voter, Drabold, Sankey and collaborators; the Introduction states what is added.