Monte Carlo Dynamics Simulations Research Articles

Synthesis, structural characterization, thermal analysis and anticorrosion properties of novel AZO-based imidazo[1,2-a]pyridine derivatives: Experimental study and theoretical approaches

In this study, innovative AZO-based imidazo[1,2-a]pyridine derivatives were synthesized using simple and efficient synthetic approach enabling large-scale production. Thermal analysis highlights their suitability for potential applications requiring resistance to moderate temperatures. These compounds were investigated as corrosion inhibitors for MS in a 1 M HCl-medium. Corrosion inhibition studies, including weight loss, PotentioDynamic polarization and electrochemical impedance spectroscopy methods, revealed that AZOs exhibited maximum corrosion inhibition efficiencies reaching 96.13 % and 94.07 %, respectively, at a concentration of 10−3 M and a temperature of 298 K, with both compounds acting as mixed-type inhibitors. The study also delved into the influence of temperature and immersion time on inhibitory performance, revealing stable inhibition within the temperature range of 298–328 K. The calculated standard free energy of adsorption was −42.06 and −40.62 kJ/mol for AZO 1 and AZO 2, respectively, suggesting mixed-type adsorption, with adsorption behavior following the Langmuir isotherm. Morphological analysis through SEM and AFM, showed a reduction in surface roughness in the presence of the inhibitors, while water contact angle measurements indicated a significant increase in surface hydrophobicity. UV–visible and FT-IR analyses confirmed the interaction of AZOs with the MS surface, suggesting the formation of a protective layer composed of adsorbed AZO molecules. Furthermore, computational investigations employing DFT, Monte Carlo, and Molecular Dynamic simulation revealed the formation of coordinate bonds between the studied inhibitors and the metal surface, confirming adsorption on the metallic surface. These findings align well with experimental outcomes, and a possible corrosion inhibition mechanism process was discussed.

Dissipative Self-Assembly of Patchy Particles under Nonequilibrium Drive: A Computational Study.

Inspired by biology and implemented using nanotechnology, the self-assembly of patchy particles has emerged as a pivotal mechanism for constructing complex structures that mimic natural systems with diverse functionalities. Here, we explore the dissipative self-assembly of patchy particles under nonequilibrium conditions, with the aim of overcoming the constraints imposed by equilibrium assembly. Utilizing extensive Monte Carlo (MC) and Molecular Dynamics (MD) simulations, we provide insight into the effects of external forces that mirror natural and chemical processes on the assembly rates and the stability of the resulting assemblies comprising 8, 10, and 13 patchy particles. Implemented by a favorable bond-promoting drive in MC or a pulsed square wave potential in MD, our simulations reveal the role these external drives play in accelerating assembly kinetics and enhancing structural stability, evidenced by a decrease in the time to first assembly and an increase in the duration the system remains in an assembled state. Through the analysis of an order parameter, entropy production, bond dynamics, and interparticle forces, we unravel the underlying mechanisms driving these advancements. We also validated our key findings by simulating a larger system of 100 patchy particles. Our comprehensive results not only shed light on the impact of external stimuli on self-assembly processes but also open a promising pathway for expanding the application by leveraging patchy particles for novel nanostructures.

Stability and dynamics of magnetic skyrmions in FM/AFM heterostructures

Magnetic skyrmions have garnered attention for their potential roles in spintronic applications, such as information carriers in computation, data storage, and nano-oscillators due to their small size, topological stability, and the requirement of small electric currents to manipulate them. Two key challenges in harnessing skyrmions are the stabilization requirement through a strong out-of-plane field, and the skyrmion Hall effect (SkHE). Here, we present a systematic model study of skyrmions in ferromagnetic/antiferromagnetic (FM/AFM) multilayer structures by employing both atomistic Monte Carlo and atomistic spin dynamics simulations. We demonstrate that skyrmions stabilized by exchange bias have superior stability to field-stabilized skyrmions due to the formation of a magnetic imprint within the AFM layer. Additionally, stacking two skyrmion hosting FM layers between two AFM layers suppresses the SkHE and enables the transport of AFM-coupled skyrmions with high velocity in the order of a few km/s. This proposed multilayer configuration could serve as a pathway to overcome existing limitations in the development of skyrmion-based devices, and the insights obtained through this study contribute significantly to the broader understanding of topological spin textures in magnetic materials. Published by the American Physical Society 2024

