Peer Revenue Publish or profit

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)

Saturday 26 September

Materials Science · Cross-listed from cond-mat.str-el

GEM: An implementation of the ghost-Gutzwiller approximation for simulating interacting quantum systems

We present GEM (Ghost Embedding Method), an open-source software package written in Python for computing equilibrium properties of strongly correlated electronic systems within the ghost-Gutzwiller approximation method. GEM provides a computationally efficient framework for studying multi-orbital lattice models. It supports zero- and finite-temperature calculations and symmetry broken phases. It is integrated with the TRIQS ecosystem, providing tools for model construction, self-consistent solution, and evaluation of physical observables. We first detail the method's theoretical formulation, then we present the software architecture, and finally we introduce some practical workflow, which also validates the implementation against established results. In particular, we illustrate the capabilities of GEM through multiorbital and finite-temperature applications and discuss its computational cost relative to more demanding quantum embedding approaches.
Read the full abstract

All authors Samuele Giuli, Tsung-Han Lee, Yong-Xin Yao, Ina Park, Harrison LaBollita, Ivan Pasqua, Nils Wentzell, Nicola Lanatà, Olivier Gingras

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

Phonon angular momentum transfer torque

Angular momentum in solids is carried by (quasi)particles such as magnons, plasmons, and phonons, and coupling between these reservoirs enables diverse hybrid phenomena. Here, we introduce the transfer of angular momentum from lattice vibrations to magnetic order through the phonon angular momentum transfer torque (PAMTT), generated by a polarized phonon bath. We identify thermal routes for generating phonon angular momentum: a phonon Edelstein accumulation in noncentrosymmetric crystals and temperature-gradient-induced PAM currents with longitudinal Seebeck-like and transverse Hall components. Using realistic order-of-magnitude estimates, we find that PAMTT could produce measurable ferromagnetic-resonance frequency shifts and, for favorable interfacial coupling, drive magnetization reversal on nanosecond timescales.
Read the full abstract

All authors Verena Brehm, Daniel A. Bustamante Lopez, Shu Zhang, Dominik Juraschek

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

Electrically switchable one-dimensional quadrupolar excitons in lateral double heterojunctions

Two neighboring lateral interfaces provide a spatial degree of freedom for controlling one-dimensional charge-transfer excitons within a single semiconductor monolayer. We investigate a type-II WS2-MoS2-WS2 double heterojunction using an effective-mass two-particle Hamiltonian with a screened Coulomb interaction. For equivalent left and right interfaces at zero electric field,inter-interface coupling produces energetically split even- and odd-parity exciton states, each with zero permanent dipole. An electric field perpendicular to the interfaces continuously converts the lower state from a quadrupolar superposition with a quadratic Stark shift into a predominantly single-interface dipolar exciton with an approximately linear shift. The spatially resolved calculation gives binding energies of approximately 104 and 100 meV, a doublet splitting of 4.6 meV, and a crossover field of 0.87 V/um for a representative 1.5-nm MoS2 strip. Projection tests show that the lowest doublet controls the response near this crossover. Strip width tunes the coupling much more strongly than the binding energy, providing geometric control of the low-field Stark sensitivity. These results establish a continuum-model route to electrically reconfigurable one-dimensional quadrupolar excitons.
Read the full abstract

Materials Science · Cross-listed from cs.LG

Phase Space Attention:A Hairer Lift Circumvents the Single-Layer Induction Obstruction

We circumvent the Sanford-Hsu-Telgarsky (SHT) single-layer induction obstruction within a linear, one-step, causal, bilinear, symplectically consistent design class on the post-RoPE substrate, by lifting attention onto a symplectic phase space, mirroring Hairer's lift of Stormer-Verlet. The lift exits the premise of the SHT counting argument rather than the bound itself. The reframing exhibits the obstruction as a filter-order gap: a one-layer bilinear score realises a $z$-transform of joint order $(0,0)$, whereas the induction discriminator requires key-side order $\geq 1$. Applying the symplectic upper shear $M_\gamma:(q,p)\mapsto(q+\gamma p,p)$ to the post-RoPE query and key streams closes it. We prove this lift is unique within the factorised subclass, exactly symplectic at operator level, and requires post-RoPE placement; and in an explicit $T_4$-only Gaussian reduction we derive a closed-form two-branch induction phase transition, held out at $r=0.9876$ with zero fitted parameters. That law is an analytically solvable limit, not a robust prediction: restoring the $T_3$ channel moves $\gamma_c$ at $d_k=64$ from 1.030 to 0.569 and removes the crossover. Deployability follows by exact derivation: the KV cache is unchanged, prefix reuse and speculative decoding are preserved, overhead is $6d$ FLOPs per token per layer, INT8 headroom grows by at most $\log_2(1+2\gamma)$ bits, fused kernels are unmodified, and no parameters are added. At 91.3M parameters a supercritical sweep locates an emergence band: induction forms 3/3 seeds at $\gamma=0.80$ in a mean of 717 steps, against 2/3 seeds and 2700 steps at $\gamma=0$. Adverse results are reported as directly: a key-only half-lift reaches 0.949 against 0.811 for the symmetric operator, so if induction accuracy is the objective, the half-lift is the better construction. Forty-one notebooks and result files ship as ancillary material.
Read the full abstract

All authors Kingsuk Maitra, Shagun Sood Morteza Hosseini, Suman Gunnala, Vikram Gupta

Materials Science

