Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical analytics”

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

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

At least 1,135 records · Page 63

Combined Neyman–Pearson chi-square: An improved approximation to the Poisson-likelihood chi-square

We describe an approximation to the widely-used Poisson-likelihood chi-square using a linear combination of Neyman’s and Pearson’s chi-squares, namely “combined Neyman–Pearson chi-square” ($X^2_{CNP}$). Through analytical derivations and toy model simulations, we show that leads to a significantly smaller bias on the best-fit model parameters compared to those using either Neyman’s or Pearson’s chi-square. When the computational cost of using the Poisson-likelihood chi-square is high, $X^2_{CNP}$ provides a good alternative given its natural connection to the covariance matrix formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.

Buchert, T.↗

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytic results for scalar-mediated Higgs boson production in association with two jets

Abstract We present compact analytic formulae for all one-loop amplitudes representing the production of a Higgs boson in association with two jets, mediated by a colour triplet scalar particle. Many of the integral coefficients present for scalar mediators are identical to the case when a massive fermion circulates in the loop, reflecting a close relationship between the two theories. The calculation is used to study Higgs boson production in association with two jets in a simplified supersymmetry (SUSY) scenario in which the dominant additional contributions arise from loops of top squarks. The results presented here facilitate an indirect search for top squarks in this channel, by a precision measurement of the corresponding cross section. However, we find that the potential for improved discrimination between the Standard Model and SUSY cases suggested by the pattern of results in the one- and two-jet samples is unlikely to be realized due to the loss in statistical power compared to an inclusive analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vhf/uhf stellar calibration error analysis.

Results of a statistical error analysis performed to determine the degree of uncertainty encountered when calibrating steerable vhf/uhf receiving antennas with the stellar calibration method. The analysis considers the propagation of precision error indices. It is shown that an antenna gain calibration by the stellar method has a one-sigma (1 sigma) uncertainty of plus or minus 0.65 dB at 1440 MHz (L-band), and plus or minus 0.8 dB (1 sigma) at 136 MHz (vhf). Somewhat increased precision can be achieved by monitoring the antenna gain-to-noise temperature ratio at a station; a worst-case uncertainty of plus or minus 0.4 dB (1 sigma) can be realized at both L-band and vhf. Finally, field test measurements of antenna gain, obtained at 136 MHz in the NASA space tracking and data acquisition network, demonstrate an uncertainty of plus or minus 1.0 dB or less, which effectively confirms the analytical result.

Taylor, R. E.↗

Three-dimensional imaging of grain boundaries via quantitative fluorescence X-ray tomography analysis

Abstract Three-dimensional visualization of material composition within multiple grains and across complex networks of grain boundaries at nanoscales can provide new insight into the structure evolution and emerging functional properties of the material for diverse applications. Here, using nanoscale scanning X-ray fluorescence tomography, coupled with an advanced self-absorption correction algorithm developed in this work, we analyze the three-dimensional gain distributions and compositions in a Ce 0.8 Gd 0.2 O 2-δ -CoFe 2 O 4 mixed ionic-electronic conductor system with high accuracy and statistical significance. Our systematic investigation reveals an additional emergent phase and uncovers highly intriguing composition stability ranges for the multiple material phases within this system. The presented visualization of composition variations across complex interfaces, supported by our quantitative composition analysis, discloses mechanistic pathways of the diverse phase transformations occurring in the material synthesis, providing insights for the optimization of transport properties in the mixed ionic-electronic conductor system.

36 MATERIALS SCIENCE↗

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum geometry embedded in unitarity of evolution: Revealing its impacts as geometric oscillation and dephasing in spin resonance and crystal bands

Quantum Hall effects provide intuitive ways of revealing the topology in crystals, i.e., each quantized “step” represents a distinct topological state. Here, we seek a counterpart for “visualizing” quantum geometry, which is a broader concept. Here we show how geometry emerges in quantum as an intrinsic consequence of unitary evolution, composing a framework compatible with quantum metric and independent of specific details or approximations, suggesting quantum geometry may have widespread applicability. Indeed, we exemplify geometric observables, such as oscillation, dephasing, in magnetic resonance or band driving scenarios. Anomalies, supported by both analytic and numerical solutions, underscore the advantages of adopting a geometric perspective, potentially yielding distinguishable experimental signatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Far-Ultraviolet Number Counts of Field Galaxies

