Atomistic simulations of calcium aluminosilicate interfaced with liquid water.

The dissolution behavior of calcium aluminosilicate based glass fibers, such as stone wool fibers, is an important consideration in mineral wool applications for both the longevity of the mineral wool products in humid environments and limiting the health impacts of released and inhaled fibers from the mineral wool product. Balancing these factors requires a molecular-level understanding of calcium aluminosilicate glass dissolution mechanisms, details that are challenging to resolve with experiment alone. Molecular dynamics simulations are a powerful tool capable of providing complementary atomistic insights regarding dissolution; however, they require force fields capable of describing not-only the calcium aluminosilicate surface structure but also the interactions relevant to dissolution phenomena. Here, a new force field capable of describing amorphous calcium aluminosilicate surfaces interfaced with liquid water is developed by fitting parameters to experimental and first principles simulation data of the relevant oxide-water interfaces, including ab initio molecular dynamics simulations performed for this work for the wüstite and periclase interfaces. Simulations of a calcium aluminosilicate surface interfaced with liquid water were used to test this new force field, suggesting moderate ingress of water into the porous glass interface. This design of the force field opens a new avenue for the further study of calcium and network-modifier dissolution phenomena in calcium aluminosilicate glasses and stone wool fibers at liquid water interfaces.

Wetting Properties of the CO2-Water-Calcite System via Molecular Simulations: Shape and Size Effects.

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO2 and the pore fluids. In this context, the interfacial tension (IFT) between CO2 and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO2-water IFT and the CO2-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO2 model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO2/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO2 atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO2. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.

Investigating the polymorphism in PR179: a combined crystal structure prediction and metadynamics study.

The polymorphism of an industrial important pigment (PR179) was studied with a combination of standard crystal structure prediction and metadynamics. The former provided a starting set of candidate polymorphs whose structural and thermal stability were then probed by metadynamics. Moreover, metadynamics allowed for exploring the free energy surface to look for other possible polymorphs that were not included in the original set. Our calculations indicate that two structures have a high structural stability and are therefore good candidates to be found in experiments. The lower energy phase of the two indeed corresponds to the known polymorph, and we suggest that the other might be a metastable polymorph not yet experimentally discovered.

Morphological and compositional evolution of the ge/si(001) surface during exposure to a si flux.

By using scanning tunneling microscopy we found that the surface reconstruction of Ge/Si(001) epilayers evolves from (M x N) to (2 x N), and eventually to (2 x 1), during exposure to a Si flux. This sequence appears to be just the inverse of that observed during the growth of Ge or SiGe alloys on Si(001). However, molecular dynamics simulations supported by ab initio calculations allow us to interpret this morphological evolution in terms of Si migration through the epilayer and complex Si-Ge intermixing below the top Ge layer.

Critical role of the surface reconstruction in the thermodynamic stability of (105) Ge pyramids on Si(001).

We show by molecular dynamics simulations on a scale comparable to experimental dimensions that a peculiar surface reconstruction of the (105) facets is responsible for the stability of Ge pyramids on Si(001). This finding is confirmed by ab initio total energy calculations for competing surface reconstructions and a very satisfactory comparison of the corresponding charge density maps to scanning tunneling microscopy measurements of the facets, both for filled and empty states.

Atomic-scale pathway of the pyramid-to-dome transition during ge growth on Si(001).

By high resolution scanning tunneling microscopy, we investigate the morphological transition from pyramid to dome islands during the growth of Ge on Si(001). We show that pyramids grow from top to bottom and that, from a critical size on, incomplete facets are formed. We demonstrate that the bunching of the steps delimiting these facets evolves into the steeper dome facets. Based on first principles and Tersoff-potential calculations, we develop a microscopic model for the onset of the morphological transition, able to reproduce closely the experimentally observed behavior.