SEARCH · Engineering Papers

Results for “Reproducibility”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Integral and Separate Effects Test Facilities To Support Water Cooled Small Modular Reactors: A Review

This study reviews previous experimental facilities and test programs relevant to water-cooled reactor system design and analysis to meet regulatory compliances. This study aims to find the best solution for designing the required experiments, obtaining necessary test data, and verifying the developed computer code/models to support the new reactor design and development while minimizing cost and time while leveraging experiences from previous facilities to minimize. Nuclear reactor licensing requires supportive design, analysis, and experimental results to ensure the safety of the full-scale prototype reactor in regular operation, as well as during postulated accident scenarios. These reactor design analyses are generally performed using system codes and other associated simulation tools that require assessment, verification and validation using an appropriate experimental dataset. Experimental facilities used for reactor system safety analysis and system code assessments are categorized into integral effects test (IET) and separate effects test (SET) facilities. Further, these IET and SET experiments and studies use geometrically scaled systems to reproduce the prototype system behavior at a reasonable cost, albeit with some scaling-related distortions. The design challenge of these model facilities is to identify and minimize scaling distortions while reproducing the most important operational phenomena in steady-state operation and in postulated accident scenarios. Lessons learned from previous experimental facilities, models, and correlations can support the development of new multipurpose, scaled, hybrid, integrated, and modular experimental facilities for advanced light water-cooled small modular reactors (SMRs). Successful operation of these facilities can significantly reduce upfront reactor development and demonstration costs and time to deployment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Benchmarking the accuracy of coverage-dependent models: adsorption and desorption of benzene on Pt (111) and Pt 3 Sn (111) from first principles

Bimetallic catalysts have demonstrated properties favorable for upgrading biofuel through catalytic hydrodeoxygenation. However, the design and optimization of such bimetallic catalysts requires the ability to construct accurate, predictive models of these systems. To generate a model that predicts the kinetic behavior of benzene adsorbed on Pt (111) and a Pt 3 Sn (111) surface alloy (Pt 3 Sn (111)), the adsorption of benzene was studied for a wide range of benzene coverages on both surfaces using density functional theory (DFT) calculations. The adsorption energy of benzene was found to correlate linearly with benzene coverage on Pt (111) and Pt 3 Sn (111); both surfaces exhibited net repulsive lateral interactions. Through an analysis of the d-band properties of the metal surface, it was determined that the coverage dependence is a consequence of the electronic interactions between benzene and the surface. The linear coverage dependence of the adsorption energy allowed us to quantify the influence of the lateral interactions on the heat of adsorption and temperature programmed desorption (TPD) spectra using a mean-field model. A comparison of our simulated TPD to experiment showed that this mean-field model adequately reproduces the desorption behavior of benzene on Pt (111) and Pt 3 Sn (111). In particular, the TPD correctly exhibits a broadening desorption peak as the initial coverage of benzene increases on Pt (111) as well as and a low temperature desorption peak on Pt 3 Sn (111). Furthermore, due to the sensitivity of the TPD peak temperature to the desorption energy, precise alignment of experimental and theoretical TPD spectra demands an accurate calculation of the adsorption energy. Therefore, an analysis of the effect of the exchange-correlation functional on TPD modeling is presented. Through this work, we show the necessity of incorporating lateral interactions into theoretical models in order to correctly predict experimental behavior.

09 BIOMASS FUELS↗

DOE OSTI2020

Microbial potentiometric sensor: A new approach to longstanding challenges

