SEARCH · Engineering Papers

Results for “shape measurement”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Using 3D Measurements from MRI to Improve VVUQ

Magnetic resonance velocimetry (MRV) is a diagnostic that can measure 3D, three-component turbulent velocity fields in arbitrarily complex flow configurations. The approach uses magnetic resonance imaging (MRI), which is a system commonly available in radiology departments or medical research centers. MRV acquires signal from hydrogen protons in water channel flows. Despite the obvious utility of the measurements for investigating a variety of flows, comparisons with computational fluid dynamic (CFD) simulations largely retain traditional metrics, comparing velocity profiles and planes at discrete regions within the flow rather than using the 3D nature of the measured field. In this effort, MRV was conducted in a simple water channel with six spanwise-centered periodic obstacles roughly shaped as cubes. At a fully turbulent Reynolds number of 15,000, the channel includes two flows: a main streamwise flow and, in the wake of the second obstacle, an injected flow oriented perpendicular to the streamwise flow. The flow geometry includes partial obstacles on the side walls and is inherently 3D because of the interaction of the two streams and the wake features, making the flow challenging to completely measure with traditional optical techniques. Aside from line profiles and planar comparisons, two 3D metrics are used to compare the data with a steady Reynolds-averaged Navier–Stokes (RANS) simulation result. A brief discussion about the comparison is provided, including comments about uncertainty.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗

DOE OSTI2025

Topology-Dependent Polymer Stretching and Scission in Solution at Extreme Shear Rates

There has been significant interest in the engineering of polymer topology to control rheology and mechanical stability in dilute solutions for applications involving extreme shear rate flows. However, methods to experimentally probe properties at relevant shear rates (≈ 10 4 –10 6 s –1 ) have remained ex situ, obscuring access to measures of polymer deformation and rheology that would otherwise provide mechanistic insight into the topology-dependent properties that control their behavior in extreme shear flows. In this study, we used novel in situ small angle neutron scattering measurements in a capillary rheometer (capillary rheo-SANS) to simultaneously measure solution viscosities and polymer deformations in high shear on a series of chemically homologous topology-defined polymers including linear, randomly branched, and star-shaped molecules. We demonstrate that differences in the onset of chain stretching and shear thinning of these polymers in dilute solution are controlled primarily by differences in their molecular relaxation time. These differences correlate with differences in chain scission inferred from ex situ measurements at more extreme shear rates. Together, the results demonstrate a direct coupling between chain deformation and scission, and suggest that the dominant effect of branching as a means to impart resilience against mechanical degradation is through differences in relaxation dynamics due to branching. We anticipate that these results will provide key insights to engineer topology-controlled polymers for rheological modification, mechanical stability, and controlled mechano-chemistry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Measurement and modeling of beam transport in the FODO line of the Spallation Neutron Source Beam Test Facility

Ongoing studies at the Spallation Neutron Source (SNS) Beam Test Facility (BTF) seek to understand and model bunch dynamics in a high-power LINAC front-end. The BTF has recently been upgraded with a reconfiguration from a U-shaped line to a Straight line. We report the current state of model benchmarking, with a focus on RMS beam sizes within the FODO line. The beam measurement is obtained via three camera/screen pairs in the FODO line. This presentation discusses the methodology and results of this measurement.

Thompson, Trent↗

DOE OSTI2023

Λ Baryon Production in ν¯µ Interactions in the MicroBooNE Detector

