Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set method”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Concepts, Challenges, and Successes in Modeling Thermodynamics of Metabolism

The modeling of the chemical reactions involved in metabolism is a daunting task. Ideally, the modeling of metabolism would use kinetic simulations, but these simulations require knowledge of the thousands of rate constants involved in the reactions. The measurement of rate constants is very labor intensive, and hence rate constants for most enzymatic reactions are not available. Consequently, constraint-based flux modeling has been the method of choice because it does not require the use of the rate constants of the law of mass action. However, this convenience also limits the predictive power of constraint-based approaches in that the law of mass action is used only as a constraint, making it difficult to predict metabolite levels or energy requirements of pathways. An alternative to both of these approaches is to model metabolism using simulations of states rather than simulations of reactions, in which the state is defined as the set of all metabolite counts or concentrations. While kinetic simulations model reactions based on the likelihood of the reaction derived from the law of mass action, states are modeled based on likelihood ratios of mass action. Both approaches provide information on the energy requirements of metabolic reactions and pathways. However, modeling states rather than reactions has the advantage that the parameters needed to model states (chemical potentials) are much easier to determine than the parameters needed to model reactions (rate constants). Herein, we discuss recent results, assumptions, and issues in using simulations of state to model metabolism.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-201-02: Application of RayStretch in Clinical Cases: A Calculation for Heterogeneity Corrections in LDR Permanent I-125 Prostate Brachytherapy

Purpose: Tissue heterogeneities and calcifications have significant impact on the dosimetry of low energy brachytherapy (BT). RayStretch is an analytical algorithm developed in our institution to incorporate heterogeneity corrections in LDR prostate brachytherapy. The aim of this work is to study its application in clinical cases by comparing its predictions with the results obtained with TG-43 and Monte Carlo (MC) simulations. Methods: A clinical implant (71 I-125 seeds, 15 needles) from a real patient was considered. On this patient, different volumes with calcifications were considered. Its properties were evaluated in three ways by i) the Treatment planning system (TPS) (TG-43), ii) a MC study using the Penelope2009 code, and iii) RayStretch. To analyse the performance of RayStretch, calcifications located in the prostate lobules covering 11% of the total prostate volume and larger calcifications located in the lobules and underneath the urethra for a total occupied volume of 30% were considered. Three mass densities (1.05, 1.20, and 1.35 g/cm3) were explored for the calcifications. Therefore, 6 different scenarios ranging from small low density calcifications to large high density ones have been discussed. Results: DVH and D90 results given by RayStretch agree within 1% with the full MC simulations. Although no effort has been done to improve RayStretch numerical performance, its present implementation is able to evaluate a clinical implant in a few seconds to the same level of accuracy as a detailed MC calculation. Conclusion: RayStretch is a robust method for heterogeneity corrections in prostate BT supported on TG-43 data. Its compatibility with commercial TPSs and its high calculation speed makes it feasible for use in clinical settings for improving treatment quality. It will allow in a second phase of this project, its use during intraoperative ultrasound planning. This study was partly supported by a fellowship grant from the Spanish Ministry of Education, by the Generalitat Valenciana under Project PROMETEOII/2013/010, by the Spanish Government under Project No. FIS2013-42156 and by the European Commission within the SeventhFramework Program through ENTERVISION (grant agreement number 264552).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

The Challenges of Air Monitoring at a Soil Remediation Project in a Suburban Setting - 20448

Air monitoring during soil excavation at a chemically and radiologically impacted site in a suburban setting presents unique challenges. Site activities were conducted in a manner that minimized the release of airborne particulates of concern, principally radionuclides of concern (ROC). Dust and air monitoring were performed to ensure the public was not being exposed to the soil contaminants. The project team collected representative air and dust data while being mindful to not impose on the lives of residents in the area. Due to the close proximity to the public, there was a high level of oversight from the client and regulatory agencies. The challenges involved with air sampling at this soil remediation site located in a suburban setting included the number of contaminants, the ranges of risks and hazards presented by the contaminants and the method of remediation, the concerns of the regulators and the public, and the perception of risk from the public. (authors)

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

Uncertainty quantification methodology for hyperbolic systems with application to blood flow in arteries

