Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical Physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

ORNL Design Engineering Library: Validated Thermoelectric Digital Twins

Radioisotope thermoelectric generators (RTGs) serve a crucial role in supplying thermal and electrical energy for reliable and long-duration power in remote and extreme environments. The thermal energy is provided by the decay of radioisotopes and is converted into electrical energy using the steady-state thermoelectric process by exploiting the Seebeck effect. Maintaining a thermal gradient through the material is necessary to drive current production. Modern predictive multi-physics tools are able to effectively describe the complex and strongly coupled phenomena needed for accurate model-based system engineering efforts needed to develop higher performing RTGs. Experimental validation of model simulations is essential to confirm the accuracy, reliability, and applicability of predictive multiphysics tools. Validation fosters confidence among users, helps meet regulatory requirements, and contributes to the ongoing improvement of simulation techniques, ultimately leading to better, safer, and more efficient designs and processes. Validated models can then be used as digital twins and can be interrogated to understand and quantify the performance gaps between theoretical and physical systems. These insights can be used to build higher performing systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spatially resolved polarization swings in the supermassive binary black hole candidate OJ 287 with first Event Horizon Telescope observations

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage.

OJ 287↗

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparing the DES-SN5YR and Pantheon+ SN cosmology analyses: investigation based on ‘evolving dark energy or supernovae systematics’?

Recent cosmological analyses measuring distances of type Ia supernovae (SNe Ia) and baryon acoustic oscillations (BAO) have all given similar hints at time-evolving dark energy. To examine whether underestimated SN Ia systematics might be driving these results, Efstathiou (2025) compared overlapping SN events between Pantheon+ and DES-SN5YR (20 per cent SNe are in common), and reported evidence for an $\sim$0.04 mag offset between the low- and high-redshift distance measurements of this subsample of events. If this offset is arbitrarily subtracted from the entire DES-SN5YR sample, the preference for evolving dark energy is reduced. In this paper, we show that this offset is mostly due to different corrections for Malmquist bias between the two samples; therefore, an object-to-object comparison can be misleading. Malmquist bias corrections differ between the two analyses for several reasons. First, DES-SN5YR used an improved model of SN Ia luminosity scatter compared to Pantheon+ but the associated scatter-model uncertainties are included in the error budget. Secondly, improvements in host mass estimates in DES-SN5YR also affected SN standardized magnitudes and their bias corrections. Thirdly, and most importantly, the selection functions of the two compilations are significantly different, hence the inferred Malmquist bias corrections. Even if the original scatter model and host properties from Pantheon+ are used instead, the evidence for evolving dark energy from CMB, DESI BAO Year 1 and DES-SN5YR is only reduced from 3.9$\sigma$ to 3.3$\sigma$, consistent with the error budget. Finally, in this investigation, we identify an underestimated systematic uncertainty related to host galaxy property uncertainties, which could increase the final DES-SN5YR error budget by 3 per cent. In conclusion, we confirm the validity of the published DES-SN5YR results.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA's Future X-ray Missions: From Constellation-X to Generation-X

Among the most important topics in modern astrophysics are the formation and evolution of supermassive black holes in concert with galaxy bulges, the nature of the dark energy equation of state, and the self-regulating symmetry imposed by both stellar and AGN feedback. All of these topics are readily addressed with observations at X-ray wavelengths. NASA's next major X-ray observatory is Constellation-X, which is being developed to perform spatially resolved high-resolution X-ray spectroscopy. Con-X will directly measure the physical properties of material near black holes' last stable orbits and the absolute element abundances and velocities of hot gas in clusters of galaxies. The Con-X mission will be described, as well as its successor, Generation-X (anticipated to fly approx.1 decade after Con-X). After describing these missions and their driving science areas, the talk will focus on areas in which Chandra observing programs may enable science with future X-ray observatories. These areas include a possible ultra-deep Chandra imaging survey as an early Universe pathfinder, a large program to spatially resolve the hot intracluster medium of massive clusters to aid dark energy measurements, and possible deep spectroscopic observations to aid in preparatory theoretical atomic physics work needed for interpreting Con-X spectra.

Hornschemeier, A.↗

Reconnection-driven Flares in M87*: Proton–Synchrotron-powered GeV Emission

