Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural settings”

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 1,297 records · Page 72

Nanoscale structural evolution and phase transformation of geopolymers: In situ SAXS/WAXS investigation under uniaxial tension at elevated temperatures

This investigation delves into the degradation mechanisms of high-density polyethylene geomembranes (PE GMXs) under a spectrum of conditions, replicating real-world scenarios within a rigorously controlled laboratory setting. The treatment protocols applied induced a notable increase in the crystallinity of the treated specimens relative to the untreated controls. Small-angle x-ray scattering (SAXS) analysis identified an initial long period (interlamellar distance) of 16.9 nm for the untreated polymer, which expanded by 19.5% at a strain of 16.7 percent. Conversely, the treated PE GMXs exhibited a more gradual elongation of the long period, with an increase of merely 10.6% at a strain of 23.3 percent. At an elevated temperature of 65°C, both samples exhibited pronounced strain hardening, with the treated PE GMXs demonstrating superior stability even at a strain of 150 percent. Wide-angle x-ray scattering (WAXS) experiments corroborated these observations, revealing that the diffraction patterns of the untreated PE remained stable up to a strain of 16.7%, whereas those of the treated PE remained distinct up to a strain of 46.1 percent. Scanning electron microscopy (SEM) images substantiated the formation of a shish–kebab structure in the treated samples. The study concludes that the geopolymer underwent oxidation and material degradation as a result of the chemical and mechanical treatments, transitioning to a more crystalline state and concomitantly losing its initial elasticity.

36 MATERIALS SCIENCE↗

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗

An active learning high-throughput microstructure calibration framework for solving inverse structure–process problems in materials informatics

Determining a process–structure–property relationship is the holy grail of materials science, where both computational prediction in the forward direction and materials design in the inverse direction are essential. Problems in materials design are often considered in the context of process–property linkage by bypassing the materials structure, or in the context of structure–property linkage as in microstructure-sensitive design problems. However, there is a lack of research effort in studying materials design problems in the context of process–structure linkage, which has a great implication in reverse engineering. In this paper, given a target microstructure, we propose an active learning high-throughput microstructure calibration framework to derive a set of processing parameters, which can produce an optimal microstructure that is statistically equivalent to the target microstructure. The proposed framework is formulated as a noisy multi-objective optimization problem, where each objective function measures a deterministic or statistical difference of the same microstructure descriptor between a candidate microstructure and a target microstructure. Furthermore, to significantly reduce the physical waiting wall-time, we enable the high-throughput feature of the microstructure calibration framework by adopting an asynchronously parallel Bayesian optimization to exploit high-performance computing resources. Case studies in additive manufacturing and grain growth are used to demonstrate the applicability of the proposed framework, where kinetic Monte Carlo (kMC) simulation is used as a forward predictive model, such that for a given target microstructure, the target processing parameters that produced this microstructure are successfully recovered.

36 MATERIALS SCIENCE↗

Reducing Qubit Requirements while Maintaining Numerical Precision for the Variational Quantum Eigensolver: A Basis-Set-Free Approach

We present a basis-set-free approach to the variational quantum eigensolver using an adaptive representation of the spatial part of molecular wave functions.Our approach directly determines system-specific representations of qubit Hamiltonians while fully omitting globally defined basis sets. In this work, we use directly determined pair-natural orbitals on the level of second-order perturbation theory. This results in compact qubit Hamiltonians with high numerical accuracy. We demonstrate initial applications with compact Hamiltonians on up to 22 qubits where conventional representation would for the same systems require 40–100 or more qubits. Here, we further demonstrate reductions in the quantum circuits through the structure ofthe pair-natural orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The gravitational lensing imprints of DES Y3 superstructures on the CMB: a matched filtering approach

Low-density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence |$\kappa$|⁠. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 (DES Y3) data set, covering approximately 4200 deg|$^2$| of the sky. We use two distinct void-finding algorithms: a 2D void-finder that operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the Marenostrum Institut de Ciències de l’Espai N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure |$A_\kappa$|⁠, the amplitude of the observed lensing signal relative to the simulation template, obtaining |$A_\kappa = 1.03 \pm 0.22$| (⁠|$4.6\sigma$| significance) for Voxel and |$A_\kappa = 1.02 \pm 0.17$| (⁠|$5.9\sigma$| significance) for 2D voids, both consistent with Lambda cold dark matter expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure |$A_\kappa = 0.87 \pm 0.15$| (⁠|$5.9\sigma$| significance). The leading source of noise in our measurements is Planck noise, implying that data from the Atacama Cosmology Telescope, South Pole Telescope and CMB-S4 will increase sensitivity and allow for more precise measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗

Dynamical Classification of Trans-Neptunian Objects Detected by the Dark Energy Survey

