Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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 1,243 records · Page 69

Quantifying neptunium oxidation states in nitric acid through spectroelectrochemistry and chemometrics

Controlled-potential in situ thin-layer spectropotentiometry was leveraged to generate visible/near-infrared (VIS/NIR) absorption spectral data sets for the development of chemometric models to quantify Np(III/IV/V/VI) oxidation states in HNO 3 . This technology would be valuable in laboratory studies and when monitoring process solutions to guide feed adjustments for radiochemical separations—the performance of which depends on oxidation state. This approach successfully isolated and stabilized Np species in pure (~99%) oxidation states without compromising solution optical properties. Multivariate curve resolution–alternating least squares models were evaluated to resolve spectral and component concentrations from a scan that sequentially produced Np(VI), Np(V), Np(IV), and Np(III) spectra with mixtures of two valences at a time. Although it provided a useful approximation, the method was not able to quantitively resolve each component likely because of rotational ambiguity. Additionally, partial least squares regression models were built from artificial and electrochemically generated VIS/NIR spectral training sets to study the effect of interionic interactions on spectral characteristics. Models built with true Bi-chemical mixtures of coexisting Np oxidation states and spectra generated from additive combinations of pure end points had similar prediction performance. This methodology can be used to directly quantify Np concentration and the ratio of Np oxidation states and other actinides in remote settings such as hot cells.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Hysteretic temperature sensitivity of wetland CH 4 fluxes explained by substrate availability and microbial activity

Methane (CH 4 ) emissions from wetlands are likely increasing and important in global climate change assessments. However, contemporary terrestrial biogeochemical model predictions of CH 4 emissions are very uncertain, at least in part due to prescribed temperature sensitivity of CH 4 production and emission. While statistically consistent apparent CH 4 emission temperature dependencies have been inferred from meta-analyses across microbial to ecosystem scales, year-round ecosystem-scale observations have contradicted that finding. Here, we show that apparent CH 4 emission temperature dependencies inferred from year-round chamber measurements exhibit substantial intra-seasonal variability, suggesting that using static temperature relations to predict CH 4 emissions is mechanistically flawed. Our model results indicate that such intra-seasonal variability is driven by substrate-mediated microbial and abiotic interactions: seasonal cycles in substrate availability favors CH 4 production later in the season, leading to hysteretic temperature sensitivity of CH 4 production and emission. Our findings demonstrate the uncertainty of inferring CH 4 emission or production rates from temperature alone and highlight the need to represent microbial and abiotic interactions in wetland biogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

An attractive model: simulating fuzzy dark matter with attractive self-interactions

