2026-09-07 | | Total: 99
Using low-temperature scanning tunneling microscopy we investigate the charge state, magnetic excitations, and Kondo features of individual Sm adatoms on graphene/Ir(111). Depending on the number and distance of their neighbors, Sm atoms can be in two discrete charge states. At certain distances, a reversible transition between these two states is induced by the electric field of the STM tip leading to concentric charge rings in the images. Only atoms in one of the two charge states exhibit magnetic excitations in d$I$/d$V$ spectra. Two such excitations are located at 35~meV and 54~meV and related to transitions from the $J = 1/2$ ground state doublet to the first crystal field split $J = 3/2$ multiplet. Together with the intra-atomic exchange excitations at higher energy, these observations indicate that Sm transfers one $6s$ electron to the substrate while it retains its gas-phase $4f$ filling. New for lanthanide adatoms, we observe a Kondo resonance. The Zeeman splitting of the Kondo peak reveals that Sm retains its large gas-phase $g$-factor. Comparison of d$I$/d$V$ spectra to cotunneling theory yields the crystal field acting on the $4f$ shell and, consequently, on the $J = 3/2$ quadruplet, confirms Sm$^+$ as the ground state, and identifies Sm$^{2+}$ as being energetically close, thereby rationalizing our observation of variable charge states.
Glyphosate is a widely used herbicide whose persistence and toxicity in aquatic and terrestrial environments demand efficient removal strategies. Here we employ GFN2 xTB calculations with implicit ALPB aqueous solvation and automated docking to investigate the adsorption of all five pH dependent ionized forms of glyphosate (G1-G5) on (10,0) single walled carbon nanotubes covalently functionalized with carboxyl groups at 0-25% coverage. Adsorption energies reveal a clear charge dependent trend: the dianionic (G4) and trianionic (G5) species exhibit the most negative binding energies over the entire functionalization range, while the protonated and neutral forms (G1, G2) bind weakly, approaching reversible adsorption at high COOH contents. The deprotonated form G5 strengthens from -2.17 eV on pristine CNT to about -5.7 eV at 10% and 25% functionalization, with a local minimum near 15-20% COOH due to steric and electrostatic crowding of adjacent groups. Decomposition of the solvation free energy shows dominant electrostatic stabilization complemented by increasingly favorable hydrogen bond contributions as carboxyl density grows. All complexes display very small HOMO-LUMO gaps (0.02-0.22 eV), indicating high electronic sensitivity to adsorption and functionalization. Overall, CNT+COOH systems can operate in both strong capture and regenerable regimes depending on glyphosate ionization state, offering a tunable platform for pH responsive remediation.
Non-reciprocal interactions are ubiquitous in active, biological, and disordered systems, generically driving them out of equilibrium. Here, we introduce a hierarchy of non-reciprocally interacting Ornstein-Uhlenbeck (NROU) models governed by a tunable non-reciprocity parameter $g$. At a special point $g=g^*$, the drift matrix becomes non-diagonalizable, realizing exceptional points (EP's) of different orders, where eigenvalues and eigenvectors simultaneously coalesce. The hierarchy encompasses non-reciprocally coupled dimers, their disordered counterparts, and a many-body chain exactly mapping onto the paradigmatic Hatano-Nelson model in the arena of non-Hermitian quantum systems. For the disordered model, we show that the distribution of the EP location $g^*$ across disorder realizations develops a universal edge singularity precisely at the clean-system EP, and is manifestly non-self-averaging. Across all models, we find that at the EP, the usual exponential relaxation of the autocorrelation and covariance functions is dressed by a polynomial-in-time prefactor whose degree is set by the order of the EP and whose detailed structure encodes the spatial architecture of the chain. At complete asymmetry, the many-body chain exhibits ``pseudo-equilibrium'': its steady-state distribution factorizes into equilibrium-like single-particle measures despite a nonzero steady-state current. Moreover, the $N$-particle interacting system decomposes into $N/2$ independent complex OU processes. Finally, using the Harada-Sasa relation, we obtain a closed-form expression for the total steady-state heat dissipation and uncover a boundary refrigeration effect, in which the boundary particles switch from acting as a hot to a cold reservoir as the non-reciprocity is tuned.
The advancement of spintronic nanoscale devices hinges on the ability to flexibly engineer magnetic spin structures in thin-film stacks with precision and control. Meeting this demand remains a challenge for stable non-collinear spin configurations and, more specifically, vertical spin spirals in thin films. Innovative methods are required for their fabrication, stabilization and control. Here, we use oblique-incidence deposition to design and stabilize vertical spin spirals at room temperature and without an external field in magnetic thin films. We induce two crossed uniaxial magnetic anisotropies at the thin film boundaries. These anisotropies are tunable in direction and strength, thus providing control over the angular range and depth profile of the resulting spin spiral. The combination of polarized neutron reflectometry and nuclear resonant scattering enables precise and direct determination of the depth-dependent spin configurations. Our results establish a single-film design approach, in which the surface anisotropies independently serve as controllable design parameters for tailoring the vertical spin-spiral profile. Potential applications include nanoscale energy-storage devices, magnetic sensors, and ferromagnetic-resonance filters, advancing all-spin-based device engineering in general.
We develop a method to extract the shear viscosity of a lipid membrane from equilibrium molecular dynamics simulations. The method characterizes the rheology of general interfacial systems embedded in three-dimensional media; we term it molecular interfacial rheology. In our simulations the planar bilayer and surrounding water are confined between solid, parallel walls. Following Onsager's regression hypothesis, membrane and water fluctuations are assumed to relax according to the coupled continuum-mechanical equations governing the confined system---which predict that the membrane transverse velocity autocorrelation function (TVACF) decays exponentially, at a rate set by the membrane and water viscosities. The measured TVACF, however, exhibits damped oscillations followed by a slowly decaying tail. We reconcile these behaviors using the Mori--Zwanzig formalism, and extract the wavevector-dependent membrane viscosity from the time-integral of the TVACF. Results from theory and simulations agree over a decade of wavevectors, and extrapolating to long wavelengths yields shear viscosities ranging from 0.064 to 0.18 pN*us/nm across two representative single-component, fluid-phase bilayers. Our results are corroborated by nonequilibrium simulations where a spatially varying in-plane body force is applied to lipid molecules, thus validating the framework of molecular interfacial rheology.
Remote epitaxy through a monolayer two-dimensional material-covered substrate establishes a crystallographic registry across the van der Waals (vdW) surface that enables the epitaxial growth, lift-off and transfer of single-crystalline films. A central belief in remote epitaxy is that the substrate facilitating the phenomenon must be a material with strong ionicity, as the interatomic electrostatic potential fluctuation in covalent and metallic materials is substantially attenuated by two-dimensional materials. Here, we show remote epitaxy is possible when the substrate is a metallic or covalently bonded material and experimentally demonstrate non-polar remote homo- and heteroepitaxy across a wide range of material systems, including both metals and semiconductors. The achieved non-polar remote interactions are designed and engineered by harnessing substrate conductivity and vicinal surface step-edge density. These findings indicate that remote epitaxy is universal and applicable to ionic, metallic, and covalent materials, expanding its capabilities and stimulating a plethora of new fundamental scientific questions about the mechanism of remote epitaxy.
We investigate how dipolar interactions between microwave-shielded polar molecules structure the self-bound droplets formed under variation of the interaction strength. We identify the transition from droplets to crystals as a finite-size first order transition. With droplet-ring states and transitional supersolid states we predict additional structure in the crystal and droplet phases, respectively. To describe this strongly correlated regime, and in particular the reconfiguration of quantum ground states, we design a variational Monte Carlo framework based on neural quantum states. It is especially suitable to describe ground states and almost degenerate states with very different configurations. Moreover, one can easily determine the superfluid fraction. Our results reveal the sequence of finite-size structures through which dipolar interactions reorganize molecular droplets into crystals.
The confluence of strong electronic correlations and Berry curvature-driven transport constitutes a largely underexplored frontier in quantum materials research, particularly in systems where electronic correlations arise from flat-band physics rather than conventional $f$-electron states. Here, we report an intrinsic Berry curvature driven anomalous Hall transport in the dilute Kondo system Ru$_2$Mn$_{0.5}$Ti$_{0.5}$Ge, a non-$f$-electron Heusler alloy hosting flat bands with van Hove singularity proximate to the Fermi level. The anomalous Hall response is strikingly non-monotonic, violating the Fermi-liquid scaling relation near the magnetic transition and is restored only at low temperature, coincident with the onset of Kondo coherence. \textit{Ab initio} calculations strongly establish that the Berry curvature originates from the Kondo hybridization induced interplay of flat bands and spin-orbit coupling mediated anticrossings at the Fermi level. Our findings demonstrate that Kondo coherence of flat bands bears a profound impact on Berry curvature associated ferroic responses, mandating a rigorous theoretical understanding of anomalous transport in the regime where reciprocal space topology and strong correlations are intrinsically intertwined.
Photons are mutually nearly noninteracting particles, so thermalized photon ensembles are commonly obtained not from direct particle-particle-interactions but rather from contact with matter, which can constitute a reservoir for the photon gas. We develop a theory model for photons in a material-filled (e.g. liquid dye) optical microcavity, with the aim to study the fluctuation properties using a free-energy description for noninteracting photons coupled to a reservoir of material electronic excitations. To begin with, we use a single mode description for the condensate. For a small relative size of the material reservoir, corresponding to the canonical regime, condensate number fluctuations are small, and the derived free energy landscape takes the usual Mexican-hat shaped form such that spontaneous symmetry breaking occurs. In contrast, for a large relative size of the reservoir, corresponding to the grand canonical regime, fluctuations become as large as the average particle number. We show that the resulting free energy landscape acquires a bowl-shaped form, with a single minimum at the origin. Thus, a macroscopic occupation of the ground state (i.e., Bose-Einstein condensation) in the absence of spontaneous symmetry breaking is expected. We also provide a model for the treatment of a photon gas trapped in a box-shaped potential with spatially distributed coupling to a reservoir. The model predicts, for example, a statistically fluctuating pattern of islands with long-range coherence, resembling transient microcondensates.
Establishing a robust and physically interpretable link between static structure and heterogeneous relaxation dynamics remains a fundamental challenge in glass physics. Here, we introduce a weighted pair-entropy descriptor based on the conventional two-body excess entropy. For this, we multiply the integrand used to calculate the excess entropy by a weight function that is directly related to the length scale of the pair correlation function. This multiplication does not affect the contribution of the short range order to the local excess entropy, but allows to take into account the structure present on intermediate distances, i.e., the medium-range order. For a canonical two-dimensional Lennard-Jones glass former, the resulting descriptor exhibits a strong correlation with particle-level dynamical propensity at long times (multiples of the alpha-relaxation time), with a maximum structure--dynamics correlation reaching about 0.9, substantially outperforming the predictive power of the conventional local pair excess entropy. These results demonstrate that incorporating a physically motivated structural length scale into entropy-based descriptors markedly enhances their predictive power while preserving physical interpretability. Our findings provide a simple and general framework for investigating structure--dynamics correlations in glass-forming systems.
We investigate reflected diffusion in a bounded domain subject to two independent, competing stopping mechanisms: an exponentially distributed bulk lifetime of rate $p$ and a surface reaction triggered when the boundary local time exceeds an independent exponential threshold of rate $q$. Denoting by $T$ the stopping time and by $L$ the acquired boundary local time at stopping, we derive their marginal distributions, joint Laplace transform, and complete hierarchy of mixed moments. These statistics are determined by the splitting probability $φ$ that surface reaction occurs before bulk decay. In particular, we establish the identity $p\expect{T}+q\expect{L}=1$ and show that the cumulative risk $pT+qL$ is exponentially distributed with unit rate. We further obtain equivalent representations of $φ$ in terms of the Robin-Laplacian and the generalized Steklov spectra. Explicit results for a three-dimensional ball reveal how competing rates $p,q$ control $φ$ and the $(T,L)$ statistics. Monte Carlo simulations test the universal cumulative-risk law.
Weyl semimetals (WSMs), with their unique topological properties and distinct electronic structure, exhibit intriguing properties when either time-reversal or inversion symmetries are broken. In this work, we consider the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between magnetic impurities in time-reversal symmetry-breaking WSMs. We derive analytical expressions for the full RKKY exchange tensor in arbitrary two-band spinful lattice systems. Our approach reveals both Heisenberg, anisotropic Ising and Dzyaloshinsky-Moriya terms, which can be calculated by energy-integrating real-space Green's functions across the entire Brillouin zone, with the band edge acting as a natural energy cutoff. We apply this framework to study a two-band tight-binding model for a time-reversal symmetry-breaking WSM that interpolates between a Weyl phase with well-separated chiral nodes and a quadratic band-touching semimetal phase. Remarkably, the spatial profile, magnitude, and anisotropic tensor structure of the exchange couplings remain persistent across the node-annihilation transition. This topological robustness reveals that short- and intermediate-range RKKY interactions are mediated by the global, Brillouin-zone-integrated quantum metric of the full valence band rather than being strictly dictated by local low-energy Berry curvature monopoles. These findings demonstrate the necessity of full-band tight-binding formulations when predicting real-space magnetic interactions, providing key insights for electric-field tuning of magnetic anisotropy and constructing realistic models of heavy-fermion and Weyl-Kondo semimetals.
We generalize the Mori-Zwanzig projection operator formalism by introducing an oblique projector using biorthogonal variables to coarse-grain the dynamics of systems far from equilibrium, such as active or living matter. Improving on previously used orthogonal projectors, the oblique projector correctly decomposes the emergent interactions of the coarse-grained particles into systematic and stochastic contributions. We show that the formalism directly connects to established results in non-equilibrium linear response theory and adiabatic perturbation theory. Based on this formalism, we develop a coarse-graining algorithm which we apply to extract the emergent non-equilibrium chiral dynamics of a passive tracer in a chiral active bath. Our theory correctly predicts the odd transport properties and a confinement-induced shift of the chiral frequency of the tracer.
Two conventional Josephson junctions connected in series can reproduce, in the static limit in which the currents through the capacitive and resistive channels are negligible, the current-phase relation of a single effective weak link with tunable transparency. Therefore, the two-junction series can be treated as a single synthetic high-transparency element. Here, we investigate to what extent this mapping remains valid under finite-frequency drive and retaining the junctions' resistive and capacitive terms. The full resistively and capacitively shunted junction equations are compared with an effective synthetic element with tunable transparency that retains the synthetic tunable-transparency current-phase relation together with effective capacitive and dissipative terms, thus reducing the two second order degree of freedom system to a single second order degree of freedom. The resulting single-element dynamics is compared with the complete two-junction system under ac excitation. The agreement is quantified through a normalized root-mean-square error between the full and effective voltage waveforms. A broad low-error region is found at low drive frequency, while pronounced deviations emerge as the drive frequency approaches the relevant plasma-frequency scale and at larger drive amplitudes. The results provide a quantitative dynamical criterion for using the reduced single-element description of a synthetic high-transparency Josephson element in superconducting circuits.
Edge-magnetoplasmon resonators have been proposed as a powerful tool to detect anyons by introducing a quantum point contact into an isolated quantum Hall system probed via radiofrequency radiation. In this paper, we study the effect of a quantum point contact embedded within an edge-magnetoplasmon resonator and how its polarization influences the propagating magnetoplasmonic mode. Combining dc and rf measurements, we unambiguously evidence the signature of both integer ($ν= 1$ and $2$) and fractional quantum Hall states ($ν= 4/3$ and $2/3$) within the radiofrequency transmission signal. Using electrostatic gating, we determine the physical parameters characterizing the electrostatic edge of an AlGaAs/GaAs based two-dimensional electron gas. We extract the dependence of the cavity perimeter with the gate voltage of the quantum point contact and fully characterize the path followed by edge magnetoplasmons in this system. Finally, we provide a geometric model in good agreement with experimental results.
The superconducting gap symmetry of high-pressure bilayer nickelates remains under debate, with weak- and strong-coupling approaches yielding different pairing tendencies. In this work, we investigate how the weak-coupling treatment of electronic states away from the Fermi level influences magnetic fluctuations and superconductivity in La$_3$Ni$_2$O$_7$. We employ a full-spectrum model based on orthonormalized projections of Kohn-Sham states onto local Ni-$e_g$ orbitals, which preserves the density-functional band structure while redistributing spectral weight over a wide energy range. Compared to a low-energy description, this approach yields enhanced interlayer spin fluctuations and a commensurate magnetic instability. Within a spin-fluctuation framework, these features favor a sign-changing $s^\pm$ superconducting state, whereas low-energy models tend to stabilize $d$-wave pairing. Our results suggest that interlayer coupling in full-energy models may play an important role in shaping the predicted pairing symmetry of bilayer nickelates.
For a one-dimensional path graph and a lifted path graph constructed from a duplication of each of its sites, we study how a reversible Markov chain can be perturbed and gradually driven into non-reversibility. The reversible Markov chain has a transition matrix that is diagonalizable and features real-valued eigenvalues and eigenvectors. The left and right eigenvectors form a biorthogonal system. We discuss in concrete examples how the transition matrix of a non-reversible Markov chain may be diagonalizable or non-diagonalizable, and it may have real eigenvalues and complex-conjugate pairs. For a number of steady states (flat, square-wave, wedge, V-shape), we compute eigenvalue spectra on both graphs and discuss the speedup that can be achieved through lifting. We develop a Green's matrix formalism, which we use to compute Kemeny times and mean first-passage times, and which provides valuable information and allows us to interpret the results for the characteristic times.
We present a systematic numerical study of the spin-1 antiferromagnetic Heisenberg model on the triangular lattice in an out-of-plane magnetic field, using Density Matrix Renormalization Group (DMRG) methods. By mapping out the quantum phase diagram as a function of the single-ion anisotropy $D_z$ and magnetic field, we identify distinct dipolar and quadrupolar spin supersolid states, characterized by spontaneous U(1) symmetry breaking with finite spin superfluid stiffness coexisting with longitudinal translational symmetry breaking. At zero field, the dipolar spin supersolid with a 'Y'-type spin configuration persists down to $D_z = 0$, whereas the quadrupolar spin supersolid prevails at large $D_z$. At intermediate fields, the phase diagram is dominated by an up-up-down phase. At high fields below saturation, a quadrupolar spin superfluid emerges in the large-$D_z$ regime, whereas a dipolar spin supersolid with a 'V'-type spin configuration dominates at small $D_z$. These phases are characterized through their order parameters and spin superfluid stiffness using calculations on various system sizes. Furthermore, the dynamical spin structure factor is obtained across the phase diagram, where characteristic spectral signatures of different phases are observed, including the gapless Goldstone mode and the roton-like minima. These features are directly accessible to inelastic neutron scattering experiments. Our results provide a theoretical understanding of the interplay between frustrations, anisotropy, and Zeeman interactions in driving distinct spin supersolid phases in the spin-1 system, which are relevant to various triangular-lattice antiferromagnets such as Na$_2$BaNi(PO$_4$)$_2$ and K$_2$Ni(SeO$_3$)$_2$.
We introduce a transparency-engineered rf-SQUID cell for Kerr-free three-wave-mixing Josephson metamaterials. This design replaces the conventional tunnel-junction element with an effective Josephson element formed by two junctions in series, yielding a non-sinusoidal energy-phase relation that can be used as a tunable nonlinear design resource. We show that the asymmetry between the two series junctions and the applied flux bias provide independent control over the local expansion of the rf-SQUID potential, enabling operating points where the leading quartic Kerr term is suppressed while the cubic nonlinearity remains finite. We derive the corresponding Kerr-free condition, identify the resulting operating ridge in parameter space, and analyze the local stability and passive-matching constraints that bound its physically accessible portion. Our results provide a compact unit-cell design principle for three-wave-mixing Josephson metamaterials and suggest a route toward Kerr-suppressed Josephson traveling-wave parametric amplifiers.
Transformation elasticity provides appropriate material distributions for elastodynamic cloaks but the required stiffness tensors generally violate the minor symmetries of Cauchy elasticity and are difficult to realise using conventional materials. Existing approaches restore these symmetries by modifying the transformed tensor, producing only an approximate cloak. In this work rather than modifying the transformed tensor we seek the best performing cloak within the class of Cauchy materials. We formulate 2D Rayleigh wave carpet cloak design as an optimisation problem governed by partial differential equations. Using a coordinate based neural-field and a differentiable finite element model solver we optimise symmetric stiffness and density fields by minimising wave field distortion. Both single frequency and broadband optimisation are considered, with the broadband model trained over multiple frequencies. Physical realisability is addressed using a database of homogenised microstructures through conditional diffusion, neural-field inverse design, and nearest-neighbour selection. FEM simulations show that the optimised Cauchy design approaches the ideal transformation-based cloak. After projection onto explicit microstructures, the homogenised representation recovers approximately 97% of the defect free reference surface-displacement magnitude, while direct FEM simulation of the fully resolved microstructured geometry recovers approximately 76%
We investigate the slowing down of dynamics in a glass-forming mixture interacting via an inverse-power-law (IPL) potential using a combination of theory and large-scale molecular dynamics simulations. We measure the static pair-correlation function, configurational entropy, inherent-structure energy, and structural relaxation time. We employ a theoretical framework to calculate the structural relaxation time $τ_α$, which is found to be in very good agreement with the simulation results. The theory identifies a local structural order which defines the cooperativity of the relaxation and brings forth a fluctuation induced parameter $ψ( T )$ and a crossover temperature $T_a$ that characterize the density and temperature dependence of the glassy dynamics. Furthermore, we determine a crossover temperature using independent dynamical and thermodynamic criteria and compare with the theoretically predicted crossover temperature $T_a$. Relaxation dynamics is shown to obey density-temperature scaling, similar to thermodynamic properties, in terms of a variable $Γ$ formed by an appropriate combination of density and temperature, characteristic of IPL interactions. Finally, we show that, when the excess thermodynamic and dynamic quantities obtained at different densities are plotted as functions of the reduced temperature $T/T_a$ (or $T_a/T$), the data collapse onto master curves with excellent agreement between theory and simulation. These scaling relations provide a unified description of the thermodynamics and dynamics in IPL systems, enabling the prediction of relaxation behavior over a wide range of densities from data at a single state point.
Heterophase interfaces are increasingly recognized as active participants in plastic deformation, yet quantitative methods for comparing their strain-accommodation capacity remain limited. Here, we present an atomistic framework for quantifying interface-mediated strain accommodation using molecular dynamics simulations. Interface broadening, roughening, and migration are extracted from atomistic trajectories and combined into an Interface Accommodation Index, with a normalized counterpart accounting for differences in initial interface structure. A weighting sensitivity analysis demonstrates that the relative ranking of interfaces is robust to the specific form of the index. The framework is demonstrated for three experimentally motivated interfaces, Al/Al3Ti, Al/Al9M2, and Al9M2/Al3Ti (M = Fe, Co, Ni), under tensile, compressive, and shear loading. Bulk simulations show that Shockley partial dislocations dominate plastic deformation, with alloy chemistry governing the transition toward mixed-character dislocation networks. Tensile loading produces the greatest interface accommodation, while shear produces comparatively limited structural evolution. The Al9M2/Al3Ti interface exhibits both the highest yield resistance and the largest accommodation response, suggesting that intermetallic-intermetallic interfaces can simultaneously sustain load and redistribute strain. These simulations quantify structural accommodation rather than ductility or fracture directly; their connection to macroscopic ductility therefore requires experimental validation. The framework provides a transferable approach for comparing interface accommodation in multiphase alloys and identifies interface chemistry and crystallography as important design variables for damage-tolerant structural materials.
We study steady-state transport in a one-dimensional Fermi-Hubbard chain coupled to particle reservoirs at the boundaries and to local dephasing baths at each site, using the time-evolving block decimation (TEBD) method to solve the Lindblad master equation. In the absence of dephasing, the current exhibits two well-separated maxima as a function of the boundary driving rate, reflecting the distinct charge and spin energy scales of the strongly correlated regime. Upon introducing dephasing, we find two distinct dephasing-induced transport-enhancement regimes, in contrast to the single enhancement previously reported for spinless fermions. Analysis of the non-equilibrium steady state in the Hamiltonian eigenbasis reveals that the two regimes originate from distinct dephasing-induced redistribution processes: the first involves redistribution within the uppermost Hubbard band, while the second involves transitions between Hubbard bands. Our results demonstrate how many-body correlations shape the interplay between coherent driving, dephasing, and quantum Zeno physics in strongly correlated open systems.
We investigated the thermoelectric effects of the Dirac electron system Ag2Te under magnetic field. Our anal-ysis based on the Boltzmann semiclassical model associated the disorder with the unconventional magnetic field responses such as linear magnetoresistance, linear Nernst effect, step-like Nernst effect, and sign change in Nernst effect. The analysis also revealed the impurity band near the Fermi energy. We simultaneously clari-fied the serious impact of the thermal Hall effect on the measurement of the Nernst effect, and we proposed the definitive solution. Our careful measurement and analysis will be the standard for the thermoelectric study under magnetic field.
The microscopic origin of charge density wave formation in quasi one dimensional Ta2NiSe7 remains actively debated, particularly regarding the relative contributions of Fermi surface nesting, electron phonon coupling, and lattice instabilities. Here, we combine temperature and orientation dependent polarized Raman spectroscopy with first principles calculations to uncover the anisotropic electronic and lattice interactions governing the CDW state in Ta2NiSe7. Heat capacity and electrical transport measurements identify an incommensurate CDW transition at 61 K. Raman spectroscopy reveals pronounced in plane anisotropy, with e ph coupling strength along intrachain b axis exceeding five times that along interchain c axis, whereas lattice anharmonicity is enhanced by threefold along c axis. First principles calculations identify Ta2 Se octahedral vibrations and Ta Se electronic states near the Fermi level as dominant channels mediating the anisotropic e ph interaction. Exceptionally strong and directional e ph coupling along b axis establishes lattice driven electronic instability as the primary mechanism underlying CDW modulation and highlights the dominance of intrachain interactions in strongly coupled Ta2NiSe7. Despite this strong coupling, the CDW remains incommensurate, indicating that lattice anharmonicity provides an additional degree of freedom. Enhanced anharmonicity along c axis suggests that anisotropic phonon phonon interactions reshape the free energy landscape and contribute to stabilizing incommensurate phase. These findings reveal a cooperative interplay between anisotropic e ph coupling and lattice anharmonicity in governing CDW formation in low dimensional quantum materials.