SEARCH · Engineering Papers

Results for “Numerical approximations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Revisiting the dark matter interpretation of excess rates in semiconductors

In light of recent results from low-threshold dark matter detectors, we revisit the possibility of a common dark matter origin for multiple excesses across numerous direct detection experiments, with a focus on the excess rates in semiconductor detectors. We explore the interpretation of the low-threshold calorimetric excess rates above 40 eV in the silicon SuperCDMS Cryogenic Phonon Detector and above 100 eV in the germanium EDELWEISS Surface detector as arising from a common but unknown origin, and demonstrate a compatible fit for the observed energy spectra in both experiments, which follow a power law of index α = 3.4 3 - 0.06 + 0.11 . Despite the intriguing scaling of the normalization of these two excess rates with approximately the square of the mass number A 2 , we argue that the possibility of common origin by dark matter scattering via nuclear recoils is strongly disfavored, even allowing for exotic condensed matter effects in an as-yet unmeasured kinematic regime, due to the unphysically large dark matter velocity required to give comparable rates in the different energy ranges of the silicon and germanium excesses. We also investigate the possibility of inelastic nuclear scattering by cosmic ray neutrons, solar neutrinos, and photons as the origin, and quantitatively disfavor all three based on known fluxes of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Simulations Supporting the Development of Northstar's Indirect Beam Parameters Monitoring System

NorthStar Medical Radioisotopes, LLC is planning to produce the important medical radioisotope molybdenum-99 (Mo-99), the parent of technetium-99m (Tc-99m), through photonuclear reactions in molybdenum-100 (Mo-100). In this approach, a target comprising multiple thin disks of enriched molybdenum metal is bombarded with a 40-MeV electron beam. Electrons impinged on the molybdenum target produce bremsstrahlung X-rays that cause the nuclear reaction. Because enriched Mo-100 is expensive, there is a desire to utilize as much beam power as possible to achieve maximum production yield and minimize the size of the target. This requirement leads to very high beam power density (and heat deposition in the target), which creates challenging requirements for the cooling of the target. The critical part of the target is the target window. It separates the high-pressure helium cooled target from the vacuum beamline and the subject of structural and thermal stress. The temperature of the target window is proportional to the energy density deposited by the beam, so it is critical to maintain the desired beam profile on the target window. The feasibility of indirectly monitoring the maximum energy density of the beam on the beam window through beam losses at the main collimator (Collimator) before the production target was verified. A model of the NorthStar beam transport line was constructed for this purpose using MAD-X and Tao/Bmad codes. Beam optics were computed for the standard operational scenario, followed by an investigation involving approximately 400 cases with parameter variations in the last tuning quadrupoles. This was done to assess the correlation between losses in the collimator and the peak energy density on the target. We developed a model to explore the potential application of Optical transition radiation (OTR) for controlling beam parameters in the NorthStar beam delivery system. This model was based on a generic formula derived from the fundamental solution of the inhomogeneous wave equation of the vector potential, and allowed us to consider various surfaces, even those with irregular or random features, using numerical integration. We applied the model to OTR generated by relativistic electrons impacting an Inconel® 718 beam window. We examined cases with different levels of the window’s surface roughness, ranging from 0.5 to 3.0 microns of root square mean (RMS) deviation. The results of the OTR simulations provided distributions of OTR photons that can be used to study the limitations of optical systems for controlling beam parameters.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Encoding trade-offs and design toolkits in quantum algorithms for discrete optimization: coloring, routing, scheduling, and other problems

Challenging combinatorial optimization problems are ubiquitous in science and engineering. Several quantum methods for optimization have recently been developed, in different settings including both exact and approximate solvers. Addressing this field of research, this manuscript has three distinct purposes. First, we present an intuitive method for synthesizing and analyzing discrete (i.e., integer-based) optimization problems, wherein the problem and corresponding algorithmic primitives are expressed using a discrete quantum intermediate representation (DQIR) that is encoding-independent. This compact representation often allows for more efficient problem compilation, automated analyses of different encoding choices, easier interpretability, more complex runtime procedures, and richer programmability, as compared to previous approaches, which we demonstrate with a number of examples. Second, we perform numerical studies comparing several qubit encodings; the results exhibit a number of preliminary trends that help guide the choice of encoding for a particular set of hardware and a particular problem and algorithm. Our study includes problems related to graph coloring, the traveling salesperson problem, factory/machine scheduling, financial portfolio rebalancing, and integer linear programming. Third, we design low-depth graph-derived partial mixers (GDPMs) up to 16-level quantum variables, demonstrating that compact (binary) encodings are more amenable to QAOA than previously understood. We expect this toolkit of programming abstractions and low-level building blocks to aid in designing quantum algorithms for discrete combinatorial problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

