A comparative study of nitrogen physisorption on different C70 crystal structures using an ab initio based potential

Quantum mechanical calculations are performed using the recently developed hybrid method for interaction energies to determine atom site Lennard-Jones potential parameters for the interactions of molecular nitrogen with C(70) molecules. This ab initio based potential is used in grand canonical Monte Carlo simulations to predict surface adsorption properties of N(2) on five known C(70) structures: rhombohedral, fcc, ideal hcp, deformed hcp, and monoclinic crystals. Because of the presence of five-membered carbon rings and the surface curvature of C(70) molecule, the Lennard-Jones potential parameters for nitrogen-carbon interactions obtained from ab initio based calculations are found to be different from that with planar graphite. The simulation results obtained from these two sets of force fields are compared and shown to differ, particularly at low coverage, where the nitrogen-carbon interactions are more important than the nitrogen-nitrogen interactions. The surface area, monolayer capacity, and isosteric heat of adsorption are calculated for various C(70) crystals and found to change appreciably as a result of the shear-induced phase transformation from hcp to rhombohedral lattice.     

Calibration of force-field dependency in free energy landscapes of peptide conformations by quantum chemical calculations

The free energy landscapes of peptide conformations were calibrated by ab initio quantum chemical calculations, after the enhanced conformational diversity search using the multicanonical molecular dynamics simulations. Three different potentials of mean force for an isolated dipeptide were individually obtained by the multicanonical molecular dynamics simulations using the conventional force fields, AMBER parm94, AMBER parm96, and CHARMm22. Each potential of mean force was then calibrated based upon the umbrella sampling algorithm from the adiabatic energy map that was calculated separately by the ab initio molecular orbital method, and all of the calibrated potentials of mean force coincided well. The calibration method was also applied to the simulations of a peptide dimer in explicit water models, and it was shown that the calibrated free energy landscapes did not depend on the force field used in the classical simulations, as far as the conformational space was sampled well. The current calibration method fuses the classical free energy calculation with the quantum chemical calculation, and it should generally make simulations for biomolecular systems much more reliable when combining with enhanced conformational sampling.     

Free energy perturbation study of water dimer dissociation kinetics

An efficient approach is described for using accurate ab initio calculations to determine the rates of elementary condensation and evaporation processes that lead to nucleation of aqueous aerosols. The feasibility of the method is demonstrated in an application to evaporation rates of water dimer at 230 K. The method, known as ABC-FEP (ab initio/classical free energy perturbation), begins with a calculation of the potential of mean force for the dissociation (evaporation) of small water clusters using a molecular dynamics (MD) simulation with a model potential. The free energy perturbation is used to calculate how changing from the model potential to a potential calculated from ab initio methods would alter the potential of mean force. The difference in free energy is the Boltzmann-weighted average of the difference between the ab initio and classical potential energies, with the average taken over a sample of configurations from the MD simulation. In principle, the method does not require a highly accurate model potential, though more accurate potentials require fewer configurations to achieve a small sampling error in the free energy perturbation step. To test the feasibility of obtaining accurate potentials of mean force from ab initio calculations at a modest number of configurations, the free energy perturbation method has been used to correct the errors when some standard models for bulk water (SPC, TIP4P, and TIP4PFQ) are applied to water dimer. To allow a thorough exploration of sampling issues, a highly accurate fit to results of accurate ab initio calculations, known as SAPT-5s, as been used a proxy for the ab initio calculations. It is shown that accurate values for a point on the potential of mean force can be obtained from any of the water models using ab initio calculations at only 50 configurations. Thus, this method allows accurate simulations of small clusters without the need to develop water models specifically for clusters.     

Multilayer adsorption of water at a rutile TiO2)(110) surface: towards a realistic modeling by molecular dynamics

We present a model combining ab initio concepts and molecular dynamics simulations for a more realistic treatment of complex adsorption processes. The energy, distance, and orientation of water molecules adsorbed on stoichiometric and reduced rutile TiO(2)(110) surfaces at 140 K are studied via constant temperature molecular dynamics simulations. From ab initio calculations relaxed atomic geometries for the surface and the most probable adsorption sites were derived. The study comprises (i) large two-dimensional surface supercells, providing a realistically low concentration of surface oxygen defects, and (ii) a water coverage sufficiently large to model the onset of the growth of a bulk phase of water on the surface. By our combined approach the influence of both, the metal oxide surface, below, and the bulk water phase, above, on the water molecules forming the interface between the TiO(2) surface and the water bulk layer is taken into account. The good agreement of calculated adsorption energies with experimental temperature programmed desorption spectra demonstrates the validity of our model.     

