Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Soft computing”

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 541 records · Page 30

Emergent disorder and mechanical memory in periodic metamaterials

Ordered mechanical systems typically have one or only a few stable rest configurations, and hence are not considered useful for encoding memory. Multistable and history-dependent responses usually emerge from quenched disorder, for example in amorphous solids or crumpled sheets. In contrast, due to geometric frustration, periodic magnetic systems can create their own disorder and espouse an extensive manifold of quasi-degenerate configurations. Inspired by the topological structure of frustrated artificial spin ices, we introduce an approach to design ordered, periodic mechanical metamaterials that exhibit an extensive set of spatially disordered states. While our design exploits the correspondence between frustration in magnetism and incompatibility in meta-mechanics, our mechanical systems encompass continuous degrees of freedom, and thus generalize their magnetic counterparts. We show how such systems exhibit non-Abelian and history-dependent responses, as their state can depend on the order in which external manipulations were applied. We demonstrate how this richness of the dynamics enables to recognize, from a static measurement of the final state, the sequence of operations that an extended system underwent. Thus, multistability and potential to perform computation emerge from geometric frustration in ordered mechanical lattices that create their own disorder.

36 MATERIALS SCIENCE↗

Far-forward neutrinos at the Large Hadron Collider

We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of charged-current tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-kT effects are approximated by a simple phenomenological model introducing a Gaussian kT-smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for D± s production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Anatomically Validated Brachial Plexus Contouring Method for Intensity Modulated Radiation Therapy Planning

Purpose: To develop contouring guidelines for the brachial plexus (BP) using anatomically validated cadaver datasets. Magnetic resonance imaging (MRI) and computed tomography (CT) were used to obtain detailed visualizations of the BP region, with the goal of achieving maximal inclusion of the actual BP in a small contoured volume while also accommodating for anatomic variations. Methods and Materials: CT and MRI were obtained for 8 cadavers positioned for intensity modulated radiation therapy. 3-dimensional reconstructions of soft tissue (from MRI) and bone (from CT) were combined to create 8 separate enhanced CT project files. Dissection of the corresponding cadavers anatomically validated the reconstructions created. Seven enhanced CT project files were then automatically fitted, separately in different regions, to obtain a single dataset of superimposed BP regions that incorporated anatomic variations. From this dataset, improved BP contouring guidelines were developed. These guidelines were then applied to the 7 original CT project files and also to 1 additional file, left out from the superimposing procedure. The percentage of BP inclusion was compared with the published guidelines. Results: The anatomic validation procedure showed a high level of conformity for the BP regions examined between the 3-dimensional reconstructions generated and the dissected counterparts. Accurate and detailed BP contouring guidelines were developed, which provided corresponding guidance for each level in a clinical dataset. An average margin of 4.7 mm around the anatomically validated BP contour is sufficient to accommodate for anatomic variations. Using the new guidelines, 100% inclusion of the BP was achieved, compared with a mean inclusion of 37.75% when published guidelines were applied. Conclusion: Improved guidelines for BP delineation were developed using combined MRI and CT imaging with validation by anatomic dissection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Fractional modeling of viscoelasticity in 3D cerebral arteries and aneurysms

We develop efficient numerical methods for fractional order PDEs, and employ them to investigate viscoelastic constitutive laws for arterial wall mechanics. Recent simulations using one-dimensional models [1] have indicated that fractional order models may offer a more powerful alternative for modeling the arterial wall response, exhibiting reduced sensitivity to parametric uncertainties compared with the integer-calculus-based models. Here, we study three-dimensional (3D) fractional PDEs that naturally model the continuous relaxation properties of soft tissue, and for the first time employ them to simulate flow structure interactions for patient-specific brain aneurysms. To deal with the high memory requirements and in order to accelerate the numerical evaluation of hereditary integrals, we employ a fast convolution method [2] that reduces the memory cost to O(log(N)) and the computational complexity to O(N log(N)). Furthermore, we combine the fast convolution with high-order backward differentiation to achieve third-order time integration accuracy. We confirm that in 3D viscoelastic simulations, the integer order models strongly depends on the relaxation parameters, while the fractional order models are less sensitive. As an application to long-time simulations in complex geometries, we also apply the method to modeling fluid-structure interaction of a 3D patient-specific compliant cerebral artery with an aneurysm. Taken together, our findings demonstrate that fractional calculus can be employed effectively in modeling complex behavior of materials in realistic 3D time-dependent problems if properly designed efficient algorithms are employed to overcome the extra memory requirements and computational complexity associated with the non-local character of fractional derivatives.

