SEARCH · Engineering Papers

Results for “BODY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2017

DISCO interacting protein 2 determines direction of axon projection under the regulation of c-Jun N-terminal kinase in the Drosophila mushroom body

Precisely controlled axon guidance for complex neuronal wiring is essential for appropriate neuronal function. c-Jun N-terminal kinase (JNK) was found to play a role in axon guidance recently as well as in cell proliferation, protection and apoptosis. In spite of many genetic and molecular studies on these biological processes regulated by JNK, how JNK regulates axon guidance accurately has not been fully explained thus far. To address this question, we use the Drosophila mushroom body (MB) as a model since the α/β axons project in two distinct directions. Here we show that DISCO interacting protein 2 (DIP2) is required for the accurate direction of axonal guidance. DIP2 expression is under the regulation of Basket (Bsk), the Drosophila homologue of JNK. We additionally found that the Bsk/DIP2 pathway is independent from the AP-1 transcriptional factor complex pathway, which is directly activated by Bsk. In conclusion, our findings revealed DIP2 as a novel effector downstream of Bsk modulating the direction of axon projection. - Highlights: • DIP2 is required for accurate direction of axon guidance in Drosophila mushroom body. • DIP2 is a downstream of JNK in the axon guidance of Drosophila mushroom body neuron. • JNK/DIP2 pathway is independent from JNK/AP-1 transcriptional factor complex pathway.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Water and heat: New constraints on the evolution of the CV chondrite parent body

This paper focuses on the post-accretion history of CV3 chondrites, through a combination of petrographic and mineralogical characterization, magnetic measurements, spectral (Raman and Infrared) and thermo-gravimetric analysis of 31 meteorites (including 7 falls, 21 Antarctic and 3 non-Antarctic finds) spanning a wide metamorphic range. We classify the 21 Antarctic chondrites and the Bukhara fall into the CV Red , CV OxA , and CV OxB subgroups. We establish quantitative parameters relevant for this sub-classification. In comparison to CV Ox , CV Red chondrites are characterized by (i) a lower abundance of matrix, (ii) a higher abundance of metal, (iii) the presence of Ni-poor sulfides. In comparison to CV OxB , CV OxA are characterized by (i) similar matrix abundance, (ii) a higher abundance of metal, (iii) the presence of metal almost exclusively under the form of awaruite, (iv) lower Ni content of sulfides, (v) lower magnetic susceptibility and saturation remanence. In this work, both CV Ox (CV OxA and CV OxB ) and CV Red experienced aqueous alteration, and contain oxyhydroxides and phyllosilicates. We show that the abundance of these hydrated secondary minerals observed today in individual CV chondrites decreases with their peak metamorphic temperature. This is interpreted either as partial dehydration of these secondary minerals or limited hydration due to the rapid exhaustion of the water reservoir during parent body thermal metamorphism. Moreover, the lower abundance of oxyhydroxides (that have a lower thermal stability than phyllosilicates and may in large part postdate the peak of thermal metamorphism) in more metamorphosed CV chondrites is interpreted as lower availability of aqueous fluids during retrograde metamorphism in these meteorites. Lastly, we show that in comparison to CV OxB , CV OxA are systematically (i) more metamorphosed, (ii) less hydrated, (iii) depleted in ferromagnetic minerals, (iv) but enriched in metal in the form of secondary awaruite. CV OxA may be thermally metamorphosed CV OxB . CV Red are significantly different from CV Ox (matrix abundances, alteration products, opaque minerals), but span the same wide metamorphic range. This could be indicative of a laterally heterogeneous CV parent body, or suggest the existence of distinct parent bodies for CV Ox and CV Red chondrites.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Impact of Pretreatment Body Mass Index on Patients With Head-and-Neck Cancer Treated With Radiation