The kinetic study of specific adsorption of phosphate species on Pt(111) in acidic solutions

The study of the adsorption/desorption mechanism of phosphate anions at Pt(111) in acidic solution of pH 4.3 and 0.8 was performed by the potential step method in order to reveal the kinetics of anion adsorption. The current-time curve due to phosphate adsorption/desorption showed various decay features, being dependent on the potential region. The rate of current decay depended on pH, being faster in a lower pH solution. Specific adsorption processes were analyzed by the Langmuir and Elovich adsorption equations and also in terms of a two-dimensional nucleation-growth mechanism in different adsorption/desorption regions. In the case of adsorption in 0.3M phosphate buffer solution of pH 4.3, random adsorption without interaction following the Langmuir adsorption, takes place at low coverage, while random adsorption with repulsive force was observed at high coverage. In the desorption process, random desorption with repulsive force takes place at high coverage, and the repulsive force disappears where random adsorption without interaction takes place at medium coverage. When the surface coverage becomes further lower, the desorption mechanism changes dramatically into a two-dimensional nucleation-growth type, suggesting that an ordered adsorbate structure is formed after a rapid discharge process of anion adsorption.     

Numerical interpretation of molecular surface field in dielectric modeling of solvation

Continuum solvent models, particularly those based on the Poisson‐Boltzmann equation (PBE), are widely used in the studies of biomolecular structures and functions. Existing PBE developments have been mainly focused on how to obtain more accurate and/or more efficient numerical potentials and energies. However to adopt the PBE models for molecular dynamics simulations, a difficulty is how to interpret dielectric boundary forces accurately and efficiently for robust dynamics simulations. This study documents the implementation and analysis of a range of standard fitting schemes, including both one‐sided and two‐sided methods with both first‐order and second‐order Taylor expansions, to calculate molecular surface electric fields to facilitate the numerical calculation of dielectric boundary forces. These efforts prompted us to develop an efficient approximated one‐dimensional method, which is to fit the surface field one dimension at a time, for biomolecular applications without much compromise in accuracy. We also developed a surface‐to‐atom force partition scheme given a level set representation of analytical molecular surfaces to facilitate their applications to molecular simulations. Testing of these fitting methods in the dielectric boundary force calculations shows that the second‐order methods, including the one‐dimensional method, consistently perform among the best in the molecular test cases. Finally, the timing analysis shows the approximated one‐dimensional method is far more efficient than standard second‐order methods in the PBE force calculations. © 2017 Wiley Periodicals, Inc.     

Quantum chemical prediction of the adsorption conformations and dynamics at HCOOH‐covered ZnO(1010) surfaces

Results from ab initio Hartree–Fock and gradient‐corrected density functional theory calculations of formic acid interactions with ZnO (101 0) surfaces are reported. Surface relaxation is found to affect equilibrium geometries and adsorption energies significantly. Large variations in adsorption energy with coverage and ordering of the adsorbates are revealed and explained in terms of strong and highly anisotropic electrostatic adsorbate–adsorbate interactions. The results are compared to published experimental and theoretical results, and differences in suggested binding geometries from the different studies are discussed. Dynamic properties of the adsorption, surface mobility, and surface reactivity are inferred from key elements of the potential energy surface obtained from the quantum chemical computations and supported by ab initio molecular dynamics simulations. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Investigation of ligand exchange reactions in aqueous uranyl carbonate complexes using computational approaches

Carbonate anion exchange reactions with water in the uranyl-carbonate and calcium-uranyl-carbonate aqueous systems have been investigated using computational methods. Classical molecular dynamics (MD) simulations with the umbrella sampling technique were employed to determine potentials of mean force for the exchange reactions of water and carbonate. The presence of calcium counter-ions is predicted to increase the stability of the uranyl-carbonate species in accordance with previous experimental observations. However, the free energy barrier to carbonate exchange with water is found to be comparable both in the presence and absence of calcium cations. Possible implications of these results for uranyl adsorption on mineral surfaces are discussed. Density functional theory (DFT) calculations were also used to confirm the trends observed in classical molecular dynamics simulations and to corroborate the validity of the potential parameters employed in the MD scheme.     

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.     

New force field for calcium binding sites in annexin-membrane complexes

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.     