Computer Science↗

Bioinspired design of flexible armor based on chiton scales

Man-made armors often rely on rigid structures for mechanical protection, which typically results in a trade-off with flexibility and maneuverability. Chitons, a group of marine mollusks, evolved scaled armors that address similar challenges. Many chiton species possess hundreds of small, mineralized scales arrayed on the soft girdle that surrounds their overlapping shell plates. Ensuring both flexibility for locomotion and protection of the underlying soft body, the scaled girdle is an excellent model for multifunctional armor design. Here we conduct a systematic study of the material composition, nanomechanical properties, three-dimensional geometry, and interspecific structural diversity of chiton girdle scales. Moreover, inspired by the tessellated organization of chiton scales, we fabricate a synthetic flexible scaled armor analogue using parametric computational modeling and multi-material 3D printing. This approach allows us to conduct a quantitative evaluation of our chiton-inspired armor to assess its orientation-dependent flexibility and protection capabilities.

Science & Technology - Other Topics↗

Metal artifact correction for x-ray computed tomography using kV and selective MV imaging

Purpose: The overall goal of this work is to improve the computed tomography (CT) image quality for patients with metal implants or fillings by completing the missing kilovoltage (kV) projection data with selectively acquired megavoltage (MV) data that do not suffer from photon starvation. When both of these imaging systems, which are available on current radiotherapy devices, are used, metal streak artifacts are avoided, and the soft-tissue contrast is restored, even for regions in which the kV data cannot contribute any information. Methods: Three image-reconstruction methods, including two filtered back-projection (FBP)-based analytic methods and one iterative method, for combining kV and MV projection data from the two on-board imaging systems of a radiotherapy device are presented in this work. The analytic reconstruction methods modify the MV data based on the information in the projection or image domains and then patch the data onto the kV projections for a FBP reconstruction. In the iterative reconstruction, the authors used dual-energy (DE) penalized weighted least-squares (PWLS) methods to simultaneously combine the kV/MV data and perform the reconstruction. Results: The authors compared kV/MV reconstructions to kV-only reconstructions using a dental phantom with fillings and a hip-implant numerical phantom. Simulation results indicated that dual-energy sinogram patch FBP and the modified dual-energy PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in the kV projections. The root-mean-square errors of soft-tissue patterns obtained using combined kV/MV data are 10–15 Hounsfield units smaller than those of the kV-only images, and the structural similarity index measure also indicates a 5%–10% improvement in the image quality. The added dose from the MV scan is much less than the dose from the kV scan if a high efficiency MV detector is assumed. Conclusions: The authors have shown that it is possible to improve the image quality of kV CTs for patients with metal implants or fillings by completing the missing kV projection data with selectively acquired MV data that do not suffer from photon starvation. Numerical simulations demonstrated that dual-energy sinogram patch FBP and a modified kV/MV PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in kV projections. Combined kV/MV images may permit the improved delineation of structures of interest in CT images for patients with metal implants or fillings.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Mobility Selective Ion Soft-Landing and Characterization Enabled Using Structures for Lossless Ion Manipulation