Dopant-modulated lattice softening drives drastic thermal conductivity reduction in \beta-FeSi2 thermoelectrics

Suppressing lattice thermal conductivity (K_lat) is pivotal for thermoelectric efficiency. While traditional strategies rely heavily on phonon scattering from mass- and size-mismatches, we demonstrate a robust K_lat suppression mechanism driven by dopant-induced lattice stiffness modulation. Through a comparative analysis of p-type (Mn) and n-type (Co, Ir) doping in the \beta-FeSi2 model system, we show that Co and Ir doping significantly reduce K_lat. Notably, Co doping achieves a ~71% reduction at 300 K even without significant mass and size contrast. By correlating transport data with neutron powder diffraction, heat capacity, and Raman spectroscopy, we reveal anomalous lattice expansion, a substantial reduction in Debye temperature, and marked vibrational redshift and broadening. These systematic changes provide strong evidence for atomic-scale lattice softening and a fundamental weakening of interatomic force constants, which synergistically lower phonon group velocities and amplify anharmonic scattering. Our findings establish lattice stiffness manipulation as a powerful strategy for thermal management, offering a distinct design pathway beyond traditional mass- and strain-fluctuation models.
Read the full abstract

All authors Cuiping Zhang, Qingyong Ren, Yangfan Cui, Chen Chen, Songbai Hu, Shengnan Dai, Chin-Wei Wang, Wanju Luo, Dexiang Gao, Bao Yuan, Junying Shen, Fan Chen, Wei Xu, Yuting Li, Mingfang Shu, Xiaoli Huang, Pengfei Qiu, Jie Ma

Materials Science

Temperature Dependent Evolution of the Electronic Structure in EuZn2As2 across the Neel Transition

Magnetoresistive materials have been tremendously important for the development of magnetic memory storage and spintronic devices. Recently, the antiferromagnetic EuX2Pn2 compounds, with X being a transition metal and Pn being a pnictogen, have seen intensive research interest due to their unusual anomalous Hall effect behavior and pronounced resistive anomaly near the Neel temperature (TN). These magnetotransport phenomena have been interpreted in the context of shortranged ferromagnetic Structuations, magnetic polaron formation, canted spin configurations, and temperature dependent metal insulator transitions in the electronic structure. Here, we report the observation of such a pronounced resistivity anomaly in EuZn2As2 near TN = 19 K. We demonstrate the suppression of this anomaly using applied magnetic fields, both in plane and out of plane. To further interpret the origin of the observed transport behavior, we studied the temperature dependent electronic structure using combined angle resolved photoemission spectroscopy (ARPES) and first principles density functional theory (DFT) calculations, which exhibits limited modifications to the bands across TN away from the Fermi energy. This lack of involvement of the electronic structure indicates a spin-scattering origin of the aforementioned transport properties, rather than a reconstruction of the Fermi surface.
Read the full abstract

All authors Milo Sprague, Anup Pradhan Sakhya, Barun Ghosh, Mazharul Islam Mondal, Arun K. Kumay, Himanshu Sheokand, Kapil Gope, Tetiana Romanova, Dariusz Kaczorowski, Arun Bansil, Madhab Neupane

Materials Science · Cross-listed from cs.RO

Active 3D weaves for load-bearing and damage-resilient locomotion

The craft of weaving, where different materials are interlaced, has tremendous potential for creating active and functional systems for use in soft robots, prosthetics, wearables, exoskeletons, and more. These textile-like systems are flexible and safe for human-machine interaction; however, this inherent flexibility limits their ability to carry loads, which is essential for many robotic functions. In this work, we introduce a general framework for integrating active materials into three-dimensional (3D) woven shells to create robotic structures that combine high axial stiffness for load bearing, low bending stiffness for efficient actuation, and system-level resilience for damage tolerance. These woven robots can be modularly assembled from 'woven corners', a fundamental unit of 3D woven structures. We use eigenvalue calculations to identify load bearing and actuation mechanisms of the 3D woven structures, and use that information to make five different robots capable of locomotion. We demonstrate that these 3D woven robots can locomote carrying loads 70 times their self-weight, and can maintain repeatable performance even after being subjected to extreme compression. This work is a pathway toward the design, manufacturing, and simulation of future 3D woven robotic systems where load bearing, high stiffness, active functional deformation, locomotion, and system-level resilience are all needed.
Read the full abstract

Friday 25 September

Materials Science

Jahn-Teller Distortion-Driven Lattice Parameter Tunability in Coherent High-Entropy Oxide Thin Films

Here, we isolate how Cu concentration and Jahn-Teller distortions control the out-of-plane lattice parameter of pseudomorphic MgCoNiCuZnO-derived high-entropy oxide films on MgO. We prepare nine compositions spanning 0-20% Cu, verify target composition and phase by X-ray fluorescence and diffraction, and grow films by pulsed laser deposition across multiple substrate temperatures. At 500 {\deg}C, the lattice parameter increases nearly linearly with Cu content; at 300{\deg}C, it decreases nonlinearly, producing a 2.6% difference between temperature regimes for equimolar MgCoNiCuZnO at our growth conditions. We attribute this change in lattice parameter to temperature-dependent coupling between Cu-driven Jahn-Teller distortions and Co valence change. While prior work identifies Co valence as the primary driver of the MgCoNiCuZnO lattice response with temperature; our results establish Cu Jahn-Teller distortion as an equally important structural control parameter.
Read the full abstract

