Computational investigation of adsorption of molecular hydrogen on lithium-doped corannulene

Density functional theory and classical molecular dynamics simulations are used to investigate the prospect of lithium-doped corannulene as adsorbent material for H(2) gas. Potential energy surface scans at the level of B3LYP/6-311G(d,p) show an enhanced interaction of molecular hydrogen with lithium-atom-doped corannulene complexes with respect to that found in undoped corannulene. MP2(FC)/6-31G(d,p) optimizations of 4H(2)-(Li(2)-C(20)H(10)) yield H(2) binding energies of -1.48 kcal/mol for the H(2)-Li interaction and -0.92 kcal/mol for the H(2)-C interaction, whereas values of -0.94 and -0.83 kcal/mol were reported (J. Phys. Chem. B 2006, 110, 7688-7694) for physisorption of H(2) on the concave and the convex side of corannulene using MP2(full)/6-31G(d), respectively. Classical molecular dynamics simulations predict hydrogen uptakes in Li-doped corannulene assemblies that are significantly enhanced with respect to that found in undoped molecules, and the hydrogen uptake ability is dependent on the concentration of lithium dopant. For the Li(6)-C(20)H(10) complex, a hydrogen uptake of 4.58 wt % at 300 K and 230 bar is obtained when the adsorbent molecules are arranged in stack configurations separated by 6.5 A, and with interlayer distances of 10 A, hydrogen uptake reaches 6.5 wt % at 300 K and 215 bar.     

Interaction of Hydrogen with MOF-5

Hydrogen storage is among the most demanding challenges in the hydrogen-based energy cycle. One proposed strategy for hydrogen storage is based on physisorption on high surface area solids such as metal-organic frameworks (MOFs). Within this class of materials, MOF-5 has been the first structure studied for hydrogen storage. The IR spectroscopy of adsorbed H2 performed at 15 K and ab initio calculations show that the adsorptive properties of this material are mainly due to dispersive interactions with the internal wall structure and to weak electrostatic forces associated with O13Zn4 clusters. Calculated and measured binding enthalpies are between 2.26 and 3.5 kJ/mol, in agreement with the H2 rotational barriers reported in the literature. A minority of binding sites with higher adsorption enthalpy (7.4 kJ/mol) is also observed. These species are probably associated with OH groups on the external surfaces present as termini of the microcrystals.     

Hydrogen adsorption on carbon-doped boron nitride nanotube

The adsorption of atomic and molecular hydrogen on carbon-doped boron nitride nanotubes is investigated within the ab initio density functional theory. The binding energy of adsorbed hydrogen on carbon-doped boron nitride nanotube is substantially increased when compared with hydrogen on nondoped nanotube. These results are in agreement with experimental results for boron nitride nanotubes (BNNT) where dangling bonds are present. The atomic hydrogen makes a chemical covalent bond with carbon substitution, while a physisorption occurs for the molecular hydrogen. For the H(2) molecule adsorbed on the top of a carbon atom in a boron site (BNNT + C(B)-H(2)), a donor defect level is present, while for the H(2) molecule adsorbed on the top of a carbon atom in a nitrogen site (BNNT + C(N)-H(2)), an acceptor defect level is present. The binding energies of H(2) molecules absorbed on carbon-doped boron nitride nanotubes are in the optimal range to work as a hydrogen storage medium.     

Fine speciation of active sites in zeolites by a CO probe: Dynamics and IR frequencies

The influence of system dynamics onto infrared spectra calculated and measured for probe molecules weakly interacting with active sites is discussed. OH stretching frequencies for Brønsted sites in Y zeolite and their shift on CO adsorption (measured by IR spectroscopy at ca. 170 K) are simulated from AIMD trajectory, based on BO DFT calculations. We present here very good agreement between experimental and computed OH frequency red‐shift on CO adsorption in dealuminated zeolite Y. It is shown that extracting nonharmonic OH frequencies from ab initio molecular dynamics trajectory provides better estimate of experimental IR results than standard vibrational analyses due to realistic computational model and better correspondence with experimental conditions.     

氢在沸石上的吸附行为

应用基于Ono-Kondo格子理论得到的通用吸附等温方程, 通过分析氢在不同温度下, 在沸石NaX、CaA、NaA和ZSM-5上的吸附数据, 确定了氢的最大单层吸附容量. 并引入维里吸附方程, 由第二维里吸附系数和圆柱孔的Lennard-Jones(12-6)势模型计算了氢与沸石微孔壁面的作用势. 结果表明, 通用吸附等温方程可较好地描述氢在沸石上的超临界吸附行为, 拟合所得的氢在沸石上的最大单层吸附容量与吸附剂相关, 而与吸附温度无关. 圆柱孔作用势模型计算所得的氢分子在沸石上的吸附作用势与吸附热相近. 氢分子间的作用力表现为吸引力.     

