SEARCH · Search NASA
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.
We present the pressure-induced polymerization of 1,4-cyclohexadiene, wherein thermal methods may access multiple reaction pathways. In contrast, light enables solely cyclobutane-based nanothreads to arise from a single viable reaction pathway.
We present a reduced-scaling auxiliary-field quantum Monte Carlo (AFQMC) framework designed for large molecular systems and ensembles, with or without coupling to optical cavities. Our approach leverages the natural block sparsity of the Cholesky decomposition (CD) of electron repulsion integrals in molecular ensembles and employs tensor hypercontraction (THC) to efficiently compress low-rank Cholesky blocks. By representing the Cholesky vectors in a mixed format, keeping high-rank blocks in block-sparse form and compressing low-rank blocks with THC, we reduce the scaling of exchange-energy evaluation from quartic to robust cubic in the number of molecular orbitals N, while lowering memory from cubic toward quadratic. Benchmark analyses on one-, two-, and three-dimensional molecular ensembles (up to ∼1,200 orbitals) show that (a) the number of nonzeros in Cholesky tensors grows linearly with system size across dimensions; (b) the average numerical rank increases sublinearly and does not saturate at these sizes; and (c) rank heterogeneity─some blocks nearly full rank and many low rank, naturally motivates the proposed mixed block sparsity and THC scheme for efficient calculation of exchange energy. In conclusion, we demonstrate that the mixed scheme yields cubic wall-time scaling with favorable prefactors and preserves AFQMC accuracy.
Direct observation of functional motions in protein structures is highly desirable for understanding how these nanomachineries of life operate at the molecular level. Because cryogenic temperatures are non-physiological and may prohibit or even alter protein structural dynamics, it is necessary to develop robust X-ray diffraction methods that enable routine data collection at room temperature. We recently reported a crystal-on-crystal device to facilitate in situ diffraction of protein crystals at room temperature devoid of any sample manipulation. Here an automated serial crystallography platform based on this crystal-on-crystal technology is presented. A hardware and software prototype has been implemented, and protocols have been established that allow users to image, recognize and rank hundreds to thousands of protein crystals grown on a chip in optical scanning mode prior to serial introduction of these crystals to an X-ray beam in a programmable and high-throughput manner. This platform has been tested extensively using fragile protein crystals. We demonstrate that with affordable sample consumption, this in situ serial crystallography technology could give rise to room-temperature protein structures of higher resolution and superior map quality for those protein crystals that encounter difficulties during freezing. This serial data collection platform is compatible with both monochromatic oscillation and Laue methods for X-ray diffraction and presents a widely applicable approach for static and dynamic crystallographic studies at room temperature.
Purpose: Determine the impact of morphologic characteristics (e.g. necrosis, vascular enhancement, and cavitation) on radiomic features from contrast enhanced CT (CE-CT) in primary lung tumors. Methods: We developed an auto-segmentation algorithm to separate lung tumors on contrast-enhanced CT into cavitation (air), necrosis, tissue, and enhancing vessels using a combination of thresholding and region-growing. An auto-segmentation algorithm was also designed to identify necrosis on FDG-PET scans. Wilcoxon rank-sum tests were used to determine if significant differences existed in radiomics features (histogram-uniformity and Laplacian-of-Gaussian average) from 249 patients, found to prognostic in previous work, based on the presence/absence of morphologic features. Feature values were also compared between the original tumor contours and contours excluding a specific morphologic feature. Comparison of necrosis segmentation on CE-CT versus FDG-PET was performed in 78 patients to assess for agreement using the concordance correlation coefficient (CCC). Results: Tumors with cavitation and enhancing vasculature had lower uniformity values (p = 0.001 and p = 0.03, respectively). Tumors with enhancing vasculature and necrosis had higher Laplacian-of-Gaussian average values (measure of “edges” within the tumor) (p < 0.001). Removing these tissue types from regions-of-interest did not drastically alter either radiomic feature value (all scenarios had R{sup 2} > 0.8). This suggests there may be interactions between morphologic characteristics and the radiomic feature value of tumor tissue. Comparison of necrosis volume and percent necrosis volume of tumor were found to have CCC values of 0.85 and 0.76, respectively between CE-CT and FDG-PET segmentation methods. Conclusions: Tumors with enhancing vasculature, necrosis, and cavitation have higher radiomic feature values that are associated with poor prognosis than tumors without these features. Removing these tissue types from quantitative assessment did not drastically impact radiomic feature values. High reproducibility of CE-CT segmented necrosis compared to FDG-PET segmented necrosis provides a reasonable validation of segmentation accuracy on CE-CT.
Purpose: We previously developed a multigene expression model of tumor radiation sensitivity index (RSI) with clinical validation in multiple independent cohorts (breast, rectal, esophageal, and head and neck patients). The purpose of this study was to assess differences between RSI scores in primary colon cancer and metastases. Methods and Materials: Patients were identified from our institutional review board–approved prospective observational protocol. A total of 704 metastatic and 1362 primary lesions were obtained from a de-identified metadata pool. RSI was calculated using the previously published rank-based algorithm. An independent cohort of 29 lung or liver colon metastases treated with 60 Gy in 5 fractions stereotactic body radiation therapy (SBRT) was used for validation. Results: The most common sites of metastases included liver (n=374; 53%), lung (n=116; 17%), and lymph nodes (n=40; 6%). Sixty percent of metastatic tumors, compared with 54% of primaries, were in the RSI radiation-resistant peak, suggesting metastatic tumors may be slightly more radiation resistant than primaries (P=.01). In contrast, when we analyzed metastases based on anatomical site, we uncovered large differences in RSI. The median RSIs for metastases in descending order of radiation resistance were ovary (0.48), abdomen (0.47), liver (0.43), brain (0.42), lung (0.32), and lymph nodes (0.31) (P<.0001). These findings were confirmed when the analysis was restricted to lesions from the same patient (n=139). In our independent cohort of treated lung and liver metastases, lung metastases had an improved local control rate compared to that in patients with liver metastases (2-year local control rate of 100% vs 73.0%, respectively; P=.026). Conclusions: Assessment of radiation sensitivity between primary and metastatic tissues of colon cancer histology revealed significant differences based on anatomical location of metastases. These initial results warrant validation in a larger clinical cohort.
Purpose: To determine whether previously identified quantitative image features (QIFs) based on {sup 18}F-fluorodeoxyglucose positron emission tomography (FDG-PET) (co-occurrence matrix energy and solidity) are able to isolate subgroups of patients who would receive a benefit or detriment from dose escalation in terms of overall survival (OS) or progression-free survival (PFS). Methods and Materials: Subgroups of a previously analyzed 225 patient cohort were generated with the use of 5-percentile increment cutoff values of disease solidity and primary tumor co-occurrence matrix energy. The subgroups were analyzed with a log-rank test to determine whether there was a difference in OS and PFS between patients treated with 60 to 70 Gy and those receiving 74 Gy. Results: In the entire patient cohort, there was no statistical difference in terms of OS or PFS between patients receiving 74 Gy and those receiving 60 to 70 Gy. It was qualitatively observed that as disease solidity and primary co-occurrence matrix energy increased, patients receiving 74 Gy had an improved OS and PFS compared with those receiving 60 to 70 Gy. The opposite trend (detriment of receiving 74 Gy) was also observed regarding low values of disease solidity and primary co-occurrence matrix energy. Conclusions: FDG-PET–based QIFs were found to be capable of isolating subgroups of patients who received a benefit or detriment from dose escalation.
We prove, using a method based on map-state duality, lower bounds on entanglement resources needed to deterministically implement a bipartite unitary using separable (SEP) operations, which include LOCC (local operations and classical communication) as a particular case. It is known that the Schmidt rank of an entangled pure state resource cannot be less than the Schmidt rank of the unitary. We prove that if these ranks are equal the resource must be uniformly (maximally) entangled: equal nonzero Schmidt coefficients. Higher rank resources can have less entanglement: we have found numerical examples of Schmidt rank 2 unitaries, which can be deterministically implemented, by either SEP or LOCC, using an entangled resource of two qutrits with less than one ebit of entanglement.
As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed employing the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.
As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed employing the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.
As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed emplying the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.
Purpose: Previous studies have demonstrated how imaging of the breast with patients lying prone using a supportive positioning device markedly facilitates longitudinal and/or multimodal image registration. In this contribution, the authors’ primary objective was to determine if there are differences in the standardized uptake value (SUV) derived from [{sup 18}F]fluorodeoxyglucose (18F-FDG) positron emission tomography (PET) in breast tumors imaged in the standard supine position and in the prone position using a specialized positioning device. Methods: A custom positioning device was constructed to allow for breast scanning in the prone position. Rigid and nonrigid phantom studies evaluated differences in prone and supine PET. Clinical studies comprised 18F-FDG-PET of 34 patients with locally advanced breast cancer imaged in the prone position (with the custom support) followed by imaging in the supine position (without the support). Mean and maximum values (SUV{sub peak} and SUV{sub max}, respectively) were obtained from tumor regions-of-interest for both positions. Prone and supine SUV were linearly corrected to account for the differences in 18F-FDG uptake time. Correlation, Bland–Altman, and nonparametric analyses were performed on uptake time-corrected and uncorrected data. Results: SUV from the rigid PET breast phantom imaged in the prone position with the support device was 1.9% lower than without the support device. In the nonrigid PET breast phantom, prone SUV with the support device was 5.0% lower than supine SUV without the support device. In patients, the median (range) difference in uptake time between prone and supine scans was 16.4 min (13.4–30.9 min), which was significantly—but not completely—reduced by the linear correction method. SUV{sub peak} and SUV{sub max} from prone versus supine scans were highly correlated, with concordance correlation coefficients of 0.91 and 0.90, respectively. Prone SUV{sub peak} and SUV{sub max} were significantly lower than supine in both original and uptake time-adjusted data across a range of index times (P < < 0.0001, Wilcoxon signed rank test). Before correcting for uptake time differences, Bland–Altman analyses revealed proportional bias between prone and supine measurements (SUV{sub peak} and SUV{sub max}) that increased with higher levels of FDG uptake. After uptake time correction, this bias was significantly reduced (P < 0.01). Significant prone-supine differences, with regard to the spatial distribution of lesions relative to isocenter, were observed between the two scan positions, but this was poorly correlated with the residual (uptake time-corrected) prone-supine SUV{sub peak} difference (P = 0.78). Conclusions: Quantitative 18F-FDG-PET/CT of the breast in the prone position is not deleteriously affected by the support device but yields SUV that is consistently lower than those obtained in the standard supine position. SUV differences between scans arising from FDG uptake time differences can be substantially reduced, but not removed entirely, with the current correction method. SUV from the two scan orientations is quantitatively different and should not be assumed equivalent or interchangeable within the same subject. These findings have clinical relevance in that they underscore the importance of patient positioning while scanning as a clinical variable that must be accounted for with longitudinal PET measurement, for example, in the assessment of treatment response.
Purpose: To investigate the feasibility of using DIR to propagate the manually contoured rectum and bladder from the 1st insertion to the new CT images on subsequent insertions and evaluate the segmentation performance. Methods: Ten cervical cancer patients, who were treated by T&R brachytherapy in 3–4 insertions, were retrospectively collected. In each insertion, rectum and bladder were manually delineated on the planning CT by a physicist and verified by a radiation oncologist. Using VelocityAI (Velocity Medical Solutions, Atlanta, GA), a rigid registration was firstly employed to match the bony structures between the first insertion and each of the following insertions, then a multi-pass B-spine DIR was carried out to further map the sub volume that encompasses rectum and bladder. The resultant deformation fields propagated contours, and dice similarity coefficient (DSC) was used to quantitatively evaluate the agreement between the propagated contours and the manually-delineated organs. For the 3rd insertion, we also evaluated if the segmentation performance could be improved by propagating the contours from the most recent insertion, i.e., the 2nd insertion. Results: On average, the contour propagation took about 1 minute. The average and standard deviation of DSC over all insertions and patients was 0.67±0.10 (range: 0.44–0.81) for rectum, and 0.78±0.07 (range: 0.63–0.87) for bladder. For the 3rd insertion, propagating contours from the 2nd insertion could improve the segmentation performance in terms of DSC from 0.63±0.10 to 0.72±0.08 for rectum, and from 0.77±0.07 to 0.79±0.06 for bladder. A Wilcoxon signed rank test indicated that the improvement was statistically significant for rectum (p = 0.004). Conclusion: The preliminary results demonstrate that deformable image registration could efficiently and accurately propagate rectum and bladder contours between CT images in different T&R brachytherapy fractions. We are incorporating the propagated contours into our learning-based method to further segment these organs.
Multiphase flow can be difficult to simulate with high accuracy due to the wide range of scales associated with various multiphase phenomena. These scales may range from the size of individual particles to the entire domain of interest. Traditionally, large scale systems can only be simulated using averaging approaches that filter out the locations of individual particles. In this work, the Euler–Lagrange method is used to simulate large-scale dense particle systems in which each individual particle is tracked. In order to accomplish this, the highly scalable spectral element code nek5000 has been extended to handle the multiple levels of multiphase coupling in these systems. These levels include what has been called one-, two-, and four-way coupling. Here, each level has been separated to detail the computational impact of each stage. A binned ghost particle algorithm has also been developed to efficiently handle the challenges of two- and four-way coupling in a parallel processing context. The algorithms and their implementations are then shown to scale to 65,536 Message Passing Interface (MPI) ranks in both the strong and weak limits. After this, validation is performed through simulation of a small-scale fluidized bed. Lastly, a large-scale fluidized bed is simulated with 65,536 MPI ranks and is able to capture the unique physics of the onset of fluidization.
The MiniFuel irradiation platform at Oak Ridge National Laboratory's High Flux Isotope Reactor (HFIR) is a flexible, high-throughput separate effects test capability. Finite element thermal models are relied upon to design MiniFuel experiments and to achieve experimental objectives. Recent reports show good agreement in the model prediction of target fuel temperatures, but as the capability of the experiments is extended to higher temperatures, the uncertainty in the model predictions must be quantified. To that end, high-impact, high-uncertainty parameters that contribute the most uncertainty to the model are identified. The uncertainty quantification was accomplished through a series of screening and sensitivity analyses. The first analysis utilizes the method of Morris to perform a computationally efficient preliminary screening that considers uncertainty in a large number of the model inputs. The most important parameters identified in the Morris screening study were then considered in a Sobol sensitivity analysis that more robustly ranks and quantifies the uncertainty associated with each parameter. From these analyses, it was determined that thermal contact conductance between components is the parameter that contributes the highest uncertainty. The estimated uncertainty of the MiniFuel model fuel temperature predictions is ±80 °C in the removable beryllium and ±40 °C in the vertical experiment facilities. In conclusion, the framework established by the series of sensitivity analyses presented herein could easily be adapted to fit the needs of accelerated fuel qualification processes.
This is the second article in a three-part series on the construction of 3-node, 9-dof membrane elements with normal-to-its-plane rotational freedoms (the so-called drilling freedoms) using parametrized variational principles. In this part, one such element is derived within the context of the assumed natural deviatoric strain (ANDES) formulation. The higher-order strains are obtained by constructing three parallel-to-sides pure-bending modes from which natural strains are obtained at the corner points and interpolated over the element. To attain rank sufficiency, an additional higher-order 'torsional' mode, corresponding to equal hierarchical rotations at each corner with all other motions precluded, is incorporated. The resulting formulation has five free parameters. When these parameters are optimized against pure bending by energy balance methods, the resulting element is found to coalesce with the optimal EFF element derived in Part I. Numerical integration as a strain filtering device is found to play a key role in this achievement.
Air stripping designs are considered to reduce the presence of volatile organic compounds in the purified water. Components of the wastewater streams are ranked by Henry's Law Constant and the suitability of air stripping in the purification of wastewater in terms of component removal is evaluated. Distillation processes are modeled in tandem with air stripping to demonstrate the potential effectiveness and utility of these methods in recycling wastewater on the Moon. Scaling factors for distillation and air stripping columns are presented to account for the difference in the lunar gravitation environment. Commercially available distillation and air stripping units which are considered suitable for Exploration Life Support are presented. The advantages to the various designs are summarized with respect to water purity levels, power consumption, and processing rates. An evaluation of reactive distillation and air stripping is presented with regards to the reduction of volatile organic compounds in the contaminated water and air. Among the methods presented, an architecture is presented for the evaluation of the simultaneous oxidation of organics in air and water. These and other designs are presented in light of potential improvements in power consumptions and air and water purities for architectures which include catalytic activity integrated into the water processor. In particular, catalytic oxidation of organics may be useful as a tool to remove contaminants that more traditional distillation and/or air stripping columns may not remove. A review of the current leading edge at the commercial level and at the research frontier in catalytically active materials is presented. Themes and directions from the engineering developments in catalyst design are presented conceptually in light of developments in the nanoscale chemistry of a variety of catalyst materials.
The purpose of this paper is to give a construction of representations of the group of currents for semisimple groups of rank greater than one. Such groups have no unitary representations in the Fock space, since the semisimple groups of this form have no nontrivial cohomology in faithful irreducible representations. Thus we first construct cohomology of the semisimple groups in nonunitary representations. The principal method is to reduce all constructions to Iwasawa subgroups (solvable subgroups of the semisimple groups), with subsequent extension to the original group. The resulting representation is realized in the so-called quasi-Poisson Hilbert space associated with natural measures on infinite-dimensional spaces.
Purpose and Objectives: To quantify, through an observer study, the reduction in metal artifacts on cone beam computed tomographic (CBCT) images using a projection-interpolation algorithm, on images containing metal artifacts from dental fillings and implants in patients treated for head and neck (H&N) cancer. Methods and Materials: An interpolation-substitution algorithm was applied to H&N CBCT images containing metal artifacts from dental fillings and implants. Image quality with respect to metal artifacts was evaluated subjectively and objectively. First, 6 independent radiation oncologists were asked to rank randomly sorted blinded images (before and after metal artifact reduction) using a 5-point rating scale (1 = severe artifacts; 5 = no artifacts). Second, the standard deviation of different regions of interest (ROI) within each image was calculated and compared with the mean rating scores. Results: The interpolation-substitution technique successfully reduced metal artifacts in 70% of the cases. From a total of 60 images from 15 H&N cancer patients undergoing image guided radiation therapy, the mean rating score on the uncorrected images was 2.3 ± 1.1, versus 3.3 ± 1.0 for the corrected images. The mean difference in ranking score between uncorrected and corrected images was 1.0 (95% confidence interval: 0.9-1.2, P<.05). The standard deviation of each ROI significantly decreased after artifact reduction (P<.01). Moreover, a negative correlation between the mean rating score for each image and the standard deviation of the oral cavity and bilateral cheeks was observed. Conclusion: The interpolation-substitution algorithm is efficient and effective for reducing metal artifacts caused by dental fillings and implants on CBCT images, as demonstrated by the statistically significant increase in observer image quality ranking and by the decrease in ROI standard deviation between uncorrected and corrected images.