The underlying hypothesis of this study is that simple potentiometric measurements between sensing electrodes and a shared reference electrode - Microbial Potentiometric Sesnor (MPS) system - can be employed in a long-term, continuous mode of operation to resolve the spatial and temporal changes in environmental systems. To address the hypothesis, (1) a conceptual description of the MPS system and its postulated principle of operation are provided; (2) the MPS system performance is documented under controlled laboratory conditions; and (3) the capabilities of the MPS system are documented under quiescent and dynamic field condition. In a laboratory setting, the variability among different MPS signals was insignificant confirming the postulated high accuracy and reproducibility of the sensor performance. It also demonstrated statistically significant correlations with dissolved oxygen (DO) and oxidation-reduction potential (ORP) sensors, while showing capabilities of operating under anoxic and anaerobic conditions. Regardless of their locations in the model wetland system, three MPS sensors functioned without interruption and cleaning for a period >2 years, and thus demonstrating long-term durability of the MPS technology. In real batch-wastewater treatment facility, the deployed MPS system signals were able to describe the organic carbon trends and correlate with each treatment phase in a cycle. Data reproducibility and reliability exceeded the expectations better describing the carbon treatment levels than the DO and ORP sensors (p < 4.4 × 10 –162 vs phase adjusted p < 3.0 × 10 –58 ). While MPS signals correlate with specific parameters that describe the local process or environments, it is more prudent to employ both the magnitude and pattern of a composite signal like the MPS signal describe the change to reflect any shift in the local environment. Finally, when compared to a baseline pattern, this change in signal magnitude and pattern reveals important information that can be employed to tailor and optimize any condition or process which involves microorganisms.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Formation trajectories of solution-processed perovskite thin films from mixed solvents

The engineering of mixed-solvent formulations and their evaporation conditions are key to reproducible perovskite coatings for high-performance photovoltaics. Here, we report a lumped-parameter evaporation model to predict the evolution of a perovskite ink liquid film over time (solvent ratio, solute concentration, and film thickness). The drying-rate model is validated via in situ film-thickness measurements, and the predicted transient liquid film state is mapped as a process path. These methods allow for the prediction of process sensitivity to local environmental factors and the understanding and visualization of a broader processing parameter space enabled through the coupling of process and ink engineering. Process maps are applied to create a new framework for scalable perovskite coating development with a goal of improving the reproducibility and transferability of perovskite fabrication. This approach is demonstrated with blade-coated FA 0.83 Cs 0.17 PbI 3 photovoltaic devices, improving the photovoltaic conversion efficiency from 17.5% ± 1.7% to 20.3% ± 0.6%.

14 SOLAR ENERGY↗

DOE OSTI2021

Radial density profile and stability of capillary discharge plasma waveguides of lengths up to 40 cm

We measured the parameter reproducibility and radial electron density profile of capillary discharge waveguides with diameters of 650 μm to 2 mm and lengths of 9 to 40 cm. To the best of the authors’ knowledge, 40 cm is the longest discharge capillary plasma waveguide to date. This length is important for ≥10 GeV electron energy gain in a single laser-driven plasma wakefield acceleration stage. Evaluation of waveguide parameter variations showed that their focusing strength was stable and reproducible to <0.2% and their average on-axis plasma electron density to <1%. These variations explain only a small fraction of laser-driven plasma wakefield acceleration electron bunch variations observed in experiments to date. Measurements of laser pulse centroid oscillations revealed that the radial channel profile rises faster than parabolic and is in excellent agreement with magnetohydrodynamic simulation results. We show that the effects of non-parabolic contributions on Gaussian pulse propagation were negligible when the pulse was approximately matched to the channel. However, they affected pulse propagation for a non-matched configuration in which the waveguide was used as a plasma telescope to change the focused laser pulse spot size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Automated Coupling of Nanodroplet Sample Preparation with Liquid Chromatography–Mass Spectrometry for High-Throughput Single-Cell Proteomics