We present a Stochastic Finite Volume - ADER (SFV-ADER) methodology for Uncertainty Quantification (UQ) in the general framework of systems of hyperbolic balance laws with uncertainty in parameters. The resulting method is weakly intrusive, meaning that deterministic solvers need minor modifications to include random contributions, and has no theoretical accuracy barrier. An illustration of the second-order version for the viscous Burgers equation with uncertain viscosity coefficient is first given in detail, under different deterministic initial conditions; the attainment of the theoretically expected convergence rate is demonstrated empirically; results and features of our scheme are discussed. We then extend the SFV-ADER method to a non-linear hyperbolic system with source terms that models one-dimensional blood flow in arteries, assuming uncertainties in a set of parameters of the problem. Results are then compared with measurements in a well-defined 1:1 experimental replica of the network of largest arteries in the human systemic circulation. For a 99%-confidence level, severe influence of parameter fluctuations is seen in pressure profiles, while secondary effects on flow rate profiles become visible, particularly in terminal branches. In conclusion, results suggest that the proposed methodology could successfully be applied to other physical problems involving hyperbolic balance laws.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Comparison of Spectroscopic Techniques for Determining the Peroxide Value of 19 Classes of Naturally Aged, Plant-Based Edible Oils

The peroxide value of edible oils is a measure of the degree of oxidation, which directly relates to the freshness of the oil sample. Several studies previously reported in the literature have paired various spectroscopic techniques with multivariate analyses to rapidly determine peroxide values using field portable and process instrumentation; those efforts presented “best-case scenarios” with oils from narrowly defined training and test sets. The purpose of this paper is to evaluate the use of near- and mid-infrared absorption and Raman scattering spectroscopies on oil samples from different oil classes, including seasonal and vendor variations, to determine which measurement technique or combination thereof is best for predicting peroxide values. Following peroxide value assays of each oil class using an established titration-based method, global and global-subset calibration models were constructed from spectroscopic data collected on the 19 oil classes used in this study. Spectra from each optical technique were used to create partial least squares regression calibration models to predict the peroxide value of unknown oil samples. A global peroxide value model based on near-infrared (8 mm optical path length) oil measurements produced the lowest RMSEP (4.9), followed by 24 mm optical path length near infrared (5.1), Raman (6.9) and 50 µm optical path length mid-infrared (7.3). However, it was determined that the Raman RMSEP resulted from chance correlations. Global peroxide value models based on low-level fusion of the NIR (8 and 24 mm optical path length) data and all infrared data produced the same RMSEP of 5.1. Global subset models, based on any of the spectroscopies and olive oil training sets from any class (pure, extra light, extra virgin), all failed to extrapolate to the non-olive oils. However, the near-infrared global subset model built on extra virgin olive oil could extrapolate to test samples from other olive oil classes. This work demonstrates the difficulty of developing a truly global method for determining peroxide value of oils.

36 MATERIALS SCIENCE↗
DOE OSTI2016

An integrated model-driven method for in-treatment upper airway motion tracking using cine MRI in head and neck radiation therapy

Purpose: For the first time, MRI-guided radiation therapy systems can acquire cine images to dynamically monitor in-treatment internal organ motion. However, the complex head and neck (H&N) structures and low-contrast/resolution of on-board cine MRI images make automatic motion tracking a very challenging task. In this study, the authors proposed an integrated model-driven method to automatically track the in-treatment motion of the H&N upper airway, a complex and highly deformable region wherein internal motion often occurs in an either voluntary or involuntary manner, from cine MRI images for the analysis of H&N motion patterns. Methods: Considering the complex H&N structures and ensuring automatic and robust upper airway motion tracking, the authors firstly built a set of linked statistical shapes (including face, face-jaw, and face-jaw-palate) using principal component analysis from clinically approved contours delineated on a set of training data. The linked statistical shapes integrate explicit landmarks and implicit shape representation. Then, a hierarchical model-fitting algorithm was developed to align the linked shapes on the first image frame of a to-be-tracked cine sequence and to localize the upper airway region. Finally, a multifeature level set contour propagation scheme was performed to identify the upper airway shape change, frame-by-frame, on the entire image sequence. The multifeature fitting energy, including the information of intensity variations, edge saliency, curve geometry, and temporal shape continuity, was minimized to capture the details of moving airway boundaries. Sagittal cine MR image sequences acquired from three H&N cancer patients were utilized to demonstrate the performance of the proposed motion tracking method. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 50 randomly selected cine image frames from each patient. The resulting average dice similarity coefficient (93.28% ± 1.46%) and margin error (0.49 ± 0.12 mm) showed good agreement between the automatic and manual results. The comparison with three other deformable model-based segmentation methods illustrated the superior shape tracking performance of the proposed method. Large interpatient variations of swallowing frequency, swallowing duration, and upper airway cross-sectional area were observed from the testing cine image sequences. Conclusions: The proposed motion tracking method can provide accurate upper airway motion tracking results, and enable automatic and quantitative identification and analysis of in-treatment H&N upper airway motion. By integrating explicit and implicit linked-shape representations within a hierarchical model-fitting process, the proposed tracking method can process complex H&N structures and low-contrast/resolution cine MRI images. Future research will focus on the improvement of method reliability, patient motion pattern analysis for providing more information on patient-specific prediction of structure displacements, and motion effects on dosimetry for better H&N motion management in radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Can fractal dimensions objectivize gastropod shell morphometrics? A case study from Lake Lugu (SW China)