In this study, while conventional ion-soft landing uses the mass-to-charge (m/z) ratio to achieve molecular selection for deposition, here we demonstrate the use of Structures for Lossless Ion Manipulation (SLIM) for mobility-based ion selection and deposition. The dynamic rerouting capabilities of SLIM were leveraged to enable the rerouting of a selected range of mobilities to a different SLIM path (rather than MS) that terminated at a deposition surface. A selected mobility range from a phosphazene ion mixture was rerouted and deposited with a current pulse (~150 pA) resembling its mobility peak. In addition, from a mixture of tetra-alkyl ammonium (TAA) ions containing chain lengths of C5–C8, selected chains (C6, C7) were collected on a surface, reconstituted into solution-phase, and subsequently analyzed with a SLIM-qToF to obtain an IMS/MS spectrum, confirming the identity of the selected species. Further, this method was used to characterize triply charged tungsten-polyoxometalate anions, PW 12 O 40 3– (WPOM). The arrival time distribution of the IMS/MS showed multiple peaks associated with the triply charged anion (PW 12 O 40 3– ), of which a selected ATD was deposited and imaged using TEM. Additionally, the identity of the deposited WPOM was ascertained using energy-dispersive (EDS) spectroscopy. Further, we present theory and computations that reveal ion landing energies, the ability to modulate the energies, and deposition spot sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High energy density soft X-ray momentum coupling to comet analogs for NEO mitigation

Here, we applied MBBAY high fluence pulsed radiation intensity driven momentum transfer analysis to calculate X-ray momentum coupling coefficients C M =(Pa s)/(J/m 2 ) for two simplified comet analog materials: i) water ice, and ii) 70% water ice and 30% distributed olivine grains. The momentum coupling coefficients (C M ) max of 50×10 –5 s/m, are about an order of magnitude greater than experimentally determined and computed MBBAY values for meteoritic materials that are analogs for asteroids. From the values for comet analog materials we infer applied energies (via momentum transfer) required to deflect an Earth crossing comet from impacting Earth by a sufficient amount (~1 cm/s) to avert collision ~a year in advance. Comet model calculations indicate for C M = 5 × 10 –4 s/m the deflection of a 2 km comet with a density 600 kg/m 3 by 1 cm/s requires an applied energy on the target surface of 5 × 10 13 J, the equivalent of 12 kT of TNT. Depending on the geometrical configuration of the interaction the explosive yield required could be an order of magnitude higher.

79 ASTRONOMY AND ASTROPHYSICS↗

Inferring colloidal interaction from scattering by machine learning

A machine learning solution for the potential inversion problem in elastic scattering is outlined. The inversion scheme consists of two major components, a generative network featuring a variational autoencoder which extracts the targeted static two-point correlation functions from experimentally measured scattering cross sections, and a Gaussian process framework which probabilistically infers the relevant structural parameters from the inverted correlation functions. Via a case study of charged colloidal suspensions, the feasibility of this approach for quantitative study of molecular interaction is critically benchmarked and its merit over existing deterministic approaches, in terms of numerical accuracy and computationally efficiency, is demonstrated.

36 MATERIALS SCIENCE↗

Operator origin of anomalous dimensions in de Sitter space

The late-time limit of the power spectrum for heavy (principal series) fields in de Sitter (dS) space yields a series of polynomial terms with complex scaling dimensions. Such scaling behavior is expected to result from an associated operator with a complex dimension. In a free theory, these complex dimensions are known to match the constraints imposed by unitarity on the space of states. Yet, perturbative corrections to the scaling behavior of operators are naively inconsistent with unitary evolution of the quantum fields in dS space. This paper demonstrates how to compute one-loop corrections to the scaling dimensions that appear in the two-point function from the field theory description in terms of local operators. We first show how to evaluate these anomalous dimensions using Mellin space, which has the feature that it naturally accommodates a scaleless regulator. We then explore the consequences for the soft de Sitter effective theory (SdSET) description that emerges in the long wavelength limit. Carefully matching between the UV and SdSET descriptions requires the introduction of novel nondynamical “operators” in the effective theory. This is not only necessary to reproduce results extracted from the Källén-Lehmann representation (that use the space of unitary states directly), but it is also required by general arguments that invoke positivity. Published by the American Physical Society 2025