Measuring Oxygen Solubility in Amorphous and Semicrystalline Polyolefins Using Test Particle Insertion: A Comparative Study of Polyethylene and Isotactic Polypropylene.

The test particle insertion method is used to study the solubility of oxygen in two commonly used polymers: polyethylene (PE) and isotactic polypropylene (iPP). Amorphous samples for both polymers were prepared by means of Monte Carlo and molecular dynamics simulations, and the oxygen solubility was measured across different temperatures. The solubility-temperature dependence for iPP proved to be nonmonotonic due to the interplay between binding and reorganizational enthalpy, while for PE, it appeared to be monotonic based on the available data in the studied temperature range. A broad comparison to experiments and simulations is included. Further oxygen insertions were also performed in semicrystalline PE and iPP samples at ambient temperature, and the obtained values were compared to a linear relationship which correlates the solubility in the purely amorphous phase with the solubility in the crystalline phase. The solubility of PE closely follows the linear relationship, while iPP exhibits some divergence. All the semicrystalline samples were previously annealed at elevated temperatures for long periods (a few μs), and a strong effect of annealing was observed on the structure and the solubility of iPP. A well-developed iPP lamellar structure emerged at longer annealing times, while PE develops that structure already in the early crystallization stages. The solubility of semicrystalline iPP samples with lamellar morphology exhibited better agreement with extrapolated solubility values of the amorphous state─the extrapolation was made using a linear relationship connecting solubility in the purely amorphous phase and solubility in mixed phases (amorphous and crystalline). Results on the correlation of the solubility with the local structural ordering are also present.

Driven Self-Assembly of Patchy Particles Overcoming Equilibrium Limitations.

Bridging biological complexity and synthetic material design, we investigate dissipative self-assembly in patchy particle systems. Utilizing Monte Carlo and Molecular Dynamics simulations, we demonstrate how external driving forces mitigate equilibrium trade-offs between assembly time and structural stability, traditionally encountered in self-assembly processes. Our findings also extend to biological-mimicking environments, where we explore the dynamics of patchy particles under crowded conditions. This comprehensive analysis offers insights into advanced material design, opening avenues for innovations in nanotechnology applications.

Removal of tricyclazole residues from water by sorption on metal organic frameworks (MOFs): A theoretical insight of the experimental data

All over the world, there are stringent regulations in place concerning the acceptable amount of pesticide residues in various food commodities, including cereals. The majority of the water tainted with pesticides comes from agricultural drainage, making it imperative that this water be cleansed of any lingering traces of the chemicals. The presence of unacceptable levels of pesticides residues is another issue hindering India's ability to export Basmati rice and other commodities. Tricyclazole, a fungicide has been found in rice in recent years, leading to the rejection of a number of bulk consignments. Therefore, three metal organic frameworks (MOFs) viz., MIL-53(Al), Fe-BTC, and ZIF-8, based on three different metal ions (aluminium, iron, and zinc) and three different organic "struts" were investigated for the removal of TRCZ from water. Interestingly, TRCZ was found to be selective towards both metal ions and linkers. The maximum adsorption capacities of TRCZ onto MIL-53(Al), Fe-BTC and ZIF-8 were 980.0, 896.7 and 904.6 μg g−1 (@1 ppm aqueous solution), respectively. Adsorption followed pseudo-second-order kinetics and best fit the Sips isotherm. The adsorption of TRCZ molecules onto the surface of the MOFs were confirmed using Fourier Transform infrared (FTIR) and X-Ray Diffraction (XRD) spectroscopy. MIL-53(Al) showed the best regeneration capacity (> 90 % efficiency after the 4th cycle). The adsorption of pesticide molecules onto MOF surfaces was studied (in vacuo and in aqueous) circumstances using Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) simulations revealing the prime competitive role of water molecules aided by π-π stacking and π-lone pair interactions during adsorption. Metal organic frameworks have the potential to replace costly commercial activated carbons in the future to treat pesticide-contaminated water.