All authors Dhiya Srikanth, Joseph Petruska, Matthew Furst, Kaylin Lamaute, Jon-Paul Maria, Saeed S. I. Almishal

Materials Science

The Young-E(3) Tensor Product Decomposition for Rotation and Permutation Equivariant Cluster Expansions

A generalization of fixed-lattice cluster expansions (CE) and atomic cluster expansions (ACE) is presented that expresses both rotation and permutation symmetries. By performing Schur–Weyl decomposition in Young subgroup-stabilized carriers, followed by joint coupling of rotation and permutation representations, arbitrary $S_N\times SO(3)$ tensor product carriers are generated. This allows us to resolve complete, orthonormal joint rotation and permutation-adapted bases for arbitrary tensor product ranks and arbitrary angular and radial tensor product content. We show that this very general Young–E(3) (YE3T) tensor product basis contains the ACE basis as a subset, corresponding to the special case where permutation symmetry character is restricted to the fully symmetric carrier. Rather than constructing an overcomplete rotation and permutation-invariant basis and reducing it a posteriori, Barthelemy et al. recently demonstrated scaling benefits by not constructing an overcomplete basis. YE3T directly produces a complete orthonormal basis without an overcomplete step and rigorously extends to permutation characters beyond permutation-symmetric features. The YE3T decomposition yields new joint irreducible rotation- and permutation-equivariant basis sets that surpass existing rotation-adapted expansions in both speed and accuracy, defining a new Pareto front for machine-learned interatomic potentials. We show that the tunability of the permutation symmetry character makes the basis useful for both atomistic and electronic-structure simulations.
Read the full abstract

Materials Science

Knowledge-Driven XRD Phase Identification via Multi-View Retrieval and Explanation

X-ray diffraction (XRD) is a experimental technique for determining the phase composition and structure of crystalline materials. However, interpreting XRD patterns is challenging, particularly in high-throughput materials discovery, where many novel materials may need to be characterized and no reference patterns are available. Consequently, machine learning is increasingly used to accelerate and automate the analysis while reducing errors associated with human interpretation. We propose a multi-decision framework for XRD phase analysis that integrates representation learning, similarity-based retrieval, and explainable decision support within a unified reference database. A convolutional autoencoder learns compact latent representations of XRD patterns that preserve structural similarity while remaining robust to variations arising from experimental noise and measurement conditions. By integrating multiple decision pathways within a shared latent space, the framework moves beyond single-label prediction toward ranked and interpretable phase analysis that mirrors expert practice. During inference, complementary decision mechanisms are applied, including latent-space classification and retrieval, explanation-guided similarity using Integrated Gradients, and composition-based similarity search. These mechanisms generate ranked candidate phase lists that are aggregated into a final prediction with an associated confidence score. Experiments on synthetic datasets demonstrate strong predictive performance, achieving 98.85\,% accuracy for crystal system classification and 95.82\,% accuracy for space group prediction on the test set, while maintaining robustness under realistic perturbations. The framework supports reliable, analyst-friendly identification of crystal phases and structures in high-throughput and exploratory materials discovery settings.
Read the full abstract

Materials Science · Cross-listed from cs.AI

CP-Agent: A Harness-Engineered Agent for Crystal Plasticity Simulation Workflows

Crystal plasticity (CP) simulations predict the mechanical behavior of polycrystalline metals, yet their routine use is hindered by the manual effort of configuring heterogeneous tools, orchestrating multi-step data pipelines, and calibrating constitutive parameters against experiments. These bottlenecks impede productivity in systematic parameter studies, motivating interest in automated workflows. This study presents CP-Agent, a harness-engineered LLM-based agent that autonomously executes complete CP modeling workflows from natural-language tasks. Operating under the ReAct paradigm, the agent reasons about tool selection and sequencing while delegating numerical search to established optimizers. The harness comprises a minimal system prompt, typed tool definitions, a dispatcher, and a safety-bounded iteration loop, encoding domain knowledge through tool schemas rather than hard-coded logic. CP-Agent is demonstrated on four case studies: calibrating four slip parameters of additively manufactured stainless steel 316L against tensile data; validating the workflow against published copper benchmarks, reproducing stress-strain and texture evolution; recovering the initial crystallographic texture of copper, where the agent correctly identifies a diffuse initial texture; and reproducing the multi-pass rolling texture evolution of a Mg-Zn-Ca alloy, where the agent chains five deformation passes and recovers the experimentally observed weakened, split basal texture. In all cases, the agent inferred the correct execution sequence from the task statement, robustly across repeated runs, and delivered physically interpretable results. This work establishes harness engineering as a systematic approach to automating CP modeling workflows while maintaining physical interpretability and auditability through visible reasoning traces.
Read the full abstract

Materials Science · Cross-listed from cond-mat.str-el

A Unified Spin-Fermion Framework for Magnetic Diversity in Chromium Monopnictides