Cohen, Timothy↗

Back-to-Back Inclusive Dijets in Deep Inelastic Scattering at Small x : Complete NLO Results and Predictions

We compute the back-to-back dijet cross section in deep inelastic scattering at small x to next-to-leading order (NLO) in the color glass condensate effective field theory. Our result can be factorized into a convolution of the Weizsäcker-Williams gluon transverse-momentum-dependent distribution function (WW gluon TMD) with a universal soft factor and an NLO coefficient function. The soft factor includes both double and single logarithms in the ratio of the relative transverse momentum P ⊥ of the dijet pair to the dijet momentum imbalance q ⊥ ; its renormalization group (RG) evolution is resummed into the Sudakov factor. Likewise, the WW TMD obeys a nonlinear RG equation in x that is kinematically constrained to satisfy both the lifetime and rapidity ordering of the projectile. Exact analytical expressions are obtained for the NLO coefficient function of transversely and longitudinally polarized photons. Our results allow for the first quantitative separation of the dynamics of Sudakov suppression from that of gluon saturation. They can be extended to other final states and provide a framework for precision tests of novel QCD many-body dynamics at the Electron-Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-G-JeP3-03: Effect of Robot Pose On Beam Blocking for Ultrasound Guided SBRT of the Prostate

Purpose: Ultrasound presents a fast, volumetric image modality for real-time tracking of abdominal organ motion. How-ever, ultrasound transducer placement during radiation therapy is challenging. Recently, approaches using robotic arms for intra-treatment ultrasound imaging have been proposed. Good and reliable imaging requires placing the transducer close to the PTV. We studied the effect of a seven degrees of freedom robot on the fea-sible beam directions. Methods: For five CyberKnife prostate treatment plans we established viewports for the transducer, i.e., points on the patient surface with a soft tissue view towards the PTV. Choosing a feasible transducer pose and using the kinematic redundancy of the KUKA LBR iiwa robot, we considered three robot poses. Poses 1 to 3 had the elbow point anterior, superior, and inferior, respectively. For each pose and each beam starting point, the pro-jections of robot and PTV were computed. We added a 20 mm margin accounting for organ / beam motion. The number of nodes for which the PTV was partially of fully blocked were established. Moreover, the cumula-tive overlap for each of the poses and the minimum overlap over all poses were computed. Results: The fully and partially blocked nodes ranged from 12% to 20% and 13% to 27%, respectively. Typically, pose 3 caused the fewest blocked nodes. The cumulative overlap ranged from 19% to 29%. Taking the minimum overlap, i.e., considering moving the robot’s elbow while maintaining the transducer pose, the cumulative over-lap was reduced to 16% to 18% and was 3% to 6% lower than for the best individual pose. Conclusion: Our results indicate that it is possible to identify feasible ultrasound transducer poses and to use the kinematic redundancy of a 7 DOF robot to minimize the impact of the imaging subsystem on the feasible beam directions for ultrasound guided and motion compensated SBRT. Research partially funded by DFG grants ER 817/1-1 and SCHL 1844/3-1.

60 APPLIED LIFE SCIENCES↗

Far-forward neutrinos at the Large Hadron Collider

We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of charged-current tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-k T effects are approximated by a simple phenomenological model introducing a Gaussian k T -smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for $D^{±}_s$ production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Table-top two-color soft X-ray laser by means of Ni-like plasmas