1. Morphometrics are fundamental for the analysis of size and shape in fossils, particularly because soft parts or DNA are rarely preserved and hard parts such as shells are commonly the only source of information. Geometric morphometrics, that is, landmark analysis, is well established for the description of shape but it exhibits a couple of shortcomings resulting from subjective choices during landmarking (number and position of landmarks) and from difficulties in resolving shape at the level of micro-sculpture. 2. With the aid of high-resolution 3D scanning technology and analyses of fractal dimensions, we test whether such shortcomings of linear and landmark morphometrics can be overcome. As a model group, we selected a clade of modern viviparid gastropods from Lake Lugu, with shells that show a high degree of sculptural variation. Linear and landmark analyses were applied to the same shells in order to establish the fractal dimensions. 3. The genetic diversity of the gastropod clade was assessed. The genetic results suggest that the gastropod clade represents a single species. The results of all morphometric methods applied are in line with the genetic results, which is that no specific morphotype could be delimited. Apart from this overall agreement, landmark and fractal dimension analyses do not correspond to each other but represent data sets with different information. Generally, the fractal dimension values quantify the roughness of the shell surface, the resolution of the 3D scans determining the level. In our approach, we captured the micro-sculpture but not the first-order sculptural elements, which explains that fractal dimension and landmark data are not in phase. 4. We can show that analyzing fractal dimensions of gastropod shells opens a window to more detailed information that can be considered in evolutionary and ecological contexts. We propose that using low-resolution 3D scans may successfully substitute landmark analyses because it overcomes the subjective landmarking. Analyses of 3D scans with higher resolution than used in this study will provide surface roughness information at the mineralogical level. We suggest that fractal dimension analyses of a combination of differently resolved 3D models will significantly improve the quality of shell morphometrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Positivity Limiters for Filtered Spectral Approximations of Linear Kinetic Transport Equations

We analyze the properties and compare the performance of several positivity limiters for spectral approximations with respect to the angular variable of linear transport equations. It is well-known that spectral methods suffer from the occurrence of (unphysical) negative spatial particle concentrations due to the fact that the underlying polynomial approximations are not always positive at the kinetic level. Positivity limiters address this defect by enforcing positivity of the polynomial approximation on a finite set of preselected points. With a proper PDE solver, they ensure positivity of the particle concentration at each step in a time integration scheme. We review several known positivity limiters proposed in other contexts and also introduce a modification for one of them. We give error estimates for the consistency of the positive approximations produced by these limiters and compare the theoretical estimates to numerical results. We then solve two benchmark problems with these limiters, make qualitative and quantitative observations about the solutions, and then compare the efficiency of the different limiters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Differential cross sections for electron impact excitation of the electronic bands of phenol

