SEARCH · Engineering Papers

Results for “Parameterization”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Perspectives on Systematic Cloud Microphysics Scheme Development With Machine Learning

Cloud microphysics—the collection of processes that govern the small‐scale formation, evolution, and interactions of liquid droplets and ice crystals in clouds and precipitation—remains a major source of uncertainty in weather and climate models. Although too small in scale to be explicitly resolved in any large‐eddy simulation, weather, or climate model, the representation of cloud microphysical processes has significant impact at the climate scale. Current microphysical schemes are limited by both parametric uncertainty, linked to uncertainty in physical parameter values, and structural uncertainty, arising from incomplete physical understanding of the processes at play or approximations made for computational efficiency. Recent advances in the application of machine learning (ML) to the physical sciences show significant potential for minimizing these limitations by leveraging high‐fidelity simulations and observations. Here we outline the challenges that must be addressed to apply ML toward cloud microphysics scheme development. This perspectives paper synthesizes recent progress in using data‐driven methods, including ML, to improve cloud microphysics parameterizations and highlights opportunities to address key uncertainties. We discuss the roles of aleatoric (irreducible, or statistical) and epistemic (reducible, or systematic) errors in contributing to microphysics parameterization uncertainty. ML can leverage observations to improve microphysical schemes via bottom‐up and top‐down constraints. Methods such as differentiable programming and ML‐enhanced sampling strategies and the creation of large scale benchmark data sets promise to bridge the gap between observations and models and to improve the consistency of cloud microphysical representation across temporal and spatial scales.

