The atomistic molecular dynamics program YASP has been parallelized for shared-memory computer architectures. Parallelization was restricted to the most CPU-time-consuming parts: neighbor-list construction, calculation of nonbonded, angle and dihedral forces, and constraints. Most of the sequential FORTRAN code was kept; parallel constructs were inserted as compiler directives using the OpenMP standard. Only in the case of the neighbor list did the data structure have to be changed. The parallel code achieves a useful speedup over the sequential version for systems of several thousand atoms and above. On an IBM Regatta p690+, the throughput increases with the number of processors up to a maximum of 12-16 processors depending on the characteristics of the simulated systems. On dual-processor Xeon systems, the speedup is about 1.7.
Molecular dynamics simulations are used to study the microscopic structure and dynamics of cations bound to cucurbit[6]uril (CB[6]) in water and in aqueous solutions of sodium, potassium, and calcium chloride. The molarities are 0.183 M for the salts and 0.0184 M for CB[6]. The cations bind only to CB[6] carbonyl oxygens. They are never found inside the CB[6] cavity. Complexes with Na(+) and K(+) mostly involve one cation, whereas with Ca(2+) single- and double-cation complexes are formed in similar proportions. The binding dynamics strongly depends on the type of cation. A smaller size or higher charge increases the residence time of a cation at a given carbonyl oxygen. When bound to CB[6], sodium and potassium cations jump mainly between nearest or second-nearest neighbors. Calcium shows no hopping dynamics. It is coordinated predominantly by one CB[6] oxygen. A few water molecules (zero to four) can occupy the CB[6] cavity, which is limited by the CB[6] oxygen faces. Their residence time is hardly influenced by sodium and potassium ions. In the case of calcium the residence time of the inner water increases notably. A simple structural model for the cation activity as "lids" over the CB[6] portal cannot, however, be identified. The slowing of the water exchange by the ions is a consequence of the generally slower dynamics in their presence and of their stable solvation shells.
A molecular dynamics model and its parametrization procedure are devised and used to study adsorption of isopropanol on platinum(111) (Pt(111)) surface in unsaturated and oversaturated coverages regimes. Static and dynamic properties of the interface between Pt(111) and liquid isopropanol are also investigated. The magnitude of the adsorption energy at unsaturated level increases at higher coverages. At the oversaturated coverage (multilayer adsorption) the adsorption energy reduces, which coincides with findings by Panja et al. in their temperature-programed desorption experiment [Surf. Sci. 395, 248 (1998)]. The density analysis showed a strong packing of molecules at the interface followed by a depletion layer and then by an oscillating density profile up to 3 nm. The distribution of individual atom types showed that the first adsorbed layer forms a hydrophobic methyl "brush." This brush then determines the distributions further from the surface. In the second layer methyl and methine groups are closer to the surface and followed by the hydroxyl groups; the third layer has exactly the inverted distribution. The alternating pattern extends up to about 2 nm from the surface. The orientational structure of molecules as a function of distance of molecules is determined by the atom distribution and surprisingly does not depend on the electrostatic or chemical interactions of isopropanol with the metal surface. However, possible formation of hydrogen bonds in the first layer is notably influenced by these interactions. The surface-adsorbate interactions influence the mobility of isopropanol molecules only in the first layer. Mobility in the higher layers is independent of these interactions.
MD studies of liquid isopropyl alcohol and melts of short poly(vinyl alcohol) (PVA) oligomers are described. The specific volume was found to depend inversely on the number N of repeat units. If the chain length is enhanced, the viscosity of the PVA melt increases and the peaks in the radial distribution function become sharper. Additional peaks that appear in melts of PVA chains are of pure intramolecular origin. The calculated radius of gyration was found to depend on the number of formula units via $N^{0.65 \pm 0.03}$. The orientation correlation functions showed that all molecular vectors of PVA melts with chain lengths N = 1, 2, 3 relax completely within a few nanoseconds. The relaxation times for the OH bond vector as obtained via the Kohlrausch‐Williams‐Watts expression showed an exponential dependence on the number of repeat units.magnified image
A Molecular dynamics simulations are used to study the structure and dynamics of the platinum(111)/poly(vinyl alcohol) (PVA) oligomers (10 monomers long) interface at 400 K. The mass density and number density distribution of separate atoms along the surface normal resemble, in general, the distributions obtained for the platinum(111)/ liquid isopropanol interface (28). The small discrepancies are dictated by the chemical bonds between monomers within PVA chains. The differences between PVA and isopropanol in the orientational structure of O-H and O-C bond vectors are large immediately at the platinum surface. The connectivity of monomers within PVA chains is also the main driving force of these changes. At longer distances, the structure of the PVA melt resembles approximately that of liquid isopropnaol. The PVA chains, which are directly adsorbed on the surface, extend into the bulk up to 2.25 nm. There exists a region where adsorbed chains mix with nonadsorbed ones and, therefore, form very likely hydrogen bonds between each other. Thus, it is expected that PVA melt, unlike liquid isopropanol, attaches relatively strongly onto the metal surface.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.