Peer Revenue Fully funded

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Condensed Matter Materials Science

cond-mat.mtrl-sci Watch Materials Science (sign in to use)

Monday 17 August

Materials Science

Toward the Ultimate Limit: Elemental Metals in One Dimension

Low-dimensional materials exhibit extraordinary properties that make them promising candidates for advanced technologies. Although they have been investigated extensively, most of the research has focused on layered two-dimensional (2D) materials. Here, inspired by recent advances in atomically thin metallenes, we further reduce dimensionality and use density-functional theory simulations to study the geometry, energetics, elasticity, and electronic structure of 40 non-magnetic one-dimensional (1D) atomic chains of elemental metals. We find that nearly all chains have a buckled ground state, nine chains are distorted, and three chains—Cd, Hg, and Sr—are semiconducting with an electronic gap. We also find that transition metals retain a substantial fraction of their 3D bulk cohesive energy even in 1D chains. We assessed chains' dynamical stabilities by molecular dynamics simulations and found that 26 of them are thermodynamically stable at 100 K. Finally, we performed chain pulling simulations to investigate the straightening dynamics of selected stable chains. Given that experimental techniques have recently reached the 1D-chain limit, our systematic study provides a foundation and timely guide to accelerating synthesis and characterization of these materials.
Read the full abstract

All authors Mohammad Bagheri, Kameyab Raza Abidi, Sushree Sarita Sahoo, Sukhbir Singh, Pekka Koskinen

Materials Science

Efficient nonequilibrium electron dynamics from first-principles: leveraging Koopmans spectral functionals and Wannier localization

We present an efficient first-principles approach for simulating the nonequilibrium electron dynamics in extended systems beyond the linear regime. The method combines Koopmans-compliant functionals, which provide an accurate quasiparticle band structures, with the real-time evolution of the electronic density matrix in a Wannier basis within the Hartree plus screened exchange (HSEX) approximation. The locality of the orbital basis enables physically motivated approximations that significantly reduce both the computational cost and memory requirements while preserving accuracy. The screened Coulomb interaction, the central ingredient of the HSEX self-energy, is computed efficiently using density-functional perturbation theory. We benchmark the approach in the linear regime against experimental spectra and reference Green's function calculations for systems featuring both weakly and strongly bound excitons. Moving to the nonlinear regime, we investigate high-harmonic generation (HHG) in silicon and lithium fluoride. While in silicon the HHG spectrum is largely governed by the quasiparticle band structure, in LiF, a material featuring strong excitonics effect, the harmonic emission is selectively enhanced at excitonic resonances, suggesting that HHG probes correlated electron-hole excitations rather than solely the quasiparticle band structure. The present framework enables fully ab-initio simulations of excitonic effects in nonlinear optical spectra at a significantly reduced computational cost compared to real-time Green's function approaches, providing an efficient route to the study of ultrafast and strong-field phenomena in solids.
Read the full abstract

All authors Giovanni Cistaro, Miguel Sá, Davide Sangalli, Antonio Picón, Nicola Colonna

Materials Science

Spin-pump-induced spin transport demonstration in a photoconductive PTCDA molecular thin film with a transparent spin current detector

We demonstrate spin-pump-induced spin transport in a photoconductive PTCDA (3,4,9,10-perylene-teracarboxylic dianhydride) molecular thin film with a transparent ITO (In2O3 + SnO2) film as a spin current detector. In a tri-layer stacking structure sample composed of ITO/PTCDA/Ni80Fe20, pure spin current is generated in the PTCDA layer by using the spin pumping driven by the ferromagnetic resonance of the Ni80Fe20 layer. The generated spin current is absorbed into the ITO layer, converted to a charge current due to the inverse spin-Hall effect of the ITO layer, and detected as an electromotive force via the ITO resistance. Also, the light irradiation effect on spin transport in PTCDA films is investigated.
Read the full abstract

Materials Science · Cross-listed from cond-mat.soft

Finite strain homogenization of periodic rod networks with application to semi-flexible biopolymers

In this work, we adopt a finite strain computational homogenization approach to characterize the response of semi-flexible biopolymer networks modeled as idealized 8- and 14-chain periodic networks. We use the geometrically exact special Cosserat rod theory to model the microscale fibers forming these 8- and 14-chain networks. This allows us to capture arbitrarily large microscale deformations. Both macroscopic strain- and stress-driven homogenization are performed to study the macroscopic uniaxial tension, compression and simple shear responses. Several phenomena unique to biopolymer networks are recovered such as strain-stiffening and volume shrinkage under uniaxial tension, softening under compression and reverse Poynting effect under simple shear. We find that nonlinearity and non-affine deformation at microscale, especially bending and buckling of microscale fibers, plays an important role in these phenomena. We obtain the postbuckled solutions of the homogenization problem using a nonlinear, imperfection-free path following approach and also check for their stability. We further compare our homogenization results with experimental data for uniaxial tension and compression of biofilament networks and find good agreement. When the fibers are replaced by helical rods in the 8-chain unit cell, we are also able to capture the enlarged stretching behaviour as shown in recently fabricated compliant metastructures.
Read the full abstract

Materials Science

First-Principles Prediction of Phonon-Mediated Infrared Optical Properties of WO$_3$ Polymorphs

