SEARCH · Engineering Papers

Results for “PRESERVATION”

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 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2017

A Real-Time Greedy-Index Dispatching Policy for using PEVs to Provide Frequency Regulation Service

This article presents a real-time greedy-index dispatching policy (GIDP) for using plug-in electric vehicles (PEVs) to provide frequency regulation services. A new service cost allocation mechanism is proposed to award PEVs based on the amount of service they provided, while considering compensations for delayed-charging and reduction of battery lifetime due to participation of the service. The GIDP transforms the optimal dispatch problem from a high-dimensional space into a one-dimensional space while preserving the solution optimality. When solving the transformed problem in real-time, the global optimality of the GIDP solution can be guaranteed by mathematically proved “indexability”. Because the GIDP index can be calculated upon the PEV’s arrival and used for the entire decision making process till its departure, the computational burden is minimized and the complexity of the aggregator dispatch process is significantly reduced. Finally, simulation results are used to evaluate the proposed GIDP, and to demonstrate the potential profitability from providing frequency regulation service by using PEVs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2015

Carbon and sulfur isotopic signatures of ancient life and environment at the microbial scale: Neoarchean shales and carbonates

An approach to coordinated, spatially resolved, in situ carbon isotope analysis of organic matter and carbonate minerals, and sulfur three- and four-isotope analysis of pyrite with an unprecedented combination of spatial resolution, precision, and accuracy is described. Organic matter and pyrite from eleven rock samples of Neoarchean drill core express nearly the entire range of δ 13 C, δ 34 S, Δ 33 S, and Δ 36 S known from the geologic record, commonly in correlation with morphology, mineralogy, and elemental composition. A new analytical approach (including a set of organic calibration standards) to account for a strong correlation between H/C and instrumental bias in SIMS δ 13 C measurement of organic matter is identified. Small (2–3 μm) organic domains in carbonate matrices are analyzed with sub-permil accuracy and precision. Separate 20- to 50-μm domains of kerogen in a single ~0.5 cm3 sample of the ~2.7 Ga Tumbiana Formation have δ 13 C = -52.3 ± 0.1‰ and -34.4 ± 0.1‰, likely preserving distinct signatures of methanotrophy and photoautotrophy. Pyrobitumen in the ~2.6 Ga Jeerinah Formation and the ~2.5 Ga Mount McRae Shale is systematically 13 C-enriched relative to co-occurring kerogen, and associations with uraniferous mineral grains suggest radiolytic alteration. A large range in sulfur isotopic compositions (including higher Δ 33 S and more extreme spatial gradients in Δ 33 S and Δ 36 S than any previously reported) are observed in correlation with morphology and associated mineralogy. Changing systematics of δ 34 S, Δ 33 S, and Δ 36 S, previously investigated at the millimeter to centimeter scale using bulk analysis, are shown to occur at the micrometer scale of individual pyrite grains. These results support the emerging view that the dampened signature of mass-independent sulfur isotope fractionation (S-MIF) associated with the Mesoarchean continued into the early Neoarchean, and that the connections between methane and sulfur metabolism affected the production and preservation of S-MIF during the first half of the planet's history.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

Biogenic manganese oxides as reservoirs of organic carbon and proteins in terrestrial and marine environments

Manganese (Mn) oxides participate in a range of interactions with organic carbon (OC) that can lead to either carbon degradation or preservation. In this work, we examine the abundance and composition of OC associated with biogenic and environmental Mn oxides to elucidate the role of Mn oxides as a reservoir for carbon and their potential for selective partitioning of particular carbon species. Mn oxides precipitated in natural brackish waters and by Mn(II)-oxidizing marine bacteria and terrestrial fungi harbor considerable levels of organic carbon (4.1-17.0 mol OC per kg mineral) compared to ferromanganese cave deposits which contain 1-2 orders of magnitude lower OC. Spectroscopic analyses indicate that the chemical composition of Mn oxide-associated OC from microbial cultures is homogeneous with bacterial Mn oxides hosting primarily proteinaceous carbon and fungal Mn oxides containing both protein- and lipopolysaccharide-like carbon. The bacterial Mn oxide-hosted proteins are involved in both Mn(II) oxidation and metal binding by these bacterial species and could be involved in the mineral nucleation process as well. By comparison, the composition of OC associated with Mn oxides formed in natural settings (brackish waters and particularly in cave ferromanganese rock coatings) is more spatially and chemically heterogeneous. Cave Mn oxide-associated organic material is enriched in aliphatic C, which together with the lower carbon concentrations, points to more extensive microbial or mineral processing of carbon in this system relative to the other systems examined in this study, and as would be expected in oligotrophic cave environments. This study highlights Mn oxides as a reservoir for carbon in varied environments. Lastly, the presence and in some cases dominance of proteinaceous carbon within the biogenic and natural Mn oxides may contribute to preferential preservation of proteins in sediments and dominance of protein-dependent metabolisms in the subsurface biosphere.

