Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “contrast transfer function”

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 73 records · Page 4

Dynamic Shaping of Grid Response of Multi-Machine Multi-Inverter Systems Through Grid-Forming IBRs

We consider the problem of controlling the frequency response of weakly-coupled multi-machine multi-inverter low-inertia power systems via grid-forming inverter-based resources (IBRs). In contrast to existing methods, our approach relies on dividing the larger system into multiple strongly-coupled subsystems, without ignoring either the underlying network or approximating the subsystem response as an aggregate harmonic mean model. Rather, through a structured clustering and recursive dynamic shaping approach, the frequency response of the overall system to load perturbations is shaped appropriately. We demonstrate the proposed approach for a three-node triangular configuration and a small-scale radial network. Furthermore, previous synchronization analysis for heterogeneous systems requires the machines to satisfy certain proportionality property. In our approach, the effective transfer functions for each cluster can be tuned by the IBRs to satisfy such property, enabling us to apply the shaping control to systems with a wider range of heterogeneous machines.

frequency-shaping control↗

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

Site-specific electronic structure of covalently linked bimetallic dyads from nitrogen K-edge x-ray absorption spectroscopy

A nitrogen K-edge x-ray absorption near-edge structure (XANES) survey is presented for tetrapyrido[3,2-a:',3'-c:3'',2''-h:2''',3'''-j]phenazine (tpphz)-bridged bimetallic assemblies that couple chromophore and catalyst transition metal complexes for light driven catalysis, as well as their individual molecular constituents. We demonstrate the high N site sensitivity of the N pre-edge XANES features, which are energetically well-separated for the phenazine bridge N atoms and for the individual metal-bound N atoms of the inner coordination sphere ligands. By comparison with the time-dependent density functional theory calculated spectra, we determine the origins of these distinguishable spectral features. We find that metal coordination generates large shifts toward higher energy for the metal-bound N atoms, with increasing shift for 3d < 4d < 5d metal bonding. This is attributed to increasing ligand-to-metal σ donation that increases the effective charge of the bound N atoms and stabilizes the N 1s core electrons. In contrast, the phenazine bridge N pre-edge peak is found at a lower energy due to stabilization of the low energy electron accepting orbital localized on the phenazine motif. While no sensitivity to ground state electronic coupling between the individual molecular subunits was observed, the spectra are sensitive to structural distortions of the tpphz bridge. Here these results demonstrate N K-edge XANES as a local probe of electronic structure in large bridging ligand motifs, able to distinctly investigate the ligand-centered orbitals involved in metal-to-ligand and ligand-to-ligand electron transfer following light absorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proteolytic processing induces a conformational switch required for antibacterial toxin delivery

Many Gram-negative bacteria use CdiA effector proteins to inhibit the growth of neighboring competitors. CdiA transfers its toxic CdiA-CT region into the periplasm of target cells, where it is released through proteolytic cleavage. The N-terminal cytoplasm-entry domain of the CdiA-CT then mediates translocation across the inner membrane to deliver the C-terminal toxin domain into the cytosol. Here, we show that proteolysis not only liberates the CdiA-CT for delivery, but is also required to activate the entry domain for membrane translocation. Translocation function depends on precise cleavage after a conserved VENN peptide sequence, and the processed ΔVENN entry domain exhibits distinct biophysical and thermodynamic properties. By contrast, imprecisely processed CdiA-CT fragments do not undergo this transition and fail to translocate to the cytoplasm. These findings suggest that CdiA-CT processing induces a critical structural switch that converts the entry domain into a membrane-translocation competent conformation.

59 BASIC BIOLOGICAL SCIENCES↗

Synthesis, Structure, and Redox Reactivity of Ni Complexes Bearing a Redox and Acid–Base Non-innocent Ligand with Ni II , Ni III , and Ni IV Formal Oxidation States