We report results from a joint theoretical and experimental investigation into electron scattering from the important organic species phenol (C{sub 6}H{sub 5}OH). Specifically, differential cross sections (DCSs) have been measured and calculated for the electron-impact excitation of the electronic states of C{sub 6}H{sub 5}OH. The measurements were carried out at energies in the range 15–40 eV, and for scattered-electron angles between 10{sup ∘} and 90{sup ∘}. The energy resolution of those experiments was typically ∼80 meV. Corresponding Schwinger multichannel method with pseudo-potentials calculations, with and without Born-closure, were also performed for a sub-set of the excited electronic-states that were accessed in the measurements. Those calculations were conducted at the static exchange plus polarisation (SEP)-level using a minimum orbital basis for single configuration interaction (MOBSCI) approach. Agreement between the measured and calculated DCSs was typically fair, although to obtain quantitative accord, the theory would need to incorporate even more channels into the MOBSCI.

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

Going Forward with the Nancy Grace Roman Space Telescope Transient Survey: Validation of Precision Forward-modeling Photometry for Undersampled Imaging

The Nancy Grace Roman Space Telescope (Roman) is an observatory for both wide-field observations and coronagraphy that is scheduled for launch in the mid-2020s. Part of the planned survey is a deep, cadenced field or fields that enable cosmological measurements with type Ia supernovae (SNe Ia). With a pixel scale of 0$_.^"$11, the Wide Field Instrument will be undersampled, presenting a difficulty for precisely subtracting the galaxy light underneath the SNe. We use simulated data to validate the ability of a forward-model code (such codes are frequently also called “scene-modeling” codes) to perform precision supernova photometry for the Roman SN survey. Our simulation includes over 760,000 image cutouts around SNe Ia or host galaxies (~10% of a full-scale survey). To have a realistic 2D distribution of underlying galaxy light, we use the VELA simulated high-resolution images of galaxies. We run each set of cutouts through our forward-modeling code which automatically measures time-dependent SN fluxes. Given our assumed inputs of a perfect model of the instrument point-spread functions and calibration, we find biases at the millimagnitude level from this method in four red filters (Y106, J129, H158, and F184), easily meeting the 0.5% Roman inter-filter calibration requirement for a cutting-edge measurement of cosmological parameters using SNe Ia. Simulated data in the bluer Z087 filter shows larger ~ 2–3 mmg biases, also meeting this requirement, but with more room for improvement. Our forward-model code has been released on Zenodo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-186: A Treatment Planning Study of Normal Tissue Sparing with Robustness Optimized IMPT, 4Pi IMRT, and VMAT for Head and Neck Cases

Purpose: We performed a retrospective dosimetric comparison study between the robustness optimized Intensity Modulated Proton Therapy (RO-IMPT), volumetric-modulated arc therapy (VMAT), and the non-coplanar 4? intensity modulated radiation therapy (IMRT). These methods represent the most advanced radiation treatment methods clinically available. We compare their dosimetric performance for head and neck cancer treatments with special focus on the OAR sparing near the tumor volumes. Methods: A total of 11 head and neck cases, which include 10 recurrent cases and one bilateral case, were selected for the study. Different dose levels were prescribed to tumor target depending on disease and location. Three treatment plans were created on commercial TPS systems for a novel noncoplanar 4π method (20 beams), VMAT, and RO-IMPT technique (maximum 4 fields). The maximum patient positioning error was set to 3 mm and the maximum proton range uncertainty was set to 3% for the robustness optimization. Line dose profiles were investigated for OARs close to tumor volumes. Results: All three techniques achieved 98% coverage of the CTV target and most photon plans had less than 110% of the hot spots. The RO-IMPT plans show superior tumor dose homogeneity than 4? and VMAT plans. Although RO-IMPT has greater R50 dose spillage to the surrounding normal tissue than 4π and VMAT, the RO-IMPT plans demonstrate better or comparable OAR (parotid, mandible, carotid, oral cavity, pharynx, and etc.) sparing for structures closely abutting tumor targets. Conclusion: The RO-IMPT’s ability of OAR sparing is benchmarked against the C-arm linac based non-coplanar 4π technique and the standard VMAT method. RO-IMPT consistently shows better or comparable OAR sparing even for tissue structures closely abutting treatment target volume. RO-IMPT further reduces treatment uncertainty associated with proton therapy and delivers robust treatment plans to both unilateral and bilateral head and neck cancer patients with desirable treatment time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-G-204-05: The Effects of CT Acquisition and Reconstruction Conditions On Computed Texture Feature Values of Lung Lesions

