Simulations of NaCl Aggregation from Solution: Solvent Determines Topography of Free Energy Landscape

Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.

Behavior of Linear and Nonlinear Dimensionality Reduction for Collective Variable Identification of Small Molecule Solution-Phase Reactions

Identifying collective variables (CVs) for chemical reactions is essential to reduce the 3N-dimensional energy landscape into lower dimensional basins and barriers of interest. However, in condensed phase processes, the nonmeaningful motions of bulk solvent often overpower the ability of dimensionality reduction methods to identify correlated motions that underpin collective variables. Yet solvent can play important indirect or direct roles in reactivity, and much can be lost through treatments that remove or dampen solvent motion. This has been amply demonstrated within principal component analysis (PCA), although less is known about the behavior of nonlinear dimensionality reduction methods, e.g., uniform manifold approximation and projection (UMAP), that have become recently utilized. The latter presents an interesting alternative to linear methods though often at the expense of interpretability. This work presents distance-attenuated projection methods of atomic coordinates that facilitate the application of both PCA and UMAP to identify collective variables in the presence of explicit solvent and further the specific identity of solvent molecules that participate in chemical reactions. The performance of both methods is examined in detail for two reactions where the explicit solvent plays very different roles within the collective variables. When applied to raw molecular dynamics data in solution, both PCA and UMAP representations are dominated by bulk solvent motions. On the other hand, when applied to data preprocessed by our attenuated projection methods, both PCA and UMAP identify the appropriate collective variables (though varying sensitivity is observed due to the presence of explicit solvent that results from the projection method). Importantly, this approach allows identification of specific solvent molecules that are relevant to the CVs and their importance.

Free energy of micellization of dodecyl phosphocholine (DPC) from molecular simulation: Hybrid PEACH-BAR method

A new method to extract the free energy of aggregation versus aggregate size from molecular simulation data is proposed and applied to a united atom model of the zwitterionic surfactant dodecyl phosphocholine in water. This system's slow dissociation rate and low critical micelle concentration (CMC of approximately 1-2 mM) make extraction of cluster free energies directly from simulation results using the "partition-enabled analysis of cluster histogram" (PEACH) method impractical. The new approach applies PEACH to a model with weakened attractions between aggregants, which allows sampling of a continuous range of cluster sizes, then recovers the free energy of aggregation under the original fully-attractive force field using the BAR free energy difference method. PEACH-BAR results are compared with free energy differences calculated via umbrella sampling, and are used to make predictions of CMC, average cluster size, and SAXS scattering profiles that are in fair agreement with experiment.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

Low-Noise Nanopore Enables In-Situ and Label-Free Tracking of a Trigger-Induced DNA Molecular Machine at the Single-Molecular Level

Solid-state nanopores have shown special high potential in a label-free molecular assay, structure identification, and target-index at the single-molecular level, even though frustrating electrical baseline noise is still one of the major factors that limit the spatial resolution and signaling reliability of solid-state nanopores, especially in small target detection. Here we develop a significant and easy-operating noise-reduction approach via mixing organic solvents with high dielectric constants into a traditional aqueous electrolyte. The strategy is generally effective for pores made of different materials, such as the most commonly used conical glass (CGN) or SiN_x. While the mechanism should be multisourced, MD simulations suggest the noise reduction may partially arise from the even ionic distribution caused by the addition of higher dielectric species. Among all solvents experimentally tested, the two with the highest dielectric constants, formamide and methylformamide, exhibit the best noise reduction effect for target detection of CGN. The power spectral density at the low-frequency limit is reduced by nearly 3 orders with the addition of 20% formamide. Our work qualifies the reliability of solid-state nanopores into much subtler scales of detection, such as dsDNAs under 100 bp. As a practical example, bare CGN is innovatively employed to perform in-situ tracking of trigger-responsive DNA machine forming oligomers.

Homogeneous Nucleation Mechanism of NaCl in Aqueous Solutions

In this study, molecular dynamic simulations are employed to investigate the homogeneous nucleation mechanism of NaCl crystal in solutions. According to the simulations, the dissolved behaviors of NaCl in water are dependent on ion concentrations. With increasing NaCl concentrations, the dissolved Na+ and Cl- ions tend to be aggregated in solutions. In combination with our recent studies, the aggregate of dissolved solutes is mainly ascribed to the hydrophobic interactions. Different from the two-step mechanism, no barrier is needed to overcome the formation of the aggregate. In comparison with the classical nucleation theory (CNT), because of the formation of solute aggregate, this lowers the barrier height of nucleation and affects the nucleation mechanism of NaCl crystal in water.

Brine rejection and hydrate formation upon freezing of NaCl aqueous solutions

Studying the freezing of saltwater on a molecular level is of fundamental importance for improving freeze desalination techniques.