ABSTRACT Fuzzy dark matter (FDM), comprised of ultralight ($m \sim 10^{-22}\,{\rm eV}$) boson particles, has received significant attention as a viable alternative to cold dark matter (CDM), as it approximates CDM on large scales (${\gtrsim}1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-α forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parametrized by an energy decay constant $f \sim 10^{15}~\rm {GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly, and its effects are initially contained to the central region of the halo.

Painter, Connor A. (ORCID:0000000235314806)↗

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]↗

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗

How Well Can CMIP6 Models Represent the Observed Influence of the Pacific and Indian Oceans on the Indian Summer Monsoon Rainfall?

This study evaluates the ability of CMIP6 climate models to simulate the observed effects of tropical Pacific and Indian Ocean sea surface temperature anomalies (SSTAs) on Indian summer monsoon rainfall (ISMR) variability. Using observational data and the large ensemble historical simulations of seven CMIP6 models from 1950 to 2014, we applied a cyclostationary linear inverse model (CS-LIM) to isolate the impacts of tropical Pacific SSTAs, Indian Ocean SSTAs and their interaction on the interannual variability of ISMR. Overall, CMIP6 models well reproduced the observed enhanced (reduced) ISMR variability from Pacific SSTAs (Indian Ocean SSTAs and the Indo-Pacific interaction), but with varying spatial patterns and magnitudes. While CESM2 and E3SM-2-0 showed the best agreement with observations for the effects of Pacific SSTAs and the Indo-Pacific interaction, respectively, CMIP6 models showed mixed results for the impacts from Indian Ocean SSTAs. Composite analysis of ISMR anomalies during the developing phases of pure and co-occurring El Niño-Southern Oscillation (ENSO) and Indian Ocean dipole (IOD) events revealed that the impacts from Pacific SSTAs were captured reasonably well by E3SM-2-0, CESM2, MIROC6, and MPI-ESM1-2-LR, while E3SM-2-0 also showed the best agreement with observations for the effects from the Indo-Pacific interaction. However, all models showed substantial biases in simulating the Indian Ocean SSTA impacts on ISMR, especially for pure El Niño events. Overall, this study provides new insights into how individual CMIP6 models simulate the isolated impacts from the tropical Pacific and Indian Oceans, which has important applications for improving ISMR predictions and interpreting ISMR future projections.

monsoon↗

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗

Analysis and Mitigation of Cascading Outages Using an Interaction Graph Addressing Transient Stability

Cascading outages of power systems pose great threats to system security and reliability, potentially leading to large-scale blackouts. For analysis and mitigation of cascading outages, this paper proposes a transient stability-incorporated interaction graph. This graph statistically quantifies the interactions among line outages and instabilities of generators, which can model propagation paths and patterns of cascading outages. Compared with an interaction graph that only models line outages, this new interaction graph provides important insights on how transient instability occurs along with cascading outages. It also offers more effective strategies for mitigating outage propagation. The proposed interaction graph can be constructed from datasets of historical or simulated cascading events. It is demonstrated on an NPCC 140-bus system with mitigation strategies.

Guo, Zhenping↗

A Novel Measurement of the Anomalous Muon Spin Precession Frequency in the Muon $g-2$ Experiment at Fermilab

The Muon $g-2$ Experiment operated at Fermi National Accelerator Laboratory (FNAL, or Fermilab) between 2018 and 2023 to produce the world's most precise measurement of the muon's \textit{anomalous magnetic moment}, $a_\mu = \frac{g_\mu - 2}{2}$, which expresses the relative deviation in the muon's $g$-factor from a baseline theoretical expectation that $g_\mu = 2$. In the Standard Model of particle physics, $g_\mu > 2$ and hence $a_\mu > 0$ by a calculable amount that depends on all possible interactions between the muon and all other fundamental particles, including any potentially undiscovered interactions beyond the Standard Model. For this reason, measurements of the electron anomaly $a_e$ and later the muon anomaly $a_\mu$ have helped guide the development of the Standard Model since the inception of quantum field theory, and the measured value of $a_\mu$ provides a valuable constraint for new hypotheses that extend the Standard Model. As of 2006, the leading measurement and Standard Model prediction for $a_\mu$ exhibited tension at the level of about three standard deviations, motivating an improved measurement at Fermilab that could test the tension more precisely. The experiment functions by storing a polarized beam of $\mu^+$ in a uniform magnetic field, which simultaneously induces circular motion and spin precession. As the stored muons undergo the Michel decay $\mu^+ \to e^+ + \nu_e + \bar{\nu}_\mu$, mediated by the parity-violating weak interaction, the rest-frame $e^+$ emission direction is correlated with the parent $\mu^+$ spin orientation. Boosting into the laboratory frame encodes this correlation in the decay $e^+$ energy, which is higher when the emission (i.e. $\mu^+$ spin direction) is aligned with the $\mu^+$ momentum, and lower when opposite. Detectors then count the rate of high-energy decay $e^+$, which modulates at the difference between the $\mu^+$ revolution and spin precession frequencies. This observed frequency, called the \textit{anomalous spin precession frequency} $\omega_a$, is directly proportional to $a_\mu$. The extraction of $\omega_a$ proceeds by fitting the time spectrum of detected $e^+$, which requires precise modeling of the $\omega_a$ oscillation as well as any perturbations from beam dynamics and detector acceptance. Using the $\omega_a$ analysis presented in this work, based on Runs 4 -- 6 of the Muon $g-2$ Experiment at Fermilab, we find that $a_\mu = 0.001\,165\,920\,738(162)$ with a relative uncertainty of 139 parts per billion.

Barrett, Tyler [Cornell U.]↗

A hybrid Penman-Monteith and machine learning model for simulating evapotranspiration and its components

Integrating physical processes with machine learning has advanced evapotranspiration (ET) simulation, yet most hybrid models fail to partition total ET into its components: soil evaporation (E) and vegetation transpiration (T). This study introduces Residual Neural Network–Penman–Monteith (RNN-PM), a novel hybrid dual-source ET model designed to overcome this limitation. The model synergizes the physically-based Penman–Monteith framework with three specialized residual neural networks trained to estimate key conductance parameters (canopy conductance, soil surface conductance, and aerodynamic conductance). Furthermore this explicit parameterization allows for the direct partitioning of total ET. Validation at National Ecological Observatory Network (NEON) flux sites using high-frequency partitioned E and T shows that RNN-PM reliably reproduces ET and the transpiration fraction (T/ET). For ET, the model achieves an average Kling–Gupta efficiency (KGE) of 0.89 and a root-mean-square error (RMSE) of 0.55 mm/day; for T/ET, the KGE is 0.87 with an RMSE of 0.06. Furthermore, RNN-PM demonstrates robust generalization, accurately simulating ET and its components well beyond the initial training dataset, even under extreme climatic conditions. This study extended the analysis by comparing the RNN-PM model with seven established dual-source ET models. The results indicate that RNN-PM outperforms both conventional machine learning models and purely physical process-based models in simulating ET components in most cases. Among the purely physical process-based dual-source models, those based on surface temperature decomposition showed improved performance as the leaf area index (LAI) decreased when evaluated against high-frequency ET component datasets. In contrast, the performance of conductance-based dual-source models declined with decreasing LAI. Although purely machine learning-based models can produce relatively accurate simulations of ET components, they often exhibit limited generalization capability, an issue that the RNN-PM model effectively overcomes. Ultimately, the RNN-PM model represents a significant advance in simulating ET components, offering a novel and scalable approach for improving the representation of land–atmosphere interactions in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Understanding the Applicability of Gas/Particle Partitioning Theory to SOA using Novel Chemically-Speciated Partitioning Measurements (Final Report)

The main objective of this research grant project was to investigate the physical-chemical processes and properties of particles and partitioning gases that influence gas/particle (G/P) partitioning and timescales for equilibration, to directly determine the conditions under which equilibrium G/P partitioning theory can be applied to secondary organic aerosol (SOA). SOA constitutes a major component of submicron particles in the atmosphere, which have important implications for air quality, human and ecosystem health, radiation balance, and climate; however SOA formation, evolution and losses are not well understood. Use of novel on-line measurements of chemically-speciated partitioning measurements of both gas-phase and particle-phase components was a central component to advancing these scientific goals. The key approach was to conduct environmental chamber studies with systems of increasing complexity to determine G/P partitioning rates and equilibrium states for chemically-speciated SOA components. The project builds on our previous work on gas-wall partitioning that must be understood in order to study G/P partitioning in chambers, and to help extrapolate the results to be atmospherically relevant. The goal was to employ measurements and modeling to systematically determine when G/P partitioning theory accurately describes the distribution of oxidized organic compounds in the gas and aerosol phases, and when and why it does not (e.g., aerosol phase-state kinetic limitations, non-ideal mixing, phase separation, oligomer formation,), for chemical and physical conditions typically found throughout the continental boundary layer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leptonic scalars and collider signatures in a UV-complete model

We study the non-standard interactions of neutrinos with light leptonic scalars ( ϕ ) in a global ( B – L )-conserved ultraviolet (UV)-complete model. The model utilizes Type-II seesaw motivated neutrino interactions with an SU(2) L -triplet scalar, along with an additional singlet in the scalar sector. This UV-completion leads to an enriched spectrum and consequently new observable signatures. We examine the low-energy lepton flavor violation constraints, as well as the perturbativity and unitarity constraints on the model parameters. Then we lay out a search strategy for the unique signature of the model resulting from the leptonic scalars at the hadron colliders via the processes H ±± → W ± W ± ϕ and H ± → W ± ϕ for both small and large leptonic Yukawa coupling cases. We find that via these associated production processes at the HL-LHC, the prospects of doubly-charged scalar H ±± can reach up to 800 (500) GeV and 1.1 (0.8) TeV at the 2 σ (5 σ ) significance for small and large Yukawa couplings, respectively. A future 100 TeV hadron collider will further increase the mass reaches up to 3.8 (2.6) TeV and 4 (2.7) TeV, at the 2 σ (5 σ ) significance, respectively. We also demonstrate that the mass of ϕ can be determined at about 10% accuracy at the LHC for the large Yukawa coupling case even though it escapes as missing energy from the detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling Phase Equilibrium of Common Sugars Glucose, Fructose, and Sucrose in Mixed Solvents

The industrial processing of sugars and sugar-containing mixtures is gaining widespread use as a form of renewable manufacturing from biomass. To aid in process modeling for design and optimization, a commonly available thermodynamic model is needed that describes the phase equilibrium of these compounds. This work compiles and compares models for solid–liquid and vapor–liquid phase equilibrium from the available data for the representative sugars glucose, fructose, and sucrose in the representative solvents water, methanol, and ethanol, including data for multisugar, multisolvent systems. Additionally, the nonrandom two-liquid (NRTL) model was chosen for these systems because of its widespread use in industry and the availability of parameters for many solvents and cosolutes. The association-NRTL (aNRTL) model was investigated as an improvement for modeling sugars, which may experience a high degree of association because of their many hydroxy groups. Both models accurately capture the data and are able to predict the behavior of multisugar systems with only solute–solvent interaction parameters. The aNRTL model shows an improvement over the baseline NRTL model that is most significant for sucrose, the component with the highest association strength, and least significant for fructose, which has the lowest association strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗

High-dimensional multivariate autoregressive model estimation of human electrophysiological data using fMRI priors

Multivariate autoregressive (MVAR) model estimation enables assessment of causal interactions in brain networks. However, accurately estimating MVAR models for high-dimensional electrophysiological recordings is challenging due to the extensive data requirements. Hence, the applicability of MVAR models for study of brain behavior over hundreds of recording sites has been very limited. Prior work has focused on different strategies for selecting a subset of important MVAR coefficients in the model to reduce the data requirements of conventional least-squares estimation algorithms. Here we propose incorporating prior information, such as resting state functional connectivity derived from functional magnetic resonance imaging, into MVAR model estimation using a weighted group least absolute shrinkage and selection operator (LASSO) regularization strategy. The proposed approach is shown to reduce data requirements by a factor of two relative to the recently proposed group LASSO method of Endemann et al (Neuroimage 254:119057, 2022) while resulting in models that are both more parsimonious and more accurate. The effectiveness of the method is demonstrated using simulation studies of physiologically realistic MVAR models derived from intracranial electroencephalography (iEEG) data. The robustness of the approach to deviations between the conditions under which the prior information and iEEG data is obtained is illustrated using models from data collected in different sleep stages. This approach allows accurate effective connectivity analyses over short time scales, facilitating investigations of causal interactions in the brain underlying perception and cognition during rapid transitions in behavioral state.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Programmable interactions and emergent geometry in an array of atom clouds

Interactions govern the flow of information and the formation of correlations between constituents of many-body quantum systems, dictating phases of matter found in nature and forms of entanglement generated in the laboratory. Typical interactions decay with distance and thus produce a network of connectivity governed by geometry—such as the crystalline structure of a material or the trapping sites of atoms in a quantum simulator. However, many envisioned applications in quantum simulation and computation require more complex coupling graphs including non-local interactions, which feature in models of information scrambling in black hole, and mappings of hard optimization problems onto frustrated classical magnets. Here we describe the realization of programmable non-local interactions in an array of atomic ensembles within an optical cavity, in which photons carry information between atomic spins. By programming the distance dependence of the interactions, we access effective geometries for which the dimensionality, topology and metric are entirely distinct from the physical geometry of the array. As examples, we engineer an antiferromagnetic triangular ladder, a Möbius strip with sign-changing interactions and a treelike geometry inspired by concepts of quantum gravity. The tree graph constitutes a toy model of holographic duality, in which the quantum system lies on the boundary of a higher-dimensional geometry that emerges from measured correlations. Overall, our work provides broader prospects for simulating frustrated magnets and topological phases, investigating quantum optimization paradigms and engineering entangled resource states for sensing and computation.

74 ATOMIC AND MOLECULAR PHYSICS↗