Purpose: To investigate the association of pretreatment body mass index (preT BMI) with outcomes of head-and-neck cancer in patients treated with radiotherapy (RT). Methods and Materials: All 1,562 patients diagnosed with head-and-neck cancer and treated with curative-intent RT to a dose of 60 Gy or higher were retrospectively studied. Body weight was measured both at entry and at the end of RT. Cancer-specific survival (CSS), overall survival (OS), locoregional control (LRC), and distant metastasis (DM) were analyzed by preT BMI (<25 kg/m{sup 2} vs. {>=}25 kg/m{sup 2}). The median follow-up was 8.6 years. Results: Patients with lower preT BMI were statistically significantly associated with poorer CSS and OS than those with higher preT BMI. There was no significant difference between preT BMI groups in terms of LRC and DM. Body weight loss (BWL) during radiation did not influence survival outcomes. However, in the group with higher preT BMI, CSS, OS, and DM-free survival of patients with less BWL during radiation were statistically longer when compared with greater BWL. Conclusion: This study demonstrates that higher preT BMI positively influenced survival outcomes for patients with head-and-neck cancer. Patients with higher preT BMI who were able to maintain their weight during radiation had significantly better survival than patients with greater BWL.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Stereotactic Ablative Body Radiation Therapy for Octogenarians With Non-Small Cell Lung Cancer

Purpose: To retrospectively investigate treatment outcomes of stereotactic ablative body radiation therapy (SABR) for octogenarians with non-small cell lung cancer (NSCLC). Methods and Materials: Between 2005 and 2012, 109 patients aged ≥80 years with T1-2N0M0 NSCLC were treated with SABR: 47 patients had histology-unproven lung cancer; 62 patients had pathologically proven NSCLC. The prescribed doses were either 50 Gy/5 fractions for peripheral tumors or 40 Gy/5 fractions for centrally located tumors. The treatment outcomes, toxicities, and the correlating factors for overall survival (OS) were evaluated. Results: The median follow-up duration after SABR was 24.2 (range, 3.0-64.6) months. Only limited toxicities were observed, except for 1 grade 5 radiation pneumonitis. The 3-year local, regional, and distant metastasis-free survival rates were 82.3%, 90.1%, and 76.8%, respectively. The OS and lung cancer-specific survival rates were 53.7% and 70.8%, respectively. Multivariate analysis revealed that medically inoperable, low body mass index, high T stage, and high C-reactive protein were the predictors for short OS. The OS for the operable octogenarians was significantly better than that for inoperable (P<.01). Conclusions: Stereotactic ablative body radiation therapy for octogenarians was feasible, with excellent OS. Multivariate analysis revealed that operability was one of the predictors for OS. For medically operable octogenarians with early-stage NSCLC, SABR should be prospectively compared with resection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Self-Management of Patient Body Position, Pose, and Motion Using Wide-Field, Real-Time Optical Measurement Feedback: Results of a Volunteer Study

Purpose: We present the results of a clinical feasibility study, performed in 10 healthy volunteers undergoing a simulated treatment over 3 sessions, to investigate the use of a wide-field visual feedback technique intended to help patients control their pose while reducing motion during radiation therapy treatment. Methods and Materials: An optical surface sensor is used to capture wide-area measurements of a subject's body surface with visualizations of these data displayed back to them in real time. In this study we hypothesize that this active feedback mechanism will enable patients to control their motion and help them maintain their setup pose and position. A capability hierarchy of 3 different level-of-detail abstractions of the measured surface data is systematically compared. Results: Use of the device enabled volunteers to increase their conformance to a reference surface, as measured by decreased variability across their body surfaces. The use of visual feedback also enabled volunteers to reduce their respiratory motion amplitude to 1.7 ± 0.6 mm compared with 2.7 ± 1.4 mm without visual feedback. Conclusions: The use of live feedback of their optically measured body surfaces enabled a set of volunteers to better manage their pose and motion when compared with free breathing. The method is suitable to be taken forward to patient studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

More Accurate Definition of Clinical Target Volume Based on the Measurement of Microscopic Extensions of the Primary Tumor Toward the Uterus Body in International Federation of Gynecology and Obstetrics Ib-IIa Squamous Cell Carcinoma of the Cervix