Single-cell proteomics can provide critical biological insight into the cellular heterogeneity that is masked by bulk-scale analysis. Here, we have developed a nanoPOTS (nanodroplet processing in one pot for trace samples) platform and demonstrated its broad applicability for single-cell proteomics. However, because of nanoliter-scale sample volumes, the nanoPOTS platform is not compatible with automated LC-MS systems, which significantly limits sample throughput and robustness. To address this challenge, we have developed a nanoPOTS autosampler allowing fully automated sample injection from nanowells to LC-MS systems. We also developed a sample drying, extraction, and loading workflow to enable reproducible and reliable sample injection. The sequential analysis of 20 samples containing 10 ng tryptic peptides demonstrated high reproducibility with correlation coefficients of >0.995 between any two samples. The nanoPOTS autosampler can provide analysis throughput of 9.6, 16, and 24 single cells per day using 120, 60, and 30 min LC gradients, respectively. As a demonstration for single-cell proteomics, the autosampler was first applied to profiling protein expression in single MCF10A cells using a label-free approach. At a throughput of 24 single cells per day, an average of 256 proteins was identified from each cell and the number was increased to 731 when the Match Between Runs algorithm of MaxQuant was used. Using a multiplexed isobaric labeling approach (TMT-11plex), ~77 single cells could be analyzed per day. We analyzed 152 cells from three acute myeloid leukemia cell lines, resulting in a total of 2558 identified proteins with 1465 proteins quantifiable (70% valid values) across the 152 cells. These data showed quantitative single-cell proteomics can cluster cells to distinct groups and reveal functionally distinct differences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Combining Mass Spectrometry of Picoliter Samples with a Multicompartment Electrodynamic Trap for Probing the Chemistry of Droplet Arrays

Single droplet levitation provides contactless access to the microphysical and chemical properties of micrometer-sized samples. Most applications of droplet levitation to chemical and biological systems use nondestructive optical techniques to probe droplet properties. To provide improved chemical specificity, we coupled a multicompartment quadrupole electrodynamic trap (QET) with single droplet mass spectrometry. Our QET continuously traps a monodisperse droplet population (tens to hundreds of droplets) and allows for the simultaneous sizing of a single droplet using its Mie scattering pattern. Single droplets are subsequently ejected into the ionization region of an ambient pressure inlet mass spectrometer. We optimized two complementary soft ionization techniques for picoliter aqueous droplets: (1) paper spray (PS) ionization and (2) thermal desorption glow discharge (TDGD) ionization. Both techniques detect oxygenated organic acids in single droplets, with signal-to-noise ratios >100 and detection limits on the order of 10 pg. Sensitivity and reproducibility across single droplets are driven by the droplet deposition location and spray stability in PS-MS and the ionization region humidity and analyte evaporation rate in TDGD-MS. Importantly, the analyte evaporation rate can control the TDGD-MS quantitative capability because high evaporation rates result in significant ion suppression. Additionally, this effect is mitigated by optimizing the vaporization temperature, droplet size range, and analyte volatility. We demonstrate quantitative and reproducible measurements with a droplet internal standard (<10% RSD) and compare the sensitivity of PS-MS and TDGD-MS. Finally, we demonstrate the application of QET-MS to the study of heterogeneous chemical kinetics with the reaction of gas phase O 3 and aqueous maleic acid droplets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effect of the Precursor Chemistry on the Crystallization of Triple Cation Mixed Halide Perovskites

Triple cation, mixed halide perovskite compositions have been reported to be more thermally stable, exhibit fewer phase impurities, and show higher power conversion efficiency and better reproducibility than single cation perovskites. In this work, we explain the formation of Cs0.05FA0.81MA0.14Pb(I0.85Br0.15)3 via a multimodal in situ study combining structural information from synchrotron grazing-incidence wide-angle X-ray scattering (GIWAXS) and optical properties from photoluminescence (PL) spectroscopy with density functional theory calculations (DFT). The focus here is on the effects of the solvent and antisolvent during crystallization. The predominantly used solvents N,N-dimethylformamide (DMF), dimethyl sulfoxide (DMSO), and the antisolvent chlorobenzene (CB) as well as the solvent-antisolvent-precursor interactions are investigated. Given the high elemental complexity and mutual interdependencies between solvent, antisolvent, and perovskite precursors, we found significant differences in the crystallization pathways. DMF-pure precursors show the formation of the DMF-containing intermediate phase and the nucleation of compositionally distinct perovskite phases, while when DMSO is added, only crystalline α- and δ-phases were found. In addition, the presence of DMSO helps the formation of α-perovskite. Coordination energy and bond order (BO) calculations support our experimental findings. Dripping of CB induces nucleation at room temperature, slows the α-phase formation rate, and appears to reduce the nucleation radius. These findings provide novel insights into solvent, antisolvent, and perovskite precursor interactions and their formation pathways. The complexity of interactions between solvents and reagents highlights the importance of understanding these effects to further improve the reproducibility and optimize processing conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Trust Not Verify? The Critical Need for Data Curation Standards in Materials Informatics