Many crystals exhibit polymorphism, undergoing atomic rearrangements that result in unit cells belonging to different symmetry groups. These structural changes directly affect lattice vibration modes and consequently influence their mid-infrared optical properties. In this study, we use tungsten trioxide (WO$_3$) as a representative system, owing to its multiple temperature-dependent phases, to study the impact of polymorphism on mid-infrared optical behaviour. Using a fully first-principles approach, we investigate three phases of WO$_3$ and evaluate their mid-infrared optical properties. Significant differences are observed among the three crystallographic phases, demonstrating the potential of this methodology as a predictive tool for materials discovery and targeted design.
Read the full abstract

Materials Science · Cross-listed from physics.optics

Transport-defined photodetection through electrically selectable nonequilibrium carrier transport

Broadband optical fields can change in both intensity and spectral distribution, but a fixed-response detector maps this evolving information onto a single electrical signal. Here we demonstrate transport-defined photodetection, in which electrical bias selects how photoexcited carriers are redistributed, escape and are collected, creating complementary response functions within one shared active region. In a GaAs/AlGaAs semiconductor ratchet, light-driven ratchet transport defines the response at 0 V, whereas the spectral evolution at -2 V is consistent with field-assisted hot-carrier transport. The overlapping states span a measured spectral range of 0.4-94.5 {\mu}m at 5 K, provide state-dependent calculated detection floors and support joint infrared operation at 30 K. Direct optical beat notes at 7.432 GHz in the mid-infrared and 17.103 GHz in the terahertz demonstrate optical-to-electrical conversion. Their common-path outputs recover an imposed spatial temperature gradient and the transient field-of-view-integrated effective radiation temperature of laser-excited graphite. These results establish post-photoexcitation transport as a function-defining design variable for semiconductor photodetectors, complementing structure-defined, field-tuned and optically encoded approaches to reconfigurable photodetection.
Read the full abstract

All authors Jiaxuan Cai, Yan Xie, Ning Yang, Weidong Chu, Peng Bai, Lianhe Li, Edmund H. Linfield, Hanbing Wang, Meng Chen, Guanchen Li, Yingxin Wang, Ziran Zhao

Materials Science

Beyond Local Berry Geometry: A First-Principles Finite-Momentum Theory of Electronic Position

Electronic position controls how a crystal polarizes and responds to an external field. In crystals, it is usually described through local changes of electronic states in momentum space. This Berry framework has reshaped modern solid-state physics, but strong fields drive electrons across a finite momentum range, where coherence between different momenta becomes part of the response. Here we establish a first-principles theory of electronic position at finite momentum that retains this missing information. We show that unequal-momentum coherence can cancel under spatial averaging and still produce polarization, forming a coherence dipole. We obtain the matrix directly from material wave functions, without model bands or fitted transition elements. In Si, the finite-momentum geometry sets a material momentum scale. Comparing this scale with the momentum change driven by the field predicts when finite-momentum physics becomes active. Crossing the scale strongly reorganizes the fifth and higher harmonics, showing that momentum-space geometry, rather than emitted photon energy alone, controls the nonlinear response. HHG is the first demonstration, but the theory applies whenever driven electrons explore a finite momentum range. It therefore extends quantum geometry beyond the local Berry limit and provides a general basis for predicting field-driven phenomena in real materials.
Read the full abstract

Sunday 16 August

Materials Science

Synthesizing like a chemist: an iterative, feedback-driven loop for materials discovery

Most computationally predicted materials are never synthesized because conventional synthesis optimization is slow, expertise-dependent, and iterative. Here we present a closed-loop framework that automates this expert workflow by placing human tacit knowledge in the loop through a large language model (LLM) that distills synthesis knowledge from the literature, high-throughput hyperspectral imaging for rapid film evaluation, and multi-objective Bayesian optimization guided by experimental feedback. In a paired optimization campaign, LLM-assisted initialization produced more Pareto-optimal samples and higher hypervolume than a Latin hypercube sampling baseline at matched trial counts, and this advantage persisted throughout iterative optimization. We demonstrate the framework by synthesizing the previously unreported perovskite-inspired compound Rb3BiI6 as thin films and validating the optimized films by optical bandgap analysis and X-ray diffraction. The framework transforms synthesis prediction from single-shot recommendation to iterative learning, providing a generalizable strategy to accelerate automated and fully autonomous experimental materials discovery.
Read the full abstract

All authors Fang Sheng, Steven B. Torrisi, Amanda Volk, Kevin Tran, Koki Nakano, Brian W. Anthony, Tonio Buonassisi

Materials Science

Crystal-structure design by agentic AI in a language of motifs

Data-driven materials discovery interpolates more reliably than it extrapolates and seldom reaches new structure types. We present MatEvolve, an agentic-AI framework designing crystals, proposing each candidate with a stated rationale and testing it. The agent reasons in an interpretable language of motifs, writing each crystal as a motif profile that describes the recurring geometric patterns—the motifs—composing it. The motif profile serves not merely as a description of a material but as the medium for material design: the agent edits the profile and constructs a crystal from the modified one, and the most promising candidates are validated by first-principles calculation. Applied to the design of rare-earth-lean permanent magnets, MatEvolve—built on the state-of-the-art language model Claude Fable~5 without fine-tuning—reaches new structural prototypes more than three times as often as generative models under an equal validation budget, at a comparable on-target-magnet rate. Beyond design, analysing the discovered crystals' human-readable profiles reveals structure–property relationships.
Read the full abstract

All authors Dinh-Khiet Le, Minh-Quyet Ha, Hong-Phuc Vu-Dinh, Takashi Miyake, Hiori Kino, Hieu-Chi Dam