A numerically exact description of ultrafast vibrational decoherence in vibration-coupled electron transfer

Broadband pump–probe spectroscopy has been widely used to measure vibrational decoherence associated with the reaction coordinate in photoinduced ultrafast vibration-coupled electron transfer (VCET) reactions. These experiments provide insight into the interplay of intramolecular coordinates along the reaction coordinate. However, a general theoretical foundation for analyzing, and even for explaining rigorously, these data is lacking. In this work, we study vibrational decoherence in a model VCET reaction using the nearly exact time-dependent density matrix renormalization group simulation method. We explore how analyzing the density matrix with quantum information measures can help elucidate the evolution of vibrational coherence in simulations of dynamics. We examine how vibrational coherence is affected by electron transfer on the timescale of approximately 100 femtoseconds. Our results suggest that electron transfer, in the nonadiabatic model, changes the vibrational equilibrium position abruptly—an example of a “quantum quench” event. This explains the concomitant vibrational decoherence. We find that abrupt vibrational decoherence can be mitigated by wavepacket motion occurring on the timescale of the electron transfer.

Science & Technology - Other Topics↗

DOE OSTI2020

Association of residential energy efficiency retrofits with indoor environmental quality, comfort, and health: A review of empirical data

This paper reviews empirical data from evaluations of the influence of residential energy efficiency retrofits on indoor environmental quality conditions and self-reported thermal comfort and health. The data were extracted from 36 studies described in 44 papers plus two reports. Nearly all reviewed studies were performed in Europe or United States. Most studies evaluated retrofits of homes with low-income occupants. Indoor radon and formaldehyde concentrations tended to increase after retrofits that did not add whole-house mechanical ventilation. Study-average indoor concentrations of nitrogen dioxide and volatile organic compounds other than formaldehyde increased and decreased with approximately equal frequency. Average indoor temperatures during winter typically increased after retrofits, usually by less than 1.5 °C. Dampness and mold, usually based on occupant's reports, almost always decreased after retrofits. Subjectively reported thermal comfort, thermal discomfort, non-asthma respiratory symptoms, general health, and mental health nearly always improved after retrofits. For asthma symptoms, the evidence of improvement slightly outweighed the evidence of worsening. There was insufficient evidence to determine whether changes in thermal comfort and health outcomes varied depending on the type of energy efficiency retrofit. The published research has numerous limitations including a lack of data from retrofits in warm-humid climates and minimal data on changes in objective health outcomes. Suggestions for future research are provided.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Growth of Jupiter: Formation in disks of gas and solids and evolution to the present epoch

The formation of Jupiter is modeled via core-nucleated accretion, and the planet's evolution is simulated up to the present epoch. Throughout the phases when the planet acquires most of its heavy-element content, the calculation of solids' accretion accounts for interactions with an evolving disk of planetesimals. The phase of growth from an embryo of a few hundred kilometers in radius until the time when the accretion of gas overtakes solids' accretion was presented by D'Angelo et al., and the same numerical methods are applied here. Those calculations followed the formation for about 4 × 10 5 years, until the epoch when the heavy-element and hydrogen/helium masses were $M_Z ≈ 7.3$ and $M_{XY} ≈ 0.15$ Earth's masses ($M_⊕$), respectively, and $\dot {M} _{XY} ≈ $$\dot {M}_Z$. In this work, the calculation is continued through the phase when $M_{XY}$ grows to equal $M_Z$ , at which age, about 2.4 × 10 6 years, the total mass of the planet is $M_p ≈ 20 M_⊕$. About 9 × 10 5 years later, $M_p$ is approximately $60 M_⊕$ and $M_Z ≈ 16 M_⊕$, three-quarters of which are delivered by planetesimals larger than 10 km in radius. Around this epoch, the contraction of the envelope dictates gas accretion rates a few times 10 –3 $M_⊕$ per year, initiating the regime of disk-limited accretion, whereby the planet can accrete all the gas provided by the disk, and its evolution is therefore tied to disk's evolution. Growth is continued by constructing simplified models of protosolar accretion disks that evolve through viscous diffusion, winds, and accretion on the planet. Jupiter's formation ends after ≈ 3.4–4.2 Myr, depending on the applied disk viscosity parameter, when nebula gas disperses. The young Jupiter is 4.5–5.5 times as voluminous as it is presently and thousands of times as luminous, ~10 –5 $L_⊙$. The heavy-element mass is ≈ 20 $M_⊕$. The evolution proceeds through the cooling and contraction phase, in isolation except for solar irradiation. After 4570 Myr, the age of the solar system, radius and luminosity of the planet are within 10$\%$ of current values, accounting also for uncertainties in the power absorbed from the Sun. During formation, and soon thereafter, the planet exhibits features, e.g., luminosity and effective temperature, that may probe aspects of the latter stages of formation, if observable. These possibly distinctive features, however, seem to disappear within a few tens of Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Impoverishing Roots Will Improve Wheat Yield and Profitability Through Increased Water and Nitrogen Use Efficiencies