Impact of aminopropyl units on physisorption in mesostructured silica is larger for CH4 than CO2 and depends on pore curvature

Mesoporous silica-based materials attract great interest for diverse applications, including capture, storage, and catalysis due to their large surface area, uniform pore size distribution, and customizable surface. The regular 2D hexagonal pore symmetry of MCM-41 makes it an ideal model. Amine groups are introduced to enhance CO2 adsorption, but the interplay of factors influencing overall performance, including pore size, remains unclear. This computational study explores the physisorption of CO2 and CH4, shedding light on the microscopic details underlying adsorption capacity and selectivity. Three aminopropyl-functionalized MCM-41 models with different pore size are compared through Grand Canonical Monte Carlo and Molecular Dynamics simulations. Aminopropyl chains affect CO2 adsorption, but its affinity for the preserved silanol groups is larger. In addition, the competitive interaction between amine and silanol groups emerged, and it is pore curvature dependent. Surprisingly, CH4 physisorption is influenced even more by surface functionalization. Results show a non-monotonic trend as a function of pore size. Local density of the alkyl chains plays the major role. Overall, the intermediate pore size demonstrates the best performance, thanks to the balance among the various effects. The study emphasizes the importance of understanding the interplay between material texture and gas interactions to optimize adsorption performance.

Effect of surface roughness on methane adsorption in shale organic nanopores from the perspective of molecular simulations

Understanding the shale gas adsorption behavior within organic nanopores is crucial for its development. The pore surfaces commonly exhibit random inherent roughness, which can significantly impact gas adsorption, while it’s still in debate about whether the roughness could enhance the gas adsorption or not. To clarify the underlying mechanisms of the effect of surface roughness, in this study, the Grand Canonical Monte Carlo and Molecular Dynamic simulations are used to investigate the methane adsorption in wrinkled graphene nanochannels with various roughness. The heterogeneity of the adsorption layer along the rough topology is described in detail. The adsorption is enhanced near the concave region and weakened near the convex region. Additionally, increased surface roughness shifts the gas molecules’ aggregation modes in the concave region from a gradual gradient to oscillatory “adsorption points”. We observe that wall roughness due to the wrinkle can enhance the adsorption amount overall, attributed to the increased surface area and the increased unit adsorption capacity due to potential energy overlap. However, the contributions of both factors differ under various pressure and roughness conditions. Then a new mathematical method to estimate the gas adsorption amounts is provided and aligns well with molecular simulation results for a wide range of pressure and roughness conditions. And a potential answer the current debate about the positive or negative impact of adsorption due to roughness is promoted. The key may lie in the dominant role of the co-existed roughness caused by the fragment basal or edge surface. Our work could shed light on the underlying effect of surface roughness on the gas adsorption in shale, and have significant implications for the reserve estimation of shale gas formations.

In-depth study of a newly synthesized imidazole derivative as an eco-friendly corrosion inhibitor for mild steel in 1 M HCl: Theoretical, electrochemical, and surface analysis perspectives

