1
|
Jobbins SA, Boulfelfel SE, Leoni S. Metashooting: a novel tool for free energy reconstruction from polymorphic phase transition mechanisms. Faraday Discuss 2018; 211:235-251. [PMID: 30039129 DOI: 10.1039/c8fd00053k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We introduce a novel scheme for the mechanistic investigation of solid-solid phase transitions, which we dub 'metashooting'. Combining transition path sampling, molecular dynamics and metadynamics, this scheme allows for both a complete mechanistic analysis and detailed mapping of the free energy surface. This is illustrated by performing metashooting calculations on the pressure-induced B4/B3 → B1 phase transition in ZnO. The resulting free energy map helps to clarify the role of intermediate configurations along this activated process and the competition between different mechanistic regimes with superior accuracy. We argue that metashooting can be efficiently applied to a broader class of activated processes.
Collapse
|
2
|
Selli D, Boulfelfel SE, Schapotschnikow P, Donadio D, Leoni S. Hierarchical thermoelectrics: crystal grain boundaries as scalable phonon scatterers. NANOSCALE 2016; 8:3729-3738. [PMID: 26815914 DOI: 10.1039/c5nr05279c] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Thermoelectric materials are strategically valuable for sustainable development, as they allow for the generation of electrical energy from wasted heat. In recent years several strategies have demonstrated some efficiency in improving thermoelectric properties. Dopants affect carrier concentration, while thermal conductivity can be influenced by alloying and nanostructuring. Features at the nanoscale positively contribute to scattering phonons, however those with long mean free paths remain difficult to alter. Here we use the concept of hierarchical nano-grains to demonstrate thermal conductivity reduction in rocksalt lead chalcogenides. We demonstrate that grains can be obtained by taking advantage of the reconstructions along the phase transition path that connects the rocksalt structure to its high-pressure form. Since grain features naturally change as a function of size, they impact thermal conductivity over different length scales. To understand this effect we use a combination of advanced molecular dynamics techniques to engineer grains and to evaluate thermal conductivity in PbSe. By affecting grain morphologies only, i.e. at constant chemistry, two distinct effects emerge: the lattice thermal conductivity is significantly lowered with respect to the perfect crystal, and its temperature dependence is markedly suppressed. This is due to an increased scattering of low-frequency phonons by grain boundaries over different size scales. Along this line we propose a viable process to produce hierarchical thermoelectric materials by applying pressure via a mechanical load or a shockwave as a novel paradigm for material design.
Collapse
Affiliation(s)
- Daniele Selli
- Cardiff University, School of Chemistry, Park Place, CF10 3AT, Cardiff, UK
| | | | | | | | | |
Collapse
|
3
|
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
|
4
|
Jansen M. The Azide/Nitrate Route, a Versatile and Prolific Procedure for the Synthesis of Alkali Oxometalates. Z Anorg Allg Chem 2012. [DOI: 10.1002/zaac.201200275] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
5
|
Boulfelfel SE, Oganov AR, Leoni S. Understanding the nature of "superhard graphite". Sci Rep 2012; 2:471. [PMID: 22745897 PMCID: PMC3384968 DOI: 10.1038/srep00471] [Citation(s) in RCA: 78] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2012] [Accepted: 04/30/2012] [Indexed: 11/09/2022] Open
Abstract
Numerous experiments showed that on cold compression graphite transforms into a new superhard and transparent allotrope. Several structures with different topologies have been proposed for this phase. While experimental data are compatible with most of these models, the only way to solve this puzzle is to find which structure is kinetically easiest to form. Using state-of-the-art molecular-dynamics transition path sampling simulations, we investigate kinetic pathways of the pressure-induced transformation of graphite to various superhard candidate structures. Unlike hitherto applied methods for elucidating nature of superhard graphite, transition path sampling realistically models nucleation events necessary for physically meaningful transformation kinetics. We demonstrate that nucleation mechanism and kinetics lead to M-carbon as the final product. W-carbon, initially competitor to M-carbon, is ruled out by phase growth. Bct-C4 structure is not expected to be produced by cold compression due to less probable nucleation and higher barrier of formation.
Collapse
|
6
|
|
7
|
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]
|
8
|
Grin Y, Nesper R. Periodic Space Partitioners (PSP) and their relations to Crystal Chemistry. ACTA ACUST UNITED AC 2011. [DOI: 10.1524/zkri.2011.1429] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Abstract
We are looking back on 30 years development of periodic space partitioners (PSP) and their relations to their periodic relatives, i.e. minimal surfaces (PMS), zero potential surfaces (P0PS), nodal surfaces (PNS), and exponential scale surfaces. Hans-Georg von Schnering and Sten Andersson have pioneered this field especially in terms of applications to crystal chemistry. This review relates the early attempts to approximate periodic minimal surfaces which established a systematic classification of all PSP in terms space group symmetry and consecutive applications in a variety of different fields. A consistent nomenclature is outlined and different methods for deriving PSP are described. Characteristic structure factor sets which solely define PNS by can be used to discriminate structure types of a given symmetry or even to determine complicated crystal structures. The concept of PSP relates space group symmetry, topology, and chemical bonding in an intriguing way and tessellations on PSP which can be generated in a straight forward way allow to predict new framework types. Through transformation of such continuous topological forms a new entry has been found for understanding and interpreting reconstructive phase transitions. Finally we indicate the importance of PSP models for soft matter science.
Collapse
Affiliation(s)
| | - Reinhard Nesper
- ETH Zürich, Laboratory of Inorganic Chemistry, Zürich, Schweiz
| |
Collapse
|
9
|
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]
|
10
|
Leoni S, Eddine Boulfelfel S, Baburin I. A Walk on the Chemical Landscape: The Role of B33 along the B1-B2 Phase Transition in RbF and NaBr. Z Anorg Allg Chem 2011. [DOI: 10.1002/zaac.201100045] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
11
|
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]
|
12
|
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
|
13
|
Nanodomain fragmentation and local rearrangements in CdSe under pressure. Proc Natl Acad Sci U S A 2008; 105:19612-6. [PMID: 19052229 DOI: 10.1073/pnas.0805235105] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Structural transformations in extended solids result from local atomic rearrangements and phase growth mechanisms. A broad class of technologically relevant properties critically depends on local structural issues connected with domain sizes, domain boundary geometries, and defects. However, a precise understanding of structural transformation mechanisms and domain formation is still an open question. Here, we demonstrate the feasibility of very detailed mechanistic investigations in real materials as a prerequisite for intelligent property control. We address the problem of domain fragmentation in bulk CdSe under pressure, jointly by molecular dynamics simulations, high-pressure experiments, and HR-TEM imaging. We show that domain fragmentation is taking place in the high-pressure regime, where nucleation events generate both zinc blende (B3) and wurtzite (B4) structural motifs and, in turn, cause the final lamellar appearance observable by high-resolution TEM. A changed nucleation pattern and a modified B3/B4 ratio represents the system's response to modified external stress conditions.
Collapse
|
14
|
Boulfelfel SE, Zahn D, Grin Y, Leoni S. Walking the path from B4- to B1-type structures in GaN. PHYSICAL REVIEW LETTERS 2007; 99:125505. [PMID: 17930518 DOI: 10.1103/physrevlett.99.125505] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2007] [Indexed: 05/25/2023]
Abstract
Molecular dynamics simulations are performed on the wurtzite-type structure (B4) to the rocksalt-type structure (B1) pressure-induced phase transition in GaN. From this, a nucleation and growth mechanism through a tetragonal metastable configuration is found. An intermediate of h-MgO type structure suggested from static calculations is ruled out. However, the pathway through the tetragonal intermediate may be altered by defect incorporation. While the overall transformation mechanism is preserved for both vacancies and Ga substitution by indium, already a 5% aluminum substitution establishes a transition route which avoids the tetragonal structure. Changes in the transformation mechanism and the resulting stabilization of the previously metastable high-pressure modification is elaborated by tracing the interplay of phase nucleation and growth and defect incorporation.
Collapse
Affiliation(s)
- S E Boulfelfel
- Max-Planck-Institut für Chemische Physik fester Stoffe, 01187 Dresden, Germany
| | | | | | | |
Collapse
|
15
|
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
|
16
|
Boulfelfel SE, Zahn D, Hochrein O, Grin Y, Leoni S. Sublattice Melting by Pressure: Superionic Conduction in Low-Dimensional Phase Interfaces. Z Anorg Allg Chem 2006. [DOI: 10.1002/zaac.200670090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
17
|
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
|
18
|
Molecular dynamics investigation of the pressure induced B1 to B2 phase transitions of RbBr. J SOLID STATE CHEM 2004. [DOI: 10.1016/j.jssc.2004.06.010] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
19
|
Zahn D, Leoni S. Nucleation and growth in pressure-induced phase transitions from molecular dynamics simulations: mechanism of the reconstructive transformation of NaCl to the CsCl-type structure. PHYSICAL REVIEW LETTERS 2004; 92:250201. [PMID: 15244985 DOI: 10.1103/physrevlett.92.250201] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2003] [Indexed: 05/24/2023]
Abstract
We perform path sampling molecular dynamics on the pressure-induced reconstructive phase transition from NaCl to CsCl type structure. Unlike the molecular dynamics simulations prior to this work our approach does not drive the process by applying elevated pressure. As a consequence, we are able to observe nucleation events that initiate the successive transformation of the crystal. The competing phases are separated by an interface exhibiting a well-defined structure that propagates through the crystal during phase transition.
Collapse
Affiliation(s)
- Dirk Zahn
- Max-Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany.
| | | |
Collapse
|