Magnetic reconnection in current layers that form intermittently in radiatively inefficient accretion flows onto black holes is a promising mechanism for particle acceleration and high-energy emission. It has been recently proposed that such layers, arising during flux eruption events, can power the rapid TeV flares observed from the core of M87. In this scenario, inverse-Compton scattering of soft radiation from the accretion flow by energetic electron–positron pairs produced near the reconnection layer was suggested as the primary emission mechanism. However, detailed calculations show that radiation from pairs alone cannot account for the GeV emission detected by the Fermi observatory. In this work, we combine analytic estimates with 3D radiative particle-in-cell simulations of pair–proton plasmas to show that the GeV emission can be naturally explained by synchrotron radiation from protons accelerated in the current sheet. Although the exact proton content of the layer is uncertain, our model remains robust across a broad range of proton-to-pair number density ratios. While protons are subdominant in number compared to pairs, our simulations demonstrate that they can be accelerated more efficiently, leading to a self-regulated steady state in which protons dominate the energy budget. Ultimately, proton synchrotron emission accounts for approximately 5%–20% of the total dissipation power. The majority is radiated as MeV photons via pair synchrotron emission, with a smaller fraction emitted as TeV photons through inverse-Compton scattering.

Active galactic nuclei↗

Measurements of Quasar Proximity Zones with the Lyα Forest of DESI Y1 Quasars

Abstract The intergalactic medium (IGM) around a quasar is shaped by its dense environment and by its excess ionizing radiation, which form a “quasar proximity zone” whose size and anisotropy depend on the quasar’s halo mass, luminosity, age, and radiation geometry. Using over 10,000 quasar pairs from the Dark Energy Spectroscopic Instrument (DESI) Year 1 data, with projected comoving separations r ⊥ < 2 h −1 Mpc, we investigate how the proximity zone of foreground quasars at z ∼ 2–3.5 affects Lyα absorption in their background quasars. The large DESI sample enables unprecedented precision in measuring this “transverse proximity” effect, allowing a detailed investigation of the signal’s dependence on the projected separation of quasar pairs and the luminosity of the foreground quasar. We find that enhanced gas clustering near quasars dominates over their ionizing effect, leading to stronger absorption on neighboring sightlines. Under the assumption that quasar ionizing luminosity is isotropic and steady, we infer the IGM overdensity profile in the vicinity of quasars, finding overdensities as high as Δ ∼ 10 at comoving distance ∼1 h −1 Mpc from the most luminous systems. Surprisingly, however, we find no significant dependence of the proximity profile on the luminosity of the foreground quasar. This lack of luminosity dependence could reflect a cancellation between higher ionizing flux and higher gas overdensity, or it could indicate that quasar emission is highly time-variable or anisotropic, so that the observed luminosity does not trace the ionizing flux on nearby sightlines.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey year 3 results: Cosmological constraints from cluster abundances, weak lensing, and galaxy clustering

Galaxy clusters provide a unique probe of the late-time cosmic structure and serve as a powerful independent test of the Λ⁢CDM model. This work presents the first set of cosmological constraints derived with ∼16,000 optically selected redMaPPer clusters across nearly 5000 deg 2 using DES year 3 datasets. Our analysis leverages a consistent modeling framework for galaxy cluster cosmology and DES-Y3 joint analyses of galaxy clustering and weak lensing (3 × 2⁢pt), ensuring direct comparability with the DES-Y3 3 × 2⁢pt analysis. Here, we obtain constraints of 𝑆 8 = 0.864 ± 0.035 and Ω m = $0.26⁢5^{+0.019}_{−0.031}$ from the cluster-based data vector. We find that cluster constraints and 3 × 2⁢pt constraints are consistent under the Λ⁢CDM model with a posterior predictive distribution (PPD) value of 0.53. The consistency between clusters and 3 × 2⁢pt provides a stringent test of Λ⁢CDM across different mass and spatial scales. Jointly analyzing clusters with 3 × 2⁢pt further improves cosmological constraints, yielding 𝑆 8 = $0.81⁢1^{+0.022}_{−0.020}$ and Ω m = $0.29⁢4^{+0.022}_{−0.033}$, a 24% improvement in the Ω m − 𝑆 8 figure of merit over 3 × 2⁢pt alone. Moreover, we find no significant deviation from the Planck CMB constraints with a probability to exceed (PTE) value of 0.6, significantly reducing previous 𝑆 8 tension claims. Finally, combining DES 3 × 2⁢pt, DES clusters, and Planck CMB places an upper limit on the sum of neutrino masses of ∑𝑚 𝜈 < 0.26 eV at 95% confidence under the Λ⁢CDM model. These results establish optically selected clusters as a key cosmological probe and pave the way for cluster-based analyses in upcoming stage-IV surveys such as LSST, Euclid, and Roman.