Lamb, Kara D. [Columbia Univ., New York, NY (Unite↗

DOE OSTI2020

Genetic circuit characterization by inferring RNA polymerase movement and ribosome usage

To perform their computational function, genetic circuits change states through a symphony of genetic parts that turn regulator expression on and off. Debugging is frustrated by an inability to characterize parts in the context of the circuit and identify the origins of failures. Here, we take snapshots of a large genetic circuit in different states: RNA-seq is used to visualize circuit function as a changing pattern of RNA polymerase (RNAP) flux along the DNA. Together with ribosome profiling, all 54 genetic parts (promoters, ribozymes, RBSs, terminators) are parameterized and used to inform a mathematical model that can predict circuit performance, dynamics, and robustness. The circuit behaves as designed; however, it is riddled with genetic errors, including cryptic sense/antisense promoters and translation, attenuation, incorrect start codons, and a failed gate. While not impacting the expected Boolean logic, they reduce the prediction accuracy and could lead to failures when the parts are used in other designs. Finally, the cellular power (RNAP and ribosome usage) required to maintain a circuit state is calculated. This work demonstrates the use of a small number of measurements to fully parameterize a regulatory circuit and quantify its impact on host.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Multi-year incubation experiments boost confidence in model projections of long-term soil carbon dynamics

Abstract Global soil organic carbon (SOC) stocks may decline with a warmer climate. However, model projections of changes in SOC due to climate warming depend on microbially-driven processes that are usually parameterized based on laboratory incubations. To assess how lab-scale incubation datasets inform model projections over decades, we optimized five microbially-relevant parameters in the Microbial-ENzyme Decomposition (MEND) model using 16 short-term glucose (6-day), 16 short-term cellulose (30-day) and 16 long-term cellulose (729-day) incubation datasets with soils from forests and grasslands across contrasting soil types. Our analysis identified consistently higher parameter estimates given the short-term versus long-term datasets. Implementing the short-term and long-term parameters, respectively, resulted in SOC loss (–8.2 ± 5.1% or –3.9 ± 2.8%), and minor SOC gain (1.8 ± 1.0%) in response to 5 °C warming, while only the latter is consistent with a meta-analysis of 149 field warming observations (1.6 ± 4.0%). Comparing multiple subsets of cellulose incubations (i.e., 6, 30, 90, 180, 360, 480 and 729-day) revealed comparable projections to the observed long-term SOC changes under warming only on 480- and 729-day. Integrating multi-year datasets of soil incubations (e.g., > 1.5 years) with microbial models can thus achieve more reasonable parameterization of key microbial processes and subsequently boost the accuracy and confidence of long-term SOC projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Global quantification of the dispersion effect with POLDER satellite data

Increased aerosols can modify the shape of the cloud Particle Size Distribution (PSD), thereby influencing the radiative properties of clouds, known as the Dispersion Effect (DE). However, a global, observation-based quantification of its impact on Aerosol-Cloud Interactions (ACI) is lacking, leading to DE being typically ignored in satellite-based estimates of ACI forcing. Here we propose a physics-based method that combines polarimetric satellite data on cloud PSD to achieve global observational quantification of DE’s impact on ACI in liquid-phase stratiform clouds. Globally, DE offsets ACI changes induced by droplet number concentration variation and liquid water path adjustment by 7% and −1.4%, respectively. Furthermore, a parameterization based on the global dataset of PSD shape parameters is developed to improve DE estimation in large-scale models. Both the quantification and parameterization enhance our understanding of DE and facilitate the inclusion of this non-negligible impact of DE on ACI in estimating aerosol climate forcing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Predicting the influence of particle size on the glass transition temperature and viscosity of secondary organic material

Abstract Atmospheric aerosols can assume liquid, amorphous semi-solid or glassy, and crystalline phase states. Particle phase state plays a critical role in understanding and predicting aerosol impacts on human health, visibility, cloud formation, and climate. Melting point depression increases with decreasing particle diameter and is predicted by the Gibbs–Thompson relationship. This work reviews existing data on the melting point depression to constrain a simple parameterization of the process. The parameter $$\xi $$ ξ describes the degree to which particle size lowers the melting point and is found to vary between 300 and 1800 K nm for a wide range of particle compositions. The parameterization is used together with existing frameworks for modeling the temperature and RH dependence of viscosity to predict the influence of particle size on the glass transition temperature and viscosity of secondary organic aerosol formed from the oxidation of $$\alpha $$ α -pinene. Literature data are broadly consistent with the predictions. The model predicts a sharp decrease in viscosity for particles less than 100 nm in diameter. It is computationally efficient and suitable for inclusion in models to evaluate the potential influence of the phase change on atmospheric processes. New experimental data of the size-dependence of particle viscosity for atmospheric aerosol mimics are needed to thoroughly validate the predictions.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Data-driven causal model discovery and personalized prediction in Alzheimer's disease

Abstract With the explosive growth of biomarker data in Alzheimer’s disease (AD) clinical trials, numerous mathematical models have been developed to characterize disease-relevant biomarker trajectories over time. While some of these models are purely empiric, others are causal, built upon various hypotheses of AD pathophysiology, a complex and incompletely understood area of research. One of the most challenging problems in computational causal modeling is using a purely data-driven approach to derive the model’s parameters and the mathematical model itself, without any prior hypothesis bias. In this paper, we develop an innovative data-driven modeling approach to build and parameterize a causal model to characterize the trajectories of AD biomarkers. This approach integrates causal model learning, population parameterization, parameter sensitivity analysis, and personalized prediction. By applying this integrated approach to a large multicenter database of AD biomarkers, the Alzheimer’s Disease Neuroimaging Initiative, several causal models for different AD stages are revealed. In addition, personalized models for each subject are calibrated and provide accurate predictions of future cognitive status.

Zheng, Haoyang (ORCID:0000000168358242)↗

DOE OSTI2024

CryoTRANS: predicting high-resolution maps of rare conformations from self-supervised trajectories in cryo-EM

Cryogenic electron microscopy (cryo-EM) has revolutionized structural biology, enabling efficient determination of structures at near-atomic resolutions. However, a common challenge arises from the severe imbalance among various conformations of vitrified particles, leading to low-resolution reconstructions in rare conformations due to a lack of particle images in these quasi-stable states. We introduce CryoTRANS, a method that predicts high-resolution maps of rare conformations by constructing a self-supervised pseudo-trajectory between density maps of varying resolutions. This trajectory is represented by an ordinary differential equation parameterized by a deep neural network, ensuring retention of detailed structures from high-resolution density maps. By leveraging a single high-resolution density map, CryoTRANS significantly improves the reconstruction of rare conformations and has been validated on four real-world datasets: alpha-2-macroglobulin, actin-binding protein complexes, SARS-CoV-2 spike glycoprotein, and the 70S ribosome. CryoTRANS can also predict high-resolution structures in cryogenic electron tomography maps using a high-resolution cryo-EM map.Cryogenic electron microscopy (cryo-EM) has revolutionized structural biology, enabling efficient determination of structures at near-atomic resolutions. However, a common challenge arises from the severe imbalance among various conformations of vitrified particles, leading to low-resolution reconstructions in rare conformations due to a lack of particle images in these quasi-stable states. We introduce CryoTRANS, a method that predicts high-resolution maps of rare conformations by constructing a self-supervised pseudo-trajectory between density maps of varying resolutions. This trajectory is represented by an ordinary differential equation parameterized by a deep neural network, ensuring retention of detailed structures from high-resolution density maps. By leveraging a single high-resolution density map, CryoTRANS significantly improves the reconstruction of rare conformations and has been validated on four real-world datasets: alpha-2-macroglobulin, actin-binding protein complexes, SARS-CoV-2 spike glycoprotein, and the 70S ribosome. CryoTRANS can also predict high-resolution structures in cryogenic electron tomography maps using a high-resolution cryo-EM map.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Atmospheric ice nucleation

Atmospheric ice nucleation is crucial for global precipitation and affects the structure, lifetime and reflectivity of clouds, thereby impacting climate. Ice nucleates in various ways from aerosol particles, termed ice-nucleating particles, over an extensive temperature and humidity range. Quantifying the kinetic and thermodynamic regimes of nucleation is necessary to relate fundamental physics to theoretically based predictions of ice formation for implementation in cloud and climate models. We review how the molecular picture of ice nucleation has advanced in recent years and consequential impacts on the interpretation and parameterization of ice nucleation. Advances include the role of interfacial free energy and pressure on ice nucleation rates, mobility regions of water that generate the critical ice nucleus, classical and non-classical pathways of nucleation, the type of ice polymorph that forms, the impact of solutes on freezing and the role of nanopores as surface features promoting ice nucleation. We also introduce currently debated and evaluated freezing parameterizations for application in model environments. Finally, we outline what we believe are the current needs for improving predictive understanding of ice nucleation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Variational approach to quantum state tomography based on maximal entropy formalism

Quantum state tomography is an integral part of quantum computation and offers the starting point for the validation of various quantum devices. One of the central tasks in the field of state tomography is to reconstruct, with high fidelity, the quantum states of a quantum system. From an experiment on a real quantum device, one can obtain the mean measurement values of different operators. With such data as input, in this report we employ the maximal entropy formalism to construct the least biased mixed quantum state that is consistent with the given set of expectation values. Even though, in principle, the reported formalism is quite general and should work for an arbitrary set of observables, in practice we shall demonstrate the efficacy of the algorithm on an informationally complete (IC) set of Hermitian operators. Such a set possesses the advantage of uniquely specifying a single quantum state from which the experimental measurements have been sampled and hence renders the rare opportunity not only to construct a least-biased quantum state but even replicate the exact state prepared experimentally within a preset tolerance. Here, the primary workhorse of the algorithm is reconstructing an energy function which we designate as the effective Hamiltonian of the system, and parameterizing it with Lagrange multipliers, according to the formalism of maximal entropy. These parameters are thereafter optimized variationally so that the reconstructed quantum state of the system converges to the true quantum state within an error threshold. To this end, we employ a parameterized quantum circuit and a hybrid quantum-classical variational algorithm to obtain such a target state, making our recipe easily implementable on a near-term quantum device.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗

DOE OSTI2025

Challenges of standard halo models in constraining galaxy properties from cosmic infrared background anisotropies

The halo model, combined with halo occupation distribution (HOD) prescriptions, is widely used to interpret cosmic infrared background (CIB) anisotropies and extract physical information about star-forming galaxies and their connection to large-scale structures. Recent CIB-specific implementations of the halo model have adopted more physical parameterizations. However, the extent to which these models can reliably recover meaningful physical parameters remains uncertain. We assessed whether the current parameterization of CIB halo models is sufficient to recover astrophysical quantities, such as star formation efficiency, η(M h , z), and halo mass at which the peak of star formation efficiency occurs, M max , when fit to mock data. We also assessed whether discrepancies arise from assumptions about galaxy emission (the HOD ingredients) or from more fundamental components in the halo model, such as bias and matter clustering. We fit the M21 CIB HOD model, implemented within the halo model framework, to mock CIB power spectra and star formation rate density (SFRD) data generated from the SIDES-Uchuu simulation, and compared the best-fit parameters to the known simulation inputs. We then repeated the analysis using a simplified version of the simulation (SSU), explicitly designed to match the HOD assumptions. A detailed comparison of model and simulation outputs was carried out to trace the origin of observed discrepancies. While the M21 HOD model provides a good fit to the mock data, it failed to recover the intrinsic parameters accurately, particularly the halo mass at which star formation efficiency peaks. This mismatch persists even when fitting data generated with the same model assumptions. We find strong agreement (within 5%) in the emission-related components (SFRD, emissivity), but observe a scale- and redshift-dependent offset exceeding 20% in the two-halo term of the CIB power spectrum. This likely arises from limitations in the treatment of halo bias and matter clustering within the linear approximation. Additionally, incorporating scatter in the SFR–halo mass relation and the spectral energy distribution (SED) templates significantly affects the shot noise (∼50%), but has only a modest impact (less than 10%) on the clustered component. These results suggest that recovering physical parameters from CIB clustering requires improvements to the cosmological ingredients of the halo model framework, such as adopting scale-dependent halo bias and nonlinear matter power spectra in addition to careful modeling of emission physics.

cosmic background radiation↗

DOE OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗

DOE OSTI2021

A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study

We investigate the viability of the phaseless finite-temperature auxiliary-field quantum Monte Carlo (ph-FT-AFQMC) method for ab initio systems using the uniform electron gas as a model. Through comparisons with exact results and FT coupled cluster theory, we find that ph-FT-AFQMC is sufficiently accurate at high to intermediate electronic densities. We show, both analytically and numerically, that the phaseless constraint at FT is fundamentally different from its zero-temperature counterpart (i.e., ph-ZT-AFQMC), and generally, one should not expect ph-FT-AFQMC to agree with ph-ZT-AFQMC in the low-temperature limit. With an efficient implementation, we are able to compare exchange-correlation energies to the existing results in the thermodynamic limit and find that the existing parameterizations are highly accurate. In particular, we found that ph-FT-AFQMC exchange-correlation energies are in better agreement with a known parameterization than is restricted path-integral MC in the regime of Θ ≤ 0.5 and rs ≤ 2, which highlights the strength of ph-FT-AFQMC.

Lee, Joonho (ORCID:0000000296671081)↗

DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗

DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗

DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

The “good,” the “bad,” and the “hidden” in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl 2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca 2+ –Cl - interactions were primarily “hidden” in the Ca-NDIS signal due to overlapping with Ca 2+ –O w and Ca 2+ –H w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models “on average,” the question “which model is the best” has not been answered but replaced by the question “which model is better for a given research.” An overall “good” model can still be inappropriate or, in some instances, “bad” and, unfortunately, produce erroneous results. Additionally, the accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