Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deconvolution”

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 253 records · Page 14

Size Effect of Local Current-Voltage Characteristics of MX 2 Nanoflakes: Local Density of States Reconstruction from Scanning Tunneling Microscopy Experiments

Local current-voltage characteristics for low-dimensional transition-metal dichalcogenides (LDTMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments, are of fundamental interest and can be useful for advanced applications. Most of the existing models either have limited applicability for complex-shaped LDTMDs (e.g., those based on the Simmons approach) or require solving of an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on the Tersoff approach). Using a serial expansion of the Tersoff formulas, we propose a flexible method to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We establish a set of key physical parameters, which characterize the tunneling current of a STM-probe–sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different numbers of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around the Fermi level, which allows us to gain insight into the local band structure of LDTMDs. The reconstructed LDOS reveal pronounced size effects for the single-layer, two-layer, and three-layer MoS 2 nanoflakes, which we relate to the low dimensionality and strong bending or corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach, allowing LDOS reconstruction, will be of critical interest for the quantitative description of STM experiments and also of use to better understand the microscopic physical aspects of the surface, strain, and bending contributions to the LDTMDs’ electronic properties.

42 ENGINEERING↗

Mixed-strategy approach to band-edge analysis and modeling in semiconductors

In semiconductor physics, the Urbach rule describes an exponential dependence of the UV-Vis absorption coefficient on photon energy at the band edge. The rule is so ubiquitously applied that a single common mechanism, variability in the local potential, is generally considered to be the cause of exponential band edges. Variability in the local potential may be caused by alloy inhomogeneity, structural defects, impurities, thermal excitation, or any other mechanism that interrupts ideal crystallinity and creates disorder. A crucial point of the Urbach rule is its development for band edges which approach linearity when plotted on a logarithmic scale. The rule was not developed in the context of nonlinear band edges with complex fine structure. Even so, this model has been applied extensively to band edges with such features. Using differential analysis and deconvolution, we develop an analytical technique capable of separating the components responsible for the fine structure observed. This approach reduces challenging signal-artifact convolution problems to a simple linear combination. Furthermore, we develop a model that can be applied to the nonlinear band-edge signals extracted. In this work, we show that in the limiting case of linearity, our model is consistent with and produces the same results as the linear Urbach model. Our model is capable of bypassing “tailing” as an indirect measure of the spread in the local potential and gives a direct measure of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Solid state core-exciton dynamics in NaCl observed by tabletop attosecond four-wave mixing spectroscopy

Nonlinear wave mixing in solids with ultrafast x-rays can provide insight into complex electronic dynamics of materials. Here, tabletop-based attosecond noncollinear four-wave mixing (FWM) spectroscopy using one extreme ultraviolet (XUV) pulse from high harmonic generation and two separately timed few-cycle near-infrared (NIR) pulses characterizes the dynamics of the Na + L 2,3 edge core-excitons in NaCl around 33.5 eV. An inhomogeneous distribution of core-excitons underlying the well-known doublet absorption of the Na+ Γ point core-exciton spectrum is deconvoluted by the resonance-enhanced nonlinear wave mixing spectroscopy. In addition, other dark excitonic states that are coupled to the XUV-allowed levels by the NIR pulses are characterized spectrally and temporally. Approximately <10fs coherence lifetimes of the core-exciton states are observed. Furthermore, the core-excitonic properties are discussed in the context of strong electron-hole exchange interactions, electron-electron correlation, and electron-phonon broadening. This investigation successfully indicates that tabletop attosecond FWM spectroscopies represent a viable technique for time-resolved solid state measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Landau-Ginzburg theory of charge density wave formation accompanying lattice and electronic long-range ordering

Here we propose an analytical Landau-Ginzburg (LG) theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include the electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the LG free energy density as a power expansion with respect to the charge density and other long-range order parameters as well as their spatial gradients and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes, and modulation periods obtained in this paper can be employed to guide the comprehensive physical explanation, deconvolution, and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Determination of β -decay feeding patterns of 88 Rb and 88 Kr using the Modular Total Absorption Spectrometer at ORNL HRIBF