Laser-produced Ni-like plasmas are known as active media for extreme ultraviolet lasing, with the flexibility to two-color lasing. Two-color laser generation is very complex at accelerator facilities. In this work, plasma lasing at the 3d{sup 9}4d{sup 1}(J = 0) → 3d{sup 9}4p{sup 1}(J = 1) (collisional-pumping process) and the 3d{sup 9}4f{sup 1}(J = 1) → 3d{sup 9}4d{sup 1}(J = 1) (photo-pumping process) transitions is studied experimentally and computationally. Several key characteristics of collisional- and photo-pumping laser, such as divergence, pointing stability, and intensity have been investigated. The measurements showed different pulse characteristics for the two lasing processes affected by plasma inhomogeneity in temperature and density. Analytical expressions of these characteristics for both collisional- and photo-pumping are derived. It is found that the plasma that maximizes the photo-pumping lasing is 20% hotter and 70% denser than the plasma that optimizes the collisional-pumping lasing. The gain of collisional pumping is ≈4 times higher than the gain for the photo-pumping. The gain lifetime is a factor of ≈5.2 larger for the monopole-pumping. Similarly, the gain thickness is a factor of ≈1.8 larger. It is also found that the gain build-up time for collisional- and photo-pumping is 0.7 ps and 0.9 ps, respectively, whereas the build-up length-scale is 11.5 μm and 6.3 μm, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

COOLING RATES FOR RELATIVISTIC ELECTRONS UNDERGOING COMPTON SCATTERING IN STRONG MAGNETIC FIELDS

For inner magnetospheric models of hard X-ray and gamma-ray emission in high-field pulsars and magnetars, resonant Compton upscattering is anticipated to be the most efficient process for generating continuum radiation. This is in part due to the proximity of a hot soft photon bath from the stellar surface to putative radiation dissipation regions in the inner magnetosphere. Moreover, because the scattering process becomes resonant at the cyclotron frequency, the effective cross section exceeds the classical Thomson value by over two orders of magnitude, thereby enhancing the efficiency of continuum production and the cooling of relativistic electrons. This paper presents computations of the electron cooling rates for this process, which are needed for resonant Compton models of non-thermal radiation from such highly magnetized pulsars. The computed rates extend previous calculations of magnetic Thomson cooling to the domain of relativistic quantum effects, sampled near and above the quantum critical magnetic field of 44.13 TG. This is the first exposition of fully relativistic, quantum magnetic Compton cooling rates for electrons, and it employs both the traditional Johnson and Lippmann cross section and a newer Sokolov and Ternov (ST) formulation of Compton scattering in strong magnetic fields. Such ST formalism is formally correct for treating spin-dependent effects that are important in the cyclotron resonance and has not been addressed before in the context of cooling by Compton scattering. The QED effects are observed to profoundly lower the rates below extrapolations of the familiar magnetic Thomson results, as expected, when recoil and Klein-Nishina reductions become important.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FORMATION OF THE ECCENTRIC-ORBIT MILLISECOND PULSAR J1903+0327 AND THE ORIGIN OF SINGLE MILLISECOND PULSARS

The millisecond pulsar (MSP) J1903+0327 is accompanied by an ordinary G dwarf star in an unusually wide (P{sub orb} {approx_equal} 95.2 days) and eccentric (e {approx_equal} 0.44) orbit. The standard model for producing MSPs fails to explain the orbital characteristics of this extraordinary binary, and alternative binary models are unable to explain the observables. We present a triple-star model for producing MSPs in relatively wide eccentric binaries with a normal (main-sequence) stellar companion. We start from a stable triple system consisting of a low-mass X-ray binary (LMXB) with an orbital period of at least 1 day, accompanied by a G dwarf in a wide and possibly eccentric orbit. Variations in the initial conditions naturally provide a satisfactory explanation for the unexplained triple component in the eclipsing soft X-ray transient 4U 2129+47 or the cataclysmic variable EC 19314-5915. The best explanation for J1903+0327, however, results from the expansion of the orbit of the LMXB, driven by the mass transfer from the evolving donor star to its neutron star companion, which causes the triple eventually to become dynamically unstable. Using numerical computations we show that, depending on the precise system configuration at the moment the triple becomes dynamically unstable, the ejection of each of the three components is possible. If the donor star of the LMXB is ejected, a system resembling J1903+0327 will result. If the neutron star is ejected, a single MSP results. This model therefore also provides a straightforward mechanism for forming a single MSP in the Galactic disk. We conclude that the Galaxy contains some 30-300 binaries with characteristics similar to J1903+0327 and about an order of magnitude fewer single MSPs produced with the proposed triple scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

