Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLASSES”

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 55 records · Page 3

Comprehensive Total Cost of Ownership Quantification for Vehicles with Different Size Classes and Powertrains

In order to accurately compare the costs of two vehicles, the total cost of ownership (TCO) should consist of all costs related to both purchasing and operating the vehicle. This TCO analysis builds on previous work to provide a comprehensive perspective of all relevant vehicle costs of ownership. In this report, we present what we believe to be the most comprehensive explicit financial analysis of the costs that will be incurred by a vehicle owner. This study considers vehicle cost and depreciation, financing, fuel costs, insurance costs, maintenance and repair costs, taxes and fees, and other operational costs to formulate a holistic total cost of ownership and operation of multiple different vehicles. For each of these cost parameters that together constitute a comprehensive TCO, extensive literature review and data analysis were performed to find representative values in order to build a holistic TCO for vehicles of all size classes. The light- and heavy-duty vehicles selected for analysis in this report are representative of those that are on the road today and expected to be available in the future. Important additive analyses in this study include systematic analysis of vehicle depreciation, in-depth examination of insurance premium costs, comprehensive maintenance and repair estimates, analysis of all relevant taxes and fees, and considerations of specific costs applicable to commercial vehicles. We find that cars depreciate faster than light trucks and that older plug-in electric vehicles have a greater depreciation rate than newer electric vehicles. Light-duty vehicle (LDV) insurance costs show comparable costs for different powertrains, and lower costs for larger size classes. Medium- and heavy-duty vehicle (MHDV) insurance costs vary significantly by vocation. Electric and electrified powertrains have lower maintenance and repair costs than internal combustion engine (ICE) powertrains for all vehicle sizes, relative to vehicle price. MHDV maintenance and repair costs depend heavily on vocation and duty cycle. LDV taxes and fees are comparable across powertrain types and size classes, though marginally higher registration fees exist for alternative fuel vehicles. MHDV fees depend on the vocation, weight rating, and state. Many electric tractor trailers would be affected by additional battery weight, reducing the available payload capacity, and this cost can be substantial. Electric vehicle charging for commercial vehicles can be time-consuming; labor rates can cause this cost to dominate TCO. With improved knowledge of each of the cost components, we calculate a lifetime TCO for comparison across vehicles of different types and attributes. For a simulated small sport utility vehicle in 2025, modeled using Autonomie, the hybrid electric vehicle (HEV) has the lowest cost, followed by the conventional ICE vehicle. For MHDV, TCO can be drastically different depending on the vocation. Long-haul vehicles typically have the lowest per-mile costs. Excluding labor costs, the class 4 delivery has a comparable TCO to the day cab. Vocational trucks, refuse trucks, and transit buses have a high per-mile cost of ownership due to maintenance and insurance. For all of these vehicles, the cost of operating the vehicle is heavily weighted by the labor of the driver, followed by the fuel costs. While the HEV begins as the lowest cost powertrain for passenger vehicles, fuel cells are forecast to reach cost parity by 2030 when hydrogen prices reach $\$ 5$/kg while battery electric vehicles (BEV) reach cost parity by 2035 at a battery cost of $\$ 98$ per usable kWh of capacity, with these two technologies being the lowest cost in 2050. For the class 8 day cab tractor, the HEV and ICE vehicle begin as the lowest cost powertrains, and the 250-mile BEV reduces in cost from the most expensive to the least expensive by 2030.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Machine Learning Framework for Conotoxin Class and Molecular Target Prediction