The Number counts of far-ultraviolet (FUV) galaxies as a function of magnitude provide a direct statistical measure of the density and evolution of star-forming galaxies. We report on the results of measurements of the rest-frame FUV number counts computed from data of several fields including the Hubble Ultra Deep Field, the Hubble Deep Field North, and the GOODS-North and -South fields. These data were obtained from the Hubble Space Telescope Solar Blind Channel of the Advance Camera for Surveys. The number counts cover an AB magnitude range from 20-29 magnitudes, covering a total area of 15.9 arcmin'. We show that the number counts are lower than those in previous studies using smaller areas. The differences in the counts are likely the result of cosmic variance; our new data cover more area and more lines of sight than the previous studies. The slope of our number counts connects well with local FUV counts and they show good agreement with recent semi-analytical models based on dark matter "merger trees".

Voyer, Elysse N.↗

The Next Generation X-Ray Galaxy Survey with eROSITA

eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e. not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50–200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect ≳15 000 galaxies (3σ significance) at 50–200 Mpc, which is ~ 100× more normal galaxies than detected in any X-ray survey to date.

Surveys↗

A Model for Atomic Precision p-Type Doping with Diborane on Si(100)-2×1

Diborane (B 2 H 6 ) is a promising molecular precursor for atomic precision p-type doping of silicon that has recently been experimentally demonstrated. In this study, we use density functional theory (DFT) calculations to determine the reaction pathway for diborane dissociating into a species that will incorporate as electrically active substitutional boron after adsorbing onto the Si(100)-2×1 surface. Our calculations indicate that diborane must overcome an energy barrier to adsorb, explaining the experimentally observed low sticking coefficient (<1 × 10 –4 at room temperature) and suggesting that heating can be used to increase the adsorption rate. Upon sticking, diborane has an ≈50% chance of splitting into two BH 3 fragments versus merely losing hydrogen to form a dimer such as B 2 H 4 . As boron dimers are likely electrically inactive, whether this latter reaction occurs is shown to be predictive of the incorporation rate. The dissociation process proceeds with significant energy barriers, necessitating the use of high temperatures for incorporation. Using the barriers calculated from DFT, we parameterize a Kinetic Monte Carlo model that predicts the incorporation statistics of boron as a function of the initial depassivation geometry, dose, and anneal temperature. Our results suggest that the dimer nature of diborane inherently limits its doping density as an acceptor precursor and furthermore that heating the boron dimers to split before exposure to silicon can lead to poor selectivity on hydrogen and halogen resists. This suggests that, while diborane works as an atomic precision acceptor precursor, other non-dimerized acceptor precursors may lead to higher incorporation rates at lower temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deacetylation and Mechanical Refining (DMR) and Deacetylation and Dilute Acid (DDA) Pretreatment of Corn Stover, Switchgrass, and a 50/50 Corn Stover/Switchgrass Blend

The most economical feedstock available to future lignocellulosic biorefineries may be a feedstock blend. Corn stover, switchgrass, and a 50:50 blend of corn stover and switchgrass were subjected to three different dilute alkali deacetylation conditions in a 90 L reactor prior to either mechanical refining (DMR) in a disc refiner followed by a Szego mill or dilute acid pretreatment (DDA) at five different conditions in a 4 L batch steam explosion reactor. The forty-five DDA and nine DMR slurries were subjected to high solids enzymatic hydrolysis at an enzyme loading of 15 mg protein per gram cellulose using Novozymes Cellic CTec3 (DDA) or a 4:1 mix of Novozymes Cellic CTec3:HTec 3 (DMR). The highest glucose yield was achieved using DMR processed corn stover from the high severity deacetylation condition (92 °C and 100 kg NaOH/ODMT) at 95.9% ± 0.8%, while the greatest monomeric glucose yields for switchgrass were between 75% to 80%. The glucose yields from saccharification of the 50:50 corn stover/switchgrass blend was approximately the average of the two individual feedstocks for both DDA and DMR processes. Statistical analysis showed that for a given feedstock type dilute alkali deacetylation severity had a significant impact (p < 0.05) on enzymatic hydrolysis yields, and for a given deacetylation severity the feedstock type had a significant effect on glucose conversion with the addition of switchgrass resulting in lower glucose yields. A blending strategy to maximize enzymatic hydrolysis sugar yields while reducing carbohydrate solubilization during the alkaline treatment step was proposed that allows deacetylation of feedstocks with different recalcitrance at optimal conditions for that feedstock, then blending prior to enzymatic hydrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Information Scrambling and Loschmidt Echo