Abbott, T. M. C. [NSF’s National Optical-Infrared ↗

Cosmological analysis of the DESI DR1 Lyα 1D power spectrum

We present the cosmological analysis of the one-dimensional Lyman-α flux power spectrum from the first data release of the Dark Energy Spectroscopic Instrument (DESI). We capture the dependence of the signal on cosmology and intergalactic medium physics using an emulator trained on a cosmological suite of hydrodynamical simulations, and we correct its predictions for the impact of astrophysical contaminants and systematics, many of these not considered in previous analyses. We employ this framework to constrain the amplitude and logarithmic slope of the linear matter power spectrum at k ★ = 0.009 km -1 s and redshift z = 3, obtaining Δ 2 ★ = 0.379 ± 0.032 and n ★ = -2.309 ± 0.019 https://github.com/igmhub/cobaya_lya_p1d. The robustness of these constraints is validated through the analysis of mocks and a large number of alternative data analysis variations, with cosmological parameters kept blinded throughout the validation process. We then combine our results with constraints from DESI BAO and temperature, polarization, and lensing measurements from Planck, ACT, and SPT-3G to set constraints on ΛCDM extensions. While our measurements do not significantly tighten the limits on the sum of neutrino masses from the combination of these probes, they sharpen the constraints on the effective number of relativistic species, N eff = 3.02 ± 0.10, the running of the spectral index, α s = 0.0014 ± 0.0041, and the running of the running, β s = -0.0006 ± 0.0048, by a factor of 1.18, 1.27, and 1.90, respectively. We conclude by outlining the improvements needed to fully reach the level of confidence implied by these uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Separation of infrared and bulk in thermal QCD

A new thermal regime of QCD, featuring decoupled scale-invariant infrared glue, has been proposed to exist both in pure-glue (N f =0) and “real-world” (N f =2+1 at physical quark masses) QCD. In this IR phase, elementary degrees of freedom flood the infrared, forming a distinct component independent from the bulk. This behavior necessitates non-analyticities in the theory. In pure-glue QCD, such non-analyticities have been shown to arise via Anderson-like mobility edges in Dirac spectra (λ IR = 0, ±λ A ≠ 0), as manifested in the dimension function d IR (λ). Here, we present the first evidence, based on lattice QCD calculation at a=0.105 fm, that this mechanism is also at work in real-world QCD, thus supporting the existence of the proposed IR regime in nature. An important aspect of our results is that, while at T = 234 MeV we find a dimensional jump between zero modes and lowest near-zero modes very close to unity (d IR = 3 to d IR ≃ 2), similar to the IR phase of pure-glue QCD, at T = 187 MeV we observe a continuous λ-dependence. This suggests that thermal states just above the chiral crossover are non-analytically (in T) connected to thermal state at T = 234 MeV, supporting the key original proposition that the transition into the IR regime occurs at a temperature strictly above the chiral crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Backup Program of the Dark Energy Spectroscopic Instrument’s Milky Way Survey

The Milky Way Backup Program (MWBP), a survey currently underway with the Dark Energy Spectroscopic Instrument (DESI) on the Nicholas U. Mayall 4 m Telescope, works at the margins of the DESI Main surveys to obtain spectra of millions of additional stars from the Gaia catalog. Efficiently utilizing times between ∼12° and 18° twilight and poor weather conditions, the MWBP extends the range of stellar sources studied to both brighter magnitudes and lower Galactic latitude and declination than the stars studied in DESI’s Main Milky Way Survey. While the MWBP prioritizes candidate giant stars selected from the Gaia catalog (using color and parallax criteria), it also includes an unbiased sample of bright stars (i.e., 11.2 ≲ G < 16 mag) as well as fainter sources (to G ≲ 19 mag). As of 2025 March 1, the survey had obtained spectra of ∼7 million stars, approximately 1.2 million of which are included in the DESI Data Release 1. The DESI spectra cover the wavelength range from 3600 to 9800 Å at a resolution λ/Δλ varying from 2000 to 5000. The full survey, when completed, will cover an area of more than 21,000 deg 2 and include approximately 10 million Gaia sources, roughly equal to the number of stellar spectra obtained through the DESI Main Survey, while only utilizing ≈9% of all DESI observing time. This paper provides an overview of the MWBP, describing the target selection, observing strategy, and an introduction to the resulting data.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of the DESI DR2 measurements of baryon acoustic oscillations from galaxies and quasars

The Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2) galaxy and quasar clustering data represents a significant expansion of data from Data Release 1 (DR1), providing improved statistical precision in baryon acoustic oscillation (BAO) constraints across multiple tracers, including bright galaxies, luminous red galaxies, emission line galaxies, and quasars. In this paper, we validate the BAO analysis of DR2. We present the results of robustness tests on the blinded DR2 data and, after unblinding, consistency checks on the unblinded DR2 data. All results are compared with those obtained from a suite of mock catalogs that replicate the selection and clustering properties of the DR2 sample. We confirm the consistency of DR2 BAO measurements with DR1 while achieving a reduction in statistical uncertainties due to the increased survey volume and completeness. The combined BAO precision, including both statistical and systematic errors, improves from ∼0.52% in DR1 to 0.30% in DR2—a factor of 1.7 gain. We assess the impact of analysis choices, including different data vectors (correlation function vs power spectrum), modeling approaches and systematics treatments, and an assumption of the Gaussian likelihood, finding that our BAO constraints are stable across these variations and assumptions with a few minor refinements to the baseline setup of the DR1 BAO analysis. We summarize a series of pre-unblinding tests that confirmed the readiness of our analysis pipeline, the final systematic errors, and the DR2 BAO analysis baseline. The successful completion of these tests led to the unblinding of the DR2 BAO measurements, ultimately leading to the DESI DR2 cosmological analysis, with their implications for the expansion history of the Universe and the nature of dark energy presented in the DESI key paper (companion paper).

