Next Article in Journal
Fly Ash-Added, Seawater-Mixed Pervious Concrete: Compressive Strength, Permeability, and Phosphorus Removal
Next Article in Special Issue
Dissolution of Portlandite in Pure Water: Part 2 Atomistic Kinetic Monte Carlo (KMC) Approach
Previous Article in Journal
Hydrazine Hydrate Induced Three-Dimensional Interconnected Porous Flower-like 3D-NiCo-SDBS-LDH Microspheres for High-Performance Supercapacitor
Previous Article in Special Issue
Dynamic Characteristics and Damage Constitutive Model of Mudstone under Impact Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach

by
Khondakar Mohammad Salah Uddin
1,*,
Mohammadreza Izadifar
2,
Neven Ukrainczyk
2,*,
Eduardus Koenders
2 and
Bernhard Middendorf
1
1
Department of Structural Materials and Construction Chemistry, University of Kassel, Mönchebergstraße 7, 34125 Kassel, Germany
2
Institute of Construction and Building Materials, Technical University of Darmstadt, Franziska-Braun-Str. 3, 64287 Darmstadt, Germany
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(4), 1404; https://doi.org/10.3390/ma15041404
Submission received: 11 January 2022 / Revised: 3 February 2022 / Accepted: 9 February 2022 / Published: 14 February 2022
(This article belongs to the Special Issue Mathematical Modeling of Building Materials)

Abstract

:
The current contribution proposes a multi-scale bridging modeling approach for the dissolution of crystals to connect the atomistic scale to the (sub-) micro-scale. This is demonstrated in the example of dissolution of portlandite, as a relatively simple benchmarking example for cementitious materials. Moreover, dissolution kinetics is also important for other industrial processes, e.g., acid gas absorption and pH control. In this work, the biased molecular dynamics (metadynamics) coupled with reactive force field is employed to calculate the reaction path as a free energy surface of calcium dissolution at 298 K in water from the different crystal facets of portlandite. It is also explained why the reactivity of the (010), (100), and (1 1 ¯ 0) crystal facet is higher compared to the (001) facet. In addition, the influence of neighboring Ca crystal sites arrangements on the atomistic dissolution rates is explained as necessary scenarios for the upscaling. The calculated rate constants of all atomistic reaction scenarios provided an input catalog ready to be used in an upscaling kinetic Monte Carlo (KMC) approach.

1. Introduction

Although the crystal dissolution (and precipitation) kinetics are of wide industrial importance, upscaling of fundamental atomistic modeling approaches still presents a major scientific challenge. The dissolution (and precipitation) of minerals phases is governed by complex physico-chemical processes starting at the atomistic scale, where the individual crystal building units (atoms or molecules) at the solid–liquid interface are dissolved and transported into the bulk solution. At mesoscale, e.g., 0.1–10 µm length scale and microsecond time scales, of interest is to predict the evolution of the crystal morphology, as a function of the transport processes and concentration effects. At the atomistic scale, molecular dynamics (MD) computations provide a powerful way for revealing how the atomistic processes of basic building units of the crystals affect the dissolution kinetics. This paper focuses on the atomistic models for the dissolution of crystals. To close the long south bridging gap between the atomistic and the meso-scale modeling approaches for cementitious materials, far-from-equilibrium conditions were considered. Therefore, the dissolution of portlandite (CH, Ca(OH)2), was selected as a relatively simple crystal representative for the dissolution/precipitation process of other more complex cementitious phases (e.g., cement clinker phases) in general. Portlandite is a key by-product of cement hydration in the amount of approximately 27% in volume, which protects the steel reinforcement from corrosion by maintaining a higher pH value [1].
Portlandite dissolution leads to carbonation that accelerates the corrosion of the reinforcement by dropping the pH value of the pore solution by approximately three units [2,3,4]. During the hydration of cement, portlandite does not only precipitated in the bulk of the hardened cement paste matrix but also fills in a more porous interfacial layer between the steel reinforcement and aggregate [5]. Besides its significance in building materials, portlandite is also used in environmental applications such as (flue) acid gas absorption in wet scrubbers, spray dryer reactors, and for pH control of wastewaters. Portlandite can also be used to reduce the carbon footprint as well as thermochemical energy storage [6,7].
Due to portlandite’s relatively simple crystal structure, we propose it as a benchmarking mineral phase in developing atomistic modeling approaches for the dissolution/precipitation process of other (more complex) cementitious minerals in general. The dissolution of cementitious minerals, especially at an atomistic level, is not fully understood yet due to the lack of experimental and modeling techniques available to reach this resolution and bridge the time and scale gaps. Within the past decade, many atomistic-level computational methods have been developed to describe cementitious materials. Density functional theory (DFT) using quantum mechanics is known to be comparatively accurate for geometry optimization [8], vibrational energy [9], band structure calculation [10], mechanical, electronic, and optical properties of the small molecule [11,12]. However, this method is highly computationally expensive and practically not applicable to a large system. On the other hand, classical force field (FF) theory is unable to explain the bond formation and bond breaking. Though Brenner FF theory is able to describe bond-breaking, the Coulomb interactions are not considered [13]. Besides, the bond-energy bond-order (BEBO) method could not explain more complex reactions. Therefore, atomistic simulation approaches using reactive force fields (ReaxFF) parameterized by quantum mechanical calculation, in combination with metadynamics (metaD), have been proposed as a potential solution to study the chemical reactions pathways with satisfactory accuracy and reasonable computing times.
ReaxFF has been developed to elucidate the interfacial reaction mechanism. It has been implemented successfully in hydrocarbons [14], polymer chemistry, metal oxides (Si/SiO2) [15,16], metal hydrides [17], and many other systems. ReaxFF calculates the molecular dynamics in femtosecond (i.e., 10−15 s) time steps for the correct integration of equations of motion. Although the ReaxFF is much more efficient compared to traditional classical force field approaches [18], sometimes it becomes computationally expensive during simulations of interesting events (e.g., atomistic transition states) that typically occurred in a longer time scale. To solve this timescale issue, metadynamics (metaD) is combined with ReaxFF, enabling a powerful coupling that can accelerate the reaction by adding a biased potential that acts on an appropriate number of degrees of freedom usually known as collective variables (CVs) (i.e., inter-atomic/molecular distance, co-ordination number, angles, etc.). During the MD simulation, the bias potential is applied as a sum of Gaussian, continuously growing and acting directly on the localized coordinates of the system [19].
The combined approach has already shown a great potential to calculate the dissolution mechanism and the reactivity of different facets of Portland cement clinkers at the atomistic scale [20]. The current contribution aims at the development of a multi-scale modeling approach that links the atomistic scale to the (sub-)micro-scale. It will open up new insights into the reaction kinetics of portlandite dissolution, and establish a benchmark approach for cementitious materials in general.
In this research (Part 1), a multistep modeling approach has been taken to achieve a clear overview of the dissolution and reactivity of different facets of Portlandite at room temperature (298 K). First, pre-simulations were run to allow for pre-hydration of (001), (100), (010) facets of portlandite, for 600 picoseconds. Then, all pre-hydrated facets were used as an initial geometry to study one-by-one dissolution mechanisms of calcium atoms by using a combination of ReaxFF and metaD.
The reactivity of different crystal facets was compared with the dissolution profile (free energy surface) of calcium. The orientation of calcium and the number of (crystal site) neighbors are different on each crystal facet and change with dissolution progress due to removal of the neighbors. This creates many possibilities which can take place and are called here dissolution scenarios. Therefore, the influence of crystal site neighbors is considered here as a critical step in upscaling the individual atomistic simulations into a meso-scale simulation. Finally, the results are used to calculate the atomistic rate constants using a transition state theory (TST). For this, the most critical possible events have to be identified and quantified to be provided as an input for upscaling using a kinetic Monte Carlo (KMC) approach (part 2 paper) and calculating the overall rate of the dissolution and morphological changes [21].