Purpose: Texture features have been investigated as a biomarker of response and malignancy. Because these features reflect local differences in density, they may be influenced by acquisition and reconstruction parameters. The purpose of this study was to investigate the effects of radiation dose level and reconstruction method on features derived from lung lesions. Methods: With IRB approval, 33 lung tumor cases were identified from clinically indicated thoracic CT scans in which the raw projection (sinogram) data were available. Based on a previously-published technique, noise was added to the raw data to simulate reduced-dose versions of each case at 25%, 10% and 3% of the original dose. Original and simulated reduced dose projection data were reconstructed with conventional and two iterative-reconstruction settings, yielding 12 combinations of dose/recon conditions. One lesion from each case was contoured. At the reference condition (full dose, conventional recon), 17 lesions were randomly selected for repeat contouring (repeatability). For each lesion at each dose/recon condition, 151 texture measures were calculated. A paired differences approach was employed to compare feature variation from repeat contours at the reference condition to the variation observed in other dose/recon conditions (reproducibility). The ratio of standard deviation of the reproducibility to repeatability was used as the variation measure for each feature. Results: The mean variation (standard deviation) across dose levels and kernel was significantly different with a ratio of 2.24 (±5.85) across texture features (p=0.01). The mean variation (standard deviation) across dose levels with conventional recon was also significantly different with 2.30 (7.11) (p=0.025). The mean variation across reconstruction settings of original dose has a trend in showing difference with 1.35 (2.60) among all features (p=0.09). Conclusion: Texture features varied considerably with variations in dose and reconstruction condition. Care should be taken to standardize these conditions when using texture as a quantitative feature. This effort supported in part by a grant from the National Cancer Institute’s Quantitative Imaging Network (QIN): U01 CA181156; The UCLA Department of Radiology has a Master Research Agreement with Siemens Healthcare; Dr. McNitt-Gray has previously received research support from Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Power increases using wind direction spatial filtering for wind farm control: Evaluation using FLORIS, modified for dynamic settings

As wind energy plays a growing role in the energy sector, new methods for controlling wind turbines and wind farms to maximize performance are garnering industry interest. A developing body of research treats the entire wind farm as a control system, with individual turbines acting as agents in a network, allowing farm-level objectives to be considered. Two promising developments in this research are wake steering control, which seeks to increase the power generated at a wind farm by directing the wakes of upstream turbines away from downstream ones, and communication-based spatial filtering, which seeks to improve the quality of information used by turbine- and farm-level controllers by combining measurements of the wind field collected at the individual turbines. The latter method has been shown to improve the estimates of wind direction at the turbines; however, the resulting potential for increased power capture warrants further investigation. With this paper, we begin to address this gap by combining wake steering with wind direction spatial filtering. To do so, we present a preliminary method for assessing the power capture of dynamic controllers using wind farm codes designed for time-averaged simulation. This allows us to generate results much more rapidly than would be possible using high-fidelity wind farm simulators and may be useful in many wind farm control design applications.

17 WIND ENERGY↗
DOE OSTI2024

Discovering invariant spatial features in electron energy loss spectroscopy images on the mesoscopic and atomic levels

Over the last two decades, Electron Energy Loss Spectroscopy (EELS) imaging with a scanning transmission electron microscope has emerged as a technique of choice for visualizing complex chemical, electronic, plasmonic, and phononic phenomena in complex materials and structures. The availability of the EELS data necessitates the development of methods to analyze multidimensional data sets with complex spatial and energy structures. Traditionally, the analysis of these data sets has been based on analysis of individual spectra, one at a time, whereas the spatial structure and correlations between individual spatial pixels containing the relevant information of the physics of underpinning processes have generally been ignored and analyzed only via the visualization as 2D maps. Here, we develop a machine learning-based approach and workflows for the analysis of spatial structures in 3D EELS data sets using a combination of dimensionality reduction and multichannel rotationally invariant variational autoencoders. This approach is illustrated for the analysis of both the plasmonic phenomena in a system of nanowires and in the core excitations in functional oxides using low loss and core-loss EELS, respectively. The code developed in this manuscript is open sourced and freely available and provided as a Jupyter notebook for the interested reader.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Use of Monte Carlo simulation to quantify CIPS risk