Materials Science

Rare-earth chalcogenide perovskites: A promising class of materials for optoelectronic applications

Rare-earth chalcogenide perovskites have attracted significant attention for optoelectronic applications due to their nontoxic composition, robust phase stability, and excellent optoelectronic properties. However, their excitonic and polaronic properties remain largely unexplored due to the high computational cost of accurate theoretical treatments. In this work, we present a comprehensive first-principles investigation of excitonic dynamics and polaronic effects in a series of III-III rare-earth chalcogenide perovskites ABX$_{3}$ (A = Y, La; B = Sc, Y; X = S, Se), along with their structural stability and optoelectronic properties, using state-of-the-art density functional theory in conjunction with many-body perturbation theory within the G$_{0}$W$_{0}$ and Bethe-Salpeter equation (BSE) frameworks. All investigated compounds satisfy the dynamical and mechanical stability criteria. They exhibit quasiparticle band gaps in the range of 2.75$-$4.47 eV, and the BSE calculations reveal strong optical absorption spanning the visible to ultraviolet regions. The computed excitonic properties indicate intermediate-to-large exciton binding energies (0.148$-$0.517 eV), moderately localized excitons, and strong electron-hole wavefunction overlap, indicative of favorable radiative recombination characteristics and enhanced light-matter interaction. Furthermore, analysis based on the Fröhlich model demonstrates intermediate-to-strong carrier-phonon coupling, with electron-phonon interactions generally stronger than hole-phonon interactions. Overall, rare-earth chalcogenide perovskites ABX$_{3}$ exhibit a compelling combination of structural stability, tunable optoelectronic properties, pronounced excitonic effects, and favorable polaronic transport, positioning them as promising lead-free materials for next-generation optoelectronic devices, including light-emitting devices and photodetectors.
Read the full abstract

Materials Science

Strain-stabilized altermagnetism and conductivity anisotropy in FeSb2

We show that FeSb2 experiences a transition from a conventional antiferromagnet to an altermagnet when tensile strain is applied. In the altermagnetic phase, the lifted Kramers degeneracy results in spin splitting up to ~0.2eV near the Fermi level even without spin-orbit coupling. The transition to the altermagnetic phase is accompanied by a dramatic change in the Fermi-surface geometry, which leads to a uniaxial conductivity anisotropy up to ~60%, much greater than those observed in typical ferromagnetic metals. Using density-functional theory and Wannier interpolated Fermi surfaces, we show that this magnetotransport behavior may function as an experimental indicator of the transition to the altermagnetic phase. These findings highlight FeSb2 as a versatile, strain-tunable platform for exploring and utilizing altermagnetic transport phenomena for spintronic devices.
Read the full abstract

All authors Masoumeh Davoudiniya, Alyssa M. Kennedy, Amy Y. Liu, Gen Yin

Materials Science

ALKEMIE Agent: an autonomous platform for computational materials design

Despite the powerful multi-scale modeling methods and high-throughput infrastructures established in the materials community, real material computation workflows remain fragmented and heavily manual, requiring researchers to constantly bridge software tools, data analysis, and intermediate decisions. This growing gap between methodological capability and practical execution highlights the need for a new kind of autonomous computational framework, one that can coordinate tools, knowledge, and workflows in a more unified and adaptive way. Here, we introduce ALKEMIE Agent, an agentic platform in which retrieval-augmented generation, a materials-computation knowledge base, registered skills, database-supported provenance, AI-assisted structure modeling, bounded task execution, tool-calling iteration, and error-diagnostic assistance are integrated within a traceable control loop. The capabilities of ALKEMIE Agent are demonstrated through applications including materials recommendation, structure modeling, phonon calculations, machine-learned interatomic potential training, LAMMPS simulations, Ab Initio Monte Carlo (AIMC) sampling, and active-learning-based materials screening. Finally, we outline the future directions and challenges for the development of agentic platforms for computational materials design.
Read the full abstract

All authors Hongfu Huang, Yuzhe Li, Ao Xu, Bo Liu, Changrui Wang, Kan Tang, Ning Yang, Shengxian Liu, Hanyu Liu, Pengpeng Zhang, Linggang Zhu, Fengkai Liu, Yichen Lu, Tong Zhao, Naihua Miao, Jian Zhou, Zhimei Sun

Materials Science

Design of altermagnetism in oxide superlattices exploiting interface effects and quantum confinement

The discovery of altermagnetism has initiated intensive research and opened new avenues for spin- tronic and transport applications. While current efforts are mostly focused on bulk materials which are typically insulating, here we propose design strategies to achieve a combination of altermag- netism and metallicity in oxide superlattices by exploiting symmetry breaking, electrostatic doping and confinement. While bulk SrCrO3 does not exhibit altermagnetism due to compensating effects between adjacent layers our density functional theory calculations with a Hubbard U parameter re- veal, that a single SrCrO3 layer confined in a (SrCrO3)1/(SrTiO3)1(001) superlattice (SL) exhibits a sizable non-relativistic spin splitting (NRSS) up to 350 meV with bulk d-wave nature due to the coexistence of orbital ordering and octahedral rotations (OORs). Since this system is insulating, we extend to SrCrO3/LaCrO3(001) SLs. In the (SrCrO3)4/(LaCrO3)4(001) SL the combination of a polar discontinuity at the interface and stronger OORs promotes metallic d-wave altermagnetism. The NRSS of up to 120 meV is contributed by the interfacial Cr d bands at the Fermi level with indications for a spin-selective Fermi surface nesting. These findings establish oxide superlattices as a promising platform to realize and explore altermagnetism for quantum transport and spintronic functionalities
Read the full abstract