The importance of data curation has been recognized in multiple areas of research; however, the discussion of this important issue is only beginning to emerge in materials science. In this Perspective, we highlight the benefits of using the standardized data curation protocols in materials science and discuss current gaps in accurate and reproducible data reporting using case studies drawn from high-impact materials science papers and well-known databases such as the Crystallography Open Database (COD) and the Cambridge Structural Database (CSD). We argue that both experimental and computational materials scientists need to embrace a culture of rigorous data curation as part of modern research data management. We propose a sample data curation pipeline for materials chemistry and illustrate its use by creating two new materials chemistry databases. Here, we hope that this perspective will serve to catalyze further discussion and promote the continuous development of rigorous data curation practices within the materials science research community. We posit that adherence to best practices of data curation will promote and enhance the reliability, reproducibility, and integrity of materials research and enable the development of reliable AI and machine learning models that critically depend on the use of quality data.

Chemical structure↗

DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗

DOE OSTI2022

Theoretical Evaluation of Metal–Ligand Bonding in Neptunium Compounds in Relation to 237 Np Mössbauer Spectroscopy

The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Learning Coarse-Grained Potentials for Binary Fluids

Even though atomistic and CG models have been used to simulate liquid nanodroplets in vapor, very few rigorous studies of the liquid-liquid interface structure are available, and most of them are limited to planar interfaces. In this work, we evaluate several existing force fields (FF)s, including two atomistic and three CG FFs, with respect to modeling the interface structure and thermodynamic properties of the water-hexane interface. Both atomistic FFs are able to quantitatively reproduce the interfacial tension and the coexisting densities of the experimentally-observed planar interface. We use the atomistic FFs to model water droplets in hexane and use these simulations to test the CG FFs. We found that the tested CG FFs cannot reproduce the interfacial tensions of planar and/or curved interfaces. Finally, we propose a new approach for learning coarse-grained (CG) potentials within the CG SDK (Shinoda-DeVane- Klein) FF framework from atomistic simulation data. We demonstrate that the new potential significantly improves the prediction of both the interfacial tension and structure of water-hexane planar and curved interfaces.

Gao, Peiyuan↗

DOE OSTI2021

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗

DOE OSTI2025

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

DPD Parameters Estimation for Simultaneously Simulating Water–Oil Interfaces and Aqueous Nonionic Surfactants

The outcome of a coarse-grained simulation within the dissipative particle dynamics framework strongly depends on the choice of the repulsive parameter between different species. Different methodologies have been used in the literature to determine these parameters toward reproducing selected experimental system properties. In this work, a systematic investigation on possible procedures for estimating the simulation parameters is conducted. We compare methods based on the Hildebrand and the Hansen solubility parameter theories, mapped into the Flory–Huggins model. We find that using the Hansen solubility parameters it is possible to achieve a high degree of coarse graining, with parameters that yield realistic values for the interfacial tension. The procedure was first applied to the water/benzene system and then validated for water/n-octane, water/1,1-dichloroethane, water/methyl cyclohexane, and water/isobutyl acetate. In all these cases, the experimental interfacial tension could be reproduced by adjusting a single correction factor. In the case of the water–benzene system, the dissipative particle dynamics parameters derived using our approach were able to simultaneously describe both the interfacial tension and micellar properties of aqueous nonionic surfactants representative of the octyl poly(ethylene oxide) C8H17O(C2H4O)mH family. We show how the parameters can be used, within the dissipative particle dynamics framework, to simulate the water/oil interface in the presence of surfactants at varying concentrations. The results show, as expected, that as the surfactant concentration increases, the interfacial tension decreases, and micelles form in bulk water.

Chemistry↗