Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

The feasibility study and characterization of a two-dimensional diode array in “magic phantom” for high dose rate brachytherapy quality assurance

Purpose: High dose rate (HDR) brachytherapy is a radiation treatment technique capable of delivering large dose rates to the tumor. Radiation is delivered using remote afterloaders to drive highly active sources (commonly {sup 192}Ir with an air KERMA strength range between 20 000 and 40 000 U, where 1 U = 1 μGy m{sup 2}/h in air) through applicators directly into the patient's prescribed region of treatment. Due to the obvious ramifications of incorrect treatment while using such an active source, it is essential that there are methods for quality assurance (QA) that can directly and accurately verify the treatment plan and the functionality of the remote afterloader. This paper describes the feasibility study of a QA system for HDR brachytherapy using a phantom based two-dimensional 11 × 11 epitaxial diode array, named “magic phantom.”Methods: The HDR brachytherapy treatment plan is translated to the phantom with two rows of 10 (20 in total) HDR source flexible catheters, arranged above and below the diode array “magic plate” (MP). Four-dimensional source tracking in each catheter is based upon a developed fast iterative algorithm, utilizing the response of the diodes in close proximity to the {sup 192}Ir source, sampled at 100 ms intervals by a fast data acquisition (DAQ) system. Using a {sup 192}Ir source in a solid water phantom, the angular response of the developed epitaxial diodes utilized in the MP and also the variation of the MP response as a function of the source-to-detector distance (SDD) were investigated. These response data are then used by an iterative algorithm for source dwelling position determination. A measurement of the average transit speed between dwell positions was performed using the diodes and a fast DAQ.Results: The angular response of the epitaxial diode showed a variation of 15% within 360°, with two flat regions above and below the detector face with less than 5% variation. For SDD distances of between 5 and 30 mm the relative response of the epitaxial diodes used in the MP is in good agreement (within 8%) with radial dose function measurements found within the TG-43 protocol, with SDD of up to 70 mm showing a 40% over response. A method for four-dimensional localization of the HDR source was developed, allowing the source dwell position to be derived within 0.50 mm of the expected position. An estimation of the average transit speed for varying step sizes was determined and was found to increase from (12.8 ± 0.3) up to (38.6 ± 0.4) cm/s for a step size of 2.5 and 50 mm, respectively.Conclusions: Our characterization of the designed QA “magic phantom” with MP in realistic HDR photon fields demonstrates the promising performance for real-time source position tracking in four dimensions and measurements of transit times. Further development of this system will allow a full suite for QA in HDR brachytherapy and analysis, and for future in vivo tracking.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Capturing the 'Representative Person' Using Exposure Equations and Input Variable Distributions - 18416

Radiological performance assessments (PAs) for radioactive waste employ the concept of future human receptors, or persons potentially exposed to radiation releases. The International Commission on Radiological Protection (ICRP) has defined the 'representative person' as a receptor, to replace the earlier concept of an 'average member of the critical group', for use in individual dose assessments and comparisons to regulatory criteria. A representative person can be assessed in PAs by defining a critical group, which is a relatively highly exposed subpopulation; and then estimating average exposures, doses, and risks for that group. Probabilistic PA models employ probability distributions to characterize uncertainty in exposure characteristics/variables at a scale aligned with contaminant transport and other aspects of the PA model. The concept of a representative person from a critical group, in the context of the first stage of an iterative approach and a probabilistic PA model, is illustrated using a receptor scenario and a set of exposure equations that incorporate probability distributions for input variables. Different scenarios may require different equations and distributions, but the approach taken is broadly generalizable. PA practitioners and regulators may benefit from such a generalizable method, as it may simplify approaches as well as facilitate reviews and comparisons between different PAs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Scalable Asynchronous Domain Decomposition Solvers