We report precise determination of ground-state feeding in the β decay of fission products is an important but challenging component in modeling reactor antineutrino flux and reactor decay heat. The Modular Total Absorption Spectrometer (MTAS) is a versatile NaI(Tl) detector array that determines the true β-decay pattern free from the pandemonium effect, including precise ground-state feeding intensities. In this paper, we report MTAS results of the β feeding intensities of 88 Rb and 88 Kr, fission products with large cumulative yields in nuclear reactors. By comparing MTAS results with previous measurements, 88 Rb provides a validation of MTAS's ability to determine ground-state feedings in β decays, while the precision of 88 Kr ground-state feeding is improved when compared with the Evaluated Nuclear Structure Data File (ENSDF). The investigation of sources that contribute to β feeding branching uncertainties in MTAS experiments is discussed in detail. Lastly, the deconvolution of 88 Rb decay spectra suggests that MTAS can distinguish an allowed β spectral shape from a first forbidden unique β spectral shape.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

String-based parametrization of nucleon GPDs at any skewness: A comparison to lattice QCD

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) valid at all skewness values. The conformal moments of the GPDs are expressed as sums of the spin-j nucleon A-form factor and the skewness-dependent spin-j nucleon D-form factor. This representation, which fulfills the polynomiality condition (due to Lorentz invariance) and does not rely on model-specific assumptions, is derived from t-channel string exchanges in anti-de Sitter spaces. The spin-j nucleon D-form factor is closely related to the spin-j nucleon A-form factor. We use the Mellin moments from empirical parton distributions to model the spin-j nucleon A-form factors. Using only five Regge slope parameters, fixed from the electromagnetic and gravitational form factors, our string-based parametrization generates accurate singlet, nonsinglet, isovector, and flavor-separated nucleon quark GPDs, along with symmetric nucleon gluon GPDs from their Mellin-Barnes integral representations. Our isovector nucleon quark GPD is in agreement with existing lattice data. Our string-based parametrization should facilitate the empirical extraction and global analysis of nucleon GPDs in exclusive processes, bypassing the deconvolution challenge.

Electron-ion collisions↗

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗

Parametrization of Generalized Parton Distributions from 𝑡-Channel String Exchange in AdS Spaces

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) that is valid for all skewness. Our approach leverages conformal moments, representing them as the sum of spin-𝑗 nucleon 𝐴-form factor and skewness-dependent spin-𝑗 nucleon 𝐷-form factor, derived from 𝑡-channel string exchange in AdS spaces consistent with Lorentz invariance and unitarity. This model-independent framework, satisfying the polynomiality condition due to Lorentz invariance, uses Mellin moments from empirical data to estimate these form factors. With just five Regge slope parameters, our method accurately produces various nucleon quark GPD types and symmetric nucleon gluon GPDs through pertinent Mellin-Barnes integrals. Our isovector nucleon quark GPD is in agreement with existing lattice data, promising to improve the empirical extraction and global analysis of nucleon GPDs in exclusive processes, by avoiding the deconvolution problem at any skewness, for the first time.

QCD phenomenology↗

Three Photon Excited Image Scanning Microscopy for in-Depth Super-Resolution Studies of Biological Samples