We demonstrate analytically and verify numerically that the out-of-time order correlator is given by the thermal average of Loschmidt echo signals. This provides a direct link between the out-of-time-order correlator—a recently suggested measure of information scrambling in quantum chaotic systems—and the Loschmidt echo, a well-appreciated familiar diagnostic that captures the dynamical aspect of chaotic behavior in the time domain, and is accessible to experimental studies.

0-dimensional systems↗

Preliminary Design of Critical Experiments Involving Commercially Available B 4 C Neutron Absorber Plates with Low-Enriched UO 2 Fuel

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) is an initiative to provide high-quality benchmark data in a standardized format for criticality safety analysts to use to validate calculational tools and nuclear data. In the last released 2022 version, the ICSBEP handbook contains over 5,000 critical experiment descriptions, results, and associated models. This paper describes the critical configurations proposed and calculations performed as part of the conceptual design phase of Integral Experiment Request (IER) 554 (IER-554). The goal of IER-554 is to add to the benchmark experiments available to the criticality safety community by designing critical experiments that can analyze how adding commercially available Boralcan plates to an assembly of low-enriched UO2 fuel rods affects the effective neutron multiplication factor (k eff ). Boralcan, a neutron absorber product developed by Rio Tinto, is made of B4C in an Al 1100 matrix and is commonly used for criticality suppression in fuel storage pools. In the current version of the ICSBEP handbook, only a few dated critical experiments involve B4C materials, and none involve Boralcan in a thin plate shape. The experiments designed as part of IER-554 are intended to be performed at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. SPRF/CX is a well-characterized assembly considered trustworthy by the benchmarking community because of the numerous high-quality evaluations with very low experimental uncertainties included in the ICSBEP handbook. Figure 1 illustrates the assembly for a particular configuration from LCT-078, one of the published benchmarks in the ICSBEP handbook. The experiments will be moderated and reflected by light water at ambient atmospheric pressure. Before the plates are inserted in the critical assembly, they will be characterized with x-ray computed tomography (XCT) to identify the sizes and distribution of the B 4 C powder particles inside the plates. The preliminary design calculations were performed with SCALE 6.3.0 using the KENO V.a sequence for the criticality calculations, the CE-TSUNAMI-3D sequence for the sensitivity and uncertainty studies, and the ENDF/B-VII.1 Continuous Energy cross section library. All the calculation results presented have 10 pcm statistical uncertainties. A set of critical experiments with Boralcan plates could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material. Experiments would also help validate k eff calculations, and the industry could use these validations to change the B loading credit limits from the US Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. This paper summarizes only part of the analysis documented in the preliminary design report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Comparison of modular analytical wake models to the Lillgrund wind plant

Efficient and accurate wake models are required for wind plant performance modeling and the suite of engineering processes that support wind plant layout, control, and monitoring. Although many analytical and engineering wake models with low computational costs have been proposed, their ability to represent the power production of large wind plants in a wide range of atmospheric conditions is not completely understood. The following validation study reviews the underlying theory for analytical wake models, outlines quality control procedures for observational data, and compares model results with observational data from the Lillgrund Wind Plant. Lillgrund makes a valuable case study for wake modeling because of its regular arrangement and the relatively close spacing of constituent wind turbines, which lead to regular and significant wake interactions within the wind plant and the development of deep array flow conditions. Formulations for the velocity deficit, wake-added turbulence, and wake superposition methods are considered in a modular sense, yielding many possible configurations to represent wind turbine wakes, of which seven are examined in detail. Velocity deficit models that account for flow conditions in the near wake are better able to reproduce power production for wind turbines in the transitional region of the wind plant, where wind turbines experience as many as five wakes from upstream turbines. Finally, in the deep array, where power production reaches asymptotic values and wake statistics become quasi-periodic, wake superposition schemes become the largest driver in error reduction; using the linear or maximum wake superposition methods can reduce the relative root mean square error by as much as 40% in the deep array.