Conotoxins are small and highly potent neurotoxic peptides derived from the venom of marine cone snails which have captured the interest of the scientific community due to their pharmacological potential. These toxins display significant sequence and structure diversity, which results in a wide range of specificities for several different ion channels and receptors. Despite the recognized importance of these compounds, our ability to determine their binding targets and toxicities remains a significant challenge. Predicting the target receptors of conotoxins, based solely on their amino acid sequence, remains a challenge due to the intricate relationships between structure, function, target specificity, and the significant conformational heterogeneity observed in conotoxins with the same primary sequence. We have previously demonstrated that the inclusion of post-translational modifications, collisional cross sections values, and other structural features, when added to the standard primary sequence features, improves the prediction accuracy of conotoxins against non-toxic and other toxic peptides across varied datasets and several different commonly used machine learning classifiers. Here, we present the effects of these features on conotoxin class and molecular target predictions, in particular, predicting conotoxins that bind to nicotinic acetylcholine receptors (nAChRs). We also demonstrate the use of the Synthetic Minority Oversampling Technique (SMOTE)-Tomek in balancing the datasets while simultaneously making the different classes more distinct by reducing the number of ambiguous samples which nearly overlap between the classes. In predicting the alpha, mu, and omega conotoxin classes, the SMOTE-Tomek PCA PLR model, using the combination of the SS and P feature sets establishes the best performance with an overall accuracy (OA) of 95.95%, with an average accuracy (AA) of 93.04%, and an f1 score of 0.959. Using this model, we obtained sensitivities of 98.98%, 89.66%, and 90.48% when predicting alpha, mu, and omega conotoxin classes, respectively. Similarly, in predicting conotoxins that bind to nAChRs, the SMOTE-Tomek PCA SVM model, which used the collisional cross sections (CCSs) and the P feature sets, demonstrated the highest performance with 91.3% OA, 91.32% AA, and an f1 score of 0.9131. The sensitivity when predicting conotoxins that bind to nAChRs is 91.46% with a 91.18% sensitivity when predicting conotoxins that do not bind to nAChRs.

59 BASIC BIOLOGICAL SCIENCES↗

Class imbalance in out-of-distribution datasets: Improving the robustness of the TextCNN for the classification of rare cancer types

In the last decade, the widespread adoption of electronic health record documentation has created huge opportunities for information mining. Natural language processing (NLP) techniques using machine and deep learning are becoming increasingly widespread for information extraction tasks from unstructured clinical notes. Disparities in performance when deploying machine learning models in the real world have recently received considerable attention. In the clinical NLP domain, the robustness of convolutional neural networks (CNNs) for classifying cancer pathology reports under natural distribution shifts remains understudied. In this research, we aim to quantify and improve the performance of the CNN for text classification on out-of-distribution (OOD) datasets resulting from the natural evolution of clinical text in pathology reports. We identified class imbalance due to different prevalence of cancer types as one of the sources of performance drop and analyzed the impact of previous methods for addressing class imbalance when deploying models in real-world domains. Our results show that our novel class-specialized ensemble technique outperforms other methods for the classification of rare cancer types in terms of macro F1 scores. We also found that traditional ensemble methods perform better in top classes, leading to higher micro F1 scores. Based on our findings, we formulate a series of recommendations for other ML practitioners on how to build robust models with extremely imbalanced datasets in biomedical NLP applications.

60 APPLIED LIFE SCIENCES↗

FY23 status report on the development of new ASME Section III, Division 5 Class B rules

This report summarizes the work done in Fiscal Year 2023 on the development of the new American Society of Mechanical Engineers Boiler and Pressure Vessel Code, Section III, Division 5, Class B rules to address the gaps identified for high temperature reactor designs. The summary of the design-by-analysis strain limit evaluation and creep-fatigue damage assessment is presented. The proposed design-by-analysis creep-fatigue damage calculation approach uses a new elastic follow-up-based Isochronous Stress Strain Curve stress relaxation procedure. This approach captures the elastic follow-up generated due to interactions of components with adjacent components, supports, and other connections in the power plant. A set of sample problems are selected to validate the proposed design-by-analysis rules for creep-fatigue damage assessment. The proposed Class B rules are evaluated against the Class A elastic design rules, and the experimental data obtained from a family of Simplified Model Test based key-feature test results. The proposed Class B creep-fatigue damage assessment methodology yields conservative design cycles estimates compared to the experimental results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interim Development of new Class B Code Case with Variable Design Lifetime and enhanced design rules to guard against cyclic structural failure mode

This report summarizes updates on the ongoing development of the new American Society of Mechanical Engineers Section III, Division 5, Class B rules to address the potential failure modes of Class B components. An initial design-by-analysis approach and supporting design rules are presented. The elastic-perfectly plastic (EPP) analysis approach has been adopted for primary stress limit and ratcheting check to prevent the component failure against the structural failure mode against primary load and strain accumulation due to cyclic loads. The creep-fatigue damage assessment uses explicitly defined elastic follow-up calculated from the stress concentration region of given component. Damage fraction calculation uses a novel coupled approach to capture the influence of the elastic follow-up and provide adequate conservatism for Class B components. The proposed Class B rules does not use the stress classification approach and uses the combined loads to assess component.

36 MATERIALS SCIENCE↗

Factorized class $\mathrm{S}$ theories and surface defects