The Cabibbo suppressed production of $\Lambda$ baryons in anti-neutrino interactions with nuclei is a rare process that is yet to be measured with a modern neutrino detector with automated reconstruction. The cross section for this process is sensitive to a number of unique nuclear effects, most notably the secondary interactions of the produced hyperon while attempting to escape from the nucleus. Other interactions within the nuclear remnant can impact the estimation of neutrino energy in oscillation measurements, and thus an accurate description of the nuclear environment is required. The strangeness violating hyperon production process is only available to anti-neutrinos. The model of this interaction is implemented into the NuWro neutrino interaction Monte Carlo simulation, and some predictions are presented, focusing on the role of nuclear effects. This model introduces a hyperon-nucleus potential, which calculations from hypernuclear theory permit to be strongly repulsive in th e case of $\Sigma$ baryons. The presence of this potential is found to sculpt the shape of the differential cross section in some variables. The MicroBooNE detector will be described, followed by a description of a measurement of the flux averaged, restricted phase space cross section of Cabibbo suppressed $\Lambda$ baryon production. A sophisticated event selection is employed, as a very large quantity of background neutrino interactions must be removed to perform the measurement with any sensitivity. This selection introduces some novel techniques such as the island finding method, and achieves a background reduction of $\sim 10^6$, with an efficiency of around 7\%. The calculation of the systematic uncertainties will be explained, including two procedures explored to handle sources of background with extremely poor simulation statistics: an in-situ constraint using data from sidebands, and a visual inspection of the data and simulation to remove the troublesome background events. The sensitivity to the $\Lambda$ baryon production cross section is calculated in the form of Bayesian posterior probability distributions, combining the systematic uncertainties with data and simulation statistical uncertainties. As a rare process, the statistical uncertainties are highly non-Gaussian, and the Bayesian approach is applied to include the full shapes of these uncertainties. Data corresponding to $2.2 \times 10^{20}$ protons on target of neutrino mode running and $4.9 \times 10^{20}$ protons on target of anti-neutrino running is analysed. When the data was unblinded, five $\Lambda$ production candidates were selected from the data, consistent with the MC simulation prediction of $5.3 \pm 1.1$ events. The final estimated cross section is $1.8^{+2.0}_{-1.6} \times 10^{-40}$cm$^2/$Ar when employing the sideband constraint procedure. A similar result of $2.0^{+2.2}_{-1.8} \times 10^{-40}$cm$^2/$Ar is obtained when performing the visual scan instead. The methods used in t his analysis are intended to be easily exported to other LArTPC detectors such as the Short Baseline Near Detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Time-resolved Rotational Velocities in the Upper Atmosphere of WASP-33 b

While steady empirical progress has been made in understanding the structure and composition of hot-planet atmospheres, direct measurements of velocity signatures, including winds, rotation, and jets, have lagged behind. Quantifying atmospheric dynamics of hot planets is critical for a complete understanding of their atmospheres, and such measurements may even illuminate other planetary properties, such as magnetic field strengths. In this manuscript we present the first detection of the Balmer lines Hα and Hβ in the atmosphere of the ultra-hot Jupiter WASP-33 b. Using atmospheric models that include the effects of atmospheric dynamics, we show that the shape of the average Balmer line transmission spectrum is consistent with rotational velocities in the planet’s thermosphere of v{sub rot}=10.1{sub −1.0}{sup +0.8} km s{sup −1}. We also measure a low-significance day- to nightside velocity shift of −4.6{sub −3.4}{sup +3.4} km s{sup −1} in the transmission spectrum, which is naturally explained by a global wind across the planet’s terminator. In a separate analysis the time-resolved velocity centroids of individual transmission spectra show unambiguous evidence of rotation, with a best-fit velocity of 10.0{sub −2.0}{sup +2.4} km s{sup −1}, consistent with the value of v {sub rot} derived from the shape of the average Balmer line transmission spectrum. Our observations and analysis confirm the power of time-resolved transmission spectra with a high signal-to-noise ratio to measure the velocity structures in exoplanet atmospheres. The high rotational and wind velocities we measure highlight the need for more detailed 3D global climate simulations of the rarefied upper atmospheres of ultra-hot gas giants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Micro-structuration effects on local magneto-transport in [Co/Pd]IrMn thin films