Multiphoton laser scanning microscopy is a powerful tool for deep imaging of thick biological samples. Image scanning microscopy (ISM) has demonstrated significant improvements in the signal-to-noise ratio in confocal laser scanning microscopy, while at the same time improving upon the effectively attainable resolution. Two-photon excitation (2PE), combined with ISM, has been shown to allow for deep tissue imaging with enhanced resolution compared to 2PE microscopy. Three-photon excitation (3PE) has enabled record imaging depth and contrast for multiphoton imaging, due to the superior suppression of out-of-focus signal generation. In this paper, we demonstrate super-resolution 3PE ISM. This is achieved using a single-photon avalanche detector array, and 1040-nm pulses for 3PE of blue fluorescence. This method enables subdiffraction limited resolution imaging of biological samples stained with blue fluorescent markers, such as mouse myocardial and spinal cord tissues stained with 4 ′ , 6 -diamidino-2-phenylindole. Deconvolution improves the resolving power further and allows for imaging with better than λ / 8 resolution with respect to the 3PE wavelength λ . With the ISM pixel reassignment procedure, we demonstrate a resolution enhancement of ∼ 1.6 laterally, compared to the resolution attained using a photomultiplier tube in a non-descanned detection arrangement, and a factor of ∼ 1.8 enhancement in axial resolution. The experimentally measured three-dimensional point spread function volume is shrunk ∼ 4.4 -fold, which is close to the theoretically expected enhancement. Published by the American Physical Society 2024

47 OTHER INSTRUMENTATION↗

BioXTAS RAW 2 : new developments for a free open-source program for small-angle scattering data reduction and analysis

BioXTAS RAW is a free open-source program for reduction, analysis and modelling of biological small-angle scattering data. Here, the new developments in RAW version 2 are described. These include improved data reduction using pyFAI ; updated automated Guinier fitting and D max finding algorithms; automated series ( e.g. size-exclusion chromatography coupled small-angle X-ray scattering or SEC-SAXS) buffer- and sample-region finding algorithms; linear and integral baseline correction for series; deconvolution of series data using regularized alternating least squares ( REGALS ); creation of electron-density reconstructions using electron density via solution scattering ( DENSS ); a comparison window showing residuals, ratios and statistical comparisons between profiles; and generation of PDF reports with summary plots and tables for all analysis. Furthermore, there is now a RAW API, which can be used without the graphical user interface (GUI), providing full access to all of the functionality found in the GUI. In addition to these new capabilities, RAW has undergone significant technical updates, such as adding Python 3 compatibility, and has entirely new documentation available both online and in the program.

97 MATHEMATICS AND COMPUTING↗

Surrogate Model for Distribution Networks Influenced by Weather

Here, we propose a method for generating reduced representations of time series and for constructing low dimensional surrogate models for time dependent calculations of power and voltage in distribution networks. We employ Fourier polynomials. The surrogate model strategy is aimed at reducing the computational cost of time dependent simulations, albeit, at the expense of fidelity. The reduced representation is achieved by identifying a small and most consequential subset of degrees of freedom. In power and voltage distribution networks dynamics that are heavily influenced by strong cyclic weather events, e.g., the hourly, diurnal and seasonal cycles, the weather/climate time series spectrum exposes these most energetic components. Once the degrees of freedom are identified their amplitudes are optimized using training data. The key challenge in using spectral methods in power network surrogates is addressing the computation of quotients. For this we propose a numerically-stable deconvolution strategy.

Weather events↗

Binary pseudo-random array standards for calibration of 3D optical surface profilers used for metrology with aspheric x-ray optics

High-accuracy surface metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate x-ray light with nanometer-scale wavelengths. The current and potential capabilities of x˗ray mirror manufacturing are limited by inherent imperfections of the integrated metrology tools. Metrology tools are currently calibrated with super-polished flat test-standard/reference mirrors. This is acceptable for fabrication of slightly curved x-ray optics. However, for even moderately curved aspherical x-ray mirrors the flat-reference calibration is not sufficiently accurate. For micro-stitching interferometry developed for surface measurements with curved x-ray mirrors, the tool aberration errors are known to be transferred into the optical surface topography of x-ray mirrors. Our approach to improving metrology is to thoroughly calibrate the measuring tool and apply the results of the calibration to deconvolution of the measured data. Here we explore the application of a recently developed technique for calibrating the instrument transfer function (ITF) of 3D optical surface profilers to metrology with significantly curved x-ray optics. The technique, based on test standards patterned with two-dimensional (2D) binary pseudo-random arrays (BPRAs), employs the unique properties of the BPRA patterns in the spatial frequency domain. The inherent 2D power spectral density of the pattern has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range and field of view of an instrument. The high efficacy of the technique has been previously demonstrated in application to metrology with flat and slightly curved optics. Here, we concentrate on development of an efficient fabrication process for production of highly randomized (HR) BPRA test standards on flat and 500-mm spherical optical substrates. We also compare and discuss the results of the ITF calibration of an interferometric microscope when using the HR BPRA standards on flat and curved substrates.