We discuss how parallel implementations of linear iterative solvers generally alternate between phases of data exchange and phases of local computation. Increasingly large problem sizes and more heterogeneous compute architectures make load balancing and the design of low latency network interconnects that are able to satisfy the communication requirements of linear solvers very challenging tasks. In particular, global communication patterns such as inner products become increasingly limiting at scale. We explore the use of asynchronous communication based on one-sided Message Passing Interface primitives in the context of domain decomposition solvers. In particular, a scalable asynchronous two-level Schwarz method is presented. We discuss practical issues encountered in the development of a scalable solver and show experimental results obtained on a state-of-the-art supercomputer system that illustrate the benefits of asynchronous solvers in load balanced as well as load imbalanced scenarios. Using the novel method, we can observe speedups of up to four times over its classical synchronous equivalent.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

TU-EF-304-07: Monte Carlo-Based Inverse Treatment Plan Optimization for Intensity Modulated Proton Therapy

Purpose: Intensity-modulated proton therapy (IMPT) is increasingly used in proton therapy. For IMPT optimization, Monte Carlo (MC) is desired for spots dose calculations because of its high accuracy, especially in cases with a high level of heterogeneity. It is also preferred in biological optimization problems due to the capability of computing quantities related to biological effects. However, MC simulation is typically too slow to be used for this purpose. Although GPU-based MC engines have become available, the achieved efficiency is still not ideal. The purpose of this work is to develop a new optimization scheme to include GPU-based MC into IMPT. Methods: A conventional approach using MC in IMPT simply calls the MC dose engine repeatedly for each spot dose calculations. However, this is not the optimal approach, because of the unnecessary computations on some spots that turned out to have very small weights after solving the optimization problem. GPU-memory writing conflict occurring at a small beam size also reduces computational efficiency. To solve these problems, we developed a new framework that iteratively performs MC dose calculations and plan optimizations. At each dose calculation step, the particles were sampled from different spots altogether with Metropolis algorithm, such that the particle number is proportional to the latest optimized spot intensity. Simultaneously transporting particles from multiple spots also mitigated the memory writing conflict problem. Results: We have validated the proposed MC-based optimization schemes in one prostate case. The total computation time of our method was ∼5–6 min on one NVIDIA GPU card, including both spot dose calculation and plan optimization, whereas a conventional method naively using the same GPU-based MC engine were ∼3 times slower. Conclusion: A fast GPU-based MC dose calculation method along with a novel optimization workflow is developed. The high efficiency makes it attractive for clinical usages.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-F-BRA-02: Electromagnetic Tracking in Brachytherapy as An Advanced Modality for Treatment Quality Assurance

Purpose: To present the use of Electromagnetic Tracking (EMT) for quality assurance in brachytherapy by means of phantom studies and to assess the clinical applicability of EMT during HDR breast brachytherapy. Methods: An EMT system was investigated to examine its suitability for clinical applications in brachytherapy. A field generator served as electromagnetic field emitter. Sensors (magnetic sensitive only), connected to a control unit, were used and their respective position and orientation inside a pre-defined measurement volume (500 mm cube length) determined. Up to three 6DoF sensors were placed on the phantom’s surface to obtain additional reference coordinates used to derive relative measured positions of a smaller 5DoF sensor inserted in the 6F catheters of the implant. The catheters were successively measured by manual displacement of the sensor at ∼40 mm/s. The measured catheter tracks, acquired multiple times at various locations (CT and treatment room), were smoothed, divided into intervals (2.5 mm dwell step size), registered (rigid Iterative Closest Point transformation) and compared against the known phantom geometry. Results: The reference coordinates were used to exclude the influence of external (e.g., respiratory-induced) motion. Precision tests in a clinical setting showed variances below 1 mm (translational) and 1° (rotational), respectively. Our method for catheter reconstruction preserved the length of the tracked catheter (within 1 mm). The measured tracking accuracy was 1±0.3 mm (maximum: 2 mm). The results are less accurate in environments potentially interfering with the magnetic field, e.g., in the vicinity of ferromagnetic table components. Conclusion: Our EMT system is able to perform reproducible and accurate catheter tracking and reconstruction. Currently, measurements of the implant geometry in HDR breast treatments are initiated. Online implant monitoring by means of EM tracking may be a first step towards advanced brachytherapy treatment QA. The study was financially supported by ELEKTA | Nucletron. There are no further disclosures and conflicts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗
DOE OSTI2013

