2026-02-11 | | Total: 34
Vibrational spectroscopy in the electron microscope can reveal phonon excitations with nanometer spatial resolution, yet routine prediction remains out of reach due to fragmented workflows requiring specialized expertise. Here we introduce PySlice, the first publicly available implementation of the Time Autocorrelation of Auxiliary Wavefunction (TACAW) method, providing an automated framework that produces momentum- and energy-resolved vibrational electron energy-loss spectra directly from atomic structures. By integrating universal machine learning interatomic potentials with TACAW, PySlice eliminates the bottleneck of per-system potential development. Users input atomic structures and obtain phonon dispersions, spectral diffraction patterns, and spectrum images through a unified workflow spanning molecular dynamics, GPU-accelerated electron scattering, and frequency-domain analysis. We outline the formulation behind the code, demonstrate its application to canonical systems in materials science, and discuss its use for advanced analysis and materials exploration. The modular Python architecture additionally supports conventional electron microscopy simulations, providing a general-purpose platform for imaging and diffraction calculations. PySlice makes vibrational spectroscopy prediction routine rather than specialized, enabling computational screening for experimental design, systematic exploration of phonon physics across materials families, and high-throughput generation of simulated data for training of future machine learning models.
We present a multiscale simulation framework to compute the current vs. voltage (I-V ) characteristics of metal/oxide/metal structures building the core of conductive bridging random access memory (CBRAM) cells and to shed light on their resistance switching properties. The approach relies on a finite element model whose input material parameters are extracted either from ab initio or from machine-learned empirical calculations. The applied techniques range from molecular dynamics and nudged elastic band to electronic and thermal quantum transport. Such an approach drastically reduces the number of fitting parameters needed and makes the resulting modeling environment more accurate than traditional ones. The developed computational framework is then applied to the investigation of an Ag/a-SiO2/Pt CBRAM, reproducing experimental data very well. Moreover, the relevance of Joule heating is assessed by considering various cell geometries. It is found that self-heating manifests itself in devices with thin conductive filaments with few-nanometer diameters and at current concentrations in the 10s-microampere range. With the proposed methodology it is now possible to explore the potential of not-yet fabricated memory cells and to reliably optimize their design.
Chemical substitution in crystalline quantum materials is a powerful way to explore the consequences of strong spin-orbit coupling on their structural and electronic properties. In this work, we present an investigation of thin films of the La$\textit{Pn}_2$ ($\textit{Pn}$~=~Sb, Bi) class of layered square-net intermetallics. We report the growth of LaBi$_2$ with a pristine layer-by-layer growth mode, classifying it as a good metal displaying superconductivity at $\sim$0.55~K. Compared to LaSb$_2$, we attribute the enhanced metallic behavior and improved growth dynamics of LaBi$_2$ to significant relativistic corrections to its electronic band structure and the resulting impact on both surface energy and intrinsic phonon scattering.
Moiré superlattices in two-dimensional (2D) materials exhibit rich quantum phenomena, but ab initio modelling of these systems remains computationally prohibitive. Existing machine learning methods for accelerating density-functional theory (DFT) can target the prediction of different quantities and often rely on the locality assumption. Here we train a Gaussian process regression SALTED model exclusively on the electron densities of small displaced bilayer structures and then extrapolate electron density prediction to the large supercells required to describe small twist angles between these bilayers. We show the necessity of long-range descriptors to yield reliable band structures and electrostatic properties of large twisted bilayer structures, when these are derived from predicted densities. We demonstrate that the choice of descriptor determines the distribution of residual density errors, which in turn affects the downstream electronic properties. We apply our models to twisted bilayer graphene, hexagonal boron nitride, and transition metal dichalcogenides, focusing on the model's capacity to predict complex phenomena, including flat band formation, bandwidth narrowing, domain-wall electric fields, and spin-orbit coupling effects. Beyond moiré materials, this approach provides a general methodology for electronic structure prediction in large-scale systems with substantial long-range phenomena related to non-local geometric information.
${\rm{EuTe}}_{4}$ is a polar charge density wave (CDW) material, with giant thermal hysteresis and non-volatile state switching under electric and optical fields, attracting great attention in recent years. However, the in-depth understanding of these anomalous phenomena remains elusive. Herein, via first-principles calculations, we reveal that the polar CDW state in ${\rm{EuTe}}_{4}$ hosts a novel layer-selective nature, wherein multiple energetically close CDW configurations coexist and exhibit low interconversion energy barriers. Monte Carlo simulations indicate that the giant thermal hysteresis in ${\rm{EuTe}}_{4}$ originates from a phase transition mainly driven by the change of configurational entropy, around which the material hosts a metastable CDW state characterized by diverse local polar configurations breaking the out-of-plane translational symmetry. The configurational composition of this metastable CDW state can be effectively controlled by electric and optical fields, thereby enabling non-volatile state switching. Our theoretical findings align well with recent experimental observations in ${\rm{EuTe}}_{4}$ and pave the way for exploring the emerging phenomena and applications of polar CDW in multilayered systems.
We compute exciton binding energies and fine-structure splittings in CdSe nanoplatelets with two zincblende geometries and one wurtzite geometry, finding that the wurtzite structure exhibits the largest bright-bright splitting due to its intrinsic in-plane anisotropy, while the zincblende structures show smaller but finite splittings arising from atomistic symmetry breaking at edges and corners. These results are obtained using a theoretical framework that we developed, which combines DFT single-particle states with screened configuration interaction, a quasi-2D dielectric screening model, and an efficient Coulomb-cutoff scheme that eliminates periodic-image interactions and enables accurate Coulomb and exchange integrals at low computational cost. This methodology provides a transferable and practical route for studying excitons in CdSe nanoplatelets and other quasi-two-dimensional nanomaterials.
We present the thermopower of EuCd2P2, a material which exhibits a large resistivity peak with significant magnetic field dependence in the temperature range of 10-25 K. In the same region we observe a highly unusual behavior of the thermopower with two sign changes and giant extrema. The overall variation of the thermopower exceeds 4 000 muV/K and takes place in an extremely narrow temperature region of less than 5 K. The anomaly is suppressed completely in a small magnetic field of 0.5 T. We discuss this observation using a simple drift-diffusion picture and taking into account that the temperature gradient inducing the thermopower voltage is accompanied by a gradient of the electrical resistivity. Our simple estimation yields the correct magnitude, shape, and field dependence of the thermopower anomaly observed in EuCd2P2. These results open a new route to giant thermopower values via gradients of electronic properties.
Twistronics is an emerging and captivating field in condensed matter physics and material science. However, accurately and efficiently calculating the electronic structures of twisted systems remains a significant challenge. To address this, we have developed MoireStudio, a universal Python-based computational package for twisted electronic structures. Its functionalities include commensurate structure search, structure generation, parameterization, and construction for tight-binding models and continuum models, and the precise incorporation of full relaxation effects. The package is applicable to arbitrary combinations of two-dimensional materials, including rectangular lattices and heterostructures. User-friendly and easy to use, MoireStudio supports parallel large-scale computations, provides visualization capabilities, and offers interfaces with third-party software. It is poised to become a convenient and powerful tool for researchers in twistronics fields.
We investigate the fracture behavior of pre-cracked triangular beam-lattices whose elements have failure stresses drawn from a Weibull distribution. Through a statistical analysis and numerical simulations, we identify and verify the existence of three distinct failure regimes: (i) disorder is effectively suppressed, (ii) disorder manifests locally near the crack tip, modifying the crack morphology, and (iii) disorder manifests globally, leading to initially diffuse failure. Our model naturally reveals the key parameters governing this behavior: the Weibull modulus, quantifying the spread in failure thresholds, and a geometric quantity termed the Slenderness Ratio. We also reproduce the disorder-induced toughening reported in previous experimental and numerical studies, further demonstrating that its manifestation depends non-monotonically on disorder. Crucially, our results indicate that this toughening cannot be simply connected to the amount of damage in the lattice, challenging interpretations that attribute increased fracture energy solely to enhanced crack tortuosity or diffuse failure. Overall, our results establish geometry as a powerful control parameter for regulating how disorder is expressed during fracture in beam-lattices, with broader implications for the disorder-induced toughening in engineered materials.
Tetragonal Ta2Se is a layered, Ta-rich chalcogenide that departs from conventional MX2 transition-metal dichalcogenides by hosting dense Ta-Ta networks capped by Se square-net layers. Here, we present a unified first-principles investigation of hydrostatic-pressure tuning in Ta2Se from 0 to 10 GPa, connecting the structural response, mechanical stability, thermophysical indicators, bonding evolution, electronic and optical behavior, lattice dynamics and superconductivity within a single framework. The derived thermophysical descriptors corroborate a pressure-stiffened lattice: density increases, Debye temperature rises, melting temperature is elevated and minimum thermal conductivity increases, whereas the Grüneisen parameter remains within a narrow window, suggesting no anomalous anharmonic softening. Bond population metrics and electron-density-difference analysis revealed a mixed metallic-covalent bonding picture dominated by a robust Ta-Ta metallic backbone, accompanied by pressure-strengthened Ta-Se hybridization. Electronic-structure calculations show persistent metallicity under compression; pressure broadens bands, reduces density of states at the Fermi level, reshapes the Fermi surface and points to a possible Lifshitz-type reconstruction without symmetry breaking. The optical response remained metallic with Drude-like low-energy behavior and pressure-tunable spectral features. The phonon dispersions exhibit no imaginary modes, confirming dynamical stability. Electron-phonon coupling calculations classify Ta2Se as a weak-coupling, phonon-mediated superconductor with Tc around 3.9 K, consistent with available experiments and establish pressure as a practical control knob for stability and superconductivity-relevant descriptors in this metal-rich layered platform.
As one of the most abundant interstitial elements, nitrogen (N) is effective in improving yield strength of metallic materials, due to interstitial solid solution strengthening. Doping N can substantially enhance the yield strength but often leads to a decreased ductility, revealing a strength-ductility trade-off phenomenon. Here, we simultaneously enhance the strength and ductility in a non-equiatomic CrMnFeCoNi high-entropy alloy via N alloying and unravel the underlying microscopic mechanisms. The N-doped alloy (1 at.% N) shows an excellent combination of higher yield strength (104% increase) and larger ductility (38% increase), with a two-stage strain hardening behavior, compared to the N-free alloy. Detailed transmission electron microscopy (TEM) analysis reveals that N-doping introduces short-range order (SRO) domains within the microstructure, leads to pronounced planar slip, and promotes the formation of nano-spaced (6-15 nm) stacking faults and deformation twins. Continuous generation and interaction of the fine-spaced SFs act as a strong barrier for dislocation movement and provide ample room for dislocation storage. The interaction of SRO with dislocations and the evolution of SFs ascribe to the first strain hardening stage, and the disordering of the SRO along with the activation of deformation twins are attributed to the second strain hardening stage. Our work shows that N-doping is effective in simultaneously improving the strength-ductility synergy and provides novel insights into alloy design with slightly elevating the SFE, and manipulating the ordered structure within the HEA.
Germanium diselenide (GeSe$_{2}$) has recently attracted substantial interest as a rare example of one-dimensional (1D) van der Waals material. Here, we investigate the photovoltaic potential of bulk-stacked GeSe$_{2}$ chains using first-principles calculations within the $GW0$ approximation and the Bethe-Salpeter equation (BSE) to capture quasiparticle and excitonic effects. The bulk GeSe$_{2}$ exhibits indirect GW band gaps of 1.92 eV (type-I) and 1.08 eV (type-II). Optical calculations show markedly stronger visible-light absorption in type-II, yielding a spectroscopically limited maximum efficiency (SLME) of ~25.6% at a 0.5 $μ$m thickness. Phonon and room-temperature ab initio molecular dynamics analyses indicate that type-II is dynamically stable, whereas type-I shows imaginary phonon modes, suggesting a propensity for structural distortion. These results identify type-II GeSe2 as a promising stable absorber for thin-film photovoltaics with enhanced flexibility compared to typical 2D vdW systems.
Controlling materials through their interactions with electromagnetic vacuum fluctuations is an emergent frontier in material engineering. Although recent experiments have demonstrated dark cavity effects for electronic material phases, like superconductivity, ferroelectricity and charge density waves, a smoking gun experiment for magnetic systems is lacking. Largely, this comes from the focus on phase transitions, where a large critical light-matter coupling is needed to observe cavity modifications. Here, we propose spiral magnets, where even a small cavity-mediated change in magnetic interactions is reflected in a change of the spiral wavelength, as a promising platform to observe cavity effects. We focus on the single-layer multiferroic NiI$_2$, interacting with electric field fluctuations from surface phonon polaritons of the paraelectric substrate SrTiO$_3$. With decreasing substrate-material distance, the ratio of nearest and third nearest neighbor exchange interactions reduces, leading to an increase of the spiral wavelength and an eventual transition into a ferromagnetic state. Our work identifies a realistic platform to observe cavity vacuum renormalization effects in magnetic systems.
Chirality governs phenomena ranging from chemical reactions to the topology of quasiparticle charge carriers. However, a direct macroscopic probe for crystal chirality remains a significant challenge, especially in time reversal symmetric systems with weak circular dichroism signal. Here, we propose the ``Chiralometer'', a mechanical detection method that probes chirality by driving angular momentum carriers out of equilibrium. Using first-principles calculations and semiclassical transport theory, we demonstrate that a temperature gradient in insulators or an electric field in metals induces uncompensated angular momentum in phonons and electrons, respectively. This imbalance generates a macroscopic mechanical torque ($τ\sim 10^{-11} N \cdot m$) well within the sensitivity of modern torque magnetometry and cantilever-based sensors. We identify robust signatures in chiral crystals such as Te, SiO$_2$, and the topological semimetal CoSi. Our work establishes mechanical torque as a fundamental order parameter for chirality, offering a transformative tool for orbitronics and chiral quantum materials.
Non-collinear Kagome antiferromagnets (AFMs) Mn3X (X = Sn, Ga, Ge, Ir, Pt) can generate an anomalous Hall effect (AHE) despite vanishing net magnetization, enabled by broken time-reversal and inversion symmetries. However, strong in-plane anisotropy has limited studies of the AFM-AHE and electronic applications to coplanar spin configurations. Non-coplanar spin textures in these systems have been realized only in low temperature spin-glass states or at interfaces with heavy metals. Here, we report an intrinsic non-coplanar spin configuration persisting up to 400 K in cubic-phase Mn3Ge, originating from coexisting symmetric and antisymmetric exchange interactions. Competing magnetic states associated with this non-coplanar spin configuration give rise to an unconventional AHE with a magnetic-field-induced sign reversal and a hump-like feature. Our findings establish a platform for non-coplanar magnetism in AFM spintronics.
Following earlier theoretical prediction, intercalated V2Se2O-family altermagnets such as RbV2Te2O and KV2Se2O have now been experimentally confirmed as d-wave altermagnets, representing the only known van der Waals layered altermagnetic systems. By combining crystal-symmetry-paired spin-momentum locking (CSML) with the layered structure, the V2Se2O-family offers a suitable platform for studying low-dimensional spintronic responses and exploring the interplay among multiple quantum degrees of freedom. To establish a concrete theoretical foundation for understanding and utilizing these materials, we investigate six representative members of the V2Se2O-family and construct a realistic tight-binding model parameterized by first-principles calculations, which is benchmarked by experimental measurements. This model accurately captures essential altermagnetic electronic properties, including CSML and noncollinear spin-conserved currents. It further incorporates strain-coupling parameters, enabling the simulation of strain-tunable responses such as the piezo-Hall effects. This realistic model allows systematic exploration of multiple degrees of freedom (like spin, valley, and layer) within a single system, and lays the groundwork for understanding their coupling with other quantum materials, such as topological insulators and superconductors, thereby advancing both the fundamental understanding and potential device applications of this novel class of layered altermagnets.
Magnetic topological materials such as Weyl and Dirac magnets exhibit unconventional electronic properties arising from the interplay between magnetic order and band topology, leading to remarkable thermomagnetic and thermoelectric effects. Here, we investigate the ANE in epitaxial thin films of the Weyl ferromagnet candidate Fe5Si3. A pronounced transverse Nernst response exceeding approximately 1.50 microvolt per kelvin is observed at room temperature, together with a giant anomalous Nernst angle of about 0.56, indicating highly efficient conversion between thermal gradients and transverse electric fields. Beyond the anomalous contribution, a sizable topological Nernst signal of approximately 0.43 microvolt per kelvin persists above room temperature, suggesting the possible presence of real-space Berry curvature associated with nontrivial spin textures. First-principles density functional theory calculations combined with symmetry analysis reveal an unconventional electronic structure in which Weyl nodal lines, nodal surfaces, and nearly flat bands coexist near the Fermi level. This rare concurrence of multiple topological band features produces a strongly enhanced and sharply energy-dependent Berry curvature, which governs both the magnitude and temperature evolution of the observed Nernst response. The close quantitative agreement between calculated anomalous Nernst conductivity and experimental results establishes the topological electronic structure as the dominant origin of the observed thermomagnetic transport, highlighting Fe5Si3 as a chemically simple, low-cost binary topological magnet for exploring both real-space and momentum-space Berry-curvature-driven thermoelectric phenomena.
Altermagnets have recently emerged as a new class of magnetic materials that combine compensated magnetic order with spin-split electronic band structures. In this work, we employ molecular beam epitaxy to grow MnTe thin films with controlled thickness on InP(111)A substrates. By performing angle-resolved photoemission spectroscopy measurements, we observe a large spin splitting of ~230 meV for bulk bands well below the Fermi level and identify surface states that cross the Fermi level. Electrical transport measurements reveal that a robust anomalous Hall (AH) effect persists down to 2 K and an AH sign reversal occurs near 175 K. By systematically tuning film thickness, growth conditions, and interfacial structure, we demonstrate that the AH response in MnTe films originates from the Berry curvature of surface states rather than from bulk bands. Our first-principles calculations reveal that this surface-state-driven AH effect is imprinted by the bulk altermagnetic order and remains unchanged for terminations with opposite Mn magnetic orientations. Our results establish a unique surface transport probe of bulk altermagnetism, demonstrate interface engineering as an effective route to generate and control the AH effect in altermagnets, and provide a unified understanding of the AH response in altermagentic MnTe films.
The quantum size effect (QSE) and strain effect are two key factors influencing the surface morphology of thin films, which can increase film surface roughness through QSE-induced thickness oscillation and strain-induced island formation, respectively. Surface roughness usually manifests in the early stages of film growth and diminishes beyond a critical thickness. In this work, we employ molecular beam epitaxy (MBE) to grow Sn(100) islands with varying thickness N on bilayer graphene-terminated 6H-SiC(0001) substrates. Scanning tunneling microscopy and spectroscopy measurements reveal an inverse surface roughness effect that highlights the interplay of QSE and misfit strain in shaping the surface morphology of Sn(100) islands. For N =< 10, the islands exhibit flat surfaces, while for N >= 26, the island surfaces become corrugated and patterned. For the intermediate range, i.e., 12 =< N =<24, both flat and patterned surfaces coexist, with the percentage coverage of the patterned surface oscillating as a function of N. By performing density functional theory calculations, we demonstrate that the unusual surface pattern evolution in our MBE-grown Sn(100) islands is a result of the interplay between QSE-induced surface roughing and tensile strain-induced smoothening effect.
Every generative model for crystalline materials harbors a critical structure size beyond which its outputs quietly become unreliable -- we call this the extrapolation frontier. Despite its direct consequences for nanomaterial design, this frontier has never been systematically measured. We introduce RADII, a radius-resolved benchmark of ${\sim}$75,000 nanoparticle structures (55-11,298 atoms) that treats radius as a continuous scaling knob to trace generation quality from in-distribution to out-of-distribution regimes under leakage-free splits. RADII provides frontier-specific diagnostics: per-radius error profiles pinpoint each architecture's scaling ceiling, surface-interior decomposition tests whether failures originate at boundaries or in bulk, and cross-metric failure sequencing reveals which aspect of structural fidelity breaks first. Benchmarking five state-of-the-art architectures, we find that: (i) all models degrade by ${\sim}13\%$ in global positional error beyond training radii, yet local bond fidelity diverges wildly across architectures -- from near-zero to over $2\times$ collapse; (ii) no two architectures share the same failure sequence, revealing the frontier as a multi-dimensional surface shaped by model family; and (iii) well-behaved models obey a power-law scaling exponent $α\approx 1/3$ whose in-distribution fit accurately predicts out-of-distribution error, making their frontiers quantitatively forecastable. These findings establish output scale as a first-class evaluation axis for geometric generative models. The dataset and code are available at https://github.com/KurbanIntelligenceLab/RADII.
The coupling between electronic excitations and lattice deformation in van der Waals ferroelectrics is governed by a competition between the electron deformation potential and the inverse piezoelectric effect. While theory predicts that piezoelectric screening should drive a polar-axis contraction in monolayer group-IV monochalcogenides, we demonstrate that in multilayer SnS, the deformation potential provides the dominant contribution, driving a polar-axis expansion even within ferroelectric domains. By correlating polarization-resolved second-harmonic generation microscopy with ultrafast reflectance spectroscopy and first-principles calculations, we resolve the anisotropic lattice response and disentangle intrinsic photostrictive strain from extrinsic thin-film interference artifacts. These results establish a microscopic hierarchy of photostrictive mechanisms and position stacking-engineered SnS as a platform for ultrafast optomechanical transduction.
In contrast to hafnium oxide (HfO2) thin films, where the appearance of switchable ferroelectric polarization can be induced by strain engineering, reliable methods of ferroelectricity control are absent in HfO2 nanoparticles. Direct experimental observations of ferroelectric hysteresis and/or ferroelectric domains, and appropriate modelling of stress-induced ferroelectric states in the nanoparticles are absent too. In this work we study the influence of chemical stresses on phase diagrams and polar properties of spherical HfO2 core-shell nanoparticles using the Landau-Ginzburg-Devonshire free energy functional with higher powers, trilinear and biquadratic couplings of polar, antipolar and nonpolar order parameters. It appeared that the ferroelectric phase can be reentrant with respect to the size of nanoparticles, because the spontaneous polarization exists in the limited range of core radii R_c, namely R_cr^min<R_c<R_cr^max. The minimal critical radius R_cr^min is mainly determined by the size dependence of the depolarization field and correlation effects. The maximal critical radius R_cr^max is mainly determined by the size dependence of chemical stresses, which are induced by the elastic defects in the shell. Analytical expressions derived for the critical radii can be generalized for hafnium-zirconium oxide nanoparticles, providing that corresponding parameters of the free energy are known from the first principles calculations.
Understanding how water interacts with graphene at the molecular level is essential for advancing nanomaterial applications in filtration, catalysis, and environmental technologies. This study establishes a quantitative baseline for assessing how structural defects, dopants, or surface functionalization may enhance water adsorption, providing insights for the rational design of graphene-based materials in water purification, sensing, and nanofluidic applications. In this work, we employed density functional theory (DFT) with the vdW-C09 functional to investigate the adsorption of a single water molecule on pristine graphene, accurately accounting for long-range dispersion forces. Three high-symmetry adsorption sites-the center of the hexagonal ring, the C-C bond, and the top site-were explored in combination with three molecular orientations: Down, H-bond, and Up configurations. The calculated adsorption energies range from -93 to -145 meV (milli-electron volts), indicating that the interaction is dominated by weak van der Waals forces characteristic of physisorption. The most stable configuration corresponds to the Down orientation above the center of the hexagonal ring, with an adsorption energy of -145 meV and an equilibrium distance of 3.27 A (angstrom), defined as the vertical separation between the oxygen atom of the water molecule and the graphene plane. These results are in close agreement with previous theoretical studies and confirm the non-reactive and hydrophobic nature of pristine graphene.
Stoichiometric materials of Eu(III) offer a promising platform for quantum memories attributable to their unique capability to display a distinctive, nondegenerate J = 0 transition, which enables precise mapping of optical quantum states into their hyperfine structure for reliable storage and retrieval on demand. However, placing Eu(III) into chiral polar structures, which are necessary for achieving narrow spectral linewidths and long optical lifetimes, is a daunting task. Here, we discover Eu2(SeO3)2(SO4)(H2O)2, a rare Eu(III) material that exhibits chiral polar symmetries encompassing both local and global structures. This unique structure is shaped by an appropriate combination of asymmetric ligands. The chirality fosters dipole-dipole interactions and J-mixing, as characterized by second-harmonic generation, photoluminescence, and magnetic susceptibility. The broken inversion symmetry is supported by the phase-matching behavior of second-harmonic generation. The J = 0 transition is observed at 578 nm with a narrow linewidth at 78 K and a microsecond-scale optical lifetime. The analysis of magnetic susceptibility data using Van Vleck theory results in an effective magnetic moment of 3.33 μB/Eu3+ and J-mixing. Heat capacity data reveal underlying phonon dynamics in the material. This study demonstrates a pathway toward realizing new stoichiometric Eu3+ compounds with potential for optically addressable quantum memory applications.
Many metal manufacturing processes involve phase change phenomena, which include melting, boiling, and vaporization. These phenomena often occur concurrently. A prototypical 1D model for understanding the phase change phenomena is the Stefan problem. There is a large body of literature discussing the analytical solution to the two-phase Stefan problem that describes only the melting or boiling of phase change materials (PCMs) with one moving interface. Density-change effects that induce additional fluid flow during phase change are generally neglected in the literature to simplify the math of the Stefan problem. In our recent work [1], we provide analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation in Cartesian coordinates. Our current work builds on our previous work to solve a more challenging problem: the three-phase Stefan problem in spherical coordinates for finite-sized particles. There are three moving interfaces in this system: the melt front, the boiling front, and the outer boundary which is in contact with the atmosphere. Although an analytical solution could not be found for this problem, we solved the governing equations using a fixed-grid sharp-interface method with second-order spatio-temporal accuracy. Using a small-time analytical solution, we predict a reasonably accurate estimate of temperature (in the three phases) and interface positions and velocities at the start of the simulation. Our numerical method is validated by reproducing the two-phase nanoparticle melting results of Font et al. [2]. Lastly, we solve the three-phase Stefan problems numerically to demonstrate the importance of kinetic energy terms during phase change of smaller (nano) particles. In contrast, these effects diminish for large particles (microns and larger).