Purpose: To more accurately define clinical target volume for cervical cancer radiation treatment planning by evaluating tumor microscopic extension toward the uterus body (METU) in International Federation of Gynecology and Obstetrics stage Ib-IIa squamous cell carcinoma of the cervix (SCCC). Patients and Methods: In this multicenter study, surgical resection specimens from 318 cases of stage Ib-IIa SCCC that underwent radical hysterectomy were included. Patients who had undergone preoperative chemotherapy, radiation, or both were excluded from this study. Microscopic extension of primary tumor toward the uterus body was measured. The association between other pathologic factors and METU was analyzed. Results: Microscopic extension toward the uterus body was not common, with only 12.3% of patients (39 of 318) demonstrating METU. The mean (±SD) distance of METU was 0.32 ± 1.079 mm (range, 0-10 mm). Lymphovascular space invasion was associated with METU distance and occurrence rate. A margin of 5 mm added to gross tumor would adequately cover 99.4% and 99% of the METU in the whole group and in patients with lymphovascular space invasion, respectively. Conclusion: According to our analysis of 318 SCCC specimens for METU, using a 5-mm gross tumor volume to clinical target volume margin in the direction of the uterus should be adequate for International Federation of Gynecology and Obstetrics stage Ib-IIa SCCC. Considering the discrepancy between imaging and pathologic methods in determining gross tumor volume extent, we recommend a safer 10-mm margin in the uterine direction as the standard for clinical practice when using MRI for contouring tumor volume.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2015

Outcomes for Spine Stereotactic Body Radiation Therapy and an Analysis of Predictors of Local Recurrence

Purpose: To investigate local control, survival outcomes, and predictors of local relapse for patients treated with spine stereotactic body radiation therapy. Methods and Materials: We reviewed the records of 332 spinal metastases consecutively treated with stereotactic body radiation therapy between 2002 and 2012. The median follow-up for all living patients was 33 months (range, 0-111 months). Endpoints were overall survival and local control (LC); recurrences were classified as either in-field or marginal. Results: The 1-year actuarial LC and overall survival rates were 88% and 64%, respectively. Patients with local relapses had poorer dosimetric coverage of the gross tumor volume (GTV) compared with patients without recurrence (minimum dose [Dmin] biologically equivalent dose [BED] 23.9 vs 35.1 Gy, P<.001; D98 BED 41.8 vs 48.1 Gy, P=.001; D95 BED 47.2 vs 50.5 Gy, P=.004). Furthermore, patients with marginal recurrences had poorer prescription coverage of the GTV (86% vs 93%, P=.01) compared with those with in-field recurrences, potentially because of more upfront spinal canal disease (78% vs 24%, P=.001). Using a Cox regression univariate analysis, patients with a GTV BED Dmin ≥33.4 Gy (median dose) (equivalent to 14 Gy in 1 fraction) had a significantly higher 1-year LC rate (94% vs 80%, P=.001) compared with patients with a lower GTV BED Dmin; this factor was the only significant variable on multivariate Cox analysis associated with LC (P=.001, hazard ratio 0.29, 95% confidence interval 0.14-0.60) and also was the only variable significant in a separate competing risk multivariate model (P=.001, hazard ratio 0.30, 95% confidence interval 0.15-0.62). Conclusions: Stereotactic body radiation therapy offers durable control for spinal metastases, but there is a subset of patients that recur locally. Patients with local relapse had significantly poorer tumor coverage, which was likely attributable to treatment planning directives that prioritized the spinal cord constraints over tumor coverage. When possible, we recommend maintaining a GTV Dmin above 14 Gy in 1 fraction and 21 Gy in 3 fractions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2013

Automatic off-body overset adaptive Cartesian mesh method based on an octree approach

