Search NASA⌕ Search

SEARCH · Search NASA

Results for “set of small ordered structures”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Deep learning simulations of the microwave sky

Here we present 500 high-resolution, full-sky millimeter-wave deep learning (DL) simulations that include lensed CMB maps and correlated foreground components. We find that these MillimeterDL simulations can reproduce a wide range of non-Gaussian summary statistics matching the input training simulations, while only being optimized to match the power spectra. The procedure we develop in this work enables the capability to mass produce independent full-sky realizations from a single expensive full-sky simulation, when ordinarily the latter would not provide enough training data. We also circumvent a common limitation of high-resolution DL simulations that they be confined to small sky areas, often due to memory or GPU issues; we do this by developing a “stitching” procedure that can recover the large-scale, high-order statistics and avoid discontinuities or repeated features in the maps. In addition, since our network takes as input a full-sky lensing convergence map, it can in principle take a full-sky lensing convergence map from any large-scale structure (LSS) simulation and generate the corresponding lensed CMB and correlated foreground components at millimeter wavelengths; this is especially useful in the current era of combining results from both CMB and LSS surveys, which require a common set of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Space–time structure of weak magnetohydrodynamic turbulence

The two-time energy spectrum of weak magnetohydrodynamic turbulence is found by applying a wave-turbulence closure to the cumulant hierarchy constructed from the dynamical equations. Solutions are facilitated via asymptotic expansions in terms of the small parameter $\varepsilon$ , describing the ratio of time scales corresponding to Alfvénic propagation and nonlinear interactions between counter-propagating Alfvén waves. The strength of nonlinearity at a given spatial scale is further quantified by an integration over all possible delta-correlated modes compliant in a given set of three-wave interactions that are associated with energy flux through the said scale. The wave-turbulence closure for the two-time spectrum uncovers a secularity occurring on a time scale of order $\varepsilon ^{-2}$ , and the asymptotic expansion for the spectrum is reordered in a manner comparable to the one-time case. It is shown that for the regime of stationary turbulence, the two-time energy spectrum exponentially decays on a lagged time scale $(\varepsilon ^2 \gamma _k^s)^{-1}$ in proportion to the strength of the associated three-wave interactions, characterized by nonlinear decorrelation frequency $\gamma _k^s$ . The scaling of the form $k_{\perp } v_0 \chi _0$ exhibited by this frequency is reminiscent of random sweeping by the outer scale with characteristic fluctuation velocity $v_0$ that is modified due to competition with Alfvénic propagation (characterized by $\chi _0$ ) at the said scale. A brief calculation of frequency broadening of the power spectrum due to nonlinear interactions is also presented.

Physics↗
DOE OSTI2021

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗
DOE OSTI2019

Renormalization group computation of likelihood functions for cosmological data sets

I show how a renormalization group (RG) method can be used to incrementally integrate the information in cosmological large-scale structure data sets (including CMB, galaxy redshift surveys, etc.). I show numerical tests for Gaussian fields, where the method allows arbitrarily close to exact computation of the likelihood function in order N time, even for problems with no symmetry, compared to N 3 for brute force linear algebra (where N is the number of data points - to be fair, methods already exist to solve the Gaussian problem in at worst N log N time, and this method will not necessarily be faster in practice). The method requires no sampling or other Monte Carlo (random) element. Nonlinearity/non-Gaussianity can be accounted for to the extent that terms generated by integrating out small scale modes can be projected onto a sufficient basis, e.g., at least in the sufficiently perturbative regime. The formulas to evaluate are straightforward and require no understanding of quantum field theory, but this paper may also serve as a pedagogical introduction to Wilsonian RG for astronomers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗
DOE OSTI2015

Homoclinic behaviors and chaotic motions of double layered viscoelastic nanoplates based on nonlocal theory and extended Melnikov method

The nonlinear dynamical equations are established for the double layered viscoelastic nanoplates (DLNP) subjected to in-plane excitation based on the nonlocal theory and von Kármán large deformation theory. The extended high dimensional homoclinic Melnikov method is employed to study the homoclinic phenomena and chaotic motions for the parametrically excited DLNP system. The criteria for the homoclinic transverse intersection for both the asynchronous and synchronous buckling cases are proposed. Lyapunov exponents and phase portraits are obtained to verify the Melnikov-type analysis. The influences of structural parameters on the transverse homoclinic orbits and homoclinic bifurcation sets are discussed for the two buckling cases. Some novel phenomena are observed in the investigation. It should be noticed that the nonlocal effect on the homoclinic behaviors and chaotic motions is quite remarkable. Hence, the small scale effect should be taken into account for homoclinic and chaotic analysis for nanostructures. It is significant that the nonlocal effect on the homoclinic phenomena for the asynchronous buckling case is quite different from that for the synchronous buckling case. Moreover, due to the van der Walls interaction between the layers, the nonlocal effect on the homoclinic behaviors and chaotic motions for high order mode is rather tiny under the asynchronous buckling condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The effect of interstellar absorption on measurements of the baryon acoustic peak in the Lyman α forest

In recent years, the autocorrelation of the hydrogen Lyman α forest has been used to observe the baryon acoustic peak at redshift 2 < z < 3.5 using tens of thousands of QSO spectra from the BOSS survey. However, the interstellar medium of the Milky Way introduces absorption lines into the spectrum of any extragalactic source. These lines, while weak and undetectable in a single BOSS spectrum, could potentially bias the cosmological signal. In order to examine this, we generate absorption line maps by stacking over a million spectra of galaxies and QSOs. Here, we find that the systematics introduced are too small to affect the current accuracy of the baryon acoustic peak, but might be relevant to future surveys such as the Dark Energy Spectroscopic Instrument (DESI). We outline a method to account for this with future data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

CMB-S4: Iterative Internal Delensing and r Constraints

Abstract The tightest constraints on the tensor-to-scalar ratio r can only be obtained after removing a substantial fraction of the lensing B -mode sample variance. The planned Cosmic Microwave Background (CMB)-S4 experiment ( cmb-s4.org ) will remove the lensing B -mode signal internally by reconstructing the gravitational lenses from high-resolution observations. We document here a first lensing reconstruction pipeline able to achieve this optimally for arbitrary sky coverage. We make it part of a map-based framework to test CMB-S4 delensing performance and its constraining power on r , including inhomogeneous noise and two non-Gaussian Galactic polarized foreground models. The framework performs component separation of the high-resolution maps, followed by the construction of lensing B -mode templates, which are then included in a parametric small-aperture map cross-spectra-based likelihood for r . We find that the lensing reconstruction and framework achieve the expected performance, compatible with the target σ ( r ) ≃ 5 · 10 −4 in the absence of a tensor signal, after an effective removal of 92%–93% of the lensing B -mode variance, depending on the simulation set. The code for the lensing reconstruction can also be used for cross-correlation studies with large-scale structures, lensing spectrum reconstruction, cluster lensing, or other CMB lensing-related purposes. As part of our tests, we also demonstrate the joint optimal reconstruction of the lensing potential with the lensing curl potential mode at second order in the density fluctuations.

Astronomy & Astrophysics↗
DOE OSTI2018

Sodium Pumping via Condensation within a Non-Wetting Porous Structure

In a sodium thermal electrochemical converter (Na-TEC), the transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is necessary to complete the thermodynamic cycle. This pumping can be typically achieved via capillary action. A unique capillary pump is explored in this work, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > π/2) stainless steel porous structure. Due to the curvature of the non-wetting liquid-vapor interface, the liquid adjacent to the interface is at a higher pressure than the vapor. This is in contrast to traditional wicks, where the liquid adjacent to the interface has a lower pressure. Conjugate heat transfer techniques are required to model the coupled momentum and thermal transport processes within this porous structure. Due to the low mass flowrates in the Na-TEC (< 0.1 mg/s), higher order friction terms are neglected and the classical Darcy’s law is applied for the non-isothermal, compressible vapor flow. At lower temperatures, the sodium vapor density becomes sufficiently small and there is a transition from classic Poiseuille flow (Kn < 0.1) to Knudsen flow (Kn > 10). The permeability of the porous structure is modeled using the Blaze-Kozeny expression, with an additional factor that accounts for this change in Kn. Initial modeling results showing expected mass flowrates, temperature profiles, and liquid-vapor interface locations are presented. The design of an experiment to validate the results from this modeling effort is also described. The main design features of the experimental set-up are presented, and the method used to measure the low sodium mass flowrates is discussed. Initial experimental results demonstrating sodium pumping for a range of temperatures (850K – 1150K) will be presented and compared against the model predictions. A number of challenges affecting the successful completion of this experiment will also be discussed. Primarily, the sodium vapor must be superheated to prevent condensation anywhere except within the porous structure. Other challenges include the requirement to keep the condensation interface above ~ 675K to prevent transition to a wetting regime, and the control of the transient heating rate so that vapor flows slower than the quasi steady-state mass flowrate, which would otherwise lead to a flooding condition that prevents pumping. Specific strategies used to address these challenges are presented, and the experimental procedures used to gather data are discussed.

14 SOLAR ENERGY↗
DOE OSTI2019

Deep Learning for Semantic Segmentation of Defects in Advanced STEM Images of Steels

Crystalline materials exhibit long-range ordered lattice unit, within which resides nonperiodic structural features called defects. These crystallographic defects play a vital role in determining the physical and mechanical properties of a wide range of material systems. While computer vision has demonstrated success in recognizing feature patterns in images with well defined contrast, automated identification of nanometer scale crystallographic defects in electron micrographs governed by complex contrast mechanism is still a challenging task. Here, building upon an advanced defect imaging mode that offers high feature clarity, we introduce DefectSegNet - a new convolutional neural network (CNN) architecture that performs semantic segmentation of three common crystallographic defects in structural alloys: dislocation lines, precipitates and voids. Results from supervised training of a small set of high-quality defect images of steels show high pixel-wise accuracy across all three types of defects: 91.60 ± 1.77% on dislocations, 93.39 ± 1.00% on precipitates, and 98.85 ± 0.56% on voids. We discuss the sources of uncertainties in CNN prediction, and the training data in terms of feature density, representation and homogeneity and their effects on deep learning performance. Further defect quantification using DefectSegNet prediction outperforms human expert average, presenting a promising new workflow for fast and statistically meaningful quantification of materials defects.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Electronic Energies Are Not Enough: An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H + GPGG Conformers

Peptides and proteins exist not as single, lowest energy structures but as ensembles of states separated by small barriers. In order to study these species we must be able to correctly identify their gas phase conformational distributions, and ion mobility spectrometry (IMS) has arisen as an experimental method for assessing the gas phase energetics of flexible peptides. Here, we present a thorough exploration and benchmarking of the low energy conformers of the small, hairpin peptide H+ GPGG with the aid of ion mobility spectrometry against a wide swath of density functionals (35 dispersion-corrected and uncorrected functionals represented by 21 unique exchange–correlation functionals) and wave function theory methods (15 total levels of theory). The three experimentally resolved IMS peaks were found to correspond to three distinct pairs of conformers, each pair composed of species differing only by the chair or boat configuration of the proline. Two of the H+ GPGG conformer pairs possess a cis configuration about the Pro–Gly1 peptide bond while the other adopts a trans configuration. While the experimental spectrum reports a higher intensity for the cis-1 conformer than the trans conformer, 13 WFT levels of theory, including a complete basis set CCSD(T) extrapolation, obtain trans to be favored in terms of the electronic energy. This same effect is seen in 14 of the 18 dispersion-corrected density functionals studied, whereas the remaining 17 functionals show more variety. Only when Gibbs free energies are considered do the WFT methods and dispersion-corrected functionals reflect the experimental distribution. CAM-B3LYP-D3BJ emerges as the best-performing density functional, matching the experimental distribution and the CCSD(T)/CBS relative energies within 6%. Further analysis reveals the trans conformer to be favored electronically, but entropically disfavored, leading to the experimental preference of the cis-1 conformer. Furthermore, these results highlight the danger in considering only the electronic energies, which is common practice in electronic structure theory predictions of conformational energy distributions. Additionally, the effects of temperature and scaling of the frequencies used in obtaining the Gibbs free energies are explored.

Energy↗
DOE OSTI2015

Detecting vortices in superconductors: Extracting one-dimensional topological singularities from a discretized complex scalar field

In type II superconductors, the dynamics of superconducting vortices determine their transport properties. In the Ginzburg-Landau theory, vortices correspond to topological defects in the complex order parameter. Extracting their precise positions and motion from discretized numerical simulation data is an important, but challenging, task. In the past, vortices have mostly been detected by analyzing the magnitude of the complex scalar field representing the order parameter and visualized by corresponding contour plots and isosurfaces. However, these methods, primarily used for small-scale simulations, blur the fine details of the vortices, scale poorly to large-scale simulations, and do not easily enable isolating and tracking individual vortices. In this paper, we present a method for exactly finding the vortex core lines from a complex order parameter field. With this method, vortices can be easily described at a resolution even finer than the mesh itself. The precise determination of the vortex cores allows the interplay of the vortices inside a model superconductor to be visualized in higher resolution than has previously been possible. Finally, by representing the field as the set of vortices, this method also massively reduces the data footprint of the simulations and provides the data structures for further analysis and feature tracking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Surface Potentials of (001), (012), (113) Hematite (α-Fe2O3) Crystal Faces in Aqueous Solution

Hematite (α-Fe2O3) is an important candidate electrode for energy system technologies such as photoelectrochemical water splitting. Conversion efficiency issues with this material are presently being addressed through nanostructuring, doping, and surface modification. However, key electrochemical properties of hematite/electrolyte interfaces remain poorly understood at a fundamental level, in particular those of crystallographically well-defined hematite faces likely present as interfacial components at the grain scale. We report a combined measurement and theory study that isolates and evaluates the equilibrium surface potentials of three nearly defect-free single crystal faces of hematite, titrated from pH 3 to 11.25. We link measured surface potentials with atomic-scale surface topology, namely the ratio and distributions of surface protonation/deprotonation site types expected from the bulk structure. The data reveal face-specific points of zero potential (PZP) relatable to points of zero net charge (PZC) that lie within a small pH window (8.35-8.85). Over the entire pH range the surface potentials show strong non-Nernstian charging at pH extremes separated by a wide central plateau in agreement with surface complexation modeling predictions, but with important face-specific distinctions. We introduce a new surface complexation model based on fitting the entire data set that depends primarily only on the proton affinities of two site types and the two associated electrical double layer capacitances. The data and model show that magnitudes of surface potential biases at the pH extremes are on the order of 100 mV, similar to the activation energy for electron hopping mobility. An energy band diagram for hematite crystallites with specific face expression and pH effects is proposed that could provide a baseline for understanding water splitting performance enhancement effects from nanostructuring, and guide morphology targets and pH for systematic improvements in efficiency.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Sodium Pumping via Condensation within a Non-Wetting Porous Structure

The transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is required to complete the thermodynamic cycle of a sodium thermal electrochemical converter (Na-TEC). A unique capillary pump is explored in this work, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > pi/2) stainless steel porous structure. Due to the curvature of the non-wetting liquid-vapor interface, the liquid adjacent to the interface is at a higher pressure than the vapor. This is in contrast to traditional wicks, where the liquid adjacent to the interface has a lower pressure. Conjugate heat transfer techniques are used to model the coupled momentum and thermal transport processes within this non-wetting porous structure. Due to the low mass flowrates in the Na-TEC (< 0.1 mg/s), higher order friction terms are neglected and the classical Darcy’s law is used for the non-isothermal, compressible vapor flow. At lower temperatures, the sodium vapor density becomes sufficiently small and there is a transition from classic Poiseuille flow (Kn < 0.1) to Knudsen flow (Kn > 10). Furthermore, transport due to binary diffusion is also considered for situations where there is a substantial presence of non-condensable argon in the sodium vapor plenum. Initial modeling results showing expected mass flowrates, temperature profiles, and liquid-vapor interface locations are presented for various physical conditions. The main design features of an experimental set-up used to demonstrate this capillary pumping is also described in detail, including the method used to measure the low sodium mass flowrates. Initial experimental results demonstrating sodium pumping for a range of temperatures (800K – 1100K) will be presented and compared against the model predictions. Finally, several of the challenges with this experiment will be discussed, and specific strategies used to address these are presented.

14 SOLAR ENERGY↗
DOE OSTI2021

Data Mining of Polymer Phase Transitions upon Temperature Changes by Small and Wide-Angle X-ray Scattering Combined with Raman Spectroscopy

The complex physical transformations of polymers upon external thermodynamic changes are related to the molecular length of the polymer and its associated multifaceted energetic balance. The understanding of subtle transitions or multistep phase transformation requires real-time phenomenological studies using a multi-technique approach that covers several length-scales and chemical states. A combination of X-ray scattering techniques with Raman spectroscopy and Differential Scanning Calorimetry was conducted to correlate the structural changes from the conformational chain to the polymer crystal and mesoscale organization. Current research applications and the experimental combination of Raman spectroscopy with simultaneous SAXS/WAXS measurements coupled to a DSC is discussed. In particular, we show that in order to obtain the maximum benefit from simultaneously obtained high-quality data sets from different techniques, one should look beyond traditional analysis techniques and instead apply multivariate analysis. Data mining strategies can be applied to develop methods to control polymer processing in an industrial context. Crystallization studies of a PVDF blend with a fluoroelastomer, known to feature complex phase transitions, were used to validate the combined approach and further analyzed by MVA.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Direct dynamics simulation of dioxetane formation and decomposition via the singlet ·O–O–CH2–CH2 · biradical: Non-RRKM dynamics

Electronic structure calculations and direct chemical dynamics simulations are used to study the formation and decomposition of dioxetane on its ground state singlet potential energy surface. The stationary points for 1O2 + C2H4, the singlet {center_dot}O-O-CH2-CH2{center_dot} biradical, the transition state (TS) connecting this biradical with dioxetane, and the two transition states and gauche {center_dot}O-CH2-CH2-O{center_dot} biradical connecting dioxetane with the formaldehyde product molecules are investigated at different levels of electronic structure theory including UB3LYP, UMP2, MRMP2, and CASSCF and a range of basis sets. The UB3LYP/6-31G* method was found to give representative energies for the reactive system and was used as a model for the simulations. UB3LYP/6-31G* direct dynamics trajectories were initiated at the TS connecting the {center_dot}O-O-CH2-CH2{center_dot} biradical and dioxetane by sampling the TS's vibrational energy levels, and rotational and reaction coordinate energies, with Boltzmann distributions at 300, 1000, and 1500 K. This corresponds to the transition state theory model for trajectories that pass the TS. The trajectories were directed randomly towards both the biradical and dioxetane. A small fraction of the trajectories directed towards the biradical recrossed the TS and formed dioxetane. The remainder formed 1O2 + C2H4 and of these {approx} 40% went directly from the TS to 1O2 + C2H4 without getting trapped and forming an intermediate in the {center_dot}O-O- CH2-CH2{center_dot} biradical potential energy minimum, a non-statistical result. The dioxetane molecules which are formed dissociate to two formaldehyde molecules with a rate constant two orders of magnitude smaller than that predicted by Rice-Ramsperger-Kassel-Marcus theory. The reaction dynamics from dioxetane to the formaldehyde molecules do not follow the intrinsic reaction coordinate or involve trapping in the gauche {center_dot}O-CH2-CH2-O{center_dot} biradical potential energy minimum. Important non-statistical dynamics are exhibited for this reactive system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Cross-nucleation between clathrate hydrate polymorphs: Assessing the role of stability, growth rate, and structure matching

Cross-nucleation is a phenomenon where a new crystal nucleates and grows upon the surface of a different polymorph. Previous studies indicate that faster growth rate of the new crystal is a necessary but not sufficient condition for cross-nucleation. The thermodynamic stability of the different polymorphs can also affect cross-nucleation by modulating the rates of crystal growth. The interplay between thermodynamic stability of the polymorphs involved, the growth rate of the crystals, and the need for creation of an interfacial transition layer that seamlessly connects the two structures has not yet been fully elucidated. Predicting cross-nucleation is particularly challenging for clathrate hydrates, for which there are sometimes several polymorphs with similar stability and for which growth rates are not known. In this work, we use molecular dynamics simulations to investigate which factor (stability, growth rate, or formation of interfacial transition layer) controls cross-nucleation between the four known Frank-Kasper clathrate hydrate polymorphs: sI, sII, TS, and HS-I. We investigate the growth and cross-nucleation of these four hydrates filled with a set of guest molecules that produce different order of stabilities for the four crystal structures. We determine that the growth rate of sII clathrate is the fastest, followed by TS, HS-I, and sI. We find that cross-nucleation into or from sII clathrates is preceded by the formation of an interfacial transition layer at the seed crystal/liquid interface because sII does not share a crystal plane with sI, HS-I, or TS. Cross-nucleation between the latter three can occur seamlessly and is determined only by their growth rates. Our results indicate that nucleation of an interfacial transition layer between non-matching polymorphs can control cross-nucleation or lack thereof under conditions of small driving force. Under conditions of sufficient supercooling clathrate hydrate polymorphs cross-nucleate into the fastest growing phase even if that new phase is less stable and does not share a common crystal plane with the initial polymorph.

36 MATERIALS SCIENCE↗