58 GEOSCIENCES↗

DOE OSTI2016

Oxygen isotope thermometry using quartz inclusions in garnet

Oxygen isotope ratios of quartz inclusions (QI) within garnet from granulite and amphibolite facies gneisses in the Adirondack Mountains, NY were analysed and used to determine metamorphic temperatures. Primary QI for eight of 12 samples have δ 18 O values significantly lower than matrix quartz (MQ). The primary QI retain δ 18 O values representative of thermal conditions during garnet crystallization, whereas the δ 18 O values of MQ were raised by diffusive exchange with other matrix minerals (e.g. mica and feldspar) during cooling. The δ 18 O differences between QI and MQ show that garnet (a mineral with slow diffusion of oxygen) can armour QI from isotopic exchange with surrounding matrix, even during slow cooling. These differences between δ 18 O in MQ and QI can further be used to test cooling rates by Fast Grain Boundary diffusion modelling. Criteria for identifying QI that preserve primary compositions and are suitable for thermometry were developed based on comparative tests. Relations between δ 18 O and inclusion size, distance of inclusion to host–garnet rim, core–rim zonation of individual inclusions, and presence or absence of petrological features (healed cracks in QI, inclusions in contact with garnet cracks lined by secondary minerals, and secondary minerals along the inclusion grain boundary) were investigated. In this study, 61% of QI preserve primary δ 18 O and 39% were associated with features that were linked to reset δ 18 O values. If δ 18 O in garnet is homogeneous and inclusions are removed, laser-fluorination δ 18 O values of bulk garnet are more precise, more accurate, and best for thermometry. Intragrain δ 18 O(Grt) profiles measured in situ by ion microprobe show no δ 18 O zonation. Almandine–rich garnet (Alm 60–75 ) from each sample was measured by laser-fluorination mass-spectrometry (LF-MS) for δ 18 O and compared with ion microprobe measurements of δ 18 O in QI for thermometry. The Δ 18 OO(Qz–Grt) values for Adirondack samples range from 2.66 to 3.24‰, corresponding to temperatures of 640–740 °C (A[Qz–Alm] = 2.71). Out of 12 samples that were used for thermometry, nine are consistent with previous estimates of peak temperature (625–800 °C) based on petrological and carbon–isotope thermometry for regional granulite and upper amphibolite facies metamorphism. The three samples that disagree with independent thermometry for peak metamorphism are from the anorthosite–mangerite–charnockite–granite suite in the central Adirondacks and yield temperatures of 640–665 °C, ~100 °C lower than previous estimates. Finally, these low temperatures could be interpreted as thermal conditions during late (post-peak) crystallization of garnet on the retrograde path.

58 GEOSCIENCES↗

DOE OSTI2020

Synchrotron x-ray fluorescence analysis reveals diagenetic alteration of fossil melanosome trace metal chemistry

A key feature of the pigment melanin is its high binding affinity for trace metal ions. In modern vertebrates trace metals associated with melanosomes, melanin-rich organelles, can show tissue-specific and taxon-specific distribution patterns. Such signals preserve in fossil melanosomes, informing on the anatomy and phylogenetic affinities of fossil vertebrates. Fossil and modern melanosomes, however, often differ in trace metal chemistry; in particular, melanosomes from fossil vertebrate eyes are depleted in Zn and enriched in Cu relative to their extant counterparts. Whether these chemical differences are biological or taphonomic in origin is unknown, limiting our ability to use melanosome trace metal chemistry to test palaeobiological hypotheses. Here, we use maturation experiments on eye melanosomes from extant vertebrates and synchrotron rapid scan-x-ray fluorescence analysis to show that thermal maturation can dramatically alter melanosome trace element chemistry. In particular, maturation of melanosomes in Cu-rich solutions results in significant depletion of Zn, probably due to low pH and competition effects with Cu. These results confirm fossil melanosome chemistry is susceptible to alteration due to variations in local chemical conditions during diagenesis. Maturation experiments can provide essential data on melanosome chemical taphonomy required for accurate interpretations of preserved chemical signatures in fossils.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