Materials Science · Cross-listed from cond-mat.supr-con

Superconducting Hydride Mg2RhH6 Experimentally Achieved at Lower Pressure

Although tremendous progress has been made in recent years in the field of polyhydride superconductors, the realization of high critical temperature superconductivity still relies on formidable high pressures. Searching for superconducting hydrides at lower pressures is of particular importance. Here we report the first experimental synthesis of the Mg2RhH6, which achieves superconductivity under a significantly reduced pressure of 30 GPa. The synthesis of Mg2RhH6 proceeds via a two step process (1) preparation of the Mg2RhH5 precursor containing hydrogen atoms stabilized by covalent bonds, followed by (2) hydrogen supplementation resulting in the filling of electrons into anti bonding orbitals above 30 GPa, which was accompanied by the structural transition from RhH5 square pyramid to RhH6 octahedron. Superconductivity is achieved at 30 GPa with a Tc of 24 K, which is further enhanced to 29 K at 53 GPa, evidenced by a sharp drop of resistivity to zero and characteristic suppression of Tc under applied magnetic fields. Our experiments prove the Mg2RhH6 superconductor to be thermodynamically stable above 30 GPa, making it the first case exhibiting a Tc of approximately 30 K at a readily accessible pressure. This study pioneers a highly promising pathway for the rational design and discovery of high temperature superconductors within phonon mediated BCS framework.
Read the full abstract

All authors Linjing Wu, Zelong Wang, Guiqi Liu, Jun Zhang, Yanfeng Ge, Yuanhao Su, Runteng Chen, Hongyu Liu, Wenmin Li, Sijia Zhang, Jingcheng Zhu, Jianfa Zhao, Zheng Deng, Shaomin Feng, Jing Song, Qingqing Liu, Xiang Li, Haozhe Liu, Panpan Kong, Xiancheng Wang, Changqing Jin

Materials Science

Uncovering the deformation mechanism of glasses during indentation through high-resolution X-ray scattering

Indentation experiments can be used to mimic real-life damage events of glasses that lead to surface flaws and thus lower practical strength. Conventional indentation studies often focus on the surface deformation after unloading. However, to understand the link between the surface deformation and structure, it is crucial to characterize the sub-surface deformation during the indentation process. The indentation-induced deformation, consisting of both elastic and plastic zones, is governed by the glass composition and structure, indentation and atmospheric conditions, and stress state. However, only a few experimental methods exist for characterizing the sub-surface indentation deformation mechanism during indentation. In this study, we use synchrotron X-ray nanoscattering to probe the deformation mechanism in situ during indentation of four types of oxide and oxynitride glasses with distinct structural features. This is done by measuring the variation in the position and intensity of the first sharp diffraction peak of the X-ray structure factor with a high spatial resolution down to ~100 nm. We find that the deformation zones of these glasses, which are characterized by the shape, size, and relative contribution between densification and shear flow under different indentation loads, vary with Poisson's ratio. Thus, our work provides new insights into the mechanical behavior of oxide glasses, contributing to the design of more damage-resistant glasses.
Read the full abstract

All authors M. Faizal Ussama Jalaludeen, Søren S. Sørensen, Johan F. S. Christensen, Anders K. R. Christensen, Sidsel Mulvad Johansen, Samraj Mollick, Yuanzheng Yue, Sharafat Ali, Sebastian Kalbfleisch, Morten M. Smedskjaer

Materials Science

Graph neural network prediction of temperature-dependent hydrogen diffusion and thermal conductivity tensors of tungsten containing helium bubbles and grain boundaries

Helium bubbles and grain boundaries in tungsten plasma-facing components alter hydrogen-isotope transport and thermal conduction by orders of magnitude, yet evaluating these transport properties for a given microstructure requires hours of molecular dynamics (MD) per configuration. We present a graph neural network surrogate that maps a tungsten atomic configuration containing helium bubbles and grain boundaries directly to the full $3\times3$ symmetric tensors of the hydrogen diffusion coefficient $D_H(T)$ and the thermal conductivity $\kappa(T)$ at arbitrary temperature. Anisotropy is captured by a rotation-equivariant tensor pooling layer; temperature enters through predicted temperature-independent parameters (an Arrhenius pair $(D_0,E_a)$, a phonon conductivity tensor, and a defect residual resistivity) expanded analytically via the Arrhenius and Wiedemann-Franz-Matthiessen relations. Training labels for 635 microstructures are generated with an embedded-atom-method potential (Green-Kubo conductivity and multi-temperature tracer diffusion); the electronic channel is calibrated against published irradiation-degradation measurements, and the pipeline is anchored to a first-principles machine-learning potential (VASP+FLARE) through paired MD calibration runs and an active-learning loop. The learned activation energies (median 0.21 eV, rising in bubble and grain-boundary structures) reproduce literature hydrogen migration barriers and trapping physics, and the equivariant pooling keeps predictions consistent across arbitrarily oriented sub-blocks. Coupled finite-element thermal-hydrogen analyses driven by the surrogate show that conductivity degradation changes predicted hydrogen permeation by a factor of 2.5 through the temperature field. The model returns both tensors in milliseconds, enabling microstructure-resolved transport input for component-scale analyses of fusion divertors.
Read the full abstract