Adsorption Activity and Molecular Dynamics Study on Anti-corrosion Mechanism of Q235 Steel

The correlation between inhibition efficiency and molecular structures of the inhibitor during hydrochloric acid corrosion of Q235 steel was studied by quantum chemical calculations and molecular dynamics(MD) simulation. The proton affinity(PA) calculations demonstrated that 2-(quinolin-2-yl)quinazolin-4(3H)-one inhibitor has the tendency to be protonated in hydrochloric acid, which was in good agreement with experimental observations. Besides, quantum chemical parameters revealed that the protonated corrosion inhibitor molecules were more easily adsorbed on Q235 steel surface and improved the corrosion resistance of steel. MD simulations were implemented to search for the adsorption behavior of this molecule on Fe (110) surface, which might be used as a convenient tool for estimating the interaction mechanism between inhibitor and iron surface.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Interaction between benzenedithiolate and gold: classical force field for chemical bonding

We have constructed a group of classical potentials based on ab initio density-functional theory (DFT) calculations to describe the chemical bonding between benzenedithiolate (BDT) molecule and gold atoms, including bond stretching, bond angle bending, and dihedral angle torsion involved at the interface between the molecule and gold clusters. Three DFT functionals, local-density approximation (LDA), PBE0, and X3LYP, have been implemented to calculate single point energies (SPE) for a large number of molecular configurations of BDT-1, 2 Au complexes. The three DFT methods yield similar bonding curves. The variations of atomic charges from Mulliken population analysis within the molecule/metal complex versus different molecular configurations have been investigated in detail. We found that, except for bonded atoms in BDT-1, 2 Au complexes, the Mulliken partial charges of other atoms in BDT are quite stable, which significantly reduces the uncertainty in partial charge selections in classical molecular simulations. Molecular-dynamics (MD) simulations are performed to investigate the structure of BDT self-assembled monolayer (SAM) and the adsorption geometry of S adatoms on Au (111) surface. We found that the bond-stretching potential is the most dominant part in chemical bonding. Whereas the local bonding geometry of BDT molecular configuration may depend on the DFT functional used, the global packing structure of BDT SAM is quite independent of DFT functional, even though the uncertainty of some force-field parameters for chemical bonding can be as large as approximately 100%. This indicates that the intermolecular interactions play a dominant role in determining the BDT SAMs global packing structure.     

A combined molecular dynamics+quantum mechanics method for investigation of dynamic effects on local surface structures

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation.     

Interaction between water molecules and zinc sulfide nanoparticles studied by temperature-programmed desorption and molecular dynamics simulations

We have investigated the bonding of water molecules to the surfaces of ZnS nanoparticles (approximately 2-3 nm sphalerite) using temperature-programmed desorption (TPD). The activation energy for water desorption was derived as a function of the surface coverage through kinetic modeling of the experimental TPD curves. The binding energy of water equals the activation energy of desorption if it is assumed that the activation energy for adsorption is nearly zero. Molecular dynamics (MD) simulations of water adsorption on 3 and 5 nm sphalerite nanoparticles provided insights into the adsorption process and water binding at the atomic level. Water binds with the ZnS nanoparticle surface mainly via formation of Zn-O bonds. As compared with bulk ZnS crystals, ZnS nanoparticles can adsorb more water molecules per unit surface area due to the greatly increased curvature, which increases the distance between adjacent adsorbed molecules. Results from both TPD and MD show that the water binding energy increases with decreasing the water surface coverage. We attribute the increase in binding energy with decreasing surface water coverage to the increasing degree of surface under-coordination as removal of water molecules proceeds. MD also suggests that the water binding energy increases with decreasing particle size due to the further distance and hence lower interaction between adsorbed water molecules on highly curved smaller particle surfaces. Results also show that the binding energy, and thus the strength of interaction of water, is highest in isolated nanoparticles, lower in nanoparticle aggregates, and lowest in bulk crystals. Given that water binding is driven by surface energy reduction, we attribute the decreased binding energy for aggregated as compared to isolated particles to the decrease in surface energy that occurs as the result of inter-particle interactions.     

Ab initio protein structure prediction with force field parameters derived from water-phase quantum chemical calculation