Surface passivation of Fe{sub 3}O{sub 4} nanoparticles with Al{sub 2}O{sub 3} via atomic layer deposition in a rotating fluidized bed reactor

Iron(II,III) oxide (Fe{sub 3}O{sub 4}) nanoparticles have shown great promise in many magnetic-related applications such as magnetic resonance imaging, hyperthermia treatment, and targeted drug delivery. Nevertheless, these nanoparticles are vulnerable to oxidation and magnetization loss under ambient conditions, and passivation is usually required for practical applications. In this work, a home-built rotating fluidized bed (RFB) atomic layer deposition (ALD) reactor was employed to form dense and uniform nanoscale Al{sub 2}O{sub 3} passivation layers on Fe{sub 3}O{sub 4} nanoparticles. The RFB reactor facilitated the precursor diffusion in the particle bed and intensified the dynamic dismantling of soft agglomerates, exposing every surface reactive site to precursor gases. With the aid of in situ mass spectroscopy, it was found that a thicker fluidization bed formed by larger amount of particles increased the residence time of precursors. The prolonged residence time allowed more thorough interactions between the particle surfaces and the precursor gas, resulting in an improvement of the precursor utilization from 78% to nearly 100%, even under a high precursor feeding rate. Uniform passivation layers around the magnetic cores were demonstrated by both transmission electron microscopy and the statistical analysis of Al mass concentrations. Individual particles were coated instead of the soft agglomerates, as was validated by the specific surface area analysis and particle size distribution. The results of thermogravimetric analysis suggested that 5 nm-thick ultrathin Al{sub 2}O{sub 3} coatings could effectively protect the Fe{sub 3}O{sub 4} nanoparticles from oxidation. The x-ray diffraction patterns also showed that the magnetic core crystallinity of such passivated nanoparticles could be well preserved under accelerated oxidation conditions. The precise thickness control via ALD maintained the saturation magnetization at 66.7 emu/g with a 5 nm-thick Al{sub 2}O{sub 3} passivation layer. This good preservation of the magnetic properties with superior oxidation resistance will be beneficial for practical magnetic-based applications.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Collaborative development of diffraction-limited beamline optical systems at US DOE light sources

An ongoing collaboration among four US Department of Energy National Laboratories has demonstrated key technology prototypes and software modeling tools required for new high-coherent flux beamline optical systems. New free electron laser (FEL) and diffraction-limited storage ring (DLSR) light sources demand wavefront preservation from source to sample to achieve and maintain optimal performance. Fine wavefront control was achieved using a novel, roomtemperature cooled mirror system called REAL (resistive element adjustable length) that combines cooling with applied, spatially variable auxiliary heating. Single-grating shearing interferometry (also called Talbot interferometry) and Hartmann wavefront sensors were developed and used for optical characterization and alignment on several beamlines, across a range of photon energies. Demonstrations of non-invasive hard x-ray wavefront sensing were performed using a thin diamond single-crystal as a beamsplitter.

Goldberg, Kenneth A.↗

DOE OSTI2026

DP-TwoLevel: two-stage gradient subspace learning for differentially private federated learning

Federated learning (FL) enables collaborative model training across distributed data sources without sharing raw data, but faces fundamental challenges in communication efficiency and privacy. Differentially private (DP) training mitigates information leakage but introduces noise that degrades model performance, especially in high-dimensional settings. We propose DP-TwoLevel, a hierarchical gradient projection method that improves utility under fixed DP constraints by exploiting low-dimensional structure in model updates. Our approach learns a two-level PCA-based representation of gradients and applies DP noise in a reduced-dimensional subspace, thereby lowering the effective noise magnitude while preserving dominant signal components. We evaluate the method across three datasets (MNIST, Fashion-MNIST, CIFAR-10) and three privacy regimes (ϵ∈0.5, 1.0, 2.0). Across nine experimental settings, DP-TwoLevel consistently outperforms DP-FedAvg, achieving an average accuracy improvement of 9.44%, with larger gains observed in lower ϵ(higher-noise) regimes (up to +22.31%). We further analyze scalability across models ranging from 100K to 1.49M parameters and identify a variance-based success criterion: performance remains strong when the projection preserves more than 75% of gradient variance, degrades in a marginal regime (65–75%), and fails below this threshold. Our results demonstrate that structure-aware dimensionality reduction can significantly improve the privacy–utility tradeoff in FL without modifying formal privacy guarantees. We also provide empirical evidence of scaling limitations for global projections and motivate per-layer extensions for larger models.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗

DOE OSTI2013

Investigation of statistical iterative reconstruction for dedicated breast CT

Purpose: Dedicated breast CT has great potential for improving the detection and diagnosis of breast cancer. Statistical iterative reconstruction (SIR) in dedicated breast CT is a promising alternative to traditional filtered backprojection (FBP). One of the difficulties in using SIR is the presence of free parameters in the algorithm that control the appearance of the resulting image. These parameters require tuning in order to achieve high quality reconstructions. In this study, the authors investigated the penalized maximum likelihood (PML) method with two commonly used types of roughness penalty functions: hyperbolic potential and anisotropic total variation (TV) norm. Reconstructed images were compared with images obtained using standard FBP. Optimal parameters for PML with the hyperbolic prior are reported for the task of detecting microcalcifications embedded in breast tissue.Methods: Computer simulations were used to acquire projections in a half-cone beam geometry. The modeled setup describes a realistic breast CT benchtop system, with an x-ray spectra produced by a point source and an a-Si, CsI:Tl flat-panel detector. A voxelized anthropomorphic breast phantom with 280 μm microcalcification spheres embedded in it was used to model attenuation properties of the uncompressed woman's breast in a pendant position. The reconstruction of 3D images was performed using the separable paraboloidal surrogates algorithm with ordered subsets. Task performance was assessed with the ideal observer detectability index to determine optimal PML parameters.Results: The authors' findings suggest that there is a preferred range of values of the roughness penalty weight and the edge preservation threshold in the penalized objective function with the hyperbolic potential, which resulted in low noise images with high contrast microcalcifications preserved. In terms of numerical observer detectability index, the PML method with optimal parameters yielded substantially improved performance (by a factor of greater than 10) compared to FBP. The hyperbolic prior was also observed to be superior to the TV norm. A few of the best-performing parameter pairs for the PML method also demonstrated superior performance for various radiation doses. In fact, using PML with certain parameter values results in better images, acquired using 2 mGy dose, than FBP-reconstructed images acquired using 6 mGy dose.Conclusions: A range of optimal free parameters for the PML algorithm with hyperbolic and TV norm-based potentials is presented for the microcalcification detection task, in dedicated breast CT. The reported values can be used as starting values of the free parameters, when SIR techniques are used for image reconstruction. Significant improvement in image quality can be achieved by using PML with optimal combination of parameters, as compared to FBP. Importantly, these results suggest improved detection of microcalcifications can be obtained by using PML with lower radiation dose to the patient, than using FBP with higher dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Quantification of organ motion based on an adaptive image-based scale invariant feature method

Purpose: The availability of corresponding landmarks in IGRT image series allows quantifying the inter and intrafractional motion of internal organs. In this study, an approach for the automatic localization of anatomical landmarks is presented, with the aim of describing the nonrigid motion of anatomo-pathological structures in radiotherapy treatments according to local image contrast.Methods: An adaptive scale invariant feature transform (SIFT) was developed from the integration of a standard 3D SIFT approach with a local image-based contrast definition. The robustness and invariance of the proposed method to shape-preserving and deformable transforms were analyzed in a CT phantom study. The application of contrast transforms to the phantom images was also tested, in order to verify the variation of the local adaptive measure in relation to the modification of image contrast. The method was also applied to a lung 4D CT dataset, relying on manual feature identification by an expert user as ground truth. The 3D residual distance between matches obtained in adaptive-SIFT was then computed to verify the internal motion quantification with respect to the expert user. Extracted corresponding features in the lungs were used as regularization landmarks in a multistage deformable image registration (DIR) mapping the inhale vs exhale phase. The residual distances between the warped manual landmarks and their reference position in the inhale phase were evaluated, in order to provide a quantitative indication of the registration performed with the three different point sets.Results: The phantom study confirmed the method invariance and robustness properties to shape-preserving and deformable transforms, showing residual matching errors below the voxel dimension. The adapted SIFT algorithm on the 4D CT dataset provided automated and accurate motion detection of peak to peak breathing motion. The proposed method resulted in reduced residual errors with respect to standard SIFT, providing a motion description comparable to expert manual identification, as confirmed by DIR.Conclusions: The application of the method to a 4D lung CT patient dataset demonstrated adaptive-SIFT potential as an automatic tool to detect landmarks for DIR regularization and internal motion quantification. Future works should include the optimization of the computational cost and the application of the method to other anatomical sites and image modalities.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Prior-based artifact correction (PBAC) in computed tomography