It is known that some theories of class S are actually factorized into multiple decoupled nontrivial four-dimensional $\mathcal{N}$ = 2 It is known that some theories of class S are actually factorized into multiple decoupled nontrivial four-dimensional $\mathcal{N}$ = 2 superconformal SU(2) QCD, describing this factorized theory as a class S theory of type A 3 on a five-punctured sphere with a twist line. Separately, we also present explicit checks that the Coulomb branch of a putative factorized class S theory has the expected product structure, in two examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of aquatic hazard concentrations for the oil spill response product class of surface washing agents using species sensitivity distributions

Surface washing agents (SWAs) are a diverse class of oil spill response products intended to facilitate removal of stranded oil from shorelines. This class of agents has high application rates relative to other categories of spill response products, but global toxicity data is generally limited to two standard test species: inland silverside and mysid shrimp. Here, in this study, we provide a framework to maximize the utility of limited toxicity data across a product class. To characterize species sensitivity to SWAs, the toxicity of three agents spanning a range of chemical and physical properties were tested in eight species. The relative sensitivity of mysids shrimp and inland silversides as surrogate test organisms was determined. Toxicity normalized species sensitivity distributions (SSDn) were used to estimate fifth centile hazard concentration (HC5) values for SWAs with limited toxicity data. Chemical toxicity distributions (CTD) of SWA HC5 values were used to compute a fifth centile chemical hazard distribution (HD5) to provide a more comprehensive assessment of hazard across a spill response product class with limited toxicity data than traditional single species or single agent approaches can give.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Identification of $\mathrm{IOMA}$-class neutralizing antibodies targeting the $\mathrm{CD4}$-binding site on the $\mathrm{HIV}$-1 envelope glycoprotein

A major goal of current HIV-1 vaccine design efforts is to induce broadly neutralizing antibodies (bNAbs). The VH1-2-derived bNAb IOMA directed to the CD4-binding site of the HIV-1 envelope glycoprotein is of interest because, unlike the better-known VH1-2-derived VRC01-class bNAbs, it does not require a rare short light chain complementarity-determining region 3 (CDRL3). Here, we describe three IOMA-class NAbs, ACS101-103, with up to 37% breadth, that share many characteristics with IOMA, including an average-length CDRL3. Cryo-electron microscopy revealed that ACS101 shares interactions with those observed with other VH1-2 and VH1-46-class bNAbs, but exhibits a unique binding mode to residues in loop D. Analysis of longitudinal sequences from the patient suggests that a transmitter/founder-virus lacking the N276 glycan might have initiated the development of these NAbs. Together these data strengthen the rationale for germline-targeting vaccination strategies to induce IOMA-class bNAbs and provide a wealth of sequence and structural information to support such strategies.

59 BASIC BIOLOGICAL SCIENCES↗

Structural insights into VRC01-class bnAb precursors with diverse light chains elicited in the IAVI G001 human vaccine trial

The development of germline-targeting vaccines represents a potentially transformative strategy to elicit broadly neutralizing antibodies (bnAbs) against HIV and other antigenically diverse pathogens. Here, we report on structural characterization of vaccine-elicited VRC01-class bnAb precursors in the IAVI G001 Phase 1 clinical trial with the eOD-GT8 60mer nanoparticle as immunogen. High-resolution X-ray structures of eOD-GT8 monomer complexed with Fabs of five VRC01-class bnAb precursors with >90% germline identity revealed a conserved mode of binding to the HIV CD4-binding site via IGHV1-2-encoded heavy chains, mirroring mature bnAb interactions. The light-chain V-gene diversity emulated VRC01 bnAbs and stabilized antigen engagement, while their conserved five-residue LCDR3 motifs prevented steric clashes. Notably, the VRC01-class bnAb precursors accommodated the N276 glycan, a key barrier in HIV Env recognition, through structural rearrangements in HCDR3 or LCDR1, despite its absence in the immunogen. Surface plasmon resonance analysis showed that 87% of elicited antibodies retained glycan binding capacity, albeit with reduced affinity. These findings validate the ability of eOD-GT8 60mer nanoparticles to prime VRC01-class bnAb precursors with native-like paratopes but with intrinsic glycan adaptability. Structural mimicry of mature bnAbs was observed even with limited somatic hypermutation, indicating that critical features are encoded in the germline repertoire. The structures highlight how germline-encoded features drive bnAb-like recognition at early stages. This work provides molecular evidence supporting germline targeting in humans and provides guidance for designing booster immunogens to shepherd affinity maturation toward broad neutralization.

Science & Technology - Other Topics↗

hi_class background evolution, initial conditions and approximation schemes