Munechika, Keiko↗

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deciphering competing elementary steps to correlate electrocatalyst chemical state with activity

The overpotential in multielectron transfer heterogeneous electrocatalysis fundamentally arises from thermodynamic and kinetic disparities among elementary steps; however, deciphering coupled and competing steps has long remained a challenge. Here, we establish an electrochemical deconvolution paradigm based on key processes in electrocatalytic reactions, such as charge accumulation, electron/proton transfer, and intermediate evolution, to resolve competing elementary steps. Taking the oxygen evolution reaction as a prototypical reaction, we design a model catalyst featuring a precise isolated cation-anion vacancy pair and track the previously elusive electrochemical behavior of lattice oxygen by disentangling interference from adsorbed oxygen intermediates. Mechanistically, the lattice oxygen oxidation pathway originates from the spontaneous, nonelectrochemical deprotonation of replenished water molecules coordinated to unsaturated cation sites. Alternating current techniques further reveal that although lattice oxygen oxidation requires a higher potential than metal oxidation, it exhibits faster kinetics, providing insight into its superior catalytic activity. These findings establish a direct experimental correlation between the initial chemical state and the catalytic activity and prove surface-confined lattice oxygen cycling. Furthermore, expanding conventional potential-current analysis into a multidimensional framework enables disentanglement of thermodynamic and kinetic contributions of key elementary steps, thereby guiding the rational optimization of various complex multielectron transfer reactions.

OER↗

Structure-based decoupling of the pro- and anti-inflammatory functions of interleukin-10

Interleukin-10 (IL-10) is an immunoregulatory cytokine with both anti-inflammatory and immunostimulatory properties and is frequently dysregulated in disease. We used a structure-based approach to deconvolute IL-10 pleiotropy by determining the structure of the IL-10 receptor (IL-10R) complex by cryo–electron microscopy at a resolution of 3.5 angstroms. The hexameric structure shows how IL-10 and IL-10Rα form a composite surface to engage the shared signaling receptor IL-10Rβ, enabling the design of partial agonists. IL-10 variants with a range of IL-10Rβ binding strengths uncovered substantial differences in response thresholds across immune cell populations, providing a means of manipulating IL-10 cell type selectivity. Some variants displayed myeloid-biased activity by suppressing macrophage activation without stimulating inflammatory CD8 + T cells, thereby uncoupling the major opposing functions of IL-10. These results provide a mechanistic blueprint for tuning the pleiotropic actions of IL-10.

60 APPLIED LIFE SCIENCES↗

A Variational Inference Approach to Inverse Problems with Gamma Hyperpriors

Hierarchical models with gamma hyperpriors provide a flexible, sparse-promoting framework to bridge L 1 and L 2 regularizations in Bayesian formulations to inverse problems. Despite the Bayesian motivation for these models, existing methodologies are limited to maximum a posteriori estimation. The potential to perform uncertainty quantification has not yet been realized. This paper introduces a variational iterative alternating scheme for hierarchical inverse problems with gamma hyperpriors. The proposed variational inference approach yields accurate reconstruction, provides meaningful uncertainty quantification, and is easy to implement. In addition, it lends itself naturally to conduct model selection for the choice of hyperparameters. Here, we illustrate the performance of our methodology in several computed examples, including a deconvolution problem and sparse identification of dynamical systems from time series data.

Bayesian shrinkage↗

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