Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “States”

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 487 records · Page 27

Probing autoionization decay lifetimes of the 4d−16ℓ core-excited states in xenon using attosecond noncollinear four-wave-mixing spectroscopy

The decay of core-excited states is a sensitive probe of autoionization dynamics and correlation effects in many-electron systems, occurring on ultrafast timescales. Xenon, with its dense manifold of autoionizing resonances that can be coupled with near-infrared light, provides a platform to investigate these processes. In this work, the autoionization decay lifetimes of 4d-16ℓ (ℓ = s, p, d, …) core-excited states in xenon atoms are probed with extreme ultraviolet (XUV) attosecond noncollinear four-wave-mixing (FWM) spectroscopy. The 4d{5/2,3/2}-16p XUV-bright states (optically dipole-allowed) exhibit decay lifetimes of ∼6 fs, which is consistent with spectator-type decay. In contrast, the 4d{5/2,3/2}-16s and 4d{5/2,3/2}-16d XUV-dark states (optically dipole forbidden) show longer decay lifetimes of ∼20 fs. Photoionization calculations confirm that all core-hole states with 4d character should decay via spectator channels in ≤6 fs, suggesting that the apparently longer dark-state decay times arise from an alternative mechanism. A few-level simulation of the FWM process shows that the inclusion of a nearby, longer-lived dark state can mimic the experimental FWM signal, suggesting population cycling with a second electronic state with non-4d character. Ab initio calculations support the presence of such multielectron excited states in the 60-70 eV range. These results demonstrate that FWM signals can encode coupled-state dynamics when probing complex systems, highlighting the importance of combining theoretical and experimental approaches to disentangle accurate core-level decay pathways and lifetimes.

Puskar, Nicolette G↗

Incommensurate and commensurate antiferromagnetic states in Ca Mn 2 As 2 and Sr Mn 2 As 2 revealed by As 75 NMR

Here, we carried out 75 As nuclear magnetic resonance (NMR) measurements on the trigonal CaMn 2 As 2 and SrMn 2 As 2 insulators exhibiting antiferromagnetic (AFM) ordered states below Néel temperatures T N = 62 and 120 K, respectively. In the paramagnetic state above T N , typical quadrupolar-split 75 As NMR spectra were observed for both systems. The 75 As quadrupolar frequency ν Q for CaMn 2 As 2 decreases with decreasing temperature, while ν Q for SrMn 2 As 2 increases, showing an opposite temperature dependence. In the AFM state, the relatively sharp and distinct 75 As NMR lines were observed in SrMn 2 As 2 and the NMR spectra were shifted to lower fields for both magnetic fields H ∥ c axis and H ∥ a b plane, suggesting that the internal fields B int at the As site produced by the Mn ordered moments are nearly perpendicular to the external magnetic field direction. No obvious distribution of B int was observed in SrMn 2 As 2 , which clearly indicates a commensurate AFM state. In sharp contrast to SrMn 2 As 2 , broad and complex NMR spectra were observed in CaMn 2 As 2 in the AFM state, which clearly shows a distribution of B int at the As site, indicating an incommensurate state. From the analysis of the characteristic shape of the observed spectra, the AFM state of CaMn 2 As 2 was determined to be a two-dimensional incommensurate state where Mn ordered moments are aligned in the ab plane. A possible origin for the different AFM states in the systems was discussed. Both CaMn 2 As 2 and SrMn 2 As 2 show very large anisotropy in the nuclear spin-lattice relaxation rate 1/T 1 in the paramagnetic state. 1/T 1 for H ∥ ab is much larger than that for H ∥ c, indicating strong anisotropic AFM spin fluctuations in both compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Searching for resonance states in Ne 22 ( p , γ ) Na 23

