Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDistributed Charge Models of Liquid Methane and Ethane for Dielectric Effects and Solvation
83
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Liquid hydrocarbons are often modeled with fixed, symmetric, atom-centered charge distributions and Lennard-Jones interaction potentials that reproduce many properties of the bulk liquid. While useful for a wide variety of applications, such models cannot capture dielectric effects important in solvation, self-assembly, and reactivity. The dielectric constants of hydrocarbons, such as methane and ethane, physically arise from electronic polarization fluctuations induced by the fluctuating liquid environment. In this work, we present non-polarizable, fixed-charge models of methane and ethane that break the charge symmetry of the molecule to create fixed molecular dipoles, the fluctuations of which reproduce the experimental dielectric constant. These models can be considered a mean-field-like approximation that can be used to include dielectric effects in large-scale molecular simulations of polar and charged molecules in liquid methane and ethane. We further demonstrate that solvation of model ionic solutes and a water molecule in these fixed-dipole models improve upon dipole-free models.
rate research
Read More
We have performed a combined experimental and theoretical study of ethane and methane at high pressures up to 120 GPa at 300 K using x-ray diffraction and Raman spectroscopy and the USPEX ab-initio evolutionary structural search algorithm, respectively. For ethane, we have determined the crystallization point, for room temperature, at 2.7 GPa and also the low pressure crystal structure (Phase A). This crystal structure is orientationally disordered (plastic phase) and deviates from the known crystal structures for ethane at low temperatures. Moreover, a pressure induced phase transition has been identified, for the first time, at 13.6 GPa to a monoclinic phase B, the structure of which is solved based on a good agreement of the experimental results and theoretical predictions. For methane, our XRD measurements are in agreement with the previously reported high-pressure structures and EOS. We have determined the equations of state of ethane and methane, which provides a solid basis for the discussion of their relative stability at high pressures.
Continuum solvation models enable electronic structure calculations of systems in liquid environments, but because of the large number of empirical parameters, they are limited to the class of systems in their fit set (typically organic molecules). Here, we derive a solvation model with no empirical parameters for the dielectric response by taking the linear response limit of a classical density functional for molecular liquids. This model directly incorporates the nonlocal dielectric response of the liquid using an angular momentum expansion, and with a single fit parameter for dispersion contributions it predicts solvation energies of neutral molecules with an RMS error of 1.3 kcal/mol in water and 0.8 kcal/mol in chloroform and carbon tetrachloride. We show that this model is more accurate for strongly polar and charged systems than previous solvation models because of the parameter-free electric response, and demonstrate its suitability for ab initio solvation, including self-consistent solvation in quantum Monte Carlo calculations.
The predictive simulation of molecular liquids requires models that are not only accurate, but computationally efficient enough to handle the large systems and long time scales required for reliable prediction of macroscopic properties. We present a new approach to the systematic approximation of the first-principles potential energy surface (PES) of molecular liquids using the GAP (Gaussian Approximation Potential) framework. The approach allows us to create potentials at several different levels of accuracy in reproducing the true PES, which allows us to test the level of quantum chemistry that is necessary to accurately predict its macroscopic properties. We test the approach by building potentials for liquid methane (CH$_4$), which is difficult to model from first principles because its behavior is dominated by weak dispersion interactions with a significant many-body component. We find that an accurate, consistent prediction of its bulk density across a wide range of temperature and pressure requires not only many-body dispersion, but also quantum nuclear effects to be modeled accurately.
Continuum solvation models enable efficient first principles calculations of chemical reactions in solution, but require extensive parametrization and fitting for each solvent and class of solute systems. Here, we examine the assumptions of continuum solvation models in detail and replace empirical terms with physical models in order to construct a minimally-empirical solvation model. Specifically, we derive solvent radii from the nonlocal dielectric response of the solvent from ab initio calculations, construct a closed-form and parameter-free weighted-density approximation for the free energy of the cavity formation, and employ a pair-potential approximation for the dispersion energy. We show that the resulting model with a single solvent-independent parameter: the electron density threshold ($n_c$), and a single solvent-dependent parameter: the dispersion scale factor ($s_6$), reproduces solvation energies of organic molecules in water, chloroform and carbon tetrachloride with RMS errors of 1.1, 0.6 and 0.5 kcal/mol respectively. We additionally show that fitting the solvent-dependent $s_6$ parameter to the solvation energy of a single non-polar molecule does not substantially increase these errors. Parametrization of this model for other solvents, therefore, requires minimal effort and is possible without extensive databases of experimental solvation free energies.
Quantum-chemical processes in liquid environments impact broad areas of science, from molecular biology to geology to electrochemistry. While density-functional theory (DFT) has enabled efficient quantum-mechanical calculations which profoundly impact understanding of atomic-scale phenomena, realistic description of the liquid remains a challenge. Here, we present an approach based on joint density-functional theory (JDFT) which addresses this challenge by leveraging the DFT approach not only for the quantum mechanics of the electrons in a solute, but also simultaneously for the statistical mechanics of the molecules in a surrounding equilibrium liquid solvent. Specifically, we develop a new universal description for the interaction of electrons with an arbitrary liquid, providing the missing link to finally transform JDFT into a practical tool for the realistic description of chemical processes in solution. This approach predicts accurate solvation free energies and surrounding atomic-scale liquid structure for molecules and surfaces in multiple solvents without refitting, all at a fraction of the computational cost of methods of comparable detail and accuracy. To demonstrate the potential impact of this method, we determine the structure of the solid/liquid interface, offering compelling agreement with more accurate (but much more computationally intensive) theories and with X-ray reflectivity measurements.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