Influence of quantum effects on the physisorption of molecular hydrogen in model carbon foams

The physisorption of molecular hydrogen in model carbon foams has been investigated from 50 K to room temperature. The study is carried out within the framework of the density functional theory for quantum liquids at finite temperatures. Calculations are performed in the grand canonical ensemble, i.e., the adsorbed fluid is assumed to be in equilibrium with an external gas of hydrogen molecules with concentrations ranging from 8×10(-4) kg m(-3) to n=71 kg m(-3). It is shown that, while strong zero-point energy effects are present even at room temperature, the adsorption isotherms exhibit only a weak dependence on the explicit incorporation of the bosonic exchange symmetry of hydrogen molecules. The increase of the average particle density prevents the deviations from the Maxwell-Boltzmann statistics to become noticeable if the system is cooled down. The volumetric storage capacity of these materials at low temperatures is about one half of the U. S. Department of Energy goal, while the gravimetric capacity is still far from the standards required by mobile applications. The relation between the microscopic structure of the hydrogen fluid and the calculated adsorption properties is also addressed.     

Modified Langmuir Isotherm for the Description of Cluster Adsorption on Surface Lyophilic Centers

An adsorption isotherm that generalizes the Langmuir equation for the cluster adsorption of several interacting molecules at one center is proposed. In the case of the adsorption of water molecules at surface hydroxyl groups on modified silica it was shown that the energy and structure characteristics of such a system, determined by the ab initio method, agree with the experimental data and lead to satisfactory agreement with the parameters of the proposed isotherm. __________ Translated from Teoreticheskaya i Eksperimental'naya Khimiya, Vol. 41, No. 5, pp. 283–289, September–October, 2005.     

Structure, energetics, and dynamics of water adsorbed on the muscovite (001) surface: a molecular dynamics simulation

Molecular dynamics (MD) computer simulations of liquid water adsorbed on the muscovite (001) surface provide a greatly increased, atomistically detailed understanding of surface-related effects on the spatial variation in the structural and orientational ordering, hydrogen bond (H-bond) organization, and local density of H2O molecules at this important model phyllosilicate surface. MD simulations at constant temperature and volume (statistical NVT ensemble) were performed for a series of model systems consisting of a two-layer muscovite slab (representing 8 crystallographic surface unit cells of the substrate) and 0 to 319 adsorbed H2O molecules, probing the atomistic structure and dynamics of surface aqueous films up to 3 nm in thickness. The results do not demonstrate a completely liquid-like behavior, as otherwise suggested from the interpretation of X-ray reflectivity measurements and earlier Monte Carlo simulations. Instead, a more structurally and orientationally restricted behavior of surface H2O molecules is observed, and this structural ordering extends to larger distances from the surface than previously expected. Even at the largest surface water coverage studied, over 20% of H2O molecules are associated with specific adsorption sites, and another 50% maintain strongly preferred orientations relative to the surface. This partially ordered structure is also different from the well-ordered 2-dimensional ice-like structure predicted by ab initio MD simulations for a system with a complete monolayer water coverage. However, consistent with these ab initio results, our simulations do predict that a full molecular monolayer surface water coverage represents a relatively stable surface structure in terms of the lowest diffusional mobility of H2O molecules along the surface. Calculated energies of water adsorption are in good agreement with available experimental data.     

Experimental and computational thermochemical study of 2- and 3-thiopheneacetic acid methyl esters

Thiophene-based compounds have widespread use in modern drug design, biodiagnostics, electronic and optoelectronic devices, and conductive polymers. The present study reports an experimental and computational thermochemical study on the relative stabilities of 2- and 3-thiopheneacetic acid methyl esters. The enthalpies of combustion and vaporization were measured by a rotating-bomb combustion calorimeter, Calvet microcalorimetry, and correlation gas chromatography, and the gas-phase enthalpies of formation at T=298.15 K were determined. Standard ab initio molecular orbital calculations at the G3 level were performed, and a theoretical study of the molecular and electronic structure of the compounds studied was carried out. Calculated enthalpies of formation, using atomization and isodesmic reactions are in very good agreement with the experimental results.     

A search for a strong physisorption site for H2 in Li-doped porous carbons