We report that more than a 60% increase in crop production is required by the 2050's to feed a growing world population. Understanding how plant functional traits and field management affect crop yields has the potential to improve agricultural productivity, minimize economic and environmental losses, and maximize food security. We explored the influence of winter wheat root characteristics and management on winter wheat growth, yield, and profit using a mechanistic and well-tested ecosystem and crop model, ecosys. We applied and further tested ecosys at an agricultural farm growing winter wheat in Ardmore, Oklahoma, United States. The model accurately predicted observed shoot carbon (R 2 = 0.95), soil moisture (R 2 = 0.67), soil temperature (R 2 = 0.91), and yield (percent error = 17%). Numerical optimization experiments were conducted to explore potential improvements of winter wheat yield and profit by modifying root characteristics, including root radius and root:shoot carbon transfer conductance, and fertilizer inputs. Our results show the potential for simultaneously improving winter wheat yields and profits. The optimum conditions were found to be in the range of root radius between 0.1 and 0.3 mm, carbon transfer conductance between 0.004 and 0.01 h -1 , and the currently applied fertilizer rate of 112 kg ha -1 . Under these conditions, improvements in yields and profits of up to approximately 25% and 110%, respectively, were modeled compared to those under baseline root traits. These improvements were achieved by impoverishing root structures, thereby increasing nutrient allocation to grains. Our results also demonstrate and motivate model structures that integrate the complex network of plant physiology, soil nutrient biogeochemistry, hydrology, and management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Here we present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A Novel Noise-Aware Classical Optimizer for Variational Quantum Algorithms

A key component of variational quantum algorithms (VQAs) is the choice of classical optimizer employed to update the parameterization of an ansatz. It is well recognized that quantum algorithms will, for the foreseeable future, necessarily be run on noisy devices with limited fidelities. Thus, the evaluation of an objective function (e.g., the guiding function in the quantum approximate optimization algorithm (QAOA) or the expectation of the electronic Hamiltonian in variational quantum eigensolver (VQE)) required by a classical optimizer is subject not only to stochastic error from estimating an expected value but also to error resulting from intermittent hardware noise. Model-based derivative-free optimization methods have emerged as popular choices of a classical optimizer in the noisy VQA setting, based on empirical studies. However, these optimization methods were not explicitly designed with the consideration of noise. In this work we adapt recent developments from the “noise-aware numerical optimization” literature to these commonly used derivative-free model-based methods. We introduce the key defining characteristics of these novel noise-aware derivative-free model-based methods that separate them from standard model-based methods. In conclusion, we study an implementation of such noise-aware derivative-free model-based methods and compare its performance on demonstrative VQA simulations to classical solvers packaged in scikit-quant.

classical optimizers↗

DOE OSTI2022

Asymptotic preserving methods for fluid electron-fluid models in the large magnetic field limit with mathematically guaranteed properties (Final Report)

The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Next-Generation Marine Energy Software Needs Assessment