The present study is concerned with the corrosion inhibition and adsorption behaviour of a series of novel imidazole derivatives. The compounds under investigation are 5,5-diphenyl-3-propyl-2-(propylthio)-3,5-dihydro-4 H-imidazol-4-one (AM3) and 3-allyl-2-(allylthio)-5,5-diphenyl-3,5-dihydro-4 H-imidazol-4-one (AM6). The objective of this study is to evaluate the efficacy of 5,5-diphenyl-3,5-dihydro-4 H-imidazol-4-one (AM6) as a corrosion inhibitor on mild steel immersed in a 1.0 M hydrochloric acid medium. This comprehensive study assesses the efficacy of these derivatives through a range of electrochemical and spectroscopic analysis techniques. Additionally, polarisation curves, electrochemical impedance spectroscopy and advanced computer simulations were employed to evaluate the efficacy and inhibition mechanism of imidazole derivatives, thereby facilitating a more profound understanding of their anticorrosive capacity. The results obtained from potentiodynamic polarisation (PDP), electrochemical frequency modulation (EFM) and electrochemical impedance spectroscopy (EIS) measurements demonstrate that the inhibition efficiency increases with increasing imidazole derivative concentration. Conversely, an inverse relationship is observed between inhibition efficiency and temperature. The thermodynamic parameters ΔG°_ads and ΔH°_ads corroborate the conclusion that the adsorption process is predominantly chemical in nature. Scanning electron microscopy (SEM), energy dispersive X-ray spectroscopy (EDX) and X-ray photoelectron spectroscopy (XPS) were employed to characterise the surface morphology. The maximum protection afforded by imidazole derivatives was 95.2 % and 93 % for compounds AM3 and AM6, respectively. Polarisation curves indicated that imidazole derivatives exhibited mixed inhibition behaviour. Surface analysis demonstrated that imidazole derivatives were effectively adsorbed onto the carbon surface, thereby significantly reducing acid damage. This finding was corroborated by DFT calculations, as well as Monte Carlo (MC) and molecular dynamics (MD) simulations. These simulations provided a comprehensive insight into the adsorption of imidazole and its protonated form onto the carbon surface, offering valuable insight into the corrosion inhibition mechanism.

Damage kinetics in high-temperature irradiated Ni crystals

High-temperature (∼800 K) ion irradiation with 380 keV Ar ions was used to modify the surface of Ni single crystals in order to generate radiation defects and mimic the effect of neutrons avoiding sample activation. The modified 800 nm thick layer was analyzed regarding structural and mechanical properties utilizing a combination of experimental and simulation techniques. The ion implantation fluence ranged from 1 × 1014 cm−2 up to 1 × 1016 cm−2, which corresponds to values of displacements per atom from 0.25 to 25, respectively. The structural characterization was performed by Rutherford Spectroscopy in Channeling mode (RBS/C) associated with scanning (SEM) and transmission electron microscopy (TEM). The experimental results were supported by Monte Carlo (McChasy) and Molecular Dynamics (MD-LAMMPS) simulations. The simulations performed using the second generation of the McChasy code show that the influence of bubbles formed inside the material is not negligible and affects the quantitative analysis of defects in the simulated spectra. A second step in damage kinetics was revealed for the highest dose (i.e., 25 dpa) as a consequence of the entanglement of dislocations and agglomeration of noble gas bubbles, which was confirmed by TEM analysis. Moreover, nanomechanical results confirm that Ar bubbles deteriorate mechanical properties such as hardness.

Thermophysical Properties and Phase Behavior of CO2 with Impurities: Insight from Molecular Simulations.

Experimentally determining thermophysical properties for various compositions commonly found in CO2 transportation systems is extremely challenging. To overcome this challenge, we performed Monte Carlo (MC) and Molecular Dynamics (MD) simulations of CO2 rich mixtures to compute thermophysical properties such as densities, thermal expansion coefficients, isothermal compressibilities, heat capacities, Joule-Thomson coefficients, speed of sound, and viscosities at temperatures of (235-313) K and pressures of (20-200) bar. We computed thermophysical properties of pure CO2 and CO2 rich mixtures with N2, Ar, H2, and CH4 as impurities of (1-10) mol % and showed good agreement with available Equations of State (EoS). We showed that impurities decrease the values of thermal expansion coefficients, isothermal compressibilities, heat capacities, and Joule-Thomson coefficients in the gas phase, while these values increase in the liquid and supercritical phases. In contrast, impurities increase the value of speed of sound in the gas phase and decrease it in the liquid and supercritical phases. We present an extensive data set of thermophysical properties for CO2 rich mixtures with various impurities, which will help to design the safe and efficient operation of CO2 transportation systems.