Isoelectronic compounds are generally expected to exhibit related electronic, structural, and magnetic behavior. Chromium monopnictides CrX (X= Sb, As, P), however, display a striking evolution from high-temperature altermagnetism in NiAs-type CrSb, through double-helical antiferromagnetism in MnP-type CrAs, to the absence of resolved long-range magnetic order in MnP-type CrP. Here, combining first-principles calculations in the non-collinear Curie-paramagnetic state, structural analysis, magnetic phase-diagram calculations, and exchange-parameter extraction, we establish a unified microscopic framework for this evolution. All three compounds share a formal high-spin $d^5$ configuration, nominally Cr$^{1+}$, coupled to itinerant states of predominantly pnictogen character. From Sb to P, chemical pressure enhances ligand itinerancy and drives the structural evolution from the high-symmetry NiAs lattice to increasingly distorted MnP-type networks. The accompanying reconstruction of competing Cr-Cr exchanges favors A-type antiferromagnetism in CrSb and double-helical order in CrAs, while placing CrP near a frustrated regime with closely competing magnetic tendencies. Our results establish chemical-pressure-driven exchange reconstruction as a general mechanism through which a common spin-fermion electronic structure can generate contrasting magnetic phases in strongly correlated materials.
Read the full abstract

Materials Science

Unconventional magnetoelastic behavior in Al rich CoFeAl Films with inverse Heusler like local order for flexible spintronics

Strain engineering of magnetic properties offers a promising route toward flexible spintronic applications. Here, we report an unconventional magnetoelastic response in Al rich Co Fe Al thin films on flexible substrates using strain dependent magneto optical Kerr effect microscopy and magnetometry measurements. While the films exhibit a conventional positive magnetostriction coefficient, consistent with standard in plane easy axis rotation under stress, their saturation magnetization increases under compressive strain and decreases under tensile strain. First principles calculations reveal that this unconventional response originates from a strain induced competition between exchange splitting and crystal field effects in an inverse Heusler like local environment created by Al enrichment. This leads to a highly sensitive, sublattice dependent magnetic state, consistent with a strain induced reconfiguration of Co and Fe moments. Our results demonstrate that local compositional tuning can fundamentally alter magnetoelastic behavior, establishing strain controlled sublattice compensation as a route toward programmable magnetic functionality in flexible spintronic systems.
Read the full abstract

All authors Rupalipriyadarsini Chhatoi, Anuroopa Behatha, Swayang Priya Mahanta, Shubhransu Sahoo, Soubhagya Dash, Bhuvneshwari Sharma, Abhisek Mishra Esita Pandey, Satadeep Bhattacharjee, Subhankar Bedanta

Materials Science

Structural prediction of B$_{18}$Y$_{2}$ cluster: A Machine-Learning-Assisted Basin-Hopping Study

The structural and optical properties of the doubly yttrium-doped boron cluster B$_{18}$Y$_2$ have been systematically investigated using density functional theory calculations. The lowest-energy structure was identified through extensive basin-hopping searches accelerated by a pre-trained MACE machine-learning potential and subsequently refined and validated at the DFT level. The resulting global-minimum structure adopts a double-ring geometry, consisting of two fused B$_9$ rings stabilized by yttrium atoms positioned above and below the boron framework. Vibrational frequency calculations confirm the dynamical stability of the optimized structure, while the calculated infrared and UV–Vis spectra provide characteristic signatures of the Y–B interactions and the electronic structure of the boron framework. These results demonstrate how yttrium doping can stabilize unusual double-ring boron architectures and highlight the effectiveness of machine-learning-assisted basin-hopping searches for exploring the complex potential-energy landscapes of doped boron clusters.
Read the full abstract

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

Fabry-Pérot Resonance Shifts as a Novel Non-destructive In Situ Surface Damage Diagnostic

Current surface damage evaluation for nanostructured electronic devices requires specialized equipment or conditions outside the designed operating environment and cannot be used for in situ evaluation of surface damage. In this work, we developed a non-destructive method to probe the surface sublimation of a semiconductor heterostructure which can be applied to other semiconductor damage mechanisms. This approach employs shifts in the resonances of an optical cavity comprising the surface and an integrated distributed Bragg reflector fabricated below the semiconductor heterostructure undergoing sublimation. We performed the demonstration of this technique under vacuum, but surface damage effects can be observed in any operational environment either in situ or ex situ.
Read the full abstract

All authors Gregory Blume, Colin Kirk, Sushil Poudel, John Hill, Alimohammed Kachwala, Joseph Grames, Marcy L. Stutzman, Matt Grau, Sylvain Marsillac

Materials Science

Retrainable physics-integrated neural differentiable modeling of sintering across material systems

Sintering is widely used to manufacture ceramics, but coupled densification and grain growth, material-dependent kinetics, and sparse measurements complicate predictive modeling and process design. We present Sinter-PiNDiff, a retrainable physics-integrated neural differentiable framework for predicting density and grain-size evolution. Two neural networks learn densification and grain-growth coefficients within coupled rate equations, while a smooth saturation factor attenuates densification near theoretical density. The same governing structure, network architecture, and training procedure were fitted independently to published data for MgO, Al-doped ZnO, and CaO-doped ThO2. Tests at held-out temperatures and compositions yielded the lowest mean error in all twelve material-metric comparisons against multilayer perceptron and residual network baselines. For MgO, Al-doped ZnO, and CaO-doped ThO2, respectively, density normalized root-mean-square errors were 14.6%, 10.8%, and 14.4%, and grain-size errors using the same metric were 8.6%, 12.1%, and 19.3%. Removing evolving density from both neural-network inputs increased density and grain-size trajectory errors in all three systems and ten of twelve aggregate errors, supporting density-dependent kinetic feedback. Deep ensembles estimated model disagreement, but empirical coverage showed that the uncertainty bands were not calibrated and did not capture all model-data discrepancies. These results establish Sinter-PiNDiff as a retrainable framework for sparse-data prediction and uncertainty-informed selection of sintering conditions.
Read the full abstract