2. Methods and Modeling Approach

The ReaxFF has been employed successfully in the cementitious system by merging two parameter sets (Si-O-H and Ca-O-H) developed by Fogarty et al. [22] and Manzano et al. [16], individually. This parameter set has already explained the absorption of water to C3S clinker phases [20]. The computer simulations were carried out by using reactive force field theory (ReaxFF) in LAAMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) platform [23].
Usually, in a force field MD simulation, the femtoseconds (10−15 s) time step is considered. However, the majority of the interesting phenomena (i.e., transition state events (TS)) occur at a much larger time scale. Therefore, the rare event calculation requires millions of time steps, which increases the computational cost. Therefore, many sampling methods have been developed to reduce the computation time.
Metadynamics is a powerful algorithm applied for the free energy calculation by accelerating the rare event sampling. In metaD, a history dependent external biased potential as a function of CVs is imposed in a controlled manner to enhance the exploration of free energy surface (FES). The potential is added along with the CVs space as a sum of Gaussians, which encourages the system to visit the unexplored configurations. After prolonging the simulation for a sufficiently long time, the biased potential is converged to reconstruct the FES as a function of that selected CVs [24].

Model Construction

The Portlandite is a trigonal layered crystal structure with space group P 3 ¯ m . Each hydroxyl group connects three Ca atoms in its (001) layer and is also surrounded by three other hydroxyl groups from the next (001) layer. In this work, a hexagonal Portlandite crystal model was used with the lattice parameter of a = b = 3.586 Å, c = 4.911 Å, α = β = 90°, and γ = 120° [25]. The fresh cleaved (001) Portlandite orthogonal periodic simulation cell (17.80 × 19.64 × 38.65) × 10−30 m3 contained 1302 atoms was built using virtual nano lab (VNL) [26] and Avogadro [26,27]. The geometry of the simulation cell was optimized using Hessian-free truncated Newton algorithm (HFTN) where the stopping tolerances for energy were 1.0 × 10−4 and force was 1.0 × 10−6 kcal mol−1 Å−1 [28]. Afterward, 6.99 × 10−27 m3 periodic cell filled with randomly distributed water molecules (density 1000 kg m−3) was added up to the optimized (001) Portlandite facet using packmole [29]. First, the simulation cell was equilibrated for 150 picoseconds (0.5 femtoseconds time steps) at 298 K and 1 atm using Nose−Hoover thermostat (NVT, canonical ensemble. Later, the equilibrated geometry was hydrated for 600 picoseconds using Nose−Hoover barostat (in NPT canonical ensemble). A periodic boundary condition was applied during the entire simulation.
The hydrated (001) facet of Portlandite (after 600 picoseconds) was taken as an input geometry to compute the dissolution mechanism of calcium using the combined metaD and ReaxFF.
The PLUMED [30,31] package as an extension of LAMMPS [23] was employed for metaD simulation. It speeds up the simulation using history-dependent biased potential. A calcium Ca-588 from (001) facet located in the center of hexagonally oriented neighboring Ca (Figure 1a) of the hydrated Portlandite facet is selected. A well-tempered metaD scheme was applied to remove calcium from the crystal facet into the pore solution. The distance between the center of mass (COM) of the crystal and Ca-588 is selected as a CV. It is worth noting that the choice of the correct CVs is crucial for getting a good physical description. The bias potential was added as a Gaussian with a height of 6.28 kJ/mol and full width at half-maximum of 0.2 × 10−10 m in every 0.02 picoseconds.
The simulation was performed for 300 picoseconds (until convergence) at 298 K with NPT ensemble and the free energy of dissolution was computed over the entire. In order to understand the influence of hexagonally oriented Ca neighbors, the identical approach was applied to compute the dissolution profile of central Ca (red) in absence of 1,2,3,4,5,6 neighbor individually being pre-deleted in a clockwise manner (Figure 1b–g). In this manner, the activation barrier for a total of seven scenarios was calculated to investigate the influence of neighbors (scenario) configurations that are expected to occur during the (longer-term) dissolution. Those scenarios are of key relevance for upscaling, e.g., to provide required inputs for kinetic Monte Carlo simulations.
A similar approach was used for other crystal facets: (100), (010), and (1 1 ¯ 0) cleavage facets of portlandite (Table A1), where after the pre-hydration simulations, the possible (most relevant) scenarios for the dissolution mechanism are simulated. The reactivities were calculated from the dissolution profile (FES) of the calcium of the individual facet. Differently to the (001) facet, where the neighbor Ca is arranged in a hexagonal pattern, the Ca in on (100), (1 1 ¯ 0), and (010) facets are arranged in a row (linear, i.e., left and right) manner). Therefore, a total of three scenarios were considered to calculate the dissolution behavior of the central Ca before and after removing 1 and 2 neighboring Ca located on both sides (marked green) in the same row (Figure 2). In addition, there is one neighbor below (direction-z) which cannot be dissolved before the central atom.
Finally, this study delivers a complete set of all (most important) events for the four different facets of portlandite, relevant for a hexagonal crystal morphology (Figure 1). Results of these simulations provide energy barriers, which are exponentially proportional to the individual rate constant of all individual events, e.g., according to transition state theory (TST). A full catalog of individual (atomistic) reaction rates provided basic input of upscaling KMC simulation for calculating the overall rate of the dissolution and morphological changes (published separately as part 2 paper to this Issue).