Solving the parquet equations for the Hubbard model beyond weak coupling

We find that imposing crossing symmetry in the iteration process considerably extends the range of convergence for solutions of the parquet equations for the Hubbard model. When crossing symmetry is not imposed, the convergence of both simple iteration and more complicated continuous loading (homotopy) methods is limited to high temperatures and weak interactions. We modify the algorithm to impose the crossing symmetry without increasing the computational complexity. We also imposed time reversal and a subset of the point group symmetries, but they did not further improve the convergence. We elaborate the details of the latency hiding scheme which can significantly improve the performance in the computational implementation. With these modifications, stable solutions for the parquet equations can be obtained by iteration more quickly even for values of the interaction that are a significant fraction of the bandwidth and for temperatures that are much smaller than the bandwidth. This may represent a crucial step towards the solution of two-particle field theories for correlated electron models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A promising limited angular computed tomography reconstruction via segmentation based regional enhancement and total variation minimization

X-ray computed tomography (CT) is a powerful and common inspection technique used for the industrial non-destructive testing. However, large-sized and heavily absorbing objects cause the formation of artifacts because of either the lack of specimen penetration in specific directions or the acquisition of data from only a limited angular range of views. Although the sparse optimization-based methods, such as the total variation (TV) minimization method, can suppress artifacts to some extent, reconstructing the images such that they converge to accurate values remains difficult because of the deficiency in continuous angular data and inconsistency in the projections. To address this problem, we use the idea of regional enhancement of the true values and suppression of the illusory artifacts outside the region to develop an efficient iterative algorithm. This algorithm is based on the combination of regional enhancement of the true values and TV minimization for the limited angular reconstruction. In this algorithm, the segmentation approach is introduced to distinguish the regions of different image knowledge and generate the support mask of the image. A new regularization term, which contains the support knowledge to enhance the true values of the image, is incorporated into the objective function. Then, the proposed optimization model is solved by variable splitting and the alternating direction method efficiently. A compensation approach is also designed to extract useful information from the initial projections and thus reduce false segmentation result and correct the segmentation support and the segmented image. The results obtained from comparing both simulation studies and real CT data set reconstructions indicate that the proposed algorithm generates a more accurate image than do the other reconstruction methods. The experimental results show that this algorithm can produce high-quality reconstructed images for the limited angular reconstruction and suppress the illusory artifacts caused by the deficiency in valid data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-F-I-54: Spatial Resolution Studies in Proton CT Using a Phase-II Prototype Head Scanner

Purpose: To characterize the modulation transfer function (MTF) of the pre-clinical (phase II) head scanner developed for proton computed tomography (pCT) by the pCT collaboration. To evaluate the spatial resolution achievable by this system. Methods: Our phase II proton CT scanner prototype consists of two silicon telescopes that track individual protons upstream and downstream from a phantom, and a 5-stage scintillation detector that measures a combination of the residual energy and range of the proton. Residual energy is converted to water equivalent path length (WEPL) of the protons in the scanned object. The set of WEPL values and associated paths of protons passing through the object over a 360° angular scan is processed by an iterative parallelizable reconstruction algorithm that runs on GP-GPU hardware. A custom edge phantom composed of water-equivalent polymer and tissue-equivalent material inserts was constructed. The phantom was first simulated in Geant4 and then built to perform experimental beam tests with 200 MeV protons at the Northwestern Medicine Chicago Proton Center. The oversampling method was used to construct radial and azimuthal edge spread functions and modulation transfer functions. The spatial resolution was defined by the 10% point of the modulation transfer function in units of lp/cm. Results: The spatial resolution of the image was found to be strongly correlated with the radial position of the insert but independent of the relative stopping power of the insert. The spatial resolution varies between roughly 4 and 6 lp/cm in both the the radial and azimuthal directions depending on the radial displacement of the edge. Conclusion: The amount of image degradation due to our detector system is small compared with the effects of multiple Coulomb scattering, pixelation of the image and the reconstruction algorithm. Improvements in reconstruction will be made in order to achieve the theoretical limits of spatial resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Imaging of cardiac perfusion of free-breathing small animals using dynamic phase-correlated micro-CT