Molecular dynamics (MD) simulations are extensively used in the study of the structures and functions of proteins. Ab initio protein structure prediction is one of the most important subjects in computational biology, and many trials have been performed using MD simulation so far. Since the results of MD simulations largely depend on the force field, reliable force field parameters are indispensable for the success of MD simulation. In this work, we have modified atom charges in a standard force field on the basis of water-phase quantum chemical calculations. The modified force field turned out appropriate for ab initio protein structure prediction by the MD simulation with the generalized Born method. Detailed analysis was performed in terms of the conformational stability of amino acid residues, the stability of secondary structure of proteins, and the accuracy for prediction of protein tertiary structure, comparing the modified force field with a standard one. The energy balance between alpha-helix and beta-sheet structures was significantly improved by the modification of charge parameters.     

O2 interaction and reactivity on a model hydroxylated rutile(110) surface

Recently several theoretical studies have examined oxygen adsorption on the clean, reduced TiO2(110) surface. However the photocatalytic behavior of TiO2 and the scavenging ability of oxygen are known to be influenced by the presence of surface hydroxyls. In this paper the chemistry of O2 on the hydroxylated TiO2 surface is investigated by means of first-principles total energy calculations and molecular dynamics (MD) simulations. The MD trajectories show a direct, spontaneous reaction between O2 and the surface hydroxyls, thus supporting the experimental hypothesis that the reaction does not necessarily pass through a chemisorbed O2 state. Following this reaction, the most stable chemisorbed intermediates are found to be peroxide species HO2 and H2O2. Although these intermediates are very stable on the short time scale of MD simulations, the energetics suggests that their further transformation is connected to a new 300 K feature observed in the experimental water temperature programmed desorption (TPD) spectrum. The participation of two less stable intermediate states, involving terminal hydroxyls and/or chemisorbed water plus oxygen adatoms, to the desorption process, is not supported by the total energy calculations. Analysis of the projected density of states, however, suggests the possibility that these intermediates have a role in completing the surface oxidation immediately before desorption.     

High-pressure in situ 129Xe NMR spectroscopy and computer simulations of breathing transitions in the metal-organic framework Ni2(2,6-ndc)2(dabco) (DUT-8(Ni))

Recently, we have described the metal-organic framework Ni(2)(2,6-ndc)(2)(dabco), denoted as DUT-8(Ni) (1) (DUT = Dresden University of Technology, 2,6-ndc = 2,6-naphthalenedicarboxylate, dabco = 1,4-diazabicyclo[2.2.2]octane). Upon adsorption of molecules such as nitrogen and xenon, this material exhibits a pronounced gate-pressure effect which is accompanied by a large change of the specific volume. Here, we describe the use of high-pressure in situ (129)Xe NMR spectroscopy, i.e., the NMR spectroscopic measurements of xenon adsorption/desorption isotherms and isobars, to characterize this effect. It appears that the pore system of DUT-8(Ni) takes up xenon until a liquid-like state is reached. Deeper insight into the interactions between the host DUT-8(Ni) and the guest atom xenon is gained from ab initio molecular dynamics (MD) simulations. van der Waals interactions are included for the first time in these calculations on a metal-organic framework compound. MD simulations allow the identification of preferred adsorption sites for xenon as well as insight into the breathing effect at a molecular scale. Grand canonical Monte Carlo (GCMC) simulations have been performed in order to simulate adsorption isotherms. Furthermore, the favorable influence of a sample pretreatment using solvent exchange and drying with supercritical CO(2) as well as the influence of repeated pore opening/closure processes, i.e., the "aging behavior" of the compound, can be visualized by (129)Xe NMR spectroscopy.     

Adsorption of polycyclic aromatic hydrocarbons at the air-water interface: molecular dynamics simulations and experimental atmospheric observations

Adsorption of benzene, naphthalene, anthracene, and phenanthrene at the aqueous surface is investigated by means of molecular dynamics simulations. Potentials of mean force, i.e., free energy profiles obtained when moving the studied molecules across an aqueous slab were evaluated. In all cases, deep surface free energy minima, corresponding to orders of magnitude of surface enhancement of the aromatic molecule, were located. This enhancement, which increases with the size of the solute, points to the importance of the aqueous surface for the chemistry of polycyclic aromatic hydrocarbons (PAHs). Supporting evidence in the atmospheric environment related to the heterogeneous chemistry of PAHs on water droplets and planar surfaces is summarized. There is good agreement between the hydration free energies computed from MD calculations and the experimentally determined values. Data pertaining to the importance of the air-water interface in the adsorption and transport of PAHs on micron sized water droplets are described. The relevant data on adsorption and reaction (ozonation and photochemical) at the air-water interface of planar surfaces and droplets are also summarized.