D.
Selli
a,
S. E.
Boulfelfel
b,
I. A.
Baburin
a,
G.
Seifert
a and
S.
Leoni
*a
aTechnische Universität Dresden, Institut für Physikalische Chemie, D-01062, Dresden, Germany. E-mail: stefano.leoni@chemie.tu-dresden.de
bStony Brook University, Department of Geosciences, New York, 11794-2100, USA
First published on 26th July 2012
Using molecular dynamics simulations and a Density Functional based Tight Binding method, the metastable germanium allotropes hR8 and cI16 are shown to interconvert by means of two sets of quasi-1D chains. The first set hosts sequences of SN2 inversions of Ge tetrahedral centers, and represents the activated step. The second set does not imply any reconstruction, but assists the first one in propagating the reconstruction. The overall process is commenced by bond nucleation, followed by chain formation and reconstruction into either structure. A novel intermediate metastable phase is visited in the process. Elementary steps of chemical reactivity are accessible due to the appropriate time and spatial resolution of the methods used. This paves the way for a chemical understanding of structure reconstruction and metastable phase formation in solid materials.
Germanium has been intensively investigated, especially in recent years.3–5 Many appealing technological properties are being discovered. While homologue Si has been the classical element for electronics, Ge is again receiving increased attention for its comparatively higher intrinsic electron mobility suitable for faster electronics, for photoluminescence and band-gap tunability, for its easier integration into semiconductor devices, and for superconducting properties.4 Squeezed under moderate pressure, germanium transforms from four-connected diamond type (cF8, sp. gr. Fdm) into denser, high-symmetry structures with increasing coordination numbers .2
From mild oxidation of precursor Zintl compounds germanium can be obtained in the allotropic clathrate II structure (cF136, sp. gr. Fdm),5 which is known in host–guest compounds6 for thermoelectric applications.7 Under pressure, cF136 transforms into β-tin type (tI4, sp. gr. I41/amd) and γ-silicon type (cI16, sp. gr. Ia
) over an intermediate, with only the moderately denser structure of trigonal symmetry, hR8 (sp. gr. R
),3b which is also accessible from another metastable phase, tP12 (sp. gr. P41212).3c The hR8 Ge allotrope comprises four-coordinated Ge only, while 2 different sites Ge1 (Wyckoff position 6c) and Ge2 (Wyckoff position 18f) are indicated by the crystallographic analysis.3b This contrasts with cF8, tI4 or cI16, which can be generated from one site only.
In the extended phase diagram of germanium, hR8 is at the gateway between denser (cF8, tI4, cI16) and open-framework (cF136) structures.3b The hR8 structure may in fact be decomposed into a Ge2 framework hosting Ge1 guest chains.3b Under compression, a structural evolution towards higher symmetric structures is common. Site splitting rarely occurs and when it does so, it is the case in exotic metals like for example Ti.8 Beyond the bare crystallographic data, we interpreted this as a manifestation of chemical reactivity and we moved on to unfolding the chemically relevant steps responsible for this differentiation. Clearly, a particular rearrangement of covalent bonds into a different pattern is controlling the formation and stability of a metastable phase here.
Using state-of-art molecular dynamics techniques combined with density-functional tight-binding (DFTB) level of theory9 for expressing interatomic forces, we have investigated the mechanism connecting hR8 to cI16, two intriguing structural patterns, and group IV (Si, Ge) characteristic assemblies (Fig. 1), which appear in the Ge phase diagram as metastable phases.
![]() | ||
Fig. 1 Crystal structures of the (a) cI16 and (b) hR8 Ge allotropes (hR8 is represented in a pseudo-cubic cell for better comparison with cI16). Their networks follow the shape of the gyroid (transparent blue), except for a set of bonds, which trespass the surface in a way that is peculiar for cI16 and hR8, respectively. Flipping these bonds formally converts cI16 into hR8. One such flip is visible in the foreground, causing inversion of a Ge tetrahedron. |
The investigation of phase transition mechanisms remains a very important but challenging task.10 This is due to the specificities of the phenomenon of phase transition, which is difficult to assess both experimentally and by theory. Recent computational methods allow for shedding light onto transformation mechanisms and reactive intermediates, though.11 A detailed atomistic map of reactive events that transform one structure into the other is the commencement of a chemical understanding of phase transitions in crystalline matter, now overdue .12
Symmetry provides, to some extent, guidance thanks to group–subgroup relationships.13 However, in a situation of framework reconstruction, relevant details of bond nucleation stay undisclosed. Heuristic matching of minimal surfaces (or approximants thereof) with crystalline matter has a long record of stimulating examples.14 The networks of cI16 and hR8 develop around the gyroid,15 a surface of intrinsic Iad symmetry. All Ge atoms are found on either side of the surface, as well as the majority of covalent bonds. Only some Ge–Ge bonds trespass the gyroid, in a way that is different and characteristic for cI16 and hR8, respectively (Fig. 1). Chemical intuition suggests that the process of their mutual structural rearrangement is “economical” and well-balanced in terms of bond-breaking and bond-forming patterns. However, how does the transformation of a framework into another one actually work, and what is the role of chemistry therein, given that covalent bonds have to be re-formed? These are the questions we address in this work.
The hR8 → cI16 phase transition was investigated with TPS. The iterations were performed within the NPT ensemble (p = 12 GPa, T = 300 K) and implemented by applying momentum modifications on selected trajectory snapshots, keeping total energy, momentum and angular momentum unchanged according to the shooting scheme.19 The pressure corresponded to the equilibrium pressure of the DFTB potential, which fulfills the condition E1 + pV1 = E2 + pV2, where the indices distinguish between hR8 and cI16. Propagating the new configuration in both directions of time provided a new trajectory that was examined for the hR8–cI16 (or cI16–hR8) process, respectively. If this condition was not met, the trajectory was rejected and the shooting process iterated. The distinction between initial and final state was made based on coordination sequences (see below). MD simulations were carried out using DFTB,9 Γ-point-only module of the CP2K package.20 The accuracy of electronic density representation was checked against DFTB k-point calculations.21 The Ge DFTB energetic score is in agreement with DFT- local-density approximation (DFT-LDA) calculations (Table S2, ESI†),22 but the equilibrium volumes are slightly overestimated.
The time step t = 0.2 fs ensured a good time reversibility. The mechanistic analysis was based on more than 50 transition pathways collected after trajectory convergence. The latter was inferred from the stability of the mechanistic features emerging during TPS MD. Different from the initial, concerted mechanistic models, the stable regime shows local events of bond flipping.
![]() | ||
Fig. 2 Configurations appearing along the transformation. (a) Ge(cI16), (b) reactive intermediate, (c) a novel metastable phase Ge(hR32) and (d) Ge(hR8). Bonds changing from one structure to the other are highlighted in blue, while red ones show stretching along the trajectory, but they are never completely broken. (c) and (d) differ only by some bonds. They are nonetheless two distinct local minima on the chemical landscape of Ge. |
Striking is the formation of two sets of Ge chains around the transition reactive intermediate (Fig. 2b). These chains, which result from the condensation of isolated bonds, are tilted with an angle of ∼60° against each other and are entangled and interconnected into a web. Blue chains are oriented along [11] (referred to as the simulation box) (Fig. 2b) and are responsible for the transformation hR8 ↔ cI16 by shortening and elongating alternate sets of bonds. Red chains (along [111]), on the contrary, are forming in the course of the transition but do not rearrange from one structure into the other. Localized “red” bonds, similarly dislocated in both limiting structures, are momentarily resonating into a chain structure.
Along the transition, an intermediate, four-connected, metastable configuration is visited (hR32) (Fig. 2c), isosymmetric with hR8 but requiring certain bonds within some of the reactive [11] chains to be flipped. Details on this novel metastable configuration, quenchable to ambient pressure based on phonon calculations and finite temperature molecular dynamics simulations (DFTB/DFT), can be found in the ESI section.†
The blue-marked dumbbells in hR8 (Fig. 2d) involve Ge1 sites only, while red ones are made by Ge2 sites. Red and blue become equivalent in cI16. Changes within blue chains are thus responsible for site differentiation. Their difference is rooted in the steps around the reactive intermediate, which we are analyzing in detail now.
The bond flip mechanism among different reactive (blue) chains as well as within a single chain is markedly asynchronous. Fig. 3 shows a close-up of the bond reorganization of a single chain as a function of time. Noteworthy is the formation of trigonal-bipyramidal geometries at Ge, suggesting sp2 + p hybridized states like they are caused by a SN2 reaction sequence mechanism, in which tetrahedral sp3-hybridized Ge atoms get inverted. Most notable is the generation of an almost regularly spaced sequence of sp2 Ge atoms. This sequence undergoes rapid rearrangements with the formation of new bonds, characteristic either of Ge(hR8) or Ge(cI16), respectively. Single bonds “nucleate” from a roughly equally spaced arrangement, triggering a sequence of bond-forming and bond-breaking events.
![]() | ||
Fig. 3 Transition mechanism for a single reaction chain (in blue in Fig. 1). Snapshots of the trajectory are shown at 0.0 (a), 880 (b), 906 (c), 960 (d), 1040 (e), 2000 (f) fs. |
A global characterization of all bond changes along a trajectory is achieved by calculating Wannier functions centers (WFCs).24 Two atoms are associated with a single WFC (within a so-called “association length”) in the presence of a chemical bond. Upon bond breaking, the WFC disappears or becomes shifted. For the characterization of the transition trajectories with respect to chemical bond changes, the motion of WFCs was monitored along the time coordinate of framework reconstruction. The shift of a bond within a chain (Fig. 3) is reflected into a jump in WFC relative position (Fig. 4). Just one snapshot every 20 fs was considered near the relaxed hR8 and cI16 phases. In the intermediate region, a finer time resolution (every 2 fs) was necessary. The displacement of Wannier centers serves here as collective order parameter.
![]() | ||
Fig. 4 Evolution of WFC relative distances, calculated for each bond and each reactive chains (four blue chains in Fig. 2b) along the time coordinate of the framework reconstruction transition. Different colors (green, black, purple and turquoise) refer to distinct chains, while distinct symbols (square, circle and triangle) refer to different bonds within a single chain. Inset: potential energy profile with hR32 as the intermediate phase. Zeroth time frame is arbitrary set. The y axis corresponds to relative displacements of WFC along the transition and is used here as an order parameter, able to distinguish between initial and final states. |
Fig. 4 shows the time evolution of WFCs within every reactive chain. Different reactive chains are distinguished by color (green, black, purple and turquoise), while bonds within a chain are marked by symbols. Clearly, the total reconstructive process is discontinuous on the timescale of the simulation. Chain bond flippings occur at different moments in the whole material, as indicated by the mutual offset of the sigmoidal curves. For each color, different symbols jump at different times, denoting non-concerted bond shifts within chains, as already highlighted in the close-up of Fig. 3. Therefore, the overall bond rearrangement process spans quite a large time window: the first bond flip occurs at ∼750 fs while the last one is at ∼1200 fs. The green curves refer to the chain of Fig. 3. As the first bond gets broken (green circle), the other bonds in the same chain start rearranging (green square, triangle), followed after some delay by bond shifts in neighboring chains (black, purple and turquoise). The larger WFC fluctuations around 850 fs (Fig. 4, green chain) and 1000 fs (green and purple chains) coincides with the formation of red chains, which is taking place just before breakage/formation of reactive (black and purple) chains, and less markedly ahead of the last chain reshuffling (∼1200 fs).
This is echoed in the changes of phonon frequencies along the transition. The stretching mode of the first flipping bond becomes negative at t = ∼750 fs, at the onset of bond nucleation (green curve and circles in Fig. 4). At ∼800 fs also, the frequency associated with black squares flips and turns negative, rapidly followed by other reactive chain modes. At 950 fs, in correspondence with the first fluctuation of the green curve (circle), the modes associated with the red, non-reactive chains become unstable and keep showing negative frequencies till 1050 fs, past the second fluctuation value of green and violet curves. The maximal wavenumbers of the red chains (−75 cm−1) are considerably smaller than the blue chain modes (−270 cm−1). Nonetheless they strongly influence phase transition by spreading the initial instability associated with single bond nucleation events over the whole material via chemical bonds.
The inset of Fig. 4 shows the potential energy profile along the reaction coordinate for a representative trajectory. For the hR8 → cI16 transformation, hR32 is initially formed, which is structurally and energetically similar to hR8, before a major reconstruction towards cI16 takes place. This allows splitting of the overall reconstruction into two “elementary” steps, hR8 → hR32, followed by hR32 → cI16. The intermediate can be formed thanks to the asynchrony of chain reconstruction, as hR8 and hR32 differ by just a few bonds. To the extent that this solid transformation can indeed be understood as a set of “molecular” steps, the Hammond postulate,27 may warrant some guidance. The hR32 intermediate is structurally and energetically closer to hR8 and is formed first. The second step involves more bond reshuffling and a larger activation barrier, as it contains the core of the reconstructions. On forming hR32 first, hR8 would be the kinetic product, cI16 the thermodynamic one.
To visualize the sp3 → (sp2 + p) → sp3 sequence of re-hybridization around Ge inside reactive chains, Wannier functions were calculated, Fig. 5. In cI16 (t = 0 fs) four-bonded Ge shows four equal sp3 Wannier functions. Around t = 880 fs, the first bond is broken, under formation of a trigonal planar center (3 Wannier functions) and of a lone pair along the chain (contracted Wannier function, associated with one Ge only). Geometry fluctuations cause the lone pair to slightly expand and contract again (t = 880–960 fs) while the intermediate arrangement is maintained. Eventually, a new bond is formed under inversion of the pristine tetrahedral arrangements at t = 960–1040 fs. This is indicated by a rapid change of spread and shape of the corresponding Wannier function along the chain axis. Minor rearrangements lead to hR8 at t = 2000 fs.
![]() | ||
Fig. 5 sp3-projected, maximally-localized Wannier functions for the Ge atom at which the transition is commenced. The evolution of bonds (transparent lobes) and lone pair (opaque lobe) are shown at different snapshots. For clarity, only the positive lobe of the MLWFs is rendered. The chain and time frames are the same as in Fig. 3. Zeroth time frame was set at cI16. |
Quasi-regularly spaced chains appeared in the intermediate regions of the phase transition reconstruct by symmetrically pairing atoms into alternating localized bonds in a sort of Bond Density Wave (BDW),28 reminiscent of a 1D Peierls distortion. Nonetheless, our analysis has shown that different modes coexist while the networks reconstruct. Particularly, no chain acts as a truly isolated 1D system, as a non-reactive “red” mode is supporting bond rearrangement. While the rearrangement of a single chain may be termed “Peierls-patterned”, a 1D chain is not acting here as the fundamental instability controlling the process.
The effect of periodic boundary conditions needs careful consideration when studying nucleation geometries. However, we believe our results to be transferable to larger systems, as the initial “nucleus” was determined to be a single bond switching, which can occur anywhere within a reactive chain. Furthermore, the non-concerted character of the bond shuffling process allows for the formation of an intermediate metastable structure (hR32), which we expect to be detectable in high-pressure experiments with an appropriate (de)compression protocol. This would represent a rare case of experimental verifiability of a solid–solid phase transition mechanism, and a major step towards controlling reactivity in the solid state.
In conclusion, by atomistic simulations we have characterized the hR8 to cI16 phase transition as a chemical isomerization reaction, which is taking place in the solid state. Sets of reactive chains are formed in the intermediate transition region. Within such reactive channels, distinct sets of bonds can rearrange from one structural pattern into the other one. Framework reconstruction is commenced by local bond nucleation, which is propagated along chains. The propagation step is assisted by a second set of chains, which is forming around the transition state only, but which is not accommodating any event of bond breakage or bond formation. In the reactive chains, bonds are reshuffled by Ge sp3 tetrahedra inversions, with the formation of intermediate trigonal-planar geometries known from SN2 type reactions.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c2ra20837g |
This journal is © The Royal Society of Chemistry 2012 |