This paper describes a method for generating adaptive structured Cartesian grids within a near-body/off-body mesh partitioning framework for the flow simulation around complex geometries. The off-body Cartesian mesh generation derives from an octree structure, assuming each octree leaf node defines a structured Cartesian block. This enables one to take into account the large scale discrepancies in terms of resolution between the different bodies involved in the simulation, with minimum memory requirements. Two different conversions from the octree to Cartesian grids are proposed: the first one generates Adaptive Mesh Refinement (AMR) type grid systems, and the second one generates abutting or minimally overlapping Cartesian grid set. We also introduce an algorithm to control the number of points at each adaptation, that automatically determines relevant values of the refinement indicator driving the grid refinement and coarsening. An application to a wing tip vortex computation assesses the capability of the method to capture accurately the flow features.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Experimental and numerical investigations of a two-body floating-point absorber wave energy converter in regular waves

This paper presents experimental and numerical studies on the hydrodynamics of a two-body floating-point absorber (FPA) wave energy converter (WEC) under both extreme and operational wave conditions. In this study, the responses of the WEC in heave, surge, and pitch were evaluated for various regular wave conditions. For extreme condition analysis, we assume the FPA system has a survival mode that locks the power-take-off (PTO) mechanism in extreme waves, and the WEC moves as a single body in this scenario. A series of Reynolds-averaged Navier-Stokes (RANS) simulations was performed for the survival condition analysis, and the results were validated with the measurements from experimental wave tank tests. For the FPA system in operational conditions, both a boundary element method (BEM) and the RANS simulations were used to analyze the motion response and power absorption performance. Additional viscous damping, primarily induced by flow separation and vortex shedding, is included in the BEM simulation as a quadratic drag force proportional to the square of body velocity. The inclusion of viscous drag improves the accuracy of BEM's prediction of the heave response and the power absorption performance of the FPA system, which agrees well with experimental data and the RANS simulation results over a broad range of incident wave periods, except near resonance in large wave height scenarios. Overall, the experimental and numerical results suggest that nonlinear effects, caused by viscous damping and interaction between waves and the FPA, significantly influence the system response and power absorption performance. Nonlinear effects were found to be particularly significant when the wave height was large and the period was near or shorter than the resonant period of the FPA system. The study is expected to be helpful for understanding the nonlinear interaction between waves and the FPA system so that the structure of the FPA can be adequately designed.

16 TIDAL AND WAVE POWER↗

DOE OSTI2016

Implementation of a volumetric modulated arc therapy treatment planning solution for kidney and adrenal stereotactic body radiation therapy

To develop a volumetric modulated arc therapy (VMAT) treatment planning solution in the treatment of primary renal cell carcinoma and oligometastatic adrenal lesions with stereotactic body radiation therapy. Single-arc VMAT plans (n = 5) were compared with clinically delivered step-and-shoot intensity-modulated radiotherapy (IMRT) with planning target volume coverage normalized between techniques. Target volume conformity, organ-at-risk (OAR) dose, treatment time, and monitor units were compared. A VMAT planning solution, created from a combination of arc settings and optimization constraints, auto-generated treatment plans in a single optimization. The treatment planning solution was evaluated on 15 consecutive patients receiving kidney and adrenal stereotactic body radiation therapy. Treatment time was reduced from 13.0 ± 2.6 to 4.0 ± 0.9 minutes for IMRT and VMAT, respectively. The VMAT planning solution generated treatment plans with increased target homogeneity, improved 95% conformity index, and a reduced maximum point dose to nearby OARs but with increased intermediate dose to distant OARs. The conformity of the 95% isodose improved from 1.32 ± 0.39 to 1.12 ± 0.05 for IMRT and VMAT treatment plans, respectively. Evaluation of the planning solution showed clinically acceptable dose distributions for 13 of 15 cases with tight conformity of the prescription isodose to the planning target volume of 1.07 ± 0.04, delivering minimal dose to OARs. The introduction of a stereotactic body radiation therapy VMAT treatment planning solution improves the efficiency of planning and delivery time, producing treatment plans of comparable or superior quality to IMRT in the case of primary renal cell carcinoma and oligometastatic adrenal lesions.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2021

EGLP-1 lowers body weight better than exendin-4 by reducing food intake and increasing basal energy expenditure in diet-induced obese mice