The outer solar system contains a large number of small bodies (known as trans-Neptunian objects or TNOs) that exhibit diverse types of dynamical behavior. The classification of bodies in this distant region into dynamical classes—subpopulations that experience similar orbital evolution—aids in our understanding of the structure and formation of the solar system. In this work, we propose an updated dynamical classification scheme for the outer solar system. This approach includes the construction of a new (automated) method for identifying mean motion resonances. We apply this algorithm to the current data set of TNOs observed by the Dark Energy Survey (DES) and present a working classification for all of the DES TNOs detected to date. Our classification scheme yields 1 inner centaur, 19 outer centaurs, 21 scattering disk objects, 47 detached TNOs, 48 securely resonant objects, 7 resonant candidates, and 97 classical belt objects. Among the scattering and detached objects, we detect 8 TNOs with semimajor axes greater than 150 au.

79 ASTRONOMY AND ASTROPHYSICS↗

Dispersion curve regression – assisted wideband local wavenumber analysis for characterizing three-dimensional (3D) profile of hidden corrosion damage

Corrosion as common damage in civil, petrochemical, nuclear, and aerospace structures affects the integrity and safety of the structures and may lead to catastrophic failures. This paper presents dispersion curve regression-assisted local wavenumber analysis method, which can analyze the time–space wavefield containing wideband information of wave-damage interaction and further extract the structural information carried by such wavefield for characterizing hidden corrosion damage in an isotropic plate. To acquire the time–space wavefield, a noncontact scanning laser Doppler vibrometer is used. In our analysis method, the acquired time–space wavefield is firstly processed to generate a local wavenumber function. The derived wavenumber function is further analyzed with a dispersion curve regression step, which searches in a set of theoretical frequency-wavenumber dispersion curves for different plate thicknesses and identifies the optimal plate thickness whose theoretical curve best matches the frequency-wavenumber relation contained in the local wavenumber function. By this means, a 3D profile including both in plane and thickness dimensions of the structure can be constructed for corrosion visualization and quantification. The experimental study demonstrates that our method can quantify the profile of hidden uniform corrosion in metal plates in terms of its in-plane shape and size as well as its out-of-plane depth in the subwavelength scale. Moreover, compared to previous Lamb wave-based corrosion inspection methods, our technique allows for noncontact 3D characterization of hidden corrosion from a far distance to the structure. Finally, we believe this work will inspire new 3D damage quantification methods that are based on wavefield analysis as well as enable potential applications for the quantification of hidden corrosion in civil, petrochemical, nuclear, and aerospace structures.

42 ENGINEERING↗

The Synthescope: A Vision for Combining Synthesis with Atomic Fabrication

Here, the scanning transmission electron microscope, a workhorse instrument in materials characterization, is being transformed into an atomic-scale material-manipulation platform. With an eye on the trajectory of recent developments and the obstacles toward progress in this field, a vision for a path toward an expanded set of capabilities and applications is provided. The microscope is reconceptualized as an instrument for fabrication and synthesis with the capability to image and characterize atomic-scale structural formation as it occurs. Further development and refinement of this approach may have substantial impact on research in microelectronics, quantum information science, and catalysis, where precise control over atomic-scale structure and chemistry of a few “active sites” can have a dramatic impact on larger-scale functionality and where developing a better understanding of atomic-scale processes can help point the way to larger-scale synthesis approaches.

36 MATERIALS SCIENCE↗

Measurement of the $\textrm{t}\overline{\textrm{t}}\textrm{H}$ and tH production rates in the H → $\textrm{b}\overline{\textrm{b}}$ decay channel using proton-proton collision data at $\sqrt{s}$ = 13 TeV

An analysis of the production of a Higgs boson (H) in association with a top quark-antiquark pair ($\textrm{t}\overline{\textrm{t}}\textrm{H}$) or a single top quark (tH) is presented. The Higgs boson decay into a bottom quark-antiquark pair (H → $\textrm{b}\overline{\textrm{b}}$) is targeted, and three different final states of the top quark decays are considered, defined by the number of leptons (electrons or muons) in the event. The analysis utilises proton-proton collision data collected at the CERN LHC with the CMS experiment at $\sqrt{s}$ = 13 TeV in 2016–2018, which correspond to an integrated luminosity of 138 fb −1 . The observed $\sqrt{s}$ production rate relative to the standard model expectation is 0.33 ± 0.26 = 0.33 ± 0.17(stat) ± 0.21(syst). Additionally, the $\textrm{t}\overline{\textrm{t}}\textrm{H}$ production rate is determined in intervals of Higgs boson transverse momentum. An upper limit at 95% confidence level is set on the tH production rate of 14.6 times the standard model prediction, with an expectation of ${19.3}_{-6.0}^{+9.2}$. Finally, constraints are derived on the strength and structure of the coupling between the Higgs boson and the top quark from simultaneous extraction of the $\textrm{t}\overline{\textrm{t}}\textrm{H}$ and tH production rates, and the results are combined with those obtained in other Higgs boson decay channels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Physics-informed machine-learning model of temperature evolution under solid phase processes