Purpose:Mouse models of cardiac diseases have proven to be a valuable tool in preclinical research. The high cardiac and respiratory rates of free breathing mice prohibit conventional in vivo cardiac perfusion studies using computed tomography even if gating methods are applied. This makes a sacrification of the animals unavoidable and only allows for the application of ex vivo methods. Methods: To overcome this issue the authors propose a low dose scan protocol and an associated reconstruction algorithm that allows for in vivo imaging of cardiac perfusion and associated processes that are retrospectively synchronized to the respiratory and cardiac motion of the animal. The scan protocol consists of repetitive injections of contrast media within several consecutive scans while the ECG, respiratory motion, and timestamp of contrast injection are recorded and synchronized to the acquired projections. The iterative reconstruction algorithm employs a six-dimensional edge-preserving filter to provide low-noise, motion artifact-free images of the animal examined using the authors' low dose scan protocol. Results: The reconstructions obtained show that the complete temporal bolus evolution can be visualized and quantified in any desired combination of cardiac and respiratory phase including reperfusion phases. The proposed reconstruction method thereby keeps the administered radiation dose at a minimum and thus reduces metabolic inference to the animal allowing for longitudinal studies. Conclusions: The authors' low dose scan protocol and phase-correlated dynamic reconstruction algorithm allow for an easy and effective way to visualize phase-correlated perfusion processes in routine laboratory studies using free-breathing mice.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Characterization of a constrained paired-view technique in iterative reconstruction for breast tomosynthesis

Purpose: The order in which the projection views are employed in the reconstruction of tomosynthesis by iterative algorithms, such as simultaneous algebraic reconstruction technique and maximum likelihood, has a strong effect on the rate of convergence, accuracy, and the edge-blurring artifacts in the reconstructed image. The purpose of this investigation was to characterize and evaluate the effects of ordering schemes on image quality for breast tomosynthesis reconstruction and to explore a new constrained paired-view technique that could provide reduction of reconstruction artifacts. In this work, the authors compared several different ordering schemes and characterized the image quality and the formation of out-of-plane artifacts. Furthermore, a new normalization method is presented. It produces more accurate reconstructions with reduced artifacts comparing to the standard method of sequential ordering.Methods: In addition to visual assessment of image quality, several indices such as the signal-difference-to-noise ratio, the artifact-spread function, and the lesion detectability (d{sup ′}) were computed to quantitatively evaluate the effect of ordering scheme. The sets of breast tomosynthesis projection images were simulated for reconstruction; one set had uniform background (white noise only) and the other two contained both anatomic background and quantum noise. Clinical breast images were also studied for comparison.Results: The authors have quantified the image quality in reconstructed slices for a range of tumor sizes. The authors’ proposed method provides better performance for all of the metrics tested (contrast, d{sup ′}, and the level of artifacts) both for the uniform phantom case and in the presence of anatomical structure.Conclusions: The paired projection normalization provides better performance in the image quality of the reconstructed slices, and results in a lower level of artifacts in the Z direction. This implies that even a relatively simple method like the “side-to-side” sequence, which pairs two symmetrical projections with equal angular distance from the central projection, would achieve better reconstructed image quality than the conventional “step-by-step” method, which uses sequential projections one after another.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Automated segmentation of dental CBCT image with prior-guided sequential random forests