Globular clusters show strong correlations between different elements, such as the well-known sodium-oxygen anticorrelation. One of the main sources of uncertainty in this anticorrelation is the 22 Ne(p,γ) 23 Na reaction rate, due to the possible influence of an unobserved resonance state at E x = 8862 keV (E r,c.m. =68 keV). The influence of two higher-lying resonance states at E x = 8894 and 9000 keV has already been ruled out by direct 22 Ne(p,γ) 23 Na measurements. To study excited states in 23 Na above the proton threshold to determine if the unconfirmed resonance states in 23 Na exist. The non-selective proton inelastic scattering reaction at low energies was used to search for excited states in 23 Na above the proton threshold. Protons scattered from various targets were momentum-analysed in the Q3D magnetic spectrograph at the Maier-Leibnitz Laboratorium, Munich, Germany. The resonance states previously reported at E x = 8862, 8894 and 9000 keV in other experiments were not observed in the present experiment at any angle. This result, combined with other non-observations of these resonance states in most other experiments, results in a strong presumption against the existence of these resonance states. Here, the previously reported resonance states at E x = 8862, 8894 and 9000 keV are unlikely to exist and should be omitted from future evaluations of the 22 Ne(p,γ) 23 Na reaction rates. Indirect studies using low-energy proton inelastic scattering are a simple and yet exceptionally powerful tool in helping to constrain astrophysical reaction rates by providing non-selective information of the excited states of nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Interplay of the Spin Density Wave and a Possible Fulde-Ferrell-Larkin-Ovchinnikov State in CeCoIn 5 in Rotating Magnetic Field

The d-wave superconductor CeCoIn 5 has been proposed as a strong candidate for supporting the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state near the low-temperature boundary of its upper critical field. Neutron diffraction, however, finds spin-density-wave (SDW) order in this part of the phase diagram for field in the a-b plane, and evidence for the SDW disappears as the applied field is rotated toward the tetragonal c axis. It is important to understand the interplay between the SDW and a possible FFLO state in CeCoIn 5 , as the mere existence of an SDW does not necessarily exclude an FFLO state. In this work, based on a model constructed on the basis of available experiments, we show that an FFLO state competes with an SDW phase. The SDW state in CeCoIn 5 is stabilized when the field is directed close to the a-b plane. When the field is rotated toward the c axis, the FFLO state emerges, and the SDW phase disappears. In the FFLO state, the nodal planes with extra quasiparticles (where the superconducting order parameter is zero) are perpendicular to the field, and in the SDW phase, the quasiparticle density of states is reduced. We test this model prediction by measuring heat transported by normal quasiparticles in the superconducting state. As a function of field, we observe a reduction of thermal conductivity for field close to the a-b plane and an enhancement of thermal conductivity when field is close to the c axis, consistent with theoretical expectations. Our modeling and experiments, therefore, indicate the existence of the FFLO state when field is parallel to the c axis.

36 MATERIALS SCIENCE↗

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings v = 0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at v = 0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023)., Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).] and find the ordered states to have lower energy at v = 0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at v = 0,±1,±2 down to ~2 K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling v = ±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Lastly, we study the effect of a C 3⁢z -rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C 3⁢z symmetry) can be stabilized at sufficiently large strain at v = 0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Constant-Depth Preparation of Matrix Product States with Adaptive Quantum Circuits

Adaptive quantum circuits, which combine local unitary gates, midcircuit measurements, and feedforward operations, have recently emerged as a promising avenue for efficient state preparation, particularly on near-term quantum devices limited to shallow-depth circuits. Matrix product states (MPS) comprise a significant class of many-body entangled states, efficiently describing the ground states of one-dimensional gapped local Hamiltonians and finding applications in a number of recent quantum algorithms. Recently, it has been shown that the Affleck-Kennedy-Lieb-Tasaki state—a paradigmatic example of an MPS—can be exactly prepared with an adaptive quantum circuit of constant depth, an impossible feat with local unitary gates alone due to its nonzero correlation length [Smith , PRX Quantum 4, 020315 (2023)]. In this work, we broaden the scope of this approach and demonstrate that a diverse class of MPS can be exactly prepared using constant-depth adaptive quantum circuits, outperforming theoretically optimal preparation with unitary circuits. We show that this class includes short- and long-ranged entangled MPS, symmetry-protected topological (SPT) and symmetry-broken states, MPS with finite Abelian, non-Abelian, and continuous symmetries, resource states for MBQC, and families of states with tunable correlation length. Moreover, we illustrate the utility of our framework for designing constant-depth sampling protocols, such as for random MPS or for generating MPS in a particular SPT phase. We present sufficient conditions for particular MPS to be preparable in constant time, with global on-site symmetry playing a pivotal role. Altogether, this work demonstrates the immense promise of adaptive quantum circuits for efficiently preparing many-body entangled states and provides explicit algorithms that outperform known protocols to prepare an essential class of states. Published by the American Physical Society 2024