3. Results and Discussion

3.1. Pre-Hydration of Portlandite

In order to understand the interaction between bulk water and Portlandite facet, the dynamics of the system have been followed for 600 picoseconds at 298 K and standard ambient pressure (101.325 kPa). The periodic boundary condition was applied to all three directions to improve the representativeness of the elementary volume cube size during simulation. Hence, when an atom passes one side of the cell, it reappears on the opposite side with the same velocity. A periodic dimension can change in size due to constant pressure boundary conditions. In addition, 0.5 fs (femtosecond) time steps were used for the entire simulation to accurately capture the fastest vibrations of hydrogen, as the smallest and lightest element in the system. Moreover, ReaxFF can be considered the smallest atoms as bonded only for sufficiently small time steps.
According to our observation, at the very initial stage, the water molecules interact with the facets of Portlandite and the charge transfer between the O atom (part of the hydroxyl group) in the first layer of (010) facet and the H ion of the water molecule is comparatively more obvious than the charge transfer between Ca and the O (water) atoms. Therefore, the electrophilic attack is much stronger than the nucleophilic attack during water adsorption on the (010) facet. Afterward, the proton transfer occurs from the hydroxyl of the surface to the inner oxygen by leaving the first oxygen-free for further reaction (hopping process) as described by Manzano et al. [32]. Pre-hydration is carried out only for 600 picoseconds.
The (010) facet (prismatic facet) of Portlandite (Figure 3d) shows higher reactivity during hydration compared to the basal facets (001). The main reason is that the atoms on the (010) facet can be the electrophilic and nucleophilic center, therefore, the charges were balanced more significantly, resulting in lower surface energy. This indicates the reason for the higher interaction between bulk water and the facet, resulting in the formation of voids (Figure 3d). Similarly, (100) and ( 1 1 ¯ 0 ) facets also exhibit such behavior at an intermediate level, indicating intermediate reactivity (Figure 3b,c).
In contrast, the (001) crystal facet indicates the lowest reactivity due to their highest surface energy and water tessellation during pre-hydration. The term ‘water tessellation’ describes the arrangement of the water molecule as a quasi-static 2D sheet barrier [33,34]. The water tessellation at the (001) facet prevents water from penetration into the crystal and the dissolution of calcium ions from the facet (Figure 3a).
Traditionally, the reactivity of surfaces was attempted to correlate with their static properties, i.e., surface energies. However, the results of dynamic pre-hydration by MD simulations at the interface between the Portlandite and water, demonstrate the lack of such a correlation. Generally, the facet properties are modified during the interaction of water. Therefore, only a correlation of the facet properties is insufficient. Hence, the reactivity of portlandite facets depends on the interfacial reaction. The free energy calculation during the dissolution of Ca from the different facets of portlandite to the solution would provide an appropriate explanation about the reactivity.