Here, a series of Ni complexes bearing a redox and acid–base noninnocent tetraamido macrocyclic ligand, H 4 -(TAML-4) {H 4 -(TAML-4) = 15,15-dimethyl-5,8,13,17-tetrahydro-5,8,13,17-tetraaza-dibenzo[a,g]cyclotridecene-6,7,14,16-tetraone}, with formal oxidation states of Ni II , Ni III , and Ni IV were synthesized and characterized structurally and spectroscopically. The X-ray crystallographic analysis of the Ni complexes revealed a square planar geometry, and the [Ni(TAML-4)] complex with the formal oxidation state of Ni IV was characterized to be [Ni III (TAML-4 •+ )] with the oxidation state of the Ni III ion and the one-electron oxidized TAML-4 ligand, TAML-4 •+ . The Ni III oxidation state and the TAML-4 radical cation ligand, TAML-4 •+ , were supported by X-ray absorption spectroscopy and density functional theory calculations. The reversible interconversions between [Ni II (TAML-4)] 2– and [Ni III (TAML-4)] − and between [Ni III (TAML-4)] − and [Ni III (TAML-4 •+ )] were demonstrated in spectroelectrochemical measurements as well as in chemical oxidation and reduction reactions. The reactivities of [Ni III (TAML-4)] − and [Ni III (TAML-4 •+ )] were then investigated in hydride transfer reactions using NADH analogs. Hydride transfer from 9,10-dihydro-10-methylacridine (AcrH 2 ) to [Ni III (TAML-4 •+ )] was found to proceed via electron transfer (ET) from AcrH 2 to [Ni III (TAML-4 •+ )] with no deuterium kinetic isotope effect (k H /k D = 1.0(2)). In contrast, hydride transfer from AcrH 2 to [NiIII(TAML-4)] − proceeded much more slowly via a concerted proton-coupled electron transfer (PCET) process with k H /k D = 7.0(5). In the latter reaction, an electron and a proton were transferred to the Ni III center and the TAML-4 ligand, respectively. The mechanisms of the ET by [Ni III (TAML-4 •+ )] and the concerted PCET by [Ni III (TAML-4)] − were ascribed to the different redox potentials of the Ni complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of Multiple Trichloroethene, cis-Dichloroethene and 1,1-Dichloroethene Degrading Propanotrophic Communities