Alkaloid extract of seed Citrullus colocynthis as novel green inhibitor for mild steel corrosion in one molar HCl acid solution: DFT and MC/MD approaches

The study was designed to explore the corrosion prevention capabilities of Citrullus colocynthis seeds alkaloid-rich extract (CSEA) on MS in a 1 M HCl environment by use of electrochemical and theoretical methods. Notably, Electrochemical Impedance Spectroscopy (EIS) and potentiodynamic polarization were used to probe the impact of immersion time and temperature. Energy-dispersive X-ray spectroscopy (EDX) and spanning electron microscopy (SEM) were used to confirm the presence of a protective layer on the substrate surface. Density functional theory (DFT) method was used to optimize the investigated species by use of B3LYP/6–31 + G(d, p) level of theory. The global and local quantum chemical reactivity descriptors were calculated to explore the inhibition corrosion efficiency and to identify the most favorable sites for electrophilic and nucleophilic attacks. Monte Carlo (MC) and molecular dynamics simulation (MDS) methods were used to study the interactions between corrosion inhibitor and metal surface. Noteworthy, results showed that CSEA exhibited an impressive inhibition efficiency, which reached 94.3% with a concentration of 2 g/L at 298 K. Potentiodynamic polarization revealed that the extract acted as a mixed-type inhibitor. Nyquist graphs showed that charge-transfer resistance and dæouble-layer capacitance both rised with increasing CSEA concentration, suggesting better inhibition efficiency. Notably, the Langmuir adsorption isotherm is well-aligned with the adsorption of inhibitor compounds. Importantly, all aforementioned theoretical methods were in agreement with the experimental findings. The outcome of the present work supported using Citrullus colocynthis seeds alkaloid-rich extract as ecofriendly agents to prevent corrosion.

Effect of molecular structure on the dynamics and viscoelasticity of vitrimers

Vitrimers are a class of dynamic polymer networks that perform bond-exchange reactions (BERs) under environmental stimuli while preserving network connectivity. The viscoelastic properties of vitrimers are not only dependent on the kinetics of BERs, but also the architectures of polymer matrix. Here, we investigate the influence of polymer architectures on the dynamics and linear viscoelasticity of vitrimers through varying the sticker (reactive bead) distribution in the precursor chain and the BERs energy barrier (Ub) by performing hybrid molecular dynamics-Monte Carlo simulations. Increasing Ub, vitrimers with block and gradient distributions exhibit considerably faster relaxations of chain ends, Rouse mode, and stress and lower zero-shear viscosity than the analogs with random and uniform distributions, in qualitative agreement with the predictions of the sticky Rouse model. In contrast, the distribution of stickers has a small influence on the relaxation of dynamic bonds regardless of the Ub value. By comparing the molecular structures of vitrimers with different distributions of stickers, we find that the molecular defects, especially the relatively longer dangling end in block and gradient distributions, impair the effect of dynamic bonds on the dynamics and viscoelasticity of vitrimers. This suggests that the dangling defects can be leveraged for further development of recyclable and healable materials with fast stress relaxation and improved flowability.

Dynamical investigation of NinAgm (n + m = 147, 309, 561) nanoalloys with core–shell orderings