3.2. Dissolution of Calcium from (001) Facet of Portlandite

Free energy surface (FES) calculation has received great interest in MD simulations for obtaining a clear impression of reaction paths, including transition states. Well-tempered metaD is an efficient method that accelerates the chemical reaction to overcome the activation barrier by. Selecting suitable CVs is the main challenge in metaD-based computation. Calculating the FES of a targeted region is the main advantage of metaD. In the current study, the ReaxFF coupled with metaD enabled the calculation of the dissolution pathways of calcium from the facets of portlandite.
Figure 4a represents the dissolution profile of calcium from (001) crystalline facet of portlandite at 298 K selecting a single CVs distance between the central Ca-588 surrounded by six Ca neighbors (hexagonally oriented) and the center of mass of the crystal. The x-axis represents the distance (reaction coordinate) in Å (10−10 m) starting from the initial state on the facet (0 Å), to the solution (final state) at the next lowest minima.
The free energy surface, from the metaD calculation, represents the movement of Ca-588 from the center of the hexagon to the pore solution required to overcome the huge activation barrier of 352.0 kJ/mol at 5.88 × 10−10 m. (Figure 4a). The total free energy change (ΔG) of +280.8 kJ/mol and high activation barrier are directed to the endergonic, i.e., thermodynamically not favorable process and less reactive facets.
Nevertheless, after removing the first and second neighbor clockwise (Figure 1b,c) the activation barrier decreases to 175.4 kJ/mol (still high); however, the dissolution process becomes thermodynamically favorable (ΔG = −23.6 kJ/mol) (Figure 4b,c). Further removal of neighboring Ca one by one reduces the activation barrier significantly for the dissolution of the same central Ca and reaches its minimum value of 25.9 kJ/mol (Table 1) when all six neighbors are missing (Figure 4d–g). The results complement the previous observation, where after 600 ps of pre-hydration, the (001) facet of Portlandite was found less reactive due to stabilizing the facet by its compact geometry and water tessellation [32]. Therefore, it is very difficult to dissolute Ca from the perfect (001) crystal facet; however, the defects, grain boundaries, and missing Ca neighbors (as demonstrated in Figure 4) could enhance its reactivity.
Furthermore, a correlation is done between the seven scenarios of Figure 1 and the resulting Free energy of activation (Table 1). The values of the activation barrier can be fitted as a function of the number of neighbors (xi) with an allometric (basic Origin) function yi = a + bxic, where a = 46.03, b = −1.43, and c = 2.98 as calibrated parameters. The R-square coefficient of 0.958 was indicated that the model follows the trend of the data satisfactorily (Figure 4h).

3.3. Dissolution of Calcium from Prismatic (100), (110), and (010) Facets of Portlandite

The free energy profile of calcium dissolution from the (010) facet shows higher reactivity than (001) in agreement with the results after 600 picoseconds of pre-hydration. The dissolution of Ca-552 from (010) facet was required to overcome the barrier of 29.9 kJ/mol at 1.71 × 10−10 m (Table 1). The total free energy change (ΔG) of −299.7 kJ/mol at 298 K indicates an exergonic reaction and thermodynamically favorable [35,36]. This explains the high reactivity and higher interaction of water, so rapidly attracted to the facet that it results in the creation of voids (Figure 3d). Moreover, the activation barrier for the dissolution of the same central calcium (Ca-552) decreased to 20.5 kJ/mol and 7.1 kJ/mol in absence of the first and second neighbor respectively which was removed one by one manner (Figure 5h,i).
Likewise, Ca dissolution (central Ca) from ( 1 1 ¯ 0 ) facet of portlandite followed the same (thermodynamically favorable) behavior as (010) facet—however, in a comparatively less reactive fashion, due to the slightly higher activation barrier (56.0 kJ/mol) in the presence of both neighbors, while lower for one neighbor missing (34.9 kJ/mol) (Figure 5d,e). Nevertheless, in absence of both neighbors, the activation barrier was found by far the lowest (0.7 kJ/mol) among all the studied scenarios during the dissolution of central Ca (Table 1).
Among the three prismatic facets of portlandite, the dissolution behavior of the (100) facet resembled the basal facet (001) in terms of facet reactivity and influence of Ca neighbors. The complete dissolution of central calcium (Ca-542) in presence of all neighbors requires 195.3 kJ/mol at 2.22 × 10−10 m which is lower than the (001) facet in a similar scenario. The total free energy change (ΔG) of +111.1 kJ/mol at 298 K represents an endergonic and thermodynamically unfavorable process (Figure 5a). The activation energy reduces to 114.6 kJ/mol when one neighbor is missing, still following a similar trend. Furthermore, in absence of both neighbors, the dissolution process become favorable (−62.3 kJ/mol) and the free energy of activation was further reduced to 70.0 kJ/mol which was higher compared to the (010) and ( 1 1 ¯ 0 ) facets even in presence of both neighbors (Figure 5b,c).
The correlation diagram represents the increase of the activation barrier with increasing the number of neighbors for all three prismatic facets (Figure 6).

3.4. Reactivity of Different Facets of Portlandite