All authors Zeping Chen, Ani Aprahamian, Khachatur V. Manukyan, Tengfei Luo

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

HOMO-LUMO Intervals from Modern Density Functionals as Fundamental Gap Estimators

Use of the Kohn-Sham eigenvalue gap (highest occupied molecular orbital to lowest unoccupied, HOMO-LUMO) as a fundamental gap (I-A, I= ionization potential, A= electron affinity) estimator is a long-standing challenge for relatively simple exchange-correlation functionals, e.g. meta-GGAs. The LAK meta-GGA [T. Lebeda et al. Phys. Rev. Lett. 133, 136402 (2024)] was developed to emphasize certain constraints in ways intended to make the LAK KS gaps close to the HSE06 hybrid functional [A.V. Krudau et al. J.Chem. Phys. 125, 224106 (2006)] gaps. The LAK functional was tested against a small atomization data set and a very small bond length set, with no molecular HOMO-LUMO gap testing. Here we test LAK against widely used, large reference data sets for molecular heats of formation, bond lengths, harmonic frequencies, I, and A and compare with r2SCAN [J.W. Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)] as well as HSE06 and PBE0, a quite different hybrid. On these tests LAK does not give meaningfully better geometries, heats of formation, HOMO-LUMO gaps, or DFT IP theorem satisfaction than r2SCAN. Compared against experimental I and A values the same holds; LAK results resemble r2SCAN results, not HSE06. Moreover, the HOMO and LUMO shifts generated by the NCAPR GGA [J. Carmona-Espindola et al., J. Chem. Phys. 157, 114109 (2022)] as derivative discontinuity estimates, correct the bare NCAPR HOMO and LUMO eigenvalues to substantially better agreement with experiment for both fundamental gap and IP theorem than do the meta-GGAs (LAK, r2SCAN) or hybrids (HSE, PBE0). We discuss possible reasons why LAK behaves differently for molecules and solids.
Read the full abstract

All authors H. Francisco, Javier Carmona-Espíndola, A. C. Cancio, José Luis Gázquez, S. B. Trickey

Materials Science

Crystallographic-orientation dependence of the early-stage oxidation of Zr single crystals: an XPS study of suboxide formation and in-depth distribution

The influence of crystallographic orientation on the early stages of the oxidation of pure Zr single crystals was studied by X ray photoelectron spectroscopy (XPS). Two samples, cut from the same a-Zr single crystal, were oxidized simultaneously at room temperature and O2 pressures of 1 X 10-8, 1 X 10-7 and 1 X 10-6 torr: one with its surface normal parallel to the c-axis (Z1, basal (0001) plane) and the other with its normal 12deg from that of a prismatic plane (Z2, near-prismatic orientation). At all three pressures the oxidation kinetics of both samples followed a three stage, logarithmic type behaviour, but the near prismatic sample (Z2) incorporated oxygen faster and to a greater extent than the basal sample (Z1) in every case. Deconvolution of the Zr 3d core level spectra resolved, in addition to metallic Zr and ZrO2, two sub-stoichiometric Zr O compounds, denoted (ZrO)a and (ZrO)b, with binding-energy shifts of 1.3 and 2.3 eV with respect to metallic Zr, consistent with sub oxides previously reported for polycrystalline Zr. Angle resolved XPS showed that ZrO2 is the outermost compound in both orientations, while the two sub-oxides are distributed nearly homogeneously through the film. Oxide thicknesses, calculated from the attenuation of the metallic Zr 3d signal, ranged from 10 to 13 A for Z1 and from 14 to 19 A for Z2 depending on the oxidation pressure, in close agreement with earlier measurements on cold worked polycrystalline Zr, whose surface texture is dominated by prismatic-oriented grains.
Read the full abstract

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

Scaling Density Functional Theory with Gaussian Splatting

Density functional theory (DFT) strikes a practical balance between accuracy and computational cost in many problems of computational chemistry and materials science. However, many DFT calculations are limited by fixed atom-centered basis sets, which dictate how accuracy and cost scale with system size. We propose Gaussian Splatting for Density Functional Theory (GS-DFT), which represents molecular orbitals as a cloud of Gaussians whose positions, shapes, and mixing coefficients are optimized jointly by gradient descent to minimize the energy without training data. Conceptually, GS-DFT is 3D Gaussian splatting with the renderer replaced by quantum mechanics. We introduce two key solver components: adaptive density fitting with screening for efficient evaluation of two-electron integrals, and a regularized differentiable orthogonalization of the molecular orbitals. Empirically, the optimized basis reaches the accuracy of the largest conventional basis sets with a fraction of the parameters, converging systematically in energy, density, and nuclear forces. At equal parameter count, it captures the stretched-bond and anion physics that fixed bases only recover with specialized basis augmentation. The resulting solver exhibits quadratic peak memory scaling in the cloud size, allowing us to simulate systems of up to 2,742 atoms (10,406 electrons) without any modifications at triple-zeta scale using a single four-GPU node.
Read the full abstract

All authors Andrés Guzmán-Cordero, Cindy Zhang, Majdi Hassan, Marta Skreta, Kirill Neklyudov, Matija Medvidović

Materials Science

Ferroelectric Switching in ZnO/Zn1-xMgxO Heterostructures: Atomistic Insights into Interfacial Coupling and Layer Architecture