Cosmological datasets have great potential to elucidate the nature of dark energy and test gravity on the largest scales available to observation. Theoretical predictions can be computed with hi_class (www.hiclass-code.net), an accurate, fast and flexible code for linear cosmology, incorporating a wide range of dark energy theories and modifications to general relativity. We introduce three new functionalities into hi_class: (1) Support for models based on covariant Lagrangians, including a constraint-preserving integration scheme for the background evolution and a series of worked-out examples: Galileon, nKGB, quintessence (monomial, tracker) and Brans-Dicke. (2) Consistent initial conditions for the scalar-field perturbations in the deep radiation era, identifying the conditions under which modified-gravity isocurvature perturbations may grow faster than adiabatic modes leading to a loss of predictivity. (3) An automated quasistatic approximation scheme allowing order-of-magnitude improvement in computing performance without sacrificing accuracy for wide classes of models. These enhancements bring the treatment of dark energy and modified gravity models to the level of detail comparable to software tools restricted to standard $\Lambda$CDM cosmologies. Finally, the hi_class code is publicly available (https://github.com/miguelzuma/hi_class_public), ready to explore current data and prepare for next-generation experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhancing dimensionality prediction in hybrid metal halides via feature engineering and class-imbalance mitigation

We present a machine learning (ML) framework for predicting the structural dimensionality of hybrid metal halides (HMHs), including organic-inorganic perovskites, using a combination of chemically-informed feature engineering and advanced class-imbalance handling techniques. This study is motivated by the small and highly imbalanced nature of experimentally available HMH datasets, which limits the applicability and reliability of conventional ML approaches. The dataset, consisting of 494 HMH structures, is highly imbalanced across dimensionality classes (0D, 1D, 2D, 3D), posing significant challenges to predictive modeling. To mitigate this limitation, the dataset was augmented to 1336 samples using the synthetic minority oversampling technique, enabling improved learning of underrepresented dimensionality classes while preserving chemically meaningful feature relationships. We developed interaction-based descriptors designed to capture coupled steric and polarity effects relevant to dimensionality prediction, which are not readily captured by standard single-parameter or composition-only descriptors. These descriptors are integrated into a multi-stage workflow combining feature selection, ensemble stacking, and performance optimization. Our approach significantly improves F1-scores for underrepresented classes, achieving robust cross-validation performance across all dimensionalities. This work demonstrates a generalizable strategy for extracting reliable and interpretable structure–dimensionality relationships from limited experimental data, enabling pre-synthesis screening of organic cations and providing a practical blueprint for small-data ML in hybrid materials systems.

36 MATERIALS SCIENCE↗

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Design by analysis rules for ASME Section III, Division 5, Class B components

The current rules for elevated temperature ASME Section III, Division 5, Class B components, other than piping, were basically adapted from the Design-by-Rule approach for Section VIII, Division 1 vessels. A goal of the proposed new rules is to explicitly account for the design life and cyclic service while recognizing the less-rigorous requirements commensurate with lesser safety consideration. Further goals are to maximize the use of modern computational technology, for example finite element analysis in conjunction with reference stress concepts, and to avoid the need for stress classification. This paper summarizes the work done to address these issues. A summary of the design-by-analysis for primary loads, strain limit evaluation, and creep-fatigue damage assessment is presented. The proposed design-by-analysis creep-fatigue damage calculation approach uses a new elastic follow-up-based Isochronous Stress Strain Curve stress relaxation procedure. A set of sample problems are selected to validate the proposed design-by-analysis rules for creep-fatigue damage assessment. The proposed Class B rules are evaluated against the Class A elastic design rules and the experimental data obtained from a family of a simplified model test-based key-feature test results. The proposed Class B creep-fatigue damage assessment methodology yields conservative design cycles estimates compared to the experimental results.

36 MATERIALS SCIENCE↗

Cultivation of Methanonezhaarchaeia, the third class of methanogens within the phylum Thermoproteota

Methane is a potent greenhouse gas, largely produced by methanogenic archaea, contributing to Earth’s dynamic climate and biogeochemical cycles. In the past decade, metagenomics revealed that lineages outside of the Euryarchaeota superphylum encode genes for methanogenesis. This was recently confirmed through the cultivation of two classes of methanogenic Thermoproteota. Thus far, all methanogens within the Thermoproteota are predicted or were shown to be methylotrophic. The only exception to this are the Nezhaarchaea, for which metagenomic predictions suggest they are CO 2 -reducing methanogens. Here, we demonstrate methanogenic activity in a third class of Thermoproteota, the Methanonezhaarchaeia. Contrary to genomic predictions for this class, we cultivated a methylotrophic species, Candidatus Methanonezhaarchaeum fastidiosum YNP3N, highlighting the importance of testing metagenomic hypotheses through experimentation. We investigate the metabolic diversity of Methanonezhaarchaeia, including metabolic modifications accompanying frequent loss of methanogenesis in this class. This highlights gaps in our understanding of the biochemistry, diversity, and evolution of thermoproteotal methanogens and their contributions to carbon cycling.

Kohtz, Anthony J. [Montana State Univ., Bozeman, M↗

Specification of Radionuclide Classes for MELCOR Molten Salt Reactor Simulations

MELCOR has been used extensively to facilitate virtual investigations into severe nuclear accidents for light-water reactors (LWRs). Non-light water reactors (non-LWRs) render some LWR-centric approaches potentially unsuitable. MELCOR has been instrumental in analyzing source terms for LWRs and has recently expanded its applicability to non-LWRs. To simplify radionuclide (RN) tracking, MELCOR currently groups elements into 17 classes, each containing representative species. This grouping, optimized for LWRs, is not appropriate for non-LWRs due to the different chemistry. This necessitates a reevaluation of radionuclide transport modeling. This report introduces a new class scheme for MELCOR tailored to MSR modeling, expanding the current 17 classes to 32 and are explained in the context of a UF 4 fueled FLiBe carrier MSR. It provides a discussion and justification for the new groupings and outlines a methodology for discovering and defining additional classes in MELCOR using a sample calculated RN inventory and a Gibbs energy minimizer (GEM).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Allometric-kinetic model predictions of concentration ratios for anthropogenic radionuclides across animal classes and food selection

Protection of the environment from radiation fundamentally relies on dose assessments for non-human biota. These dose assessments use measured or predicted concentrations of radionuclides in soil or water combined with Concentration Ratios (CRs) to estimate tissue concentrations in animals and plants, yet there is a paucity of CR data relative to the vast number of potential taxa and radioactive contaminants in the environment and their taxon-specific ecosystems. Because there are many taxa each having very different behaviors and biology and many possible bioavailable radionuclides, CRs have the potential to vary by orders-of-magnitude, as often seen in published data. Given the diversity of taxa, the International Commission on Radiological Protection (ICRP) has categorized non-human biota into 12 representative organisms (ROs), while the U.S. Department of Energy (DOE) uses the non-taxon specific categories of terrestrial, riparian, and aquatic animals. The question we examine here is: do these categorizations adequately “represent” the broad diversity of animals? To explore this question, here we analyze an Allometric-Kinetic (A-K) model to calculate radionuclide-specific CRs for common animal classes which are then further subcategorized into herbivores, omnivores, carnivores, and invertebrate detritivores. Comparisons in CRs among animal classes exhibited only small differences, but there was order of magnitude differences between herbivores, carnivores, and especially detritivores, for many radionuclides of interest. These findings suggest that the ICRP ROs and the DOE categories are reasonable, but their accuracy could be improved by including sub-categories related to animal dietary ecology. Finally, comparisons of A-K model predicted CR values to published CRs show order-of-magnitude variations, providing justification for additional studies of assimilation across radionuclides and animal classes.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Visualization of conformational changes and membrane remodeling leading to genome delivery by viral class-II fusion machinery

Chikungunya virus (CHIKV) is a human pathogen that delivers its genome to the host cell cytoplasm through endocytic low pH-activated membrane fusion mediated by class-II fusion proteins. Though structures of prefusion, icosahedral CHIKV are available, structural characterization of virion interaction with membranes has been limited. Here, we have used cryo-electron tomography to visualize CHIKV’s complete membrane fusion pathway, identifying key intermediary glycoprotein conformations coupled to membrane remodeling events. Using sub-tomogram averaging, we elucidate features of the low pH-exposed virion, nucleocapsid and full-length E1-glycoprotein’s post-fusion structure. Contrary to class-I fusion systems, CHIKV achieves membrane apposition by protrusion of extended E1-glycoprotein homotrimers into the target membrane. The fusion process also features a large hemifusion diaphragm that transitions to a wide pore for intact nucleocapsid delivery. Our analyses provide comprehensive ultrastructural insights into the class-II virus fusion system function and direct mechanistic characterization of the fundamental process of protein-mediated membrane fusion.

59 BASIC BIOLOGICAL SCIENCES↗

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