Considering the free energy calculation for the different facets of Portlandite in different scenarios, both (001) and (100) facets were not reactive especially in the perfect condition (no missing neighbors); however, they become reactive in absence of neighbors. In contrast, (010) and ( 1 1 ¯ 0 ) were found reactive in all cases where the reactivity increases with the number of missing neighbors. Therefore, neighbors played a significant role in the reactivity of the dissolution mechanism.
Finally, the thermodynamic properties for the dissolution of calcium at the perfect crystalline facets, the reactivity order for the different facets of portlandite are shown in Figure 7.

3.5. Upscaling the Dissolution Rate for Different Facets of Portlandite

To upscale the atomistic simulations towards much larger timescales and meso-scopic—i.e., sub-micrometer system-sizes, kinetic Monte Carlo (KMC) simulations are applied. The principal idea of a KMC simulation is in bridging the timescale gap by coarse-graining the time evolution and focusing on discrete rare events using Markovian state-to-state dynamics [37]. The calculated activation barrier obtained here by MD simulations (using ReaxFF coupled with MetaD) is tabulated for all scenarios for the dissolution of calcium from the different crystalline facets (Table 1). This table is considered to be the most critical and tedious step for the KMC upscaling model. Results of this paper are provided as input data for KMC simulations in Part 2 publication [38], where the rate of calcium dissolution for the individual scenarios was calculated using transition state theory [39]. Finally, KMC can predict the overall mesoscopic (far-from equilibrium) rate of dissolution and morphological change for any user-defined initial crystal size and morphology.

4. Conclusions

The primary objective of this research was to achieve a deeper understanding of the dissolution mechanism of portlandite; therefore, different crystal faces were considered. The upscaling approach was established for the atomistic scale simulation using ReaxFF coupled with the metaD to feed the sub-micro KMC modeling approach. Simulations identified the key dissolution reaction events from possible scenarios depending on nearest neighbors Ca configurations, and the different crystal facets of Portlandite.
The results show that the (010) and ( 1 1 ¯ 0 ) facets (prismatic facet) were found most reactive and the dissolution of calcium was thermodynamically favorable. Besides, the lower the number of (dissolved) Ca-neighbors, the lower the activation barrier—i.e., higher atomistic event dissolution rate. In contrast, (001) and (100) were shown the lowest reactivity, and the dissolution was unfavorable at the perfect facet. However, the dissolution rate increased by decreasing the number of (dissolved) Ca-neighbors, which explains that the crystal defects, grain boundaries, and absent neighbors increase the dissolution rate.
The (001) facet was found to be the least reactive facet, where the water tessellation maintained the atomic arrangement of the crystal facets, preventing Ca dissolution. Whereas high interaction of water and lower activation barrier on the dissolution profile of Ca indicated the (010) facet to be the most reactive one.
Finally, the calculated dissolution rate catalog of the most important scenarios is obtained to feed KMC simulations for upscaling computations to predict the mesoscopic dissolution rates and evolution of the crystal morphology.

Author Contributions

Conceptualization, K.M.S.U., N.U. (MD scenarios for KMC upscaling); Methodology, investigation, visualization, K.M.S.U.; Validation, K.M.S.U. and N.U.; Writing—original draft preparation, K.M.S.U.; Writing—review and editing, N.U., B.M., M.I. and E.K.; Supervision, N.U., B.M. and E.K.; Funding acquisition, N.U., E.K. and B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Research Foundation (DFG) organization under project number 455605608 titled: ‘Elucidation of initial cement dissolution mechanism by the gap-bridging multiscale modeling approach (CEM-bridge)’.

Acknowledgments

The authors would like to express their gratitude to Nikolas Luke, M. Sc. for his technical support in running MD simulation successfully in the high-performance Linux cluster at the University of Kassel. The authors also acknowledge support from the German Research Foundation (DFG) and the Open Access Publishing Fund of Technische Universität Darmstadt.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MDMolecular Dynamics
CHPortlandite (Ca(OH)2)
KMCKinetic Monte Carlo
ReaxFFReactive Force Field
LAMMPS Large-Scale Atomic/Molecular Massively Parallel Simulator
metaDMetadynamics
CVsCollective Variables
TSTransition State
VNLVirtual Nano Lab
FESFree Energy Surfaces (kJ/mol)
HFTNHessian-Free Truncated Newton Algorithm
NVTNose-Hoover Thermostat
NPTNose-Hoover Pressure Barostat
Δ G * Free Energy of Activation Barrier (kJ/mol)
Δ G Free Energy Change (kJ/mol)

Appendix A

Table A1. Crystallographic data of Portlandite consisting of different crystal facets at 298 K.
Table A1. Crystallographic data of Portlandite consisting of different crystal facets at 298 K.
Crystal Facet of PortlanditeThe Dimension of the Sumulation Cell
3/10−30 m3)
No. of Atoms in the Simulation Cell
(001)17.67, 19.94, 40.751307
(100)20.75, 19.94, 38.421431
(1 1 ¯ 0)20.75, 26.64, 46.412186
(010)17.67, 20.75, 34.941155