Ferroelectric switching in heterostructures couples composition, layer topology, temperature, and interfacial boundary conditions. ReaxFF molecular dynamics isolates these variables in ZnO/ Zn1- xMgxO/ZnO and Zn1-xMgxO/ZnO/ Zn1-xMgxO stacks. Within pristine, initially single-domain models, coupling to switchable Zn1-xMgxO reduces the applied field required to reverse ZnO by up to fivefold. Temperature generally lowers the coercive field, whereas Mg concentration produces a nonmonotonic response. At equal ZnO and Zn1-xMgxO (ZMO) proportions, structures with ZnO at the center switch at lower fields than those with ZMO at the center at all four temperatures examined, demonstrating a topology-dependent response. Layer-resolved trajectories reveal topology-dependent switching sequences with direction-dependent redistribution of normal stress near the heterointerfaces, consistent with a stress-assisted cooperative pathway. Limiting MgO-containing structures exhibit sequential multilevel switching or low-polarity trapping, depending on thickness and temperature. Fixed-charge atomistic simulations complement previous continuum descriptions by resolving structural, energetic, and local stress evolution under a common applied field. TEM and STEM-EDS observations provide experimental structural context for the modeled architectures. Together, the results establish layer topology and interfacial mechanical confinement as design variables for wurtzite ferroelectric heterostructures.
Read the full abstract

All authors Alireza Sepehrinezhad, Ali Mohammadi Dinani, Ece Gunay, Elizabeth C. Dickey, Jon-Paul Maria, Susan Trolier-McKinstry, Adri C.T. van Duin

Materials Science · Cross-listed from physics.optics

High-Yield, 150 ps Polymer Scintillators via Synergistic Purcell and High-Z Enhancement

Developing plastic scintillators that simultaneously achieve high light yield (LY), ultrafast timing, and scalable processability remains a formidable challenge in radiation detection. Here, we report a solution-processable hybrid scintillator that integrates high-Z sensitization with Purcell-enhanced radiative recombination within a conjugated-polymer platform. By embedding HfO2 nanoparticles and resonant Ag-SiO2 plasmonic nanoantennas into a poly(9,9-dioctylfluorene) matrix, we exploit the intrinsic sub-nanosecond kinetics of conjugated polymers while overcoming their characteristic low stopping power. The HfO2 nanocrystals amplify energy deposition and suppress interchain aggregation, preserving high photoluminescence quantum yields in the solid state. Simultaneously, the nanoantennas induce strong plasmon-exciton coupling, accelerating radiative decay to ~150 ps. This synergistic architecture yields LY=6.1 photons/keV, comparable to standard commercial plastics with nanosecond lifetimes, yet over an order of magnitude higher than current sub-nanosecond state-of-the-art materials. This performance represents one of the most favourable combinations of efficiency and timing resolution reported to date. Our versatile strategy provides a framework for next-generation organic scintillators, merging record-breaking temporal response with high efficiency for advanced radiation-sensing architectures.
Read the full abstract

All authors Xiaohe Zhou, Hiba H. Karakkal, Chenger Wang, Andrea Fratelli, Francesco Bruni, Matteo L. Zaffalon, Leonardo Poletti, Saptarshi Chakraborty, Emanuele Mazzola, Roberto Lorenzi, Francesco Carulli, Francesca Rossi, Francesco Meinardi, Luca Gironi, Sergio Brovelli

Materials Science

Dielectric response of the antiferromagnetic multiferroics with cycloidal equilibrium for the polarization related to the noncollinear spins

Linear dynamics of the small amplitude perturbations in the antiferromagnetic multiferroics is analytically considered. Equilibrium state in the considered regime is assumed to be cycloidal structure both for the antiferromagnetic vector and the nonzero magnetization vector (presence of weak ferromagnetism in the system is included). The electric polarization of spin origin related to the "noncollinear" spins is considered and its influence on the spin dynamics is included. The dynamical electric susceptibility tensor is obtained to demonstrate response of the multiferroics on the alternating electric field. Existence of two polarization projections is found, response of each of them on two projections of the electric field perpendicular to the cycloid direction is presented. Different roles of the in plane and out of plane electric perturbations is demonstrated analytically. Essential role of the low frequency gapless spin wave in the dielectric response is described. Contributions of the spiral structure and the weak ferromagnetism in the dielectric response are specified.
Read the full abstract

Materials Science

Polar discontinuity screening by charge disproportionation in ferroelectric-nickelate superlattices

We use density functional theory (DFT) to demonstrate rearrangement of charge disproportionation as a new mechanism for screening polar discontinuities at heterointerfaces. We focus on SmNiO$_3$, a III-III perovskite characterized by an insulating ground state that originates from disproportionation of the electrons at the Ni sites. When interfaced with ferroelectric, II-IV perovskite BaTiO$_3$, this system can display a polar discontinuity whose magnitude is determined by the nature of the interface (NiO$_2$-BaO or SmO-TiO$_2$) and the orientation of the BaTiO$_3$ spontaneous polarization. In the orientation maximizing the interface charge, SmNiO$_3$ compensates by transferring 1 $e$/f.u. between the two interfaces and modifying the disproportionation sequence from Ni$^{4+}$/Ni$^{2+}$ to Ni$^{2+}$/Ni$^{2+}$ at the SmO-TiO$_2$ interface and to Ni$^{4+}$/Ni$^{4+}$ at the NiO$_2$-BaO interface. While this mechanism is potentially universal in all perovskites exhibiting charge disproportionation, our formal analysis of BiNiO$_3$/BaTiO$_3$ superlattices suggests that larger polar discontinuities may need more complex disproportionation patterns for screening.
Read the full abstract