ABSTRACT The structures and dynamical properties of core–shell bimetallic Ni-Ag nanoalloys varying with different sizes and compositions have been studied using Monte Carlo (MC) and Molecular Dynamic (MD) simulations. We have considered the compositions in which the size of the core increases while the total number of atoms is fixed. In this sense, two (Ni13Ag134, Ni55Ag92), three (Ni13Ag296, Ni55Ag254, Ni147Ag162) and four (Ni13Ag548, Ni55Ag506, Ni147Ag414 and Ni309Ag252) compositions were considered for 147, 309 and 561 atoms, respectively. Highly symmetric Mackay icosahedral structures with centred symmetric cores appear for these specific sizes and compositions. Also, smaller Ni atoms tend to occupy the core and Ag atoms prefer to segregate to the surface of the nanoalloy due to its lower surface and cohesive energy. Then, the lowest energy structures obtained by Basin Hopping MC simulations were used as initial configurations for melting simulations. The transitions between different chemical ordering patterns with increasing temperature are possible in these systems while they are still in the solid state. Although there are clear differences in the melting process of the compositions with increasing size of the core, for all cases, surface melting occurs indicating that the Ag shell melts before the inner Ni core.

Adsorption and Diffusion Properties of Gas in Nanopores of Kerogen: Insights from Grand Canonical Monte Carlo and Molecular Dynamics Simulations

Investigating the adsorption and diffusion processes of shale gas within the nanopores of kerogen is essential for comprehending the presence of shale gas in organic matter of shale. In this study, an organic nanoporous structure was constructed based on the unit structure of Longmaxi shale kerogen. Grand canonical Monte Carlo and molecular dynamics simulation methods were employed to explore the adsorption and diffusion mechanisms of pure CH4, CO2, and N2, as well as their binary mixtures with varying mole fractions. The results revealed that the physical adsorption characteristics of CH4, CO2, and N2 gases on kerogen adhered to the Langmuir adsorption law. The quantity of adsorbed gas molecules increased with rising pressure but decreased with increasing temperature. The variation in the heat of adsorption was also analyzed. Under identical temperature and pressure conditions, the adsorption of CH4 increased with higher mole fractions of CH4, whereas it decreased with greater mole fractions of CO2 and N2. Notably, CO2 molecules exhibited a robust interaction with kerogen molecules compared to the adsorption properties of CH4 and N2. Furthermore, the self-diffusion coefficient of gas within kerogen nanopores gradually decreased with increasing pressure or decreasing temperature. The diffusion capacity of gas molecules followed the descending order N2 > CH4 > CO2 under the same pressure and temperature conditions.

Modeling segregated solutes in plastically deformed alloys using coupled molecular dynamics-Monte Carlo simulations

A microscopic understanding of the complex solute-defect interaction is pivotal for optimizing the alloy’s macroscopic mechanical properties. Simulating solute segregation in a plastically deformed crystalline system at atomic resolution remains challenging. The objective is to efficiently model and predict a physically informed segregated solute distribution rather than simulating a series of diffusion kinetics. To address this objective, we coupled molecular dynamics (MD) and Monte Carlo (MC) methods using a novel method based on virtual atoms technique. We applied our MD-MC coupling approach to model off-lattice carbon (C) solute segregation in nanoindented Fe-C samples containing complex dislocation networks. Our coupling framework yielded the final configuration through efficient parallelization and localized energy computations, showing C Cottrell atmospheres near dislocations. Different initial C concentrations resulted in a consistent trend of C atoms migrating from less crystalline distortion to high crystalline distortion regions. Besides unraveling the strong spatial correlation between local C concentration and defect regions, our results revealed two crucial aspects of solute segregation preferences: (1) defect energetics hierarchy and (2) tensile strain fields near dislocations. The proposed approach is generic and can be applied to other material systems as well.

Memristive feature and mechanism induced by laser-doping in defect-free 2D semiconductor materials