The mechanism of hydrogen absorption between two coronene molecules has been studied by first principle calculations. Examination of different sites for H(2) molecule confirmed the classical picture of physisorption. We have also considered molecular hydrogen adsorption in a charged carbon structure achieved by doping with lithium at a density corresponding to the intercalate compound LiC(6). We have performed different types of calculations [Hartree-Fock and density functional theory (DFT)] for various atomic basis sets using CRYSTAL98, GAUSSIAN98, and DMOL3 codes. B3LYP-DFT (B3LYP-three-parameter functional of Backe, Lee, Yang and Parr) energy minimization calculations unravel that there is a stable adsorption site for molecular hydrogen in Li-doped sp(2) carbon structure. These calculations also give an insight into the atomic configurations of interlayer species (H(2) and Li) as the interlayer spacing increases. It can be shown that large changes occur in the positions and electronic properties of interlayer species. Hydrogen molecule does not show any tendency for dissociation and adopts a position in the interlayer void that is deeply related to that of lithium ions. We have evidenced a rather large charge transfer from lithium and capping hydrogen species on neighboring slab carbon atoms that induce the stabilization of molecular hydrogen. We have also found that rotating one carbon layer with respect to the other one (at constant interlayer distance) does not change the adsorption energy to a large extent. The best adsorption site is about five times deeper than the physisorption site found in the undoped case and occurs at an interlayer separation of 5.5+/-0.5 A. The corresponding atomic configuration consists in a hydrogen molecule standing (nearly) perpendicular to the plane surface surrounded by the three lithium ions in a configuration close to that of the LiC(6) intercalation compound.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Electrostatic model calculations on multiple adsorption at NaCl surfaces

A previously proposed electrostatic model for physisorption at ionic solids is extended to multiple adsorption of small molecules. A new algorithm is developed to avoid interpenetration of the interacting systems. The geometry optimization procedure is described. Ab initio calculations are used for the application of the method to the adsorption of CO and CO₂ at NaCl(100) surfaces simulated by Na₂₅Cl₂₅ clusters. The orientation of the adsorbate molecules in dependence on the cumulative atomic multipole moments (CAMMs) of the cluster atoms is discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 685–693, 1997     

Enthalpy and entropy effects in hydrogen adsorption on carbon nanotubes

Interaction energies and entropies associated with hydrogen adsorption on the inner and outer surfaces of zigzag single-wall carbon nanotubes (SWCNT) of various diameters are analyzed by means of molecular mechanics, density functional theory, and ab initio calculations. For a single molecule the strongest interaction, which is 3.5 greater than that with the planar graphite sheet, is found inside a (8,0) nanotube. Adsorption on the outer surfaces is weaker than that on graphite. Due to the steric considerations, both processes are accompanied by an extremely strong decline in entropy. Absence of specific adsorption sites and weak attractive interaction between hydrogen molecules within carbon nanotubes results in their close packing at low temperatures. Using the calculated geometric and thermodynamic parameters in Langmuir isotherms we predict the adsorption capacity of SWCNTs at room temperature to be smaller than 1 wt % even at 100 bar.     

分子模拟预测含配位键不饱和Cu的有机骨架(MOFs)材料的储氢性能

本文根据RI-UMP2方法、TZVPP基组和BSSE校正计算得到的能量数据,推导了描述氢分子与含有不饱和配位键的Cu的金属羧酸配合物相互作用的分子力学力场. 用巨正则系综蒙特卡洛模拟(GCMC)计算了氢在含四方形配位Cu的MOFs材料上的吸附等温线. 通过对比CuBTC的实验吸附数据发现,虽然理论预测含有四方形配位Cu的MOFs具有比含有四面体配位Zn的MOFs与氢分子更强的相互作用,实验合成得到的材料尚未反映这一区别. 以现有的CuBTC实验数据为参照消除理论计算和实验测定的系统偏差,预测了3种含有四方形配位Cu的MOFs材料以CuBTC作为参照的储氢能力.     

Thermochemistry of biphenylcarboxylic and dicarboxylic acids. A combined experimental and theoretical study

The standard molar enthalpies of combustion and sublimation of 2- and 4-biphenylcarboxylic acid, 2,2'- and 4,4'-biphenyldicarboxylic acid were measured and the gas-phase enthalpies of formation, at T= 298.15 K, were determined. Ab initio calculations were performed and a theoretical study on molecular structure of all the biphenyl acid isomers has been carried out. Calculated enthalpies of formation using appropriate isodesmic reactions are compared with experimental values, and a good agreement is observed. Estimates of enthalpies of formation for the isomers, which were not studied experimentally, are presented. All the acids containing at least one ortho COOH are comparatively less stable than their isomers having just meta or para COOH group(s).     