Purpose: Cone-beam computed tomography (CBCT) is an increasingly utilized imaging modality for the diagnosis and treatment planning of the patients with craniomaxillofacial (CMF) deformities. Accurate segmentation of CBCT image is an essential step to generate 3D models for the diagnosis and treatment planning of the patients with CMF deformities. However, due to the image artifacts caused by beam hardening, imaging noise, inhomogeneity, truncation, and maximal intercuspation, it is difficult to segment the CBCT. Methods: In this paper, the authors present a new automatic segmentation method to address these problems. Specifically, the authors first employ a majority voting method to estimate the initial segmentation probability maps of both mandible and maxilla based on multiple aligned expert-segmented CBCT images. These probability maps provide an important prior guidance for CBCT segmentation. The authors then extract both the appearance features from CBCTs and the context features from the initial probability maps to train the first-layer of random forest classifier that can select discriminative features for segmentation. Based on the first-layer of trained classifier, the probability maps are updated, which will be employed to further train the next layer of random forest classifier. By iteratively training the subsequent random forest classifier using both the original CBCT features and the updated segmentation probability maps, a sequence of classifiers can be derived for accurate segmentation of CBCT images. Results: Segmentation results on CBCTs of 30 subjects were both quantitatively and qualitatively validated based on manually labeled ground truth. The average Dice ratios of mandible and maxilla by the authors’ method were 0.94 and 0.91, respectively, which are significantly better than the state-of-the-art method based on sparse representation (p-value < 0.001). Conclusions: The authors have developed and validated a novel fully automated method for CBCT segmentation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2011

Resolving all-order method convergence problems for atomic physics applications

The development of the relativistic all-order method where all single, double, and partial triple excitations of the Dirac-Hartree-Fock wave function are included to all orders of perturbation theory led to many important results for the study of fundamental symmetries, development of atomic clocks, ultracold atom physics, and others, as well as provided recommended values of many atomic properties critically evaluated for their accuracy for a large number of monovalent systems. This approach requires iterative solutions of the linearized coupled-cluster equations leading to convergence issues in some cases where correlation corrections are particularly large or lead to an oscillating pattern. Moreover, these issues also lead to similar problems in the configuration-interaction (CI)+all-order method for many-particle systems. In this work, we have resolved most of the known convergence problems by applying two different convergence stabilizer methods, namely, reduced linear equation and direct inversion of iterative subspace. Examples are presented for B, Al, Zn{sup +}, and Yb{sup +}. Solving these convergence problems greatly expands the number of atomic species that can be treated with the all-order methods and is anticipated to facilitate many interesting future applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

A general framework and review of scatter correction methods in x-ray cone-beam computerized tomography. Part 1: Scatter compensation approaches

Since scattered radiation in cone-beam volume CT implies severe degradation of CT images by quantification errors, artifacts, and noise increase, scatter suppression is one of the main issues related to image quality in CBCT imaging. The aim of this review is to structurize the variety of scatter suppression methods, to analyze the common structure, and to develop a general framework for scatter correction procedures. In general, scatter suppression combines hardware techniques of scatter rejection and software methods of scatter correction. The authors emphasize that scatter correction procedures consist of the main components scatter estimation (by measurement or mathematical modeling) and scatter compensation (deterministic or statistical methods). The framework comprises most scatter correction approaches and its validity also goes beyond transmission CT. Before the advent of cone-beam CT, a lot of papers on scatter correction approaches in x-ray radiography, mammography, emission tomography, and in Megavolt CT had been published. The opportunity to avail from research in those other fields of medical imaging has not yet been sufficiently exploited. Therefore additional references are included when ever it seems pertinent. Scatter estimation and scatter compensation are typically intertwined in iterative procedures. It makes sense to recognize iterative approaches in the light of the concept of self-consistency. The importance of incorporating scatter compensation approaches into a statistical framework for noise minimization has to be underscored. Signal and noise propagation analysis is presented. A main result is the preservation of differential-signal-to-noise-ratio (dSNR) in CT projection data by ideal scatter correction. The objective of scatter compensation methods is the restoration of quantitative accuracy and a balance between low-contrast restoration and noise reduction. In a synopsis section, the different deterministic and statistical methods are discussed with respect to their properties and applications. The current paper is focused on scatter compensation algorithms. The multitude of scatter estimation models will be dealt with in a separate paper.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Evaluating iterative reconstruction performance in computed tomography