We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.

36 MATERIALS SCIENCE↗

Viral Inhibitors Derived from Macroalgae, Microalgae, and Cyanobacteria: A review of antiviral potential throughout pathogenesis

Viruses are abiotic obligate parasites utilizing complex mechanisms to hijack cellular machinery and reproduce, causing multiple harmful effects in the process. Viruses represent a growing global health concern; at the time of writing, COVID-19 has killed at least two million people around the world and devastated global economies. Lingering concern regarding the virus’ prevalence yet hampers return to normalcy. While catastrophic in and of itself, COVID-19 further heralds in a new era of human-disease interaction characterized by the emergence of novel viruses from natural sources with heretofore unseen frequency. Due to deforestation, population growth, and climate change, we are encountering more viruses that can infect larger groups of people with greater ease and increasingly severe outcomes. The devastation of COVID-19 and forecasts of future human/disease interactions call for a creative reconsideration of global response to infectious disease. There is an urgent need for accessible, cost-effective antiviral (AV) drugs that can be mass-produced and widely distributed to large populations. Development of AV drugs should be informed by a thorough understanding of viral structure and function as well as human biology. To maximize efficacy, minimize cost, and reduce development of drug-resistance, these drugs would ideally operate through a varied set of mechanisms at multiple stages throughout the course of infection. Due to their abundance and diversity, natural compounds are ideal for such comprehensive therapeutic interventions. Promising sources of such drugs are found throughout nature; especially remarkable are the algae, a polyphyletic grouping of phototrophs that produce diverse bioactive compounds. While not much literature has been published on the subject, studies have shown that these compounds exert antiviral effects at different stages of viral pathogenesis. In this review, we follow the course of viral infection in the human body and evaluate the AV effects of algae-derived compounds at each stage. Specifically, we examine the AV activities of algae-derived compounds at the entry of viruses into the body, transport through the body via the lymph and blood, infection of target cells, and immune response. We discuss what is known about algae-derived compounds that may interfere with the infection pathways of SARS-CoV-2; and review which algae are promising sources for AV agents or AV precursors that, with further investigation, may yield life-saving drugs due to their diversity of mechanisms and exceptional pharmaceutical potential.

Algae, antiviral, SARS-CoV-2, seaweed, cyanobacter↗

LAVA 1.0: A general-purpose python toolkit for calculation of material properties with LAMMPS and VASP

Here, we introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion. These functionalities are demonstrated for different interatomic potentials and DFT calculations in Al.

36 MATERIALS SCIENCE↗

A risk-informed assessment tool to support review of license amendment requests for spent nuclear fuel dry cask storage systems

We report that while risk-informed safety evaluations are common methodologies used by the U.S. Nuclear Regulatory Commission (NRC), their application in the field of commercial spent nuclear fuel (SNF) dry storage has received limited attention. For instance, license amendment requests (LARs) for SNF dry cask storage systems (DCSSs) are evaluated according to U.S. NRC standard review plans (SRPs). However, risk-informed LAR review strategies could improve the understanding of critical system modifications, and lead to more predictable, efficient, and consistent LAR review processes. A methodology to develop a risk-informed LAR review tool is presented. Although focused on DCSSs, this methodology can be used to develop similar tools for a broad spectrum of safety-relevant systems. The end product (i.e., the tool) includes a tree diagram to help visualize the review risk significances assigned to a predefined set of system modifications and supporting rationale documentation. Instructions support correct tool application and provide the user with a path to incorporate currently unevaluated modifications in the tool structure. The demonstration of the methodology leads to the conclusion that risk-informed reviews of DCSS LARs are possible, although the tool precision could benefit from additional risk information. Further, the results indicate that many typical LARs involve low-risk or medium-risk modifications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗

Broadening the scope of high structural dimensionality nanomaterials using pyridine-based curcuminoids

We present a new heteroditopic ligand (3pyCCMoid) that contains the typical skeleton of a curcuminoid (CCMoid) decorated with two 3-pyridyl groups. The coordination of 3pyCCMoid with ZnII centres results in a set of novel coordination polymers (CPs) that display different architectures and dimensionalities (from 1D to 3D). Our work analyses how synthetic methods and slight changes in the reaction conditions affect the formation of the final materials. Great efforts have been devoted toward understanding the coordination entities that provide high dimensional systems, with emphasis on the characterization of 2D materials, including analyses of different types of substrates, stability and exfoliation in water. Here, we foresee the great use of CCMoids in the field of CPs and emphasize 3pyCCMoid as a new-born linker.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