Adsorption mechanism of water molecules surrounding Au nanoparticles of different sizes

Molecular dynamic simulation is used to investigate the adsorption mechanism of water molecules surrounding Au nanoparticles with different sizes. Our results show that the adsorption mechanism of the water molecules in the first water shell will be influenced by the size of the Au nanoparticle. For the larger Au nanoparticles, the hydrogen bonding of water molecules adsorbed on the surface of the Au nanoparticles are arranged in a two-dimensional structure, while those adsorbed on the edge of the surface of the Au nanoparticles are arranged in a three-dimensional structure. However, in the case of the smallest Au nanoparticle, the hydrogen bonding of the water molecules on the first adsorbed layer are arranged only in a three-dimensional structure. The arrangement of the water molecules in the first water shell can be determined by orientation order parameter. The water molecules that adsorb on the larger Au nanoparticles tend to arrange in an irregular arrangement, while those adsorbed on the smallest Au nanoparticle tend to arrange a regular arrangement. Interestingly, the water molecules adsorbed on the smallest nanoparticle are arranged in a bulklike structure in the first shell.     

DFT studies of the interactions of a graphene layer with small water aggregates

We have investigated the structure, adsorption, electronic states, and charge transfer of small water aggregates on the surface of a graphene layer using density functional theory. Our calculations were focused on water adsorbates containing up to five water molecules interacting with one and both sides of a perfect freestanding sheet. Different orientations of the aggregates with respect to the graphene sites were considered. The results show that the adsorption energy of one water molecule is primarily determined by its orientation, although it is also strongly dependent on the implemented functional scheme. Despite its intrinsic difficulties with dispersion interactions, the Perdew and Wang's exchange-correlation functional may be a viable alternative to investigate the adsorption of large molecular aggregates on a graphene surface. Although water physisorption is expected to occur in the regime of droplets, we found no induced impurity states close to the Fermi level of graphene interacting with small water clusters. In order to investigate the donor/acceptor tendency of the water clusters on graphene, we have performed a Bader charge analysis. Considering the charge transfer mechanism, we have noticed that it should preferentially occur from water to graphene only when the oxygen atom is pointing toward the surface. Otherwise, and in the case of larger adsorbed clusters, charge transfers systematically occur from graphene to water.     

Calculation of enthalpies of hydrogenation of hydrocarbons

Sets of hydrogen molecule equivalents have been developed which permit the calculation of hydrogenation of different types of carbon-carbon bonds from ab initio total energies (3-21G and 6-31G* basis sets, and, to a more limited extent, for MP2/6-31G* data) of reactants and products. The calculated enthalpies of hydrogenation are in good agreement with experiment for unstrained molecules, with average errors on the order of 2 kcal/mol. The 6-31G* equivalents allow the enthalpies for strained molecules to be calculated accurately, but the 3-21G equivalents do not. The equivalents for both basis sets have been tested by calculating the enthalpies of hydrogenation of carbon-carbon bonds in nitrogen- and oxygen-containing organic molecules, free radicals, and classical carbocations. The results are in good agreement with experiment in most cases.     

含羧酸镁(钙)官能团的多孔芳香骨架材料储氢性能的预测

用MP2方法,TZVPP基组以及基组重叠误差(BSSE)校正计算了氢分子与修饰在多孔芳香骨架(PAF)上的羧酸镁、羧酸钙官能团的相互作用,并建立了描述这一相互作用的分子力学力场.在此基础上用巨正则系综蒙特卡洛(GCMC)模拟预测了氢气在该种新型PAF材料上的吸附等温线.量子化学计算结果表明,每个羧酸镁、羧酸钙官能团分别可以提供13、14个氢分子吸附位点,与每个氢分子的平均结合能在8kJ·mol-1左右.通过比较不同温度和压力下材料的绝对吸附量和超额吸附量发现,在PAF骨架中引入羧酸镁、羧酸钙官能团可以显著提高材料的综合储氢性能,达到并超过了美国能源部提出的2015年储氢标准.同时该工作还揭示了氢吸附量与材料的表面积、空腔体积和分子作用强度间的复杂关系.     

First-principle study of adsorption of hydrogen on Ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H2 molecules, respectively). Additionally, a molecular adsorption state of H2 above the Ti atom is observed for the first time and is attributed to the polarization of the H2 molecule by the Ti cation. Our results parallel recent findings for H2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.