Purpose: Image quality in computed tomography (CT) often suffers from artifacts which may reduce the diagnostic value of the image. In many cases, these artifacts result from missing or corrupt regions in the projection data, e.g., in the case of metal, truncation, and limited angle artifacts. The authors propose a generalized correction method for different kinds of artifacts resulting from missing or corrupt data by making use of available prior knowledge to perform data completion. Methods: The proposed prior-based artifact correction (PBAC) method requires prior knowledge in form of a planning CT of the same patient or in form of a CT scan of a different patient showing the same body region. In both cases, the prior image is registered to the patient image using a deformable transformation. The registered prior is forward projected and data completion of the patient projections is performed using smooth sinogram inpainting. The obtained projection data are used to reconstruct the corrected image. Results: The authors investigate metal and truncation artifacts in patient data sets acquired with a clinical CT and limited angle artifacts in an anthropomorphic head phantom data set acquired with a gantry-based flat detector CT device. In all cases, the corrected images obtained by PBAC are nearly artifact-free. Compared to conventional correction methods, PBAC achieves better artifact suppression while preserving the patient-specific anatomy at the same time. Further, the authors show that prominent anatomical details in the prior image seem to have only minor impact on the correction result. Conclusions: The results show that PBAC has the potential to effectively correct for metal, truncation, and limited angle artifacts if adequate prior data are available. Since the proposed method makes use of a generalized algorithm, PBAC may also be applicable to other artifacts resulting from missing or corrupt data.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Metallic artifact mitigation and organ-constrained tissue assignment for Monte Carlo calculations of permanent implant lung brachytherapy

Purpose: To investigate methods of generating accurate patient-specific computational phantoms for the Monte Carlo calculation of lung brachytherapy patient dose distributions. Methods: Four metallic artifact mitigation methods are applied to six lung brachytherapy patient computed tomography (CT) images: simple threshold replacement (STR) identifies high CT values in the vicinity of the seeds and replaces them with estimated true values; fan beam virtual sinogram replaces artifact-affected values in a virtual sinogram and performs a filtered back-projection to generate a corrected image; 3D median filter replaces voxel values that differ from the median value in a region of interest surrounding the voxel and then applies a second filter to reduce noise; and a combination of fan beam virtual sinogram and STR. Computational phantoms are generated from artifact-corrected and uncorrected images using several tissue assignment schemes: both lung-contour constrained and unconstrained global schemes are considered. Voxel mass densities are assigned based on voxel CT number or using the nominal tissue mass densities. Dose distributions are calculated using the EGSnrc user-code BrachyDose for{sup 125}I, {sup 103}Pd, and {sup 131}Cs seeds and are compared directly as well as through dose volume histograms and dose metrics for target volumes surrounding surgical sutures. Results: Metallic artifact mitigation techniques vary in ability to reduce artifacts while preserving tissue detail. Notably, images corrected with the fan beam virtual sinogram have reduced artifacts but residual artifacts near sources remain requiring additional use of STR; the 3D median filter removes artifacts but simultaneously removes detail in lung and bone. Doses vary considerably between computational phantoms with the largest differences arising from artifact-affected voxels assigned to bone in the vicinity of the seeds. Consequently, when metallic artifact reduction and constrained tissue assignment within lung contours are employed in generated phantoms, this erroneous assignment is reduced, generally resulting in higher doses. Lung-constrained tissue assignment also results in increased doses in regions of interest due to a reduction in the erroneous assignment of adipose to voxels within lung contours. Differences in dose metrics calculated for different computational phantoms are sensitive to radionuclide photon spectra with the largest differences for{sup 103}Pd seeds and smallest but still considerable differences for {sup 131}Cs seeds. Conclusions: Despite producing differences in CT images, dose metrics calculated using the STR, fan beam + STR, and 3D median filter techniques produce similar dose metrics. Results suggest that the accuracy of dose distributions for permanent implant lung brachytherapy is improved by applying lung-constrained tissue assignment schemes to metallic artifact corrected images.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Adaptive nonlocal means filtering based on local noise level for CT denoising