17 WIND ENERGY↗

CODEX weak lensing mass catalogue and implications on the mass–richness relation

The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at 0.35 < z < 0.65. It was selected from ROSAT data in the 10 000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX programme on cluster mass measurement, including a reanalysis of CFHTLS Wide data, with 25 individual lensing-constrained cluster masses. We employ lensfit shape measurement and perform a conservative colour–space selection and weighting of background galaxies. Using the combination of shape noise and an analytic covariance for intrinsic variations of cluster profiles at fixed mass due to large-scale structure, miscentring, and variations in concentration and ellipticity, we determine the likelihood of the observed shear signal as a function of true mass for each cluster. We combine 25 individual cluster mass likelihoods in a Bayesian hierarchical scheme with the inclusion of optical and X-ray selection functions to derive constraints on the slope α, normalization β, and scatter σln λ|μ of our richness–mass scaling relation model in log-space: ${\langle {\rm In}\,\, \lambda\!\!\mid\!\!\mu\rangle = \alpha\mu + \beta,} $ with μ = ln (M 200c /M piv ), and M piv = 10 14.81 M ⊙ . We find a slope $\alpha = 0.49^{+0.20}_{-0.15}$, normalization $\exp (\beta) = 84.0^{+9.2}_{-14.8}$, and $\sigma _{\ln \lambda | \mu } = 0.17^{+0.13}_{-0.09}$ using CFHT richness estimates. In comparison to other weak lensing richness–mass relations, we find the normalization of the richness statistically agreeing with the normalization of other scaling relations from a broad redshift range (0.0 < z < 0.65) and with different cluster selection (X-ray, Sunyaev–Zeldovich, and optical).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mineralogy of the Solar System

The coming decade will witnesses the first sample return missions from solar system bodies in 30 years. These samples will all be very small, some missions return only a few milligrams of total mass. Fortunately, the capability of modem methods to characterize ultra-small samples is well established from analysis of interplanetary dust particles (IDPs), interstellar grains recovered from meteorites, and other materials requiring ultra-sensitive analytical capabilities. Powerful analytical techniques are available that require, under favorable circumstances, single particles of only a few nanograms for entire suites of fairly comprehensive characterizations. A returned sample of greater than 1,000 particles with total mass of just one microgram permits comprehensive quantitative geochemical measurements that are impractical to can-y out in situ by flight instruments. With the Galileo flybys of Gaspra and Ida, it is now recognized that even very small airless bodies have indeed developed a particulate regolith. Acquiring a sample of the bulk regolith, a simple sampling strategy, provides two critical pieces of information about the body. Regolith samples are excellent bulk samples since they normally contain all the key components of the local environment, albeit in particulate form. Furthermore, since this fine fraction dominates remote measurements, regolith samples also provide information about surface alteration processes and are a key link to remote sensing of other bodies. Studies indicate that a statistically significant number of nanogram-sized particles should be able to characterize the regolith of a primitive asteroid, although the presence of larger components within even primitive meteorites (e.g.. Murchison), e.g. chondrules, CAI, large crystal fragments, etc., points out the limitations of using data obtained from nanogram-sized samples to characterize entire primitive asteroids. However, most important asteroidal geological processes have left their mark on the matrix, since this is the finest-grained portion and therefore most sensitive to chemical and physical changes. Thus, the following information can be learned from this fine grain size fraction alone: (1) mineral paragenesis; (2) regolith processes, (3) bulk composition; (4) conditions of thermal and aqueous alteration (if any); (5) relationships to planets, comets, meteorites (via isotopic analyses, including oxygen; (6) abundance of water and hydrated material; (7) abundance of organics; (8) history of volatile mobility, (9) presence and origin of presolar and/or interstellar material.

Zolensky, Michael E.↗