All authors Edith Simmen, Jonathan Amara, Philippe Ghosez, Nicola A. Spaldin

Materials Science

A quantitative phase-field model for grain boundary trapping and diffusion of hydrogen

Grain boundaries influence hydrogen transport in polycrystalline metals by acting as both trapping sites and interconnected diffusion pathways. Conventional interpretations of thermal desorption spectroscopy (TDS) and permeation experiments often treat traps as isolated defects, neglecting grain boundary connectivity and potentially misinterpreting experimental data. Here, we develop a quantitative phase-field model for hydrogen diffusion, grain boundary trapping, and grain boundary-assisted transport. The formulation preserves the physical grain boundary volume independently of the numerical interface thickness and accounts for anisotropic diffusion parallel and perpendicular to grain boundaries. Benchmark simulations verify quantitative behavior under interface upscaling. The model shows that grain boundary diffusion can significantly shift TDS peaks, making Kissinger-type analyses unreliable when grain boundary transport is significant. Effective diffusion coefficients are also shown to depend strongly on grain size, trapping free energy, temperature, and grain boundary mobility. These results highlight the coupled role of trapping thermodynamics, bulk diffusion, and grain-boundary transport.
Read the full abstract

Materials Science

High-throughput computational discovery of CuI-based ternary $p$-type transparent conductors

P-type transparent conducting materials (TCMs) remain scarce, limiting progress in transparent electronics and tandem photovoltaics. Here, we perform a high-throughput computational search across CuI-based ternary X–Cu–I compounds using the minima hopping method combined with density functional theory. Filtering for thermodynamic stability, optical transparency, and light hole effective mass, we identify 58 candidate p-type TCMs, of which 49 are previously unreported. Twenty are p-type degenerate semiconductors with intrinsic hole carriers, including Cu$_2$FI$_2$ (calculated gap 3.75 eV, hole effective mass $0.262\,m_0$). We further use this dataset to revisit the chemical modulation of the valence band (CMVB) design rules, originally formulated for Cu-based delafossite oxides. Four generalizations emerge: linear Cu coordination is not required, as tetrahedral coordination (CN = 4) is compatible with p-type TCM character; the mechanism extends to halide-based systems with I $5p$ anion states; the Cu $d$ state energy center $\epsilon_d$ is a quantitative descriptor for valence band dispersion; and Cu undercoordination in monovalent-anion hosts drives intrinsic p-type carrier generation through stoichiometric charge balance. These results establish a generalized computational framework for the design of p-type halide TCMs.
Read the full abstract

Materials Science

Multiscale computational study of the dielectric response of semi-crystalline polyethylene with chemical defects

Polymers are widely used as functional insulating materials, but predicting their dielectric behavior is challenging because multiple mechanisms, acting across different spatial and temporal scales, contribute to their response. In polyethylene (PE), one of the most common polymers, electronic, atomic, and mesoscale dynamics must all be considered. In this work, semicrystalline models of polyethylene were developed and investigated using a multiscale approach. Quantum simulations were employed to describe electronic and vibrational contributions to the dielectric response, while classical molecular dynamics simulations captured the behavior of polymer chains at room temperature and frequencies down to the 10-100 MHz range. The study focuses on the impact of radio-oxidation defects on the dielectric properties of polyethylene. Quantum calculations reveal how electronic and vibrational contributions depend on the local atomic environment surrounding these defects. Comparison with molecular dynamics results highlights the influence of temperature on the static dielectric constant. Structural analyses further assess the effect of each defect type on PE crystallinity. Analysis of dipolar correlation functions shows that different defects interact with one another and with the polymer matrix in distinct ways, affecting permittivity and dielectric loss peaks. For all defect types, defect-defect interactions provide the largest contribution to dielectric response. While most oxidized groups exhibit positive coupling with the PE matrix, alcohol defects display a negative cross-correlation, partially offsetting their impact on the dielectric properties. Our results also show that ketone groups produce the largest dielectric loss between the defects studied.
Read the full abstract

Materials Science

First-principles calculations of the thermodynamic factor of multicomponent alloys

Computing non-dilute diffusion coefficients from Onsager transport coefficients of a multicomponent alloy requires the thermodynamic factor, which measures the curvature of the free energy with respect to composition. As diffusion in metals is mediated by vacancies, the curvature that matters is not that of the vacancy-free alloy, which is straightforward to compute, but that of an alloy carrying a dilute concentration of vacancies. The resulting matrix is nearly singular, and composition fluctuations from which it is obtained converge slowly in Monte Carlo simulations. Here we develop two routes that avoid sampling vacancy composition fluctuations. The first constructs the thermodynamic factor from two quantities, the free-energy curvature of the vacancy-free alloy and the vacancy concentration. The second truncates the semi-grand canonical partition function at a single vacancy and recovers the full thermodynamic factor of the vacancy-containing alloy. The same separation shows that the largest eigenvalue of the diffusion matrix is the vacancy tracer diffusion coefficient for any number of components and any degree of non-ideality, while the remaining eigenvalues scale with the vacancy concentration. Across binary through quinary alloys of the Hf–Mo–Nb–Ti–Zr system, both routes reproduce eigenvalues and eigenvectors of diffusion matrices obtained from fully converged Monte Carlo simulations. The single-vacancy expansion reaches this accuracy with roughly one tenth of the sampling effort. A scaling relation extends both routes to local vacancy chemical potentials away from equilibrium, without additional simulations. These results allow thermodynamic descriptions from atomistic models or CALPHAD assessments to be used directly in mesoscale simulations of mass transport.
Read the full abstract

