2025-11-17 | | Total: 68
Interoperability and cross-validation remains a significant challenge in the computational materials discovery community. In this context, we introduce a common input/output standard designed for internal translation by various workflow managers (AiiDA, PerQueue, Pipeline Pilot, and SimStack) to produce results in a unified schema. This standard aims to enable engine-agnostic workflow execution across multiple density functional theory (DFT) codes, including CASTEP, GPAW, Quantum ESPRESSO, and VASP. As a demonstration, we have implemented a workflow to calculate the open-circuit voltage across several battery cathode materials using the proposed universal input/output schema. We analyze and resolve the challenges of reconciling energetics computed by different DFT engines and document the code-specific idiosyncrasies that make straightforward comparisons difficult. Motivated by these challenges, we outline general design principles for robust automated DFT workflows. This work represents a practical step towards more reproducible and interoperable workflows for high-throughput materials screening, while highlighting challenges of aligning electronic properties, especially for non-pristine structures.
Spins in semiconductor quantum dots offer a scalable approach to quantum computing; however, precise control and efficient readout of large quantum dot arrays remain challenging, mainly due to the hyperdimensional voltage space required for tuning multiple gates per dot. To automate this process, large datasets are required for testing and training autotuning algorithms. To address the demand for such large datasets, we introduce RF-Squad, a physics-based simulator designed to realistically replicate radiofrequency (RF) reflectometry measurements of quantum dot arrays, with the ability to go beyond the Constant Interaction Model (CIM) and simulate physical phenomena such as tunnel coupling, tunnel rates, and quantum confinement. Implemented in JAX, an accelerated linear algebra library, RF-Squad achieves high computational speed, enabling the simulation of a 100x100 pixel charge stability diagram of a double quantum dot (DQD) in 52.1 $\pm$0.2 milliseconds at the CIM level. Using optimization algorithms, combined with it's layered architecture, RF-Squad allows users to balance physical accuracy with computational speed, scaling from simple to highly detailed models.
Nonlinear processes can offer pathways to next-generation sensors and frequency mixing devices to overcome modern imaging, detection, and communication challenges. In this article, we report on strong second and third-order nonlinear optical responses in visible and terahertz (THz) light in single crystals of the noble metal dichalcogenide PdTe2. We find that buried conduction and valence topological surface states of PdTe2 lead to resonant enhancement of optical second-harmonic generation. On the other hand, although the nonlinear responses obtained with THz excitation are not close to this resonance, they can be clearly observed in reflection geometry, even in the presence of broadband excitation, where optical filters are not necessary to observe the odd-order higher harmonic gain. By carefully considering the radiative photocurrent framework of stimulated THz emission, we are able to extract fingerprints of both second- and third-order processes in the THz regime, and show that PdTe2 is a promising material candidate for radio frequency rectification, frequency mixing, and beam focusing.
The manipulation of magnetization in ferromagnetic metals (FMs) through orbital torque (OT) has emerged as a promising route for energy-efficient magnetic devices without relying on heavy metals. While Ti and Cu are among the most extensively studied light nonmagnetic metals (NMs) for OT devices, theoretical calculations of the resulting torque have remained limited. Here, we present a systematic and quantitative theoretical study of current-induced torques in Ti/FM and Cu/FM (FM = Co, Ni) bilayers using semi-realistic tight-binding models derived from ab initio electronic structures. We find that the torque in Ti/FM is larger for Ni than for Co, but this trend does not necessarily hold in Cu/FM, revealing that the FM dependence of OT is not universal but varies with the orbital current source. Moreover, the dependence of OT on NM thickness clearly indicates its NM bulk origin in both Ti- and Cu-based systems. Notwithstanding, the quantitative characteristics of OT cannot be explained by a simplified picture based on the individual bulk properties of the NM or FM layers. These results provide microscopic insight and practical guidance for designing light-metal-based orbitronic devices.
Angle-resolved photoemission spectroscopy (ARPES) has recently emerged as a direct probe of excitonic correlations in two-dimensional semiconductors, resolving their dispersion and dynamics in energy-momentum space, including dark exciton states inaccessible to optical techniques. However, the ARPES fingerprint of charged excitons (trions), which plays a key role in all doped and gated 2D material systems, has remained unknown so far. We present a first theoretical analysis of trion signatures in monolayer transition-metal dichalcogenides, highlighting how the additional charge carrier modifies the spectral position and shape relative to neutral excitons in ARPES spectra. Interestingly, we further predict that mass-imbalanced trions yield a characteristic double-peak structure, clearly separated in energy and line shape from neutral excitons. The predicted temperature dependence of these features offers guidance for experimental investigations aimed at identifying trionic states, thereby establishing a framework for ARPES studies of many-body Coulomb complexes in doped two-dimensional semiconductors.
The search for materials with intrinsically low thermal conductivity ($κ_L$) is critical for energy applications, yet conventional descriptors often fail to capture the complex interplay between bonding and lattice dynamics. Here, first-principles calculations are used to contrast the thermal transport in covalent zincblende (zb) and metavalent rocksalt (rs) BeO. We find that the metavalent bonding in rs-BeO enhances lattice anharmonicity, activating multi-phonon scattering channels and suppressing phonon transport. This results in an ultralow $κ_L$ of 24 W m$^{-1}$ K$^{-1}$ at 300 K, starkly contrasting with the zb phase (357 W m$^{-1}$ K$^{-1}$). Accurately modeling such strongly anharmonic systems requires explicit inclusion of temperature-dependent phonon renormalization and four-phonon scattering. These contributions, negligible in zb-BeO, are essential for high-precision calculations of the severely suppressed $κ_L$ in rs-BeO. Finally, we identify three key indicators to guide the discovery of metavalently bonded, incipient-metallic materials: (i) an NaCl-type crystal structure, (ii) large Grüneisen parameters ($\textgreater$2), and (iii) a breakdown of the Lyddane-Sachs-Teller relation. These findings provide microscopic insight into thermal transport suppression by metavalent bonding and offer a predictive framework for identifying promising thermoelectrics and phase-change materials.
We propose a microscopic, weak-coupling mechanism by which generic Chern bands relax toward ideal bands. We consider coupling interacting electrons to a Caldeira-Leggett like Ohmic bosonic bath. Using the Born-Markov approximation, Slater determinant states of a Chern band under Hartree-Fock approximation evolve toward Slater determinant states corresponding to an ideal Chern band. We validate our proposal by performing numerical simulation of a massive Dirac model, showing that the Berry curvature and quantum metric indeed co-evolve to saturate the trace condition. Our proposal provides a concrete dissipative route to realize ideal Chern bands, a fundamental building block for the stabilization of fractional Chern insulators.
Nanoscale topologically non-trivial magnetization configurations generate significant interest due to both the fundamental properties of their knotted structures and their potential applications in ultra-efficient computing devices. While such textures have been widely studied in two dimensions, three-dimensional (3D) systems can yield more complex configurations, resulting in richer topologies and dynamic behaviors. However, reliably nucleating these 3D textures has proven challenging, and so far, 3D configurations such as vortex rings and hopfions can often only be observed forming spontaneously in relatively uncontrolled manners. Here, we demonstrate that through the 3D nanopatterning of chiral single crystal helimagnets into nano-tori, the controlled formation of a magnetic double helix can be achieved. This surface-localized topological state is stabilized by the interplay of intrinsic exchange interactions of the single crystal with the extrinsic emergent effects of the patterned geometry. These double helices host magnetic defects akin to supercoiling in circular DNA and climbing vines. We expect this study to serve as a foundation for future research combining single crystal systems with 3D nanopatterning, offering a new degree of control over emergent phenomena in nanoscale magnets and wider quantum material systems.
The non-Hermitian skin effect (NHSE), characterized by the extensive localization of bulk modes at the boundaries, has attracted significant attention as a hallmark feature of non-Hermitian topology. This localization invalidates the conventional Bloch band theory, necessitating an analysis under open boundary conditions even in the thermodynamic limit. The Amoeba formulation addresses this challenge by computing the spectral potential rather than the spectrum itself. Based on the (strong) Szegö limit theorem and its topological generalization, this approach reduces the evaluation of the potential to an optimization problem involving the Ronkin function. However, while the generalized Szegö limit theorem is formally applicable in arbitrary dimensions, its implementation is limited to single-band systems, and its applicability to multiband systems remains unclear. In this paper, we establish the Wiener-Hopf factorization (WHF) of the non-Bloch Hamiltonian as a powerful framework, providing a unified and rigorous foundation for Amoeba analysis in multiband systems. By combining the WHF with Hermitian doubling, we first elucidate the applicability criteria for the generalized Szegö limit theorem in multiband systems. We then show that the WHF provides the natural mathematical origin for the symmetry-decomposed Ronkin function in symmetry class AII$^\dagger$, leading to a rigorous proof of the generalized Szegö limit theorem for these systems and opening a path toward systematic generalizations to other symmetry classes.
We develop a Fokker-Planck approach to describe the dynamics of staggered magnetization and thermal fluctuations in a two-dimensional antiferromagnetic system with uniaxial anisotropy. Beginning with a classical model for the antiferromagnetic system, we incorporate a Landau-Lifshitz-Gilbert equation augmented by Langevin fields to account for thermal fluctuations, and we derive the Fokker-Planck equation governing the probability distribution function of the spin configuration. Employing the mean-field approximation, we derive the equations of motion for the spin polarization and the two-time spin-spin correlation functions. The methodology is applied to the study of spin-wave dynamics and to the formulation of a phenomenological model for resistance fluctuations in two-dimensional antiferromagnetic semiconductors.
The Eigenstate Thermalization Hypothesis (ETH) has been established as a cornerstone for understanding thermalization in quantum many-body systems. Recently, there has been growing interest in the full ETH, which extends the framework of the conventional ETH and postulates a smooth function to describe the multi-point correlations among matrix elements. Within this framework, free cumulants play a central role, and most previous studies have primarily focused on closed systems. In this paper, we extend the analysis to a system-bath setup, considering both an idealized case with a random-matrix bath and a more realistic scenario where the bath is modeled as a defect Ising chain. In both cases, we uncover a universal scaling of microcanonical free cumulants of system observables with respect to the interaction strength. Furthermore we establish a connection between this scaling behavior and the thermalization dynamics of the thermal free cumulants of corresponding observables.
Room temperature superconductivity remains elusive, and hydrogen-base compounds despite remarkable transition temperatures(Tc) typically require extreme pressures that hinder application. To accelerate discovery under moderate pressures, an interpretable framework based on symbolic regression is developed to predict Tc in hydrogen-based superconductors. A key descriptor is an integrated density of states (IDOS) within 1 eV of the Fermi level (EF), which exhibits greater robustness than conventional single-point DOS features. The resulting analytic model links electronic-structure characteristics to superconducting performance, achieves high accuracy (RMSEtrain = 20.15 K), and generalizes well to external datasets. By relying solely on electronic structure calculations, the approach greatly accelerates materials screening. Guided by this model, four hydrogen-based candidates are identified and validated via calculation: Na2GaCuH6 with Tc =42.04 K at ambient pressure (exceeding MgB2), and NaCaH12, NaSrH12, and KSrH12 with Tc up to 162.35 K, 86.32 K, and 55.13 K at 100 GPa, 25 GPa, and 25 GPa, respectively. Beyond rapid screening, the interpretable form clarifies how hydrogen-projected electronic weight near EF and related features govern Tc in hydrides, offering a mechanism-aware route to stabilize high-Tc phases at reduced pressures.
Multiple Andreev reflections in voltage-biased Josephson junctions give rise to the subharmonic gap structure in the conductance, which is widely used to characterize transport properties and estimate the induced gap in the junctions. Here we theoretically investigate the evolution of the subharmonic gap structure in spinful Josephson junctions. Spin mixing introduced by the spin-orbit coupling opens avoided crossings in the dispersion relation of the leads, which, as we demonstrate, subsequently results in pronounced multiple Andreev reflection features in the conductance traces. We analyze how these features evolve under an external magnetic field and explain that their visibility in conductance is governed by the spin polarization of the bands.
We theoretically investigate orbital accumulation driven by chiral phonons via orbital-dependent electron-lattice coupling. We derive a formula for the orbital accumulation induced by classical lattice dynamics or nonequilibrium phonons, using the Berry curvature, linear response theory, and the nonequilibrium Green's function method. We show that chiral phonons primarily couple to orbital quadrupole moments and that static orbital dipole accumulation can be generated in the second-order of lattice displacement. Our study provides a useful method for generating orbital accumulation without using spin-orbit interaction.
A fundamental characteristic of a superconducting state is the coherence length $ξ$. Multicomponent superconductors, particularly ones breaking multiple symmetries, are characterized by multiple coherence lengths. Here we show that even, nominally $single$-component superconductors under certain conditions are characterized by multiple coherence lengths. We consider nearest-neighbor pairing interactions on a square lattice that leads to $s$-wave and $d$-wave representations of link superconducting order parameter. We show that even if the subdominant order parameter is completely suppressed in the ground state, it results in multiple correlation lengths with nontrivial hierarchy, resulting in important physical consequences in inhomogeneous solutions. Under certain conditions, this leads to type-1.5 superconductivity, where magnetic field penetration length falls between two coherence lengths, leading to vortex clustering in an external magnetic field.
We derive the full spectrum of decorated Cayley trees that constitute tree analogs of selected two-dimensional Euclidean lattices; namely of the Lieb, the double Lieb, the kagome, and the star lattice. The common feature of these Euclidean lattices is that their nearest-neighbor models give rise to flat energy bands interpretable through compact localized states. We find that the tree analogs exhibit similar flat or nearly flat energy bands at the corresponding energies. Interestingly, such flat bands in the decorated Cayley trees acquire an interpretation that is absent in their Euclidean counterparts: as edge states localized to the inner or the outer boundary of the tree branches. In particular, we establish an exact correspondence between the Lieb-Cayley tree and an ensemble of one-dimensional Su-Schrieffer-Heeger chains, which maps topological edge states on one side of the chains to flat-band states localized in the bulk of the tree, furnishing the flat energy band with a topological stability. Similar mapping to topological edge states or to states bound to edge defects in one-dimensional chains is shown for flat-band states in all the considered tree decorations. We finally show that the persistence of exact flat bands on infinite decorated trees (i.e., Bethe lattices) arises naturally from a covering interpretation of tree graphs. Our findings reveal a rich landscape of flat-band and topological phenomena in non-Euclidean systems, where geometry alone can generate and stabilize unconventional quantum states.
This work presents electronic and magnetic properties of selected members in the cubic Laves phase series Nd1-xPrxCoNi and Ce0.25Pr0.75CoNi, together with the corresponding binary compositions (NdCo2, NdNi2, PrCo2, PrNi2, CeCo2, CeNi2), using soft x-ray absorption spectroscopy, x-ray magnetic circular dichroism (XMCD), density-functional theory, and crystal field multiplet calculations. All transition-metal moments saturate below 1 T, while the rare-earth moments do not saturate even at 5 T, consistent with van Vleck paramagnetic contributions and crystal field suppression. While the sum rules are widely used to extract element-specific magnetic moments from XMCD, we show that for 3d transition metals, their application requires accurate estimates of the number of unoccupied 3d states. We observe a finite magnetic moment on Ni, challenging the common assumption of its nonmagnetic character in Laves phases. The orbital magnetic moments were determined using the spin rules, while the spin moments were estimated from single-ion values from multiplet calculations, due to the invalidity of the spin sum rule for light rare-earth elements. The magnetic moments of Nd and Pr are found to be suppressed relative to their free-ion values, with multiplet theory indicating that this is due to crystal field effects. Our results confirm that Nd and Pr maintain localized 4f3 and 4f2 configurations, respectively, and that their element-specific magnetic moments are robust to rare-earth substitution. Ce, on the other hand, exhibits a tunable mixed-valent ground state with both magnetic 4f1 and nonmagnetic 4f0 components. The relative fraction of these states varies with the electronegativity of the surrounding 3d transition metals, revealing a pathway to tune Ce magnetism via composition. This work establishes a framework for accurately interpreting XMCD in light rare-earth-based intermetallics.
We study the collective behavior of a polar flock in an inhomogeneous environment in two-dimensions. The inhomogeneity is modelled by introducing regions at random locations on the substrate with higher noise but accessible for the flock to move. Hence inside such regions the particles orientation get randomised. Such inhomogeneities are different from the physical disorder, which obstructs the space for the incoming particles. The study focuses on how the phase behavior of polar flock changes by tuning the packing fraction of inhomogeneity. As packing fraction increases, the system crosses over from long-range to quasi long range order and ultimately to a disordered phase, while the order disorder transition for flocking changes from discontinuous to continuous. The resultant phase behavior of polar flock patterns here is comparable to that exhibited in the presence of physical disorder.
In this paper we attempt to construct the topological theory of superfluid helium $4$ in the framework of the (rigid) $U(1)$ model in which the initial $U(1)$ group is destroyed with appearance of (topologically nontrivial) domains separated by domain walls treated as step voltages between domains (e.g. with neighboring topological numbers). This can explain the superfluid properties in a helium $4$ specimen as well as the appearance of topologically nontrivial vortices therein.
Electron accumulation in interatomic regions is a fundamental quantum phenomenon dictating chemical bonding and material properties, yet its origin remains elusive across disciplines. Here, we report a counterintuitive quantum effect -- potential-barrier affinity (PBA) -- revealed by solving the Schrödinger equation for a crystalline potential. PBA effect drives significant interatomic electron accumulation when electron energy exceeds the barrier maximum. This effect essentially determines interatomic electron density patterns, governing microstructures and properties of condensed matters. Our theory overturns the traditional wisdom that the interstitial electron localization in electride requires potential-well constraints or hybrid orbitals, and it serves the fundamental mechanism underlying the formation of conventional chemical bonds. This work delivers a paradigm shift in understanding electron distribution and establishes a theoretical foundation for the microscopic design of material properties.
Systems of active particles can show a large variety of collective behavior. In theory, two aspects determine the collective behavior: the model at the particle level and the parameter regime. While many studies consider a single model and study its parameter regime, here, we focus on the former aspect. Motivated by experiments that study dilute suspensions of Chlamydomonas reinhardtii in a self-generated oxygen gradient, we compare various models with external field-dependent motility to understand how the collective behavior changes between models. We vary the particle-particle interaction from no interactions to steric interactions, the particle shape from round disks to dumbbells, the self-propulsion mechanism from constant speed to rocking motion, and the particle's center of mass from the geometric center to off-center. We find that changes in the model of the active agents can lead to similar statistics in the dilute regime and different collective behavior in the dense regime. We conclude that active particle models do not easily generalize for different real active agents, but instead require a clear understanding of the agents' microscopic properties.
Continuous-Time Recurrent Neural Networks (CTRNNs) have been widely used for their capacity to model complex temporal behaviour. However, their internal dynamics often remain difficult to interpret. In this paper, we propose a new class of CTRNNs based on Hopfield-like associative memories with asymmetric couplings. This model combines the expressive power of associative memories with a tractable mathematical formalism to characterize fluctuations in nonequilibrium dynamics. We show that this mathematical description allows us to directly compute the evolution of its macroscopic observables (the encoded features), as well as the instantaneous entropy and entropy dissipation of the system, thereby offering a bridge between dynamical systems descriptions of low-dimensional observables and the statistical mechanics of large nonequilibrium networks. Our results suggest that these nonequilibrium associative CTRNNs can serve as more interpretable models for complex sequence-encoding networks.
The electronic structure of semiconducting 2D materials such as transition metal dichalcogenides (TMDs) is known to be tunable by its environment, from simple external fields applied with electrical contacts up to complex van der Waals heterostructure assemblies. However, conventional alloying from reference binary TMD compounds to composition-controlled ternary alloys also offers unexplored opportunities. In this work, we use nano-angle resolved photoemission spectroscopy (nano-ARPES) and density functional theory (DFT) calculations to study the structural and electronic properties of different alloy compositions of bulk WS2(1-x)Se2x. Our results demonstrate the continuous variation of the band structure and the progressive evolution of the valence band splitting at the K points from 420 to 520 meV in bulk WS2(1-x)Se2x. We also carried out scanning tunneling microscopy (STM) measurements and DFT to understand the possible S or Se substitutions variants in WS2(1-x)Se2x alloys, with different local atomic configurations. Our work opens up perspectives for the fine control of the band dispersion in van der Waals materials and demonstrate how the band structure can be tuned in bulk TMDs. The collected information can serve as a reference for future applications.
We report a comprehensive investigation of the quasi-one-dimensional spin-chain compound Ca3CoIrO6 (CCIO) using a combination of structural, magnetic, thermodynamic, transport, Raman, and dielectric measurements. Temperature-dependent neutron powder diffraction confirms the rhombohedral R-3c structure down to 5 K without any structural phase transition. DC magnetization, ac susceptibility, and relaxation measurements reveal a gradual evolution from a high-temperature paramagnetic-like state to a partially disordered antiferromagnetic (PDA) state below 100 K, accompanied by slow cluster-like spin dynamics followed by a freezing transition near 30 K. Isothermal magnetic hysteresis M(H) loops demonstrate partial chain freezing, while robust exchange bias is observed in field-cooled protocols, highlighting the interplay between PDA ordering and frozen spins. Resistivity and specific heat data indicate strong coupling between spin and charge degrees of freedom, accompanied by activated transport behavior. Raman spectroscopy identifies pronounced anomalies in phonon frequencies and linewidths across multiple magnetic regimes, reflecting strong spin-lattice coupling. Polarization-electric field (P-E) measurements reveal temperature-dependent crossovers from linear dielectric to weakly hysteretic behavior, consistent with short-range polar correlations driven by spin-lattice interactions. These findings establish CCIO as a prototypical quasi-one-dimensional frustrated spin-chain system where geometric frustration, spin-orbit coupling, and low-dimensionality generate field-tunable PDA order, glassy spin dynamics, exchange bias, and magnetodielectric coupling. These results provide new insights into frustration-driven phases in low-dimensional oxides and point towards potential multifunctional applications based on intrinsic magnetodielectric and exchange bias phenomena.
The optical properties of the superconducting K$_{0.8}$Fe$_{1.7}$(Se$_{0.73}$S$_{0.27}$)$_2$ single crystals with a critical temperature $T_c\approx 26$ K have been measured in the {\it ab} plane in a wide frequency range using both infrared Fourier-transform spectroscopy and spectroscopic ellipsometry at temperatures of 4--300 K. The normal-state reflectance of K$_{0.8}$Fe$_{1.7}$(Se$_{0.73}$S$_{0.27}$)$_2$ is analyzed using a Drude-Lorentz model with one Drude component. The temperature dependences of the plasma frequency, optical conductivity, scattering rate, and dc resistivity of the Drude contribution in the normal state are presented. In the superconducting state, we observe a signature of the superconducting gap opening at $2Δ$(5~K) = 11.8~meV. An abrupt decrease in the low-frequency dielectric permittivity $\varepsilon _1(ω)$ at $T < T_c$ also evidences the formation of the superconducting condensate. The superconducting plasma frequency $ω_{pl,s} = (213\pm 5)$~cm$^{-1}$ and the magnetic penetration depth $λ=(7.5\pm 0.2)$~$μ$m at $T=5$~K are determined.