LL

L. Lin

info

Please Note

7 records found

Journal article (2017) - Sayed Jamali, Thijs Vlugt, Li-Chiang Lin
Reverse osmosis constitutes a large portion of currently operating commercial water desalination systems. Employing membranes with large water fluxes while maintaining high salt rejection is of central importance in decreasing the associated energy consumption and costs. The ultrathin-film nature of zeolite nanosheets and their versatile pore structures provides great opportunities in desalination. To push forward the development of zeolite nanosheets for water desalination, nonequilibrium molecular dynamics simulations were carried out to systematically study zeolites as RO membranes and establish fundamental structure-performance relationships. We have identified that zeolite nanosheets can achieve a high salt rejection rate close to 100% while allowing nearly 2 orders of magnitude higher water permeability than currently available membranes. Moreover, the effects of the pore density, inclusion of cages, and free energy barrier on water permeability and salt rejection are unraveled, leading to important insights toward the rational design of novel zeolite membranes. ...
Journal article (2017) - Tim M. Becker, J. Heinen, David Dubbeldam, Li Chiang Lin, Thijs J H Vlugt
(Graph Presented) The family of M-MOF-74, with M = Co, Cr, Cu, Fe, Mg, Mn, Ni, Ti, V, and Zn, provides opportunities for numerous energy related gas separation applications. The pore structure of M-MOF-74 exhibits a high internal surface area and an exceptionally large adsorption capacity. The chemical environment of the adsorbate molecule in M-MOF-74 can be tuned by exchanging the metal ion incorporated in the structure. To optimize materials for a given separation process, insights into how the choice of the metal ion affects the interaction strength with adsorbate molecules and how to model these interactions are essential. Here, we quantitatively highlight the importance of polarization by comparing the proposed polarizable force field to orbital interaction energies from DFT calculations. Adsorption isotherms and heats of adsorption are computed for CO2, CH4, and their mixtures in M-MOF-74 with all 10 metal ions. The results are compared to experimental data, and to previous simulation results using nonpolarizable force fields derived from quantum mechanics. To the best of our knowledge, the developed polarizable force field is the only one so far trying to cover such a large set of possible metal ions. For the majority of metal ions, our simulations are in good agreement with experiments, demonstrating the effectiveness of our polarizable potential and the transferability of the adopted approach. ...
Journal article (2016) - Zhaochuan Fan, Li-Chiang Lin, Wim Buijs, Thijs Vlugt, MA van Huis
Cation exchange is a powerful tool for the synthesis of nanostructures such as core–shell nanocrystals, however, the underlying mechanism is poorly understood. Interactions of cations with ligands and solvent molecules are systematically ignored in simulations. Here, we introduce the concept of pseudoligands to incorporate cation-ligand-solvent interactions in molecular dynamics. This leads to excellent agreement with experimental data on cation
exchange of PbS nanocrystals, whereby Pb ions are partially replaced by Cd ions from solution. The temperature and the ligand-type control the exchange rate and equilibrium composition of cations in the nanocrystal. Our simulations reveal that Pb ions are kicked out by exchanged Cd interstitials and migrate through interstitial sites, aided by local relaxations at core–shell interfaces and point defects. We also predict that high-pressure conditions facilitate strongly enhanced cation exchange reactions at elevated temperatures. Our approach is easily extendable to other semiconductor compounds and to other families of
nanocrystals. ...
Journal article (2016) - David Cohen-Tanugi, Li Chiang Lin, Jeffrey C. Grossman
While single-layer nanoporous graphene (NPG) has shown promise as a reverse osmosis (RO) desalination membrane, multilayer graphene membranes can be synthesized more economically than the single-layer material. In this work, we build upon the knowledge gained to date toward single-layer graphene to explore how multilayer NPG might serve as a RO membrane in water desalination using classical molecular dynamic simulations. We show that, while multilayer NPG exhibits similarly promising desalination properties to single-layer membranes, their separation performance can be designed by manipulating various configurational variables in the multilayer case. This work establishes an atomic-level understanding of the effects of additional NPG layers, layer separation, and pore alignment on desalination performance, providing useful guidelines for the design of multilayer NPG membranes. ...
Journal article (2016) - Amber Janda, Bess Vlaisavljevich, Li Chiang Lin, Berend Smit, Alexis T. Bell
The effects of zeolite structure on the kinetics of n-butane monomolecular cracking and dehydrogenation are investigated for eight zeolites differing in the topology of channels and cages. Monte Carlo simulations are used to calculate enthalpy and entropy changes for adsorption (ΔHads-H+ and ΔSads-H+) of gas-phase alkanes onto Brønsted protons. These parameters are used to extract intrinsic activation enthalpies (ΔHint ), entropies (ΔSint ), and rate coefficients (kint) from measured data. As ΔSads-H+ decreases (i.e., as confinement increases), ΔHint and ΔSint for terminal cracking and dehydrogenation decrease for a given channel topology. These results, together with positive values observed for ΔSint , indicate that the transition states for these reactions resemble products. For central cracking (an earlier transition state), ΔHint is relatively constant, while ΔSint increases as ΔSads-H+ decreases because less entropy is lost upon protonation of the alkane. Concurrently, selectivities to terminal cracking and dehydrogenation decrease relative to central cracking because ΔSint decreases for the former reactions. Depending on channel topology, changes in the measured rate coefficients (kapp) with confinement are driven by changes in kint or by changes in the adsorption equilibrium constant (Kads-H+). Values of ΔSint and ΔHint are positively correlated, consistent with weaker interactions between the zeolite and transition state and with the greater freedom of movement of product fragments within more spacious pores. These results differ from earlier reports that ΔHint and ΔSint are structure-insensitive and that kapp is dominated by Kads-H+. They also suggest that ΔSads-H+ is a meaningful descriptor of confinement for zeolites having similar channel topologies. ...
In this modeling study, the uses of nitrogen (77.3 K), probe molecule of choice for decades, and argon, opted as alternative in the 2015 IUPAC report on adsorptive characterization, as probe molecules for geometric surface area determination are compared. Graphene sheets possessing slit-shaped pores with varying size (width) are chosen as model porous solids, and different methods for the determination of specific surface areas are investigated. The BET method, which is the most commonly applied analysis, is compared to the Langmuir and relatively recently proposed ESW (excess sorption work) method. We show that either using argon or nitrogen as adsorptive, the physical meaningfulness of adsorption-derived surface areas highly depends on the pore size. When less than two full layers of adsorbate molecules can be formed within slitlike pores of a graphitic material (Dpore < 5.8 Å for Ar/N2), adsorption-derived surface areas are about half that of the geometric surface area. Between two and four layers (6.8 < Dpore < 12.8 Å), adsorption surface areas can be significantly larger (up to 75%) than the geometric surface area because monolayer−multilayer formation and pore filling cannot be distinguished. For four or more layers of adsorbate molecules (Dpore > 12.8 Å), adsorption-derived surface areas are comparable to their geometrically accessible counterparts. Note that for the Langmuir method this only holds if pore-filling effects are excluded during determination. This occurs in activated carbon materials as well. In the literature, this indistinguishability issue has been largely overlooked, and erroneous claims of materials with extremely large surface areas have been made. Both the BET and Langmuir areas, for Dpore > 12.8 Å, correspond to geometric surface areas, whereas the ESW method yields significantly lower values. For the 6.8 Å < Dpore < 12.8 Å range, all methods erroneously overestimate the specific surface area. For the energetically homogeneous graphene sheets, differences between argon and nitrogen for the assessment of surface areas are minor. ...