2025-06-16 | | Total: 25
Recent advances in Foundation Models for Materials Science are poised to revolutionize the discovery, manufacture, and design of novel materials with tailored properties and responses. Although great strides have been made, successes have been restricted to materials classes where multi-million sample data repositories can be readily curated (e.g., atomistic structures). Unfortunately, for many structural and functional materials (e.g., mesoscale structured metal alloys), such datasets are too costly or prohibitive to construct; instead, datasets are limited to very few examples. To address this challenge, we introduce a novel machine learning approach for learning from hyper-sparse, complex spatial data in scientific domains. Our core contribution is a physics-driven data augmentation scheme that leverages an ensemble of local generative models, trained on as few as five experimental observations, and coordinates them through a novel diversity curation strategy to generate a large-scale, physically diverse dataset. We utilize this framework to construct PolyMicros, the first Foundation Model for polycrystalline materials (a structural material class important across a broad range of industrial and scientific applications). We demonstrate the utility of PolyMicros by zero-shot solving several long standing challenges related to accelerating 3D experimental microscopy. Finally, we make both our models and datasets openly available to the community.
Wood's increasing role as a structural resource in sustainable materials selection demands accurate characterization of its mechanical behavior. Its performance arises from a hierarchical structure, where the dominant load-bearing component is the S2 layer of tracheid cell walls-a thick, fiber-reinforced composite of cellulose microfibrils embedded in hemicelluloses and lignin. Due to the small dimensions and anisotropic nature of the S2 layer, mechanical testing presents significant challenges, particularly in producing homogeneous stress and strain fields. In this study, we apply micropillar compression (MPC) combined with digital image correlation (DIC) to Norway spruce tracheids, enabling direct and model-free strain measurements at the cell wall scale. Micropillars were oriented at different microfibril angles (MFAs), confirming the expected dependence of stiffness and yield stress on ultrastructural alignment, with higher stiffness and yield stress at low MFAs. For these under compression fibril-aligned kink bands occurred, while shear related failure was observed at higher angles. A parameter study on the acceleration voltage of the Scanning Electron Microscope revealed that electron beam exposure significantly degrades pillar integrity, which could explain data scatter and mechanical underestimation in earlier MPC studies. By controlling imaging protocols and using DIC-based strain measurements, we report the highest direct measurements of wood cell wall stiffness to date-up to 42 GPa for MFA=0{\deg}-closer matching micromechanical model predictions compared to previous results. Findings are compared with Finite Element Method-based displacement corrections to establish a robust protocol for probing soft, anisotropic biological composites' mechanical behavior while clarifying longstanding inconsistencies in reported results of wood MPC measurements.
In passive sky-facing radiative cooling, wavelength-selective thermal emitters in the atmospheric transparency window of 8-13 $\mu$m can reach lower temperatures compared to broadband emitters, but broadband emitters always have higher cooling power when the emitter is warmer than the ambient. Here, we propose a temperature-tunable thermal emitter that switches between a wavelength-selective state -- with high emissivity only in the atmospheric transparency window of 8-13 $\mu$m -- and a broadband-emissive state with high emissivity in the 3-25 $\mu$m range, thus maintaining high cooling potential across all temperatures. We also propose a realization of such a temperature-tunable emitter using the phase transition of vanadium dioxide (VO$_2$), which can be tuned to the ambient temperature using a combination of doping and defect engineering.
A DC arc with a meltable metal anode in a near-atmospheric pressure hydrocarbon gas is an emerging method for producing single-walled carbon nanotubes (SWCNTs). In these systems, evaporation of the molten metal anode determines the formation of catalyst seed particles needed for SWCNT growth, and therefore, should be monitored, controlled, and optimized. Evaluating the anode ablation rate by weighing the anode before and after a synthesis run is unfeasible due to anode carburization in the hydrocarbon atmosphere. To overcome this, we implemented a high-speed, 2D, 2-color pyrometry for reliable temperature measurements of the molten anode in a DC arc. The obtained temperature fields were used to calculate the anode ablation rates. Results showed the importance of resolving the arc and molten pool dynamics, as well as addressing the issue of reflections. Furthermore, significant changes in ablation rates were revealed upon addition of CH4, which must be considered when scaling up the production of SWCNTs.
Precession of a converged beam during acquisition of a 4D-STEM dataset improves strain, orientation, and phase mapping accuracy by averaging over continuous angles of illumination. Precession experiments usually rely on integrated systems, where automatic alignments lead to fast, high-quality results. The dependence of these experiments on specific hardware and software is evident even when switching to non-integrated detectors on a precession tool, as experimental set-up becomes challenging and time-consuming. Here, we introduce multi-angle precession electron diffraction (MAPED): a method to perform precession electron diffraction by collecting sequential 4D-STEM scans at different incident beam tilts. The multiple diffraction datasets are averaged together post-acquisition, resulting in a single dataset that minimizes the impact of the curvature and orientation of the Ewald sphere relative to the crystal under study. Our results demonstrate that even four additional tilts improved measurement of material properties, namely strain and orientation, as compared to single-tilt 4D-STEM experiments. We show the versatility and flexibility of our few-tilt precession approach with data collected on a number of microscopes with different hardware configurations and a variety of detectors.
Indium gallium nitride (InGaN) quantum well (QW) micro- and nano-scale light-emitting diodes (LEDs) are vital to the next-generation optoelectronic applications such as ultrahigh-speed optical interconnects and near-eye augmented/virtual reality displays. However, scaling to nanoscale dimensions introduces several challenges, including an increased nonradiative recombination from surface states and defect-related emission degradation. In this work, we present a high-performance, strain-relaxed nanoscale LED array, fabricated using a top-down nanostructuring approach applied to a high-quality, all-InGaN (quantum well/barrier) LED structure grown by plasma-assisted molecular beam epitaxy (PAMBE) on bulk GaN substrates. Combining an ultra-low threading dislocation density bulk GaN substrate and effective sidewall passivation with atomic layer deposition (ALD) results in excellent diode rectification characteristics and extremely low leakage currents. Electroluminescence measurements reveal bright blue emission with a peak external quantum efficiency (EQE) of 0.46% at an elevated injection current density of ~1.2 kA/cm2.
We investigate the structural and electronic conditions conducive to superconductivity in layered nickelates using density functional theory with Hubbard corrections (DFT+$U$). For both the bilayer and 1-3 polymorphs of La$_3$Ni$_2$O$_7$, we find that the in-plane Ni-O-Ni bond angles under pressure strongly correlate with the experimentally observed superconducting transition temperature ($T_c$) dome, and may serve as a reasonable proxy. Under compressive strain, the bond angles straighten, peaking near 2\% strain-consistent with experimental reports of superconductivity in strained bilayer thin films. However, the bond angles at this strain are more bent than those achieved under hydrostatic pressure, correlating with a lower $T_c$. We show that increasing the number of NiO$_2$ layers, as in La$_4$Ni$_3$O$_{10}$, or substituting heavier rare-earth elements (e.g., Pr) raises the pressure required to reach the structural configuration associated with superconductivity. Our results indicate that these systems require higher external pressure to achieve in-plane bond straightening. Varying the on-site Coulomb interaction $U$ reveals that stronger electronic correlations delay the structural transition and favor high-spin states. This suggests that moderate correlation strength may be optimal for superconductivity, with stronger correlation preventing the formation of favorable bond geometries. Electronic structure analysis shows that the Ni $e_g$ orbitals dominate near the Fermi level and shift downward with pressure, enhancing Ni-O hybridization. These results highlight how pressure and strain tune structural features that may be essential for engineering high-$T_c$ phases in nickelate superconductors.
Magnetic hopfions are non-collinear spin textures that are characterized by an integer topological invariant, called Hopf index. The three-dimensional magnetic solitons can be thought of as a tube with a twisted magnetization that has been closed at both ends to form a torus. The tube consists of a magnetic whirl called in-plane skyrmion or bimeron. Although hopfions have been observed by microscopy techniques, their detection remains challenging as they lack an electronic hallmark so far. Here we predict a three-dimensional orbital Hall effect caused by hopfion textures: When an electric field is applied, the hopfion generates a transverse current of orbital angular momentum. The effect arises due to the local emergent field that gives rise to in-plane and out-of-plane orbital Hall conductivities. This orbital Hall response can be seen as a hallmark of hopfions and allows us to distinguish them from other textures, like skyrmioniums, that look similar in real-space microscopy experiments. While the two-dimensional topological invariant of a skyrmion determines its topological Hall transport, the unique three-dimensional topological orbital Hall effect can be identified with the three-dimensional topological invariant that is the Hopf index. Our results make hopfions attractive for spin-orbitronic applications because their orbital signatures allow for their detection in devices and give rise to large orbital torques.
We report theoretical best estimates of vertical transition energies (VTEs) for a large number of excited states and molecules: the \textsc{quest} database. This database includes 1489 \emph{aug}-cc-pVTZ VTEs (731 singlets, 233 doublets, 461 triplets, and 64 quartets) for both valence and Rydberg transitions occurring in molecules containing from 1 to 16 non-hydrogen atoms. \textsc{Quest} also includes a significant list of VTEs for states characterized by a partial or genuine double-excitation character, known to be particularly challenging for many computational methods. The vast majority of the reported values are deemed chemically-accurate, that is, are within $\pm0.05$ eV of the FCI/\emph{aug}-cc-pVTZ estimate. This allows for a balanced assessment of the performance of popular excited-state methodologies. We report the results of such benchmarks for various single- and multi-reference wavefunction approaches, and provide extensive supporting information allowing testing of other models. All corresponding data associated with the \textsc{quest} database, along with analysis tools, can be found in the associated \textsc{GitHub} repository at the following URL: https://github.com/pfloos/QUESTDB.
Harnessing spin currents to control magnon dynamics enables new functionalities in magnonic devices. Here, we demonstrate current-controlled magnon-magnon coupling between cavity and boundary modes in an ultrathin film of Bi-doped yttrium iron garnet (BiYIG). Cavity modes emerge in a BiYIG region between two Pt nanostripes, where interfacial anisotropy modifies the magnon dispersion. These modes hybridize with boundary magnons confined within the Pt-capped BiYIG, resulting in an anticrossing gap. Modeling based on dipole-exchange spin-wave dispersion accurately reproduces the observed modes and their hybridization. Spin current injection via the spin Hall effect in a Pt nanostripe disrupts the cavity boundary conditions and suppresses both cavity modes and hybridization upon driving the system beyond the damping compensation threshold. Furthermore, tuning the microwave power applied to a microstrip antenna enables controlled detuning of the anticrossing gap. Our findings provide a platform for exploring spin current-magnon interactions and designing on-chip reconfigurable magnonic devices.
We report the detection of pure orbital currents generated by both coherent and thermal magnons in the magnetic insulator Bi-doped yttrium iron garnet (BiYIG). The pumping of orbital and spin currents is jointly investigated in nano-devices made of naturally oxidized Cu, pure Cu, Pt, and Cr. The absence of charge conduction in BiYIG and the negligible spin-to-charge conversion of oxidized Cu allows us to disambiguate the orbital current contribution. Comparative measurements on YIG and BiYIG show that the origin of the orbital pumping in BiYIG/oxidized Cu is the dynamics of the orbital magnetization in the magnetic insulator. In Cr, the pumping signal is dominated by the negative spin Hall effect rather than the positive orbital Hall effect, indicating that orbital currents represent a minority of the total angular momentum current pumped from the magnetic insulator. Our results also evidence that improving the interfacial transparency significantly enhances pumping efficiencies not only for spin, but also for orbital currents.
The thermal properties of solids under nanoscale confinement are currently not understood at the atomic level. Recent numerical studies have highlighted the presence of a minimum in the thermal conductivity as a function of thickness for ultrathin films at a thickness about 1-2 nm, which cannot be described by existing theories. We develop a theoretical description of thin films which predicts a new physical law for heat transfer at the nanoscale. In particular, due to the strong redistribution of phonon momentum states in reciprocal space (with a transition from a spherical Debye surface to a different homotopy group $\mathbb{Z}$ at strong confinement), the low-energy phonon density of states no longer follows Debye's law but rather a cubic law with frequency, which then crosses over to Debye's law at a crossover frequency proportional to the average speed of sound of the material and inversely proportional to the film thickness. Concomitantly, this implies that the phonon population becomes dominated by low-energy phonons as confinement increases, which then leads to a higher thermal conductivity under extreme confinement. The theory is able to reproduce the thermal conductivity minimum in recent molecular simulations data for ultrathin silicon and provides useful guidelines as to tune the minimum position based on the mechanical properties of the material.
Laser-induced inertial cavitation (LIC)-where microscale vapor bubbles nucleate due to a focused high-energy pulsed laser and then violently collapse under surrounding high local pressures-offers a unique opportunity to investigate soft biological material mechanics at extremely high strain rates (>1000 1/s). Traditional rheological tools are often limited in these regimes by loading speed, resolution, or invasiveness. Here we introduce novel machine learning (ML) based microrheological frameworks that leverage LIC to characterize the viscoelastic properties of biological materials at ultra-high strain rates. We utilize ultra-high-speed imaging to capture time-resolved bubble radius dynamics during LIC events in various soft viscoelastic materials. These bubble radius versus time measurements are then analyzed using a newly developed Bubble Dynamics Transformer (BDT), a neural network trained on physics-based simulation data. The BDT accurately infers material viscoelastic parameters, eliminating the need for iterative fitting or complex inversion processes. This enables fast, accurate, and non-contact characterization of soft materials under extreme loading conditions, with significant implications for biomedical applications and materials science.
In this article, we provide a brief perspective on recent developments in the study of shear thickening in dense suspensions. We give a rapid overview of the state of the art and discuss current models aiming to describe this particular rheology. Although most of the experiments and simulation studies are conducted in "ideal" flows, where the sample is confined without an open boundary condition, we have decided to highlight more realistic flow conditions. We further provide an overview on how to relate the recently proposed constitutive models to these more practical flow conditions like pipe flow or flow down an incline.
Mixed ionic-electronic conductor (MIEC) exhibit both high ionic and electronic conductivity to improve the battery performance. In this work, we investigate the mechanism and stability of transport channels in recently developed our MIEC material, amorphous Ti-doped lithium phosphorus sulfide (LPS), using molecular dynamics (MD) simulations with a 99\% accurate machine-learning force field (MLFF) trained on \textit{ab-initio} MD data. The achieved MLFF helps efficient large-scale MD simulations on LPS with three Ti concentrations (10\%, 20\%, and 30\%) and six temperatures (25$^\mathrm{o}$C to 225$^\mathrm{o}$C) to calculate ionic conductivity, activation energy, Li-ion transport mechanism, and configurational entropy. Results show that ionic conductivities and activation energies are consistent with our recent experimental values. Moreover, Li-ion transport occurs via free-volume diffusion facilitated by the formation of disordered Li-S polyhedra. Enhanced stability of transport channels at 10\% and 20\% Ti doping, compared to 0\% and 30\%, is observed through the analysis of the vibrational and configurational entropy of these disordered Li-S polyhedra. Overall, this study highlights the utility of MLFF-based large-scale MD simulations in explaining the transport mechanism and its stability of Li-ion in Ti doped LPS electrolyte with significant computational efficiency.
Chalcogenide perovskites, particularly II-IV ABS3 compounds, are a promising class of materials for optoelectronic applications. However, these materials frequently exhibit instability in two respects: 1) a preference for structures containing one-dimensional edge- or face-sharing octahedral networks instead of the three-dimensional corner-sharing perovskite framework (polymorphic instability), and (2) a tendency to decompose into competing compositions (hull instability). We evaluate the stability of 81 ABS3 compounds using Density Functional Theory, finding that only BaZrS3 and BaHfS3 are both polymorphically and hull stable, with the NH4CdCl3-type structure the preferred polymorph for 77% of these compounds. Comparison with existing tolerance factor models demonstrates that these approaches work well for known perovskites but overpredict stability for compositions without published experimental results. Polymorphic stability analysis reveals that perovskite structures are stabilized by strong B-S bonding interactions, while needle structures exhibit minimal B-S covalency, suggesting that electrostatic rather than covalent interactions drive the preference for edge-sharing motifs. Hull stability analysis comparing ABS3 to ABO3 analogues reveals a weaker inductive effect in sulfides as a possible explanation for the scarcity of sulfides compared with oxides. The relative instability of ABS3 compounds is further supported by experimental synthesis attempts. These findings provide fundamental insights into the origins of instability in chalcogenide perovskites and highlight the challenges in expanding this promising materials class beyond the few materials that have been reported to date.
In crystals with broken time-reversal symmetry, zone-center phonons can acquire a finite angular momentum via velocity-dependent forces on the nuclei. Despite having the same order of magnitude as the electron spin angular momentum, the phonon angular momentum can be hard to detect because the frequency splitting is small. Here, by developing a theory of lattice magneto-optical effects in reflection and transmission, we show that infrared magnetic circular dichroism is a sensitive probe of zone-center phonon chirality. We evaluate the infrared magneto-optical Faraday, Kerr, and circular-dichroism spectra of CrI$_3$ from time-dependent density-functional theory in the adiabatic local-density approximation. We find sizeable circular dichroism from the infrared-active E$_u$ mode at $\approx 214$ cm$^{-1}$, even though the calculated splitting is only 0.22 cm$^{-1}$.
Crystalline materials can form different structural arrangements (i.e. polymorphs) with the same chemical composition, exhibiting distinct physical properties depending on how they were synthesized or the conditions under which they operate. For example, carbon can exist as graphite (soft, conductive) or diamond (hard, insulating). Computational methods that can predict these polymorphs are vital in materials science, which help understand stability relationships, guide synthesis efforts, and discover new materials with desired properties without extensive trial-and-error experimentation. However, effective crystal structure prediction (CSP) algorithms for inorganic polymorph structures remain limited. We propose ParetoCSP2, a multi-objective genetic algorithm for polymorphism CSP that incorporates an adaptive space group diversity control technique, preventing over-representation of any single space group in the population guided by a neural network interatomic potential. Using an improved population initialization method and performing iterative structure relaxation, ParetoCSP2 not only alleviates premature convergence but also achieves improved convergence speed. Our results show that ParetoCSP2 achieves excellent performance in polymorphism prediction, including a nearly perfect space group and structural similarity accuracy for formulas with two polymorphs but with the same number of unit cell atoms. Evaluated on a benchmark dataset, it outperforms baseline algorithms by factors of 2.46-8.62 for these accuracies and improves by 44.8\%-87.04\% across key performance metrics for regular CSP. Our source code is freely available at https://github.com/usccolumbia/ParetoCSP2.
We present a Semiempirical Pseudopotential Method for accurately computing the band structures and Bloch states of monolayer transition-metal dichalcogenides (TMDCs), including MoS2, MoSe2, WS2, and WSe2. Our approach combines local and non-local pseudopotentials, carefully fitted to replicate fully self-consistent density-functional theory results, while using only a minimal set of empirical parameters. By expressing the total potential as a sum of a few separable components, we achieve both accuracy and computational efficiency. The resulting pseudopotentials are transferable to more complex systems such as bilayer, trilayer, and moire superlattices, offering a reliable foundation for large-scale simulations. When integrated with artificial intelligence techniques, this method provides a powerful tool for the design and exploration of TMDC-based nanodevices for next-generation optoelectronic applications.
Uni-traveling-carrier photodiodes (UTC-PDs), which utilize only electrons as the active carriers, have become indispensable in high-speed optoelectronics due to their unique capabilities, such as high saturation power and broad bandwidth. However, extending the operating wavelengths into the visible region for wider applications is challenging due to the lack of suitable wide-bandgap III-V semiconductor combinations with the necessary band alignment and lattice matching. Here, we show that a UTC-PD based on a van der Waals heterojunction composed of a 2D transition metal dichalcogenide, molybdenum disulfide (MoS2), as a photoabsorption layer and a gallium nitride (GaN) film as a carrier transport layer, offers a solution to this challenge. The fast vertical carrier transport across the heterointerface is enabled by the direct epitaxial growth of a MoS2 layer on a GaN film. Our device demonstrates a frequency response in the several-GHz range with a quantum efficiency on the order of 1% throughout the entire visible spectrum, highlighting the promise for high-speed visible optoelectronics.
The development of $p$-type oxide semiconductors remains impeded by the inherently low-lying valence-band maximum (VBM) dominated by O-2$p$ states. A prevailing approach to mitigate this limitation is to elevate the VBM by introducing cation states that hybridize with O-2$p$ orbitals or lie energetically above the O-2$p$ level. Nevertheless, the $p$-type oxides reported to date exhibit limited hole mobilities. To expand the search space, it is essential to accurately understand the intrinsic difficulty of introducing holes into O-2$p$ orbitals. Accordingly, we evaluated 845 oxides to identify those in which holes can be doped into O-2$p$ bands. Our high-throughput screening revealed CaCdO$_2$ as the only promising exemplar, in which the VBM is slightly hybridized with deep-lying Cd-3$d$ states. Our screening suggests that hole doping into O-2$p$ orbitals is extremely difficult and thus reinforces the effectiveness of the traditional ``VBM-raising strategy.''
The catalytic activity of metal nanoparticles toward nucleation of single-walled carbon nanotubes (SWCNTs) is fundamental to achieving structure-controlled growth using catalytic chemical vapor deposition (CVD). Despite the success of oxidized catalysts in SWCNT growth, there is a lack of understanding regarding how oxygen influences the catalysts and the nucleation process. Quantum chemical molecular dynamics (MD) simulations employing density functional tight binding (DFTB) demonstrate that the kinetics of carbon nucleation on an iron nanoparticle catalyst can be tuned via oxygen loading. Increasing the oxygen content in the catalyst leads to activation of surface-bound carbon species and enhanced carbon chain growth due to respective weakening and strengthening of the C-C and Fe-C bonding. This is due to oxygen modulating the electronic structure of the iron catalyst, with the Fermi level of the catalyst increasing proportionally with oxygen content until the iron:oxygen stoichiometry reaches parity. The increase in Fe 3d states near the Fermi level also promotes the donation of electron density into unoccupied C 2p states, activating C-C bonds which in turn facilitates carbon chain growth and slows carbon ring condensation.
An impact of thin metallic ZrN layers on Si and sapphire wafers on substrate temperature during MBE growth of GaN nanowires is studied. Using nucleation kinetics of GaN as a sensitive probe we show that a thin ZrN layer strongly increases the substrate temperature, which significantly affects the dimensions and density of the nanowires. To quantify the effect we developed a technique of optical pyrometer calibration that allows reliable determination of emissivity, and thus precise measurement of temperature of substrates with unknown optical parameters, such as ZrN buffers of various thicknesses. Our results show that emissivity of ZrN-coated Si and sapphire wafers differs significantly from the bulk ZrN and increases drastically for films thinner than ~100 nm. Simple calculations indicate that ignoring the influence of the thin film may lead to huge errors in temperature readings and consequently to losing the growth control. Then, we show that we can compensate for the impact of ZrN buffer on substrate temperature and grow identical nanowire arrays on Si substrates with and without ZrN layers. Finally, having identical arrays of GaN nanowires we used X-ray diffraction to compare nanowire arrangements on Si and ZrN/Si substrates with a thin SiN nucleation layer.
The theoretical picture of thermal field-emission (TFE) from semiconductors has been limited to 1D and 2D models. This can be attributed to the complex and interdependent phenomena that is involved in TFE from semiconductors which makes the calculations cumbersome. Such limitations result in a partial understanding of the underlying physics of semiconducting surfaces under high electrical fields, which requires the addition of the temporal dimension (4D) to yield a realistic model. Here we develop a 3D model of TFE from semiconductors that can take arbitrary geometries and doping levels. Our model successfully reproduces the characteristic saturation plateau of some semiconductors, as well as its dependence in temperature. The model is found to be in good agreement with experimental data from ntype Germanium at a qualitative level. We propose this model as a platform for future extensions into the full 4D framework, incorporating temporal dynamics for a more complete and predictive description of thermal-field emission from semiconductors.
Nonadiabatic Thouless pumping of electrons is studied within the framework of topological Floquet engineering, particularly focused on how atomic lattice dynamics affect the emergent Floquet topological phase in trans-polyacetylene under the driving electric field. As similarly done in the earlier work [Zhou and Kanai, J. Phys. Chem. Lett., 12, 4496 (2021)], the real-time time-dependent density functional theory and Ehrenfest dynamics simulations were used to investigate the extent to which the number of pumped charges remains equal to the topological invariant, the winding number, when the temperature effect of ions and the dynamical coupling of electrons and ions are taken into account. Our theoretical work shows that the Floquet topological phase remains intact but the condition on the driving field necessary for observing the topological phase becomes more limiting.