Materials Science · Cross-listed from cond-mat.mes-hall

Observation of nonlocal ferron-drag thermoelectricity

The Peltier effect induces a heat current when a charge current passes through a conductor. Since a charge current is conserved at the junction between two different conductors, the difference between the heat flowing in both conductors for the same charge current leads to heating or cooling of the interface, providing an operating mechanism of solid-state heat pumps. Here, we report observation of heat absorption and release signals even in a junction-free, homogeneous metal when placed in proximity to a ferroelectric insulator. Our experiments using active thermographic imaging techniques confirm the prediction of the ferron-drag effect, i.e., the nonlocal excitation of ferrons, the collective excitation of the ferroelectric order, by conduction electrons in the adjacent metal. We reveal the electric-polarization-direction dependence of the temperature change signals and their unexpected increase with the metal thickness beyond the charge screening length, uncovering additional electron-phonon-ferron interactions in the metal/ferroelectric hybrid structure. The discovery of crosstalk between metals and ferroelectrics via remote ferrons could become both a nuisance and an opportunity for highly integrated circuits with ferroelectric barrier materials and revolutionize the design architecture of thermoelectric devices.
Read the full abstract

All authors Takuma Itoh, Takamasa Hirai, Ping Tang, Hossein Sepehri-Amin, Ryo Iguchi, Yusuke Kozuka, Takao Shimizu, Soshi Akita, Shunsuke Mori, Gerrit E. W. Bauer, Ken-ichi Uchida

Materials Science · Cross-listed from cond-mat.dis-nn

A General-purpose Solver of Fourier Neural Swarm Operator Towards Accurate and Efficient Mechanical Modeling of Ultra Large Composite Materials

Composite media with complex microstructures exhibit highly tailorable mechanical properties but remain challenging to model efficiently and accurately. Conventional homogenization often oversimplifies microstructural effects, whereas multiscale approaches typically require costly coupling across spatial and temporal scales. To address these limitations, we propose a two-scale neural-swarm framework for large-scale mechanical modeling of heterogeneous composites. At the local scale, the mechanical characteristics of representative microstructural features are encoded into building-block Fourier neural operators (FNOs) using level-set representations. At the global scale, these pretrained FNOs are assembled into an FNO swarm according to the spatial distribution of microstructural constituents. A coarse-mesh finite element model is employed to provide global physical guidance, while Schwarz iteration is used to synchronize neighboring FNOs and enforce consistency across shared interfaces. The proposed framework is validated through nonlinear simulations of SiC-Al composites with diverse microstructural configurations. Compared with nonlinear finite element analysis, the FNO-swarm method achieves comparable accuracy while reducing computational cost by orders of magnitude. For an extreme dual-property SiC-Al composite containing more than a billion nodal points, the proposed approach predicts the mechanical response within approximately one hour, demonstrating exceptional scalability. Furthermore, the framework naturally accommodates arbitrary Dirichlet boundary conditions and complex domain geometries. The proposed neural-swarm strategy provides a robust and scalable paradigm for large-scale mechanics simulations, reconciling the longstanding trade-off between computational efficiency and physical fidelity in heterogeneous materials modeling.
Read the full abstract

Materials Science

Intrinsic Wannier Functions for Hamiltonian downfolding

Downfolding ab initio material band structure into a low-energy subspace spanned by orbitals of specified atomic character, a procedure known as Wannier downfolding, is a common task in the simulation of complex materials. We introduce the Intrinsic Wannier Function (IWF) method to Wannierize bands with given atomic character. The method is non-iterative and requires only a single dimensionless parameter to disentangle bands. In benchmarks on silicon, graphene, and the three-band model of a mercury cuprate, we show that Intrinsic Wannier Functions provide high quality downfolded band structures compared to those from standard approaches such as Maximally Localized Wannier Functions and the Selected Columns of the Density Matrix method. Further, their straightforward implementation and robustness positions Intrinsic Wannier Functions as a general and useful tool for Wannier downfolding in materials electronic structure and in high-throughput applications.
Read the full abstract

Materials Science

Individual Vanadium Dopants Form Deep In-Gap States in Monolayer WS2

Point defects in atomically thin materials have a strong impact on physical properties and those that induce in-gap states are advantageous for quantum information science and engineering (QISE). However, dopant engineering consisting of well-controlled synthesis and robust identification of in-gap states is challenging. In this work, we addressed this challenge by first using finely tuned chemical vapor deposition to incorporate vanadium dopants into a monolayer WS2 (V-WS2). Next, we utilized a suite of scanned probe microscopy techniques to identify and characterize individual dopants. The latter included conductive atomic force microscopy (cAFM), low temperature scanning tunneling microscopy and spectroscopy (STM/STS), and scanning transmission electron microscopy and unambiguously revealed that vanadium dopants form deep in-gap states 0.35 eV above the valence band maximum in V-WS2. Our experimental results are well supported by first principles calculations and taken together demonstrate that V-WS2 is a promising platform for QISE applications.
Read the full abstract

All authors Tianhui Zhu, Carlos A. Gonzalez, Shihao Tu, Søren Tornøe, Ivan Pelayo, Dong-Rong Wu, Zhehao Ge, Hem Prasad Bhusal, Kenji Watanabe, Takashi Taniguchi, Nobuhiko P. Kobayashi, Yuan Ping, Jairo Velasco Jr., Aiming Yan