Memristors as non-volatile memory devices have gained numerous attentions owing to their advantages in storage, in-memory computing, synaptic applications, etc. In recent years, two-dimensional (2D) materials with moderate defects have been discovered to exist memristive feature. However, it is very difficult to obtain moderate defect degree in 2D materials, and studied on modulation means and mechanism becomes urgent and essential. In this work, we realized memristive feature with a bipolar switching and a configurable on/off ratio in a two-terminal MoS2 device (on/off ratio ~100), for the first time, from absent to present using laser-modulation to few-layer defect-free MoS2 (about 10 layers), and its retention time in both high resistance state and low resistance state can reach 2 × 104 s. The mechanism of the laser-induced memristive feature has been cleared by dynamic Monte Carlo simulations and first-principles calculations. Furthermore, we verified the universality of the laser-modulation by investigating other 2D materials of TMDs. Our work will open a route to modulate and optimize the performance of 2D semiconductor memristive devices.

Physisorption or chemisorption: Insight from AI computing model based on DFT, MC/MD-simulation for prediction of MOF-based inhibitor adsorption on Cu in brine solution

This study investigates the dynamics of cerium metal organic framework (Ce-MOF) adsorption on the Cu surface in brine solution. This research integrates computational modelling based on Density Functional Theory (DFT), Monte Carlo (MC), and Molecular Dynamics (MD) simulations for enhanced predictive modelling. Additionally, by the utilization of the Dmol3 module within Biovia materials studio software, DFT calculations /geometric configurations through the Double Numerical plus d-functions (DND) were achieved. Notably, the results demonstrate a robust negative adsorption energy (−363.54 kcal/mol), indicating a highly stable interaction leading to the formation of a protective film on the metal surface. The study also proposes a corrosion inhibition mechanism involving physical and chemical barriers. The findings not only contribute to fundamental corrosion science but also propose innovative strategies for designing MOF-based inhibitors tailored to specific environmental conditions. This research introduces a novel perspective by combining theoretical calculations, and simulation models for predicting inhibitor behavior, fostering advancements in corrosion protection strategies.

Molecular Dynamics Simulation of Methane Adsorption and Diffusion

To study the adsorption and diffusion mechanism of low-rank coal at the microscopic level, samples of Fukang low-rank coal were collected, and the elemental composition, carbon type distribution and functional group type of the Fukang low-rank coal structure were determined by elemental analysis (Ea), Fourier-transform interferometric radiometer (FTIR), X-ray photoelectron spectroscopy (XPS) and 13C nuclear magnetic resonance (13C NMR) experiments to construct a 2D molecular structure of the coal and a 3D macromolecular structure model. The adsorption and diffusion characteristics of methane were researched by giant regular Monte Carlo (GCMC) and molecular dynamics (MD) simulation methods. The results showed that the excess adsorption amount of methane increased and then decreased with the increase in pressure. The diffusion of methane showed two stages with increasing pressure: a sharp decrease in the diffusion coefficient from 0.5 to 5.0 MPa and a slow decrease in the diffusion coefficient from 5.0 to 15.0 MPa. The lower the pressure, the larger the effective radius of the CH4 and C atoms, and the higher the temperature, the more pronounced the diffusion and the larger the effective radius.

Refinement of Kelvin equation for condensation of water in CSH slits based on molecular simulation

The adsorption and condensation of water molecules in the pores of cementitious materials play a crucial role in the mechanical performance and durability of concrete. The Kelvin equation is generally employed to describe the water condensation in porous materials. However, it requires modification for condensation behavior in nano and micro scales of pores. In the present work, various widths of slit models for several cementitious minerals were established. The adsorption isotherms in slit pores were calculated using the hybrid of Grand Canonical Monte Carlo and Molecular Dynamics simulations, obtaining the relationship between the critical pressures for condensation and pore sizes. Then, a modified Kelvin equation was proposed, with discussion on a new reasonable method for determining the real thickness of adsorbed water film. The results demonstrate that condensation will occur in the slits within the critical width range, which promotes the surface adsorption to achieve complete adsorption. Therefore, the location of the interface between the adsorbed water film and condensed water bulk is fixed in the same mineral during condensation. The modified Kelvin equation shows higher accuracy than traditional Kelvin equations at the nano and micro scale. This work promotes a deeper understanding of the interaction mechanism of water with cementitious minerals and inspires more fundamental realization in the design of concrete durability.