Crud Induced Power Shift (CIPS), caused by the accumulation of boron compounds in corrosion product deposits on high power fuel assemblies, is an important operational issue for Pressurized Water Reactors (PWRs). The Boron-induced Offset Analysis (BOA) computer code has been developed by EPRI in collaboration with Westinghouse to help utility engineers avoid onset of CIPS. The code, used with an established procedure, predicts the maximum boron mass that will be accumulated in the upcoming cycle for a given core design and plant chemistry and operating program. The calculated boron mass above a given threshold predicts the risk of CIPS to the upcoming cycle. BOA has been successful in correctly predicting the risk for onset or lack of CIPS in many cycles. However, like most codes modeling complex physical phenomena, it is subject to uncertainties in several of the inputs used to describe these phenomena. The current BOA procedure results in a single 'best estimate' prediction, which has been judged to provide engineers with a practical way to assess CIPS risk. As such, however, there is a probability that the actual plant CIPS performance may be underestimated. This paper describes a method to assess the probability (with an associated confidence level) of no CIPS, mild to moderate CIPS, or severe CIPS using BOA for a given core design. This is accomplished through the use of a Monte Carlo simulation. The simulation is applied to several plant examples, and the results are evaluated in terms of the likelihood of mild or severe CIPS. A suggested set of criteria for decision making regarding CIPS risk is presented in the conclusion. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Accurate Binding Energies for Lithium Polysulfides and Assessment of Density Functionals for Lithium–Sulfur Battery Research

Lithium–sulfur batteries have high theoretical energy density, but a better knowledge of their intimate structural details will be helpful in improving their conductivity and long-term cycling behavior. In order to identify the stationary configurations of lithium polysulfides (Li 2 S n , 2 ≤ n ≤ 8) formed in the charging and discharging processes of the lithium–sulfur batteries, ab initio molecular dynamics was employed to sample the configuration space of Li 2 S n , followed by optimization of structures by CCSD(T)-F12b/aug-cc-pVDZ. Here, using the optimized stationary points, we have created the LiSAE38 benchmark database of atomization energies (AEs) of 38 lithium polysulfide isomers by using the higher-level WMS and W3X-L methods. In addition, the performances of 39 density functionals have been assessed against the benchmark AEs and the relative stabilities of Li 2 S n isomers. On the basis of the assessments with the def2-QZVP basis set, the PW6B95, B97-1, B3LYP-D3, TPSS, and DSD-PBEP86 density functionals are the most accurate for the AEs, whereas the mPW2-PLYP-D, DSD-PBEP86, PW6B95, HSE06, and PBEQIDH functionals are most accurate for the relative energies. Local functionals are of special interest because of their lower cost (faster timings in terms of computer resources) for large systems; among the tested local functionals, MN15-L and revM06-L are the most accurate for the calculations of relative stabilities.

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

The Ground State Electronic Energy of Benzene

We report here on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground-state energy of the benzene molecule in a standard correlation-consistent basis set of double-ζ quality. As a broad international endeavor, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around -863 mE H . However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE H ), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

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

From minimum-viable-products to full models: a step-wise development of diagnostic forward models in support of design, analysis and modelling on the ST40 tokamak

Like most magnetic confined fusion experiments, the ST40 tokamak started off with a small subset of diagnostics and gradually increased the diagnostic set to include more complex and comprehensive systems. To make the most of each operational phase, forward models of various diagnostics are used and developed to aid design, provide consistency-checks during commissioning, test analysis methods, and build workflows to constrain high-level parameters to inform interpretation, theory and modelling. For new models and new analysis workflows, minimum-viable-products are released early, and their complexity is increased in a step-wise manner, facilitating the support of all programme phases on multiple parallel applications, while enabling learning opportunities and feedback loops. In this contribution we review the philosophy, scope and architecture of the framework under development. We discuss the details of some forward models, with examples on how they are used to aid diagnostic design, to investigate analysis methodologies through synthetic data, and how they are embedded in experimental analysis workflows. We compare previously published experimental results with new, more advanced analysis workflows employing more recent, detailed models and new diagnostic data, providing confirmation of the published material from the 2021–22 experimental campaign.

integrated data analysis↗