Aerobic cometabolism offers a viable strategy for the remediation of chlorinated solvent plumes at oxic sites where anaerobic approaches are limited. In this study, propane-enriched mixed cultures (derived from agricultural soils and an impacted site sediment) which previously degraded 1,4-dioxane, were evaluated for their capacity to also degrade trichloroethene (TCE), cis-1,2-dichloroethene (cDCE), and 1,1-dichloroethene (1,1-DCE) over successive transfers. Sustained biodegradation of TCE and cDCE was observed across multiple enrichments, and cultures enriched on one compound generally degraded the other. In contrast, 1,1-DCE biodegradation was restricted to a subset of cultures and removal times increased over transfers. Further, 1,1-DCE removal was absent at elevated concentrations, both trends consistent with inhibitory or toxic effects. Whole genome sequencing analyses revealed pronounced substrate-dependent selection of microbial communities, with cDCE-degrading cultures being dominated by Mycobacterium and Mycolicibacterium, whereas TCE-degrading cultures were dominated by Rhodococcus. Rhodococcus metagenome-assembled genomes (MAGs) in the TCE degrading cultures classified as R. opacus or R. wratislaviensis. 1,1-DCE degrading cultures were dominated by Pseudonocardia, although the associated MAGs contained a truncated propane monooxygenase alpha subunit. Functional gene analysis identified both group 5 (prmABCD) and putative group 6 propane monooxygenases. The following KBase narratives contain the quality controlled reads, MAGs (fasta assemblies) and the prokka annotations for each assembly TCE Site 1A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254918) TCE Soil 2A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254919) TCE Soil T3 A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254920) TCE Soil T4 A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254921) cDCE Site 1A 1B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254915) cDCE Soils T2 A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254927) cDCE Soil T3 A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254928) cDCE Soil 4A and B Propanotrophic MAGs (https://narrative.kbase.us/narrative/254942) 1,1-DCE T2 T3 Propanotrophic MAGs (https://narrative.kbase.us/narrative/254903)

59 BASIC BIOLOGICAL SCIENCES↗

Embedding Reverse Electron Transfer Between Stably Bare Cu Nanoparticles and Cation‐Vacancy CuWO 4

Cu nanoparticles (NPs) have attracted widespread attention in electronics, energy, and catalysis. However, conventionally synthesized Cu NPs face some challenges such as surface passivation and agglomeration in applications, which impairs their functionalities in the physicochemical properties. Here, the issues above by engineering an embedded interface of stably bare Cu NPs on the cation-vacancy CuWO 4 support is addressed, which induces the strong metal-support interactions and reverse electron transfer. Various atomic-scale analyses directly demonstrate the unique electronic structure of the embedded Cu NPs with negative charge and anion oxygen protective layer, which mitigates the typical degradation pathways such as oxidation in ambient air, high-temperature agglomeration, and CO poisoning adsorption. Kinetics and in situ spectroscopic studies unveil that the embedded electron-enriched Cu NPs follow the typical Eley-Rideal mechanism in CO oxidation, contrasting the Langmuir-Hinshelwood mechanism on the traditional Cu NPs. This mechanistic shift is driven by the Coulombic repulsion in anion oxygen layer, enabling its direct reaction with gaseous CO to form the easily desorbed monodentate carbonate.

Cu nanoparticles↗

Transferring a Molecular Foundation Model for Polymer Property Predictions

Transformer-based large language models have remarkable potential to accelerate design optimization for applications such as drug development and material discovery. Self-supervised pretraining of transformer models requires large-scale data sets, which are often sparsely populated in topical areas such as polymer science. Further, state-of-the-art approaches for polymers conduct data augmentation to generate additional samples but unavoidably incur extra computational costs. In contrast, large-scale open-source data sets are available for small molecules and provide a potential solution to data scarcity through transfer learning. In this work, we show that using transformers pretrained on small molecules and fine-tuned on polymer properties achieves comparable accuracy to those trained on augmented polymer data sets for a series of benchmark prediction tasks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural architecture of TolQ-TolR inner membrane protein complex from opportunistic pathogen Acinetobacter baumannii

Gram-negative bacteria harness the proton motive force (PMF) within their inner membrane (IM) to uphold cell envelope integrity, an indispensable aspect for both division and survival. The IM TolQ-TolR complex is the essential part of the Tol-Pal system, serving as a conduit for PMF energy transfer to the outer membrane. Here we present cryo–electron microscopy reconstructions ofAcinetobacter baumanniiTolQ in apo and TolR-bound forms at atomic resolution. The apo TolQ configuration manifests as a symmetric pentameric pore, featuring a transmembrane funnel leading toward a cytoplasmic chamber. In contrast, the TolQ-TolR complex assumes a proton nonpermeable stance, characterized by the TolQ pentamer’s flexure to accommodate the TolR dimer, where two protomers undergo a translation-based relationship. Our structure-guided analysis and simulations support the rotor-stator mechanism of action, wherein the rotation of the TolQ pentamer harmonizes with the TolR protomers’ interplay. These findings broaden our mechanistic comprehension of molecular stator units empowering critical functions within the Gram-negative bacterial cell envelope.

Science & Technology - Other Topics↗

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗

Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. Here, we leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T 1/2 ) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T 1/2 values. We then compare ML T 1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T 1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T 1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T 1/2 values with a Pearson’s r of 0.82. In contrast, DFA-predicted T 1/2 values have large errors and show no correlation to estimated T 1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plasmon Dynamics Driven by Aggregation of Tris(2,2′-bipyridine)ruthenium(II)-Functionalized Gold Nanoparticles Probed by XANES and Transient Absorption Spectroscopy

