2026-05-29 | | Total: 40
We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.
Conjugated polymers exhibit unique spin-dependent phenomena arising from weak yet critical hyperfine interactions. Understanding these spin effects, particularly the spin-dependent formation and decay of correlated spin pairs, is important for advancing both organic electronics and polymer-based spintronics. Intrinsic magnetic-field responses such as magnetoresistance have primarily been investigated in diode architectures, where electrons and holes coexist. However, such systems are less suitable for probing bipolaron formation in unipolar transport, and the relationship between polymer structure and bipolaron formation in lightly doped polymers remains unclear. Here, we systematically investigate intrinsic magnetoresistance in representative conjugated polymers using field-effect transistors and observe a generally positive magnetoresistance. First-principles simulations reveal that bipolarons preferentially form on short conjugated segments associated with amorphous regions. Moreover, comparisons across these polymers show that enhanced charge trapping correlates with stronger magnetoresistance, implying promoted bipolaron formation. Bipolaron-incorporated energylevel-alignment modeling near metal/polymer interfaces suggests that charge traps can increase the bipolaron density.
Message-passing neural networks (MPNNs) are widely used for molecular property prediction, but their deployment as monolithic architectures makes it difficult to identify how specific message-passing operators affect performance. We present an operator-level factorial benchmark that decomposes 2D molecular MPNNs into the three families of message-seed initialization, node-edge fusion, and node update operators. The resulting 84 configurations are benchmarked on ten MoleculeNet datasets under a shared experimental setup and statistical analysis protocol. Across this controlled design, performance variation is associated primarily with message construction rather than update complexity. Message-seed initialization shows significant family-level effects for both regression and classification, node-edge fusion shows a significant family-level effect for regression with descriptive advantages for concatenation-based mixing, and the update family shows no statistically supported effect for either endpoint family. A representation probe into the Quinethazone molecule further demonstrates that concatenation-based mixing can better differentiate chemically distinct heteroatoms and withstand oversmoothing than Hadamard gating. Representative configurations selected separately for classification and regression recover competitive performance relative to established molecular graph neural network (GNN) baselines, ranking numerically best on eight of ten benchmark datasets. These empirical results are interpreted through concise mechanistic analyses of representative node-edge fusion and update operators. Our findings provide empirical design heuristics for molecular MPNNs by turning model design from a search over monolithic architectures into a targeted assessment of where and how chemical information enters the message-passing pipeline.
Ternary transition-metal nitrides with layered crystal structures host anisotropic bonding and reduced dimensionality that may enable unconventional electronic and magnetic behavior. Yet, synthesis of such nitride thin films remains challenging because reactive sputtering often favors metastable rocksalt-derived structures. We report the composition-dependent synthesis, structure, and properties of layered FeMN2 (M = W, Mo) thin films with triangular Fe sublattices, prepared by reactive sputtering and post-deposition NH3 annealing. Using synchrotron grazing-incidence wide-angle X-ray scattering (GIWAXS) and X-ray absorption spectroscopy (XAS), we investigate the evolution of phase formation, crystallographic texture, and local Fe coordination across composition. Both layered phase form over broad composition ranges. FeWN2 maintains high phase purity across compositions, consistent with cation substitution, whereas FeMoN2 only exhibits good phase purity at Fe-poor compositions. Azimuthal GIWAXS analysis shows that Fe-rich films in both systems exhibit out-of-plane fiber texture, which progressively evolves toward predominantly in-plane orientation near stoichiometry in FeWN2, while FeMoN2 develops a more randomly oriented polycrystalline microstructure. Electrical measurements reveal relatively low and composition-insensitive resistivity in FeWN2 (~1 mΩ.cm), whereas FeMoN2 exhibits a pronounced resistivity maximum near nominal stoichiometry. Preliminary room-temperature magnetization measurements on FeWN2 further reveal weak ferromagnetic-like behavior in Fe-poor films, while stoichiometric compositions remain predominantly paramagnetic. These results demonstrate fundamentally different structural accommodation mechanisms in FeWN2 and FeMoN2, and highlight the strong coupling between composition, microstructure, and electronic/magnetic properties in layered nitride thin films.
Noncollinear (NC) magnetism and spin-orbit coupling (SOC) are indispensable for predictive ab initio materials simulations with pronounced relativistic effects and magnetic frustration, yet they significantly increase the cost of cubic-scaling density functional theory (DFT) by introducing complex 2-component wavefunctions per electron and consequently much larger eigenproblems. We present a GPU-centric high-performance framework for NC-SOC DFT that combines: (i) algorithmic advances for solving finite-element (FE) discretized DFT equations; (ii) residual-based Chebyshev filtered subspace iteration (R-ChFSI), tolerant to inexact matrix-vector products, for the resulting sparse generalized eigenproblem; (iii) a matrix-free strategy for accelerating FE Poisson solver; (iv) R-ChFSI-enabled mixed-precision computation with block floating-point compressed MPI communication at compression ratios over 4x, preserving double-precision robustness while reducing compute and data movement costs; and (v) a communication efficient band-partitioning algorithm to improve scalability. Numerical results demonstrate improved time-to-solution and excellent scaling on exascale architectures, enabling fully relativistic pseudopotential DFT simulations of up to 100,000 electrons.
We investigate the propagation of longitudinal elastic waves in one-dimensional periodic composite rods composed of alternating segments of a shape-memory alloy (NiTiCu) and a polymer spacer (Parylene C). In the thin-rod regime, the longitudinal phase velocity reduces to $c=\sqrt{E/ρ}$, which coincides with the regime in which the elastic modulus of NiTiCu has been measured directly through its acoustic response across the martensitic transformation. Using the standard transfer-matrix method along the heating and cooling branches of the transformation separately, we compute the Bloch band structure of the infinite periodic system and the transmission spectrum of finite composite rods. Because the elastic modulus of NiTiCu follows different paths upon heating and cooling, the same external temperature within the transformation interval corresponds to two different phase fractions and, consequently, to two different phononic spectra. The resulting hysteresis of the underlying material is thus transferred to the collective acoustic response of the periodic structure: stop-band edges trace closed loops in the temperature--frequency plane, and the transmission coefficient of a finite rod at a fixed temperature depends on the previous thermal history. We further show that the geometric filling fraction of the active segment provides a complementary tuning mechanism, modifying the width of the spectral hysteresis loops and the position of specific gap closures independently of temperature. These results illustrate how a first-order structural phase transition with intrinsic thermal hysteresis manifests itself in the dispersion relation of a periodic elastic medium.
Generating and characterizing uniform and staggered spin polarization in antiferromagnets is one of the key challenges for antiferromagnetic spintronic technology. Here, we perform perturbative theory, group-theoretical symmetry analysis, low energy and ab initio simulations to propose that the magnetic precession near the equilibrium magnetic axis could generate finite uniform and staggered spin polarization at the opposite magnetic sublattices (referring to total magnetic and Néel vector generation) in a single AFM semiconductors. This response does not require the heterojunction setup and could eliminate the lattice mismatch issues at the junction. Through scrutinizing all symmetrically-protected vanishing magnetic moment groups and especially focusing on parity-time (PT ) invariant groups, we identify the symmetry constraints that describe the staggered spin accumulation responses, and disclose their fieldlike and dampinglike characters. This unravels a hidden spin accumulation mode in AFM semiconductors. Furthermore, we simulate such an effect using a perturbative approach and suggest that electric gate field and Floquet light-dressing can effectively manipulate these responses.
Alloy-based perovskite solar cells offer tunable properties and improved stability, but their complexity has impeded accurate modeling, hindering development. We present a machine-learning (ML) accelerated atomistic modeling approach for the phase stability of (Cs/FA)Pb(Br/I)3 and (Cs/FA)Sn(Br/I)3 perovskites, with FA being formamidinium. To make such quaternary alloys tractable, we adopt a two-level ML strategy, combining 1) graph neural network interatomic potentials trained on density functional theory data for efficient structure relaxations with 2) secondary ML models for direct energy prediction from unrelaxed structures. These models enable computations of free energy landscapes across compositions and phases, capturing alloy disorder and FA molecular orientations. Our results reveal narrower stable composition regions for the Sn-based system compared to its Pb-based counterpart, limiting options for compositional engineering. Maximum stability occurs at high I content, and no stabilization is observed near the center of the composition space. Our results guide the design of stable perovskites.
Machine learning interatomic potentials (MLIPs) have transformed materials discovery by leveraging graph neural networks (GNNs) to predict material properties with near density functional theory (DFT) accuracy. While large-scale pretrained foundation models offer transferable baseline representations, they frequently struggle to generalise to out-of-distribution (OOD) target systems -- a common challenge in modelling complex or chemically diverse materials. Fine-tuning is the standard remedy, but the high cost of generating DFT-labelled configurations confines adaptation to data-scarce regimes, where over-parameterised GNNs amplify overfitting and degrade target-domain performance. To address this, we propose a prototype-based alignment approach for data-efficient fine-tuning of MLIPs. Our method identifies local structural similarities between the source and target domains by grouping atoms with analogous chemical environments based on their latent representations. Each target-domain atom's energy contribution is aligned to its source-domain prototype, introducing an inductive bias that anchors fine-tuned representations to the pretrained structure, encouraging effective reuse of learned interactions and improving generalisation without restrictive assumptions on the target chemistry. We evaluate our method on the rMD17 benchmark using equivariant MACE and invariant SchNet across varying data budgets, and extend evaluation to the MACE-OFF foundation models on the SPICE dataset. Our approach consistently improves predictive accuracy in the low-data regime, reducing energy MAE by up to 18% over standard fine-tuning baselines.
Pnictogen-based solar absorbers have gained prominence as promising nontoxic and stable alternatives to lead-halide perovskites (LHPs), but are severely limited by carrier localization, preventing their performance from approaching those of LHPs. Recent efforts have uncovered routes to overcome carrier localization, but these early efforts only considered intrinsic factors. Herein, we push beyond these limited early efforts, examining the role of defects, not only on cold carriers but also hot carriers. Focusing on the structurally one-dimensional pnictogen chalcohalide BiSBr, we find that whilst this material intrinsically does not exhibit carrier localization, vacancies introduced during synthesis or post-treatment lead to pronounced extrinsic self-trapping via the formation of defect-bound hot polarons-excited charge-carriers strongly coupled to local defect-induced vibrational modes. These above-gap defect states divert hot carriers from cooling to the band edge, thus depleting the mobile carrier population. Our findings establish the key role of defect-bound hot polarons in mediating extrinsic localization and offer new mechanistic insights into the interplay between defects, lattice coupling, and excited-state charge-carrier transport, which are critical to designing efficient perovskite-inspired solar absorbers.
Nontrivial topological spin-textures, such as skyrmions, merons, bimerons, and skyrmioniums, are envisioned as robust building blocks for future memory and logic devices. Controllable transformations between these states require a quantum-mechanical description of electronic degrees of freedom and atomic-scale insight beyond existing phenomenological models. Here, we report an atomic-scale investigation of topological phase transitions and their protection in the two-dimensional van der Waals ferromagnet Fe$_3$GeTe$_2$ (FGT) using a combined experimental-theoretical approach. Synchrotron-based Fourier transform holography directly images labyrinth domains, isolated skyrmions, mixed labyrinth-skyrmion phases, and skyrmion bags with high spatial resolution. We compare these observations to simulations based on an electronic lattice Hamiltonian that captures both metallicity and relativistic spin-orbit coupling in FGT. By systematically exploring a broad range of temperatures and magnetic fields, we map the mechanisms governing topological transitions and their stability. This sequential-integrated experimental-theoretical framework advances understanding of spin-texture interactions and enables precise control of external tuning parameters. Our results establish a platform for creating, stabilizing, and manipulating topological states, paving the way for engineered spin-texture transitions in next-generation spintronic technologies.
The efficiency of optical cooling of nuclear spins in a CdTe/CdZnTe quantum well is investigated as a function of an external magnetic field. Our results confirm that there is indeed an optimal external magnetic field for optical cooling. We associate it with the kinetic local field $B_{KL}$ defined by the heating rate of the spin-spin reservoir due to the fluctuations of the hyperfine interaction. We also propose an experimental technique for measuring $B_{KL}$. For our sample we find that $B_{KL}=1.0\pm0.4$ G and it is independent of the electron polarization and pump power. The measured values of the kinetic local fields are in good agreement with a theoretical calculation $B_{KL} = 0.7$ G, taking into account indirect spin-spin interactions of Cd and Te nuclear spins and their considerably different hyperfine constants. The hyperfine constants of the magnetic isotopes of Cd and Te in CdTe are estimated.
The hydrogen evolution reaction (HER) is central to sustainable hydrogen production, and nitrogen coordinated dual atom catalysts (DACs) offer a promising route to noble metal activity at low cost. Yet their vast compositional and coordination design space remains underexplored, as density functional theory (DFT) screening at scale is prohibitive. Here, we map the HER landscape of graphene supported TM2@Nx-Gr DACs, screening 23 transition metals across 20 nitrogen coordination motifs using a machine learning potential (MLP) benchmarked against DFT. Intermediate coordination (2N to 4N) consistently yields near-optimal ΔGH*, with Ti2@2Na, Mn2@2Na, Fe2@2Na, Cu2@2Na, Rh2@2Na, Zr2@2Na, Zr2@2Nb, Zr2@2Nc, Nb2@2Nc, Zr2@2Nd, Mn2@2Ne, Mn2@2Nf, Ti2@3Na, Au2@3Na, Fe2@3Na, Pd2@3Nb, Rh2@3Nc, Rh2@3Nd, Au2@3Nd, V2@4Na, Ti2@4Nb, Pd2@4Nb, Ti2@4Nc, Cr2@4Nd, Ni2@4Nd, Cu2@4Nd emerging as standout, synthesizable candidates, most exhibiting metallic or narrow gap (<0.25 eV) character. The MLP reaches near-DFT accuracy, with a mean absolute error of 80 meV for Gibbs binding free energies at orders of magnitude lower computational cost, establishing MLP driven screening as a practical engine for next-generation catalyst discovery.
Hybrid structure of two zigzag-edge graphene nanoribbons with a fused porphine ring in between, results in two distinct nearly degenerate ground states: a semiconducting antiferromagnetic state and a conducting ferromagnetic state with unequal and opposite Fermi velocities of majority and minority spins, the former having slightly higher stability. Such ground states result from the broken bipartite symmetry induced by the porphine ring. The incorporation of different transition metal atoms in the porphine cavity reduces their energy difference but keeps their electronic properties mostly unchanged. The splitting of the $d$-orbitals in the distorted square-planar ligand field of porphine produces a high spin ground state that breaks the global time reversal symmetry ($\mathcal{T}$). The opposite Fermi velocities of the majority and minority spins in the ferromagnetic ground state and lower sensitivity of the conducting majority spin channel to the edge disorder, make this class of quasi-one-dimensional hybrid structures promising for dual spin-filtering device applications.
Spin-orbit coupling (SOC) plays an important role in determining the structural and electronic properties of recently proposed two-dimensional planar pentagonal materials. In this work, density functional theory calculations are employed to investigate SOC effects in p-MS$_{2}$ systems (M = Si, Ge, and Pb). Our results indicate that the p-SiS$_{2}$ structure is likely unstable, except for p-GeS$_{2}$ and p-PbS$_{2}$. A detailed j-resolved (total angular momentum) orbital analysis reveals that SOC enhances electronic localization, leading to a slight structural contraction and a reconstruction of electronic states near the Fermi level, this effect becoming stronger for heavier M atoms. While p-GeS$_{2}$ remains metallic, SOC drives a metal-semiconductor transition in p-PbS$_{2}$ and opening a quasi-direct band gap of about 0.475 eV. In addition, the conduction band minimum state of p-PbS$_{2}$ exhibits pronounced anisotropy along the S-S bonds. These findings provide insight into SOC-driven structural and electronic reconstruction in planar pentagonal chalcogenides p-MS$_{2}$ and suggest that p-PbS$_{2}$ may be a promising candidate for gas-sensing applications.
Doping in two-dimensional materials has emerged as an effective tool for modulating their electronic properties and thereby enabling their multifunctional applications. In this work, we present a first-principles study on induced effective magnetic moment and metallicity in MoS$_2$ monolayer by substitutional doping of group-5 transition metal (TM) elements -- V, Nb and Ta. From our study, we observe that the V doping induces half-metallicity, whereas metallic characteristics are observed in the case of Nb and Ta doping. Moreover, V and Ta-doped MoS$_2$ monolayers are observed to show total induced magnetic moments of 0.922 and 0.624 $μ_{\rm B}$, respectively. Importantly, the combined effects of strong spin-orbit coupling (SOC), broken inversion symmetry, and structural asymmetry is observed to lead to a permanent valley polarization in the V- and Ta-MoS$_2$ systems. In particular, we observed a valley polarization of 121 and 21 meVs for V and Ta-doped MoS$_2$, respectively. Furthermore, an enhanced piezoelectric coefficient for the doped systems is observed compared to pristine MoS$_2$. Notably, the simultaneous presence of half-metallicity, substantial valley polarization, and enhanced piezoelectricity in V-doped MoS$_2$ establishes this system as a promising multifunctional platform for next-generation spintronic, valleytronic, and piezoelectric nanodevices. Overall, our findings provide fundamental insights into engineering coupled spin-valley-mechanical degrees of freedom in two-dimensional materials for advanced quantum and nanoelectronic applications.
A multi-scale finite-temperature micromagnetic model is presented, based on the Landau-Lifshitz equation and the Bernoulli differential equation. This model accurately reproduces classic Maxwell magnetostatics of paramagnets for high temperatures and accurately reproduces standard micromagnetics described by the conventional Landau-Lifshitz model in ferromagnets. The Landau-Lifshitz-Bernoulli (LLBe) model can, by design, directly couple atomic-scale simulations with micromagnetics and output consistent predictions of bulk magnetic properties at finite temperatures, from below to above the material's Curie temperature. The LLBe model is validated against established solvers: MUMAX3 for zero-temperature micromagnetics, and FEMCE for high-temperature classic magnetostatics. We present an application of the LLBe model by simulating Heat-Assisted magnetic recording on a thin magnetic track with local heating, demonstrating the multi-scale finite-temperature capabilities of the LLBe.
Since its original prediction 25 years ago, room-temperature out-of-plane ferroelectricity in compressively strained SrTiO3 remains an ongoing pursuit. In this work, we investigate the structural, electrical and electromechanical properties of highly strained epitaxial SrTiO3 capacitors with rare earth nickelate electrodes. The SrTiO3 layers experience compressive strains up to -3% and exhibit pronounced tetragonality, comparable to that of bulk PbTiO3. Variable-temperature electrical measurements and room-temperature piezoresponse force microscopy reveal butterfly-shaped capacitance-voltage hysteresis and domain-like electromechanical response typical of ferroelectric materials. However, the overall behavior is inconsistent with a stable ferroelectric state. We therefore propose an alternative mechanism for the observed polarity in our samples based on spatially inhomogeneous internal fields. Our first-principles calculations show that such fields may arise from charge discontinuities between the formally charged NdNiO3 layers and charge-neutral SrTiO3 layers.
Multicomponent transition metal carbides are promising for extreme-environment applications, but identifying compositions that are both synthesizable and hard remains challenging. We fine-tune the MACE machine-learned interatomic potential on approximately 28,000 density functional theory calculations spanning the composition space of groups 4-6 transition metals and carbon to predict the thermodynamic stability and elastic properties of multicomponent carbides. The fine-tuned model achieves formation energy errors of ~ 10 meV/atom for thermodynamically relevant structures with only 20% of the training data. We screen over 1500 equiatomic compositions across rocksalt, hexagonal, and hcp prototypes, combining free energy models with elasticity-based hardness surrogates. Synthesizability predictions at 1500°C agree well with experimental reports for both single-phase and multiphase carbides. The group number of the constituent metals governs both stability and hardness. Free energy contributions from short-range order are small, typically a few meV/atom, indicating that a perfectly disordered solid solution provides a reasonable approximation for high-throughput screening. For compositions mixing group 4/5 and group 6 metals, we identify a new family of stacking-ordered phases with formation energies well below those of disordered rocksalt or hexagonal structures. DFT calculations corroborate these predictions and suggest that stacking-ordered phases should be experimentally accessible in multicomponent carbides. This study provides a framework for screening synthesizable multicomponent materials with target properties, identifies promising carbide compositions across the full nine-component space, and reveals a new class of stacking-ordered carbides accessible only in multicomponent compositions.
Predictive modeling of the infrared dielectric function in polar materials is crucial for thermal management and infrared devices design. While the Green-Kubo molecular dynamics (MD) framework provides a nonperturbative route to compute dielectric responses from dipole fluctuations, yet it commonly relies on the fixed-charge approximation that neglects dynamic charge redistribution during atomic motion. Here, we employ a machine-learning neuroevolution potential with dynamic charges (qNEP) combined with Green-Kubo MD to investigate the dynamic charge effect on the infrared dielectric response of rutile TiO$_2$, a material with large Born effective charges (BEC). Our results show that dynamic charge effects become increasingly important at elevated temperatures and are essential for accurately predicting longitudinal optical phonon features and infrared reflectance. This work establishes that accurate prediction of infrared optical properties in polar materials under thermal excitation requires explicit treatment of dynamic charge evolution.
The kagome lattice, renowned for hosting topological band structures and rich magnetic behaviors, offers an exceptional setting to investigate unconventional transport in magnetic topological systems. Controlling the polarity of the anomalous Nernst effect (ANE) is crucial for designing flexible thermoelectric devices, such as thermopiles, where the ability to switch the thermoelectric voltage sign can dramatically enhance energy conversion efficiency and output. Here, we demonstrate such a bipolar ANE in the kagome magnet YCo3, driven by a temperature-induced topological Lifshitz transition. With a Curie temperature TC~225 K, sizable anomalous Hall and Nernst effects emerge below TC. Supported by the first-principles calculations, the AHE and ANE are suggested to be dominated by the intrinsic mechanism. Furthermore, the intrinsic anomalous Hall conductivity exhibits a piecewise-linear dependence on magnetization, with an abrupt slope change near 100 K, consistent with the Karplus-Luttinger mechanism. Concurrently, the anomalous Nernst coefficient SAyx reverses its sign around the same temperature, realizing the crucial bipolarity. These anomalies could be interpreted as a topological Lifshitz transition, enabled by the evolution of Co moments that could shift the Fermi level relative to Weyl nodes. Our work reveals YCo3 as a prototypical kagome magnet where temperature and magnetism directly govern both Weyl node topology and the bipolar ANE, opening a pathway to magnetically control thermoelectric output in topological quantum materials.
Metal-organic frameworks (MOFs) are excellent candidates for water harvesting due to their tunable pore environments, which can be precisely engineered to capture and release water in arid conditions. Integrating artificial intelligence (AI) into MOF discovery can further accelerate the design of high-performance sorbents by identifying structural features that enhance atmospheric water harvesting (AWH), stability, and cycling efficiency. In this Perspective, we examine key MOF design principles, including cooperative adsorption, operational relative humidity (RH), uptake capacity, hysteresis, and scalability. We highlight recent design advancements such as multivariate strategies and long-arm linker extension, and examine how these principles tune pore capacity and hydrophilicity, while preserving stability and crystallinity. Furthermore, we discuss how AI, large language models (LLMs), and data mining can accelerate the discovery process through predictive synthesis, inverse design, and elucidating synthesis-structure-property relationships for the next generation of MOF water harvesters.
Memristive devices have gained significant attention for their potential in next-generation non-volatile memory and neuromorphic computing architectures. Among emerging candidates, transition metal oxides have proven particularly promising. While the switching mechanism in Ta/HfO$_2$/Pt devices was long attributed solely to oxygen vacancy based filaments, recent experimental evidence suggests a more complex dual-regime: the diffusion of metal cations also contributes to the formation of a conductive bridge. However, the precise atomistic mechanisms governing this metal cation migration remain poorly understood. Additionally, the role of defects such as oxygen vacancies present in the transition metal oxide in determining the final filament size and shape is also not well understood. Here, we employ molecular dynamics (MD) simulations with dynamic charge transfer to provide a detailed analysis of the atomistic mechanisms governing the co-formation of Ta-cation and oxygen-deficient filaments. We clearly show how varying the initial oxygen vacancy concentrations and spatial configurations within the HfO$_2$ matrix influences the final morphology and dimensions of the conductive filament. The switching is governed by field-driven formation and rupture of a hybrid Ta-cation-rich, oxygen-deficient filament in HfO$_2$. Our simulations closely match experiment, validating the model as a robust framework for understanding switching in oxide memristors and guiding designs that reduce cycle-to-cycle and device-to-device variability -- key barriers to high-performance devices.
Calcium batteries (CBs) are an attractive post-Li-ion technology, offering the appeal of Ca's natural abundance and high volumetric energy density. However, practical realization of CBs remains limited by the scarcity of positive electrode (cathode) materials that support reversible Ca$^{2+}$ (de)intercalation under electrochemical conditions. To address this challenge, we screen the materials project (MP) database for novel host structures that can intercalate Ca using geometry- and chemistry-based design principles. Specifically, we employ the Voronoi polyhedral volume as a descriptor of site compatibility for hosting Ca in potential frameworks. Further, we down-select candidate structures progressively through criteria including charge neutrality, absence of non-Ca mobile cations, thermodynamic (meta)stability, average voltage, and Ca migration barriers ($E_m$) using foundational machine-learning (ML) models. Subsequently, we validate the ensemble-ML-predicted $E_m$ in a subset of the final candidates using density functional theory based nudged elastic band calculations. Overall, from an initial pool of 52,945 MP structures, our workflow identifies 37 promising Ca cathode candidates, several of which exhibit favorable combinations of thermodynamic (meta)stability, voltage, and Ca-mobility, marking them as strong candidates for synthesis and electrochemical characterization. Particularly, we identify two Ca cathode candidates with markedly low Ca$^{2+}~E_m$ (CaSc$_2$V$_2$O$_8$ and CaVSO$_4$F$_3$), and four cathode candidates with thermodynamically stable charged states (Ca$_3$(CoO$_2$)$_4$, Ca$_3$Mn$_4$(TeO$_6$)$_2$, CaVF$_5$, and CaVSO$_4$F$_3$). Beyond identifying Ca-cathodes, our work establishes geometry-based descriptors and ML-based workflows as transferable methods for high-throughput screening, enabling the rapid discovery of novel materials for battery and other applications.
A Kondo insulator (KI) is a prototypical example of a highly entangled phase of matter, where many-body interactions between local moments and delocalized electrons engender the non-magnetic insulating ground state. Conventionally, the local moments arise from atomic multiplet states with a narrow bandwidth, limiting Kondo coherence to low temperatures. Here, we realize a new paradigm for constructing the KI state with hybridized molecular orbitals in FeSb2. Resonant inelastic X-ray scattering (RIXS) at the Fe L-edge reveals distinct signatures of band-like continuum states and localized states. Comparisons with first-principles calculations establish a mixed-configuration ground state with hybridized Fe d-Sb p molecular orbitals as basis states. By systematically investigating the RIXS momentum, temperature, and doping dependences, we find propagating collective modes commensurate with many-body charge and spin excitations. Our results pave the way for understanding the emerging class of unconventional d electron insulators and engineering high temperature Kondo many-body states.