Materials Science

Hydrogen (deuterium) dynamics and thermal stability in ion-irradiated platinum-hydride thin films synthesized at low temperature

Hydrogen (H) and deuterium (D) interactions with transition metals play a central role in heterogeneous catalysis and hydrogen-related technologies. While H-Pt surface interactions have been extensively studied, direct investigations of hydrogen incorporation and transport in Pt remain limited due to its low solubility. Here, we study H(D) incorporation and desorption dynamics in metastable $PtH(D)_x$ thin films prepared by low-energy ion irradiation, enabling hydrogen loading far above equilibrium concentrations. Nuclear reaction analysis (NRA) reveals a nonuniform hydrogen depth profile with two accumulation regions: the subsurface and the film-substrate interface. Thermal desorption spectroscopy (TDS) exhibits two desorption peaks near 190 and 230 K, consistent with hydrogen release from these sites. Resistance relaxation measurements, analyzed within a two-parallel-channel conduction model, indicate different relaxation kinetics for subsurface and near-interface hydrogen. Arrhenius analysis reveals two thermally activated processes for $PtH_x$ with an average hydrogen concentration of $x = 0.15$, with activation energies of $130 \pm 18$ meV (subsurface) and $164 \pm 26$ meV (near interface). Above 140 K, D exhibits slower relaxation rates with activation energies of $117 \pm 8$ and $121 \pm 7$ meV for $PtD_x$ prepared under the same implantation dose. Within experimental uncertainty, the activation barriers remain comparable, while the prefactors are reduced significantly for D, indicating isotope-dependent attempt frequencies and zero-point energy effects. TDS simulations based on the Polanyi-Wigner formalism reproduce the experimental desorption spectra by resolving subsurface and near-interface contributions, in agreement with the NRA profile. These findings provide insight into hydrogen kinetics in $PtH_x$ for Pt-based catalysis, sensing, and hydrogen-metal interactions.
Read the full abstract

All authors S. S. Das, T. Ozawa, Y. Komatsu, R. Shimizu, T. Hitosugi, K. Fukutani

Saturday 15 August

Materials Science · Cross-listed from math-ph

Disentangling Measurement Protocol from Material Physics in Isotypic and Anisotypic SiC/Si Heterojunctions: A Metrology-Driven Study for Betavoltaic Energy Conversion

We compare isotypic (n-SiC/n-Si) and anisotypic (n-SiC/p-Si) heterojunctions using Kelvin probe force microscopy, scanning tunnelling spectroscopy, precision electrical measurements with a Keithley 2450 SourceMeter, and a 0.1 farad supercapacitor charge-storage test. Two-wire probing shows an enhanced low-bias current in the isotypic structures; four-wire Kelvin sensing removes the voltage-drop contribution of external leads and contacts and collapses this current to a much lower bulk and interface limited level. The high density of interface states detected by KPFM and STS, up to five times ten to the twelfth per square centimetre per electronvolt in dislocation-rich regions, provides a plausible origin for the enhanced local response. Anisotypic junctions exhibit a genuine diode-like current-voltage characteristic and an open-circuit voltage of about 1.2 volts under beta irradiation. In the capacitor test with the patented charge-pump converter, the anisotypic structure reaches 0.25 volts versus 0.05 volts for the isotypic one after 30 minutes, giving an energy ratio of about 25. We formulate a measurement-interpretation framework based on four-wire Kelvin sensing, generator polarity with cathode to LO, and high-impedance output-off state. Our results show that anisotypic SiC/Si heterojunctions are the preferred platform for betavoltaic conversion when measured correctly.
Read the full abstract

Materials Science · Cross-listed from physics.comp-ph

Scalable dynamical inference of phase-field fracture from sparse and partial measurements

Evolving crack fields in structural health monitoring and fracture assessment must often be inferred from sparse mechanical measurements rather than dense full-field observations. We develop CNN2D–ConvGRU, a convolutional-recurrent framework for measurement-conditioned reconstruction of time-dependent phase-field brittle fracture. At each load step, the model maps a fixed-length history of phase and displacement fields and sparse current-step displacement measurements to the current full-field state. New measurements are assimilated during sequential deployment, making the framework a state-inference surrogate rather than an autonomous time integrator. It reproduces crack paths, damage evolution, and bulk displacement response, with the largest errors near propagating crack tips and steep displacement gradients and some drift at late stages. Without retraining, weights learned on a $256 \times 256$ raster are evaluated on a $512 \times 512$ raster of the same physical domain and finite-element discretization using a proportionally refined measurement grid. This empirical raster-and-sensing transfer preserves the principal damage topology and global damage evolution, although fine-scale displacement errors increase near crack tips. Comparisons with alternative spatial and temporal architectures show that CNN2D–ConvGRU offers a favorable balance between reconstruction accuracy and computational cost. Relative to repeated finite-element solutions, sequential reconstruction achieves mean speedups of $175\times$ on CPU and $253\times$ on GPU. These results demonstrate efficient full-field fracture-state reconstruction from sparse observations while retaining the spatial structure and history dependence of phase-field fracture.
Read the full abstract

Materials Science

Intrinsic Defects in Amorphous Optical Coatings of TiO$_2$-doped GeO$_2$ for Gravitational-wave Detectors