Energy conversion dynamics is critical for advancing next-generation photovoltaics, optoelectronics, and light-harvesting technologies. Noble metal plasmonic nanoparticles play a pivotal role as nanoscale electromagnetic confinement structures, driving photon-induced chemical reactions. In this study, we explore the effects of [Ru(bpy) 3 ] 2+ functionalization and aggregation on citrate-capped gold nanoparticles (AuNPs) of 40 and 100 nm diameters, focusing on molecule-plasmon interactions and their influence on electronic and energy dissipation properties. X-ray absorption near-edge spectroscopy (XANES) revealed that [Ru(bpy) 3 ] 2+ functionalization induces controlled aggregation without altering the oxidation state of gold. A more pronounced white-line intensity is observed in 40 nm AuNPs, consistent with greater s–p–d hybridization and a higher density of surface states, likely influenced by both nanoparticle size and aggregation. Transient absorption (TA) spectroscopy highlights faster electron–phonon relaxation dynamics in aggregated 40 nm nanoparticles, which is attributed to increased electron delocalization and more efficient coupling to the phonon bath. In contrast, 100 nm nanoparticles exhibit minimal changes due to a lower degree of aggregation. Interestingly, we observe that enhanced electron–phonon coupling in aggregated nanoparticles coincides with a slowing of electron–electron scattering. These observations suggest a competitive interplay between the two relaxation pathways, where enhanced energy transfer to the lattice in aggregated systems can suppress electronic thermalization. In conclusion, these findings underscore the critical role of nanoparticle size, aggregation, and molecule–surface interactions in modulating plasmonic dynamics and excited-state lifetimes and further provide valuable insights into designing tailored plasmonic systems with transformative potential for sensing, catalysis, and energy conversion.

36 MATERIALS SCIENCE↗

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

Molecular Geometry and Electronic Structure of Copper Corroles

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree–Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S o ) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T o ), which is approximately 5 kcal/mol higher in energy compared to the S o state for both the bare and substituted copper corroles. Notably, the planarity of the T o geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S o and T o geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D 3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Examination of the sensitivity of quasifree reactions to details of the bound-state overlap functions

It is often stated that heavy-ion nucleon knockout reactions are mostly sensitive to the tails of the bound-state wave functions. In contrast, (p,2p) and (p,pn) reactions are known to access information on the full overlap functions within the nucleus. We analyze the oxygen isotopic chain and explore the differences between single-particle wave functions generated with potential models, used in the experimental analysis of knockout reactions, and ab initio computations from self-consistent Green's function theory. Contrary to common belief, we find that not only the tail of the overlap functions, but also their internal part is assessed in both reaction mechanisms, which are crucial to yield accurately determined spectroscopic information.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Isolation, bonding and reactivity of a monomeric stibine oxide

Abstract In contrast to phosphine oxides and arsine oxides, which are common and exist as stable monomeric species featuring the corresponding pnictoryl functional group (Pn=O/Pn + –O − ; Pn = P, As), stibine oxides are generally polymeric, and the properties of the unperturbed stiboryl group (Sb=O/Sb + –O − ) remain unexplored. We now report the isolation of the monomeric stibine oxide, Dipp 3 SbO (where Dipp = 2,6-diisopropylphenyl). Spectroscopic, crystallographic and computational studies provide insight into the nature of the Sb=O/Sb + –O − bond. Moreover, isolation of Dipp 3 SbO allows the chemistry of the stiboryl group to be explored. Here we show that Dipp 3 SbO can act as a Brønsted base, a hydrogen-bond acceptor and a transition-metal ligand, in addition engaging in 1,2-addition, O-for-F 2 exchange and O-atom transfer. In all cases, the reactivity of Dipp 3 SbO differed from that of the lighter congeners Dipp 3 AsO and Dipp 3 PO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational models of direct and indirect X‐ray breast imaging detectors for in silico trials

Abstract Background To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed. Purpose The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems. Methods We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open‐source GPU‐accelerated Monte Carlo x‐ray imaging simulation code (MC‐GPU). We build on the previous version of the MC‐GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)‐approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x‐ray photon till its first interaction within the x‐ray detector. On the other hand, the signal spread in the x‐ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron–hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth‐dependent spread and collection of optical photons and relies on a pre‐computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x‐ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data. Results We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors. Conclusions We have presented computational models of three commercially available FDA‐approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC‐GPU codes that can be used to replicate three systems are available in open source format through GitHub.

Sengupta, Aunnasha↗