We measured the local magneto-transport (MT) signal with an out-of-plane magnetic field, including magneto-resistance (MR) and Extraordinary Hall effect (EHE), in exchange-biased [Co/Pd]IrMn thin multilayers that are micro-structured with a 100 μm window. We found that when measured locally around the window, the MT signal deviate from the expected behavior. We studied possible causes, including film micro-structuration, electrical contact geometry as well as magnetic field angular tilt. We found that tilting the magnetic field direction with respect to the normal direction does not significantly affect the MT signal, whereas the positioning and geometry of the contacts seem to highly affect the MT signal. For comparison purposes, we carried these MT measurements using the Van-der-Pauw method on a set of four microscopic contacts directly surrounding the window, and on another set of micro-contacts located outside the window, as well as a set of four contacts positioned several millimeters away of each other at the corners of the wafer. If the contacts are sufficiently far apart, the EHE and MR signals have the expected shape and are not significantly affected by the presence of the window. If, on the other hand, the contacts are micro-positioned, the shape of the EHE signal is drastically deformed, and may be modeled as a mix of the standard EHE and MR signals measured on the outer contacts. Furthermore, if the micro-contacts are located directly around the window, the deformation is amplified, and the weight of the MR signal in the mix is further increased by about 40 %. This suggests that the electron path in the Hall geometry is disturbed by both the proximity of the electrodes and by the presence of the window, which both contribute to the deformation for about two-third and one third, respectively.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Shaping dark photon spectral distortions

The cosmic microwave background (CMB) spectrum is an extraordinary tool for exploring physics beyond the Standard Model. The exquisite precision of its measurement makes it particularly sensitive to small effects caused by hidden sector interactions. In particular, CMB spectral distortions can unveil the existence of dark photons which are kinetically coupled to the standard photon. In this work, we use the COBE-FIRAS dataset to derive accurate and robust limits on photon-to-dark-photon oscillations for a large range of dark photon masses, from 10$^{−10}$ to 10$^{−4}$ eV. We consider in detail the redshift dependence of the bounds, computing CMB distortions due to photon injection/removal using a Green’s function method. Our treatment improves on previous results, which had set limits studying energy injection/removal into baryons rather than photon injection/removal, or ignoring the redshift evolution of distortions. The difference between our treatment and previous ones is particularly noticeable in the predicted spectral shape of the distortions, a smoking gun signature for photon-to-dark-photon oscillations. The characterization of the spectral shape is crucial for future CMB missions, which could improve the present sensitivity by orders of magnitude, exploring regions of the dark photon parameter space that are otherwise difficult to access [inline-graphic not available: see fulltext].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Accelerating Large-Format Metal Additive Manufacturing: How Controls R&D is Driving Speed, Scale, and Efficiency

This article highlights work at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility to develop closed-loop, feedback control for laser-wire based Directed Energy Deposition, a form of metal Big Area Additive Manufacturing (m-BAAM), a process being developed in partnership with GKN Aerospace specifically for the production of Ti-6Al-4V pre-forms for aerospace components. A large-scale structural demonstrator component is presented as a case-study in which not just control, but the entire 3D printing workflow for m-BAAM is discussed in detail, including design principles for large-format metal AM, toolpath generation, parameter development, process control, and system operation, as well as post-print net-shape geometric analysis and finish machining. In terms of control, a multi-sensor approach has been utilized to measure both layer height and melt pool size, and multiple modes of closed-loop control have been developed to manipulate process parameters (laser power, print speed, deposition rate) to control these variables. Layer height control and melt pool size control have yielded excellent local (intralayer) and global (component-level) geometry control, and the impact of melt pool size control in particular on thermal gradients and material properties is the subject of continuing research. Further, these modes of control have allowed the process to advance to higher deposition rates (exceeding 7.5 lb/hr), larger parts (1-meter scale), shorter build times, and higher overall efficiency. The control modes are examined individually, highlighting their development, demonstration, and lessons learned, and it is shown how they operate concurrently to enable the printing of a large-scale, near net shape Ti-6Al-4V component.

Gibson, Brian↗

DOE OSTI2022

Joint Determination of Reactor Antineutrino Spectra from U 235 and Pu 239 Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of 235U and 239Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of 235U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant 235U and 239Pu isotopes and improves the uncertainty of the 235U spectral shape to about 3%. The 235U and 239Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the 235U and 239Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Weak Gravitational Lensing of Low Surface Brightness Galaxies in the DES Year 3 Catalog