It is well known that GLP-1 activates GLP-1R to reduce body weight by inhibiting eating. GLP-1 is cleaved by the neutral endopeptidase (NEP) 24.11 into a pentapeptide GLP-1 (32–36) amide, which increases basal energy expenditure and inhibits weight gain in obese mice. It is well known that GLP-1 analogs can reduce weight by suppressing eating. However, there are few reports of reducing weight through the dual effects of inhibiting eating and increasing basic energy. Here, we report the peptide EGLP-1, a GLP-1 analogue, which can reduce food intake and increase basal energy expenditure. In C2C12 myotubes, EGLP-1 can increase both phosphorylation of acetyl CoA carboxylase (ACC) and the ratio between phosphorylation of ACC and the total expression of ACC (pACC/ACC). In diet-induced obese mice, EGLP-1 is more effective than exendin-4 in reducing body weight, reducing fat mass and improving hepatic steatosis. At the same time, EGLP-1 can improve hyperglycemia, reduce food intake, and improve insulin resistance, just like exendin-4. In addition, EGLP-1, not exendin-4, can improve physiological parameters associated with lipid metabolism and increase oxygen consumption by increasing uncoupling proteins 3 (UCP3) expression and pACC/ACC ratio in skeletal muscle. Taken together, this data showed that EGLP-1 is able to reduce body weight by reducing food intake and increasing basal energy expenditure, suggesting it may be more effective in treating diabetic and non-diabetic overweight or obese people than pure GLP-1R agonist exendin-4.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2012

Dynamical separation of spherical bodies in supersonic flow

An experimental and computational investigation of the unsteady separation behaviour of two spheres in a highly supersonic flow is carried out. The spherical bodies, initially touching, are released with negligible relative velocity, an arrangement representing the idealized binary fragmentation of a meteoritic body in the atmosphere. In experiments performed in a Mach-4 Ludwieg tube, nylon spheres are initially suspended in the test section by weak threads and, following detachment of the threads by the arrival of the flow, fly freely according to the aerodynamic forces experienced. The resulting sphere motions and unsteady flow structures are recorded using high-speed shadowgraphy. The qualitative separation behaviour and the final lateral velocity of the smaller sphere are found to vary strongly with both the radius ratio and the initial alignment angle of the two spheres. More disparate radii and initial configurations in which the smaller sphere centre lies downstream of that of the larger sphere both increase the tendency for the smaller sphere to be entrained within the flow region bounded by the bow-shock of the larger body, rather than expelled from this region. At a critical angle for a given radius ratio (or a critical radius ratio for a given angle), transition from entrainment to expulsion occurs; at this critical value, the final lateral velocity is close to maximum due to the surfing effect noted by Laurence & Deiterding (J. Fluid Mech., vol. 676, 2011, pp. 396-431) at hypersonic Mach numbers. A high-precision tracking algorithm is used to provide quantitative comparisons between experiments and high-resolution inviscid numerical simulations, with generally favourable agreement.

Laurence, Stuart J↗

DOE OSTI2022

High Accuracy Transition Metal Effective Cores for the Many-Body Diffusion Monte Carlo Method

Practical applications of the real-space diffusion Monte Carlo (DMC) method require the removal of core electrons, where currently localization approximations of semilocal potentials are generally used in the projector. Accurate calculations of complex solids and large molecules demand minimizing the impact of approximated atomic cores. Prior works have shown that the errors from such approximations can be sizable in both finite and periodic systems. In this work, we show that a class of differential pseudopotentials, known as pseudo-Hamiltonians, can be constructed for the 3d transition metal atoms, entirely removing the need for any localization scheme in the DMC projector. As a proof of principle, we demonstrate the approach for the case of Co. In order to minimize errors in the pseudo-Hamiltonian at the many-body level, we generalize the recently proposed correlation-consistent pseudopotential generation scheme to successively close semilocal representations of the differential potentials. Our generation scheme successfully produces potentials tailored specifically for real space projector quantum Monte Carlo methods with low error at the many-body level, i.e., with many-body scattering properties very close to relativistic all-electron results. In particular, we show that the agreement with respect to atomic and molecular quantities reach chemical accuracy in many cases-on par with the most accurate semilocal pseudopotentials available. Further, our pseudo-Hamiltonian generation scheme utilizes standard quantum chemistry codes designed only to work with semilocal pseudopotentials, enabling straightforward generation of pseudo-Hamiltonians for additional elements in future works.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