An enhanced SOCP-based method for feeder load balancing using the multi-terminal soft open point in active distribution networks

The integration of distributed generators (DGs) exacerbates the feeder power flow fluctuation and load unbalanced condition in active distribution networks (ADNs). The unbalanced feeder load causes inefficient use of network assets and network congestion during system operation. The flexible interconnection based on the multi-terminal soft open point (SOP) significantly benefits the operation of ADNs. The multi-terminal SOP, which is a controllable power electronic device installed to replace the normally open point, provides accurate active and reactive power flow control to enable the flexible connection of feeders. An enhanced SOCP-based method for feeder load balancing using the multi-terminal SOP is proposed in this paper. Furthermore, by regulating the operation of the multi-terminal SOP, the proposed method can mitigate the unbalanced condition of feeder load and simultaneously reduce the power losses of ADNs. Then, the original non-convex model is converted into a second-order cone programming (SOCP) model using convex relaxation. In order to tighten the SOCP relaxation and improve the computation efficiency, an enhanced SOCP-based approach is developed to solve the proposed model. Finally, case studies are performed on the modified IEEE 33-node system to verify the effectiveness and efficiency of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sparse-view proton computed tomography using modulated proton beams

Purpose: Proton imaging that uses a modulated proton beam and an intensity detector allows a relatively fast image acquisition compared to the imaging approach based on a trajectory tracking detector. In addition, it requires a relatively simple implementation in a conventional proton therapy equipment. The model of geometric straight ray assumed in conventional computed tomography (CT) image reconstruction is however challenged by multiple-Coulomb scattering and energy straggling in the proton imaging. Radiation dose to the patient is another important issue that has to be taken care of for practical applications. In this work, the authors have investigated iterative image reconstructions after a deconvolution of the sparsely view-sampled data to address these issues in proton CT. Methods: Proton projection images were acquired using the modulated proton beams and the EBT2 film as an intensity detector. Four electron-density cylinders representing normal soft tissues and bone were used as imaged object and scanned at 40 views that are equally separated over 360°. Digitized film images were converted to water-equivalent thickness by use of an empirically derived conversion curve. For improving the image quality, a deconvolution-based image deblurring with an empirically acquired point spread function was employed. They have implemented iterative image reconstruction algorithms such as adaptive steepest descent-projection onto convex sets (ASD-POCS), superiorization method–projection onto convex sets (SM-POCS), superiorization method–expectation maximization (SM-EM), and expectation maximization-total variation minimization (EM-TV). Performance of the four image reconstruction algorithms was analyzed and compared quantitatively via contrast-to-noise ratio (CNR) and root-mean-square-error (RMSE). Results: Objects of higher electron density have been reconstructed more accurately than those of lower density objects. The bone, for example, has been reconstructed within 1% error. EM-based algorithms produced an increased image noise and RMSE as the iteration reaches about 20, while the POCS-based algorithms showed a monotonic convergence with iterations. The ASD-POCS algorithm outperformed the others in terms of CNR, RMSE, and the accuracy of the reconstructed relative stopping power in the region of lung and soft tissues. Conclusions: The four iterative algorithms, i.e., ASD-POCS, SM-POCS, SM-EM, and EM-TV, have been developed and applied for proton CT image reconstruction. Although it still seems that the images need to be improved for practical applications to the treatment planning, proton CT imaging by use of the modulated beams in sparse-view sampling has demonstrated its feasibility.

60 APPLIED LIFE SCIENCES↗