Search NASA⌕ Search

SEARCH · Search NASA

Results for “empirical solvent polarity”

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.

Molecular Lignin Solubility and Structure in Organic Solvents

Lignin, a polymer found in the secondary plant cell wall of terrestrial plants, is the single largest source of renewable aromatics and has attracted considerable attention as a feedstock for potential industrial use. However, the secondary plant cell wall is a crowded environment, and lignin in its native form interacts with other biomass components within a larger network. Application of some organic solvents is known to liberate lignin from this network and creates lignin-rich streams suitable for conversion into target products. Through molecular-scale lignin simulation, we analyze how diverse lignin polymers change their structure in response to varying organic solvent environments. We quantify the relationship between solvent polarity and lignin polymer extension, observing maximum polymer expansion and solvation for solvents with polarity near those of dimethyl sulfoxide. From our observations at the nanoscale, increasing syringyl content within lignin polymers reduces the expansion of the polymer in organic solvent environments and decreases the free energy difference compared to aqueous solvent environments, thereby reducing solubility for high syringyl lignin polymers. The conformational transition rates between lignin polymer shapes increased through a combination of the solvent diffusion constant and polymer extension. The molecular simulations indicate that there is likely no single optimal organic solvent for lignin. Different solvent mixtures have optimal or near-optimal properties in solubilizing lignin polymers, thereby disrupting interactions with other biopolymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