Purpose: To develop and evaluate an image-domain noise reduction method based on a modified nonlocal means (NLM) algorithm that is adaptive to local noise level of CT images and to implement this method in a time frame consistent with clinical workflow. Methods: A computationally efficient technique for local noise estimation directly from CT images was developed. A forward projection, based on a 2D fan-beam approximation, was used to generate the projection data, with a noise model incorporating the effects of the bowtie filter and automatic exposure control. The noise propagation from projection data to images was analytically derived. The analytical noise map was validated using repeated scans of a phantom. A 3D NLM denoising algorithm was modified to adapt its denoising strength locally based on this noise map. The performance of this adaptive NLM filter was evaluated in phantom studies in terms of in-plane and cross-plane high-contrast spatial resolution, noise power spectrum (NPS), subjective low-contrast spatial resolution using the American College of Radiology (ACR) accreditation phantom, and objective low-contrast spatial resolution using a channelized Hotelling model observer (CHO). Graphical processing units (GPU) implementation of this noise map calculation and the adaptive NLM filtering were developed to meet demands of clinical workflow. Adaptive NLM was piloted on lower dose scans in clinical practice. Results: The local noise level estimation matches the noise distribution determined from multiple repetitive scans of a phantom, demonstrated by small variations in the ratio map between the analytical noise map and the one calculated from repeated scans. The phantom studies demonstrated that the adaptive NLM filter can reduce noise substantially without degrading the high-contrast spatial resolution, as illustrated by modulation transfer function and slice sensitivity profile results. The NPS results show that adaptive NLM denoising preserves the shape and peak frequency of the noise power spectrum better than commercial smoothing kernels, and indicate that the spatial resolution at low contrast levels is not significantly degraded. Both the subjective evaluation using the ACR phantom and the objective evaluation on a low-contrast detection task using a CHO model observer demonstrate an improvement on low-contrast performance. The GPU implementation can process and transfer 300 slice images within 5 min. On patient data, the adaptive NLM algorithm provides more effective denoising of CT data throughout a volume than standard NLM, and may allow significant lowering of radiation dose. After a two week pilot study of lower dose CT urography and CT enterography exams, both GI and GU radiology groups elected to proceed with permanent implementation of adaptive NLM in their GI and GU CT practices. Conclusions: This work describes and validates a computationally efficient technique for noise map estimation directly from CT images, and an adaptive NLM filtering based on this noise map, on phantom and patient data. Both the noise map calculation and the adaptive NLM filtering can be performed in times that allow integration with clinical workflow. The adaptive NLM algorithm provides effective denoising of CT data throughout a volume, and may allow significant lowering of radiation dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Monte Carlo study for designing a dedicated “D”-shaped collimator used in the external beam radiotherapy of retinoblastoma patients

Purpose: Retinoblastoma is the most common intraocular malignancy in the early childhood. Patients treated with external beam radiotherapy respond very well to the treatment. However, owing to the genotype of children suffering hereditary retinoblastoma, the risk of secondary radio-induced malignancies is high. The University Hospital of Essen has successfully treated these patients on a daily basis during nearly 30 years using a dedicated “D”-shaped collimator. The use of this collimator that delivers a highly conformed small radiation field, gives very good results in the control of the primary tumor as well as in preserving visual function, while it avoids the devastating side effects of deformation of midface bones. The purpose of the present paper is to propose a modified version of the “D”-shaped collimator that reduces even further the irradiation field with the scope to reduce as well the risk of radio-induced secondary malignancies. Concurrently, the new dedicated “D”-shaped collimator must be easier to build and at the same time produces dose distributions that only differ on the field size with respect to the dose distributions obtained by the current collimator in use. The scope of the former requirement is to facilitate the employment of the authors' irradiation technique both at the authors' and at other hospitals. The fulfillment of the latter allows the authors to continue using the clinical experience gained in more than 30 years. Methods: The Monte Carlo codePENELOPE was used to study the effect that the different structural elements of the dedicated “D”-shaped collimator have on the absorbed dose distribution. To perform this study, the radiation transport through a Varian Clinac 2100 C/D operating at 6 MV was simulated in order to tally phase-space files which were then used as radiation sources to simulate the considered collimators and the subsequent dose distributions. With the knowledge gained in that study, a new, simpler, “D”-shaped collimator is proposed. Results: The proposed collimator delivers a dose distribution which is 2.4 cm wide along the inferior-superior direction of the eyeball. This width is 0.3 cm narrower than that of the dose distribution obtained with the collimator currently in clinical use. The other relevant characteristics of the dose distribution obtained with the new collimator, namely, depth doses at clinically relevant positions, penumbrae width, and shape of the lateral profiles, are statistically compatible with the results obtained for the collimator currently in use. Conclusions: The smaller field size delivered by the proposed collimator still fully covers the planning target volume with at least 95% of the maximum dose at a depth of 2 cm and provides a safety margin of 0.2 cm, so ensuring an adequate treatment while reducing the irradiated volume.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Practical considerations for maximizing heat production in a novel thermobrachytherapy seed prototype

