1
|
Karalis K, Zahn D, Prasianakis NI, Niceno B, Churakov SV. Deciphering the molecular mechanism of water boiling at heterogeneous interfaces. Sci Rep 2021; 11:19858. [PMID: 34615926 PMCID: PMC8494797 DOI: 10.1038/s41598-021-99229-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 09/17/2021] [Indexed: 11/23/2022] Open
Abstract
Water boiling control evolution of natural geothermal systems is widely exploited in industrial processes due to the unique non-linear thermophysical behavior. Even though the properties of water both in the liquid and gas state have been extensively studied experimentally and by numerical simulations, there is still a fundamental knowledge gap in understanding the mechanism of the heterogeneous nucleate boiling controlling evaporation and condensation. In this study, the molecular mechanism of bubble nucleation at the hydrophilic and hydrophobic solid-water interface was determined by performing unbiased molecular dynamics simulations using the transition path sampling scheme. Analyzing the liquid to vapor transition path, the initiation of small void cavities (vapor bubbles nuclei) and their subsequent merging mechanism, leading to successively growing vacuum domains (vapor phase), has been elucidated. The molecular mechanism and the boiling nucleation sites' location are strongly dependent on the solid surface hydrophobicity and hydrophilicity. Then simulations reveal the impact of the surface functionality on the adsorbed thin water molecules film structuring and the location of high probability nucleation sites. Our findings provide molecular-scale insights into the computational aided design of new novel materials for more efficient heat removal and rationalizing the damage mechanisms.
Collapse
Affiliation(s)
| | - Dirk Zahn
- Lehrstuhl für Theoretische Chemie/Computer Chemie Centrum, Friedrich-Alexander Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Nikolaos I Prasianakis
- Laboratory for Waste Management (LES), Paul Scherrer Institute, 5232, Villigen, Switzerland
| | - Bojan Niceno
- Laboratory of Scientific Computing and Modelling (LSM), Paul Scherrer Institute, 5232, Villigen, Switzerland
| | - Sergey V Churakov
- Institute of Geological Sciences, University of Bern, 3012, Bern, Switzerland.
- Laboratory for Waste Management (LES), Paul Scherrer Institute, 5232, Villigen, Switzerland.
| |
Collapse
|
2
|
Zahn D. Exploring the Mechanisms of Reactions in Solution from Transition Path Sampling Molecular Dynamics Simulations. J Chem Theory Comput 2015; 2:107-14. [PMID: 26626385 DOI: 10.1021/ct0501755] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recent advances in molecular dynamics simulations of rare reaction events and aggregation processes are reviewed. Therein the central focus is dedicated to employing the transition path sampling method to study reactions in solution. We describe systematic approaches for generating initial transition pathways and efficient strategies for computationally feasible exploration of further transition routes. The unprejudiced study of reaction mechanisms is illustrated for reactions in aqueous solution and other complex systems. Transition path sampling allows very detailed investigation of solvent effects. Apart from stabilization of reactant, transition, or product state ensembles, this also includes the role of the solvent as a heat bath and as a putative reaction partner. The latter issue is of particular importance for reactions in aqueous solutions, which involve proton-transfer steps that may be assisted by water molecules via the Grotthuss mechanism.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Strasse 40, 01187 Dresden, Germany
| |
Collapse
|
3
|
Ru Q, Fadda HM, Li C, Paul D, Khaw PT, Brocchini S, Zloh M. Molecular dynamic simulations of ocular tablet dissolution. J Chem Inf Model 2013; 53:3000-8. [PMID: 24073784 DOI: 10.1021/ci4002348] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Small tablets for implantation into the subconjunctival space in the eye are being developed to inhibit scarring after glaucoma filtration surgery (GFS). There is a need to evaluate drug dissolution at the molecular level to determine how the chemical structure of the active may correlate with dissolution in the nonsink conditions of the conjunctival space. We conducted molecular dynamics simulations to study the dissolution process of tablets derived from two drugs that can inhibit fibrosis after GFS, 5-fluorouracil (5-FU) and the matrix metalloprotease inhibitor (MMPi), ilomastat. The dissolution was simulated in the presence of simple point charge (SPC) water molecules, and the liquid turnover of the aqueous humor in the subconjunctival space was simulated by removal of the dissolved drug molecules at regular intervals and replacement by new water molecules. At the end of the simulation, the total molecular solvent accessible surface area of 5-FU tablets increased by 60 times more than that of ilomastat as a result of tablet swelling and release of molecules into solution. The tablet dissolution pattern shown in our molecular dynamic simulations tends to correlate with experimental release profiles. This work indicates that a series of molecular dynamic simulations can be used to predict the influence of the molecular properties of a drug on its dissolution profile and could be useful during preformulation where sufficient amounts of the drug are not always available to perform dissolution studies.
Collapse
Affiliation(s)
- Qian Ru
- UCL School of Pharmacy , 29/39 Brunswick Square, London WC1N 1AX, United Kingdom
| | | | | | | | | | | | | |
Collapse
|
4
|
Zahn D. Nucleation mechanism and kinetics of the perovskite to post-perovskite transition of MgSiO3 under extreme conditions. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.04.058] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
5
|
Zahn D, Anwar J. Collective displacements in a molecular crystal polymorphic transformation. RSC Adv 2013. [DOI: 10.1039/c3ra40653a] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
6
|
Zahn D. Tackling time-reversibility in transition path sampling molecular dynamics simulations. MOLECULAR SIMULATION 2012. [DOI: 10.1080/08927022.2011.614241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
7
|
Zahn D. Modeling martensic transformations in crystalline solids: validity and redesign of geometric approaches. ACTA ACUST UNITED AC 2011. [DOI: 10.1524/zkri.2011.1368] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
8
|
Zahn D, Tlatlik H. Atomistic In Situ Investigation of the Morphogenesis of Grains during Pressure-Induced Phase Transitions: Molecular Dynamics Simulations of the B1-B2 Transformation of RbCl. Chemistry 2010; 16:13385-9. [DOI: 10.1002/chem.201001481] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
9
|
Rowley CN, Woo TK. New shooting algorithms for transition path sampling: Centering moves and varied-perturbation sizes for improved sampling. J Chem Phys 2009; 131:234102. [DOI: 10.1063/1.3274203] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Frenzel J, Seifert G, Zahn D. Surface Effects in the Pressure-Induced Structural Transformation of a ZnO Nanorod. Z Anorg Allg Chem 2009. [DOI: 10.1002/zaac.200900260] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
11
|
Zahn D. Minimum energy pathways of brittle and ductile deformation/fracture processes. J Chem Phys 2008; 128:184707. [DOI: 10.1063/1.2918280] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Zahn D. Atomistic Mechanisms of Phase Separation and Formation of Solid Solutions: Model Studies of NaCl, NaCl−NaF, and Na(Cl1-xBrx) Crystallization from the Melt. J Phys Chem B 2007; 111:5249-53. [PMID: 17439270 DOI: 10.1021/jp068766u] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The crystallization of sodium chloride from its melt and mixtures with other sodium halides is investigated by means of transition path sampling molecular dynamics simulations. From this we explore the nucleation mechanisms of both the solidification and the melting process at the atomistic level of detail. By incorporation of impurities the nucleation picture of the eutectic mixtures changes considerably. Doping the NaCl crystal with fluoride ions, we observed the substitutional defects to act as favored nucleation centers for the melting transition. This phenomenon plays a critical role during the solidification process of NaCl-NaF melts of low NaF concentration and is demonstrated to account for the segregation of fluoride ions. While NaCl-NaF corresponds to a eutectic system, we also investigated NaCl-NaBr mixtures. The bromide ions were observed to behave very similarly to chloride ions. As a consequence, no phase separation occurs and Na(Cl1-xBrx) solid solutions are formed. At the example of these two prototypes we demonstrate the study of the atomistic mechanisms related to phase separation processes and solid solution formation during the nucleation and growth of crystals from multinary melts.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Strasse 40, 01187 Dresden, Germany.
| |
Collapse
|
13
|
Rowley CN, Woo TK. Generation of initial trajectories for transition path sampling of chemical reactions with ab initio molecular dynamics. J Chem Phys 2007; 126:024110. [PMID: 17228946 DOI: 10.1063/1.2424712] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Transition path sampling is an innovative method for focusing a molecular dynamics simulation on a reactive event. Although transition path sampling methods can generate an ensemble of reactive trajectories, an initial reactive trajectory must be generated by some other means. In this paper, the authors have evaluated three methods for generating initial reactive trajectories for transition path sampling with ab initio molecular dynamics. The authors have tested each of these methods on a set of chemical reactions involving the breaking and making of covalent bonds: the 1,2-hydrogen elimination in the borane-ammonia adduct, a tautomerization, and the Claisen rearrangement. The first method is to initiate trajectories from the potential energy transition state, which was effective for all reactions in the test set. Assigning atomic velocities found using normal mode analysis greatly improved the success of this method. The second method uses a high temperature molecular dynamics simulation and then iteratively reduces the total energy of the simulation until a low temperature reactive trajectory is found. This was effective in generating a low temperature trajectory from an initial trajectory run at 3000 K of the tautomerization reaction, although it failed for the other two. The third uses an orbital based bias potential to find a reactive trajectory and uses this trajectory to initiate an unbiased trajectory. The authors found that a highest occupied molecular orbital-lowest unoccupied molecular orbital bias could be used to find a reactive trajectory for the Claisen rearrangement, although it failed for the other two reactions. These techniques will help make it practical to use transition path sampling to study chemical reaction mechanisms that involve bond breaking and forming.
Collapse
Affiliation(s)
- Christopher N Rowley
- Centre for Catalysis Research and Innovation, Department of Chemistry, University of Ottawa, 10 Marie-Curie, Ottawa, Ontario, K1N 6N5, Canada
| | | |
Collapse
|
14
|
Zahn D. Competing Evaporation and Condensation Processes during the Boiling of Methane. J Phys Chem B 2006; 110:19601-4. [PMID: 17004826 DOI: 10.1021/jp062221i] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The atomistic mechanism of the boiling of methane is explored from molecular dynamics simulations. The liquid --> vapor transition is initiated by local density fluctuations resulting in a nanometer-sized domain that exhibits both liquid and vapor characteristics. Though the rates of evaporation and condensation events increase dramatically in this area, the overall balance exhibits only a marginal net rate of evaporation. Growth of the precritical domain leads to the nucleation of a vapor phase in which isolated methane molecules are confined by a liquid-vapor interface. After crossing the transition state, the system experiences progressive destabilization of the liquid phase and the evaporation processes clearly outnumber the condensation events.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Strasse 40, 01187 Dresden, Germany.
| |
Collapse
|
15
|
Zahn D, Leoni S. Mechanisms and Nucleation Characteristics of the Pressure-Induced B1−B2 Transition in Potassium Halides: A Question of Ion Hardness and Softness. J Phys Chem B 2006; 110:10873-7. [PMID: 16771339 DOI: 10.1021/jp061538f] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The transformation of potassium bromide from the B1 to the high-pressure B2 structure type is investigated by means of molecular dynamics simulations and compared to previous studies of KF and KCl. The underlying simulation scheme is based on the transition path sampling approach, which allows an unbiased investigation of the phase transition and offers a unique perspective for studying the involved mechanisms at the atomistic level of detail. Our analysis reveals identical mechanisms for the overall transition in KF, KCl, and KBr, but rather dissimilar characteristics of the nucleation and growth of phases. The transformation of KCl may be initiated by both K+ and Cl- ion displacement, exhibiting no preference for either species. However, for KF and KBr, we identified a clear favoring of column-wise F- and K+ displacement, respectively. Such tendencies have important implications on the morphogenesis of the phase nuclei and account for the observation of short-ranged coexisting nucleation centers, resulting in the formation of nanosized twin domains separated by mirror planes on completion of the transition. On the basis of a systematic study of potassium halides, we present a conclusive explanation for the observed nucleation characteristics, which is expected to be of general relevance to pressure-induced phase transitions in ionic compounds.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Str. 40, 01187 Dresden, Germany.
| | | |
Collapse
|
16
|
Zahn D. Unprejudiced identification of reaction mechanisms from biased transition path sampling. J Chem Phys 2005; 123:044104. [PMID: 16095343 DOI: 10.1063/1.1950672] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for unprejudiced investigation of reaction mechanisms from molecular-dynamics simulations is presented. It combines the transition path sampling approach with a biasing strategy, which (a) allows optimization of transition paths crossing an energy minimum of the transition state surface. The bias is then used to (b) find reaction pathways, which follow different mechanistic routes. In the first step the manifold of similar trajectories that correspond to the same mechanism is reduced to a single characteristic dynamical path. Our method then allows a systematic search for further reaction mechanisms and the related energy barriers. It is illustrated at the example of a single particle in a two-dimensional potential and of the rather complex process of the pressure-induced insertion of a helium atom into a C60 buckyball molecule.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Strasse 40, 01187 Dresden, Germany.
| |
Collapse
|
17
|
Abstract
Insight into the boiling of water is obtained from molecular dynamics simulations. The process is initiated by the spontaneous formation of small vacuum cavities in liquid water. By themselves, these defects are very short lived. If, however, several cavities occur at close distances, they are likely to merge into larger vacuum holes. At the liquid-vapor interfaces, single or small groups of water molecules tend to leave the liquid surface. Once the system is propagated beyond the transition state, these evaporation events outnumber the competing reintegration into the hydrogen-bonded network.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institut für Chemische Physik fester Stoffe, Nöthnitzer Strasse 40, 01187 Dresden, Germany.
| |
Collapse
|