Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThermal ignition revisited with two-dimensional molecular dynamics: role of fluctuations in activated collisions
49
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The problem of thermal ignition in a homogeneous gas is revisited from a molecular dynamics perspective. A two-dimensional model is adopted, which assumes reactive disks of type A and B in a fixed area that react to form type C products if an activation threshold for impact is surpassed. Such a reaction liberates kinetic energy to the product particles, representative of the heat release. The results for the ignition delay are compared with those obtained from the continuum description assuming local thermodynamic equilibrium, in order to assess the role played by molecular fluctuations. Results show two regimes of non-equilibrium ignition whereby ignition occurs at different times as compared to that from the continuum description. The first regime is at low activation energies, where the ignition time is found to be higher than that expected from theory for all values of heat release, in agreement with predictions from Prigogine and Xhrouet who attribute this departure to non-equilibrium effects. Results suggest the ignition is spatially homogeneous in this regime. The second regime occurs at high activation energies and sufficiently large heat release values. In this regime, ignition times are found to be dependent on domain size, with larger domains yielding shorter ignition delays than expected. Results for larger systems agree with the expectations by Prigogine and Mahieu, who predict a non-equilibrium reaction rate larger than expected for a homogeneous system in equilibrium. Results yield a large variance for ignition times under these conditions, suggesting a departure from homogeneous combustion. The results obtained are in qualitative agreement with experimental observations of auto-ignition at relatively low temperatures, where hot-spot ignition and associated ignition delays lower than predicted are generally observed.
rate research
Read More
The present study addresses the role of molecular non-equilibrium effects in thermal ignition problems. We consider a single binary reaction of the form A+B -> C+C. Molecular dynamics calculations were performed for activation energies ranging between RT and 7.5RT and heat release of 2.5RT and 10RT. The evolution of up to 10,000 particles was calculated as the system undergoes a thermal ignition at constant volume. Ensemble averages of 100 calculations for each parameter set permitted to determine the ignition delay, along with a measure of the stochasticity of the process. A well behaved convergence to large system sizes is also demonstrated. The ignition delay calculations were compared with those obtained at the continuum level using rates derived from kinetic theory: the standard rate assuming that the distribution of the speed of the particles is the Maxwell-Boltzmann distribution, and the perturbed rates by Prigogine and Xhrouet [1] for an isothermal system, and Prigogine and Mahieu [2] for an energy releasing reaction, obtained by the Chapman-Enskog perturbation procedure. The molecular results were found in very good agreement with the latter at low temperatures, confirming that non-equilibrium effects promote the formation of energetic particles, that serve as seeds for subsequent reaction events: i.e., hot spots. This effect was found to lower the ignition delay by up to 30%. At high temperatures, the ignition delay obtained from the standard equilibrium rate was found to be up to 60% longer than the molecular calculations. This effect is due to the rapidity of the reactive collisions that do not allow the system to equilibrate. For this regime, none of the perturbation solutions obtained by the Chapman-Enskog procedure were valid. This study thus shows the importance of non-equilibrium effects in thermal ignition problems, for most temperatures of practical interest.
Molecular Dynamics studies of chemical processes in solution are of great value in a wide spectrum of applications, which range from nano-technology to pharmaceutical chemistry. However, these calculations are affected by severe finite-size effects, such as the solution being depleted as the chemical process proceeds, which influence the outcome of the simulations. To overcome these limitations, one must allow the system to exchange molecules with a macroscopic reservoir, thus sampling a Grand-Canonical ensemble. Despite the fact that different remedies have been proposed, this still represents a key challenge in molecular simulations. In the present work we propose the Constant Chemical Potential Molecular Dynamics (C$mu$MD) method, which introduces an external force that controls the environment of the chemical process of interest. This external force, drawing molecules from a finite reservoir, maintains the chemical potential constant in the region where the process takes place. We have applied the C$mu$MD method to the paradigmatic case of urea crystallization in aqueous solution. As a result, we have been able to study crystal growth dynamics under constant supersaturation conditions, and to extract growth rates and free-energy barriers.
A comprehensive microscopic understanding of ambient liquid water is a major challenge for $ab$ $initio$ simulations as it simultaneously requires an accurate quantum mechanical description of the underlying potential energy surface (PES) as well as extensive sampling of configuration space. Due to the presence of light atoms (e.g., H or D), nuclear quantum fluctuations lead to observable changes in the structural properties of liquid water (e.g., isotope effects), and therefore provide yet another challenge for $ab$ $initio$ approaches. In this work, we demonstrate that the combination of dispersion-inclusive hybrid density functional theory (DFT), the Feynman discretized path-integral (PI) approach, and machine learning (ML) constitutes a versatile $ab$ $initio$ based framework that enables extensive sampling of both thermal and nuclear quantum fluctuations on a quite accurate underlying PES. In particular, we employ the recently developed deep potential molecular dynamics (DPMD) model---a neural-network representation of the $ab$ $initio$ PES---in conjunction with a PI approach based on the generalized Langevin equation (PIGLET) to investigate how isotope effects influence the structural properties of ambient liquid H$_2$O and D$_2$O. Through a detailed analysis of the interference differential cross sections as well as several radial and angular distribution functions, we demonstrate that this approach can furnish a semi-quantitative prediction of these subtle isotope effects.
We used molecular dynamics simulations to predict the steady state crystal shape of naphthalene grown from ethanol solution. The simulations were performed at constant supersaturation by utilizing a recently proposed algorithm [Perego et al., J. Chem. Phys., 142, 2015, 144113]. To bring the crystal growth within the timescale of a molecular dynamics simulation we applied Well-Tempered Metadynamics with a spatially constrained collective variable, which focuses the sampling on the growing layer. We estimated that the resulting steady state crystal shape corresponds to a rhombic prism, which is in line with experiments. Further, we observed that at the investigated supersaturations, the ${00bar{1}}$ face grows in a two step two dimensional nucleation mechanism while the considerably faster growing faces ${1bar{1}0}$ and ${20bar{1}}$ grow new layers with a one step two dimensional nucleation mechanism.
We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