Purpose: Iterative reconstruction (IR) offers notable advantages in computed tomography (CT). However, its performance characterization is complicated by its potentially nonlinear behavior, impacting performance in terms of specific tasks. This study aimed to evaluate the performance of IR with both task-specific and task-generic strategies. Methods: The performance of IR in CT was mathematically assessed with an observer model that predicted the detection accuracy in terms of the detectability index (d′). d′ was calculated based on the properties of the image noise and resolution, the observer, and the detection task. The characterizations of image noise and resolution were extended to accommodate the nonlinearity of IR. A library of tasks was mathematically modeled at a range of sizes (radius 1–4 mm), contrast levels (10–100 HU), and edge profiles (sharp and soft). Unique d′ values were calculated for each task with respect to five radiation exposure levels (volume CT dose index, CTDI{sub vol}: 3.4–64.8 mGy) and four reconstruction algorithms (filtered backprojection reconstruction, FBP; iterative reconstruction in imaging space, IRIS; and sinogram affirmed iterative reconstruction with strengths of 3 and 5, SAFIRE3 and SAFIRE5; all provided by Siemens Healthcare, Forchheim, Germany). The d′ values were translated into the areas under the receiver operating characteristic curve (AUC) to represent human observer performance. For each task and reconstruction algorithm, a threshold dose was derived as the minimum dose required to achieve a threshold AUC of 0.9. A task-specific dose reduction potential of IR was calculated as the difference between the threshold doses for IR and FBP. A task-generic comparison was further made between IR and FBP in terms of the percent of all tasks yielding an AUC higher than the threshold. Results: IR required less dose than FBP to achieve the threshold AUC. In general, SAFIRE5 showed the most significant dose reduction potentials (11–54 mGy, 77%–84%), followed by SAFIRE3 (7–36 mGy, 50%–61%) and IRIS (6–26 mGy, 37%–50%). The dose reduction potentials highly depended on task size and task contrast, with tasks of lower contrasts and smaller sizes, i.e., more challenging tasks, indicating higher dose reductions. Softer edge profile showed higher dose reduction potentials with SAFIRE3 and SAFIRE5, but not with IRIS. The task-generic comparison between IR and FBP demonstrated the overall superiority of IR performance, as IR allowed a larger percent of tasks to exceed the threshold AUC: IRIS, 8%–12%; SAFIRE3, 10%–16%; and SAFIRE5, 20%–33%. The improvement with IR was generally more pronounced at lower dose levels. Conclusions: Expanding beyond traditional contrast and noise based assessments of IR, we performed both task-specific and task-generic evaluations of IR performance. The task-specific evaluation demonstrated the dependency of IR’s dose reduction potential on task attributes, which can be employed to optimize IR for clinical indications with specific range of size and contrast. The task-generic evaluation demonstrated IR’s overall superiority over FBP in terms of the range of tasks exceeding a threshold performance level, which can be employed for general comparisons between algorithms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

The current capacity of computers makes it possible to perform simulations of small systems with portable, explicit-solvent potentials achieving high degree of accuracy. However, simplified models must be employed to exploit the behavior of large systems or to perform systematic scans of smaller systems. While powerful algorithms are available to facilitate the sampling of the conformational space, successful applications of such models are hindered by the availability of simple enough potentials able to satisfactorily reproduce known properties of the system. We develop an interatomic potential to account for a number of properties of proteins in a computationally economic way. The potential is defined within an all-atom, implicit solvent model by contact functions between the different atom types. The associated numerical values can be optimized by an iterative Monte Carlo scheme on any available experimental data, provided that they are expressible as thermal averages of some conformational properties. We test this model on three different proteins, for which we also perform a scan of all possible point mutations with explicit conformational sampling. The resulting models, optimized solely on a subset of native distances, not only reproduce the native conformations within a few Angstroms from the experimental ones, but show the cooperative transition between native and denatured state and correctly predict the measured free-energy changes associated with point mutations. Moreover, differently from other structure-based models, our method leaves a residual degree of frustration, which is known to be present in protein molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Technical Note: Evaluation of a 160-mm/256-row CT scanner for whole-heart quantitative myocardial perfusion imaging