We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at $g-i\ge 0.60$ and $g-i\le 0.60$, respectively. We use the Y3 \metacal\ shape catalog as the source sample, with a number density of {\color{red}$5.59 \mathrm{gal}/\mathrm{arcmin}^2$}. We measure the tangential shear around the lens galaxies across angular scales of $0.25-400$ arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. W e fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at $7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}$. We place a 95\% upper bound on the subhalo mass at $2.8\times 10^{11} M_\mathrm{\odot}$. We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR \citep{Moster2010}. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Tungsten erosion and divertor leakage from the DIII-D SAS-VW tungsten-coated divertor in experiments with neon gas seeding

Collector probes have been used to examine tungsten divertor leakage in a variety of scenarios with low-Z impurity seeding during operation with the new tungsten-coated SAS-VW divertor in DIII-D. Measurements of tungsten deposition on collector probes inserted into the far Scrape-off-Layer (SOL) are used to deduce how efficiently tungsten leaks out of the closed, V-shaped divertor after it is eroded from the target surfaces. Qualitative differences in the tungsten deposition patterns across the collector probes provide clear experimental evidence that the SOL conditions depend on the low-Z impurity seeding conditions. These measurements show that in scenarios where neon gas is injected into the plasma, the tungsten divertor leakage and SOL transport depend on the poloidal location from which the neon is injected. In particular, neon injection from the Inner Midplane and Outer Midplane appear to each result in higher divertor leakage by a factor of 2 to 3 compared to cases with neon injection from either the SOL Crown or from the SAS-VW divertor itself.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Mass and density of individual frozen hydrometeors

A new precipitation sensor, the Differential Emissivity Imaging Disdrometer (DEID), is used to provide the first continuous measurements of the mass, diameter, and density of individual hydrometeors. The DEID consists of an infrared camera pointed at a heated aluminum plate. It exploits the contrasting thermal emissivity of water and metal to determine individual particle mass by assuming that energy is conserved during the transfer of heat from the plate to the particle during evaporation. Particle density is determined from a combination of particle mass and morphology. A Multi-Angle Snowflake Camera (MASC) was deployed alongside the DEID to provide refined imagery of particle size and shape. Broad consistency is found between derived mass–diameter and density–diameter relationships and those obtained in prior studies. However, DEID measurements show a generally weaker dependence with size for hydrometeor density and a stronger dependence for aggregate snowflake mass.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

The effects of low- Z shielding on uranium isotope discrimination using the time-emission profiles of long-lived delayed neutrons

Characteristic radiation signals can often be obscured or distorted by relatively small amounts of shielding, as is the case for highly-enriched uranium (HEU). Methods that can detect and identify special nuclear material (SNM) even in the presence of shielding are therefore of high value in many nuclear security and nonproliferation applications. It has been previously shown that the buildup and decay time profiles of long-lived delayed neutron emission from induced fission are unique and can be used to successfully discriminate between uranium isotopes and infer enrichment without the need for a calibration standard. Because the long-lived delayed neutron precursors exhibit decay times on the order of tens of seconds, their time profiles are not sensitive to the delays associated with scattering or diffusion in shielding materials, which occur on much shorter time scales. The use of capture-based neutron detectors can also mitigate the effects of shielding, because unlike recoil-based detectors, they do not have a low-energy detection threshold, so that neutrons that lose their energy in the shielding are still detectable. In this work, the effects of low-Z shielding on the discrimination of uranium isotopes using delayed neutron time profile measurements are investigated for the first time. Furthermore, the delayed neutron buildup and decay time profiles of polyethylene-shielded HEU and depleted uranium objects are measured using two different capture-gated detector designs. The recorded profiles show good agreement with the shapes predicted from six-group delayed neutron models even with the shielding in place, and the enrichment of the HEU and DU objects was estimated accurately to within 5%. While this method shows potential for discriminating and characterizing extreme cases of enrichment when a moderate amount of shielding is present, simulation results suggest that significant limitations exist for intermediate enrichment levels due to distortions to the delayed neutron time profile caused by neutron reflection in the shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Interplay of Quantum Size Effect and Tensile Strain on Surface Morphology of β-Sn(100) Islands