The increased laser power of future gravitational-wave detectors will require mirror coatings with optical absorption below 0.1 ppm per mirror. TiO$_2$-doped GeO$_2$, currently the best high-index material for reducing room-temperature coating thermal noise, still exhibits ppm-level absorption even after extrinsic contamination is minimized. Using ab-initio simulations and absorption measurements, we identify oxygen-deficient Ti-rich environments as the origin of this residual absorption. We find that ordinary structural disorder in the amorphous network can localize electronic states but does not produce defects capable of absorbing 1064-nm light. In contrast, oxygen vacancies in compact Ti-rich environments create localized Ti$^{3+}$–Ti$^{3+}$-like polaron-pair or mixed Ti-polaron states states with transitions near 1064 nm. Photothermal measurements show increased absorption after dry/inert annealing, supporting the formation of these reduction-sensitive defects. These results show that the residual absorption is not an intrinsic limitation of TiO$_2$-doped GeO$_2$, but a process-dependent defect that may be mitigated through control of oxygen stoichiometry during deposition and annealing.
Read the full abstract

All authors K. Prasai, K. Lee, R. Bassiri, A. Davenport, D. A. Drabold, M. M. Fejer, A. Markosyan, C.S. Menoni, S. Tait

Materials Science

What makes a useful molecular model of biochar? A community roadmap

Biochars are disordered carbonaceous materials produced by biomass pyrolysis, with applications spanning soil amendment, water remediation, carbon storage, and functional materials. Although they share structural features with other disordered carbons such as coal, kerogen, and activated carbons, the questions posed to biochar models are distinct, and no single model can answer all of them equally well. Model usefulness must be defined relative to a specific question and validated against independent experimental observables. This community roadmap, arising from a CECAM workshop, critically maps current molecular approaches: experimentally guided top-down reconstruction, mimetic bottom-up simulation, and hybrid methods. We argue that first-generation models have been more successful than is often acknowledged, provided they are built at sufficient length scale and with explicit control over microporosity and bulk chemistry. Structural and equilibrium interfacial properties are increasingly tractable with classical force fields, whereas dynamic and reactive behaviours require selective use of reactive methods within multiscale workflows. A parallel, largely unaddressed gap concerns the mineral and ash components of biochar, and the changes the material undergoes during ageing in soil. We identify seven open questions current models cannot yet answer reliably, and five community priorities: force field benchmarking, open model and data repositories, shared classification and metadata standards, ensemble validation, and training in reproducible practice. Across these, sustained interaction with experimentalists is essential to ground models in real observables and document where they fail. Progress will be accelerated by adapting transferable methods from coal, kerogen, and clay-organic matter frameworks rather than repeating trial-and-error development.
Read the full abstract

All authors Valentina Sierra-Jimenez, Jonathan P. Mathews, Luca Bellucci, Edo Boek, Carla de Tomas, Manuel Garcia-Perez, Stef Ghysels, Paola Giudicianni, Corinna Maria Grottola, Kelly Anne Hawboldt, Robert L. Johnson, Fenna B.E. Kolff, Jean-Marc Leyssale, Diego Liberati, Francisco J. Martin-Martinez, Jacob W. Martin, Ondřej Mašek, Mohammad Mezbah Ul Hoque, Audrey Ngambia, Amaël Obliger, Frederik Ossler, Muhammad Riaz, John M. Tobin, Xiaolei Zhang, Valentina Erastova

Materials Science

Barium Hexaferrite Thin Films as a Scalable Magnetic-Insulator Platform for Proximity-Engineered Spintronics

Rare-earth iron garnets, such as yttrium iron garnet (YIG) and thulium iron garnet (TmIG), are the benchmark magnetic insulators for spintronic and magnonic devices, but achieving usable perpendicular magnetic anisotropy (PMA) in these materials typically relies on substrate strain- engineering, requiring careful lattice-matching and specific growth conditions that constrain ma- terial accessibility. Here we establish sputter grown barium hexaferrite (BaFe12O19, BaM) as a magnetic-insulator alternative with strong intrinsic perpendicular anisotropy, requiring no strain engineering. X-ray diffraction, transmission electron microscopy and Raman spectroscopy confirm stoichiometric films with atomically smooth surfaces, while first-principles calculations corroborate a robust ferrimagnetic ground state. The films exhibit square out-of-plane hysteresis with a coercive field of nearly 0.1 T. Unlike rare-earth garnets, the perpendicular anisotropy in BaM is intrinsic to its magnetoplumbite crystal structure, arising independent of highly ordered strain. Interfaced with Pt and with exfoliated BiSbTeSe2 (BSTS), BaM induces proximity induced anomalous Hall trans- port, confirming efficient interfacial exchange coupling, while the BSTS/BaM heterostructure shows an additional Hall contribution suggestive of non-collinear interfacial spin textures. These results position BaM thin films as a scalable magnetic-insulator platform for spintronic and topological heterostructure devices beyond the constraints of garnet chemistry.
Read the full abstract

All authors Shyam Sundar Poriah, Sanjana D. S., Agrim Sharma, Sreelakshmi M. Nair, Pankaj Bhardwaj, Laxmipriya Nanda, Aryaman Das, Jagadish Rajendran, R. S. Patel, Manish Jain, Dhavala Suri

Materials Science

Claude-SpinDynamics: a cross-platform, dual-precision CPU/GPU micromagnetic simulator with native mumax3 script compatibility