79 ASTRONOMY AND ASTROPHYSICS↗

Parametric amplification of spin-motion coupling in three-dimensional trapped-ion crystals

Three-dimensional (3D) crystals offer a route to scaling up trapped-ion systems for quantum sensing and quantum simulation applications; however, engineering coherent spin-motion couplings and effective spin-spin interactions in large crystals poses technical challenges associated with decoherence and prolonged timescales to generate appreciable entanglement. Here, we explore the possibility of speeding up these interactions in 3D crystals via parametric amplification. For this purpose, we derive a general Hamiltonian for the parametric amplification of spin-motion coupling that is broadly applicable to normal modes with motion transverse to or along the spatial extent of the crystal. Unlike in lower-dimensional crystals, we find that the ability to faithfully (uniformly) amplify the spin-spin interactions in 3D crystals depends on the physical implementation of the spin-motion coupling. We consider the light-shift gate, and the so-called phase-insensitive and phase-sensitive Mølmer-Sørensen (MS) gates, and we find that only the phase-sensitive MS gate can be faithfully amplified in general 3D crystals. We discuss a situation where nonuniform amplification can be advantageous. We also reconsider the effect of counter-rotating terms on parametric amplification and find that they are not as detrimental as previous studies suggest.

Optical parametric oscillators & amplifiers↗

Celestial Quantum Error Correction. Part II. From qudits to celestial CFT

A holographic CFT description of asymptotically flat spacetimes inherits vacuum degeneracies and IR divergences from its gravitational dual. We devise a Quantum Error Correcting (QEC) framework to encode both effects as correctable fluctuations on the CFT dual. The framework is physically motivated by embedding a chain of qudits in the so-called Klein spacetime and then taking a continuum 𝑁 → ∞ limit. At finite 𝑁 the qudit chain 1) enjoys a discrete version of celestial symmetries and 2) supports a Gottesman-Kitaev-Preskill (GKP) code. The limit results in hard states with quantized BMS hair in the celestial torus forming the logical subspace, robust under errors induced by soft radiation. Technically, the construction leverages the recently studied 𝑤 1+∞ hierarchy of soft currents and its realization from a sigma model in twistor space.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey Year 6 results: Redshift calibration of the MagLim++ lens sample

In this work, we derive and calibrate the redshift distribution of the MagLim++ lens galaxy sample used in the Dark Energy Survey Year 6 (DES Y6) 3 x 2pt cosmology analysis. The 3 x 2pt analysis combines galaxy clustering from the lens galaxy sample and weak gravitational lensing. The redshift distributions are inferred using the SOMPZ method - a Self-Organizing Map framework that combines deep-field multi-band photometry, wide-field data, and a synthetic source injection ( B alrog) catalog. Key improvements over the DES Year 3 (Y3) calibration include a noise-weighted SOM metric, an expanded Balrog catalogue, and an improved scheme for propagating systematic uncertainties, which allows us to generate O(10 8 ) redshift realizations that collectively span the dominant sources of uncertainty. These realizations are then combined with independent clustering-redshift measurements via importance sampling. The resulting calibration achieves typical uncertainties on the mean redshift of 1-2%, corresponding to a 20-30% average reduction relative to DES Y3. We compress the n(z) uncertainties into a small number of orthogonal modes for use in cosmological inference. Marginalizing over these modes leads to only a minor degradation in cosmological constraints. Here, this analysis establishes the MagLim++ sample as a robust lens sample for precision cosmology with DES Y6 and provides a scalable framework for future surveys.