The quantum size effect (QSE) and strain effect are two key factors influencing the surface morphology of thin films, which can increase film surface roughness through QSE-induced thickness oscillation and strain-induced island formation, respectively. Surface roughness usually manifests in the early stages of film growth and diminishes beyond a critical thickness. In this work, we employ molecular beam epitaxy (MBE) to grow β-Sn(100) islands with varying thickness N on bilayer graphene-terminated 6H-SiC(0001) substrates. Scanning tunneling microscopy and spectroscopy measurements reveal an inverse surface roughness effect that highlights the interplay of QSE and misfit strain in shaping the surface morphology of β-Sn(100) islands. For N ≤ 10, the islands exhibit flat surfaces, while for N ≥ 26, the island surfaces become corrugated and patterned. For the intermediate range, i.e., 12 ≤ N ≤ 24, both flat and patterned surfaces coexist, with the percentage coverage of the patterned surface oscillating as a function of N. By performing density functional theory calculations, we demonstrate that the unusual surface pattern evolution in our MBE-grown β-Sn(100) islands is a result of the interplay between QSE-induced surface roughing and tensile strain-induced smoothening effect.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

End block dynamics in unentangled polymers by dielectric spectroscopy

Dielectric spectroscopy measures the dynamics of polymer melts over a broad frequency range. Developing a theory for the spectral shape can extend the analysis of dielectric spectra beyond determining relaxation times from the peak maxima and adds physical meaning to shape parameters determined with empirical fit functions. Toward this goal, we use the experimental results on unentangled poly(isoprene), and unentangled poly(butylene oxide), polymer melts, to test whether the concept of end blocks could be one reason for the Rouse model deviating from experimental data. These end blocks have been suggested by simulations and neutron spin echo spectroscopy and are a consequence of the monomeric friction coefficient depending on the position of the bead in the chain. The concept of an end block is an approximation which partitions the chain in a middle and two end blocks to avoid overparameterization by a continuous position dependent change of the friction parameter. Analysis of dielectric spectra shows that the deviations of the calculated from the experimental normal mode cannot be related to the end block relaxation. However, the results do not contradict an end block hiding below the segmental relaxation peak. Finally, it seems that the results are compatible with an end block being the specific part of the sub-Rouse chain interpretation close to the chain ends.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Collinear antiferromagnetic order in a quasi-two-dimensional triangular lattice compound DyNiAl 4 Ge 2

The two-dimensional magnetic triangular lattice, with rare-earth local moments coupled through long-range interactions and hosting strong spin-orbit coupling, serves as an ideal platform for exploring novel quantum phenomena. Here, we report the physical properties of DyNiAl 4 Ge 2 single crystals through comprehensive structural, magnetic, heat capacity, electrical transport, and neutron powder diffraction measurements. DyNiAl 4 Ge 2 crystallizes in a trigonal lattice with the space group $R\bar{3}m$, where Dy atoms form two-dimensional triangular lattice layers within the 𝑎𝑏 plane and stack along the 𝑐 axis. Magnetic susceptibility measurement reveals a magnetic phase transition at 𝑇 𝑁 = 9.0 K, further confirmed by the 𝜆-shaped peak in the heat capacity curve. The obtained temperature-magnetic field phase diagram is relatively simple compared to other isostructural compounds, suggesting the absence of magnetic frustration. This is further evidenced by neutron powder diffraction, which shows a collinear antiferromagnetic ground state with a propagation vector 𝒌 = (0, 0, 1.5), implying the absence of magnetic frustration in DyNiAl 4 Ge 2 . In this magnetic structure, Dy moments are aligned ferromagnetically within each triangular layer and arranged antiferromagnetically between adjacent layers along the 𝑐 axis. In conclusion, our results establish DyNiAl 4 Ge 2 as a model system for understanding the interplay between RKKY interaction and geometrical frustration in this family of compounds.

Wu, HengHeng [Chinese Academy of Sciences (CAS), S↗

DOE OSTI2020

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2–3 × 10 −4 (Δ E / E ). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (Δ E / E = 3 × 10 −3 ) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

59 BASIC BIOLOGICAL SCIENCES↗