Purpose: A combination of hyperthermia and radiation in the treatment of cancer has been proven to provide better tumor control than radiation administered as a monomodality, without an increase in complications or serious toxicities. Moreover, concurrent administration of hyperthermia and radiation displays synergistic enhancement, resulting in greater tumor cell killing than hyperthermia and radiation delivered separately. The authors have designed a new thermobrachytherapy (TB) seed, which serves as a source of both radiation and heat for concurrent brachytherapy and hyperthermia treatments when implanted in solid tumors. This innovative seed, similar in size and geometry to conventional seeds, will have self-regulating thermal properties. Methods: The new seed's geometry is based on the standard BEST Model 2301{sup 125}I seed, resulting in very similar dosimetric properties. The TB seed generates heat when placed in an oscillating magnetic field via induction heating of a ferromagnetic Ni–Cu alloy core that replaces the tungsten radiographic marker of the standard Model 2301. The alloy composition is selected to undergo a Curie transition near 50 °C, drastically decreasing power production at higher temperatures and providing for temperature self-regulation. Here, the authors present experimental studies of the magnetic properties of Ni–Cu alloy material, the visibility of TB seeds in radiographic imaging, and the ability of seed prototypes to uniformly heat tissue to a desirable temperature. Moreover, analyses are presented of magnetic shielding and thermal expansion of the TB seed, as well as matching of radiation dose to temperature distributions for a short interseed distance in a given treatment volume. Results: Annealing the Ni–Cu alloy has a significant effect on its magnetization properties, increasing the sharpness of the Curie transition. The TB seed preserves the radiographic properties of the BEST 2301 seed in both plain x rays and CT images, and a preliminary experiment demonstrates thermal self-regulation and adequate heating of a tissue-mimicking phantom by seed prototypes. The effect of self-shielding of the seed against the external magnetic field is small, and only minor thermal stress is induced in heating of the seeds from room temperature to well above the seed operating temperature. With proper selection of magnetic field parameters, the thermal dose distribution of an arrangement of TB and hyperthermia-only seeds may be made to match with its radiation dose distribution. Conclusions: The presented analyses address several practical considerations for manufacturing of the proposed TB seeds and identify critical issues for the prototype implementation. The authors’ preliminary experiments demonstrate close agreement with the modeling results, confirming the feasibility of combining sources of heat and radiation into a single thermobrachytherapy seed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Iterative image-domain decomposition for dual-energy CT

Purpose: Dual energy CT (DECT) imaging plays an important role in advanced imaging applications due to its capability of material decomposition. Direct decomposition via matrix inversion suffers from significant degradation of image signal-to-noise ratios, which reduces clinical values of DECT. Existing denoising algorithms achieve suboptimal performance since they suppress image noise either before or after the decomposition and do not fully explore the noise statistical properties of the decomposition process. In this work, the authors propose an iterative image-domain decomposition method for noise suppression in DECT, using the full variance-covariance matrix of the decomposed images. Methods: The proposed algorithm is formulated in the form of least-square estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, the authors include the inverse of the estimated variance-covariance matrix of the decomposed images as the penalty weight in the least-square term. The regularization term enforces the image smoothness by calculating the square sum of neighboring pixel value differences. To retain the boundary sharpness of the decomposed images, the authors detect the edges in the CT images before decomposition. These edge pixels have small weights in the calculation of the regularization term. Distinct from the existing denoising algorithms applied on the images before or after decomposition, the method has an iterative process for noise suppression, with decomposition performed in each iteration. The authors implement the proposed algorithm using a standard conjugate gradient algorithm. The method performance is evaluated using an evaluation phantom (Catphan©600) and an anthropomorphic head phantom. The results are compared with those generated using direct matrix inversion with no noise suppression, a denoising method applied on the decomposed images, and an existing algorithm with similar formulation as the proposed method but with an edge-preserving regularization term. Results: On the Catphan phantom, the method maintains the same spatial resolution on the decomposed images as that of the CT images before decomposition (8 pairs/cm) while significantly reducing their noise standard deviation. Compared to that obtained by the direct matrix inversion, the noise standard deviation in the images decomposed by the proposed algorithm is reduced by over 98%. Without considering the noise correlation properties in the formulation, the denoising scheme degrades the spatial resolution to 6 pairs/cm for the same level of noise suppression. Compared to the edge-preserving algorithm, the method achieves better low-contrast detectability. A quantitative study is performed on the contrast-rod slice of Catphan phantom. The proposed method achieves lower electron density measurement error as compared to that by the direct matrix inversion, and significantly reduces the error variation by over 97%. On the head phantom, the method reduces the noise standard deviation of decomposed images by over 97% without blurring the sinus structures. Conclusions: The authors propose an iterative image-domain decomposition method for DECT. The method combines noise suppression and material decomposition into an iterative process and achieves both goals simultaneously. By exploring the full variance-covariance properties of the decomposed images and utilizing the edge predetection, the proposed algorithm shows superior performance on noise suppression with high image spatial resolution and low-contrast detectability.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Combined iterative reconstruction and image-domain decomposition for dual energy CT using total-variation regularization