References

  1. Beaudoin, J.J.; Catinaud, S.; Marchand, J. Volume stability of calcium hydroxide in aggressive solutions. Cem. Concr. Res. 2001, 31, 149–151. [Google Scholar] [CrossRef]
  2. Lea, F.M. Lea’s Chemistry of Cement and Concrete, 4th ed.; Elsevier Butterworth-Heinemann: Oxford, UK, 2004; ISBN 0750662565. [Google Scholar]
  3. Taylor, H.F.W. Cement Chemistry; Thomas Telford: London, UK, 1997; pp. 271–272. [Google Scholar]
  4. Glass, G.K.; Reddy, B.; Buenfeld, N.R. Corrosion inhibition in concrete arising from its acid neutralisation capacity. Corros. Sci. 2000, 42, 1587–1598. [Google Scholar] [CrossRef]
  5. Moorehead, D.R. Cementation by the carbonation of hydrated lime. Cem. Concr. Res. 1986, 16, 700–708. [Google Scholar] [CrossRef]
  6. Couling, D.J.; Das, U.; Green, W.H. Analysis of Hydroxide Sorbents for CO2 Capture from Warm Syngas. Ind. Eng. Chem. Res. 2012, 51, 13473–13481. [Google Scholar] [CrossRef]
  7. Nikulshina, V.; Gálvez, M.E.; Steinfeld, A. Kinetic analysis of the carbonation reactions for the capture of CO2 from air via the Ca(OH)2–CaCO3–CaO solar thermochemical cycle. Chem. Eng. J. 2007, 129, 75–83. [Google Scholar] [CrossRef]
  8. Izadifar, M.; Dolado, J.S.; Thissen, P.; Ayuela, A. Interactions between Reduced Graphene Oxide with Monomers of (Calcium) Silicate Hydrates: A First-Principles Study. Nanomaterials 2021, 11, 2248. [Google Scholar] [CrossRef] [PubMed]
  9. Izadifar, M.; Königer, F.; Gerdes, A.; Wöll, C.; Thissen, P. Correlation between Composition and Mechanical Properties of Calcium Silicate Hydrates Identified by Infrared Spectroscopy and Density Functional Theory. J. Phys. Chem. C 2019, 123, 10868–10873. [Google Scholar] [CrossRef]
  10. Ren, K.; Zheng, R.; Yu, J.; Sun, Q.; Li, J. Band Bending Mechanism in CdO/Arsenene Heterostructure: A Potential Direct Z-Scheme Photocatalyst. Front. Chem. 2021, 9, 788813. [Google Scholar] [CrossRef] [PubMed]
  11. Ren, K.; Zheng, R.; Xu, P.; Cheng, D.; Huo, W.; Yu, J.; Zhang, Z.; Sun, Q. Electronic and Optical Properties of Atomic-Scale Heterostructure Based on MXene and MN (M = Al, Ga): A DFT Investigation. Nanomaterials 2021, 11, 2236. [Google Scholar] [CrossRef] [PubMed]
  12. Ren, K.; Shu, H.; Huo, W.; Cui, Z.; Yu, J.; Xu, Y. Mechanical, electronic and optical properties of a novel B2P6 monolayer: Ultrahigh carrier mobility and strong optical absorption. Phys. Chem. Chem. Phys. 2021, 23, 24915–24921. [Google Scholar] [CrossRef] [PubMed]
  13. Brenner, D.W. Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films. Phys. Rev. B 1990, 42, 9458–9471. [Google Scholar] [CrossRef]
  14. Chenoweth, K.; van Duin, A.C.T.; Goddard, W.A. ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. J. Phys. Chem. A 2008, 112, 1040–1053. [Google Scholar] [CrossRef] [Green Version]
  15. Van Duin, A.C.T.; Strachan, A.; Stewman, S.; Zhang, Q.; Xu, X.; Goddard, W.A. ReaxFF SiO Reactive Force Field for Silicon and Silicon Oxide Systems. J. Phys. Chem. A 2003, 107, 3803–3811. [Google Scholar] [CrossRef] [Green Version]
  16. Manzano, H.; Pellenq, R.J.; Ulm, F.J.; Buehler, M.J.; van Duin, A.C. Hydration of calcium oxide surface predicted by reactive force field molecular dynamics. Langmuir 2012, 28, 4187–4197. [Google Scholar] [CrossRef]
  17. Cheung, S.; Deng, W.Q.; van Duin, A.C.; Goddard, W.A. ReaxFF(MgH) reactive force field for magnesium hydride systems. J. Phys. Chem. A 2005, 109, 851–859. [Google Scholar] [CrossRef] [Green Version]
  18. Senftle, T.P.; Hong, S.; Islam, M.M.; Kylasa, S.B.; Zheng, Y.; Shin, Y.K.; Junkermeier, C.; Engel-Herbert, R.; Janik, M.J.; Aktulga, H.M.; et al. The ReaxFF reactive force-field: Development, applications and future directions. NPJ Comput. Mater. 2016, 2, 9396. [Google Scholar] [CrossRef]
  19. Barducci, A.; Bonomi, M.; Parrinello, M. Metadynamics. WIREs Comput. Mol. Sci. 2011, 1, 826–843. [Google Scholar] [CrossRef]
  20. Salah Uddin, K.M.; Middendorf, B. Reactivity of Different Crystalline Surfaces of C3S During Early Hydration by the Atomistic Approach. Materials 2019, 12, 1514. [Google Scholar] [CrossRef] [Green Version]
  21. Martin, P.; Gaitero, J.J.; Dolado, J.S.; Manzano, H. New Kinetic Monte Carlo Model to Study the Dissolution of Quartz. ACS Earth Space Chem. 2021, 5, 516–524. [Google Scholar] [CrossRef]
  22. Fogarty, J.C.; Aktulga, H.M.; Grama, A.Y.; van Duin, A.C.; Pandit, S.A. A reactive molecular dynamics simulation of the silica-water interface. J. Chem. Phys. 2010, 132, 174704. [Google Scholar] [CrossRef]
  23. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  24. Laio, A.; Gervasio, F.L. Metadynamics: A method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science. J. Phys. Chem. A 2008, 71, 126601. [Google Scholar] [CrossRef]
  25. Busing, W.R.; Levy, H.A. Neutron Diffraction Study of Calcium Hydroxide. J. Chem. Phys. 1957, 26, 563–568. [Google Scholar] [CrossRef]
  26. Schneider, J.; Hamaekers, J.; Chill, S.T.; Smidstrup, S.; Bulin, J.; Thesen, R.; Blom, A.; Stokbro, K. ATK-ForceField: A new generation molecular dynamics software package. Model. Simul. Mater. Sci. Eng. 2017, 25, 85007. [Google Scholar] [CrossRef] [Green Version]
  27. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeersch, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [Green Version]
  28. Tuckerman, M.E.; Alejandre, J.; López-Rendón, R.; Jochim, A.L.; Martyna, G.J. A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal–isobaric ensemble. J. Phys. A Math. Gen. 2006, 39, 5629–5651. [Google Scholar] [CrossRef]
  29. Martínez, L.; Andrade, R.; Birgin, E.G.; Martínez, J.M. PACKMOL: A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157–2164. [Google Scholar] [CrossRef]
  30. Bonomi, M.; Branduardi, D.; Bussi, G.; Camilloni, C.; Provasi, D.; Raiteri, P.; Donadio, D.; Marinelli, F.; Pietrucci, F.; Broglia, R.A.; et al. PLUMED: A portable plugin for free-energy calculations with molecular dynamics. Comput. Phys. Commun. 2009, 180, 1961–1972. [Google Scholar] [CrossRef] [Green Version]
  31. Tribello, G.A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED 2: New feathers for an old bird. Comput. Phys. Commun. 2014, 185, 604–613. [Google Scholar] [CrossRef] [Green Version]
  32. Manzano, H.; Durgun, E.; López-Arbeloa, I.; Grossman, J.C. Insight on Tricalcium Silicate Hydration and Dissolution Mechanism from Molecular Simulations. ACS Appl. Mater. Interfaces 2015, 7, 14726–14733. [Google Scholar] [CrossRef]
  33. Miranda, P.B.; Xu, L.; Shen, Y.R.; Salmeron, M. Icelike Water Monolayer Adsorbed on Mica at Room Temperature. Phys. Rev. Lett. 1998, 81, 5876–5879. [Google Scholar] [CrossRef]
  34. Ostroverkhov, V.; Waychunas, G.A.; Shen, Y.R. New information on water interfacial structure revealed by phase-sensitive surface spectroscopy. Phys. Rev. Lett. 2005, 94, 46102. [Google Scholar] [CrossRef] [Green Version]
  35. Salah Uddin, K.M.; Middendorf, B. Atomistic modeling of early Hydration of C3S. In Computational Modelling of Concrete and Concrete Structures, 1st ed.; CRC Press/Balkema: Boca Raton, FL, USA, 2018; pp. 103–107. [Google Scholar]
  36. Salah Uddin, K.M.; Middendorf, B. Insight into the reactivity of different crystalline phases of C3S during early hydration by the atomistic approach. In Proceedings of the 20th International Conference on Building Materials, Bauhaus-Universität, Wimar, Germany, 12–14 September 2018; Volume 1, pp. 731–736. [Google Scholar]
  37. Chatterjee, A.; Vlachos, D.G. An overview of spatial microscopic and accelerated kinetic Monte Carlo methods. J. Comput.-Aided Mater. Des. 2007, 14, 253–308. [Google Scholar] [CrossRef]
  38. Izadifar, M.; Ukrainczyk, N.; Salah Uddin, K.M.; Middendorf, B.; Koenders, E. Dissolution of Portlandite in Pure Water: Part 2 Atomistic Kinetic Monte Carlo (KMC) Approach. Materials 2022, 15, 1442. [Google Scholar] [CrossRef]
  39. Keith, J.L.; King, M.C. The development of transition-state theory. J. Phys. Chem. A 1983, 87, 2657–2664. [Google Scholar]