Smith, Kevin C. (ORCID:0000000223971518)↗

Numerical Circuit Synthesis and Compilation for Multi-State Preparation

Near-term quantum computers have significant error rates and short coherence times, so compilation of circuits to be as short as possible is essential. Two types of compilation problems are typically considered: circuits to prepare a given state from a fixed input state, called 'state preparation'; and circuits to implement a given unitary operation, for example by 'unitary synthesis'. In this paper we solve a more general problem: the transformation of a set of m states to another set of m states, which we call 'multi-state preparation'. State preparation and unitary synthesis are special cases; for state preparation, m=1, while for unitary synthesis, m is the dimension of the full Hilbert space. We generate and optimize circuits for multi-state preparation numerically. In cases where a top-down approach based on matrix decompositions is also possible, our method finds circuits with substantially (up to 40 %) fewer two-qubit gates. We discuss possible applications, including efficient preparation of macroscopic superposition ('cat') states and synthesis of quantum channels.

Szasz, Aaron↗

A FFT Method for the Quasiclassical Selection of Initial Ro-Vibrational States of Triatomic Molecules

This paper describes the use of an exact fast Fourier transform (FFT) method to prepare specified vibrational-rotational states of triatomic molecules. The method determines the Fourier coefficients needed to describe the coordinates and momenta of a vibrating-rotating triatomic molecule. Once the Fourier coefficients of a particular state are determined, it is possible to easily generate as many random sets of initial cartesian coordinates and momenta as desired. All the members of each set will correspond to the particular vibrational-rotational state selected. For example, in the case of the ground vibrational state of a non-rotating water molecule, the calculated actions of 100 sets of initial conditions produced actions within 0.001 h(bar) of the specified quantization values and energies within 5 cm(sup -1) of the semiclassical eigenvalue. The numerical procedure is straightforward for states in which all the fundamental frequencies are independent. However for states for which the fundamental frequencies become commensurate (resonance states), there are additional complications. In these cases it is necessary to determine a new set of "fundamental" frequencies and to modify the quantization conditions. Once these adjustments are made, good results are obtained for resonance states. The major problems are in labelling the large number of Fourier coefficients and the presence of regions of chaotic motion. Results are presented for the vibrational states of H2O and HCN and the ro-vibrational states of H2O.

Eaker, Charles W.↗

Alternating and Gaussian Fermionic Isometric Tensor Network States

Isometric tensor networks in two dimensions enable efficient and accurate study of quantum many-body states, yet the effect of the isometric restriction on the represented quantum states is not fully understood. We address this question in two main contributions. First, we introduce an improved variant of isometric tensor network states (isoTNS) in two dimensions, where the isometric arrows on the columns of the network alternate between pointing upward and downward; hence the name alternating isometric tensor network states. Second, we introduce a numerical tool—the isometric Gaussian fermionic TNS (isoGfTNS)—that incorporates isometric constraints into the framework of Gaussian fermionic tensor network states. We demonstrate in numerous ways that alternating isoTNSs represent many-body ground states of two-dimensional quantum systems significantly better than the original isoTNSs. First, we show that the entanglement in an isoTNS is mediated along the isometric arrows and that alternating isoTNSs mediate entanglement more efficiently than conventional isoTNSs. Second, alternating isoTNSs correspond to a deeper, and thus more representative, sequential-circuit construction of depth 𝒪⁢(𝐿𝑥 ⋅𝐿𝑦) compared to the original isoTNSs of depth 𝒪⁢(𝐿𝑥 +𝐿𝑦). Third, using the Gaussian framework and gradient-based energy minimization, we provide numerical evidence of better bond-dimension scaling and variational energy of alternating isoGfTNSs for ground states of various free-fermionic models, including the Fermi surface, the band insulator, and the 𝑝𝑥 +𝑖⁢𝑝𝑦 mean-field superconductor. Finally, benchmarking on the transverse-field Ising model, we demonstrate that an alternating isoTNS provides substantially improved performance and stability relative to the original isoTNS for the ground-state search algorithm in interacting systems.

