Adsorption of hydrogen in covalent organic frameworks using expanded Wang–Landau simulations

Accelerated convergence via adiabatic sampling for adsorption and desorption processes

Under isothermal conditions, phase transitions occur through a nucleation event when conditions are sufficiently close to coexistence. The formation of a nucleus of the new phase requires the system to overcome a free energy barrier of formation, whose height rapidly rises as supersaturation decreases. This phenomenon occurs both in the bulk and under confinement and leads to a very slow kinetics for the transition, ultimately resulting in hysteresis, where the system can remain in a metastable state for a long time. This has broad implications, for instance, when using simulations to predict phase diagrams or screen porous materials for gas storage applications. Here, we leverage simulations in an adiabatic statistical ensemble, known as adiabatic grand-isochoric ensemble (μ, V, L) ensemble, to reach equilibrium states with a greater efficiency than its isothermal counterpart, i.e., simulations in the grand-canonical ensemble. For the bulk, we show that at low supersaturation, isothermal simulations converge slowly, while adiabatic simulations exhibit a fast convergence over a wide range of supersaturation. We then focus on adsorption and desorption processes in nanoporous materials, assess the reliability of (μ, V, L) simulations on the adsorption of argon in IRMOF-1, and demonstrate the efficiency of adiabatic simulations to predict efficiently the equilibrium loading during the adsorption and desorption of argon in MCM-41, a system that exhibits significant hysteresis. We provide quantitative measures of the increased rate of convergence when using adiabatic simulations. Adiabatic simulations explore a wide temperature range, leading to a more efficient exploration of the configuration space.

Design and applications of three dimensional covalent organic frameworks

We summarize in this review the current state-of-the-art development of three dimensional covalent organic frameworks.

Classical and quantum many-body effects on the critical properties and thermodynamic regularities of silicon

Using molecular simulation, we determine the critical properties of Si as well as the loci for several remarkable thermodynamic contours spanning the supercritical region of the phase diagram. We consider a classical three-body potential as well as a quantum (tight-binding) many-body model, and determine the loci for the ideality contours, including the Zeno line and the H line of ideal enthalpy. The two strategies (classical or quantum) lead to strongly asymmetric binodals and to critical properties in good agreement with each other. The Zeno and H lines are found to remain linear over a wide temperature interval, despite the changes in electronic structure undergone by the fluid along these contours. We also show that the classical and quantum model yield markedly different results for the parameters defining the H line, the exponents for the power-laws underlying the line of minima for the isothermal enthalpy and for the density required to achieve ideal behavior, most notably for the enthalpy.

Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. IV. Performance of many-body force fields and tight-binding schemes for the fluid phases of silicon

We extend Expanded Wang-Landau (EWL) simulations beyond classical systems and develop the EWL method for systems modeled with a tight-binding Hamiltonian. We then apply the method to determine the partition function and thus all thermodynamic properties, including the Gibbs free energy and entropy, of the fluid phases of Si. We compare the results from quantum many-body (QMB) tight binding models, which explicitly calculate the overlap between the atomic orbitals of neighboring atoms, to those obtained with classical many-body (CMB) force fields, which allow to recover the tetrahedral organization in condensed phases of Si through, e.g., a repulsive 3-body term that favors the ideal tetrahedral angle. Along the vapor-liquid coexistence, between 3000 K and 6000 K, the densities for the two coexisting phases are found to vary significantly (by 5 orders of magnitude for the vapor and by up to 25% for the liquid) and to provide a stringent test of the models. Transitions from vapor to liquid are predicted to occur for chemical potentials that are 10%-15% higher for CMB models than for QMB models, and a ranking of the force fields is provided by comparing the predictions for the vapor pressure to the experimental data. QMB models also reveal the formation of a gap in the electronic density of states of the coexisting liquid at high temperatures. Subjecting Si to a nanoscopic confinement has a dramatic effect on the phase diagram with, e.g. at 6000 K, a decrease in liquid densities by about 50% for both CMB and QMB models and an increase in vapor densities between 90% (CMB) and 170% (QMB). The results presented here provide a full picture of the impact of the strategy (CMB or QMB) chosen to model many-body effects on the thermodynamic properties of the fluid phases of Si.

Impact of Friedel oscillations on vapor-liquid equilibria and supercritical properties in two and three dimensions

We determine the impact of the Friedel oscillations on the phase behavior, critical properties, and thermodynamic contours in films [two dimensions (2D)] and bulk phases [three dimensions (3D)]. Using expanded Wang-Landau simulations, we calculate the grand-canonical partition function and, in turn, the thermodynamic properties of systems modeled with a linear combination of the Lennard-Jones and Dzugutov potentials, weighted by a parameter X (0<X<1). Varying X allows us to control the height of the first Friedel oscillation and to provide a complete characterization of the effect of the metal-like character in the potential on the thermodynamic properties over a wide range of conditions. For 3D systems, we are able to show that the critical parameters exhibit a linear dependence on X and that the loci for the thermodynamic state points, for which the system shows the same compressibility factor or enthalpy as an ideal gas, are two straight lines spanning the subcritical and supercritical regions of the phase diagram for all X values. Reducing the dimensionality to 2D results in a loss of impact of the Friedel oscillation on the critical properties, as evidenced by the virtually constant critical density across the range of X values. Furthermore, our results establish that the straightness of the two ideality lines is retained in 2D and is independent from the height of the first Friedel oscillation in the potential.

Scaling Laws and Critical Properties for fcc and hcp Metals

The determination of the critical parameters of metals has remained particularly challenging both experimentally, because of the very large temperatures involved, and theoretically, because of the many-body interactions that take place in metals. Moreover, experiments have shown that these systems exhibit an unusually strong asymmetry of their binodal. Recent theoretical work has led to new similarity laws, based on the calculation of the Zeno line and of the underlying Boyle parameters, which provided results for the critical properties of atomic and molecular systems in excellent agreement with experiments. Using the recently developed expanded Wang-Landau (EWL) simulation method, we evaluate the grand-canonical partition function, over a wide range of conditions, for 11 fcc and hcp metals (Ag, Al, Au, Be, Cu, Ir, Ni, Pb, Pd, Pt, and Rh), modeled with a many-body interaction potential. This allows us to calculate the binodal, Zeno line, and Boyle parameters and, in turn, obtain the critical properties for these systems. We also propose two scaling laws for the enthalpy and entropy of vaporization, and identify critical exponents of 0.4 and 1.22 for these two laws, respectively.