Over the past decade the marine energy industry has continued to grow and evolve, with new concepts and technologies constantly being pursued. Additionally, the field of computing is vastly different today than it was five or ten years ago. By utilizing advanced software and hardware architectures, like graphics processing units as well as parallelization and high-performance computing resources, software can produce higher quality outputs and a higher volume of outputs. These software and hardware resources can enable the marine energy community to exploit computational advancements from other research fields, which can include machine learning, differentiable programming, and controls co-design. Better integration of existing software and development of potential new software is necessary to take advantage of trends in modern computing and respond to the current and future needs of the marine energy community. In order to better understand the existing marine energy software landscape and industry needs, DOE's Water Power Technologies Office (WPTO) tasked Sandia National Laboratories and the National Renewable Energy Laboratory to update the needs assessment by identifying existing software gaps and software needs, and assisting WPTO in planning the next wave of marine energy software development. The proposed effort involved cataloguing and analyzing the available data on existing software related to marine energy. The marine energy software landscape has vastly changed in the last ten years. There are now nearly 230 different software packages utilized by the marine energy sector, compared to a decade ago when the Cardinal Engineering survey identified approximately 40 software packages. In 2012, the marine energy software landscape was captured in two tables, whereas the current marine energy software landscape required development of a software database to collect and categorize software.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗

DOE OSTI2021

Approaching the self-dual point of the sinh-Gordon model

One of the most striking but mysterious properties of the sinh-Gordon model (ShG) is the b → 1/b self-duality of its S-matrix, of which there is no trace in its Lagrangian formulation. Here b is the coupling appearing in the model’s eponymous hyperbolic cosine present in its Lagrangian, cosh(bΦ). In this paper we develop truncated spectrum methods (TSMs) for studying the sinh-Gordon model at a finite volume as we vary the coupling constant. We obtain the expected results for b << 1 and intermediate values of b, but as the self-dual point b = 1 is approached, the basic application of the TSM to the ShG breaks down. We find that the TSM gives results with a strong cutoff Ec dependence, which disappears according only to a very slow power law in Ec. Standard renormalization group strategies — whether they be numerical or analytic — also fail to improve upon matters here. We thus explore three strategies to address the basic limitations of the TSM in the vicinity of b = 1. In the first, we focus on the small-volume spectrum. We attempt to understand how much of the physics of the ShG is encoded in the zero mode part of its Hamiltonian, in essence how ‘quantum mechanical’ vs ‘quantum field theoretic’ the problem is. In the second, we identify the divergencies present in perturbation theory and perform their resummation using a supra-Borel approximate. In the third approach, we use the exact form factors of the model to treat the ShG at one value of b as a perturbation of a ShG at a different coupling. In the light of this work, we argue that the strong coupling phase b > 1 of the Lagrangian formulation of model may be different from what is naïvely inferred from its S-matrix. In particular, we present an argument that the theory is massless for b > 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Spectrally accelerated edge and scrape-off layer gyrokinetic turbulence simulations

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case (Oliveira et al., 2022 [26]) and compare these results with previously validated grid simulations (Ulbl et al., 2023 [25]). The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. As a consequence, the spectral approach reduces the computational cost (CPUh) by at least an order of magnitude, of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

Gyrokinetic↗

DOE OSTI2023

Single-photon absorption and emission from a natural photosynthetic complex

Photosynthesis is generally assumed to be initiated by a single photon from the Sun, which, as a weak light source, delivers at most a few tens of photons per nanometre squared per second within a chlorophyll absorption band1. Yet much experimental and theoretical work over the past 40 years has explored the events during photosynthesis subsequent to absorption of light from intense, ultrashort laser pulses. Here, we use single photons to excite under ambient conditions the light-harvesting 2 (LH2) complex of the purple bacterium Rhodobacter sphaeroides, comprising B800 and B850 rings that contain 9 and 18 bacteriochlorophyll molecules, respectively. Excitation of the B800 ring leads to electronic energy transfer to the B850 ring in approximately 0.7 ps, followed by rapid B850-to-B850 energy transfer on an approximately 100-fs timescale and light emission at 850–875 nm. Using a heralded single photon source20,21 along with coincidence counting, we establish time correlation functions for B800 excitation and B850 fluorescence emission and demonstrate that both events involve single photons. We also find that the probability distribution of the number of heralds per detected fluorescence photon supports the view that a single photon can upon absorption drive the subsequent energy transfer and fluorescence emission and hence, by extension, the primary charge separation of photosynthesis. An analytical stochastic model and a Monte Carlo numerical model capture the data, further confirming that absorption of single photons is correlated with emission of single photons in a natural light-harvesting complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