All authors Damien K. J. Lee, Shashank Saxena, Anton Van der Ven, Anirudh Raju Natarajan

Materials Science

Electron paramagnetic resonance of Dy${^3+}$-doped CaWO$_4$: spin Hamiltonian, crystal-field analysis, and linear electric field effect

We report an electron paramagnetic resonance (EPR) study of Dy${^3+}$ in CaWO$_4$ single crystals spanning dopant concentrations from 80~ppb to 247~ppm. The $g$ factors of the ground Kramers doublet at the $S_4$ point-group site, $g_\parallel = 7.22 \pm 0.01$ and $g_\perp = 5.45 \pm 0.01$, and the hyperfine constants of \ce{^{161}Dy} and $^{163}$Dy are determined with an order-of-magnitude improvement in precision over the only previous report. A crystal-field analysis constrained jointly by published optical levels and by the measured $g$ tensor identifies the ground doublet as the lower branch of an anticrossing between the near-degenerate $\|{\pm}11/2>$ and $\|{\mp}13/2>$ states, and explains why optical data alone left the $g$ values unconstrained. \add{The same parameter set, tested on the excited multiplets of the $^6H$ term, reproduces the ${}^6H_{13/2}$ crystal-field levels better than the optical parameters themselves.} The angular dependence of the linewidth in the $ab$ plane reveals broadening by random internal electric fields through the linear electric field effect. Calibrating these fields with the residual Er$^{3+}$ present in the same crystals yields the first electric-field coupling parameters of Dy$^3+$, $(B_{31}^2+B_{36}^2)^{1/2} = (34\pm5)\times10^{-6}$~cm/V, three times the Er$^{3+}$ value and the largest reported for a rare-earth ion in this host.The residual linewidth and the concentration-dependent lineshape asymmetry follow the quadratic field scaling expected for a second-order Stark shift, which escapes the inversion-site cancellation constraining the linear effect. Both couplings trace back to the $\sim 20$~cm$^{-1}$ gap to the first opposite-symmetry doublet, a consequence of the dense crystal-field structure of the $J=15/2$ manifold.
Read the full abstract

All authors Achuthan Manoj Kumar, Larissa Aboudem-Joumessi, Remy Dassonneville, Adrien Savoyant, Karolina Waszowska, Pengrui Jiao, Patrice Bertet, Philippe Goldner, Sylvain Bertaina

Materials Science

Current-voltage characteristics and resistive switching in epitaxial La$_{0.67}$Sr$_{0.33}$MnO$_3$/SrMnO$_3$/La$_{0.67}$Sr$_{0.33}$MnO$_3$ multilayer

In present study, we investigate the structural properties and current-voltage ($I$-$V$) characteristics of an epitaxial multilayer composed of La$_{0.67}$Sr$_{0.33}$MnO$_3$ (25 nm)/SrMnO$_3$ (70 nm)/La$_{0.67}$Sr$_{0.33}$MnO$_3$ (25 nm), grown on different substrates i.e., SrTiO$_3$(100), LaAlO$_3$(100), Si (100). The used substrates not only have different chemical compositions in line with the film materials, also they have different lattice parameters which would tune the material chemistry and lattice strain at the interface significantly. Our $I$-$V$ measurements reveal distinct electrical behaviors depending on the substrate. The multilayers grown on oxide SrTiO$_3$(100) and LaAlO$_3$(100) substrates exhibit regular $I$-$V$ with slight nonlinearity at low applied voltages. In contrast, $I$-$V$ in multilayer with Si(100) substrate exhibits large asymmetry and notable resistive switching (RS) behavior at room temperature, where a resistance ratio around 10 between the high resistance state (HRS) and low resistance state (LRS) is observed. While the observed $I$-$V$ behavior is sensitive to substrate-induced interfacial disorder and formation of trap bands within bandgap of SrMnO$_3$, the present results have potential applications in future memory and memristive devices.
Read the full abstract

All authors A. G. A. Rahman, R. K. Patel, Chandrani Nath, A. K. Pramanik

Materials Science

Modern triad for magnetic material science: geometric shape, electronic bands and spin textures

To enable novel concepts in nanotechnologies, it is insightful to design materials with qualitatively new performances originating from recent discoveries in fundamental material science and physics. Curvilinear magnetism emerged as a tool to design chiral and anisotropic responses of materials at the nanoscale using the effects of the geometric shape and topology of the object. This concept is distinct yet complementary to material screening, which is primarily based on the optimization of intrinsic microscopic properties to modify magnetic textures for specific applications. Contemporary curvilinear magnetism focuses on mean-field micromagnetics, while the explicit contribution of the lattice structure and electronic degrees of freedom remains beyond consideration. In this review, we will emphasize on the effects which appear when this coupling is accounted for and that have resulted in novel concepts for the design of “metageometric” materials. Further development of this field will provide appealing technological prospects for fundamental understanding of curvilinear magnetic nanostructures and their applications in energy efficient and scalable nanoelectronics.
Read the full abstract

All authors Denys Makarov, Oleksandr Pylypovskyi, Rui Xu, Carmine Ortix