Quantitative spintronics increasingly depends on a handful of GPU micromagnetic codes, all of which require NVIDIA hardware and single precision throughout, leaving researchers without such hardware unable to run even the standard validation problems. We report Claude-SpinDynamics (Claude-SD), a new open source micromagnetic simulator with a cross platform C++20 core (Windows and Linux) and a Python interface that closes this gap: a complete CPU build, validated by the same test suite as the GPU path, runs every unit test and uMAG standard problem with no accelerator at all, alongside GPU builds offering both single and double precision, a choice of two demagnetization FFT backends, and natively implemented spin-orbit, spin-transfer, and Zhang-Li torques, Dzyaloshinskii-Moriya interaction, and percell materials.Claude-SD natively interprets mumax3's .mx3 scripting language, so existing community scripts run unmodified; under matched conditions the two codes agree cell-by cell to single-precision round-off, and, together with mumax+ and OOMMF, to within 2% on the uMAG dynamic-switching standard problem, with MuMax-CO agreeing to mumax3 to float32 round-off on the same problem. Benchmarked head-to-head against these three codes, Claude-SD's single-precision build is the fastest solver on small and two-dimensional problems and remains competitive at the largest grid sizes, while its double-precision and dual-FFT-backend paths are unmatched among GPU micromagnetic codes. A GPU replicabatching extension further advances an entire ensemble of finite-temperature trajectories in a single kernel launch per step, giving one to two orders of magnitude of throughput over a per-trial loop while reproducing single-trajectory results to numerical round-off. The complete source is openly licensed and distributed with runnable example notebooks and documentation for independent reproduction.
Read the full abstract

Materials Science

Linking Electronic Bonding and Short-range Order to Strength in $\alpha$-Titanium Alloys: A First-Principles Study

The development of accurate strength prediction models for titanium alloys is critical for advanced materials design. This study systematically examines how the mechanical properties of $\alpha$-Ti are affected by substitutional (X = Al, V, Mo) and interstitial (Y = H, C, N, O) alloying elements, with a focus on electronic bonding. Using density functional theory (DFT), we uncover the short-range ordering (SRO) of substitutional atoms and quantify their influence on the electronic bonding and mechanical behavior. The primary novelty of this work lies in developing a predictive model for tensile strength that goes beyond traditional empirical approaches. To quantify the contributions of individual solutes to strengthening, we use physically grounded quantum-chemical descriptors, such as the Integrated Crystal Orbital Hamilton Population (ICOHP), which is a direct measure of bond strength derived from first-principles calculations. The resulting formula quantitatively predicts the tensile strength of a wide range of $\alpha$-Ti alloys, demonstrating a significant advancement in the computational design of high-performance structural materials.
Read the full abstract

Materials Science · Cross-listed from physics.optics

Driven Time Crystal in Low-Symmetry ENZ Conductors

In recent years, epsilon-near-zero (ENZ) materials have attracted a great deal of attention in nonlinear optics, as they combine field enhancement with strong, ultrafast nonlinearities. In particular, transparent conducting oxides, such as ITO, have emerged as a promising class of materials and have been extensively exploited to achieve temporal optical responses varying on the femtosecond scale using an optical pump. Most of the solutions discussed so far in the literature rely on effective $\chi^{(3)}$ modulations, wherein the dominant material response is controlled by the envelope of the optical pump. Here, it is shown that low-symmetry conductors can provide an interesting alternative to transparent conducting oxides and a more natural implementation of time-crystalline behavior in optical systems with optical-cycle modulation. I demonstrate that ENZ confinement, combined with the strong anomalous-velocity nonlinearity of low-symmetry conductors, enables a subwavelength nanoparticle to develop a time-crystalline response under optical pumping. For sufficiently strong pumping, this response can overcome dissipative losses and lead to parametric amplification. Furthermore, the pump can strongly tailor the scattering and extinction of a weak probe and, in extreme cases, render the extinction negative. In this regime, the driven nanoparticle effectively amplifies the probe beam.
Read the full abstract

Materials Science

Broadband phonon-velocity suppression and a finite anisotropic crossover in twisted bilayer SnSe

Moiré superlattices reshape lattice dynamics without altering chemical composition, yet how crystal anisotropy modifies this control remains unclear. We combine density-functional-theory (DFT)-calibrated lattice-dynamical calculations with angle-matched untwisted controls to study puckered bilayer SnSe across seven commensurate twist angles ($3.18^\circ$–$8.77^\circ$). At 300 K, twisting suppresses the band-path heat-capacity-weighted mean-square group velocity to 2.6–8.4% of the control values; the suppression spans a broad frequency range rather than a few soft branches. The velocity response crosses over between $4.78^\circ$ and $3.82^\circ$ into a regime where the relaxed stacking textures and frequency-resolved velocity profiles become self-similar, with the normalized mean-square velocity ratio spanning only 11.1% of its mean across the three smallest angles—a finite anisotropic crossover, not a singular-angle condition. Direct DFT–MACE force-constant agreement ($r=0.996$), uniform $4\times4\times1$ stability scans, and acoustic-sum-rule and path-density tests support the trend. The equilibrium trend is defined by six structures after excluding one relaxation-sensitive case. These results extend phonon twistronics to low-symmetry layered materials and identify crystal anisotropy as a key determinant of finite-angle phonon crossover behavior.
Read the full abstract

All authors Peng Kang, Wei Yin, Da Wan, Shulin Bai, Sirui Fan, Qi Zou, Hongfeng Li, Xiao Xiang, Zhen Li, Yu Liu, Lei Zheng, Li-Dong Zhao