Wu, Yantao [Chinese Academy of Sciences, Beijing (↗

Cheaper and more noise-resilient quantum state preparation using eigenvector continuation

Subspace methods are powerful, noise-resilient methods that can effectively prepare ground states on quantum computers. The challenge is to get a subspace with a small condition number that spans the states of interest using minimal quantum resources. In this work, we will use eigenvector continuation to build a subspace from the low-lying states of a set of Hamiltonians. The basis vectors are prepared using truncated versions of standard state preparation methods such as imaginary-time evolution (ITE), adiabatic state preparation (ASP), and variational quantum eigensolver. By using these truncated methods combined with eigenvector continuation, we can directly improve upon them, obtaining more accurate ground-state energies at a reduced cost. We use several spin systems to demonstrate convergence even when methods like ITE and ASP fail, such as ASP in the presence of level crossings and ITE with vanishing energy gaps. We also showcase the noise resilience of this approach beyond the gains already made by having shallower quantum circuits. Furthermore, our findings suggest that eigenvector continuation can be used to improve existing state preparation methods in the near term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gromov ground state in phase space engineering for fusion energy

Phase space engineering by rf waves plays important roles in both thermal D-T fusion and nonthermal advanced fuel fusion, but not all phase space manipulation is allowed; certain fundamental limits exist. In addition to Liouville's theorem, which requires the manipulation to be volume preserving, Gromov's nonsqueezing theorem imposes another constraint. Here, the Gardner ground state is defined as the ground state accessible by smooth volume-preserving maps. However, the extra Gromov constraint should produce a higher-energy ground state. An example of a Gardner ground state forbidden by Gromov's nonsqueezing theorem is given. The challenge question is “What is the Gromov ground state, i.e., the lowest energy state accessible by smooth symplectic maps?” This is a difficult problem. As a simplification, we conjecture that the linear Gromov ground state problem is solvable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE

Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accurate energies for ππ* excited states via exchange scaling: the XS-CASSCF method

The state-averaged complete-active space self-consistent field method (SA-CASSCF) is a widely employed electronic structure method used for studying photochemistry and dynamics owing to its ability to provide a reliable description even of complicated cases while still retaining computational efficiency. However, SA-CASSCF suffers from one Achilles heel, related to the description of ionic ππ* excited states, whose energy is often overestimated by 1–2 eV. In light of this challenge, we present the XS-CASSCF method, a new approach based on the idea of exchange scaling (XS) that screens the involved energy terms to improve the excitation energies of singlet ionic ππ* states. First, we illustrate the power of the XS-CASSCF method using hexatriene and para-quinodimethane as examples, showing that it corrects the targeted ionic states while leaving the other states largely unaffected, giving root-mean-square errors (RMSE) below 0.2 eV for the four lowest states in both cases. Subsequently, XS-CASSCF vertical excitation energies are tested against theoretical best estimates for a set of 11 molecules and 56 excited states. XS-CASSCF performs exceptionally well for the ππ* states of hydrocarbons, reducing the RMSE over 21 excitation energies from 0.96 to 0.27 eV. In the challenging subset of molecules with heteroatoms and a larger number of ππ* and nπ* states, we find that improvements can also be obtained, albeit not as pronounced. We conclude with an outlook into more realistic molecular materials focusing on their singlet–triplet (S 1 /T 1 ) gaps, finding that significant improvements can be obtained along the whole range of S 1 /T 1 gaps studied, going from 0.1 eV to more than 1.5 eV. Owing to notable improvements across significant classes of molecules combined with its conceptual simplicity, we believe that XS-CASSCF is a promising addition to the electronic structure toolbox, serving both as a standalone electronic structure method and as a starting point for further correlated treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

State Machine Modeling of the Space Launch System Solid Rocket Boosters

The Space Launch System is a Shuttle-derived heavy-lift vehicle currently in development to serve as NASA's premiere launch vehicle for space exploration. The Space Launch System is a multistage rocket with two Solid Rocket Boosters and multiple payloads, including the Multi-Purpose Crew Vehicle. Planned Space Launch System destinations include near-Earth asteroids, the Moon, Mars, and Lagrange points. The Space Launch System is a complex system with many subsystems, requiring considerable systems engineering and integration. To this end, state machine analysis offers a method to support engineering and operational e orts, identify and avert undesirable or potentially hazardous system states, and evaluate system requirements. Finite State Machines model a system as a finite number of states, with transitions between states controlled by state-based and event-based logic. State machines are a useful tool for understanding complex system behaviors and evaluating "what-if" scenarios. This work contributes to a state machine model of the Space Launch System developed at NASA Ames Research Center. The Space Launch System Solid Rocket Booster avionics and ignition subsystems are modeled using MATLAB/Stateflow software. This model is integrated into a larger model of Space Launch System avionics used for verification and validation of Space Launch System operating procedures and design requirements. This includes testing both nominal and o -nominal system states and command sequences.

model-based systems engineering↗

Vibrational Quantum-State-Controlled Reactivity in the O 2 + + C 3 H 4 Reaction

Quantum-state-controlled reactivity is a long-standing goal in the field of physical chemistry. In this work, we explore the vibrational-state-dependent behavior of the ion–molecule reaction between O 2 + in distinct vibrational states and two isomers of C 3 H 4 , allene (H 2 C 3 H 2 ) and propyne (H 3 C 3 H). While most products are formed regardless of the vibrational state of O 2 + , the branching ratios are influenced by vibrational excitation, and a new product, C 2 O + , appears exclusively in the excited-state reactions. This selective formation of C 2 O + demonstrates that vibrational excitation can effectively activate a reaction pathway, providing direct evidence of quantum-state control in the reactivity. These results represent an important step toward the goal of quantum-state-controlled chemistry in molecular systems.

Chemical reactions↗

Deciphering the Conflict Between Ion and Electron Percolating Networks in Solid-State Battery Cathodes

High-energy- and power-density solid state batteries require an optimal cathode composition and microstructural arrangement of cathode active material, solid-state electrolyte, conductive carbon, and binder to simultaneously support lithium-ion transport, electron conduction, and storage capacity. The ion and electron conducting phases in solid-state cathodes counteract each other's percolating networks as their mass ratios increase or decrease relative to each other. Here, we investigate targeted mass ratio variations of argyrodite solid electrolyte and two different types of conductive carbon (particles and fibers) in composite LiNi0.8Mn0.1Co0.1O2 (NMC811) solid-state cathodes to ascertain the ionic-electronic tradeoffs in cathode performance. Through ionic and electronic conductivity measurements on composite cathodes, as well as rate-testing and cycling performance in full cells, it is shown that the conductive carbon fibers form a percolative electronic network within the composite at a lower mass ratio (3-5 wt%) than particulate carbon (>5 wt%). The threshold to achieve electronic percolation coincides with higher accessible capacity in the cathode as the active material particles become electronically connected. However, carbon loadings beyond this percolation threshold lead to increased ion transport resistance, arising from disruptions to ionic conduction pathways and degraded contact at the interface between the electrolyte and active materials. Imaging, spectroscopy, and physics-based models quantitatively describe the relationship between carbon and electrolyte compositions and the cell's capacity and rate performance through percolation theory. This work demonstrates the importance of quantitatively understanding percolating networks in solid-state cells and that strategic engineering of conductive carbon morphologies can further increase the energy- and power-density of solid-state cells.

25 ENERGY STORAGE↗

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