dark energy↗

Atmospheric ice nucleation

Atmospheric ice nucleation is crucial for global precipitation and affects the structure, lifetime and reflectivity of clouds, thereby impacting climate. Ice nucleates in various ways from aerosol particles, termed ice-nucleating particles, over an extensive temperature and humidity range. Quantifying the kinetic and thermodynamic regimes of nucleation is necessary to relate fundamental physics to theoretically based predictions of ice formation for implementation in cloud and climate models. We review how the molecular picture of ice nucleation has advanced in recent years and consequential impacts on the interpretation and parameterization of ice nucleation. Advances include the role of interfacial free energy and pressure on ice nucleation rates, mobility regions of water that generate the critical ice nucleus, classical and non-classical pathways of nucleation, the type of ice polymorph that forms, the impact of solutes on freezing and the role of nanopores as surface features promoting ice nucleation. We also introduce currently debated and evaluated freezing parameterizations for application in model environments. Finally, we outline what we believe are the current needs for improving predictive understanding of ice nucleation.

54 ENVIRONMENTAL SCIENCES↗

Dark Energy Survey Year 6 results: Clustering redshifts and importance sampling of self-organized-maps 𝑛⁡(𝑧) realizations for 3 × 2 ⁢pt samples

This work is part of a series establishing the redshift framework for the 3 × 2 ⁢pt analysis of the Dark Energy Survey Year 6 (DES Y6). For DES Y6, photometric redshift distributions are estimated using self-organizing maps (SOMs), calibrated with spectroscopic and many-band photometric data. To overcome limitations from color-redshift degeneracies and incomplete spectroscopic coverage, we enhance this approach by incorporating clustering-based redshift constraints (clustering-z, or WZ) from angular cross-correlations with BOSS and eBOSS galaxies and eBOSS quasar samples. We define a WZ likelihood and apply importance sampling to a large ensemble of SOM-derived 𝑛⁡(𝑧) realizations, selecting those consistent with the clustering measurements to produce a posterior sample for each lens and source bin. The analysis uses angular scales corresponding to 1.5–5 Mpc to optimize signal-to-noise ratio while mitigating modeling uncertainties and marginalizes over redshift-dependent galaxy bias and other systematics informed by the N-body simulation CARDINAL . While a sparser spectroscopic reference sample limits WZ constraining power at 𝑧 >1.1, particularly for source bins, we demonstrate that combining SOM with WZ improves redshift accuracy and enhances the overall cosmological constraining power of DES Y6. As a result, we estimate an improvement in 𝑆 8 of approximately 10% for cosmic shear and 3 ×2⁢pt analysis, primarily due to the WZ calibration of the source samples.

Cosmological parameters↗

Charge correlations and magnetoelastic coupling in intercalated transition metal dichalcogenides

The large van der Waals gap in transition metal dichalcogenides (TMDs) offers an avenue to tune the ground state of 2D materials through the intercalation of magnetic atoms. Here, we investigate the charge correlations in Fe 1/3 ⁢TaS 2 , Co 1/3 ⁢TaS 2 , and Fe 0.35 ⁢NbS 2 by combining angle-resolved photoemission spectroscopy (ARPES), x-ray scattering, magnetometry, and density functional theory (DFT). We find that, while short-range charge fluctuations develop in Ta-based compounds, Fe 0.35 ⁢NbS 2 exhibits long-range charge order which is strongly coupled with magnetic order and tunable by external magnetic field. Our electronic structure analysis reveals that intercalation reconstructs the Fermi surface via charge transfer and band renormalization, yet does not generate the nesting conditions compatible with the observed ordering vectors. Complementary phonon calculations further exclude a conventional electron-phonon origin of charge order. Together, these results establish magnetoelastic coupling as the dominant mechanism behind charge ordering in Fe 0.35 ⁢NbS 2 and highlight the contrasting role of Nb and Ta hosts in stabilizing correlated ground states in intercalated TMDs.

36 MATERIALS SCIENCE↗