Figure 1. Dissolution scenarios of hexagonally oriented central Ca (red) from (001) facet of Portlandite (a) in the presence and (bg) absence of 1,2,3,4,5,6 neighbors (green) removed one by one (depicted in a clockwise manner).
Figure 1. Dissolution scenarios of hexagonally oriented central Ca (red) from (001) facet of Portlandite (a) in the presence and (bg) absence of 1,2,3,4,5,6 neighbors (green) removed one by one (depicted in a clockwise manner).
Materials 15 01404 g001
Figure 2. Dissolution scenarios of central Ca (violet) in absence of first and second nearest neighbors (green) removed one by one from: (a) (100), (b) (1 1 ¯ 0) and (c) (010) facet of Portlandite.
Figure 2. Dissolution scenarios of central Ca (violet) in absence of first and second nearest neighbors (green) removed one by one from: (a) (100), (b) (1 1 ¯ 0) and (c) (010) facet of Portlandite.
Materials 15 01404 g002
Figure 3. Comparison of the reactivity between the different facets of Portlandite: (a) (001), (b) (100), (c) ( 1 1 ¯ 0 ), and (d) (010). Snapshots after pre-hydration simulations for 600 picoseconds at 298 K.
Figure 3. Comparison of the reactivity between the different facets of Portlandite: (a) (001), (b) (100), (c) ( 1 1 ¯ 0 ), and (d) (010). Snapshots after pre-hydration simulations for 600 picoseconds at 298 K.
Materials 15 01404 g003
Figure 4. Representative dissolution profile (free energy surface) of the hexagonally oriented central Ca (red Ca in Figure 1) from (001) crystalline facet of Portlandite in different scenarios: the presence of all six Ca neighbors (a) and absence of the different number of 1,2,3,4,5,6 (bg) hexagonally oriented neighbors. (h) The correlation diagram of free energy of activation of central Ca of (001) facet and hexagonally oriented number of Ca neighbors at 298 K.
Figure 4. Representative dissolution profile (free energy surface) of the hexagonally oriented central Ca (red Ca in Figure 1) from (001) crystalline facet of Portlandite in different scenarios: the presence of all six Ca neighbors (a) and absence of the different number of 1,2,3,4,5,6 (bg) hexagonally oriented neighbors. (h) The correlation diagram of free energy of activation of central Ca of (001) facet and hexagonally oriented number of Ca neighbors at 298 K.
Materials 15 01404 g004
Figure 5. Characteristic dissolution profile (free energy surface) of the central Ca (violet Ca in the Figure 2) (ac), from (100); (df); from ( 1 1 ¯ 0 ); and (gi); from (010); crystalline facets of Portlandite in the different scenarios: before and after removal of one and two neighbors.
Figure 5. Characteristic dissolution profile (free energy surface) of the central Ca (violet Ca in the Figure 2) (ac), from (100); (df); from ( 1 1 ¯ 0 ); and (gi); from (010); crystalline facets of Portlandite in the different scenarios: before and after removal of one and two neighbors.
Materials 15 01404 g005
Figure 6. Correlation diagram of free energy of activation of central Ca (violet Ca in Figure 2) of (100), ( 1 1 ¯ 0 ), and (010) crystalline facets of Portlandite as a function of number of Ca neighbors at 298 K.
Figure 6. Correlation diagram of free energy of activation of central Ca (violet Ca in Figure 2) of (100), ( 1 1 ¯ 0 ), and (010) crystalline facets of Portlandite as a function of number of Ca neighbors at 298 K.
Materials 15 01404 g006
Figure 7. Overall reactivity order of different facets of Portlandite comparing dissolution profile of central calcium in presence of all neighbors (at 298 K).
Figure 7. Overall reactivity order of different facets of Portlandite comparing dissolution profile of central calcium in presence of all neighbors (at 298 K).
Materials 15 01404 g007
Table 1. Free energy change of different facets of Portlandite during the dissolution of calcium at standard ambient temperature (298 K) and pressure (101.325 kPa).
Table 1. Free energy change of different facets of Portlandite during the dissolution of calcium at standard ambient temperature (298 K) and pressure (101.325 kPa).
Portlandite Crystal PlanePossible Scenarios of Dissolution of Central CaFree Energy of Activation ( Δ G * )
kJ/mol
Free Energy Change ( Δ G )
kJ/mol
Thermodynamic Properties
(001)In presence of 6 neighborsΔ352.0+280.8endergonic
In absence of 1 neighbor199.1+117.7endergonic
In absence of 2 neighbors175.4−23.6exergonic
In absence of 3 neighbors56.1−180.5exergonic
In absence of 4 neighbors55.8−44.3Exergonic
In absence of 5 neighbors54.9−130.0Exergonic
In absence of 6 neighbors25.9−147.6Exergonic
(100)In presence of 2 neighbors195.3+111.1Endergonic
In absence of 1 neighbor114.6+5.6Endergonic
In absence of 2 neighbors70.0−62.3Exergonic
( 1 1 ¯ 0 )In presence of 2 neighbors56.0−92.2Exergonic
In absence of 1 neighbor34.9−59.4Exergonic
In absence of 2 neighbors0.7−211.9Exergonic
(010)In presence of 2 neighbors29.9−299.7Exergonic
In absence of 1 neighbor20.5−237.9Exergonic
In absence of 2 neighbors7.1−126.0Exergonic
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salah Uddin, K.M.; Izadifar, M.; Ukrainczyk, N.; Koenders, E.; Middendorf, B. Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach. Materials 2022, 15, 1404. https://doi.org/10.3390/ma15041404

AMA Style

Salah Uddin KM, Izadifar M, Ukrainczyk N, Koenders E, Middendorf B. Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach. Materials. 2022; 15(4):1404. https://doi.org/10.3390/ma15041404

Chicago/Turabian Style

Salah Uddin, Khondakar Mohammad, Mohammadreza Izadifar, Neven Ukrainczyk, Eduardus Koenders, and Bernhard Middendorf. 2022. "Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach" Materials 15, no. 4: 1404. https://doi.org/10.3390/ma15041404

APA Style

Salah Uddin, K. M., Izadifar, M., Ukrainczyk, N., Koenders, E., & Middendorf, B. (2022). Dissolution of Portlandite in Pure Water: Part 1 Molecular Dynamics (MD) Approach. Materials, 15(4), 1404. https://doi.org/10.3390/ma15041404

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop