Andrew G. Stack

Bio: Andrew G. Stack is an academic researcher from Oak Ridge National Laboratory. The author has contributed to research in topics: Nucleation & Dissolution. The author has an hindex of 26, co-authored 77 publications receiving 1859 citations. Previous affiliations of Andrew G. Stack include Georgia Institute of Technology & Lawrence Livermore National Laboratory.

Accurate Rates of the Complex Mechanisms for Growth and Dissolution of Minerals Using a Combination of Rare-Event Theories

Abstract: Mineral growth and dissolution are often treated as occurring via a single reversible process that governs the rate of reaction. We show that multiple distinct intermediate states can occur during both growth and dissolution. Specifically, we used metadynamics, a method for efficiently exploring the free-energy landscape of a system, coupled to umbrella sampling and reactive flux calculations to examine the mechanism and rates of attachment and detachment of a barium ion onto a stepped barite (BaSO4) surface. The activation energies calculated for the rate-limiting reactions, which are different for attachment and detachment, precisely match those measured experimentally during both growth and dissolution. These results can potentially explain anomalous non-steady-state mineral reaction rates observed experimentally and will enable the design of more efficient growth inhibitors and facilitate an understanding of the effect of impurities.

Growth Rate of Calcite Steps As a Function of Aqueous Calcium-to-Carbonate Ratio: Independent Attachment and Detachment of Calcium and Carbonate Ions

Abstract: Growth rates of monolayer-height steps on the {1014} calcite surface have been measured as a function of the aqueous calcium-to-carbonate ratio. The maximum growth rates of the two common crystallographic orientations were found to deviate from the ideal stoichiometric ratio of 1:1, and dissolution features were observed under supersaturated solutions containing high calcium-to-carbonate ratios. To explain these phenomena, a theory is applied that treats the rates of attachment and detachment of aqueous calcium and carbonate ions separately. The resultant attachment rate constants are 1−3 orders of magnitude smaller than the water exchange rate of the constituent aqueous ions, suggesting that ligand-exchange processes may directly drive attachment. The broader implication is that the saturation state alone is not adequate to fully describe the rates of the multiple, independent reactions that occur on mineral surfaces under these conditions.

CO2 sorption to subsingle hydration layer montmorillonite clay studied by excess sorption and neutron diffraction measurements.

Abstract: Geologic storage of CO2 requires that the caprock sealing the storage rock is highly impermeable to CO2. Swelling clays, which are important components of caprocks, may interact with CO2 leading to volume change and potentially impacting the seal quality. The interactions of supercritical (sc) CO2 with Na saturated montmorillonite clay containing a subsingle layer of water in the interlayer region have been studied by sorption and neutron diffraction techniques. The excess sorption isotherms show maxima at bulk CO2 densities of ≈0.15 g/cm3, followed by an approximately linear decrease of excess sorption to zero and negative values with increasing CO2 bulk density. Neutron diffraction experiments on the same clay sample measured interlayer spacing and composition. The results show that limited amounts of CO2 are sorbed into the interlayer region, leading to depression of the interlayer peak intensity and an increase of the d(001) spacing by ca. 0.5 A. The density of CO2 in the clay pores is relatively stab...

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

The ReaxFF reactive force-field: development, applications and future directions

Abstract: The reactive force-field (ReaxFF) interatomic potential is a powerful computational tool for exploring, developing and optimizing material properties. Methods based on the principles of quantum mechanics (QM), while offering valuable theoretical guidance at the electronic level, are often too computationally intense for simulations that consider the full dynamic evolution of a system. Alternatively, empirical interatomic potentials that are based on classical principles require significantly fewer computational resources, which enables simulations to better describe dynamic processes over longer timeframes and on larger scales. Such methods, however, typically require a predefined connectivity between atoms, precluding simulations that involve reactive events. The ReaxFF method was developed to help bridge this gap. Approaching the gap from the classical side, ReaxFF casts the empirical interatomic potential within a bond-order formalism, thus implicitly describing chemical bonding without expensive QM calculations. This article provides an overview of the development, application, and future directions of the ReaxFF method.

Observation of Single Colloidal Platinum Nanocrystal Growth Trajectories

Abstract: It is conventionally assumed that the growth of monodisperse colloidal nanocrystals requires a temporally discrete nucleation followed by monomer attachment onto the existing nuclei. However, recent studies have reported violations of this classical growth model, and have suggested that inter-particle interactions are also involved during the growth. Mechanisms of nanocrystal growth still remain controversial. Using in situ transmission electron microscopy, we show that platinum nanocrystals can grow either by monomer attachment from solution onto the existing particles or by coalescence between the particles. Surprisingly, an initially broad size distribution of the nanocrystals can spontaneously narrow. We suggest that nanocrystals take different pathways of growth based on their size- and morphology-dependent internal energies. These observations are expected to be highly relevant for other nanocrystal systems.