The Importance of Three-Body Interactions in Molecular Dynamics Simulations of Water with the Fragment Molecular Orbital Method

The analytic first derivative with respect to nuclear coordinates is formulated and implemented in the framework of the three-body fragment molecular orbital (FMO) method. The gradient has been derived and implemented for restricted Hartree-Fock, second-order Møller-Plesset perturbation, and density functional theories. The importance of the three-body fully analytic gradient is illustrated through the failure of the two-body FMO method during molecular dynamics simulations of a small water cluster. The parallel implementation of the fragment molecular orbital method, its parallel efficiency, and its scalability on the Blue Gene/Q architecture up to 262,144 CPU cores, are also discussed.

Pruitt, Spencer R.↗

DOE OSTI2026

Hierarchical Truncations for Many-Body Expansion Potentials

In this work, a new strategy to truncate high-order terms in the many-body expansion (MBE) is proposed. This new approach, which we call a hierarchical many-body expansion (HMBE), is based on a hierarchical partition of the system into multitier fragments and can in principle be applied to any large molecular system. Numerical tests on a series of (H 2 O) 64 structures are presented, demonstrating satisfactory relative energies between the clusters and binding energies of individual clusters compared with full-cluster calculations, with significantly fewer high-order terms computed than conventional MBE. The hierarchical truncation can be augmented by certain many-body terms for fragments at the interface between the partitions (called “Schengen terms”) to further improve accuracy. This work establishes the HMBE scheme as a promising framework to model very large systems (e.g., proteins), which are naturally built on a hierarchical structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Energy-Screened Many-Body Expansion: A Practical Yet Accurate Fragmentation Method for Quantum Chemistry

We introduce an implementation of the truncated many-body expansion, MBE(n), in which the n-body corrections are screened using the effective fragment potential force field, and only those that exceed a specified energy threshold are computed at a quantum-mechanical level of theory. This energy-screened MBE(n) approach is tested at the n = 3 level for a sequence of water clusters, (H2O) N=6–34 . A threshold of 0.25 kJ/mol eliminates more than 80% of the subsystem electronic structure calculations and is even more efficacious in that respect than is distance-based screening. Even so, the energy-screened MBE(3) method is faithful to a full-system quantum chemistry calculation to within 1–2 kJ/mol/monomer, even in good quality basis sets such as aug-cc-pVTZ. These errors can be reduced by means of a two-layer approach that involves a Hartree–Fock calculation for the entire cluster. Such a correction proves to be necessary in order to obtain accurate relative energies for conformational isomers of (H2O)20, but the cost of a full-system Hartree–Fock calculation remains smaller than the cost of three-body subsystem calculations at correlated levels of theory. At the level of second-order Møller–Plesset perturbation theory (MP2), a screened MBE(3) calculation plus a full-system Hartree–Fock calculation is less expensive than a full-system MP2 calculation starting at N = 12 water molecules. This is true even if all MBE(3) subsystem calculations are performed on a single 40-core compute node, i.e., without significant parallelization. Energy-screened MBE(n) thus provides a fragment-based method that is accurate, stable in large basis sets, and low in cost, even when the latter is measured in aggregate computer time.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Data-Driven Many-Body Models for Molecular Fluids: CO 2 /H 2 O Mixtures as a Case Study

In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO2) and water (H2O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born–Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO 2 –CO 2 and CO 2 –H 2 O dimers and the energetics of small clusters, as well as to correctly reproduce various properties in both gas and liquid phases. Furthermore, building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