Purpose: The authors investigated the performance of a recently introduced 160-mm/256-row CT system for low dose quantitative myocardial perfusion (MP) imaging of the whole heart. This platform is equipped with a gantry capable of rotating at 280 ms per full cycle, a second generation of adaptive statistical iterative reconstruction (ASiR-V) to correct for image noise arising from low tube voltage potential/tube current dynamic scanning, and image reconstruction algorithms to tackle beam-hardening, cone-beam, and partial-scan effects. Methods: Phantom studies were performed to investigate the effectiveness of image noise and artifact reduction with a GE Healthcare Revolution CT system for three acquisition protocols used in quantitative CT MP imaging: 100, 120, and 140 kVp/25 mAs. The heart chambers of an anthropomorphic chest phantom were filled with iodinated contrast solution at different concentrations (contrast levels) to simulate the circulation of contrast through the heart in quantitative CT MP imaging. To evaluate beam-hardening correction, the phantom was scanned at each contrast level to measure the changes in CT number (in Hounsfield unit or HU) in the water-filled region surrounding the heart chambers with respect to baseline. To evaluate cone-beam artifact correction, differences in mean water HU between the central and peripheral slices were compared. Partial-scan artifact correction was evaluated from the fluctuation of mean water HU in successive partial scans. To evaluate image noise reduction, a small hollow region adjacent to the heart chambers was filled with diluted contrast, and contrast-to-noise ratio in the region before and after noise correction with ASiR-V was compared. The quality of MP maps acquired with the CT system was also evaluated in porcine CT MP studies. Myocardial infarct was induced in a farm pig from a transient occlusion of the distal left anterior descending (LAD) artery with a catheter-based interventional procedure. MP maps were generated from the dynamic contrast-enhanced (DCE) heart images taken at baseline and three weeks after the ischemic insult. Results: Their results showed that the phantom and animal images acquired with the CT platform were minimally affected by image noise and artifacts. For the beam-hardening phantom study, changes in water HU in the wall surrounding the heart chambers greatly reduced from >±30 to ≤ ± 5 HU at all kVp settings except one region at 100 kVp (7 HU). For the cone-beam phantom study, differences in mean water HU from the central slice were less than 5 HU at two peripheral slices with each 4 cm away from the central slice. These findings were reproducible in the pig DCE images at two peripheral slices that were 6 cm away from the central slice. For the partial-scan phantom study, standard deviations of the mean water HU in 10 successive partial scans were less than 5 HU at the central slice. Similar observations were made in the pig DCE images at two peripheral slices with each 6 cm away from the central slice. For the image noise phantom study, CNRs in the ASiR-V images were statistically higher (p < 0.05) than the non-ASiR-V images at all kVp settings. MP maps generated from the porcine DCE images were in excellent quality, with the ischemia in the LAD territory clearly seen in the three orthogonal views. Conclusions: The study demonstrates that this CT system can provide accurate and reproducible CT numbers during cardiac gated acquisitions across a wide axial field of view. This CT number fidelity will enable this imaging tool to assess contrast enhancement, potentially providing valuable added information beyond anatomic evaluation of coronary stenoses. Furthermore, their results collectively suggested that the 100 kVp/25 mAs protocol run on this CT system provides sufficient image accuracy at a low radiation dose (<3 mSv) for whole-heart quantitative CT MP imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Iterative Joint Ptychography-Tomography with Total Variation Regularization

In order to determine the 3D structure of a thick sample, researchers have recently combined ptychography (for high resolution) and tomography (for 3D imaging) in a single experiment. 2-step methods are usually adopted for reconstruction, where the ptychography and tomography problems are often solved independently. In this paper, we provide a novel model and ADMM-based algorithm to jointly solve the ptychography-tomography problem iteratively, also employing total variation regularization. The proposed method permits large scan stepsizes for the ptychography experiment, requiring less measurements and being more robust to noise with respect to other strategies, while achieving higher reconstruction quality results.

2-step methods↗