Purpose: Dual-energy CT (DECT) is being increasingly used for its capability of material decomposition and energy-selective imaging. A generic problem of DECT, however, is that the decomposition process is unstable in the sense that the relative magnitude of decomposed signals is reduced due to signal cancellation while the image noise is accumulating from the two CT images of independent scans. Direct image decomposition, therefore, leads to severe degradation of signal-to-noise ratio on the resultant images. Existing noise suppression techniques are typically implemented in DECT with the procedures of reconstruction and decomposition performed independently, which do not explore the statistical properties of decomposed images during the reconstruction for noise reduction. In this work, the authors propose an iterative approach that combines the reconstruction and the signal decomposition procedures to minimize the DECT image noise without noticeable loss of resolution. Methods: The proposed algorithm is formulated as an optimization problem, which balances the data fidelity and total variation of decomposed images in one framework, and the decomposition step is carried out iteratively together with reconstruction. The noise in the CT images from the proposed algorithm becomes well correlated even though the noise of the raw projections is independent on the two CT scans. Due to this feature, the proposed algorithm avoids noise accumulation during the decomposition process. The authors evaluate the method performance on noise suppression and spatial resolution using phantom studies and compare the algorithm with conventional denoising approaches as well as combined iterative reconstruction methods with different forms of regularization. Results: On the Catphan©600 phantom, the proposed method outperforms the existing denoising methods on preserving spatial resolution at the same level of noise suppression, i.e., a reduction of noise standard deviation by one order of magnitude. This improvement is mainly attributed to the high noise correlation in the CT images reconstructed by the proposed algorithm. Iterative reconstruction using different regularization, including quadratic orq-generalized Gaussian Markov random field regularization, achieves similar noise suppression from high noise correlation. However, the proposed TV regularization obtains a better edge preserving performance. Studies of electron density measurement also show that our method reduces the average estimation error from 9.5% to 7.1%. On the anthropomorphic head phantom, the proposed method suppresses the noise standard deviation of the decomposed images by a factor of ∼14 without blurring the fine structures in the sinus area. Conclusions: The authors propose a practical method for DECT imaging reconstruction, which combines the image reconstruction and material decomposition into one optimization framework. Compared to the existing approaches, our method achieves a superior performance on DECT imaging with respect to decomposition accuracy, noise reduction, and spatial resolution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2014

Fast dose kernel interpolation using Fourier transform with application to permanent prostate brachytherapy dosimetry

Purpose: Boyer and Mok proposed a fast calculation method employing the Fourier transform (FT), for which calculation time is independent of the number of seeds but seed placement is restricted to calculation grid points. Here an interpolation method is described enabling unrestricted seed placement while preserving the computational efficiency of the original method. Methods: The Iodine-125 seed dose kernel was sampled and selected values were modified to optimize interpolation accuracy for clinically relevant doses. For each seed, the kernel was shifted to the nearest grid point via convolution with a unit impulse, implemented in the Fourier domain. The remaining fractional shift was performed using a piecewise third-order Lagrange filter. Results: Implementation of the interpolation method greatly improved FT-based dose calculation accuracy. The dose distribution was accurate to within 2% beyond 3 mm from each seed. Isodose contours were indistinguishable from explicit TG-43 calculation. Dose-volume metric errors were negligible. Computation time for the FT interpolation method was essentially the same as Boyer's method. Conclusions: A FT interpolation method for permanent prostate brachytherapy TG-43 dose calculation was developed which expands upon Boyer's original method and enables unrestricted seed placement. The proposed method substantially improves the clinically relevant dose accuracy with negligible additional computation cost, preserving the efficiency of the original method.

61 RADIATION PROTECTION AND DOSIMETRY↗